Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solid dynamics”

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

A stress-based fracture model for reacting metal ejecta

The evolution of reacting metal ejecta continues to be a topic of interest at the forefront of metals in reactive and extreme environments. Ejecta are small particles formed when the surface of a metal undergoes Richtmyer–Meshkov instability from a strong shock. Experiments have shown that in the case where ejecta are in ambient conditions that induce a reaction, the ejecta behave irregularly. The ejecta temperature rises and then plateaus, and the acceleration profile shows unexpected jumps. These variations are assumed to be related to the exothermic heat release and particle mass loss caused by the reaction. To explain this phenomenon, efforts to model this in simulations have increased. While current models can capture many of these physical processes, they currently assign a constant reaction shell thickness with little physical reasoning. This work remedies this problem by assigning a dynamic physically informed shell thickness to the reacting particles, using solid analysis. The shell thickness of the particles impacts the rate of change of reacted mass in the system, as well as the rate at which the particles react. The model is based on a simple stress–strain relationship and gives a dynamic assignment for when the reacting particle should begin to fracture. We compare our model to the previous computational and simulation data to analyze the effects of different model parameters.

42 ENGINEERING↗
DOE OSTI2025

Dynamic Carrier Modulation via Nonlinear Acoustoelectric Transport in van der Waals Heterostructures

Dynamically manipulating carriers in van der Waals heterostructures could enable solid-state quantum simulators with tunable lattice parameters. A key requirement is the formation of deep potential wells to reliably trap excitations. Here, we report the observation of nonlinear acoustoelectric transport and dynamic carrier modulation in boron nitride-encapsulated graphene devices coupled to intense surface acoustic waves (SAWs) on LiNbO 3 substrates. SAWs generate strong acoustoelectric current densities ( J AE ), transitioning from linear to nonlinear regimes with increasing SAW intensity. In the nonlinear regime, periodic carrier (electrons, holes, or their mixtures) stripes emerge. Using counter-propagating SAWs, we create standing SAWs (SSAWs) to dynamically manipulate charge distributions without static gates. The saturation of J AE , attenuation transitions, and tunable resistance peaks confirms strong carrier localization. Finally, these results establish SAWs as a powerful tool for controlling carrier dynamics in two-dimensional (2D) materials, paving the way for the development of time-dependent quantum systems and acoustic lattices for quantum simulation.

2D materials↗
NASA NTRS2019

In Situ Observation of Growth Dynamics in DECLIC Directional Solidification Insert Onboard ISS: DSI-R Flight Campaign

