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 307 records · Page 17

Electromechanical Dynamics Simulations of Superconducting LSM Rocket Launcher System in Attractive-Mode

Further feasibility study on a superconducting linear synchronous motor (LSM) rocket launcher system is presented on the basis of dynamic simulations of electric power, efficiency and power factor as well as the ascending motions of the launcher and rocket. The advantages of attractive-mode operation are found from comparison with repulsive-mode operation. It is made clear that the LSM rocket launcher system, of which the long-stator is divided optimally into 60 sections according to launcher speeds, can obtain high efficiency and power factor.

Yoshida, Kinjiro↗

Modeling bicarbonate formation in an alkaline solution with multi-level quantum mechanics/molecular dynamics simulations

Understanding carbonate speciation and how it may be modulated is essential for the advancement of carbon dioxide (CO 2 ) capture and storage technologies, which often rely on the transformation of CO 2 into carbonate, e.g. via the formation of carbonate minerals. To date, few atomic-level, quantum-mechanics-based simulations have been carried out to characterize how carbonic acid (H 2 CO 3 ) and bicarbonate ($HCO^{-}_{3}$) form in aqueous solution, and how pH affects this process. Recently, Martirez and Carter utilized rare-event sampling density functional theory molecular dynamics simulations in combination with multi-level embedded correlated wavefunction theory, thus accounting for both solvent dynamics and electron correlation accurately, to elucidate the mechanism of H 2 CO 3 formation in neutral solution (J. Am. Chem. Soc., 145, 12561, 2023). Here, we perform a complementary simulation using the same method to map out the energetics of $HCO^{-}_{3}$ formation from dissolved CO 2 in basic solution. We find that, as in H 2 CO 3 formation, including water dynamics is important to obtain an accurate prediction of the energetics for the aforementioned reaction. Furthermore, only with MD did we identify the correct pathway for the reaction, in which water – not hydroxide – acts as the initial nucleophile and only at the transition state does it lose a proton.

74 ATOMIC AND MOLECULAR PHYSICS↗

Molecular dynamics simulation of zirconium tungstate amorphization and the amorphous-crystalline interface

Classical molecular dynamics (MD) simulations were performed to provide a conceptual understanding of the amorphous-crystalline interface for a candidate negative thermal expansion (NTE) material, ZrW 2 O 8 . Simulations of pressure-induced amorphization at 300 K indicate that an amorphous phase forms at pressures of 10 GPa and greater, and this phase persists when the pressure is subsequently decreased to 1 bar. However, the crystalline phase is recovered when the slightly distorted 5 GPa phase is relaxed to 1 bar. Simulations were also performed on a two-phase model consisting of the high-pressure amorphous phase in direct contact with the crystalline phase. Upon equilibration at 300 K and 1 bar, the crystalline phase remains unchanged beyond a thin layer of disrupted structure at the crystalline-amorphous interface. In conclusion, the differences in local atomic structure at the interface are quantified from the simulation trajectories.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Molecular Dynamics Simulation of the Migration of Tilt Grain Boundaries in Ni and Ni{sub 3}Al

The migration of the 〈100〉 and 〈111〉 tilt boundaries in Ni and the Ni{sub 3}Al intermetallic compound is studied by molecular dynamics simulation. The low-angle 〈100〉 boundaries in Ni and Ni{sub 3}Al are shown to migrate much more slowly (by approximately two times at a temperature of 1700 K) that the 〈111〉 boundaries, which is related to different migration mechanisms of the 〈100〉 and 〈111〉 boundaries. The migration of the 〈100〉 boundaries is found to occur by the splitting of paired grain-boundary dislocations followed by a change of partner dislocations, which proceeds via the slip of split dislocations. The 〈111〉 tilt boundaries move according to a combined action of the following two mechanisms: the mechanism described above and the mechanism consisting in joint slip of paired grain-boundary dislocations, which have common slip planes in contrast to the grain-boundary dislocations in 〈100〉 boundaries. Since the joint slip of paired dislocations has a relatively low activation energy, the mobility of the 〈111〉 tilt boundaries is higher than that of the 〈100〉 boundaries. Our molecular dynamics studies demonstrate that the migration rate of analogous boundaries in Ni{sub 3}Al is significantly lower than in Ni (by approximately three times at a temperature of 1700 K), which is partly caused by the additional energy consumed for the formation of a disordered region behind a migrating boundary in Ni{sub 3}Al. Because of a relatively low boundary mobility in Ni{sub 3}Al, the contribution of diffusion atomic displacements during boundary migration is higher as compared to Ni.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Pressure-induced phase transformations in Nb 2 O 5 from ab initio molecular dynamics simulations

