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,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Structural evolution in Au- and Pd-based metallic glass forming liquids and the case for improved molecular dynamics force fields

The results of a combined experimental and computational investigation of the structural evolution of Au 81 Si 19 , Pd 82 Si 18 , and Pd 77 Cu 6 Si 17 metallic glass forming liquids are presented. Electrostatically levitated metallic liquids are prepared, and synchrotron x-ray scattering studies are combined with embedded atom method molecular dynamics simulations to probe the distribution of relevant structural units. Metal–metalloid based metallic glass forming systems are an extremely important class of materials with varied glass forming ability and mechanical processibility. High quality experimental x-ray scattering data are in poor agreement with the data from the molecular dynamics simulations, demonstrating the need for improved interatomic potentials. The first peak in the x-ray static structure factor in Pd 77 Cu 6 Si 17 displays evidence for a Curie–Weiss type behavior but also a peak in the effective Curie temperature. In conclusion, a proposed order parameter distinguishing glass forming ability, 1/(S(T,q 1 )−1, shows a peak in the effective Curie temperature near a crossover temperature established by the behavior of the viscosity, T A .

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

Free Surface Relaxations of Star-Shaped Polymer Films

Here, the surface relaxation dynamics of supported star-shaped polymer thin films are shown to be slower than the bulk, persisting up to temperatures at least 50 degrees above the bulk glass transition temperature T$bulk\atop{g}$. This behavior, exhibited by star-shaped polystyrenes (SPSs) with functionality f = 8 arms and molecular weights per arm M arm < M e (M e is the entanglement molecular weight), is shown by molecular dynamics simulations to be associated with a preferential localization of these macromolecules at the free surface. This new phenomenon is in notable contrast to that of linear chain polymer thin film systems where the surface relaxations are enhanced in relation to the bulk; this enhancement persists only for a limited temperature range above the bulk T$bulk\atop{g}$. Finally, evidence of the slow surface dynamics, compared to the bulk, for temperatures well above T g and at length and time scales not associated with the glass transition has not previously been reported for polymers.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Beyond nothingness in the formation and functional relevance of voids in polymer films

Abstract Voids—the nothingness—broadly exist within nanomaterials and impact properties ranging from catalysis to mechanical response. However, understanding nanovoids is challenging due to lack of imaging methods with the needed penetration depth and spatial resolution. Here, we integrate electron tomography, morphometry, graph theory and coarse-grained molecular dynamics simulation to study the formation of interconnected nanovoids in polymer films and their impacts on permeance and nanomechanical behaviour. Using polyamide membranes for molecular separation as a representative system, three-dimensional electron tomography at nanometre resolution reveals nanovoid formation from coalescence of oligomers, supported by coarse-grained molecular dynamics simulations. Void analysis provides otherwise inaccessible inputs for accurate fittings of methanol permeance for polyamide membranes. Three-dimensional structural graphs accounting for the tortuous nanovoids within, measure higher apparent moduli with polyamide membranes of higher graph rigidity. Our study elucidates the significance of nanovoids beyond the nothingness, impacting the synthesis‒morphology‒function relationships of complex nanomaterials.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Effects of electronic excitation on cascade dynamics in nickel–iron and nickel–palladium systems

Using molecular dynamics simulations and the two-temperature model, we provide in this paper a comparison of the surviving damage from single ion irradiation events in nickel-based alloys, for cascades with and without taking into account the effects of the electronic excitations. We find that including the electronic effects impacts the amount of the resulting damage and the production of isolated defects. Finally, irradiation of nickel–palladium systems results in larger numbers of defects compared to nickel–iron systems, with similar numbers of isolated defects. We additionally investigate the mass effect on the two-temperature model in molecular dynamics simulations of cascades.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Distance-dependent dielectric constant at the calcite/electrolyte interface: Implication for surface complexation modeling

Hypothesis: The electrical double layer formed at the mineral/electrolyte interface is often modeled using mean-field approaches based on a continuum description of the solvent whose dielectric constant is assumed to decrease monotonically with decreasing distance to the surface. In contrast, molecular simulations show that the solvent polarizability oscillates near the surface similar to the water density profile - as shown previously, for example, by Bonthuis et al. (D.J. Bonthuis, S. Gekle, R.R. Netz, Dielectric Profile of Interfacial Water and its Effect on Double-Layer Capacitance, Phys Rev Lett 107(16) (2011) 166102). Here, we showed that molecular and mesoscale pictures agree by spatially averaging the dielectric constant obtained from molecular dynamics simulations over the distances relevant to the mean-field representation. In addition, the values of capacitances used to describe the electrical double layer in Surface Complexation Models (SCMs) of the mineral/electrolyte interface can be estimated using molecularly informed spatially averaged dielectric constants and positions of hydration layers. Experiments: First, we used molecular dynamics simulations to model the calcite 101¯4/electrolyte interface. Next, by using atomistic trajectories, we calculated the distance-dependent static dielectric constant and water density in the direction normal to the. Finally, we applied spatial compartmentalization consistent with the model of parallel-plate capacitors connected in series to estimate SCM capacitances. Findings: Computationally expensive simulations are required to determine the dielectric constant profile of interfacial water near the mineral surface. On the other hand, water density profiles are readily assessable from much shorter simulation trajectories. Our simulations confirmed that dielectric and water density oscillations at the interface are correlated. Here, we parametrized linear regression models to estimate the dielectric constant directly from the local water density. This is a significant computational shortcut compared to slowly converging calculations relying on total dipole moment fluctuations. The amplitude of the interfacial dielectric constant oscillation can exceed the dielectric constant of the bulk water, suggesting an ice-like frozen state, but only if there are no electrolyte ions. The interfacial accumulation of electrolyte ions causes a decrease in the dielectric constant due to the reduction of water density and re-orientation of water dipoles in ion hydration shells. Finally, we show how to use the computed dielectric properties to estimate SCM's capacitances.

58 GEOSCIENCES↗
DOE OSTI2016

Predicting kinetic nanocrystal shapes through multi-scale theory and simulation: Polyvinylpyrrolidone-mediated growth of Ag nanocrystals

In the shape-controlled synthesis of colloidal Ag nanocrystals, structure-directing agents, particularly polyvinylpyrrolidone (PVP), are known to be a key additive in making nanostructures with well-defined shapes. Although many Ag nanocrystals have been successfully synthesized using PVP, the mechanism by which PVP actuates shape control remains elusive. Here, we present a multi-scale theoretical framework for kinetic Wulff shape predictions that accounts for the chemical environment, which we used to probe the kinetic influence of the adsorbed PVP film. Within this framework, we use umbrella-sampling molecular dynamics simulations to calculate the potential of mean force and diffusion coefficient profiles of Ag atom deposition onto Ag(100) and Ag(111) in ethylene glycol solution with surface-adsorbed PVP. We use these profiles to calculate the mean-first passage times and implement extensive Brownian dynamics simulations, which allows the kinetic effects to be quantitatively evaluated. Our results show that PVP films can regulate the flux of Ag atoms to be greater towards Ag(111) than Ag(100). PVP’s preferential binding towards Ag(100) over Ag(111) gives PVP its flux-regulating capabilities through the lower free-energy barrier of Ag atoms to cross the lower-density PVP film on Ag(111) and enhanced Ag trapping by the extended PVP film on Ag(111). Under kinetic control, {100}-faceted nanocrystals will be formed when the Ag flux is greater towards Ag(111). Furthermore, the predicted kinetic Wulff shapes are in agreement with the analogous experimental system.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Interactions between displacement cascade and dislocation and their influences on Peierls stress in Fe-20Cr-25Ni alloys

The interactions between dislocations and displacement cascades are investigated using molecular dynamics simulations. The primary knock-on atoms (PKAs) are placed at different distances to edge and screw dislocation lines. Specifically, an fcc Fe-20Cr-25Ni system is studied; the alloy is the based alloy for Alloy 709. The number of defects is calculated using Wigner-Seitz analysis in order to illustrate the effect of dislocations on cascades. The displacement cascades in systems containing dislocations tend to generate more total defects compared with bulk systems, but eventually lead to less surviving defects in the materials matrix due to interactions with dislocations. The changes in atomic structures of the dislocations after interacting with cascades are analyzed to understand how the primary damage affects dislocations. The displacement cascade potentially causes dislocation climb in edge dislocations and cross-slip in screw dislocations, which could serve as additional mechanisms contributing to the radiation-induced hardening (RIH) compared with conventional RIH models. To reveal this, Peierls stresses are calculated before and after cascades using molecular static simulations. This study provides critical information of how defect production and dislocations are correlated, especially when total dose is high, which needs to be taken into account in upper scale models, such as mean-field rate theory.

Lu, Zizhe↗
DOE OSTI2018

A comparative study on the oxidation of two-dimensional Ti 3 C 2 MXene structures in different environments

The oxidation of two-dimensional Ti 3 C 2 MXene structures has been recognized as a promising method for the formation of hybrid structures of carbon supported nanotitania. Here, we studied the oxidation of Ti 3 C 2 MXene structures using reactive force field (ReaxFF) based molecular dynamics simulations. To investigate the effect of oxidizing agent, we used three different environments including dry air (oxygen molecules), wet air (oxygen and water molecules) and hydrogen peroxide (H 2 O 2 molecules). Oxidation simulations were performed at different temperatures of 1000, 1500, 2000, 2500 and 3000 K. We found that by controlling the temperature, carbon supported titania can be formed by diffusion of Ti atoms to the surface of MXene structures during the oxidation. These results were confirmed by the increase in the average bond orders of C–C and Ti–O and decrease in the average bond orders of Ti–C. Moreover, it was revealed that by increasing the temperature, the oxidation rate increases, and the rate depends on the oxidant according to the following order: H 2 O 2 > wet air > dry air. This was validated by heating the MXene in wet air and dry air followed by examining the products using X-ray diffraction and Raman spectroscopy. To investigate the effect of the presence of the oxidant, the MXene structure was also heated under vacuum, and it was found that in this case the MXene structure was converted to cubic TiC. Furthermore, it was indicated that during the heating process of the MXene structure in a vacuum, the average bond orders for C–C, Ti–C and Ti–O do not change significantly due to the lack of oxidants in the environment, leaving room for only topotactic transformation of atoms during the heating. In conclusion, the formation of the TiC structure by heating the MXene in an inert environment was confirmed experimentally as well.

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

Application of Pinniped Vibrissae to Aeropropulsion

Vibrissae of Phoca Vitulina (Harbor Seal) and Mirounga Angustirostris (Elephant Seal) possessundulations along their length. Harbor Seal Vibrissae were shown to reduce vortex induced vibrations and reduce dragcompared to appropriately scaled cylinders and ellipses. Samples of Harbor Seal vibrissae, Elephant Seal vibrissae andCalifornia Sea Lion vibrissae were collected from the Marine Mammal Center in California. CT scanning, microscopy and3D scanning techniques were utilized to characterize the whiskers. Computational fluid dynamics simulations of thewhiskers were carried out to compare them to an ellipse and a cylinder. Leading edge parameters from the whiskerswere used to create a 3D profile based on a modern power turbine blade. The NASA SW-2 facility was used to performwind tunnel cascade testing on the 'Seal Blades'. Computational Fluid Dynamics simulations were used to studyincidence angles from -37 to +10 degrees on the aerodynamic performance of the Seal Blade. The tests and simulationswere conducted at a Reynolds number of 100,000. The Seal Blades showed consistent performance improvements overthe baseline configuration. It was determined that a fuel burn reduction of approximately 5 could be achieved for a fixedwing aircraft. Noise reduction potential is also explored

vortex induced vibration↗
NASA NTRS2013

Conductance of Ion Channels - Theory vs. Experiment

Transmembrane ion channels mediate a number of essential physiological processes in a cell ranging from regulating osmotic pressure to transmission of neural signals. Kinetics and selectivity of ion transport is of critical importance to a cell and, not surprisingly, it is a subject of numerous experimental and theoretical studies. In this presentation we will analyze in detail computer simulations of two simple channels from fungi - antiamoebin and trichotoxin. Each of these channels is made of an alpha-helical bundle of small, nongenomically synthesized peptides containing a number of rare amino acids and exhibits strong antimicrobial activity. We will focus on calculating ionic conductance defined as the ratio of ionic current through the channel to applied voltage. From molecular dynamics simulations, conductance can be calculated in at least two ways, each involving different approximations. Specifically, the current, given as the number of charges transferred through the channel per unit of time, can be obtained from the number of events in which ions cross the channel during the simulation. This method works well for large currents (high conductance values and/or applied voltages). If the number of crossing events is small, reliable estimates of current are difficult to achieve. Alternatively, conductance can be estimated assuming that ion transport can be well approximated as diffusion in the external potential given by the free energy profile. Then, the current can be calculated by solving the one-dimensional diffusion equation in this external potential and applied voltage (the generalized Nernst-Planck equation). To do so three ingredients are needed: the free energy profile, the position-dependent diffusion coefficient and the diffusive flux of ions into the channel. All these quantities can be obtained from molecular dynamics simulations. An important advantage of this method is that it can be used equally well to estimating large and small currents. In addition, once the free energy profile becomes available the full current-voltage dependence can be readily obtained. For both channels we carried out calculations using both approaches. We also tested the main assumptions underlying the diffusive model, such as uncorrelated nature of individual crossing events and Fickian diffusion. The accuracy and consistency of different methods will be discussed. Finally we will discuss how comparisons between calculated and measured ionic conductance and selectivity of transport can be used for determining structural models of the channels.

Pohorille, Andrew↗
DOE OSTI2018

Developing an interatomic potential for martensitic phase transformations in zirconium by machine learning

Abstract Atomic simulations provide an effective means to understand the underlying physics of structural phase transformations. However, this remains a challenge for certain allotropic metals due to the failure of classical interatomic potentials to represent the multitude of bonding. Based on machine-learning (ML) techniques, we develop a hybrid method in which interatomic potentials describing martensitic transformations can be learned with a high degree of fidelity from ab initio molecular dynamics simulations (AIMD). Using zirconium as a model system, for which an adequate semiempirical potential describing the phase transformation process is lacking, we demonstrate the feasibility and effectiveness of our approach. Specifically, the ML-AIMD interatomic potential correctly captures the energetics and structural transformation properties of zirconium as compared to experimental and density-functional data for phonons, elastic constants, as well as stacking fault energies. Molecular dynamics simulations successfully reproduce the transformation mechanisms and reasonably map out the pressure–temperature phase diagram of zirconium.

Zong, Hongxiang↗
DOE OSTI2018

Microscopic Behaviors of Tri-n-Butyl Phosphate, n-Dodecane and Their Mixtures at Air/Liquid and Liquid/Liquid Interfaces: An AMBER Polarizable Force Field Study

In solvent extraction processes for recovering metal ions from used nuclear fuel, as well as other industrial applications, a better understanding of the metal complex phase transfer phenomenon would greatly aid ligand design and process optimization. We have approached this challenge by utilizing classical molecular dynamics simulations technique to gain visual appreciation of the vapor/liquid and liquid/liquid interface between tri-n-butyl phosphate (TBP) and n-dodecane with air and water. In this study, we successfully reparameterized polarizable force fields for TBP and n-dodecane that accurately reproduced several of their thermophysical properties such as: density, heat of vaporization, and dipole moment. Our models were able to predict the surface and interfacial tension of different systems when compared to experimental results also performed by us. Through this study, we gained numerous atomistic understanding of the behaviors of TBP and n-dodecane at the interface against air and water, useful in further computational studies of such systems. Lastly, our studies indicate that the initial configuration of a simulation may have a large effect on the final result.

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

Application of Pinniped Vibrissae to Aeropropulsion

Vibrissae of Phoca Vitulina (Harbor Seal) and Mirounga Angustirostris (Elephant Seal) possess undulations along their length. Harbor Seal Vibrissae were shown to reduce vortex induced vibrations and reduce drag compared to appropriately scaled cylinders and ellipses. Samples of Harbor Seal vibrissae, Elephant Seal vibrissae and California Sea Lion vibrissae were collected from the Marine Mammal Center in California. CT scanning, microscopy and 3D scanning techniques were utilized to characterize the whiskers. Computational fluid dynamics simulations of the whiskers were carried out to compare them to an ellipse and a cylinder. Leading edge parameters from the whiskers were used to create a 3D profile based on a modern power turbine blade. The NASA SW-2 facility was used to perform wind tunnel cascade testing on the 'Seal Blades'. Computational Fluid Dynamics simulations were used to study incidence angles from -37 to +10 degrees on the aerodynamic performance of the Seal Blade. The tests and simulations were conducted at a Reynolds number of 100,000. The Seal Blades showed consistent performance improvements over the baseline configuration. It was determined that a fuel burn reduction of approximately 5 could be achieved for a fixed wing aircraft. Noise reduction potential is also explored.

wakes↗
DOE OSTI2016

Local structure, composition, and crystallization mechanism of a model two-phase “composite nanoglass”

We report a detailed study of the local composition and structure of a model, bi-phasic nanoglass with nominal stoichiometry Cu{sub 55}Nb{sub 45}. Three dimensional atom probe data suggest a nanoscale-phase-separated glassy structure having well defined Cu-rich and Nb-rich regions with a characteristic length scale of ≈3 nm. However, extended x-ray absorption fine structure analysis indicates subtle differences in the local environments of Cu and Nb. While the Cu atoms displayed a strong tendency to cluster and negligible structural order beyond the first coordination shell, the Nb atoms had a larger fraction of unlike neighbors (higher chemical order) and a distinctly better-ordered structural environment (higher topological order). This provides the first experimental indication that metallic glass formation may occur due to frustration arising from the competition between chemical ordering and clustering. These observations are complemented by classical as well as ab initio molecular dynamics simulations. Our study indicates that these nanoscale phase-separated glasses are quite distinct from the single phase nanoglasses (studied by Gleiter and others) in the following three respects: (i) they contain at least two structurally and compositionally distinct, nanodispersed, glassy phases, (ii) these phases are separated by comparatively sharp inter-phase boundaries, and (iii) thermally induced crystallization occurs via a complex, multi-step mechanism. Such materials, therefore, appear to constitute a new class of disordered systems that may be called a composite nanoglass.

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

The thermodynamic and kinetic interactions of He interstitial clusters with bubbles in W

Due to its enviable properties, tungsten is a leading candidate plasma facing material in nuclear fusion reactors. However, like many other metals, tungsten is known to be affected by the high doses of helium atoms incoming from the plasma. Indeed, the implanted interstitial helium atoms cluster together and, upon reaching a critical cluster size, convert into substitutional nanoscale He bubbles. These bubbles then grow by absorbing further interstitial clusters from the matrix. This process can lead to deleterious changes in microstructure, degradation of mechanical properties, and contamination of the plasma. In order to better understand the growth process, we use traditional and accelerated molecular dynamics simulations to investigate the interactions between interstitial He clusters and pre-existing bubbles. These interactions are characterized in terms of thermodynamics and kinetics. We show that the proximity of the bubble leads to an enhancement of the trap mutation rate and, consequently, to the nucleation of satellite bubbles in the neighborhood of existing ones. We also uncover a number of mechanisms that can lead to the subsequent annihilation of such satellite nanobubbles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

NSTXU Diagnostic Disruption Dynamic Loading Represented by Response Spectra

This article presents the results of transient dynamic simulations of loads due to disruption eddy currents on the NSTXU vacuum vessel. Dynamic loading at diagnostic mounting locations is expressed as response spectra derived from the time history results of the dynamic structural simulations of a variety of disruption scenarios. The disruption simulations draw on a history of the project assessments of worst case disruptions for specific components. Major efforts to assess disruption loading have included the vacuum vessel which is the major structural support for the machine, as well as the passive plates (PPs), high harmonic fast wave (HHFW) antenna, and centerstack casing. Each one of these efforts included transient electromagnetic simulations producing time-dependent eddy current Lorentz loads (and in some cases halo loads) which then were applied to time-dependent structural dynamic analyses intended to obtain the proper dynamic amplification factors. In some instances, the EM model and structural model were identical allowing direct transfer of EM forces to the structural model. In other cases, the EM and structural model were not identical and the vector potential (VP) transfer method was used. The results files from these analyses were available (or re-run) to post process in ANSYS Classic time history postprocessor. In conclusion, the ANSYS command is used to create response spectra from time history data at desired points on the vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Communication: Dynamical and structural analyses of solid hydrogen under vapor pressure

Nuclear quantum effects play a dominant role in determining the phase diagram of H{sub 2}. With a recently developed quantum molecular dynamics simulation method, we examine dynamical and structural characters of solid H{sub 2} under vapor pressure, demonstrating the difference from liquid and high-pressure solid H{sub 2}. While stable hexagonal close-packed lattice structures are reproduced with reasonable lattice phonon frequencies, the most stable adjacent configuration exhibits a zigzag structure, in contrast with the T-shape liquid configuration. The periodic angular distributions of H{sub 2} molecules indicate that molecules are not a completely free rotor in the vapor-pressure solid reflecting asymmetric potentials from surrounding molecules on adjacent lattice sites. Discrete jumps of librational and H–H vibrational frequencies as well as H–H bond length caused by structural rearrangements under vapor pressure effectively discriminate the liquid and solid phases. The obtained dynamical and structural information of the vapor-pressure H{sub 2} solid will be useful in monitoring thermodynamic states of condensed hydrogens.

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

The hTERT core promoter forms three parallel G-quadruplexes

The structure of the 68 nt sequence with G-quadruplex forming potential within the hTERT promoter is disputed. One model features a structure with three stacked parallel G-quadruplex units, while another features an unusual duplex hairpin structure adjoined to two stacked parallel and antiparallel quadruplexes. We report here the results of an integrated structural biology study designed to distinguish between these possibilities. As part of our study, we designed a sequence with an optimized hairpin structure and show that its biophysical and biochemical properties are inconsistent with the structure formed by the hTERT wild-type sequence. By using circular dichroism, thermal denaturation, nuclear magnetic resonance spectroscopy, analytical ultracentrifugation, small-angle X-ray scattering, molecular dynamics simulations and a DNase I cleavage assay we found that the wild type hTERT core promoter folds into a stacked, three-parallel G-quadruplex structure. The hairpin structure is inconsistent with all of our experimental data obtained with the wild-type sequence. All-atom models for both structures were constructed using molecular dynamics simulations. These models accurately predicted the experimental hydrodynamic properties measured for each structure. We found with certainty that the wild-type hTERT promoter sequence does not form a hairpin structure in solution, but rather folds into a compact stacked three-G-quadruplex conformation.

59 BASIC BIOLOGICAL SCIENCES↗