Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic Simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

A method for data handling numerical results in parallel OpenFOAM simulations

Parallel computational fluid dynamics simulations produce vast amount of numerical result data. This paper introduces a method for reducing the size of the data by replaying the interprocessor traffic. The results are recovered only in certain regions of interest configured by the user. A known test case is used for several mesh partitioning scenarios using the OpenFOAM toolkit{sup ®}[1]. The space savings obtained with classic algorithms remain constant for more than 60 Gb of floating point data. Our method is most efficient on large simulation meshes and is much better suited for compressing large scale simulation results than the regular algorithms.

42 ENGINEERING↗
DOE OSTI2015

Constraint methods that accelerate free-energy simulations of biomolecules

Atomistic molecular dynamics simulations of biomolecules are critical for generating narratives about biological mechanisms. The power of atomistic simulations is that these are physics-based methods that satisfy Boltzmann’s law, so they can be used to compute populations, dynamics, and mechanisms. But physical simulations are computationally intensive and do not scale well to the sizes of many important biomolecules. One way to speed up physical simulations is by coarse-graining the potential function. Another way is to harness structural knowledge, often by imposing spring-like restraints. But harnessing external knowledge in physical simulations is problematic because knowledge, data, or hunches have errors, noise, and combinatoric uncertainties. Here, we review recent principled methods for imposing restraints to speed up physics-based molecular simulations that promise to scale to larger biomolecules and motions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2011

Dynamic heat capacity of the east model and of a bead-spring polymer model.

In this report we have presented a brief review of the glass transition and one means of characterizing glassy materials: linear and nonlinear thermodynamic oscillatory experiments to extract the dynamic heat capacity. We have applied these methods to the east model (a variation of the Ising model for glass forming systems) and a simple polymeric system via molecular dynamics simulation, and our results match what is seen in experiment. For the east model, since the dynamics are so simple, a mathematical model is developed that matches the simulated dynamics. For the polymeric system, since the system is a simulation, we can instantaneously 'quench' the system - removing all vibrational energy - to separate the vibrational dynamics from dynamics associated with particle rearrangements. This shows that the long-time glassy dynamics are due entirely to the particle rearrangements, i.e. basin jumping on the potential energy landscape. Finally, we present an extension of linear dynamic heat capacity to the nonlinear regime.

36 MATERIALS SCIENCE↗
DOE OSTI2012

A FIRST ESTIMATE OF RADIO HALO STATISTICS FROM LARGE-SCALE COSMOLOGICAL SIMULATION

We present a first estimate based on a cosmological gas dynamics simulation of galaxy cluster radio halo counts to be expected in forthcoming low-frequency radio surveys. Our estimate is based on a FLASH simulation of the {Lambda}CDM model for which we have assigned radio power to clusters via a model that relates radio emissivity to cluster magnetic field strength, intracluster turbulence, and density. We vary several free parameters of this model and find that radio halo number counts vary by up to a factor of two for average magnetic fields ranging from 0.2 to 3.1 {mu}G. However, we predict significantly fewer low-frequency radio halos than expected from previous semi-analytic estimates, although this discrepancy could be explained by frequency-dependent radio halo probabilities as predicted in reacceleration models. We find that upcoming surveys will have difficulty in distinguishing models because of large uncertainties and low number counts. Additionally, according to our modeling we find that expected number counts can be degenerate with both reacceleration and hadronic secondary models of cosmic-ray generation. We find that relations between radio power and mass and X-ray luminosity may be used to distinguish models, and by building mock radio sky maps we demonstrate that surveys such as LOFAR may have sufficient resolution and sensitivity to break this model degeneracy by imaging many individual clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Adaptive two-regime method: Application to front propagation

