Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamics 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 541 records · Page 30

Computational Fluid Dynamics Simulation of Methane Slosh and Drain Experiments

NASA possesses a wealth of historical cryogenic experiments that provide valuable insights into the design, troubleshooting, and understanding involved in the complex fluid and thermodynamics of managing cryogenic propellants. One approach to leveraging these historical datasets is by simulating these experiments to validate the accuracy of simulation environments and constituent models. This study focused on simulating a selection of the K site test series for both static and sloshing pressurized liquid methane draining experiments conducted at NASA in the 1970’s, utilizing computational fluid dynamics. The simulations were performed in the ANSYS FLUENT environment using the Volume of Fluid (VOF) numerical approach. A k-omega turbulence model was used with interfacial turbulence damping, and accurately predicted the amount of pressurant needed to maintain the required tank pressure throughout the static drain. The static simulation predicted the temperature stratification in the ullage observed at the end of the drain. During the methane expulsion with sloshing test, many features were successfully captured using the k-omega turbulence model with interfacial turbulence damping included. The rate of phase change and liquid temperature was overpredicted compared to the experimental measurements. The overprediction may be attributed to uncertainties in the vessel geometry, methane pressurant temperature and composition, and methodological differences in how the sloshing frequency was adjusted during the expulsion.

Cryogenic Propellants↗

Probing High-Pressure Structural Evolution in Polyurea with In Situ Energy-Dispersive X-ray Diffraction and Molecular Dynamics Simulations

Polyurea, an elastomer with a phase-segregated microstructure, has been proven as an effective coating in defense applications. To gain a more complete understanding of the high-pressure atomic-level morphology of these phases and to validate molecular dynamics (MD) simulations, multi-angle energy-dispersive X-ray diffraction experiments were performed in situ up to pressures of ~6 GPa at room temperature. Structure factors were obtained and compared to MD simulations with an average error of less than 5% between major peak positions. The first sharp diffraction peak shifted from 4.56 Å to lower d-spacing with pressure, indicating compression between hard segments. This was further supported by the behavior of a peak at ~3.86 Å from the pair distribution function (PDF), suspected to represent π-stacking and separation between soft segments. Compression within the hard segments themselves is minimal as low-r peaks in the PDF are not greatly affected by pressure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mott Transition in a Metallic Liquid: Gutzwiller Molecular Dynamics Simulations

We present here a formulation of quantum molecular dynamics that includes electron correlation effects via the Gutzwiller method. Our new scheme enables the study of the dynamical behavior of atoms and molecules with strong electron interactions. The Gutzwiller approach goes beyond the conventional mean-field treatment of the intra-atomic electron repulsion and captures crucial correlation effects such as band narrowing and electron localization. Finally, we use Gutzwiller quantum molecular dynamics to investigate the Mott transition in the liquid phase of a single-band metal and uncover intriguing structural and transport properties of the atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗

Integrating CryoEM and Molecular dynamic simulations [Slide]

López et al. analyze the structure and dynamics of lipid nanodiscs, both empty and in complex with pro-apoptotic BAX protein. Information from the computational derived structures prompted changes in image processing strategy, which led to quality improvements in Cryo-EM reconstructions. Together, these findings provide novel insights on the activity of BAXv.

59 BASIC BIOLOGICAL SCIENCES↗

From Electronic Structure to Ion Transport: Photoelectron Spectroscopy and Molecular Dynamics Simulations Reveal the Role of Anions in Lithium Battery Electrolytes

Electrolyte anions are pivotal for lithium battery performance, yet their fundamental electronic structural properties are not well understood. In this work, we employ a combination of negative-ion photoelectron spectroscopy (NIPES), ab initio calculations, and molecular dynamics (MD) simulations to investigate the electronic structures of three representative electrolyte anions. This multiscale approach enables us to elucidate how their intrinsic electronic properties govern anion–solvent interactions in gas-phase clusters, as well as lithium-ion (Li + ) solvation structures and ion transport behavior in the condensed phase. NIPES reveals that difluoro(oxalato)borate (DFOB – ), bis(fluorosulfonyl)imide (FSI – ), and bis(oxalato)borate (BOB – ) all exhibit high electron binding energies, with vertical/adiabatic detachment energies increasing from DFOB – (6.09/5.70 eV) to FSI – (6.80/6.10 eV) to BOB – (6.82/6.40 eV), correlating with enhanced oxidation stability. Ab initio calculations reveal that DFOB – /FSI – –solvent complexes bind Li + ∼ 10 kcal/mol stronger than BOB – series, aligning with the strength of a Li + –anion model. DFOB – exhibits pronounced charge localization on both oxygen and fluorine atoms, enabling their involvement in Li + coordination. In contrast, fluorine atoms in FSI – are largely electron-depleted and remain excluded from direct Li + binding. MD simulations further demonstrate that LiDFOB and LiFSI systems exhibit Li + diffusion coefficients three and five times higher than those of LiBOB across four common solvents. Notably, LiFSI salt in acetonitrile (AN) exhibits the fastest Li + diffusion among 12 electrolyte systems, highlighting the synergistic effect of FSI – and AN in promoting ion mobility. In conclusion, these findings provide a molecular-level understanding of the critical roles of anion and its microsolvation in optimizing Li + diffusion dynamics, once again emphasizing the positioning of FSI – and DFOB – as prime candidates for next-generation electrolytes.

