Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Thermal exchange-correlation functionals: Capturing quantum electron behavior in warm, dense plasmas

We summarize and give perspective upon recent progress in developing non-empirical constraint-based thermal (i.e., free energy) exchange-correlation (XC) density functionals essential for accurate description of the quantum behavior of electrons in warm, dense plasmas. After delineating the critical role of ground-state functionals for zero-temperature, time-dependent DFT, we outline the underpinnings of local density approximation, generalized gradient approximation (GGA), and meta-GGA XC free-energy functionals. Two basic thermalization principles for upgrading ground-state XC functionals to successful thermal ones are emphasized. Then, we turn to a long-standing challenge, assessment of the accuracy of well-founded functionals. Unlike the ground state, there are a few exact results for large T and P. An exception is path integral Monte Carlo (PIMC) data for dense H/D and He plasmas. For those, we did ab initio molecular dynamics simulations under selected thermodynamic conditions employing five thermal XC functionals: two approximate thermal GGAs, fully thermal GGA, an approximate meta-GGA, and fully thermal meta-GGA. Comparisons with the PIMC data show that functionals thermalized by augmenting a non-thermal functional with a lower-level thermal contribution are inferior to functionals with thermal XC and spatial inhomogeneity effects taken into account at the same level of refinement. We believe this and similar evidence should be convincing to the high-energy density physics community of the necessity of use of proper thermal XC functionals in simulation studies of finite-temperature quantum effects in warm, dense plasmas.

Ab-initio molecular dynamics↗
DOE OSTI2025

Impact of Interfacial Structure on Heterogeneous Nucleation of Amorphous Carbonates

For this work, classical molecular dynamics simulations were performed to provide physical insight into the impact of interfacial structure on the heterogeneous nucleation of amorphous calcium carbonate (ACC, CaCO 3 ·H 2 O) and amorphous magnesium carbonate (AMC, MgCO 3 ·H 2 O) by using α-quartz as a model substrate. Interfacial structure and energies were computed for ACC and AMC in contact with the (100), (001), and (101) α-quartz surfaces. The simulations showed α-quartz surfaces drew water molecules out of the carbonate nuclei to form a partial hydration layer. The formation of a partial hydration layer and its disruption to the ACC/AMC structure meant the α-quartz–ACC/AMC interfaces were not energetically favored relative to separate α-quartz–water and ACC/AMC–water interfaces and, thus, homogeneous ACC/AMC nucleation was favored over heterogeneous nucleation. The CMD simulations hence provided an atomic-level explanation for a reported nonclassical growth mechanism whereby carbonate minerals grow via homogeneous nucleation and subsequent surface attachment of amorphous intermediates.

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

Unexpected inverse correlations and cooperativity in ion-pair phase transfer

Liquid/liquid extraction is one of the most widely used separation and purification methods, where a forefront of research is the study of transport mechanisms for solute partitioning and the relationships that these have to solution structure at the phase boundary. To date, organized surface features that include protrusions, water-fingers, and molecular hinges have been reported. Many of these equilibrium studies have focused upon small-molecule transport – yet the extent to which the complexity of the solute, and the competition between different solutes, influence transport mechanisms have not been explored. Here we report molecular dynamics simulations that demonstrate that a metal salt (LiNO 3 ) can be transported via a protrusion mechanism that is remarkably similar to that reported for H 2 O by tri-butyl phosphate (TBP), a process that involves dimeric assemblies. Yet the LiNO 3 out-competes H 2 O for a bridging position between the extracting TBP dimer, which in-turn changes the preferred transport pathway of H 2 O. Examining the electrolyte concentration dependence on ion-pair transport unexpectedly reveals an inverse correlation with the extracting surfactant concentration. As [LiNO 3 ] increases, surface adsorbed TBP becomes a limiting reactant in correlation with an increased negative surface charge induced by excess interfacial NO 3 – , however the rate of transport is enhanced. Within the highly dynamic interfacial environment, we hypothesize that this unique cooperative effect may be due to perturbed surface organization that either decreases the energy of formation of transporting protrusion motifs or makes it easier for these self-assembled species to disengage from the surface.

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

