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,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Understanding irradiation growth through atomistic simulations: Defect diffusion and clustering in α-zirconium and the influence of alloying elements - 2016-0093

Irradiation-induced structural changes of α-zirconium alloys and in particular the effect of iron were investigated by molecular dynamics simulations using embedded atom potentials derived from first-principles calculations. The simulations revealed that at temperatures between 500 and 600 K self-interstitial atoms (SIAs) diffuse rapidly in a cooperative movement, preferably parallel to basal planes (a directions; ), forming nano-clusters with an extension in and . Vacancies diffuse more slowly than SIAs and remain isolated for a longer period of time. Nano-clusters associated with SIAs cause a pronounced overall expansion in a directions, as well as local strains. Under compressive strain in the c direction, vacancy diffusivity increases in the c direction. In contrast, the diffusivity of SIAs increases in the c direction under a tensile strain in the c direction. SIA nano-clusters are highly mobile within basal planes. Vacancy clusters grow by merging, leading to a contraction in the a direction, compensating for the expansion caused by SIA nano-clusters and possibly contributing to the plateau in growth after the initial rapid expansion. At the onset of breakaway growth, possibly due to stress buildup, the vacancy nano-clusters can condense into c loops, thereby diminishing the compensation effect. The alloying elements iron, nickel, chromium, and niobium liberated from secondary phase particles under irradiation or already in solution are attracted to vacancies and SIAs and are found inside vacancy and SIA loops. The interaction of alloying elements with defect clusters is discussed, with a particular focus on iron. Iron has been found to promote cluster formation in zirconium, and the structures of zirconium-iron clusters have been analyzed. Tin is repelled by SIA clusters and only weakly attracted by vacancies. Niobium impedes the diffusion of SIAs (and therefore may increase annihilation rates with nearby vacancies) and does not destabilize vacancy or SIA clusters. Ab initio calculations of the dimensional and elastic coefficients of the intermetallic phases occurring in secondary phase particles, such as Zr{sub 2}Fe and Zr{sub 3}Fe, are presented, allowing an assessment of local strains in a zirconium matrix. Thus, novel results from extended molecular dynamics simulations provide new insights and contribute to a deeper understanding of the complex mechanisms causing irradiation-induced dimensional changes and the breakaway growth of zirconium alloys. (authors)

36 MATERIALS SCIENCE↗
NASA NTRS2010

Computational Fluid Dynamics (CFD) Simulations of Jet Mixing in Tanks of Different Scales

For long-duration in-space storage of cryogenic propellants, an axial jet mixer is one concept for controlling tank pressure and reducing thermal stratification. Extensive ground-test data from the 1960s to the present exist for tank diameters of 10 ft or less. The design of axial jet mixers for tanks on the order of 30 ft diameter, such as those planned for the Ares V Earth Departure Stage (EDS) LH2 tank, will require scaling of available experimental data from much smaller tanks, as well designing for microgravity effects. This study will assess the ability for Computational Fluid Dynamics (CFD) to handle a change of scale of this magnitude by performing simulations of existing ground-based axial jet mixing experiments at two tank sizes differing by a factor of ten. Simulations of several axial jet configurations for an Ares V scale EDS LH2 tank during low Earth orbit (LEO) coast are evaluated and selected results are also presented. Data from jet mixing experiments performed in the 1960s by General Dynamics with water at two tank sizes (1 and 10 ft diameter) are used to evaluate CFD accuracy. Jet nozzle diameters ranged from 0.032 to 0.25 in. for the 1 ft diameter tank experiments and from 0.625 to 0.875 in. for the 10 ft diameter tank experiments. Thermally stratified layers were created in both tanks prior to turning on the jet mixer. Jet mixer efficiency was determined by monitoring the temperatures on thermocouple rakes in the tanks to time when the stratified layer was mixed out. Dye was frequently injected into the stratified tank and its penetration recorded. There were no velocities or turbulence quantities available in the experimental data. A commercially available, time accurate, multi-dimensional CFD code with free surface tracking (FLOW-3D from Flow Science, Inc.) is used for the simulations presented. Comparisons are made between computed temperatures at various axial locations in the tank at different times and those observed experimentally. The affect of various modeling parameters on the agreement obtained are assessed.

