SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

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

At least 721 records · Page 40

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Si–Cl 2 –Ar + Atomic Layer Etching Window: A Fundamental Study Using Molecular Dynamics Simulations and a Reduced Order Model

Silicon (Si) atomic layer etching (ALE) by alternating exposure to chlorine gas (Cl 2 ) and argon ions (Ar + ) is studied by using molecular dynamics (MD) simulations and a reduced order model (ROM). Here, the purpose of this study is to elucidate the properties of the ALE window, a range of ion energies where the amount of Si etched over a series of cycles is nonzero and nearly independent of ion energy. Experimental studies of the Si–Cl 2 –Ar + ALE system report contradictory results related to the ALE window’s ion energy range. Both MD simulations and the ROM show that there is an ALE window present from approximately 15 to 20 eV for normal incidence argon ions. The Si–Cl 2 –Ar + system, therefore, exhibits a narrow ALE window. The amount of Si etched per cycle is less than one atomic layer because of the higher etch yield of Cl atoms relative to atomic Si and silicon chlorides. A modified version of the ROM with an artificially increased Si physical sputtering threshold energy expands the ALE window, illustrating the importance of the difference in chemical and physical sputtering threshold energies in the ALE window energy range. The ROM is also used to examine the dependence of the EPC on the Ar + ion fluence.

energy↗

DOE OSTI2018

Structure–Property of Lithium–Sulfur Nanoparticles via Molecular Dynamics Simulation

Lithium-sulfur (Li-S) batteries offer higher energy densities than most reported lithium-ion batteries. However, our understanding of Li-S battery is still largely unknown at the level of the nanoscale. Here, the structural properties of Li-S materials were investigated via molecular dynamics (MD) simulations using the ReaxFF force field. everal Li-S nanoparticles with different Li:S composition ratios (2:1 and 2:8) and various structures are studied. Our MD simulations show that, among the four structures we constructed for Li 2 S 8 nanoparticles, the core-shell structure is the most thermodynamically stable one during the charging (delithiation) process. In contrast to bulk crystal Li 2 S, we find the presence of mixed lithium sulfide and polysulfide species are common features for these Li-S (Li 2 S, Li 2 S 8 ) nanoparticles. The complex distribution of these sulfide and polysulfide speciation are dictated by both stoichiometry and local atomic structures in the nanoparticle. Finally, these findings will provide insight into further development of functionalized lithium-sulfur cathodes.

25 ENERGY STORAGE↗

NASA NTRS2016

Assessment of the Impacts of ACLS on the ISS Life Support System Using Dynamic Simulations in V-HAB

The Advanced Closed Loop System (ACLS) is currently under development by Airbus Defense and Space and is slated for launch to the International Space Station (ISS) in 2017. The addition of new hardware into an already complex system such as the ISS life support system (LSS) always poses operational risks. It is therefore important to understand the impacts ACLS will have on the existing systems to ensure smooth operations for the ISS. This analysis can be done by using dynamic computer simulations and one possible tool for such a simulation is the Virtual Habitat (V-HAB). Based on MATLAB, V-HAB has been under development at the Institute of Astronautics of the Technical University of Munich (TUM) since 2004 and in the past has been successfully used to simulate the ISS life support systems. The existing V-HAB ISS simulation model treated the interior volume of the space station as one large, ideally-stirred container. This model was improved to allow the calculation of the atmospheric composition inside individual modules of the ISS by splitting it into twelve distinct volumes. The virtual volumes are connected by a simulation of the inter-module ventilation flows. This allows for a combined simulation of the LSS hardware and the atmospheric composition aboard the ISS. A dynamic model of ACLS is added to the ISS Simulation and several different operating modes for both ACLS and the existing ISS life support systems are studied and the impacts of ACLS on the rest of the system are determined. The results suggest that the US, Russian and ACLS CO2 systems can operate at the same time without impeding each other. Furthermore, based on the results of this analysis, the US and ACLS Sabatier systems can be operated in parallel as well to a achieve a very low CO2 concentration in the cabin atmosphere.

Putz, Daniel↗

NASA NTRS2016

Assessment of the Impacts of ACLS on the ISS Life Support System using Dynamic Simulations in V-HAB