Here, the equations of state for the H-, L- and T-Nb 2 O 5 polymorphs were investigated with density functional theory. The shock Hugoniot curves were predicted from ab initio molecular dynamics simulations along 300–2500 K isotherms. Enthalpies derived from simulations at 300 K show that the H → L transition occurs at 6.2 GPa, in agreement with the measured value of ~5.9 GPa, and is more favorable than the H → T transition. Pressure-induced amorphization (PIA) of L-Nb 2 O 5 is in the range ~18-25 GPa, consistent with experiments reporting PIA above 21.4 GPa. The simulated PIA of T-Nb 2 O 5 between ~16 and 26 GPa reproduces PIA observed above 19.2 GPa.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Path integral Monte Carlo and density functional molecular dynamics simulations of warm dense MgSiO 3

In order to provide a comprehensive theoretical description of MgSiO 3 at extreme conditions, we combine results from path integral Monte Carlo and density functional molecular dynamics simulations and generate a consistent equation of state for this material. We consider a wide range of temperature and density conditions from 10 4 to 10 8 K and from 0.321 to 64.2 g/cm -3 (0.1- to 20-fold the ambient density). We study how the L and K shell electrons are ionized with increasing temperature and pressure. We derive the shock Hugoniot curve and compare with experimental results. Our Hugoniot curve is in good agreement with the experiments, and we predict a broad compression maximum that is dominated by the K shell ionization of all three nuclei while the peak compression ratio of 4.70 is obtained when the Si and Mg nuclei are ionized. To conclude, we analyze the heat capacity and structural properties of the liquid.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quantum dynamics simulation of the advection-diffusion equation

The advection-diffusion equation is simulated via several quantum algorithms. Three formulations are considered: (1) Trotterization, (2) variational quantum time evolution (VarQTE), and (3) adaptive variational quantum dynamics simulation (AVQDS). These schemes were originally developed for the Hamiltonian simulation of many-body quantum systems. The finite-difference discretized operator of the transport equation is formulated as a Hamiltonian and solved without the need for ancillary qubits. Computations are conducted on a quantum simulator (IBM Qiskit Aer) and a superconducting quantum hardware (IBM Fez). The former emulates the latter without the noise. The actual hardware implementation experiences significant noise. The results of the quantum simulator are compared with data from direct numerical simulation (DNS) with infidelities of the order 10 −5 . In the quantum simulator, Trotterization is observed to have the lowest infidelity and is suitable for fault-tolerant computation. The AVQDS algorithm requires the lowest gate count and circuit depth. The VarQTE algorithm is the next best in terms of gate counts, but the number of its optimization variables is directly proportional to the number of qubits. Due to current hardware limitations, Trotterization cannot be implemented, as it has an overwhelmingly large number of operations. Meanwhile, AVQDS and VarQTE can be executed at the hardware level. These algorithms present a new paradigm for computational transport phenomena on quantum computers.