The Adaptive Two-Regime Method (ATRM) is developed for hybrid (multiscale) stochastic simulation of reaction-diffusion problems. It efficiently couples detailed Brownian dynamics simulations with coarser lattice-based models. The ATRM is a generalization of the previously developed Two-Regime Method [Flegg et al., J. R. Soc., Interface 9, 859 (2012)] to multiscale problems which require a dynamic selection of regions where detailed Brownian dynamics simulation is used. Typical applications include a front propagation or spatio-temporal oscillations. In this paper, the ATRM is used for an in-depth study of front propagation in a stochastic reaction-diffusion system which has its mean-field model given in terms of the Fisher equation [R. Fisher, Ann. Eugen. 7, 355 (1937)]. It exhibits a travelling reaction front which is sensitive to stochastic fluctuations at the leading edge of the wavefront. Previous studies into stochastic effects on the Fisher wave propagation speed have focused on lattice-based models, but there has been limited progress using off-lattice (Brownian dynamics) models, which suffer due to their high computational cost, particularly at the high molecular numbers that are necessary to approach the Fisher mean-field model. By modelling only the wavefront itself with the off-lattice model, it is shown that the ATRM leads to the same Fisher wave results as purely off-lattice models, but at a fraction of the computational cost. The error analysis of the ATRM is also presented for a morphogen gradient model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Molecular simulation study of dynamical properties of room temperature ionic liquids with carbon pieces

Room temperature ionic liquids (RTILs) with dispersed carbon pieces exhibit distinctive physiochemical properties. In order to explore the molecular mechanism, RTILs/carbon pieces mixture we investigated it by molecular dynamics (MD) simulation in this work. Rigid and flexible carbon pieces in the form of graphene with different thicknesses and carbon nanotubes in different sizes were dispersed in a representative RTIL 1-butyl-3-methyl-imidazolium dicyanamide ([Bmim][DCA]). Our study demonstrated that the diffusion coefficients of RTILs in the presence of flexible carbons are similar to those of bulk RTILs at varying temperatures, which is in contrast to the decreased diffusion of RTILs in the presence of rigid carbons. In addition, interfacial ion number density at rigid carbon surfaces was higher than that at flexible ones, which is correlated with the accessible external surface area of carbon pieces. The life time of cation-anion pair in the presence of carbon pieces also exhibited a dependence on carbon flexibility. RTILs with dispersed rigid carbon pieces showed longer ion pair life time than those with flexible ones, in consistence with the observation in diffusion coefficients. Furthermore, this work highlights the necessity of including the carbon flexibility when performing MD simulation of RTILs in the presence of dispersed carbon pieces in order to obtain the reliable dynamical and interfacial structural properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Subdiffusive High-Pressure Hydrogen Gas Dynamics in Elastomers

We report elastomeric rubber materials serve a vital role as sealing materials in the hydrogen storage and transport infrastructure. With applications including O-rings and hose liners, these components are exposed to pressurized hydrogen at a range of temperatures, cycling rates, and pressure extremes. High-pressure exposure and subsequent rapid decompression often lead to cavitation and stress-induced damage of the elastomer due to localization of the hydrogen gas. Here, we use all-atom classical molecular dynamics simulations to assess the impact of compositional variations on gas diffusion within the commonly used elastomer ethylene–propylene–diene monomer (EPDM). With the aim to build a predictive understanding of precursors to cavitation and to motivate material formulations that are less sensitive to hydrogen-induced failure, we perform systematic simulations of gas dynamics in EPDM as a function of temperature, gas concentration, and cross-link density. Our simulations reveal anomalous, subdiffusive hydrogen motion at pressure and intermediate times. We identify two groups of gas with different mobilities: one group exhibiting high mobility and one group exhibiting low mobility due to their motion being impeded by the polymer. With decreasing temperatures, the low-mobility group shows increased gas localization, the necessary precursor for cavitation damage in these materials. At lower temperatures, increasing cross-link density led to greater hydrogen gas mobility and a lower fraction of caged hydrogen, indicating that increasing cross-link density may reduce precursors to cavitation. Finally, we use a two-state kinetic model to determine the energetics associated with transitions between these two mobility states.

36 MATERIALS SCIENCE↗
DOE OSTI2022

NPP Simulators for Coupled Thermal and Electric Power Dispatch