The study of solidification microstructure formation is of utmost importance for materials design and processing, as solid-liquid interface patterns largely govern mechanical and physical properties. Pattern selection occurs under dynamic conditions of growth in which the initial morphological instability evolves nonlinearly and undergoes a reorganization process. The dynamic and nonlinear nature of this instability renders in situ observation of the interface an invaluable tool to gain knowledge on the time-evolution of the interface pattern. Transparent organic analogs, which solidify like metallic alloys, allow direct visualization of interface dynamics. Extensive ground-based studies of both metallic and organic bulk samples have established the presence of significant convection during solidification processes that alters the formation of interfacial microstructures. A reduced-gravity environment is therefore mandatory for fluid flow elimination in bulk samples. In the framework of the CNES project MISOL3D (MIcrostrutures de SOLidification 3D) and the NASA projects DSIP (Dynamical Selection of 3D Interface Patterns), SPADES (SPAtiotemporal Evolution of three-dimensional DEndritic array Structures) and CAMUS (ComputAtional Studies of MicrostrUcture Formation During Alloy Solidification in Microgravity), we participated in the development of the Directional Solidification Insert (DSI) of the DEvice for the study of Critical Liquids and Crystallization (DECLIC). The DECLIC-DSI is dedicated to in situ and real-time characterization of solid-liquid interface patterns during directional solidification of transparent alloys in diffusive transport regime. Between April 2010 and March 2011, the first ISS campaign (DSI) explored the entire range of microstructures resulting in unprecedented observations. A second campaign (DSI-R), performed between October 2017 and December 2018, in which the insert contained an alloy of higher solute concentration, allowed to complete the benchmark database. The increase of solute concentration resulted in well-developed dendritic patterns at lower velocities (lower interface curvature and larger tip radius). The microstructure resulting from dendritic growth is dominant in metallurgy so that it is fundamental to understand the mechanisms of its formation. The main aims of this experimental campaign are to understand: the mechanisms of the cell to dendrite transition, the fundamental mechanisms of sidebranching formation, the dependence of dendrite tip shapes on growth conditions, the interaction of primary array and secondary sidebranches, and the influence of subgrain boundaries on the spatiotemporal organization of the array structure. Preparation, analysis and interpretation of the experiments performed onboard ISS are considerably enhanced by experiments performed on ground using thin-samples (Pr. Trivedi’s group) and phase-field simulations of microstructure formation in a diffuse growth regime (Pr. Karma’s group). In this summary, we will present an initial assessment of the results obtained during the DSI-R campaign. Nomenclature Vp: pulling velocity G: thermal gradient L: solidified length rtip: tip curvature radius Acronyms/Abbreviations CADMOS: Centre d'Aide au Développement des activités en Micro-pesanteur et des Opérations Spatiales CNES: Centre National d’Études Spatiales DECLIC: Device for the study of Critical LIquids and Crystallization DSI: Directional Solidification Insert DSI-R: Directional Solidification Insert - Refurbish ISS: International Space Station NASA: National Aeronautics and Space Administration PF: phase-field SCN: Succinonitrile 3D: three-dimensional

Mota, F.L.↗
DOE OSTI2024

Exploring Li-Ion Transport Properties of Li 3 TiCl 6 : A Machine Learning Molecular Dynamics Study

We performed large-scale molecular dynamics simulations based on a machine-learning force field (MLFF) to investigate the Li-ion transport mechanism in cation-disordered Li 3 TiCl 6 cathode at six different temperatures, ranging from 25°C to 100°C. In this work, deep neural network method and data generated by ab − initio molecular dynamics (AIMD) simulations were deployed to build a high-fidelity MLFF. Radial distribution functions, Li-ion mean square displacements (MSD), diffusion coefficients, ionic conductivity, activation energy, and crystallographic direction-dependent migration barriers were calculated and compared with corresponding AIMD and experimental data to benchmark the accuracy of the MLFF. From MSD analysis, we captured both the self and distinct parts of Li-ion dynamics. The latter reveals that the Li-ions are involved in anti-correlation motion that was rarely reported for solid-state materials. Similarly, the self and distinct parts of Li-ion dynamics were used to determine Haven’s ratio to describe the Li-ion transport mechanism in Li 3 TiCl 6 . Obtained trajectory from molecular dynamics infers that the Li-ion transportation is mainly through interstitial hopping which was confirmed by intra- and inter-layer Li-ion displacement with respect to simulation time. Ionic conductivity (1.06 mS/cm) and activation energy (0.29eV) calculated by our simulation are highly comparable with that of experimental values. Overall, the combination of machine-learning methods and AIMD simulations explains the intricate electrochemical properties of the Li 3 TiCl 6 cathode with remarkably reduced computational time. Thus, our work strongly suggests that the deep neural network-based MLFF could be a promising method for large-scale complex materials.