Alipanah, Hirad [Univ. of Pittsburgh, PA (United S↗

Topological Effects in Isolated Poly[ n ]catenanes: Molecular Dynamics Simulations and Rouse Mode Analysis

Poly[n]catenanes are mechanically interlocked polymers consisting of interlocking ring molecules. We report that over the years, researchers have speculated that the permanent topological interactions within the poly[n]catenane backbone could lead to unique dynamical behaviors. To investigate these unusual polymers, molecular dynamics simulations of isolated poly[n]catenanes have been conducted, along with a Rouse mode analysis. Owing to the mechanical bonds within the molecule, the dynamics of poly[n]catenanes at short length scales are significantly slowed and the distribution of relaxation times is broadened; these same behaviors have been observed in melts of linear polymers and are associated with entanglement. Lastly, despite these entanglement-like effects, at large length scales poly[n]catenanes do not relax much slower than isolated linear polymers and are less strongly impacted by increased segmental stiffness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular Dynamics Simulation on the Effect of MgO/Al{sub 2}O{sub 3} Ratio on Structure and Properties of Blast Furnace Slag Under Different Basicity Conditions

The SiO{sub 2}-Al{sub 2}O{sub 3}-CaO-MgO is the basic structural system of blast furnace slag and the composition directly affects the performance of the slag. Molecular dynamics simulations were carried out to analyze the local structure, structural unit, bond angle, transport properties, and enthalpy of the slag with the increase of MgO/Al{sub 2}O{sub 3} mass ratio under different basicity conditions. It was found that the change of MgO/Al{sub 2}O{sub 3} ratio does not affect the short-range order structure, but it will reduce the overall stability of the network structure. Through the analysis of the structural unit of the slag, it was found that the polymerization degree of the system decreases with the increase of MgO/Al{sub 2}O{sub 3} ratio, indicating that the complex network structure of the system is partially depolymerized. Besides, the self-diffusion coefficients of each ion were concluded and the magnitudes were observed to be in the following order: Mg{sup 2+}>Ca{sup 2+}>Al{sup 3+}>O{sup 2-}>Si{sup 4+}. The diffusivity of the slag increases with the increase of MgO/Al{sub 2}O{sub 3} ratio and the corresponding viscosity decreases. And the enthalpy of the slag is increased with the increase of MgO/Al{sub 2}O{sub 3} ratio or basicity, therefore the fuel consumption should also be properly adjusted during the operation of blast furnace to ensure the target temperature.

36 MATERIALS SCIENCE↗

Initial position optimization in molecular dynamics simulations for a Coulomb system

A new algorithm for molecular dynamics (MD) simulations is developed to optimize plasma particle distributions at given initial temperatures. By combining velocity scaling and reassignment, the method effectively eliminates the initial rise and oscillation in temperatures observed with randomly distributed positions. These rises and oscillations are undesired numerical artifacts observed in conventional plasma MD simulations, arising from unoptimized particle positions. The algorithm demonstrates temperature relaxation without initial rises or oscillations, as well as precise flow velocity relaxation, enabling accurate measurement of relaxation times. The code is accelerated using graphics processing units for parallel processing, enhancing the study of plasma dynamics. The proposed method for distributing physically valid particles in MD simulations enables accurate studies of intrinsic collision processes in plasmas, including the dynamics of strongly coupled plasmas, plasma–wave interactions, and transport phenomena in magnetized plasmas. The paper concludes with a discussion of potential applications and future enhancements to the algorithm.

Jo, Jawon (ORCID:0009000924193285)↗

Initial hydration process of calcium silicates in Portland cement: A comprehensive comparison from molecular dynamics simulations

As the main components of Portland cement, calcium silicates show substantial differences in their hydration reactivity which have not been fully explained. A comprehensive comparison of the initial hydration process of calcium silicates, namely dicalcium silicate (C{sub 2}S) and tricalcium silicate (C{sub 3}S), was conducted using molecular dynamics simulations. The initial hydration process was divided into three stages using cut-off times of 0.001 ns and 3 ns. The hydration of M{sub 3}-C{sub 3}S (010) was more evident than that of β-C{sub 2}S (100), supported by the hydroxylation degree, radial distribution function, atomic density profile, etc. The coordination number of the surface Ca atoms might be the underlying reason for such a difference. Interactions between cement surfaces and water molecules were mainly characterised by solid OH bonding and Ca-water O bonding. Dissolution of Ca atoms was observed, although quite scarce, while no dissolution of Si atoms was observed.

36 MATERIALS SCIENCE↗

Predicting CaO-(MgO)-Al2O3-SiO2 glass reactivity in alkaline environments from force field molecular dynamics simulations

In this investigation, force field-based molecular dynamics (MD) simulations have been employed to generate detailed structural representations for a range of amorphous quaternary CaO-MgO-Al{sub 2}O{sub 3}-SiO{sub 2} (CMAS) and ternary CaO-Al{sub 2}O{sub 3}-SiO{sub 2} (CAS) glasses. Comparison of the simulation results with select experimental X-ray and neutron total scattering and literature data reveals that the MD-generated structures have captured the key structural features of these CMAS and CAS glasses. Based on the MD-generated structural representations, we have developed two structural descriptors, specifically (i) average metal oxide dissociation energy (AMODE) and (ii) average self-diffusion coefficient (ASDC) of all the atoms at melting. Both structural descriptors are seen to more accurately predict the relative glass reactivity than the commonly used degree of depolymerization parameter, especially for the eight synthetic CAS glasses that span a wide compositional range. Hence these descriptors hold great promise for predicting CMAS and CAS glass reactivity in alkaline environments from compositional information.

36 MATERIALS SCIENCE↗

Determination of the structure of lithium niobosilicate glasses by molecular dynamics simulation with a new Nb-O potential

A partial charge empirical potential for the Nb-O pair has been developed based on existing parameters of other pairs in a Teter potential to enable molecular dynamics simulations of both lithium niobate crystal and lithium niobosilicate glasses. The developed potential is capable of describing structural features of lithium niobosilicate glasses in a wide composition range, including niobium coordination number and bond length, density, bond angle distribution, polyhedral distribution and fraction of bridging oxygen. Furthermore, the results obtained using this new potential show good agreement with experimental data of density and structure from Raman spectroscopy of lithium niobosilicate glasses and X-ray absorption spectroscopy near Nb K-edge.

36 MATERIALS SCIENCE↗

Determination of bimetallic architectures in nanometer-scale catalysts by combining molecular dynamics simulations with x-ray absorption spectroscopy

In this study, we present an approach for the determination of an atomic structure of small bimetallic nanoparticles by combining extended X-ray absorption fine structure spectroscopy and classical molecular dynamics simulations based on the Sutton-Chen potential. The proposed approach is illustrated in the example of PdAu nanoparticles with ca 100 atoms and narrow size and compositional distributions. Using a direct modeling approach and no adjustable parameters, we were able to reproduce the size and shape of nanoparticles as well as the intra-particle distributions of atoms and metal mixing ratios and to explore the influence of these parameters on the local structure and dynamics in nanoparticles.

74 ATOMIC AND MOLECULAR PHYSICS↗

Machine Learning‐Augmented Molecular Dynamics Simulations (MD) Reveal Insights Into the Disconnect Between Affinity and Activation of ZTP Riboswitch Ligands

Abstract The challenge of targeting RNA with small molecules necessitates a better understanding of RNA–ligand interaction mechanisms. However, the dynamic nature of nucleic acids, their ligand‐induced stabilization, and how conformational changes influence gene expression pose significant difficulties for experimental investigation. This work employs a combination of computational and experimental methods to address these challenges. By integrating structure‐informed design, crystallography, and machine learning‐augmented all‐atom molecular dynamics simulations (MD), we synthesized, biophysically and biochemically characterized, and studied the dissociation of a library of small molecule activators of the 5‐aminoimidazole–4–carboxamide ribonucleotide triphosphate (ZTP) riboswitch, a ligand‐binding RNA motif that regulates bacterial gene expression. We uncovered key interaction mechanisms, revealing valuable insights into the role of ligand binding kinetics on riboswitch activation. Further, we established that ligand on‐rates determine activation potency as opposed to binding affinity and elucidated RNA structural differences, which provide mechanistic insights into the interplay of RNA structure on riboswitch activation.

Chemistry↗

Thermal conductivity of MoS 2 monolayers from molecular dynamics simulations

Quantification of lattice thermal conductivity of two-dimensional semiconductors like MoS 2 is necessary for the design of electronic and thermoelectric devices, but direct experimental measurements on free-standing samples is challenging. Molecular dynamics simulations, with appropriate corrections, can provide a reference value for thermal conductivity for these material systems. Here, we construct a new empirical forcefield of the Stillinger-Weber form, parameterized to phonon dispersion relations, lattice constants and elastic moduli and we use it to compute a material-intrinsic thermal conductivity of 38.1 W/m-K at room temperature and estimate a maximum thermal conductivity of 85.4 W/m-K at T = 200 K. We also identify that phonon scattering by the large isotopic mass distribution of Mo and S contributes a significant correction (>45%) to the thermal conductivity at low temperatures.

36 MATERIALS SCIENCE↗

Molecular dynamics simulation of phase competition in terbium

The competition among multiple solid phases determines the final microstructures of a material. Such competition can originate at the very beginning of the solidification process. We report the results of molecular dynamics simulation of the phase competition between the hexagonal close-packed (hcp), face-centered cubic (fcc), and body-centered cubic (bcc) phases during the solidification of pure Tb. We found that the liquid supercooled below the hcp melting temperature has both bcc and hcp/fcc nuclei, but only the bcc nuclei grow such that the liquid always solidifies into the bcc phase, even at temperatures where the hcp phase is more stable. The hcp phase can only form in the last liquid droplet or at the bcc grain boundaries. Depending on the bcc grain orientations, the hcp phase jammed between the bcc grains either completely disappears or slowly grows via a solid-state massive transformation mechanism. Once the hcp phase becomes large enough, the stresses associated with its appearance can trigger a martensitic transformation. Yet, not the entire bcc phase is consumed by the martensitic transformation and the remaining bcc phase is transformed into the hcp phase via the solid-state massive transformation mechanism. If the supercooling is too large, the nucleation becomes almost barrier free and the liquid solidifies into a structure consisting of ultra-fine hcp and bcc grains after which the bcc phase quickly disappears.

36 MATERIALS SCIENCE↗