Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic Simulation”

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

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

At least 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Six-degree-of-freedom aircraft simulation with mixed-data structure using the applied dynamics simulation language, ADSIM

A realistic simulation of an aircraft in the flight using the AD 100 digital computer is presented. The implementation of three model features is specifically discussed: (1) a large aerodynamic data base (130,00 function values) which is evaluated using function interpolation to obtain the aerodynamic coefficients; (2) an option to trim the aircraft in longitudinal flight; and (3) a flight control system which includes a digital controller. Since the model includes a digital controller the simulation implements not only continuous time equations but also discrete time equations, thus the model has a mixed-data structure.

Savaglio, Clare↗
DOE OSTI2014

Lenard-Balescu calculations and classical molecular dynamics simulations of electrical and thermal conductivities of hydrogen plasmas

Here, we present a discussion of kinetic theory treatments of linear electrical and thermal transport in hydrogen plasmas, for a regime of interest to inertial confinement fusion applications. In order to assess the accuracy of one of the more involved of these approaches, classical Lenard-Balescu theory, we perform classical molecular dynamics simulations of hydrogen plasmas using 2-body quantum statistical potentials and compute both electrical and thermal conductivity from out particle trajectories using the Kubo approach. Our classical Lenard-Balescu results employing the identical statistical potentials agree well with the simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

Carbon Nanotube Based Molecular Electronics and Motors: A View from Classical and Quantum Dynamics Simulations

The tubular forms of fullerenes popularly known as carbon nanotubes are experimentally produced as single-, multiwall, and rope configurations. The nanotubes and nanoropes have shown to exhibit unusual mechanical and electronic properties. The single wall nanotubes exhibit both semiconducting and metallic behavior. In short undefected lengths they are the known strongest fibers which are unbreakable even when bent in half. Grown in ropes their tensile strength is approximately 100 times greater than steel at only one sixth the weight. Employing large scale classical and quantum molecular dynamics simulations we will explore the use of carbon nanotubes and carbon nanotube junctions in 2-, 3-, and 4-point molecular electronic device components, dynamic strength characterization for compressive, bending and torsional strains, and chemical functionalization for possible use in a nanoscale molecular motor. The above is an unclassified material produced for non-competitive basic research in the nanotechnology area.

Srivastava, Deepak↗
DOE OSTI2023

Molecular Dynamics Simulations of Calcite Fracture in Water

Calcite (CaCO 3 ) is one of the most common minerals in geologic and engineered systems. It is often in contact with aqueous solutions, causing chemically assisted fracture that is critical to understanding the stability of subsurface systems and manmade structures. Calcite fracture was evaluated with reactive molecular dynamics simulations, including the impacts of crack tip geometry (notch), the presence of water, and surface hydroxyl groups. Chemo-mechanical weakening was assessed by comparing the loads where fracture began to propagate. Our analyses show that in the presence of a notch, the load at which crack growth begins is lower, compared to the effect of water or surface hydroxyls. Additionally, the breaking of two adjacent Ca–O bonds is the kinetic limitation for crack initiation, since transiently broken bonds can reform, not resulting in crack growth. In aqueous environments, fresh (not hydroxylated) calcite surfaces exhibited water strengthening. Manual addition of H + and/or OH – species on the (104) calcite surface resulted in chemo-mechanical weakening of calcite by 9%. Achieving full hydroxylation of the calcite surface was thermodynamically and kinetically limited, with only 0.17–0.01 OH/nm 2 surface hydroxylation observed on the (104) surface at the end of the simulations. In conclusion, the limited reactivity of pure water with the calcite surface restricts the chemo-mechanical effects and suggests that reactions between physiosorbed water and localized structural defects may be dominating the chemo-mechanical process in the studies where water weakening has been reported.

Calcite↗
DOE OSTI2023

Molecular dynamics simulations of fission gas xenon (Xe) diffusion at UO 2 grain-boundaries (Rev.1)

