Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer simulations”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Computer simulation of reaction and transport of core-scale serpentinization of Fe-bearing olivine and evolution of geological hydrogen

We present reactive flow and transport simulations of core-scale serpentinization of forsteritic, Fe-bearing olivine under hydrothermal conditions. The model captures fluid flow, evolving porosity–permeability, and redox-controlled H 2 production in a cylindrical core over 55 days at 245 °C and 37 bar. Key processes include olivine dissolution, precipitation of serpentine and magnetite, and oxidation of Fe2+ to Fe 3+ . Here, results show strong coupling among flow velocity, alteration front propagation, and spatially heterogeneous H 2 generation. Elevated H2 concentrations align with inlet-localized magnetite precipitation, consistent with experiments. Secondary mineral formation reduces porosity and permeability, altering transport pathways. Reactive flow at the fluid–solid interface governs H 2 generation rates and distribution. Dissolved SiO2 promotes serpentine and talc formation by suppressing brucite, while bicarbonate extends reaction duration and moderates surface complexation effects.

Geologic hydrogen↗
DOE OSTI2018

Computational simulations of [email protected] and [email protected] [email protected] nanotubes

Highlights: • [email protected] nanotubes are easily formed than [email protected] nanotubes. • Armchair [email protected] nanotubes present type-II behavior, while zigzag [email protected] nanotubes present type-I behavior. • The [email protected] had the piezoelectricity response improved until 10% of elongation. - Abstract: The structural, electronic, and mechanical properties of armchair and zigzag chiralities double-walled and [email protected] [email protected] and [email protected] nanotubes were investigated by periodic DFT/B3LYP calculations with an all-electron basis set. For both chiralities, [email protected] presents minor strain and deposition energies, which predict that this nanotube can be easier formed and the GaN is the most favorable substrate (the core) than ZnO. On the other hand, the zigzag [email protected] did not exhibit the major piezoelectric response, which is three times smaller than the [email protected] nanotube, showing that the compression of the [email protected] nanotube length is not favorable to this property. However, the piezoelectricity can be improved when the zigzag [email protected] is under elongation. The elastic constants showed that the [email protected] nanotubes are more rigid than the homogenous nanotubes and present higher piezoelectric constants. In addition, the projected DOS shows that [email protected] has a type-II interface and [email protected] has a type-I interface. Based on the results obtained from our theoretical models, the nanotubes have great potential for the experimental development of new electronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Hydrogen Dynamics in Supercritical Water Probed by Neutron Scattering and Computer Simulations

In this work, an investigation of supercritical water is presented combining inelastic and deep inelastic neutron scattering experiments and molecular dynamics simulations based on a machine-learned potential of ab initio quality. The local hydrogen dynamics is investigated at 250 bar and in the temperature range of 553–823 K, covering the evolution from subcritical liquid to supercritical gas-like water. The evolution of libration, bending, and stretching motions in the vibrational density of states is studied, analyzing the spectral features by a mode decomposition. Moreover, the hydrogen nuclear momentum distribution is measured, and its anisotropy is probed experimentally. It is shown that hydrogen bonds survive up to the higher temperatures investigated, and we discuss our results in the framework of the coupling between intramolecular modes and intermolecular librations. Results show that the local potential affecting hydrogen becomes less anisotropic within the molecular plane in the supercritical phase, and we attribute this result to the presence of more distorted hydrogen bonds.

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

High temperature breakdown of the Stokes-Einstein relation in a computer simulated Cu-Zr melt

Transport properties and the Stokes-Einstein (SE) relation in liquid Cu{sub 8}Zr{sub 3} are studied by molecular dynamics simulation with a modified embedded atom potential. The critical temperature T{sub c} of mode coupling theory (MCT) is derived as 930 K from the self-diffusion coefficient D and viscosity η. The SE relation breaks down around T{sub SE} = 1900 K, which is far above T{sub c}. At temperatures below T{sub SE}, the product of D and η fluctuates around a constant value, similar to the prediction of MCT near T{sub c}. The influence of the microscopic atomic motion on macroscopic properties is investigated by analyzing the time dependent liquid structure and the self-hole filling process. The self-holes for the two components are preferentially filled by atoms of the same component. The self-hole filling dynamics explains the different breakdown behaviors of the SE relation in Zr-rich liquid CuZr{sub 2} compared to Cu-rich Cu{sub 8}Zr{sub 3}. At T{sub SE}, a kink is found in the temperature dependence of both partial and total coordination numbers for the three atomic pair combinations and of the typical time of self-hole filling. This indicates a strong correlation between liquid structure, atomic dynamics, and the breakdown of SE relation. The previously suggested usefulness of the parameter d(D{sub 1}/D{sub 2})/dT to predict T{sub SE} is confirmed. Additionally we propose a viscosity criterion to predict T{sub SE} in the absence of diffusion data.

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