25 ENERGY STORAGE↗

Molecular Dynamics Simulations of Austenite-Martensite Interfaces in NiTi Shape Memory Alloys

The unique properties of shape memory alloys (SMAs) arise from a reversible martensitic transformation. The nucleation and migration of austenite-martensite interfaces are the key to understanding the SMA properties. Molecular dynamics (MD) simulations can provide important atomic-scale information about these aspects of the transformation, but their time scales prevent the interface formation under near-equilibrium conditions relevant to experiment. We present a new MD methodology which allows for the natural formation of energetically preferred austenite-martensite interfaces under near-equilibrium conditions. Our simulation demonstrates that the interfaces in NiTi are semi-coherent, composed of a series of terrace planes and structural disconnections, and they migrate rapidly through single crystals with only a small thermodynamic driving force. In bi-crystals and polycrystals, the migration of these same interfaces is significantly impeded by grain boundaries and stored elastic energy. This behavior can result in SMA hysteresis via several mechanisms associated with nucleation and non-elastic strain accommodation.

Gabriel Plummer↗

Dynamic simulation of powder packing structure for powder bed additive manufacturing

Here, powder packing structure is a critical parameter of powder bed-based additive manufacturing (AM). Experimental characterization of powder is typically limited to measuring bulk properties, whereas many numerical models of AM powder packing are based on geometrical consideration without accounting for particle-to-particle interactions. In the present paper, the powder packing dynamics is simulated using a discrete element method (DEM)-based model that solves the mechanical contact forces and moments between individual particles. As DEM uses explicit time integration, a main challenge in modeling dynamics of metallic powder packing is the need for extremely fine time increment size (e.g., in the order of 1 ns for a 10-μm-diameter particle). The effect of mass scaling, employed for speeding up the calculation, on the simulation results is examined in a test case of powder particles packed inside a box container. The calculated packing density for two different particle size distributions is validated against independent literature data for laser powder bed AM with AISI 316L stainless steel powder. The sensitivity of key input parameters (e.g., friction coefficient) is further evaluated in this test case. The powder packing model is then applied to a practical situation of binder jet AM involving rolling of multiple layers of IN718 powder particles onto a powder bed, for which the calculated packing density is also validated with independent literature data.

42 ENGINEERING↗

Polariton relaxation under vibrational strong coupling: Comparing cavity molecular dynamics simulations against Fermi’s golden rule rate

Under vibrational strong coupling (VSC), the formation of molecular polaritons may significantly modify the photo-induced or thermal properties of molecules. In an effort to understand these intriguing modifications, both experimental and theoretical studies have focused on the ultrafast dynamics of vibrational polaritons. Here, following our recent work [Li et al., J. Chem. Phys. 154, 094124 (2021)], we systematically study the mechanism of polariton relaxation for liquid CO 2 under a weak external pumping. Classical cavity molecular dynamics (CavMD) simulations confirm that polariton relaxation results from the combined effects of (i) cavity loss through the photonic component and (ii) dephasing of the bright-mode component to vibrational dark modes as mediated by intermolecular interactions. The latter polaritonic dephasing rate is proportional to the product of the weight of the bright mode in the polariton wave function and the spectral overlap between the polariton and dark modes. Both these factors are sensitive to parameters such as the Rabi splitting and cavity mode detuning. Compared to a Fermi’s golden rule calculation based on a tight-binding harmonic model, CavMD yields a similar parameter dependence for the upper polariton relaxation lifetime but sometimes a modest disagreement for the lower polariton. We suggest that this disagreement results from polariton-enhanced molecular nonlinear absorption due to molecular anharmonicity, which is not included in our analytical model. Finally, we also summarize recent progress on probing nonreactive VSC dynamics with CavMD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Iron surface corrosion in supercritical CO2 at atomic scale investigated by molecular dynamics simulations