The Light Water Reactor Sustainability (LWRS) program within the United States Department of Energy supports extending the operation of the U.S. commercial nuclear power plant (NPP) fleet. Within the LWRS program, the Flexible Plant Operation and Generation (FPOG) Pathway works to diversify the revenue streams of light water reactors (LWRs) by opening opportunities for the co-generation of non-electric products in addition to supplying electrical power to the grid. Recent events have added greater motivation to these efforts. For example, the recent Inflation Reduction Act (IRA) passed by the U.S. federal government offers substantial tax incentives for producing clean hydrogen, the technology readiness level of dispatchable and high-efficiency hydrogen production has dramatically increased in a short time, and societal response to world climate change is driving a transition away from fossil fuels. Producing hydrogen with maximum efficiency using nuclear power requires dispatching both electrical and thermal power from the nuclear plant to the hydrogen plant, so testing concepts of operations for combined electrical and thermal power dispatch (TPD) from an NNP to a hydrogen plant is of interest. This report documents achievement of the Light Water Reactor Sustainability (LWRS) program milestone “Install and demonstrate a vendor-developed simulator on the Human Systems Simulation Laboratory (HSS) for dispatch of LWR electrical power to a close-coupled electrolysis plant” with a due date of Dec. 22, 2022. Several factors provide motivation for this effort. Coupling the power generation deck of a nuclear power plant to a hydrogen production facility introduces new possibilities for operational transients that must be addressed. In particular, the performance of the integrated system during startup and shutdown of the hydrogen production facility, as well as offnormal conditions, need to be evaluated to ensure there are no adverse effects on the operation of the existing NPP. The concept of operations involving the NPP, the hydrogen plant, and the electric power grid must be tested using NPP simulators and operating procedures that have been modified for TPD operations. These tests must also include dynamic simulations of the coupled tertiary thermal and electric loads as well as coordinated activities with NPP operators, tertiary load operators and grid power coordinators. The report summarizes progress in developing and testing full-scope NPP simulators at the HSSL, including a generic BWR simulator from GSE Systems, Inc. and generic PWR simulator from Westinghouse. In the case of the TPD-GBWR Simulator from GSE Systems, Inc., a BWR is thermally coupled to a high temperature electrolysis (HTE) plant that produces hydrogen and oxygen from de-ionized water. The hydrogen plant is not explicitly simulated but only included as a transient heat sink. A thermal power dispatch (TPD) system transfers heat between the steam systems at the BWR and the hydrogen plant. Operational results from two versions of the modified simulator are presented. The first version uses synthetic oil as a heat transfer fluid in a closed delivery heat loop (DHL) that generates steam at the hydrogen plant. The second version uses steam as the heat transfer fluid in a delivery steam line (DSL) to provide steam to the hydrogen plant. For both versions, the estimated thermal power delivery distance is approximately one kilometer. The amount of thermal power dispatched in the simulators is 15% of the total reactor thermal power such that the simulators provide a tool to study the feasibility of coupling a BWR to industrial processes that benefit from a combination electrical and thermal power dispatch. Ongoing work within a CRADA is also developing a full-scope PWR simulator provided by Westinghouse for both thermal and electric power coupling. This simulator is based on a PWR plant with two three-loop Westinghouse reactors. Westinghouse PWRs are sufficiently similar that a simulator of a three-loop reactor is an appropriate representation for two-loop and four-loop PWR reactors. The three-loop simulator will initially be modified for close-coupling to a 100 MW HTE hydrogen production plant that will require approximately 25 MW of thermal power while operating at its maximum rated capacity. The simulator testing will include full coupling to dynamic simulations of a hydrogen production plant and a representative bulk electric grid. The simulator provided by Westinghouse is similar to the GPWR simulator that INL has already obtained from GSE Systems but has a few important added benefits. First, the Westinghouse simulator is based on digital controls and has additional screens that can be called up to show parameter trends to assist operators in decision-making. The Westinghouse simulator also has upgrades to the controls and hardware representations, such as valve actuators, that make it more realistic and flexible in terms of accurately sim

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2017

Operator splitting method for simulation of dynamic flows in natural gas pipeline networks