Breisacher, Kevin↗
DOE OSTI2016

Structural crossover in a supercooled metallic liquid and the link to a liquid-to-liquid phase transition

Time-resolved synchrotron measurements were carried out to capture the structure evolution of an electrostatically levitated metallic-glass-forming liquid during free cooling. The experimental data shows a crossover in the liquid structure at ~1000 K, about 115 K below the melting temperature and 150 K above the crystallization temperature. The structure change is characterized by a dramatic growth in the extended-range order below the crossover temperature. Molecular dynamics simulations have identified that the growth of the extended-range order was due to an increased correlation between solute atoms. These results provide structural evidence for a liquid-to-liquid-phase-transition in the supercooled metallic liquid.

Glass transitions↗
DOE OSTI2016

Shock induced phase transition of water: Molecular dynamics investigation

Molecular dynamics simulations were carried out using numerous force potentials to investigate the shock induced phenomenon of pure bulk liquid water. Partial phase transition was observed at single shock velocity of 4.0 km/s without requirement of any external nucleators. Change in thermodynamic variables along with radial distribution function plots and spectral analysis revealed for the first time in the literature, within the context of molecular dynamic simulations, the thermodynamic pathway leading to formation of ice VII from liquid water on shock loading. The study also revealed information for the first time in the literature about the statistical time-frame after passage of shock in which ice VII formation can be observed and variations in degree of crystallinity of the sample over the entire simulation time of 100 ns.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Damage initiation and evolution in Al-Si layered microstructures under shock loading conditions at atomic scales

Designing materials with microstructural features like multiphase interfaces have shown significant promise for usage in the next generation defense and nuclear applications. The dynamic response of these interfaces in extreme environments is observed to vary with the deformability of individual phases, which could subsequently alter the favored damage nucleation sites related to spall failure. Here, we investigate the role of spacing of interfaces in a nanocrystalline layered Aluminum-Silicon system on the shock wave propagation behavior, and microstructural and defect evolution using classical molecular dynamics simulations. The simulations are carried out with different hypothetical Al/Si microstructures including variations in the distribution of Si grains as layers in a nanocrystalline Al matrix. The molecular dynamics simulations suggest that spall failure is preferentially initiated at Al/Si interfaces. In addition, for the same system volume and same concentration of Si, a higher number of layers is marked by an increase in shock wave velocity and reduced resistance to spall failure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Structural determinants of protocadherin-15 mechanics and function in hearing and balance perception

The vertebrate inner ear, responsible for hearing and balance, is able to sense minute mechanical stimuli originating from an extraordinarily broad range of sound frequencies and intensities or from head movements. Integral to these processes is the tip-link protein complex, which conveys force to open the inner-ear transduction channels that mediate sensory perception. Protocadherin-15 and cadherin-23, two atypically large cadherins with 11 and 27 extracellular cadherin (EC) repeats, are involved in deafness and balance disorders and assemble as parallel homodimers that interact to form the tip link. Here we report the X-ray crystal structure of a protocadherin-15 + cadherin-23 heterotetrameric complex at 2.9-Å resolution, depicting a parallel homodimer of protocadherin-15 EC1-3 molecules forming an antiparallel complex with two cadherin-23 EC1-2 molecules. In addition, we report structures for 10 protocadherin-15 fragments used to build complete high-resolution models of the monomeric protocadherin-15 ectodomain. Molecular dynamics simulations and validated crystal contacts are used to propose models for the complete extracellular protocadherin-15 parallel homodimer and the tip-link bond. Steered molecular dynamics simulations of these models suggest conditions in which a structurally diverse and multimodal protocadherin-15 ectodomain can act as a stiff or soft gating spring. These results reveal the structural determinants of tip-link–mediated inner-ear sensory perception and elucidate protocadherin-15’s structural and adhesive properties relevant in disease.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Real-Time Time-Dependent Density Functional Theory for Simulating Nonequilibrium Electron Dynamics