Data-Driven Approach for Controlled Icosahedral Boron- Rich Compound Growth

This final technical report summarizes the research accomplishments and research highlights at the end of the funding period. This project aimed to leverage existing and new computational data produced from first-principles and molecular dynamics simulations to understand the thermodynamic, mechanical, and electronic properties of icosahedral boron compounds. The goal is to achieve targeted material properties by controlling the synthesis routes of these boron-rich compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Colloidal Clusters and Networks Formed by Oppositely Charged Nanoparticles with Varying Stiffnesses

Colloidal clusters and gels are ubiquitous in science and technology. Particle softness has a strong effect on interparticle interactions; however, our understanding of the role of this factor in the formation of colloidal clusters and gels is only beginning to evolve. Here, we report the results of experimental and simulation studies of the impact of particle softness on the assembly of clusters and networks from mixtures of oppositely charged polymer nanoparticles (NPs). Experiments were performed below or above the polymer glass transition temperature, at which the interaction potential and adhesive forces between the NPs were significantly varied. Hard NPs assembled in fractal clusters that subsequently organized in a kinetically arrested colloidal gel, while soft NPs formed dense precipitating aggregates, due to the NP deformation and the decreased interparticle distance. Importantly, interactions of hard and soft NPs led to the formation of discrete precipitating NP aggregates at a relatively low volume fraction of soft NPs. A phenomenological model was developed for interactions of oppositely charged NPs with varying softnesses. Furthermore, the experimental results were in agreement with molecular dynamics simulations based on the model. This work provides insight on interparticle interactions before, during, and after the formation of hard–hard, hard–soft, and soft–soft contacts and has impact for numerous applications of reversible colloidal gels, including their use as inks for additive manufacturing.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Rate of Para-Orthohydrogen Conversion in Cryogenic Vortex Tube

Raising efficiency of the cooling process for cryogenic hydrogen and minimizing hydrogen boil-off during storage and transportation of liquid hydrogen are critically important factors for widening implementation of hydrogen-based clean energy systems. One novel approach involves utilization of relatively simple catalyzed vortex tubes, where hydrogen in the primarily para- nuclear spin isomer form is converted into the ortho- form while consuming a significant amount of heat. To design such systems effectively, better understanding is required about catalyst-assisted para-orthohydrogen conversion rates in high-speed vortical flows. In this study, a computational fluid dynamics simulation has been set up to model an experimental system with cryogenic hydrogen vortex tubes. Mesh-verification and a validation study has been conducted first for a non-catalyzed tube. Then, the rate coefficient of the para-orthohydrogen conversion of cryogenic hydrogen has been determined by matching numerical results with experimental data available for a catalyzed vortex tube. Furthermore, the presented information can help design and optimize cryogenic hydrogen cooling within vortex tubes.

08 HYDROGEN↗
DOE OSTI2021

The RIT1 C-terminus associates with lipid bilayers via charge complementarity

RIT1 is a member of the Ras superfamily of small GTPases involved in regulation of cellular signaling. Mutations to RIT1 are involved in cancer and developmental disorders. Like many Ras subfamily members, RIT1 is localized to the plasma membrane. However, RIT1 lacks the C-terminal prenylation that helps many other subfamily members adhere to cellular membranes. For this study, we used molecular dynamics simulations to examine the mechanisms by which the C-terminal peptide (CTP) of RIT1 associates with lipid bilayers. We show that the CTP is unstructured and that its membrane interactions depend on lipid composition. While a 12-residue region of the CTP binds strongly to anionic bilayers containing phosphatidylserine lipids, the CTP termini fray from the membrane allowing for accommodation of the RIT1 globular domain at the membrane-water interface.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Radiative and Nonradiative Recombinations in Organic Radical Emitters: The Effect of Guest–Host Interactions