Here, we develop an operator splitting method to simulate flows of isothermal compressible natural gas over transmission pipelines. The method solves a system of nonlinear hyperbolic partial differential equations (PDEs) of hydrodynamic type for mass flow and pressure on a metric graph, where turbulent losses of momentum are modeled by phenomenological Darcy-Weisbach friction. Mass flow balance is maintained through the boundary conditions at the network nodes, where natural gas is injected or withdrawn from the system. Gas flow through the network is controlled by compressors boosting pressure at the inlet of the adjoint pipe. Our operator splitting numerical scheme is unconditionally stable and it is second order accurate in space and time. The scheme is explicit, and it is formulated to work with general networks with loops. We test the scheme over range of regimes and network configurations, also comparing its performance with performance of two other state of the art implicit schemes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Body-centered orthorhombic C 16 : A novel topological node-line semimetal

We identify by ab initio calculations a novel topological semimetal carbon phase in all-sp 2 bonding networks with a 16-atom body-centered orthorhombic unit cell, termed bco-C 16 . Total-energy calculations show that bco-C 16 is comparable to solid fcc-C 60 in energetic stability, and phonon and molecular dynamics simulations confirm its dynamical stability. This all-sp 2 carbon allotrope can be regarded as a three-dimensional modification of graphite, and its simulated x-ray diffraction (XRD) pattern matches well a previously unexplained diffraction peak in measured XRD spectra of detonation and chimney soot, indicating its presence in the specimen. Electronic band structure calculations reveal that bco-C 16 is a topological node-line semimetal with a single nodal ring. Lastly, these findings establish a novel carbon phase with intriguing structural and electronic properties of fundamental significance and practical interest.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Impact of pH on Aqueous-Phase Phenol Hydrogenation Catalyzed by Carbon-Supported Pt and Rh

In aqueous phase, the rates of hydrogenation of aromatic substrates such as phenol on Pt/C and Rh/C are influenced by varying activity of hydronium ions. Decreasing the pH from 8 to 1 increases the rate of hydrogenation of phenol on Pt at 20 bar H 2 and 80 °C by 15-fold. This increase is attributed to weakening of the hydrogen binding energy (HBE) on the metal surface with decreasing pH. A weaker HBE at lower pH is also predicted by ab initio molecular dynamics simulations, providing atomistic insight into the impact of electrolyte ion distribution and interfacial solvent reorganization on HBE. The lower HBE results in a decrease in the activation energy for addition of adsorbed H from the metal to the adsorbed organic (with a Brønsted–Evans–Polanyi slope of ~1). In conclusion, the kinetic model derived accounts also for the lack of pH dependence at low hydrogen coverages (at 1 bar H 2 on Pt or up to 70 bar H 2 on Rh), when the weakening of the HBE decreases the hydrogen coverage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

The Potosi Reservoir Model 2013c, Property Modeling Update