The explicit real-time propagation approach for time-dependent density functional theory (RT-TDDFT) has increasingly become a popular first-principles computational method for modeling various time-dependent electronic properties of complex chemical systems. Here, in this Perspective, we provide a nontechnical discussion of how this first-principles simulation approach has been used to gain novel physical insights into nonequilibrium electron dynamics phenomena in recent years. Following a concise overview of the RT-TDDFT methodology from a practical standpoint, we discuss our recent studies on the electronic stopping of DNA in water and the Floquet topological phase as examples. Our discussion focuses on how RT-TDDFT simulations played a unique role in deriving new scientific understandings. We then discuss existing challenges and some new advances at the frontier of RT-TDDFT method development for studying increasingly complex dynamic phenomena and systems.

Hamiltonians↗
DOE OSTI2016

An implicit divalent counterion force field for RNA molecular dynamics

How to properly account for polyvalent counterions in a molecular dynamics simulation of polyelectrolytes such as nucleic acids remains an open question. Not only do counterions such as Mg{sup 2+} screen electrostatic interactions, they also produce attractive intrachain interactions that stabilize secondary and tertiary structures. Here, we show how a simple force field derived from a recently reported implicit counterion model can be integrated into a molecular dynamics simulation for RNAs to realistically reproduce key structural details of both single-stranded and base-paired RNA constructs. This divalent counterion model is computationally efficient. It works with existing atomistic force fields, or coarse-grained models may be tuned to work with it. We provide optimized parameters for a coarse-grained RNA model that takes advantage of this new counterion force field. Using the new model, we illustrate how the structural flexibility of RNA two-way junctions is modified under different salt conditions.

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

Understanding supercooling mechanism in sodium sulfate decahydrate phase-change material

Salt hydrate-based phase-change materials are considered promising for future heat storage applications in residential heating/cooling systems. Smooth phase transition from the liquid to solid phase and vice versa is essential for effective heat exchanger; however, supercooling in salt hydrates delays the onset of liquid–solid phase transition. We investigate the molecular level mechanism of supercooling in sodium sulfate decahydrate (SSD). SSD is a complex salt hydrate whose properties are governed by electrostatic forces that include pure Coulombic interactions as well as hydrogen bonds. Experimentally, we examine the importance of a nucleator in reducing supercooling temperatures. We investigated the effect of various mass concentrations of a borax nucleator on a decrease of supercooling temperatures. Molecular dynamics simulation techniques are used to obtain a basic understanding of supercooling in SSD. We observe that by introducing borax as a nucleator, there is a decrease in the supercooling temperature before nucleation. Our molecular dynamics simulations show that long-range electrostatics between sodium and sulfate ion pairs and that with polar water molecules is responsible for delayed nucleation in SSD that results in supercooling, and also, dynamics of charged molecules slows down. The lack of crystallization leads to amorphous structures in supercooled SSD.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mean force kinetic theory applied to self-diffusion in supercritical Lennard-Jones fluids

A new kinetic-theory-based calculation of the self-diffusion coefficient for dense supercritical Lennard-Jones fluids is presented. The mean force kinetic theory, which was recently developed for transport in dense plasmas, is applied for the calculation of diffusion in dense neutral fluids. The calculation only requires the pair distribution function, a quantity that is readily calculable from equilibrium statistical mechanics for many systems, including the Lennard-Jones fluid. The self-diffusion coefficients are compared with calculations from molecular dynamics simulations, and good agreement at high density is demonstrated, even in the vicinity of the solid–fluid coexistence line. A comparison of different kinetic models with molecular dynamics simulations demonstrates that the transport coefficients have important contributions due to particle interaction via a potential of mean force and local correlations, which increase the collision rate. The new calculations compare well to those from free-volume theory and overcome a limitation of this theory that prevents its use in systems that interact via long range monotonic potentials. It is expected that this approach will also apply to other systems, including neutral-plasma and neutral-electrolyte mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Modulating electronic coupling at the quantum dot/molecule interface by wavefunction engineering