Radical-carrying organic molecules have received significant attention to bypass the issue related to harvesting triplet excitons in current light-emitting materials. While the computational efforts conducted so far have treated these radical emitters as isolated entities, in actual devices, they are embedded in a host matrix and subject to emitter–host interactions. In this study, by combining molecular dynamics simulations and density functional theory calculations, the impact of the host matrix on the optoelectronic performance of radical emitters is evaluated, taking as a representative example the (4-ncarbazolyl-2,6-dichlorophenyl)bis(2,4,6-trichlorophenyl)-methyl (TTM-3NCz) radical emitter dispersed in a 4,4-bis(carbazol-9-yl)biphenyl (CBP) host. A morphological analysis shows that steric effects around the radical centers, carried by the TTM electron-poor moieties of the emitters, disfavor π–π interactions with the host molecules, which leads to random intermolecular orientations around the TTM moieties. The 3NCz electron-rich moieties of the emitters, however, have much lesser spatial hindrance for intermolecular π–π stacking, which modulates the structural and electronic properties of the emitters in the host matrix. The influence of dynamic and static disorders on the radiative and nonradiative recombination processes is also investigated and it is found that the rates of nonradiative recombination are small, which opens the way to 100% internal quantum efficiency for the doublet-based emission process.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nucleation Rate Theory for Coordination Number: Elucidating Water-Mediated Formation of a Zigzag Na 2 SO 4 Morphology

Predicting and controlling nanostructure formation during nucleation can pave the way to synthesizing novel energy materials via crystallization. However, such control over nucleation and crystallization remains challenging due to an inadequate understanding of critical factors that govern evolving atomistic structures and dynamics. Here, in this paper, we utilize coordination number as a reaction coordinate and rate theory to investigate how sodium sulfate, commonly known as a phase-change energy material, nucleates in a supersaturated aqueous solution. In conjunction with ab initio and force field-based molecular dynamics simulation, the rate theoretical analysis reveals that sodium sulfate from an initially dissolved metastable state transits to a heterogeneous mixture of prenucleated clusters and finally to a large cylindrical zigzag morphology. Measurements of Raman spectra and their ab initio modeling confirm that this nucleated morphology contains a few waters for every sulfate. Rate processes such as solvent exchange and desolvation exhibit high sensitivity to the evolving prenucleation/nucleation structures, providing a means to distinguish between critical nucleation precursors. Desolvation and forming the first-shell interionic coordination structure via monomer-by-monomer addition around sulfates are found to explain the formation of large nuclei. Thus, a detailed understanding of the step-by-step structure formation across scales has been achieved. This can be leveraged to predict nucleation-related structures and dynamics and potentially control the synthesis of novel phase-change materials for energy applications.

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

Thermally frustrated phase transition at high pressure in B 2-ordered FeV

X-ray diffraction measurements of equiatomic B2-ordered FeV were performed in a diamond-anvil cell at room temperature at several pressure points up to 80 GPa that showed the cubic phase to be stable with no indication of structural phase transitions. Density functional theory at 0 K predicts Fermi surface nesting, an electronic topological transition, and a phonon dynamical instability within the experimentally investigated pressure range. Nevertheless, the instability is absent in phonon dispersion curves extracted from ab initio molecular dynamics simulations below the critical volume at temperatures as low as 10 K, indicating that thermal atomic displacements can frustrate the phase transition by renormalizing the phonon dispersion curves. Ferrimagnetism is critical for the stability of the cubic phase at low temperature, but thermal atomic displacements are enough to support the structure at and above the Néel temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Insights on the dissolution of water in an albite melt at high pressures and temperatures from a direct structural analysis