Selvaraj, Selva Chandrasekaran (ORCID:000000029023↗
DOE OSTI2024

eReaxFF force field development for BaZr 0.8 Y 0.2 O 3-δ solid oxide electrolysis cells applications

The use of solid-oxide materials in electrocatalysis applications, especially in hydrogen-evolution reactions, is promising. However, further improvements are warranted to overcome the fundamental bottlenecks to enhancing the performance of solid-oxide electrolysis cells (SOECs), which is directly linked to the more-refined fundamental understanding of complex physical and chemical phenomena and mass exchanges that take place at the surfaces and in the bulk of electrocatalysis materials. Here, we developed an eReaxFF force field for barium zirconate doped with 20 mol% of yttrium, BaZr 0.8 Y 0.2 O 3-δ (BZY20) to enable a systematic, large-length-scale, and longer-timescale atomistic simulation of solid-oxide electrocatalysis for hydrogen generation. All parameters for the eReaxFF were optimized to reproduce quantum-mechanical (QM) calculations on relevant condensed phase and cluster systems describing oxygen vacancies, vacancy migrations, electron localization, water adsorption, water splitting, and hydrogen generation on the surfaces of the BZY20 solid oxide. Using the developed force field, we performed both zero-voltage (excess electrons absent) and non-zero-voltage (excess electrons present) molecular dynamics simulations to observe water adsorption, water splitting, proton migration, oxygen-vacancy migrations, and eventual hydrogen-production reactions. Based on investigations offered in the present study, we conclude that the eReaxFF force field-based approach can enable computationally efficient simulations for electron conductivity, electron leakage, and other non-zero-voltage effects on the solid oxide materials using the explicit-electron concept. Moreover, we demonstrate how the eReaxFF force field-based atomistic-simulation approach can enhance our understanding of processes in SOEC applications and potentially other renewable-energy applications.

08 HYDROGEN↗
DOE OSTI2026

Exploring sodium dynamics in the dilute oxygen regime of mixed oxy-sulfide NaPSiSO glasses

All-solid-state‑sodium batteries have great potential to lower the cost of grid-scale energy storage systems. However, they require a low cost, high conductivity solid electrolyte to become a reality. The electrochemical properties of mixed oxy-sulfide glasses NaPSiSO exceed those of either the pure sulfide or the pure oxide glass alone. In this work, we have tested models of sodium ion dynamics in the NaPSiSO series: zNa 2 S+(1-z)[(1-y)[(1-x)SiS 2 + xPS 5/2 ] + yNaPO 3 ] where 0.57 ≤ z ≤ 0.588, 0 ≤ x ≤ 0.35, and 0.15 ≤ y ≤ 0.35 using 23 Na nuclear magnetic resonance to directly probe the local fields and their fluctuations at the nuclear sites. We determined the distribution of activation energies in these materials arising from distributions of local environments and provided quantitative information about the connectivity of the ionic pathways through which the sodium ions move. We found that the sodium ion conduction dynamics were described as an ionic percolation through a Gaussian distribution of energy barriers and that the composition dependence of the ionic conductivity in the dilute oxygen regime is well-described by the Nernst-Einstein relation modified for percolation.

25 ENERGY STORAGE↗
DOE OSTI2026

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Using Plate Finite Elements for Modeling Fillets in Design, Optimization, and Dynamic Analysis

A methodology has been developed that allows the use of plate elements instead of numerically inefficient solid elements for modeling structures with 90 degree fillets. The technique uses plate bridges with pseudo Young's modulus (Eb) and thickness (tb) values placed between the tangent points of the fillets. These parameters are obtained by solving two nonlinear simultaneous equations in terms of the independent variables rlt and twallt. These equations are generated by equating the rotation at the tangent point of a bridge system with that of a fillet, where both rotations are derived using beam theory. Accurate surface fits of the solutions are also presented to provide the user with closed-form equations for the parameters. The methodology was verified on the subcomponent level and with a representative filleted structure, where the technique yielded a plate model exhibiting a level of accuracy better than or equal to a high-fidelity solid model and with a 90-percent reduction in the number of DOFs. The application of this method for parametric design studies, optimization, and dynamic analysis should prove extremely beneficial for the finite element practitioner. Although the method does not attempt to produce accurate stresses in the filleted region, it can also be used to obtain stresses elsewhere in the structure for preliminary analysis. A future avenue of study is to extend the theory developed here to other fillet geometries, including fillet angles other than 90 and multifaceted intersections.

Brown, A. M.↗
DOE OSTI2021

Molecular dynamics simulations of energy accommodation between gases and polymers for ultra-low thermal conductivity insulation

Determining the energy accommodation between gases and solids is essential to developing porous thermal insulation materials with ultra-low effective thermal conductivity that reduce energy use, greenhouse gas emissions, and fossil fuel consumption. The energy accommodation coefficients of most gases, however, have been rarely studied, especially with respect to solids that have relatively high thermal resistivity, e.g., polymers. Here, by using all-atom nonequilibrium molecular dynamics simulations, we reveal the accommodation coefficients of He, Ar, N 2 , and O 2 with polymers, mainly polystyrene. We find that their values are around 0.51, 0.72, 0.79, and 0.90, respectively, suggesting a critical reexamination of the commonly used theoretical maximum value of 1. We have also conducted experiments and validated the value for air, which is about 0.81. Such a change in accommodation coefficients can lead to a reduction of about 70%, 50%, 35%, and 20% in the thermal conductivity of He, Ar, N 2 , and O 2 gases in nano pores (below 100 nm) or at low pressures (below 1 millibar). With these new accommodation coefficients, we find that in a 10 nm pore with ambient pressure at 300 K, the gas thermal conductivity of He, Ar, N 2 , and O 2 in porous polystyrene can be as low as 9.7 × 10 –4 , 3.4 × 10 –4 , 7.3 × 10 –4 , and 8.5 × 10 –4 W·m –1 ·K –1 , respectively, which are two to three orders of magnitude lower than their bulk values, promising higher thermal resistivity of insulation materials. This work reveals the fundamental energy exchange between gases and polymers, providing important guidance for designing high-performance thermal insulation materials for various applications.

36 MATERIALS SCIENCE↗
NASA NTRS2002

Current Progress of a Finite Element Computational Fluid Dynamics Prediction of Flutter for the AeroStructures Test Wing

This progress report focuses on the use of the STructural Analysis RoutineS suite program, SOLIDS, input for the AeroStructures Test Wing. The AeroStructures Test Wing project as a whole is described. The use of the SOLIDS code to find the mode shapes of a structure is discussed. The frequencies, and the structural dynamics to which they relate are examined. The results of the CFD predictions are compared to experimental data from a Ground Vibration Test.

Arena, Andrew S., Jr.↗
NASA NTRS1992

Delta II Geotail -- 1st Stage and Solid Motor Booster Erection

The Geotail mission's goal was to investigate the structure and dynamics of the geomagnetic tail that extends on the nightside of the Earth. The launch date was July 24, 1992. This video shows the Delta II on the pad, being prepared for the launch. The first stage and the solid motor booster are shown being moved into place on the rocket.

Source record↗
DOE OSTI2026

Time-Evolved Hetero-Alkali Interphases Enable Long-Life Sulfide-Based Anode-Free Solid-State Batteries

Sulfide-based anode-free solid-state batteries (AFSSBs) offer compelling advantages in terms of energy density and safety, yet their practical implementation is severely hindered by undesirable interfacial reactions between sulfide solid electrolytes (SEs) and freshly plated lithium (Li), as well as non-uniform Li plating/stripping behavior. Herein, an effective interfacial stabilization strategy by incorporating sodium bis(fluorosulfonyl)imide (NaFSI) additive into the Li5.4PS4.4Cl1.6 (LPSC) is investigated. Unlike conventional Li-based additives that form static passivation layers, NaFSI introduces a transient hetero-alkali chemistry that kinetically governs interphase evolution during fresh Li plating. NaFSI induces a timesequenced interphase evolution: an initial NaF/LiF-rich layer that suppresses early sulfide reduction, followed by a LiF/Li3N-rich layer that optimizes Li⁺ transport during repeated anode-free cycling. This evolved robust and fast ion conducting layer mitigates interfacial impedance growth, enhances Li + transport kinetics, and suppresses localized Li growth and filamentary shorting. As a result, the anode-free full cell with NaFSI modified LPSC as the interlayer exhibits an excellent cycling stability over 500 cycles at 0.2 C with a capacity retention of 77.6%, whereas the cell with bare LPSC suffers from rapid capacity decay after 100 cycles, retaining only 32.1% of its initial capacity. This work establishes dynamic heteroalkali additive chemistry as a general strategy to kinetically program solid-solid interphases, guiding the interface design in anode-free solid-state batteries.

25 ENERGY STORAGE↗
DOE OSTI2020

Atomistic simulations of the thermal conductivity of liquids

We present a method based on sinusoidal approach to equilibrium molecular dynamics (SAEMD) to compute the thermal conductivity of liquids Similar to nonequilibrium molecular dynamics, and unlike equilibrium simulations based on the Green-Kubo formalism, the method only requires the calculation of forces and total energies. The evaluation of heat fluxes and energy densities is not necessary, thus offering the promise of efficiently implementing first principles simulations based on density functional theory or deep molecular dynamics. Our approach is a generalization of SAEMD for solids, where the thermal conductivity is computed in the steady state, instead of a transient regime, thus properly taking into account diffusive terms in the heat equation. We present results for liquid water at ambient conditions and under pressure and discuss simulation requirements to obtain converged values of the thermal conductivity as a function of size and simulation time.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Unraveling the convoluted and dynamic interphasial mechanisms on Li metal anode

Accurate understanding of the chemistry of solid-electrolyte interphase (SEI) is key to developing new electrolytes for high-energy batteries using lithium metal (Li-0) anodes(1). SEI is generally believed to be formed by the reactions between Li-0 and electrolyte(2,3). However, our new study shows this is not the whole story. Through synchrotron-based X-ray diffraction and pair distribution function analysis, we reveal a much more convoluted formation mechanism of SEI, which receives considerable contributions from electrolyte, cathode, moisture and native surface species on Li-0, with highly dynamic nature during cycling. Using isotope labelling, we traced the origin of LiH to electrolyte solvent, moisture and a new source: the native surface species (LiOH) on pristine Li-0. When lithium accessibility is very limited as in the case of anode-free cells, LiOH develops into plate-shaped large crystals during cycling. Alternatively, when the lithium source is abundant, as in the case of Li||NMC811 cells, LiOH reacts with Li-0 to form LiH and Li2O. While the desired anion-derived LiF-rich SEI is typically found in the concentrated electrolytes or their derivatives, we found it can also be formed in low-concentration electrolyte via the crosstalk effect, emphasizing the importance of formation cycle protocol and opening up opportunities for low-cost electrolyte development.

Polzin, Bryant J.↗
DOE OSTI2022

An index to characterize gas‐solid and solid‐solid mixing from average volume fraction fields

Abstract A mixing index based on solid volume fraction fields is developed for gas‐solid flows. Conventional mixing indices are based on particle realizations of granular mixing and are applicable to experimental data or discrete element method simulations. However, these indices cannot be used as‐is for multifluid models, and an index for characterizing mixing in gas‐solid flows from continuous fields is needed. The performance of the new mixing index is tested in two applications. The first is a 3D simulation of the mixing of biomass and sand in a fluidized bed reactor, and the second is a 2D simulation of binary particle segregation in a fluidized bed. The simulations are performed using OpenFOAM®. The mixing index is used to quantify gas‐solid mixing using solid volume fractions and solid‐solid mixing using solid fractions. The formulation of conventional mixing indices is extended to be used with solid volume fractions fields, and methods for performance improvement are presented.

09 BIOMASS FUELS↗
DOE OSTI2023

Probing the Morphology and Mobility of Amines in Porous Silica CO 2 Sorbents by 1 H T 1 – T 2 Relaxation Correlation NMR

Molecular and oligomeric amines supported in porous oxide supports are a promising class of CO 2 sorbent materials studied for CO 2 removal from diverse streams such as flue gas and ambient air. Among the various amines investigated, low molecular weight, hyperbranched poly(ethyleneimine) (PEI), and tetraethylenepentamine (TEPA) are among the most extensively studied. While macroscopic structure–performance relationships relating the support structure, amine loading, and other factors affecting CO 2 sorption capacities and kinetics have been developed, structural and dynamic information about the organic amine phase in the porous support is less plentiful. The structure and mobility of amines impregnated in the pores of porous supports directly impact gas sorption, as the accessibility of amine sites in the pores directly relates to amine distribution in the pores and overall pore filling as well as the dynamics of the amine chains. Here, we prepare a family of mesoporous silica SBA-15 materials containing varying loadings of oligomeric (PEI) and molecular (TEPA) amines. 1 H T 1 –T 2 relaxation correlation solid-state NMR experiments are used to characterize the structural and dynamic properties of the confined amines. Both TEPA and PEI are shown to form multiple different domains in the pores, each with distinguishable dynamic properties. TEPA and PEI form more rigid layers around the silica support walls at lower organic loading fractions, characterized by lower mobilities, followed by the formation of more mobile domains less engaged in pore wall interactions at higher loadings. TEPA shows faster mobilities than PEI because of its lower molecular weight. TEPA also appears to more easily transfer between domains within the pores, leading to generally faster CO 2 uptake rates with higher sorption capacities, while PEI located closer to the pore walls remained much less mobile and is thus less engaged in CO 2 capture.

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

Mesoscale Modeling of Hydrogels Under Frictional Shear Stress

Hydrogels are three-dimensional networks of hydrophilic polymers often used as a simplified model of hydrated biological materials, from cartilaginous joints to the ocular tear film. However, the lubrication mechanisms of hydrogels remain poorly understood, partly due to their complex polymeric structure, which creates blurred interfaces during sliding that are challenging to study experimentally. In this study, we employ dissipative particle dynamics (DPD) to investigate the frictional behavior of a polymeric hydrogel network sliding against a solid wall in an explicit viscous solvent. This computational approach enables us to model hydrodynamic interactions and mesoscale polymer dynamics, capturing key aspects of hydrogel friction. Our simulations reveal that hydrogel friction is governed by the interplay between polymer relaxation and viscous shear, characterized by the Weissenberg number (Wi). At low Wi, friction coefficient remain nearly constant, dominated by polymer relaxation. However, at higher Wi, friction is dominated by viscous drag within a near-wall solvent layer, leading to a linear increase in friction coefficient with Wi. Furthermore, our results demonstrate an inverse relationship between the friction coefficient and the applied normal load, consistent with experimental observations. This work provides new insights into the fundamental tribological properties of hydrogels, shedding light on the micromechanics of hydrogel friction. Improving our understanding of hydrogel structure and dynamics under friction advances our knowledge of the mechanisms regulating biological lubrication in health and disease.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nanoscale Structure and Dynamics in Geochemical Systems

Neutron scattering is a powerful tool to elucidate the structure and dynamics of systems that are important to geochemists, including ion association in complex aqueous solutions, solvent-exchange reactions at mineral–water interfaces, and reaction and transport of fluids in nanoporous materials. This article focusses on three techniques: neutron diffraction, which can reveal the atomic-level structure of aqueous solutions and solids; quasi-elastic neutron scattering, which measures the diffusional dynamics at mineral–water interfaces; and small-angle neutron scattering, which can show how properties of nanoporous systems change during gas, liquid, and solute imbibition and reaction. Furthermore, the usefulness and applicability of the experimental results are extended by rigorous comparison to computational simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