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 1,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

The Nature of Hydrogen Adsorption on Platinum in the Aqueous Phase

H2 adsorption properties on Pt in aqueous phase was determined by a reaction kinetic measurement of gas H2 reacting with liquid D2O into liquid HDO, gas HD and D2, and a DFT based ab inito molecular dynamic simulation of H2 adsorption on Pt (111) in the presence and absence of water. Both approaches show that the H2 adsorption on Pt in water is much weaker than that in the gas phase. The presence of interfacial water destabilizes the adsorbed H on Pt via decreasing the adsorption heat (hydrogen binding energy) by 19 - 22 kJ·molH2-1 and inducing an extra entropy loss of 50 - 70 J·mol H2-1·K-1. Upon H adsorption on Pt, the closest water layer near Pt surface shifts to a larger distance from the surface and re-orientates to a more ordered structure; in turn, the water layer limits the translation mobility of adsorbed H on Pt. The influence of interfacial water property on H2 adsorption is also reflected by the decrease of adsorption heat with the decrease of local pH provided by the support. Research was funded by the Laboratory Directed Research and Development program at Pacific Northwest National Laboratory.

Yang, Guoju↗
DOE OSTI2017

Facet-Dependent Thermal Instability in LiCoO2

Thermal runaways triggered by the oxygen release from oxide cathode materials pose a major safety concern for widespread application of lithium ion batteries. Utilizing in situ aberration-corrected scanning transmission electron microscopy (STEM) and electron energy loss spectroscopy (EELS) at high temperatures, we show that oxygen release from LixCoO2 cathode crystals is occurring at the surface of particles. We correlated this local oxygen evolution from the LixCoO2 structure with local phase transitions spanning from layered to spinel and then to rock salt structure upon exposure to elevated temperatures. Ab initio molecular dynamics simulations (AIMD) results show that oxygen release is highly dependent on LixCoO2 facet orientation. While the [001] facets are stable at 300 degrees C, oxygen release is observed from the [012] and [104] facets, where under-coordinated oxygen atoms from the delithiated structures can combine and eventually evolve as O-2. The novel understanding that emerges from the present study piovides in-depth insights into the thermal runaway mechanism of Li-ion batteries and can assist the design and fabrication of cathode crystals with the most thermally stable facets.

ab initio molecular dynamics↗
DOE OSTI2021

Three-dimensional atomic packing in amorphous solids with liquid-like structure

Liquids and solids are two fundamental states of matter. However, our understanding of their three-dimensional atomic structure is mostly based on physical models. Here we use atomic electron tomography to experimentally determine the three-dimensional atomic positions of monatomic amorphous solids, namely a Ta thin film and two Pd nanoparticles. We observe that pentagonal bipyramids are the most abundant atomic motifs in these amorphous materials. Instead of forming icosahedra, the majority of pentagonal bipyramids arrange into pentagonal bipyramid networks with medium-range order. Molecular dynamics simulations further reveal that pentagonal bipyramid networks are prevalent in monatomic metallic liquids, which rapidly grow in size and form more icosahedra during the quench from the liquid to the glass state. Furthermore, these results expand our understanding of the atomic structures of amorphous solids and will encourage future studies on amorphous–crystalline phase and glass transitions in non-crystalline materials with three-dimensional atomic resolution.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shape-preserving amorphous-to-crystalline transformation of CaCO3 revealed by in situ TEM

Organisms use inorganic ions and macromolecules to regulate crystallization from amorphous- precursors, endowing natural biominerals with complex morphologies and enhanced properties. The mechanisms by which modifiers enable these shape preserving transformations are poorly understood. We used in situ liquid-phase transmission electron microscopy to follow the evolution from amorphous calcium carbonate to calcite in the presence of additives. A combination of contrast analysis and infrared spectroscopy shows that Mg ions, which are widely present in seawater and biological fluids, alter the transformation pathway in a concentration-dependent manner. The ions bring excess (structural) water into the amorphous bulk so that a direct transformation is triggered by dehydration in the absence of morphological changes. Molecular dynamics simulations suggest Mg-incorporated water induces structural fluctuations, allowing transformation without the need to nucleate a separate crystal. Thus, the obtained calcite retains the original morphology of the amorphous state, biomimetically achieving the morphological control of crystals seen in biominerals