The water dissolution mechanism in silicate melts under high pressures is not well understood. Here we present the first direct structure investigation of a water-saturated albite melt to monitor the interactions between water and the network structure of silicate melt at the molecular level. In situ high-energy X-ray diffraction was carried out on the NaAlSi 3 O 8 -H 2 O system at 800 °C and 300 MPa, at the Advanced Photon Source synchrotron facility. The analysis of the X-ray diffraction data was augmented with classical Molecular Dynamics simulations of a hydrous albite melt, incorporating accurate water-based interactions. The results show that metal–oxygen bond breaking at the bridging sites occurs overwhelmingly at the Si site upon reaction with H 2 O, with subsequent Si–OH bond formation and negligible Al–OH formation. Furthermore, we see no evidence for the dissociation of the Al 3+ ion from the network structure upon breaking of the Si–O bond in the hydrous albite melt. The results also indicate that the Na + ion is an active participant in the modifications of the silicate network structure of the albite melt upon water dissolution at high P–T conditions. We do not find evidence for the Na + ion dissociating from the network structure upon depolymerization and subsequent formation of NaOH complexes. Instead, our results show that the Na + ion persists as a structure modifier with a shift away from Na–BO bonding to an increase in the extent of Na-NBO bonding, in parallel with pronounced depolymerization of the network. Our MD simulations show that the Si–O and Al–O bond lengths are expanded by about 6% in the hydrous albite melt compared to those of the dry melt at high P–T conditions. The changes in the network silicate structure of a hydrous albite melt at high pressure and temperature, as revealed in this study, must be considered in the advancement of water dissolution models of hydrous granitic (or alkali aluminosilicate) melts.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nucleation of Dislocations in 3.9 nm Nanocrystals at High Pressure

As circuitry approaches single nanometer length scales, it has become important to predict the stability of single nanometer-sized metals. The behavior of metals at larger scales can be predicted based on the behavior of dislocations, but it is unclear if dislocations can form and be sustained at single nanometer dimensions. Here, we report the formation of dislocations within individual 3.9 nm Au nanocrystals under nonhydrostatic pressure in a diamond anvil cell. We used a combination of x-ray diffraction, optical absorbance spectroscopy, and molecular dynamics simulation to characterize the defects that are formed, which were found to be surface-nucleated partial dislocations. These results indicate that dislocations are still active at single nanometer length scales and can lead to permanent plasticity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Development of a Computational Framework for Multiphysics Multiphase Species Tracking using NEAMS Tools

This report implements a high-fidelity multiphysics modeling framework using the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program tools to track isotopic species in Molten Salt Reactors (MSRs), with a specific focus on the 91-depletion chain within the Molten Salt Reactor Experiment (MSRE). The model integrates neutronics, thermal-hydraulics, depletion, and thermochemistry to simulate the production, transport, and phase transitions of isotopes under steady-state and transient conditions. The main findings reveal that isotopes such as bromine-91 largely remain in the liquid phase, while others, including krypton-91and yttrium-91, transition to the gas phase, significantly influencing the reactor’s radiological source term. The study also shows that during transients, like a reactivity insertion transient, rapid void formation and the expansion of the liquid-gas interface led to substantial transfers of dissolved isotopes into the gas phase, altering isotope distribution and largely increasing the source term in the off-gas system. Additionally, the research highlights that short-lived isotopes dominate the initial off-gas response during transients, while longer-lived isotopes determine the equilibrium state, underscoring the necessity of dynamic simulations for accurate species tracking and reactor safety analysis. The developed methodology will be applied in the future to the tracking of a larger number of species and introduce other species tracking mechanisms, such as deposition and plating.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Proton Collective Quantum Tunneling Induces Anomalous Thermal Conductivity of Ice under Pressure

Proton tunneling is believed to be nonlocal in ice, but its range has been shown to be limited to only a few molecules. Here, we measured the thermal conductivity of ice under pressure up to 50 GPa and found it increases with pressure until 20 GPa but decreases at higher pressures. We attribute this nonmonotonic thermal conductivity to the collective tunneling of protons at high pressures, supported by large-scale quantum molecular dynamics simulations. The collective tunneling loops span several picoseconds in time and are as large as nanometers in space, which match the phonon periods and wavelengths, leading to strong phonon scattering at high pressures. Finally, our results show direct evidence of global quantum motion existing in high-pressure ice and provide a new perspective to understanding the coupling between phonon propagation and atomic tunneling.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Three-state harmonic models for photoinduced charge transfer