The Advanced Closed Loop System (ACLS) is currently under development by Airbus Defense and Space and is slated for launch to the International Space Station (ISS) in 2017. The addition of new hardware into an already complex system such as the ISS life support system (LSS) always poses operational risks. It is therefore important to understand the impacts ACLS will have on the existing systems to ensure smooth operations for the ISS. This analysis can be done by using dynamic computer simulations and one possible tool for such a simulation is Virtual Habitat (V-HAB). Based on Matlab (Registered Trademark) V-HAB has been under development at the Institute of Astronautics of the Technical University Munich (TUM) since 2006 and in the past has been successfully used to simulate the ISS life support systems. The existing V-HAB ISS simulation model treated the interior volume of the space station as one large ideally-stirred container. This model was improved to allow the calculation of the atmospheric composition inside the individual modules of the ISS by splitting it into ten distinct volumes. The virtual volumes are connected by a simulation of the inter-module ventilation flows. This allows for a combined simulation of the LSS hardware and the atmospheric composition aboard the ISS. A dynamic model of ACLS is added to the ISS simulation and different operating modes for both ACLS and the existing ISS life support systems are studied to determine the impacts of ACLS on the rest of the system. The results suggest that the US, Russian and ACLS CO2 systems can operate at the same time without impeding each other. Furthermore, based on the results of this analysis, the US and ACLS Sabatier systems can be operated in parallel as well to achieve the highest possible CO2 recycling together with a low CO2 concentration.

Puetz, Daniel↗

DOE OSTI2019

Amphiphilic peptide binding on crystalline vs. amorphous silica from molecular dynamics simulations

The leucine-lysine amphiphilic peptide LKα14 has been used to study fundamental driving forces in processes such as peptide-surface binding and biomineralization. Here, we employ molecular dynamics (MD) simulations in tandem with replica exchange metadynamics to probe the binding mechanism and thermodynamics of LKα14 on silica. We also investigate the effect that the nature of the silica surface – crystalline vs. amorphous, has on the binding properties and peptide-surface conformations. We find that water adsorbs differently on both surfaces; it forms a denser interfacial layer on the crystalline surface, compared to the amorphous surface. This causes the peptide to bind more strongly on the amorphous surface than the crystalline surface. Cluster analysis shows that the peptide adopts a helical conformation at both surfaces, with a greater distribution of states on the crystalline surface. Peptide binding is primarily through lysine interactions, in line with prior experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Early stage oxynitridation process of Si(001) surface by NO gas: Reactive molecular dynamics simulation study

Initial stage of oxynitridation process of Si substrate is of crucial importance in fabricating the ultrathin gate dielectric layer of high quality in advanced MOSFET devices. The oxynitridation reaction on a relaxed Si(001) surface is investigated via reactive molecular dynamics (MD) simulation. A total of 1120 events of a single nitric oxide (NO) molecule reaction at temperatures ranging from 300 to 1000 K are statistically analyzed. The observed reaction kinetics are consistent with the previous experimental or calculation results, which show the viability of the reactive MD technique to study the NO dissociation reaction on Si. We suggest the reaction pathway for NO dissociation that is characterized by the inter-dimer bridge of a NO molecule as the intermediate state prior to NO dissociation. Although the energy of the inter-dimer bridge is higher than that of the intra-dimer one, our suggestion is supported by the ab initio nudged elastic band calculations showing that the energy barrier for the inter-dimer bridge formation is much lower. The growth mechanism of an ultrathin Si oxynitride layer is also investigated via consecutive NO reactions simulation. The simulation reveals the mechanism of self-limiting reaction at low temperature and the time evolution of the depth profile of N and O atoms depending on the process temperature, which would guide to optimize the oxynitridation process condition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

Molecular Dynamics Simulations of Uranyl and Uranyl Carbonate Adsorption at Alumino-silicate Surfaces

Adsorption at mineral surfaces is a critical factor controlling the mobility of uranium(VI) in aqueous environments. Therefore, molecular dynamics (MD) simulations were performed to investigate uranyl(VI) adsorption onto two neutral alumino-silicate surfaces, namely the orthoclase (001) surface and the octahedral aluminum sheet of the kaolinite (001) surface. Although uranyl preferentially adsorbed as a bi-dentate innersphere complex on both surfaces, the free energy of adsorption at the orthoclase surface (-15 kcal mol-1) was significantly more favorable than that at the kaolinite surface (-3 kcal mol-1), which was attributed to differences in surface functional groups and to the ability of the orthoclase surface to dissolve a surface potassium ion upon uranyl adsorption. The structures of the adsorbed complexes compared favorably with X-ray absorption spectroscopy results. Simulations of the adsorption of uranyl complexes with up to three carbonate ligands revealed that uranyl complexes coordinated to up to 2 carbonate ions are stable on the orthoclase surface whereas uranyl carbonate surface complexes are unfavored at the kaolinite surface. Combining the MD-derived equilibrium adsorption constants for orthoclase with aqueous equilibrium constants for uranyl carbonate species indicates the presence of adsorbed uranium complexes with one or two carbonates in alkaline conditions, in support of current uranium(VI) surface complexation models.