In this work, we use wavefunction engineering by varying the size of Quantum Dots (QDs) and tuning the delocalization (or diffuseness) of frontier orbitals of an acceptor molecule to modulate charge transfer dynamics at the QD/molecule interface. For this purpose, we apply our recently developed bulk-adjusted linear combination of atomic orbitals (BA-LCAO) approach for nanostructures and a density functional theory (DFT) for the acceptor molecules. These electronic structure calculations, combined with extensive molecular dynamics simulations using a fragmented molecular mechanics (FraMM) force field, reveal intimate details of charge transfer across the QD/Acceptor interface. For the spherical wurtzite-(CdSe) 201 and (CdSe) 693 nanostructures, as model QDs with respective 2.8 and 4.1 nm diameters, and anthraquinone-2,3-dicarboxylic acid and its derivatives with the 7-OH, 7-OF, 10-BH, and 10-CH 2 substituents, as model molecular acceptors, we find that (1) both the electron donating and withdrawing groups greatly enhance hole transfer by means of diffusing the acceptor HOMO; (2) electron transfer is affected only by the electron donating groups; (3) solvent effects are largely negligible for the orbital overlaps, and (4) consistent with spatial confinement theories, the electron density of the smaller QD penetrates farther into the vacuum than the corresponding density of the larger QD leading to stronger coupling with the acceptor. These findings suggest that (a) one can effectively control charge transfer across the QD/molecule interface by either changing the size of the QD or by tuning diffuseness of frontier orbitals of the acceptor molecule and (b) the combination of the recently developed BA-LCAO approach for QDs with a DFT for the acceptor molecules, facilitated by the use of the FraMM force field and extensive molecular dynamics simulations, provide qualitatively accurate description of charge transfer dynamics at the QD/acceptor interface.

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

Generalized 3D beam dynamics model for industrial traveling wave linacs design and simulations

Beam dynamics simulations in traveling wave (TW) accelerating structures with significant beam loading is a challenging problem. Some codes are capable of calculating the TW electromagnetic fields, and some can track particles through such fields, but most cannot treat both self-consistently. A few commercial codes can model the physics correctly at the expense of many processor-hour computations to obtain converged self-consistent solutions. However, simple, accurate equations of motion for intensive beam dynamics in TW accelerating structure analysis have been previously obtained by Masunov and subsequently implemented in the Hellweg code, which was developed at the Moscow Engineering-Physics Institute. Hellweg is based on equations that allow fast simulations of beam dynamics while taking into account such effects as beam loading, space charge, and external magnetic fields. In this paper, we describe in detail some recent improvements to the Hellweg physics kernel. We have generalized the Masunov results into a 3D set of equations of motion, which include all spatial components of the radiofrequency (RF) and external magnetic fields. We have also improved the Lapostolle space charge model to the general 3D ellipsoid form for any dimensions’ ratio, consideration for the particles outside the beam core, and the fields from the neighboring bunches that can exist in the real machine. These modifications allow approaching the Hellweg accuracy to the self-consistent commercial codes while keeping the simulation time short, which is essential during the linac design and optimization stage. The implementation of these new capabilities in Hellweg is carefully benchmarked against other codes and analytical calculations. The code is freely available with an open source license.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Coordinate transformation methodology for simulating quasistatic elastoplastic solids

Molecular dynamics simulations frequently employ periodic boundary conditions where the positions of the periodic images are manipulated in order to apply deformation to the material sample. For example, Lees-Edwards conditions use moving periodic images to apply simple shear. Here, we examine the problem of precisely comparing this type of simulation to continuum solid mechanics. We employ a hypo-elastoplastic mechanical model, and develop a projection method to enforce quasi-static equilibrium. We introduce a simulation framework that uses a fixed Cartesian computational grid on a reference domain, and imposes deformation via a time-dependent coordinate transformation to the physical domain. As a test case for our method, we consider the evolution of shear bands in a bulk metallic glass using the shear transformation zone theory of amorphous plasticity. We examine the growth of shear bands in simple shear and pure shear conditions as a function of the initial preparation of the bulk metallic glass.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Extended asymmetric hot region formation due to shockwave interactions following void collapse in shocked high explosive