crystallization, biomineralization, calcium carbon↗
DOE OSTI2016

Can We Describe Graphene Confined Water Structures as Overlapping of Approaching Graphene–Water Interfacial Structures?

We investigate the microscopic mechanisms of the overlap of interfacial structures in confined fluids and attempt to answer the question whether the confined structures can be predicted from the original density profiles of individual solid–fluid interfaces. For that purpose, we perform (globally) isobaric–isothermal (locally, grand canonical) molecular dynamics simulations to extract not only the axial distribution functions of the water sites for the uncoupled graphene–water interfaces, but also those corresponding to the confined aqueous environments over the interplate range $8 ≤ h (Å) ≤ 28$ typically at ambient conditions. We have tested two (i.e., an arithmetic and a geometric) superposition approximations for the singlet density of confined water between flat graphene plates, as well as for a combination of flat and corrugated graphene plates. The outcome of this study suggests that the answer to the title’s question is a “yes”, provided that the interplate distance h is large enough to avoid fluid geometric packing frustration.

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

Gel phase formation in dilute triblock copolyelectrolyte complexes

Assembly of oppositely charged triblock copolyelectrolytes into phase-separated gels at low polymer concentrations (<1% by mass) has been observed in scattering experiments and molecular dynamics simulations. Here we show that in contrast to uncharged, amphiphilic block copolymers that form discrete micelles at low concentrations and enter a phase of strongly interacting micelles in a gradual manner with increasing concentration, the formation of a dilute phase of individual micelles is prevented in polyelectrolyte complexation-driven assembly of triblock copolyelectrolytes. Gel phases form and phase separate almost instantaneously on solvation of the copolymers. Furthermore, molecular models of self-assembly demonstrate the presence of oligo-chain aggregates in early stages of copolyelectrolyte assembly, at experimentally unobservable polymer concentrations. Our discoveries contribute to the fundamental understanding of the structure and pathways of complexation-driven assemblies, and raise intriguing prospects for gel formation at extraordinarily low concentrations, with applications in tissue engineering, agriculture, water purification and theranostics.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Prebiotic hydrocarbon synthesis in impacting reduced astrophysical icy mixtures

We present results of prebiotic organic synthesis in shock-compressed reducing mixtures of simple ices from quantum molecular dynamics simulations extended to close to chemical equilibrium timescales. Given the relative abundance of carbon in reduced forms in astrophysical ices as well as the tendency of these mixtures to form complex hydrocarbons under the presence of external stimuli, it is possible that cometary impacts on a planetary surface could have yielded a larger array of prebiotic organic compounds than previously investigated. We find that the high pressures and temperatures due to shock compression yield a large assortment of carbon- and nitrogen-bonded extended structures that are highly reactive with short molecular lifetimes. Expansion and cooling causes these materials to break apart and form a wide variety of stable, potentially life-building compounds, including long-chain linear and branched hydrocarbons, large heterocyclic compounds, and a variety of different amines and exotic amino acids. Lastly, our results help provide a bottom-up understanding of hydrocarbon impact synthesis on the early Earth and its role in producing life-building molecules from simple starting materials.

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

Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes

Understanding the impact of extractant functionalization on metal-binding energetics in liquid-liquid extraction is essential to guide the development of better separation processes. Traditionally, computational extractant design uses electronic structure calculations on metal-ligand clusters to determine the metal-binding energy of the lowest energy state. Although highly accurate, this approach does not account for all of the relevant physics encountered under experimental conditions. Such methodologies often neglect entropic contributions such as temperature effects and ligand flexibility, in addition to approximating solvent-extractant interactions with implicit solvent models. In this study, we use classical molecular dynamics simulations with an advanced sampling method, metadynamics, to map out extractant molecule conformational free energies in the condensed phase. Here we generate the complete conformational landscape in solution for a family of bidentate malonamide-based extractants with different functionalizations of the headgroup and the side chains. In particular, we show how such alkyl functionalization reshapes the free energy landscape, affecting the free energy penalty of organizing the extractant into the cis-like metal-binding conformation from the trans-like conformation of the free extractant in solution. Specifically, functionalizing alkyl tails to the center of the headgroup has a greater influence on increasing molecular rigidity and disfavoring the binding conformation than functionalizing side chains. These findings are consistent with trends in metal-binding energetics based on experimentally reported distribution ratios. We also consider a different bidentate extractant molecule, carbamoylmethylphosphine oxide, and show how the choice of solvent can further reshape the conformational energetic landscape. This study demonstrates the feasibility of using molecular dynamics simulations with advanced sampling techniques to investigate extractant conformational energetics in solution, which, more broadly, will enable extractant design that accounts for entropic effects and explicit solvation.

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

Electronic structure of aqueous solutions: Bridging the gap between theory and experiments

Predicting the electronic properties of aqueous liquids has been a long-standing challenge for quantum mechanical methods. However, it is a crucial step in understanding and predicting the key role played by aqueous solutions and electrolytes in a wide variety of emerging energy and environmental technologies, including battery and photoelectrochemical cell design. Here, we propose an efficient and accurate approach to predict the electronic properties of aqueous solutions, on the basis of the combination of first-principles methods and experimental validation using state-of-the-art spectroscopic measurements. We present results of the photoelectron spectra of a broad range of solvated ions, showing that first-principles molecular dynamics simulations and electronic structure calculations using dielectric hybrid functionals provide a quantitative description of the electronic properties of the solvent and solutes, including excitation energies. The proposed computational framework is general and applicable to other liquids, thereby offering great promise in understanding and engineering solutions and liquid electrolytes for a variety of important energy technologies.

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

Experimental and Numerical Investigation of a Miniature Additively Manufactured Vortex Tube

Ranque–Hilsch vortex tubes are simple devices that can produce a cooling effect using compressed air. A key advantage of vortex tubes is the lack of moving solid parts; however, their efficiencies are relatively low. The present study focuses on the development of a miniature variable-diameter tube using additive manufacturing. A metal-based 3D printing technique was utilized to fabricate this vortex tube monolithically. Computational fluid dynamics simulations employing software STAR-CCM + with a compressible Reynolds-Averaged Navier–Stokes (RANS) approach and the elliptic-blending lag k-epsilon turbulence model have been applied to model thermofluid processes inside the vortex tube, to good agreement with the experiment. In conclusion, a temperature decrease of 13.3 °C and a cooling power of approximately 4 W were experimentally achieved with a pressure ratio of 4 in the air at normal conditions. This result shows promise for the goal of utilizing additive manufacturing to design and build complex-geometry vortex tubes intended for use with cryogenic fluids.

08 HYDROGEN↗
DOE OSTI2020

Facile diamond synthesis from lower diamondoids

Carbon-based nanomaterials have exceptional properties that make them attractive for a variety of technological applications. Here, we report on the use of diamondoids (diamond-like, saturated hydrocarbons) as promising precursors for laser-induced high-pressure, high-temperature diamond synthesis. The lowest pressure and temperature $(P-T)$ conditions that yielded diamond were 12 GPa (at ~2000 K) and 900 K (at ~20 GPa), respectively. This represents a substantially reduced transformation barrier compared with diamond synthesis from conventional (hydro)carbon allotropes, owing to the similarities in the structure and full sp 3 hybridization of diamondoids and bulk diamond. At 20 GPa, diamondoid-to-diamond conversion occurs rapidly within <19 μs. Molecular dynamics simulations indicate that once dehydrogenated, the remaining diamondoid carbon cages reconstruct themselves into diamond-like structures at high $P-T$. This study is the first successful mapping of the $P-T$ conditions and onset timing of the diamondoid-to-diamond conversion and elucidates the physical and chemical factors that facilitate diamond synthesis.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Effects of carotenoids on lipid bilayers