Environmental Molecular Sciences Laboratory↗

DOE OSTI2021

Investigations of water/oxide interfaces by molecular dynamics simulations

Water/oxide interfaces are ubiquitous on earth and show significant influence on many chemical processes. For example, understanding water and solute adsorption as well as catalytic water splitting can help build better fuel cells and solar cells to overcome our looming energy crisis; the interaction between biomolecules and water/oxide interfaces is one hypothesis to explain the origin of life. However, knowledge in this area is still limited due to the difficulty of studying water/solid interfaces. As a result, research using increasingly sophisticated experimental techniques and computational simulations has been carried out in recent years. Although it is difficult for experimental techniques to provide detailed microscopic structural information, molecular dynamics (MD) simulations have satisfactory performance. In this report, we discuss classical and ab initio MD simulations of water/oxide interfaces. Generally, we are interested in the following questions: How do solid surfaces perturb interfacial water structure? How do interfacial water molecules and adsorbed solutes affect solid surfaces and how do interfacial environments affect solvent and solute behavior? Finally, we discuss progress in the application of neural network potential based MD simulations, which offer a promising future because this approach has already enabled ab initio level accuracy for very large systems and long trajectories.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

Spin-lattice dynamics simulation of external field effect on magnetic order of ferromagnetic iron

Modeling of field-induced magnetization in ferromagnetic materials has been an active topic in the last dozen years, yet a dynamic treatment of distance-dependent exchange integral has been lacking. In view of that, we employ spin-lattice dynamics (SLD) simulations to study the external field effect on magnetic order of ferromagnetic iron. Our results show that an external field can increase the inflection point of the temperature. Also the model provides a better description of the effect of spin correlation in response to an external field than the mean-field theory. An external field has a more prominent effect on the long range magnetic order than on the short range counterpart. Furthermore, an external field allows the magnon dispersion curves and the uniform precession modes to exhibit magnetic order variation from their temperature dependence.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Overcoming the Time Limitation in Molecular Dynamics Simulation of Crystal Nucleation: A Persistent-Embryo Approach

The crystal nucleation from liquid in most cases is too rare to be accessed within the limited time scales of the conventional molecular dynamics (MD) simulation. Here, we developed a “persistent embryo” method to facilitate crystal nucleation in MD simulations by preventing small crystal embryos from melting using external spring forces. We applied this method to the pure Ni case for a moderate undercooling where no nucleation can be observed in the conventional MD simulation, and obtained nucleation rate in good agreement with the experimental data. Moreover, the method is applied to simulate an even more sluggish event: the nucleation of the B2 phase in a strong glass-forming Cu-Zr alloy. The nucleation rate was found to be 8 orders of magnitude smaller than Ni at the same undercooling, which well explains the good glass formability of the alloy. In conclusion, our work opens a new avenue to study solidification under realistic experimental conditions via atomistic computer simulation.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Modelling dislocation-graphene interactions in a BCC Fe matrix by molecular dynamics simulations and gradient plasticity theory

We report graphene nanosheets (GNS) can enhance the strength and ductility of metal-based composites as they can obstruct the propagation of dislocations. The present article employs Molecular Dynamics (MD) simulations to investigate dislocation-GNS interaction mechanisms and possible influencing factors, including the number of GNS layers, the thickness of the metallic amorphous layer and the C - C bond strength. The results indicated that the shear strength of the metal/GNS interface and the bending stiffness of GNS determined the ability of GNS to block dislocation transmission. A physically based phenomenological parameter that can capture such dislocation-GNS interactions is the mechanical interface energy that has been put forth within gradient plasticity. By fitting the theoretical expressions to the simulation data, it was possible to obtain estimates for the mechanical interface energy for the GNS. It was found that increasing the GNS layers and adding an amorphous layer resulted in a strengthening in the stress–strain response and increased the value of this interfacial parameter. This indicates that the mechanical interfacial energy can be a unified measure for capturing and tuning the strength of various interfaces such as grain boundaries, GNS, amorphous-crystalline interface and bimetal interfaces.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Structural features of ISG borosilicate nuclear waste glasses revealed from high-energy X-ray diffraction and molecular dynamics simulations