In both continuum hydrodynamics simulations and also multimillion atom reactive molecular dynamics simulations of shockwave propagation in single crystal pentaerythritol tetranitrate (PETN) containing a cylindrical void, we observed the formation of an initial radially symmetric hot spot. By extending the simulation time to the nanosecond scale, however, we observed the transformation of the small symmetric hot spot into a longitudinally asymmetric hot region extending over a much larger volume. Performing reactive molecular dynamics shock simulations using the reactive force field (ReaxFF) as implemented in the LAMMPS molecular dynamics package, we showed that the longitudinally asymmetric hot region was formed by coalescence of the primary radially symmetric hot spot with a secondary triangular hot zone. We showed that the triangular hot zone coincided with a double-shocked region where the primary planar shockwave was overtaken by a secondary cylindrical shockwave. The secondary cylindrical shockwave originated in void collapse after the primary planar shockwave had passed over the void. A similar phenomenon was observed in continuum hydrodynamics shock simulations using the CTH hydrodynamics package. Furthermore, the formation and growth of extended asymmetric hot regions on nanosecond timescales has important implications for shock initiation thresholds in energetic materials.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Simulating Neutron Stars with Solid Quark Cores: Rotations, Oscillations, and Binary Mergers [Slides]

Neutron stars feature complexity from GR, QM, particle physics and nuclear physics effects which cannot be investigated in the lab up to the densities in their core. GWs from binary neutron star mergers are on the rise to provide more insight into the physics inside NSs. There is a possibility that involves crystalline, color-superconducting quark matter, with implications towards altering the GW signal if NSs harbor this exotic matter. There are no 3D dynamical simulations of solid quark matter inside Neutron stars. We perform the first 3D dynamical simulations of solid quark matter inside NSs. These are performed in the Newtonian Approximation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Biochemical and physiological flexibility accompanies reduced cellulose biosynthesis in Brachypodium cesa1 S830N

Here, we present a study into the mechanisms of primary cell wall cellulose formation in grasses, using the model cereal grass Brachypodium distachyon. The exon found adjacent to the BdCESA1 glycosyltransferase QXXRW motif was targeted using Targeting Induced Local Lesions in Genomes (TILLING) and sequencing candidate amplicons in multiple parallel reactions (SCAMPRing) leading to the identification of the Bdcesa1 S830N allele. Plants carrying this missense mutation exhibited a significant reduction in crystalline cellulose content in tissues that rely on the primary cell wall for biomechanical support. However, Bdcesa1 S830N plants failed to exhibit the predicted reduction in plant height. In a mechanism unavailable to eudicotyledons, B. distachyon plants homozygous for the Bdcesa1 S830N allele appear to overcome the loss of internode expansion anatomically by increasing the number of nodes along the stem. Stem biomechanics were resultantly compromised in Bdcesa1 S830N . The Bdcesa1 S830N missense mutation did not interfere with BdCESA1 gene expression. However, molecular dynamic simulations of the CELLULOSE SYNTHASE A (CESA) structure with modelled membrane interactions illustrated that Bdcesa1 S830N exhibited structural changes in the translated gene product responsible for reduced cellulose biosynthesis. Molecular dynamic simulations showed that substituting S830N resulted in a stabilizing shift in the flexibility of the class specific region arm of the core catalytic domain of CESA, revealing the importance of this motion to protein function.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2002

Simulating the Dynamics of Earth's Core: Using NCCS Supercomputers Speeds Calculations

If one wanted to study Earth's core directly, one would have to drill through about 1,800 miles of solid rock to reach liquid core-keeping the tunnel from collapsing under pressures that are more than 1 million atmospheres and then sink an instrument package to the bottom that could operate at 8,000 F with 10,000 tons of force crushing every square inch of its surface. Even then, several of these tunnels would probably be needed to obtain enough data. Faced with difficult or impossible tasks such as these, scientists use other available sources of information - such as seismology, mineralogy, geomagnetism, geodesy, and, above all, physical principles - to derive a model of the core and, study it by running computer simulations. One NASA researcher is doing just that on NCCS computers. Physicist and applied mathematician Weijia Kuang, of the Space Geodesy Branch, and his collaborators at Goddard have what he calls the,"second - ever" working, usable, self-consistent, fully dynamic, three-dimensional geodynamic model (see "The Geodynamic Theory"). Kuang runs his model simulations on the supercomputers at the NCCS. He and Jeremy Bloxham, of Harvard University, developed the original version, written in Fortran 77, in 1996.

Source record↗