Carotenoids have been found to be important in improving the integrity of biomembranes in eukaryotes. However, the molecular details of how carotenoids modulate the physical properties of biomembranes are unknown. To this end, we have conducted a series of molecular dynamics simulations of different biologically-relevant membranes in the presence of carotenoids. The carotenoid effect on the membrane was found to be specific to the identity of the carotenoid and the composition of the membrane itself. Therefore, different classes of carotenoids produce a different effect on the membrane, and different membrane phases are affected differently by carotenoids. It is apparent from our data that carotenoids do trigger the bilayer to become thinner. The mechanism by which this occurs depends on two competing factors, the ability of the lipid tails of opposing monolayers to either (1) compress or (2) interdigitate as the bilayer condenses. Indeed, carotenoids directly influence the physical properties via these two mechanisms, thus compacting the bilayer. Furthermore, the degree to which these competing mechanisms are utilized depends on the bilayer phase and the carotenoid identity.

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

Crystal structure of the human NK 1 tachykinin receptor

The NK 1 tachykinin G-protein–coupled receptor (GPCR) binds substance P, the first neuropeptide to be discovered in mammals. Through activation of NK 1 R, substance P modulates a wide variety of physiological and disease processes including nociception, inflammation, and depression. Human NK 1 R (hNK 1 R) modulators have shown promise in clinical trials for migraine, depression, and emesis. However, the only currently approved drugs targeting hNK 1 R are inhibitors for chemotherapy-induced nausea and vomiting (CINV). To better understand the molecular basis of ligand recognition and selectivity, we solved the crystal structure of hNK 1 R bound to the inhibitor L760735, a close analog of the drug aprepitant. Our crystal structure reveals the basis for antagonist interaction in the deep and narrow orthosteric pocket of the receptor. We used our structure as a template for computational docking and molecular-dynamics simulations to dissect the energetic importance of binding pocket interactions and model the binding of aprepitant. Furthermore, the structure of hNK 1 R is a valuable tool in the further development of tachykinin receptor modulators for multiple clinical applications.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2010

Control Activity in Support of NASA Turbine Based Combined Cycle (TBCC) Research

Control research for a Turbine Based Combined Cycle (TBCC) propulsion system is the current focus of the Hypersonic Guidance, Navigation, and Control (GN&C) discipline team. The ongoing work at the NASA Glenn Research Center (GRC) supports the Hypersonic GN&C effort in developing tools to aid the design of control algorithms to manage a TBCC airbreathing propulsion system during a critical operating period. The critical operating period being addressed in this paper is the span when the propulsion system transitions from one cycle to another, referred to as mode transition. One such tool, that is a basic need for control system design activities, is computational models (hereto forth referred to as models) of the propulsion system. The models of interest for designing and testing controllers are Control Development Models (CDMs) and Control Validation Models (CVMs). CDMs and CVMs are needed for each of the following propulsion system elements: inlet, turbine engine, ram/scram dual-mode combustor, and nozzle. This paper presents an overall architecture for a TBCC propulsion system model that includes all of the propulsion system elements. Efforts are under way, focusing on one of the propulsion system elements, to develop CDMs and CVMs for a TBCC propulsion system inlet. The TBCC inlet aerodynamic design being modeled is that of the Combined-Cycle Engine (CCE) Testbed. The CCE Testbed is a large-scale model of an aerodynamic design that was verified in a small-scale screening experiment. The modeling approach includes employing existing state-of-the-art simulation codes, developing new dynamic simulations, and performing system identification experiments on the hardware in the NASA GRC 10 by10-Foot Supersonic Wind Tunnel. The developed CDMs and CVMs will be available for control studies prior to hardware buildup. The system identification experiments on the CCE Testbed will characterize the necessary dynamics to be represented in CDMs for control design. These system identification models will also be the reference models to validate the CDM and CVM models. Validated models will give value to the tools used to develop the models.