As part of a larger project co-funded by the United States Department of Energy (US DOE) to evaluate the potential of formations within the Cambro-Ordovician strata above the Mt. Simon as potential targets for carbon sequestration in the Illinois and Michigan Basins, the Illinois Clean Coal Institute (ICCI) requested Schlumberger to evaluate the potential injectivity and carbon dioxide (CO2) plume size of the Cambrian Potosi Formation. The evaluation of this formation was accomplished using wireline data, core data, pressure data, and seismic data from this project as well as two other separately funded projects: the US DOE-funded Illinois Basin–Decatur Project (IBDP) being conducted by the Midwest Geological Sequestration Consortium (MGSC) in Macon County, Illinois, and the Illinois Industrial Carbon Capture and Sequestration (ICCS) project funded through the American Recovery and Reinvestment Act. In 2010, technical performance evaluations on the Cambrian Potosi Formation were performed through reservoir modeling. The data included formation tops from mud logs, well logs from the Verification Well #1 (VW1) and the Injection Well (CCS1), structural and stratigraphic formation from three dimensional (3D) seismic data, and field data from several waste water injection wells for Potosi Formation. The intention was for 2.2 million tons per annum (2 million tonnes per annum [MTPA]) of CO2 to be injected for 20 years. In the Task Error! Reference source not found., the 2010 Potosi heterogeneous model (referred to as the "Potosi Dynamic Model 2010") was re-run using a new injection scenario of 3.5 million tons per annum (3.2 MTPA) for 30 years. The extent of the Potosi Dynamic Model 2010, however, appeared too small for the new injection target. The models size was insufficient to accommodate the evolution of the plume. The new model, Potosi Dynamic Model 2013a, was built by extending the Potosi Dynamic Model 2010 grid to 30 by 30 mi (48 by 48 km), while preserving all property modeling workflows and layering. This model was retained as the base case. In the preceding Task [1], the Potosi reservoir model was updated to take into account the new data from the Verification Well #2 (VW2) which was drilled in 2012. The porosity and permeability modeling was revised to take into account the log data from the new well. Revisions of the 2010 modeling assumptions were also done on relative permeability, capillary pressures, formation water salinity, and the maximum allowable well bottomhole pressure. Dynamic simulations were run using the injection target of 3.5 million tons per annum (3.2 MTPA) for 30 years. This dynamic model was named Potosi Dynamic Model 2013b. In this Task, a new property modeling workflow was applied, where seismic inversion data guided the porosity mapping and geobody extraction. The static reservoir model was fully guided by PorosityCube interpretations and derivations coupled with petrophysical logs from three wells. The two main assumptions are: porosity features in the PorosityCube that correlate with lost circulation zones represent vugular zones, and that these vugular zones are laterally continuous. Extrapolation was done carefully to populate the vugular facies and their corresponding properties outside the seismic footprint up to the boundary of the 30 by 30 mi (48 by 48 km) model. Dynamic simulations were also run using the injection target of 3.5 million tons per annum (3.2 MTPA) for 30 years. This new dynamic model was named Potosi Dynamic Model 2013c. Reservoir simulation with the latest model gives a cumulative injection of 43 million tons (39 MT) in 30 years with a single well, which corresponds to 40% of the injection target. The injection rate is approx. 3.2 MTPA in the first six months as the well is injecting into the surrounding vugs, and declines rapidly to 1.8 million tons per annum (1.6 MTPA) in year 3 once the surrounding vugs are full and the CO2 start to reach the matrix. After, the injection rate declines gradually to 1.2 million tons per annum (1.1 MTPA) in year 18 and stays constant. This implies that a minimum of three (3) wells could be required in the Potosi to reach the injection target. The injectivity evaluated in this Task was higher compared to the preceding Task, since the current facies modeling (guided by the porosity map from the seismic inversion) indicated a higher density of vugs within the vugular zones. 5 As the CO2 follows the paths where vugs interconnection exists, a reasonably large and irregular plume extent was created. After 30 years of injection, the plume extends 13.7 mi (22 km) in E-W and 9.7 mi (16 km) in N-S directions. After injection finishes, the plume continues to migrate laterally, mainly driven by the remaining pressure gradient. After 60 years post-injection, the plume extends 14.2 mi (22.8 km) in E-W and 10 mi (16 km) in N-S directions, and remains constant as the remaining pressure gradient has become very low. Should the targeted cumulative injection of 106 million tons (96 MT) be achieved; a much larger plume extent could be expected. The increase of reservoir pressure at the end of injection is approximately 1,200 psia (8,274 kPa) around the injector and gradually decreases away from the well. The reservoir pressure increase is less than 10 psia (69 kPa) beyond 14 mi (23 km) away from injector. Should the targeted cumulative injection of 106 million tons (96 MT) be achieved; a much larger areal pressure increase could be expected. The reservoir pressure declines rapidly during the first 30 years post injection and the initial reservoir pressure is nearly restored after 100 years post-injection. The present evaluation is mainly associated with uncertainty on the vugs permeability and interconnectivity. The use of porosity mapping from seismic inversion might have reduced the uncertainty on the lateral vugs body distributions. However, major uncertainties on the Potosi vugs permeability remains. Therefore, injection test and pressure interference test among the wells could be considered to evaluate the local vugs permeability, extent, and interconnectivity. Facies modeling within the Potosi has yet to be thoroughly addressed. The carbonates during the time of deposition are believed to be regionally extensive. However, it may be worth delineating the reservoir with other regional wells or modern day analogues to understand the extent of the Potosi. More specifically, the model could incorporate lateral changes or trends if deemed necessary to represent facies transition. Data acquisitions to characterize the fracture pressure gradient, the formation water properties, the relative permeability, and the capillary pressure could also be considered in order to allow a more rigorous evaluation of the Potosi storage performance. A simulation using several injectors could also be considered to determine the required number of wells and appropriate spacing to achieve the injection target while taking into account the pressure interference.