The diffusivity of fission gas xenon (Xe) at UO 2 grain-boundaries is one of the most important parameters in mechanistic modeling of fission gas diffusion in UO 2 based nuclear fuels. In this report, we use molecular dynamics simulations to investigate the Xe diffusivity in UO 2 grain-boundaries, employing the many-body potential developed by Cooper, Rushton and Grimes for UO 2 . Three different types of grain-boundaries are investigated, twist Σ5, tilt Σ5, and a random grain-boundary. Diffusion activation energies in the range of 0.39 – 1.46 eV are obtained for the Xe diffusivity. Comparison to results for the uranium vacancy diffusivity from MD simulations employing the same methodology suggests a weak to moderate attractive Xe-uranium vacancy binding energy depending on the grain-boundary type.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Investigating structure and dynamics of unentangled poly(dimethyl- co -diphenyl)siloxane via molecular dynamics simulation

Polysiloxane is one of the most important polymeric materials in technological use. Polydimethylsiloxane displays glass-like mechanical properties at low temperatures. Incorporation of phenyl siloxane, via copolymerization for example, improves not only the low-temperature elasticity but also enhances its performance over a wide range of temperatures. Copolymerization with the phenyl component can significantly change the microscopic properties of polysiloxanes, such as chain dynamics and relaxation. However, despite much work in the literature, the influence of such changes is still not clearly understood. In this work, we systematically study the structure and dynamics of random poly(dimethyl- co -diphenyl)siloxane via atomistic molecular dynamics simulations. As the molar ratio Φ of the diphenyl component increases, we find that the size of the linear copolymer chain expands. At the same time, the chain-diffusivity slows down by over an order of magnitudes. In conclusion, the reduced diffusivity appears to be a result of a complex interplay between the structural and dynamic changes induced by phenyl substitution.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Self-Healing Nanomaterials: Multimillion-Atom Reactive Molecular Dynamics Simulations

Organometal halide perovskites are attracting great attention as promising material for solar cells because of their high power conversion efficiency. The high performance has been attributed to the existence of free charge carriers and their large diffusion lengths, but the nature of carrier transport at the atomistic level remains elusive. Here, nonadiabatic quantum molecular dynamics simulations elucidate the mechanisms underlying the excellent free-carrier transport in CH 3 NH 3 PbI 3 . Pb and I sublattices act as disjunct pathways for rapid and balanced transport of photoexcited electrons and holes, respectively, while minimizing efficiency-degrading charge recombination. On the other hand, CH 3 NH 3 sublattice quickly screens out electrostatic electron-hole attraction to generate free carriers within 1 ps. Together this nano-architecture lets photoexcited electrons and holes dissociate instantaneously and travel far away to be harvested before dissipated as heat. As a result, this work provides much needed structure-property relationships and time-resolved information that potentially lead to rational design of efficient solar cells.

14 SOLAR ENERGY↗
DOE OSTI2015

Solvent controlled ion association in structured copolymers: Molecular dynamics simulations in dilute solutions

Tailoring the nature of individual segments within ion containing block co-polymers is one critical design tool to achieve desired properties. The local structure including the size and distribution of the ionic blocks, as well as the long range correlations, are crucial for their transport ability. Here, we present molecular dynamics simulations on the effects of varying the concentrations of the ionizable groups on the conformations of pentablock ionomer that consist of a center block of ionic sulfonated styrene tethered to polyethylene and terminated by a bulky substituted styrene in dilute solutions. Sulfonation fractions f (0 ≤ f ≤ 0.55), spanning the range from ionomer to polyelectrolytes, were studied. Results for the equilibrium conformation of the chains in water and a 1:1 mixture of cyclohexane and heptane are compared to that in implicit poor solvents with dielectric constants ε = 1.0 and 77.73. In water, the pentablock collapses with the sulfonated groups on the outer surface. As f increases, the ionic, center block increasingly segregates from the hydrophobic regions. In the 1:1 mixture of cyclohexane and heptane, the flexible blocks swell, while the center ionic block collapses for f > 0. For f = 0, all blocks swell. In both implicit poor solvents, the pentablock collapses into a nearly spherical shape for all f . In conclusion, the sodium counterions disperse widely throughout the simulation cell for both water and ε = 77.73, whereas for ε = 1.0 and mixture of cyclohexane and heptane, the counterions largely condense onto the collapsed pentablock.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Transient β-hairpin formation in α-synuclein monomer revealed by coarse-grained molecular dynamics simulation