Vitrification is a widely accepted method to immobilize nuclear waste. Detailed structural information is critical to understand the physical and chemical behaviors, including the long-term chemical durability of these glasses that are to be stored in geological storage sites. High-energy X-ray diffraction studies are used to obtain accurate structural information of the International Simple Glass (ISG), a six-component borosilicate model nuclear waste glass with a composition (mol%) of 60.2SiO 2 -16.0B 2 O 3 -12.6Na 2 O-3.8Al 2 O 3 -5.7CaO-1.7ZrO 2 . Classical molecular dynamics (MD) simulations with recently developed effective potentials are utilized to generate structural models of the ISG glasses to provide interpretation of the high-energy X-ray structural data and additional short and medium range structural details. The atomic structure model generated from MD simulations shows an excellent agreement with high-energy X-ray diffraction, measured for the first time for ISG, with an R x value of 5.5%. Systematic investigations of the medium-range structural features of ISG are performed, built upon recent study of the short-range structure. In particular, the glass former cation oxygen polyhedral analysis show that [SiO 4 ]-[SiO 4 ] and [BO 4 ]-[SiO 4 ] connections are the most abundant while [ZrO 6 ] octohedra having higher probability to connect to another [SiO 4 ], [BO 4 ] and [BO 3 ]. Moreover, we have studied the ISG surface structure using MD simulations. Noticeable differences in terms of both chemical compositions and structural features between the bulk and surface are observed. There exist increased small-sized rings, under coordinated Si species, high 3-coordinated boron content (86.5%) on the surface. Very importantly, sodium ions are found to be enriched on the glass surface, consistent with recent experimental results. Furthermore, both bulk and surface structural information are discussed in terms of the dissolution and corrosion mechanisms of ISG and related nuclear waste glasses.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Coarse-grained molecular dynamics simulation of the interface behaviour and self-assembly of CTAB cationic surfactants

Here, we study the behaviour at interfaces and the micelle self-assembly of a cationic surfactant (CTAB) by Molecular Dynamics (MD) simulations of coarse-grained models. We consider both the standard (with explicit water) Martini force field and the implicit solvent version of the Martini force field (Dry Martini). First, we study the behaviour of CTAB at a water/vacuum interface, at a water/organic solvent interface and in a pre-assembled CTAB micelle using both standard and Dry Martini and all-atomic simulations. Our results indicate that there are significant quantitative differences between the predictions of the two models. Interestingly, implicit solvent simulations with Dry Martini show good quantitative agreement with all-atomic MD simulations, better than explicit solvent Martini MD simulations. The computational efficiency of the Martini and Dry Martini models allowed us to study the self-assembly of CTAB in a large system with many micelles. We observe the self-assembly of CTAB into micelles and also the exchange of CTAB molecules between micelles by events such as micelle fusion and fission which are difficult to observe in all-atomic MD simulations due to the time and length scales involved. Under the studied conditions, both Martini models predict a rather different self-assembly behaviour. The standard Martini model predicts a final equilibrium state with spherical micelles with an average size of ≈70 CTAB molecules. In contrast, the Dry Martini model predicts the formation of large tubular micelles with ≈330 CTAB molecules. Compared with experiments, standard Martini and Dry Martini underestimate and overestimate, respectively, the micelle size.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

A classical molecular dynamics simulation method for the formation of “dry” gels from boro-aluminosilicate glass structures

In contact with water, glass transforms into amorphous and porous structures called gels. A simulation method based on classical molecular dynamics is proposed here to mimic “dry” gels forming from initial oxide glass structures. Six glass compositions were investigated. Two behaviours were evidenced depending on the initial glass composition, and in particular on the quantity of elements removed. If a large quantity of soluble elements (B, Na) was removed, it induced an increase in the average pore size within the gels, and the time needed to stabilise the gel structure increased because more local atomic rearrangements occurred. The gel network displayed a higher proportion of Si-Q4 at the expense of Si-Q3 and a lower average ring size compared to the glass network, irrespective of the glass composition. Surface effects were also highlighted in the dry gels, such as the presence of 3-coordinated Al and a decrease in the average angle Si-O-Si and Al-O-Al. Our findings will be compared to both wet gels and experimental data in further studies, to help find the best procedure to simulate such structures.

36 MATERIALS SCIENCE↗

DOE OSTI2014