Understanding the corrosion behavior of steels in supercritical carbon dioxide (S-CO2) is essential for ensuring the safe application of S-CO2 as a heat-transfer fluid in high-temperature energy systems, including advanced nuclear reactors. In this work, molecular dynamics (MD) simulations using ReaxFF potential are performed to explore the atomic-scale corrosion mechanisms of body-centered cubic iron (BCC-Fe) in S-CO2. The results show that CO2 molecules in S-CO2 decompose at the Fe surface, generating free C and O atoms that form Fe-C and Fe-O bonds and subsequently produce oxides and carbides. Concurrently, Fe atoms dissolve from the surface and diffuse into the S-CO2 region, resulting in interdiffusion of Fe, C and O atoms at the interface. The corrosion-layer thickness calculations show that high pressure and temperature induced by S-CO2 have stronger effects than surface orientation on the corrosion process. In addition, surface Fe atoms undergo substantial displacement under S-CO2 exposure, further accelerating corrosion. When a radiation-induced void is introduced near the Fe surface, the corrosion is enhanced. The void-matrix interface expands the reaction surface area and simultaneously induces corrosion reactions inside the bulk, resulting in a deeper penetration of C and O and thicker corrosion layers. All these results indicate that high-temperature, high-pressure and radiation-induced voids can seriously affect the corrosion of Fe in S-CO2, and must be considered to better use S-CO2 in nuclear facilities.

Li, Wenhua↗

EXAFS and molecular dynamics simulation studies of Cu-Zr metallic glass: Short-to-medium range order and glass forming ability

Highlights: • Atomic structures of Cu{sub 100−x}Zr{sub x} MGs were studied by EXAFS and MD simulations. • The change of GFA with composition is not dominated by the SRO and MROs. • An apparent reduction of “soft region” at the local high-GFA compositions. - Abstract: Atomic structures of Cu{sub 100−x}Zr{sub x} (x = 41.18–66.67) metallic glasses were investigated by extended X-ray absorption fine structure (EXAFS) spectroscopy and molecular dynamics (MD) simulations. It is demonstrated that both the fraction and average volume of the typical atomic clusters, characterized by Voronoi cells, are of compositional dependence. With the increasing of Zr content, the fraction of Cu and Zr centered Voronoi clusters with coordination number lower than 11 and 14 respectively increases, while that of clusters with larger coordination numbers reduces. Among the major Voronoi clusters, the Cu-centered icosahedra are distinctive. Their average volume is far smaller than that of the other Cu-centered Voronoi cells with the same coordination number. More importantly, they prefer to interlink into pentagon-rich icosahedral super-clusters. However, compared with their neighbors, the glasses with x = 44 and x = 50 exhibit no other particularity in Cu-centered icosahedra than a slightly enhanced interlink at x = 50 although glass forming ability (GFA) reaches a local maximum at these two compositions, indicating that there are other factors dominating the GFA change with composition. A significant reduction in free volume for the loosely packed regions takes place at x = 44 and x = 50, due to which the corresponding Cu-Zr metallic glasses have the maximum microhardness.

36 MATERIALS SCIENCE↗

Predicting initial dissolution rates using structural features from molecular dynamics simulations

Predicting chemical durability of glass materials is important for various applications from daily life such as drink glass and kitchen ware to advanced technologies such as nuclear waste disposal and biomedicine. In this work, we explored prediction of initial dissolution rate through structural features from molecular dynamics (MD) simulations for a wide range of glass compositions (total 28) including borosilicate and aluminosilicate glasses, ZrO 2 -containing and V 2 O 5 -containing boroaluminosilicate glasses. The initial dissolution rates (r 0 ) measured experimentally at 90 °C with varying solution conditions were correlated with structural features (e.g., polyhedral linkages and non-bridging oxygen species) obtained from MD simulations, either from this study or from literature. Since hydrolysis of the glass network through breaking of the network former linkages (e.g., Si-O-Si, Si-O-Al, etc.) is a critical step of network glass dissolution, the statistics of these linkages obtained from MD were correlated to r0 through linear regression, where the coefficient of determination (R 2 ) and root mean square error are found to be 0.949 and 0.681, respectively. This model was compared and discussed with existing models developed by various approaches including machine learning, the kinetic rate equation, topological constraint theory, and other descriptors from MD simulations. The discussion provides insights on future model improvements to predict glass dissolution. In addition, the impact of V 2 O 5 on the glass dissolution was examined in detail, implicating that the impact is not the same across all glass compositions and test conditions.