Parkinson’s disease, originating from the intrinsically disordered peptide α-synuclein, is a common neurodegenerative disorder that affects more than 5% of the population above age 85. It remains unclear how α-synuclein monomers undergo conformational changes leading to aggregation and formation of fibrils characteristic for the disease. In the present study, we perform molecular dynamics simulations (over 180 μs in aggregated time) using a hybrid-resolution model, Proteins with Atomic details in Coarse-grained Environment (PACE), to characterize in atomic detail structural ensembles of wild type and mutant monomeric α-synuclein in aqueous solution. The simulations reproduce structural properties of α-synuclein characterized in experiments, such as secondary structure content, long-range contacts, chemical shifts, and {sup 3}J(H{sub N}H{sub C{sub α}})-coupling constants. Most notably, the simulations reveal that a short fragment encompassing region 38-53, adjacent to the non-amyloid-β component region, exhibits a high probability of forming a β-hairpin; this fragment, when isolated from the remainder of α-synuclein, fluctuates frequently into its β-hairpin conformation. Two disease-prone mutations, namely, A30P and A53T, significantly accelerate the formation of a β-hairpin in the stated fragment. We conclude that the formation of a β-hairpin in region 38-53 is a key event during α-synuclein aggregation. We predict further that the G47V mutation impedes the formation of a turn in the β-hairpin and slows down β-hairpin formation, thereby retarding α-synuclein aggregation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1996

Two-Dimensional Dynamic Simulation of a Continuous Foil Bearing

In this paper, the two dimensional(radial and circumferential) transient Navier-Stokes equations are used to solve the hydrodynamic problem in conjunction with the time dependent motion of the journal, and the deformable, spring supported foil. The elastic deformation of the foil and its supports are simulated by a finite element model. The time-dependent Navier-Stokes formulation is used to solve for the interaction between the fluid lubricant, the motion of the journal and the deformable foil boundary. The steady state, the quasi-transient and the full transient dynamic simulation of the foil-fluid journal interaction are examined on a comparative basis. For the steady state simulation, the fluid lubricant pressures are evaluated for a particular journal position, by means of an iterative scheme until convergence is achieved in both the fluid pressures and the corresponding foil deformation. For the quasi-transient case, the transient motion of the journal is calculated using a numerical integration scheme for the velocity and displacement of the journal. The deformation of the foil is evaluated through numerical iteration in feedback mode with the fluid film pressure generated by the journal motion until convergence at every time step is achieved. For the full transient simulation, a parallel real-time integration scheme is used to evaluate simultaneously the new journal position and the new deformed shape of the foil at each time step. The pressure of the fluid lubricant is iterated jointly with the corresponding journal position and the deformed foil geometry until convergence is achieved. A variable time-stepping Newmark-Beta integration procedure is used to evaluate the transient dynamics at each time step of the bearing.

Braun, M. Jack↗
DOE OSTI2021

Evaluation of Force Fields for Molecular Dynamics Simulations of Platinum in Bulk and Nanoparticle Forms