Does fluoride disrupt hydrogen bond network in cationic lipid bilayer? Time-dependent fluorescence shift of Laurdan and molecular dynamics simulations

Time-dependent fluorescence shift (TDFS) of Laurdan embedded in phospholipid bilayers reports on hydration and mobility of the phospholipid acylgroups. Exchange of H{sub 2}O with D{sub 2}O prolongs the lifetime of lipid-water and lipid-water-lipid interactions, which is reflected in a significantly slower TDFS kinetics. Combining TDFS measurements in H{sub 2}O and D{sub 2}O hydrated bilayers with atomistic molecular dynamics (MD) simulations provides a unique tool for characterization of the hydrogen bonding at the acylgroup level of lipid bilayers. In this work, we use this approach to study the influence of fluoride anions on the properties of cationic bilayers composed of trimethylammonium-propane (DOTAP). The results obtained for DOTAP are confronted with those for neutral phosphatidylcholine (DOPC) bilayers. Both in DOTAP and DOPC H{sub 2}O/D{sub 2}O exchange prolongs hydrogen-bonding lifetime and does not disturb bilayer structure. These results are confirmed by MD simulations. TDFS experiments show, however, that for DOTAP this effect is cancelled in the presence of fluoride ions. We interpret these results as evidence that strongly hydrated fluoride is able to steal water molecules that bridge lipid carbonyls. Consequently, when attracted to DOTAP bilayer, fluoride disrupts the local hydrogen-bonding network, and the differences in TDFS kinetics between H{sub 2}O and D{sub 2}O hydrated bilayers are no longer observed. A distinct behavior of fluoride is also evidenced by MD simulations, which show different lipid-ion binding for Cl{sup −} and F{sup −}.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

High-temperature chromium diffusion in austenitic stainless steel: Ab initio molecular dynamics simulations

Chromium self-diffusion through stainless steel (SS) matrix and along grain boundaries is an important mechanism controlling SS structural materials corrosion. For this work, Cr diffusion in austenitic SS was simulated using canonical ab initio molecular dynamics with realistic models of type-316 SS bulk, with and without Cr vacancies, and a low-energy Σ3 twin boundary typically observed at active corrosion sites. Cr self-diffusion coefficients at 750 and 850 °C calculated using Einstein’s diffusion equation are 4.2 × 10 –6 and 8.1 × 10 –6 Å 2 ps –1 in pristine bulk, 3.8 × 10 –3 and 5.5 × 10 –3 Å 2 ps –1 in bulk including Cr vacancies, and 9.5 × 10 –2 and 1.0 × 10 –1 Å 2 ps –1 at a Σ3[1 1 1]60° twin boundary.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Performance and Feature Improvements in Parareal-based Power System Dynamic Simulation

In recent years, a novel Parareal-based approach has been developed for fast transient simulations of large power system interconnections. Parareal belongs to the class of Parallel-in-time algorithms for solution of systems of differential-algebraic equations in parallel over an interval of time. The selection of a reasonably fast and accurate coarse solution is crucial to improve the performance of Parareal algorithm. Semi-analytical solution methods are one promising approach to achieve this goal. They have been investigated, and some preliminary results are presented here. In addition, Parareal-based simulator has been expanded to enable co-simulation with OpenDSS, a widely used open-source distribution system simulator. Preserving the parallel nature of the Parareal approach and taking advantage of the parallel capabilities of the latest versions of OpenDSS, each distribution system can be solved in their entirety on different processors in parallel within the main Parareal simulator. This paper also presents the structure of the transmission and distribution co-simulation and some results with different dynamic models of inverter-based resources in the distribution systems.

Park, Byungkwon↗

DOE OSTI2019

An Overview of Recent Standard and Accelerated Molecular Dynamics Simulations of Helium Behavior in Tungsten

One of the most critical challenges for the successful adoption of nuclear fusion power corresponds to plasma-facing materials. Due to its favorable properties in this context (low sputtering yield, high thermal conductivity, high melting point, among others), tungsten is a leading candidate material. Nevertheless, tungsten is affected by the plasma and fusion byproducts. Irradiation by helium nuclei, in particular, strongly modifies the surface structure by a synergy of processes, whose origin is the nucleation and growth of helium bubbles. In this review, we present recent advances in the understanding of helium effects in tungsten from a simulational approach based on accelerated molecular dynamics, which emphasizes the use of realistic parameters, as are expected in experimental and operational fusion power conditions.

36 MATERIALS SCIENCE↗