36 MATERIALS SCIENCE↗

Ab initio study of the structure and properties of amorphous silicon hydride from accelerated molecular dynamics simulations

This paper presents a large-scale ab initio simulation study of amorphous silicon hydride (a-Si 1-x H x ) with an emphasis on the structure and properties of the material across a range of hydrogen concentration by combining accelerated molecular dynamics (MD) simulations with first-principles density-functional calculations. The accelerated MD scheme relied on classical metadynamics, which enabled the development of 2500+ high-quality structural models of a-Si 1-x H x , with system sizes ranging from 150 to 6000 atoms and hydrogen concentrations vary from 6 to 20 at. %. The resulting amorphous networks were found to be completely free from any coordination defects and that they all exhibited a pristine band-gap in their electronic spectrum. The microstructural properties of hydrogen distributions were examined with an emphasis on the presence of isolated and clustered environments of hydrogen atoms. The results were compared with experimental data obtained from X-ray diffraction, infrared spectroscopy and nuclear magnetic resonance studies.

36 MATERIALS SCIENCE↗

Ab initio study of the structure and properties of amorphous silicon hydride from accelerated molecular dynamics simulations

This paper presents a large-scale ab initio simulation study of amorphous silicon hydride (a-Si 1-x H x ) with an emphasis on the structure and properties of the material across a range of hydrogen concentration by combining accelerated molecular dynamics (MD) simulations with first-principles density-functional calculations. The accelerated MD scheme relied on classical metadynamics, which enabled the development of 2600+ high-quality structural models of a-Si 1-x H x , with system sizes ranging from 150 to 6,000 atoms and hydrogen concentrations vary from 6 to 20 at. %. The resulting amorphous networks were found to be completely free from any coordination defects and that they all exhibited a pristine band-gap in their electronic spectrum. The microstructural properties of hydrogen distributions were examined with great emphasis on the presence of isolated and clustered environments of hydrogen atoms. The results were compared with a suite of experimental data obtained from x-ray diffraction, infrared spectroscopy, spectroscopic ellipsometry and nuclear magnetic resonance studies.

36 MATERIALS SCIENCE↗

Stability of Calcium Ion Battery Electrolytes: Predictions from Ab Initio Molecular Dynamics Simulations

Multivalent batteries, such as magnesium-ion, calcium-ion, and zinc-ion batteries, have attracted significant attention as next-generation electrochemical energy storage devices to complement conventional lithium-ion batteries (LIBs). Among them, calcium-ion batteries (CIBs) are the least explored due to the difficult reversible Ca deposition-dissolution. In this work, we examined the stability of four different Ca salts with weakly coordinating anions and three different solvents commonly employed in existing battery technologies to identify suitable candidates for CIBs. By employing Born-Oppenheimer molecular dynamics (BOMD) simulations on salt-Ca and solvent-Ca interfaces, we find that the tetraglyme solvent and carborane salt are promising candidates for CIBs. Due to the strong reducing nature of the calcium surface, the other salts and solvents readily decompose. Further, we explain the microscopic mechanisms of salt/solvent decomposition on the Ca surface using time-dependent projected density of states, time-dependent charge-transfer plots, and climbing-image nudged elastic band calculations. Collectively, this work presents the first mechanistic assessment of the dynamical stability of candidate salts and solvents on a Ca surface using BOMD simulations, and provides a predictive path toward designing stable electrolytes for CIBs.

Born-Oppenheimer molecular dynamics↗

Nonadiabatic molecular dynamics simulation of C2H22+ in a strong laser field*

We investigate the alignment dependence of the strong laser dissociation dynamics of molecule C 2 H 2 2 + in the frame of real-time and real-space time-dependent density function theory coupled with nonadiabatic quantum molecular dynamics (TDDFT-MD) simulation. This work is based on a recent experiment study “ultrafast electron diffraction imaging of bond breaking in di-ionized acetylene” [Wolter et al , Science 354 , 308–312 (2016)]. Our simulations are in excellent agreement with the experimental data and the analysis confirms that the alignment dependence of the proton dissociation dynamics comes from the electron response of the driving laser pulse. Our results validate the ability of the TDDFT-MD method to reveal the underlying mechanism of experimentally observed and control molecular dissociation dynamics.