Adushita, Yasmin↗
DOE OSTI2023

ReaxFF reactive force field model enables accurate prediction of physiochemical and mechanical properties of crystalline and amorphous shape‐memory polyurethane

Abstract We investigate reactive force field (ReaxFF) prediction of physiochemical and mechanical properties, and stimuli response behavior of a shape‐memory polyurethane (SMPU) rigid segment. We used SMPUs as a platform to link the rigid/hard segment domain crystallinity with properties that can also be experimentally measured. Specifically, we developed models for crystalline and amorphous forms of a common rigid segment of SMPUs, for example, 4,4′‐diphenylmethane diisocyanate (MDI) with n‐butanediol (BDO), denoted as MDI‐BDO. ReaxFF molecular dynamics simulations with equilibrium dynamics are performed by controlling mass, temperature, and pressure or volume to predict the structural and physiochemical properties of crystalline and amorphous rigid segments of SMPU systems, followed by simulating their XRD and FTIR patterns, respectively. Non‐equilibrium ReaxFF simulations are implemented by using uniaxial box deformation technique to study response‐behavior of crystalline SMPU under tensile loading and during stress‐relaxation. The model and simulation scheme presented in this study demonstrates a robust theoretical pathway toward application‐dictated, on‐demand design of high‐performance SMPUs in applications such as soft robotics.

Afroz, Mohammad M.↗
DOE OSTI2024

Simulating water dynamics related to pedogenesis across space and time: Implications for four-dimensional digital soil mapping

Digital soil mapping (DSM) relies on machine-learning and geostatistics to represent soil property observations across space. DSM techniques are powerful but often empirical, being limited to the quality and density of point samples. Water dynamics are closely related to soil variability, and the physics that govern water movement are well known. Hydrological properties can hence be simulated by physical models through space and time, unveiling key characteristics about soils. We propose the use of hydrologic models to map soils across the surface (2D), depth (1D), and time (1D)–which provides a 4D approach to digital soil mapping (4DSM). The Distributed Hydrology Soil Vegetation Model (DHSVM) was applied to a watershed currently under pasture. Moisture sensors and wells were installed at different depths in the watershed on summit, sideslope and toeslope positions to validate the model. DHSVM simulations of soil moisture distribution and depth to saturation were performed during the hydrological year (October 2008-September 2009). Clusters of similar pixels based on soil moisture values were determined using Dynamic Time Warping (DTW) to align temporal data and K-means. Clustering was performed both seasonally and for the entire year. Temporal patterns simulated by DHSVM matched measurements given by moisture sensors and wells. Seasonal clusters differed from the annual cluster. Distinct clusters were observed for each season and with depth, showing that spatiotemporal soil variability is lost when statically assessing soils. Spatiotemporal clusters corroborated field observations of fragipan occurrence not explicitly spatially mapped by Soil Survey Geographic Database (SSURGO). If a connection can be made between water and soils, static and dynamic soil variability can be predicted using physically based hydrologic models. Hydrologic models can benefit soil mapping by enabling reliable 4D simulation of water dynamics, which are fundamental to soil variability and soil classification and directly relate to biological, physical and chemical soil processes not captured by typical soil sampling protocols.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Cluster Dynamics Modeling with Bubble Nucleation, Growth and Coalescence