A widely used strategy for simulating the charge transfer between donor and acceptor electronic states in an all-atom anharmonic condensedphase system is based on invoking linear response theory to describe the system in terms of an effective spin-boson model Hamiltonian. Extending this strategy to photoinduced charge transfer processes requires also taking into consideration the ground electronic state in addition to the excited donor and acceptor electronic states. In this paper, we revisit the problem of describing such nonequilibrium processes in terms of an effective three-state harmonic model. We do so within the framework of nonequilibrium Fermi’s golden rule (NE-FGR) in the context of photoinduced charge transfer in the carotenoid–porphyrin–C 60 (CPC 60 ) molecular triad dissolved in explicit tetrahydrofuran (THF). To this end, we consider different ways for obtaining a three-state harmonic model from the equilibrium autocorrelation functions of the donor–acceptor, donor–ground, and acceptor–ground energy gaps, as obtained from all-atom molecular dynamics simulations of the CPC 60 /THF system. The quantum-mechanically exact time-dependent NE-FGR rate coefficients for two different charge transfer processes in two different triad conformations are then calculated using the effective three-state model Hamiltonians as well as a hierarchy of more approximate expressions that lead to the instantaneous Marcus theory limit. Our results show that the photoinduced charge transfer in CPC 60 /THF can be described accurately by the effective harmonic three-state models and that nuclear quantum effects are small in this system.

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

Regimes of evaporation and mixing behaviors of nanodroplets at transcritical conditions

The objective of this paper is to examine the fundamental mechanisms responsible for the transition between subcritical evaporation and supercritical dense-fluid-mixing in the absence of convection effects, specifically focusing on the liquid–vapor interfacial dynamics. To isolate the dynamics of this transition process, we characterize the different physical behaviors exhibited by an -dodecane nanoscale droplet placed in different nitrogen ambient conditions across the fuel’s critical point. We employ a continuum-based interface-resolving diffuse-interface method to explore the underlying phase-exchange mechanisms that bring about such distinct dynamics. Following the comparison against molecular dynamics simulations and experiments of evaporating droplets and experimental data for vapor–liquid equilibria, a parametric study at various ambient conditions and droplet sizing is performed to identify four regimes of evaporation/mixing behaviors: sub- and supercritical droplet evaporation, and sub- and supercritical dense-fluid-mixing. It is shown that the distinction in the phase-exchange mechanisms in these four regimes are brought about by the different thermodynamic phases the droplet center can exhibit during the evaporation/mixing process: subcritical liquid, supercritical liquid-like, subcritical gaseous, and supercritical gas-like, respectively. It is shown that the subcritical dense-fluid-mixing behavior is a direct result of nanoconfinement of the liquid–vapor interfacial structure and thus is not present for large droplet sizes. Finally, the present study also shows that the supercritical phase-exchange dynamics can follow two different pathways: supercritical droplet-like evaporation and supercritical dense-fluid-mixing. Furthermore, promoting the early transition to supercritical dense-fluid-mixing can significantly expedite the phase-exchange process through the disintegration of the liquid-like droplet core.

42 ENGINEERING↗
DOE OSTI2020

Pore-scale simulation of gas flow in microscopic permeable media with complex geometries

Lattice Boltzmann method (LBM) is an efficient tool to perform direct numerical simulation of gas flow in micro-size pores. Despite great advances in LBM and its use in describing gas flow in micro-size channels in slip and transitional flow regimes, exploring formulations that are able to capture flow behavior in domains with complex boundary geometries remains a challenging task. As a result, the impact of complexities in pore structures on gas flow may not be fully explored using LBM. In this study, we propose to use the combined bounce-back and Maxwellian diffusive reflection scheme to capture the slip velocity with a new set of slip coefficients. Optimal values of slip coefficients are determined using a gradient based method, where the results from the linearized Boltzmann equation is used as the reference solution. The proposed formulation is further validated against predictions from molecular dynamics simulation in the presence of complex geometries by introducing obstacles of different shapes in straight channels. LBM with the proposed boundary treatment is then used to investigate gas flow in a synthetic microscopic permeable medium. Here, results indicate that gas exhibits different flow configurations as Knudsen number varies, and apparent gas permeability appears to have up to a second-order dependency on the reciprocal mean pressure in slip and transitional flow regimes.

03 NATURAL GAS↗