Here, understanding the size- and shape-dependent properties of platinum nanoparticles is critical for enabling the design of nanoparticle-based applications with optimal and potentially tunable functionality. Toward this goal, we evaluated nine different empirical potentials with the purpose of accurately modeling faceted platinum nanoparticles using molecular dynamics simulation. First, the potentials were evaluated by computing bulk and surface properties - surface energy, lattice constant, stiffness constants, and the equation of state - and comparing these to prior experimental measurements and quantum mechanics calculations. Then, the potentials were assessed in terms of the stability of cubic and icosahedral nanoparticles with faces in the {100} and {111} planes, respectively. Although none of the force fields predicts all the evaluated properties with perfect accuracy, one potential - the embedded atom method formalism with a specific parameter set - was identified as best able to model platinum in both bulk and nanoparticle forms.

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

Molecular Dynamics Simulation Study of the Early Stages of Nucleation of Iron Oxyhydroxide Nanoparticles in Aqueous Solutions

Nucleation is a fundamental step in crystal growth. Of environmental and materials relevance are reactions that lead to nucleation of iron oxyhydroxides in aqueous solutions. These reactions are difficult to study experimentally due to their rapid kinetics. Here, we used classical molecular dynamics simulations to investigate nucleation of iron hydroxide/oxyhydroxide nanoparticles in aqueous solutions. Results show that in a solution containing ferric ions and hydroxyl groups, iron–hydroxyl molecular clusters form by merging ferric monomers, dimers, and other oligomers, driven by strong affinity of ferric ions to hydroxyls. When deprotonation reactions are not considered in the simulations, these clusters aggregate to form small iron hydroxide nanocrystals with a six-membered ring-like layered structure allomeric to gibbsite. By comparison, in a solution containing iron chloride and sodium hydroxide, the presence of chlorine drives cluster assembly along a different direction to form long molecular chains (rather than rings) composed of Fe–O octahedra linked by edge sharing. Further, in chlorine-free solutions, when deprotonation reactions are considered, the simulations predict ultimate formation of amorphous iron oxyhydroxide nanoparticles with local atomic structure similar to that of ferrihydrite nanoparticles. Overall, our simulation results reveal that nucleation of iron oxyhydroxide nanoparticles proceeds via a cluster aggregation-based nonclassical pathway.

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

Assessing Melting and Solid–Solid Transition Properties of Choline Chloride via Molecular Dynamics Simulations

Choline chloride (ChCl) is used extensively as a hydrogen bond donor in deep eutectic solvents (DESs). However, determining its melting properties experimentally is challenging due to decomposition upon melting, leading to widely varying literature values. Accurate melting properties are crucial for understanding the solid–liquid phase behavior of ChCl-containing DESs. Here, we employ molecular dynamics simulations to compute the phase transitions of ChCl, testing a variety of atomistic force fields. We find that the results are sensitive to the choice of force field, but a melting temperature of 627 K and a melting enthalpy of 7.8 kJ/mol seem most reasonable, in good agreement with some literature values. Furthermore, we suggest these as the likely melting properties of ChCl, though the results are tentative due to limited experimental data for the liquid ChCl phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2017

Plausible Compositions of the Seven TRAPPIST-1 Planets Using Long-Term Dynamical Simulations

TRAPPIST-1 is a nearby ultra-cool dwarf that is host to a remarkable planetary system consisting of seven transiting planets. The orbital properties and radii of the planets have been well constrained, and recently, the masses of the inner six planets have been measured with additional ground- and space-based photometric observations. Large uncertainties in these mass measurements have prevented a robust analysis of the planetary compositions. Here, we perform many thousands of N-body dynamical simulations with planet properties perturbed from the observed values and identify those that are stable for millions of years. This allows us to identify self-consistent orbital solutions that can be used in future studies. From our range of dynamical masses, we find that most of the planets are consistent with an Earth-like composition, where TRAPPIST-1f is likely to have a volatile-rich envelope.

Quarles, B.↗
NASA NTRS2017

Plausible Compositions of the Seven TRAPPIST-1 Planets Using Long-Term Dynamical Simulations