Evaluation and optimization of the performance of frame geometries for lithium-ion battery application by computer simulation

Tailoring battery geometries is essential for many applications, as geometry influences the delivered capacity value. Two geometries, frame and conventional, have been studied and, for a given scan rate of 330C, the square frame shows a capacity value of 305,52 Ahm{sup −2}, which is 527 times higher than the one for the conventional geometry for a constant the area of all components.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Computational simulations of solvation force and squeezing out of dodecane chain molecules in an atomic force microscope

Understanding the squeeze out behaviors of liquid films at nanometer scale in an atomic force microscope (AFM) has been a significant interest since the 1990s. We carry out all-atom static-mode AFM simulations in a liquid-vapor molecular dynamics ensemble to investigate the solvation force oscillation and squeeze out mechanisms of a confined linear dodecane fluid between a gold AFM tip and a mica substrate. Solvation force oscillations are found to be associated with the layering transition of the liquid film and unstable jumps of the AFM tip. Detailed structural analyses and molecular animations show that the local permeation of chain molecules and the squeeze out of molecules near the edge of contact promote the layering transition under compression. The confinement-induced slow down dynamics is manifested by the decrease in diffusivity and increase in rotational relaxation times. Furthermore, the persistent diffusive behavior of dodecane chain molecules even in the single-monolayer film is attributed to the chain sliding motions in the film due to the substantial vacancy space and thermal fluctuations.

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

The interplay between chromosome stability and cell cycle control explored through gene–gene interaction and computational simulation

Chromosome stability models are usually qualitative models derived from molecular-genetic mechanisms for DNA repair, DNA synthesis, and cell division. While qualitative models are informative, they are also challenging to reformulate as precise quantitative models. In this report we explore how (A) laboratory experiments, (B) quantitative simulation, and (C) seriation algorithms can inform models of chromosome stability. Laboratory experiments were used to identify 19 genes that when over-expressed cause chromosome instability in the yeast Saccharomyces cerevisiae . To better understand the molecular mechanisms by which these genes act, we explored their genetic interactions with 18 deletion mutations known to cause chromosome instability. Quantitative simulations based on a mathematical model of the cell cycle were used to predict the consequences of several genetic interactions. These simulations lead us to suspect that the chromosome instability genes cause cell-cycle perturbations. Cell-cycle involvement was confirmed using a seriation algorithm, which was used to analyze the genetic interaction matrix to reveal an underlying cyclical pattern. The seriation algorithm searched over 10 14 possible arrangements of rows and columns to find one optimal arrangement, which correctly reflects events during cell cycle phases. To conclude, we illustrate how the molecular mechanisms behind these cell cycle events are consistent with established molecular interaction maps.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Phonon second harmonic generation in NaBr studied by inelastic neutron scattering and computer simulation

The phenomenon of second harmonic generation (SHG) was found for phonons in anharmonic NaBr by inelastic neutron scattering. The temperature dependence of this phonon SHG was measured from 300 K to 650 K. At 300 K the second harmonic (SH) is seen as a high-energy branch around 33 meV, nearly independent of $\overrightarrow{Q}$. The temperature effective potential (TDEP) method and classical molecular dynamics (MD) simulation with machine learning interatomic potential were able to reproduce the SH, and showed that SHG occurs with the flat transverse optical (TO) phonon branch. A classical model of a nonlinear medium explains the intensity and lifetime of the SH, compared to those of the TO modes. Also successful was a quantum model based on the Heisenberg-Langevin equation for interacting phonons coupled to a thermal bath, which also predicts a spectral distribution of the SH. In conclusion, the measured temperature dependence of the intensity of the second harmonic showed that it follows the Planck distribution of a one-phonon quasiparticle, and not two TO phonons.

36 MATERIALS SCIENCE↗