Physics↗

Shock compression of crystalline TeO 2 to the high-pressure fluid regime: Insights from ab initio molecular dynamics simulations

The shock response of fully-dense and porous crystalline tellurium dioxide (TeO 2 ⁠) to the high-pressure and high-temperature fluid regime was investigated within the framework of density functional theory with Mermin’s generalization to finite temperatures. The principal and porous shock Hugoniot curves were predicted from canonical ab initio molecular dynamics (AIMD) simulations, with the phase space sampled along isotherms up to 80 000 K, for densities ranging from ρ = 3 to 17 g/cm 3 . The polymorphs investigated are α-TeO 2 paratellurite (⁠P4 1 2 1 2), TeO 2 cotunnite (⁠Pnma⁠), and TeO 2 post-cotunnite (⁠P2 1 /m⁠). Based on the discontinuity found in the calculated U s – u p slope of TeO 2 post-cotunnite at a shock velocity of U s ≃ 8.35km/s and a particle velocity of u p ≃ 3.64 km/s, the shock melting temperature and pressure are predicted to be ≃ 6500 K and ≃ 170 GPa. Results from the AIMD simulations are in line with the static compression data of TeO 2 paratellurite and cotunnite, and with the recent shock Hugoniot data for single-crystal α- TeO 2 for pressures up to 85 GPa, obtained using the inclined-mirror method and the velocity interferometer system for any reflector combined with powder gun and two-stage light-gas gun.

74 ATOMIC AND MOLECULAR PHYSICS↗

Adsorption mechanisms of cesium at calcium-silicate-hydrate surfaces using molecular dynamics simulations

We report understanding the mechanisms of interactions and immobilization of radionuclides in cement systems at the molecular level is key to developing safe and novel solutions for the long-term storage of nuclear wastes. The molecular level adsorption mechanisms of cesium (Cs + ) ions onto tobermorite 9 Å, tobermorite 14 Å, and jennite – three crystalline structural analogues for calcium-silicate-hydrates, the main components in cement-based materials, were investigated using molecular dynamics (MD) simulations. Convergence monitoring of the simulations, using the root mean square displacement with average correlation analysis, indicated that MD trajectories of 10–15 ns were needed given the size of the selected computational cell to properly capture the dynamic process of adsorption. Cs + ions adsorbed as inner-sphere complexes at the tobermorite surfaces, while they showed lower affinity for jennite. The charge and structure of the surface layer strongly influenced the adsorption mechanisms of Cs + ions. Strong association with the hexagonal cavities of the silica sites were seen at the tetrahedral SiO 4 surface of the tobermorite structures. Cs + ion adsorption was attributed to co-ion adsorption at the octahedral CaO 6 surface of tobermorite 9 Å, cation exchange at that of tobermorite 14 Å, and electrostatic interaction with the surface hydroxyl groups for jennite. Distribution coefficients of Cs + ions onto cement pastes derived from molecular dynamics by taking the surface structure and dynamics of Cs + ions directly into account showed a range of overlap with experimental results.

42 ENGINEERING↗

Ballistic-diffusive phonon transport in cellulose nanocrystals by ReaxFF molecular dynamics simulations

Cellulose nanocrystals (CNCs), a ubiquitous nano-sized natural biopolymer, have drawn considerable attention in recent years due to their excellent mechanical and chemical properties. However, thermal properties of CNCs, which are critical for their potential applications as structural or functional materials, are less investigated. In this work, their thermal properties were systematically studied by molecular dynamics (MD) simulations in terms of the effects of polymorphs, size, temperature, and strain. It was found that the thermal conductivities depend on the polymorphs of CNCs, but the thermal conductance of unit chain of each polymorph is similar. Moreover, the strong dependence of the thermal conductivities of CNCs on their length and cross-sectional area indicates the existence of remarkable ballistic-diffusive phonon transport. Single chain CNCs can achieve a thermal conductivity of ~6 W/(m∙K), implying their potential as high thermal conductive polymer. The mechanical strain can enhance the alignment of cellulose skeleton, which reduces phonon scattering and thus increases thermal conductivity. Possible underlying mechanism of the varying thermal conductivity was further discussed and attributed to the shift of phonon frequencies and varied degree of orientation. This study reveals the sophisticated ballistic-diffusive phonon transport in CNCs, and pave path to future design of CNC-based materials with desirable thermal properties.

42 ENGINEERING↗