The topic of this communication pertains to defect formation in irradiated solids such as plasma-facing tungsten submitted to helium implantation in fusion reactor com- ponents, and nuclear fuel (metal and oxides) submitted to volatile ssion product generation in nuclear reactors. The purpose of this progress report is to describe ef- forts towards addressing the prediction of long-time evolution of defects via continuum cluster dynamics simulation. The di culties are twofold. First, realistic, long-time dynamics in reactor conditions leads to a non-dilute di usion regime which is not accommodated by the prevailing dilute, stressless cluster dynamics theory. Second, long-time dynamics calls for a large set of species (ideally an in nite set) to capture all possible emerging defects, and this represents a computational bottleneck. Extensions beyond the dilute limit is a signi cant undertaking since no model has been advanced to extend cluster dynamics to non-dilute, deformable conditions. Here our proposed approach to model the non-dilute limit is to monitor the appearance of a spatially localized void volume fraction in the solid matrix with a bell shape pro le and insert an explicit geometrical bubble onto the support of the bell function. The newly cre- ated internal moving boundary provides the means to account for the interfacial ux of mobile species into the bubble, and the growth of bubbles allows for coalescence phenomena which captures highly non-dilute interactions. We present a preliminary interfacial kinematic model with associated interfacial di usion transport to follow the evolution of the bubble in any number of spatial dimensions and any number of bubbles, which can be further extended to include a deformation theory. Finally we comment on a computational front-tracking method to be used in conjunction with conventional cluster dynamics simulations in the non-dilute model proposed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

RF Gap Transformation in PARMILA

The PARMILA (Phase And Radial Motion in Ion Linear Accelerators) code is used to design ion linacs and simulate the beam dynamics therein. The original code was developed in the 1960’s and later the Los Alamos Accelerator Code Group developed a PC version that has been used worldwide. Recently, in an effort to modernize the code and use it as the basis for a high-performance beam dynamics simulator with the LANSCE linac, the details of the code have come under additional scrutiny. PARMILA is a “z code” that uses various transformations to approximately represent the effect of the beam line element on the beam. In the case of the RF gap, it treats the particle dynamics inside the gap as a combination of drift + kick at the mid-plane + drift. Recently, we uncovered an undocumented keyword “off-exact” in PARMILA that enables the code to use the full modified-Bessel functions instead of low-order approximations in the RF gap transformation. This change improved the accuracy of the transformation for off-axis particles. However, a discrepancy between the analytical formula and the code was found for the transformation of the transverse coordinates. The purpose of this tech note is to document the derivation of the gap transformation and the related DTL cell length calculation used in the PARMILA code.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Population dynamics, information transfer, and spatial organization in a chemical reaction network under spatial confinement and crowding conditions

Here, we investigate, via Brownian dynamics simulations, the reaction dynamics of a generic, nonlinear chemical network under spatial confinement and crowding conditions. In detail, the Willamowski-Rossler chemical reaction system has been “extended” and considered as a prototype reaction-diffusion system. These results are potentially relevant to a number of open problems in biophysics and biochemistry, such as the synthesis of primitive cellular units (protocells) and the definition of their role in the chemical origin of life and the characterization of vesicle-mediated drug delivery processes. More generally, the computational approach presented in this work makes the case for the use of spatial stochastic simulation methods for the study of biochemical networks in vivo where the “well-mixed” approximation is invalid and both thermal and intrinsic fluctuations linked to the possible presence of molecular species in low number copies cannot be averaged out.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Anomalous fast atomic dynamics in bulk metallic glasses

It is believed that the strong glass formers tend to be dense packing and have rather slow atomic dynamics. Here, we report that the anomalous fast atomic dynamics exist in an Au 55 Cu 20 Ag 5 Si 20 metallic glass (MG) with good glass forming ability (GFA) by combining X-ray photon correlation spectroscopy measurements with ab initio molecular dynamics simulations. The fast dynamics are mainly originated from the increased mobility of Cu atoms by Ag addition. Ag atoms are inclined to locate at Cu-centered <0, 2, 8, 1> motifs, which slow or accelerate the dynamics of shell Cu atoms, depending on whether Ag atoms appear at the first shells or not, leading to the formation of structural heterogeneity at the atomic level. The results enrich the understanding of atomic dynamics in Cue-Ag-containing MGs and will trigger more studies on the dynamics-GFA relationship in MGs.

74 ATOMIC AND MOLECULAR PHYSICS↗