Stueber, Thomas J.↗
DOE OSTI2023

Gaussian approximation potentials for accurate thermal properties of two-dimensional materials

Two-dimensional materials (2DMs) continue to attract a lot of attention, particularly for their extreme flexibility and superior thermal properties. Molecular dynamics simulations are among the most powerful methods for computing these properties, but their reliability depends on the accuracy of interatomic interactions. While first principles approaches provide the most accurate description of interatomic forces, they are computationally expensive. In contrast, classical force fields are computationally efficient, but have limited accuracy in interatomic force description. Machine learning interatomic potentials, such as Gaussian Approximation Potentials, trained on density functional theory (DFT) calculations offer a compromise by providing both accurate estimation and computational efficiency. Here, in this work, we present a systematic procedure to develop Gaussian approximation potentials for selected 2DMs, graphene, buckled silicene, and h-XN (X = B, Al, and Ga, as binary compounds) structures. We validate our approach through calculations that require various levels of accuracy in interatomic interactions. The calculated phonon dispersion curves and lattice thermal conductivity, obtained through harmonic and anharmonic force constants (including fourth order) are in excellent agreement with DFT results. HIPHIVE calculations, in which the generated GAP potentials were used to compute higher-order force constants instead of DFT, demonstrated the first-principles level accuracy of the potentials for interatomic force description. Molecular dynamics simulations based on phonon density of states calculations, which agree closely with DFT-based calculations, also show the success of the generated potentials in high-temperature simulations.

2D Materials↗
DOE OSTI2016

Design of specially adapted reactive coordinates to economically compute potential and kinetic energy operators including geometry relaxation

Quantum dynamics simulations require prior knowledge of the potential energy surface as well as the kinetic energy operator. Typically, they are evaluated in a low-dimensional subspace of the full configuration space of the molecule as its dimensionality increases proportional to the number of atoms. This entails the challenge to find the most suitable subspace. We present an approach to design specially adapted reactive coordinates spanning this subspace. In addition to the essential geometric changes, these coordinates take into account the relaxation of the non-reactive coordinates without the necessity of performing geometry optimizations at each grid point. The method is demonstrated for an ultrafast photoinduced bond cleavage in a commonly used organic precursor for the generation of electrophiles. The potential energy surfaces for the reaction as well as the Wilson G-matrix as part of the kinetic energy operator are shown for a complex chemical reaction, both including the relaxation of the non-reactive coordinates on equal footing. A microscopic interpretation of the shape of the G-matrix elements allows to analyze the impact of the non-reactive coordinates on the kinetic energy operator. Additionally, we compare quantum dynamics simulations with and without the relaxation of the non-reactive coordinates included in the kinetic energy operator to demonstrate its influence.

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

Nanometer-Scale Correlations in Aqueous Salt Solutions

The intermediate-range (1–4 nm) correlation of cations, an- ions, and water in aqueous alkaline earth salt solutions is measured using synchrotron X-ray diffraction. We differentiate from the entire solution structure factor, S(Q), a separate region at low Q (<1 °A-1) containing local diffraction maxima (pre-peaks) that indicate nanometer-scale oscillatory behavior. These features are quantitatively reproduced with classical molecular dynamics simulations. At high concentrations, the pre-peaks emerge from correlations arising from existence of small quasi close-packed lattice- like structures comprised of cation hydration spheres. We also analyze the concentration dependence of the pre-peak and discuss the overall results in light of the rich literature dealing with intermediate-range correlations underlying universal phenomena in concentrated electrolytes. The authors were supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences. Pacific Northwest National Laboratory (PNNL) is operated by Battelle for the U.S. Department of Energy under Contract No. DE-AC05-76RL01830. Computing resources were generously allocated by PNNL's Institutional Computing program. DOE Contract No. DEAC02- 06CH11357 supports operation of the Advanced Photon Source at Argonne National Laboratory.

Chemical physics, nucleation, Madelung potential↗