TRAPPIST-1 is a nearby ultra-cool dwarf that is host to a remarkable planetary system consisting of seven transiting planets. The orbital properties and radii of the planets have been well constrained, and recently, the masses of the inner six planets have been measured with additional ground- and space-based photometric observations. Large uncertainties in these mass measurements have prevented a robust analysis of the planetary compositions. Here, we perform many thousands of N-body dynamical simulations with planet properties perturbed from the observed values and identify those that are stable for millions of years. This allows us to identify self-consistent orbital solutions that can be used in future studies. From our range of dynamical masses, we find that most of the planets are consistent with an Earth-like composition, where TRAPPIST-1f is likely to have a volatile-rich envelope.

B. Quarles↗
NASA NTRS2017

Plausible Compositions of the Seven TRAPPIST-1 Planets Using Long-Term Dynamical Simulations

TRAPPIST-1 is a nearby ultra-cool dwarf that is host to a remarkable planetary system consisting of seven transiting planets. The orbital properties and radii of the planets have been well constrained, and recently, the masses of the inner six planets have been measured with additional ground- and space-based photometric observations. Large uncertainties in these mass measurements have prevented a robust analysis of the planetary compositions. Here, we perform many thousands of N-body dynamical simulations with planet properties perturbed from the observed values and identify those that are stable for millions of years. This allows us to identify self-consistent orbital solutions that can be used in future studies. From our range of dynamical masses, we find that most of the planets are consistent with an Earth-like composition, where TRAPPIST-1f is likely to have a volatile-rich envelope.

B. Quarles↗
DOE OSTI2018

Predicting the Conformational Variability of Abl Tyrosine Kinase using Molecular Dynamics Simulations and Markov State Models

Understanding protein conformational variability remains a challenge in drug discovery. The issue arises in protein kinases, whose multiple conformational states can affect the binding of small-molecule inhibitors. To overcome this challenge, we propose a comprehensive computational framework based on Markov state models (MSMs). Our framework integrates the information from explicit-solvent molecular dynamics simulations to accurately rank-order the accessible conformational variants of a target protein. We tested the methodology using Abl kinase with a reference and blind-test set. Only half of the Abl conformational variants discovered by our approach are present in the disclosed X-ray structures. We report the approach successfully identified a protein conformational state not previously observed in public structures but evident in a retrospective analysis of Lilly in-house structures: the X-ray structure of Abl with WHI-P154. Using a MSM-derived model, the free energy landscape and kinetic profile of Abl was analyzed in detail highlighting opportunities for targeting the unique metastable states.

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

Assessment of the impact of reactor residence time distribution on non-equilibrium product selectivity of polypropylene pyrolysis using reactive molecular dynamics simulations

Mass residence time distribution (RTD) is considered to be an important factor controlling the product selectivity in the pyrolysis of biomass and plastic wastes along with the pyrolysis chemistry. However, due to the complex pyrolysis chemistry of biomass and plastic waste, the coupling between the reaction chemistry, RTD, and product selectivity is challenging to understand. Here, we introduce a reaction molecular dynamics-based method to examine pyrolysis chemistry and species timescales to assess the impact of RTD on product selectivity and yield. To validate this method, reactive molecular dynamics simulations were conducted for polypropylene pyrolysis and its non-equilibrium product selectivity using 6 different RTDs. We find that the RTD and the reaction chemistry control the peak non-equilibrium product concentrations. The peak monomer (C 3 H 6 ) concentration during pyrolysis can be increased by up to 25 % by using a narrow RTD in the case of polypropylene pyrolysis. We also find that product selectivity is strongly affected by the average residence time and RTD. This coupling between the reaction chemistry, RTD, and product selectivity highlights the need to understand detailed reaction chemistry to control RTD and optimize non-equilibrium product selectivity during polymer and biomass pyrolysis. The present method provides a new way to design RTD for reactors to reach maximized product selectivity of plastic waste and biomass.

33 ADVANCED PROPULSION SYSTEMS↗