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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Benchmark of the LAMMPS code on CRESCO and SUMMIT HPC systems

In this work we present the benchmark results of a numerical model where the solid Deuterium grows under gas pressure. The two-phase molecular Deuterium systems (gas + solid) is modeled using classical molecular dynamics approach and the LAMMPS code. The performance between different modern supercomputers are analyzed and compared. In addition, the numerical simulations reach high accuracy to produce results that match with the data collected through laboratory experiments.

Lupo Pasini, Massimiliano↗
DOE OSTI2023

Degradation Modeling of SOEC Systems

This presentation was delivered at the 2023 AICHE Annual Meeting. This presentation focuses on dynamic modeling of physical degradation of solid oxide cell systems and failure analysis.

Le, Quang↗
DOE OSTI2024

High-Fidelity Arc-Discharge Model for Hydrogen-Plasma-Smelting-Reduction of Iron Ore

Electrification and use of renewable hydrogen is currently a necessity for decarbonizing the iron-and-steel industry. In this regard, hydrogen plasma smelting reduction (HPSR) is a novel pathway that is being explored for reduction of iron ore. HPSR provides several decarbonization merits compared to conventional blast furnaces. Firstly, the use of renewable hydrogen drastically reduces the CO2 emissions compared to the use of coke. Secondly, renewable electricity in the form of a thermal plasma for making reactive hydrogen species (radicals, ions) are more efficient at reducing iron ore compared to neutral H2. Thirdly, a molten product compatible with downstream processes is obtained from the intense heat transfer from the plasma. However, the scale-up of this technology requires fundamental exploration of hydrogen plasma dynamics and its interaction with complex solid material that include phase changing iron-ore and slag. In this work, we present a first principles continuum scale model for thermal plasmas in Ar/H2 gas mixtures typically used for HPSR. The thermal plasma governing equations for mass, momentum and energy with Lorentz force and Joule heating source terms are solved along with electromagnetic equations for electrostatic and magnetic vector potential. Our solver will be based on Pele, a suite of reacting flow solvers designed for advanced scientific computing architectures (Henry De Frahan et al., Proceedings of SIAM Parallel Processing, 13-25, 2024), and will utilize adaptive mesh generation for enhanced resolutions at locations of intense physicochemical interactions. This study will present the impact of Ar to H2 ratios on excited/dissociated hydrogen species concentrations, plasma temperature and conductivity along with the impact of outgassed species (water, metal vapor, O, OH radicals) from ore surface on gas phase chemistry. Furthermore, the heat and species flux to the surface will be quantified as a function of applied voltages in a transferred arc configuration.

hydrogen plasma↗
DOE OSTI2022

Ab Initio Phase Diagram of Chromium to 2.5 TPa

Chromium possesses remarkable physical properties such as hardness and corrosion resistance. Chromium is also a very important geophysical material as it is assumed that lighter Cr isotopes were dissolved in the Earth’s molten core during the planet’s formation, which makes Cr one of the main constituents of the Earth’s core. Unfortunately, Cr has remained one of the least studied 3d transition metals. In a very recent combined experimental and theoretical study (Anzellini et al., Scientific Reports, 2022), the equation of state and melting curve of chromium were studied to 150 GPa, and it was determined that the ambient body-centered cubic (bcc) phase of crystalline Cr remains stable in the whole pressure range considered. However, the importance of the knowledge of the physical properties of Cr, specifically its phase diagram, necessitates further study of Cr to higher pressure. In this work, using a suite of ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries, we obtain the theoretical phase diagram of Cr to 2.5 TPa. We calculate the melting curves of the two solid phases that are present on its phase diagram, namely, the lower-pressure bcc and the higher-pressure hexagonal close-packed (hcp) ones, and obtain the equation for the bcc-hcp solid–solid phase transition boundary. We also obtain the thermal equations of state of both bcc-Cr and hcp-Cr, which are in excellent agreement with both experimental data and QMD simulations. We argue that 2180 K as the value of the ambient melting point of Cr which is offered by several public web resources (“Wikipedia,” “WebElements,” “It’s Elemental,” etc.) is most likely incorrect and should be replaced with 2135 K, found in most experimental studies as well as in the present theoretical work.

equation of state↗
NASA NTRS2004

Thermodynamic Modeling of a Solid Oxide Fuel Cell to Couple with an Existing Gas Turbine Engine Model

The Controls and Dynamics Technology Branch at NASA Glenn Research Center are interested in combining a solid oxide fuel cell (SOFC) to operate in conjunction with a gas turbine engine. A detailed engine model currently exists in the Matlab/Simulink environment. The idea is to incorporate a SOFC model within the turbine engine simulation and observe the hybrid system's performance. The fuel cell will be heated to its appropriate operating condition by the engine s combustor. Once the fuel cell is operating at its steady-state temperature, the gas burner will back down slowly until the engine is fully operating on the hot gases exhausted from the SOFC. The SOFC code is based on a steady-state model developed by The U.S. Department of Energy (DOE). In its current form, the DOE SOFC model exists in Microsoft Excel and uses Visual Basics to create an I-V (current-voltage) profile. For the project's application, the main issue with this model is that the gas path flow and fuel flow temperatures are used as input parameters instead of outputs. The objective is to create a SOFC model based on the DOE model that inputs the fuel cells flow rates and outputs temperature of the flow streams; therefore, creating a temperature profile as a function of fuel flow rate. This will be done by applying the First Law of Thermodynamics for a flow system to the fuel cell. Validation of this model will be done in two procedures. First, for a given flow rate the exit stream temperature will be calculated and compared to DOE SOFC temperature as a point comparison. Next, an I-V curve and temperature curve will be generated where the I-V curve will be compared with the DOE SOFC I-V curve. Matching I-V curves will suggest validation of the temperature curve because voltage is a function of temperature. Once the temperature profile is created and validated, the model will then be placed into the turbine engine simulation for system analysis.

Brinson, Thomas E.↗
DOE OSTI2021

Improved modeling of the solid-to-plasma transition of polystyrene ablator for laser direct-drive inertial confinement fusion hydrocodes

We report that the target performance of laser direct-drive inertial confinement fusion (ICF) can be limited by the development of hydrodynamic instabilities resulting from the nonhomegeneous laser absorption at the target surface, i.e., the laser imprint on the ablator. To understand and describe the formation of these instabilities, the early ablator evolution during the laser irradiation should be considered. In this work, an improved modeling of the solid-to-plasma transition of a polystyrene ablator for laser direct-drive ICF is proposed. This model is devoted to be implemented in hydrocodes dedicated to ICF which generally assume an initial plasma state. The present approach consists of the two-temperature model coupled to the electron, ion and neutral dynamics including the chemical fragmentation of polystyrene. The solid-to-plasma transition is shown to significantly influence the temporal evolution of both free electron density and temperatures, which can lead to different shock formation and propagation compared with an initial plasma state. The influence of the solid-to-plasma transition on the shock dynamics is evidenced by considering the scaling law of the pressure with respect to the laser intensity. The ablator transition is shown to modify the scaling law exponent compared with an initial plasma state.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Examining normal modes as fundamental heat carriers in amorphous solids: The case of amorphous silicon

Normal mode decomposition of atomic vibrations has been used to provide microscopic understanding of thermal transport in amorphous solids for decades. In normal mode methods, it is naturally assumed that atoms vibrate around their equilibrium positions, and that individual normal modes are the fundamental vibrational excitations transporting heat. With the abundance of predictions from normal mode methods and experimental measurements now available, we carefully analyze these calculations in amorphous silicon, a model amorphous solid. We find a number of discrepancies, suggesting that treating individual normal modes as fundamental heat carriers may not be accurate in amorphous solids. Furthermore, our classical and ab initio molecular dynamics simulations of amorphous silicon demonstrate a large degree of atomic diffusion, especially at high temperatures, leading to the conclusion that thermal transport in amorphous solids could be better described starting from the perspectives of liquid physics rather than from crystalline solids.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Semi-coupled flow and thermal analysis of the field joint during rapid pressurization of the redesigned solid rocket motor

A transient, semi-coupled, multi-dimensional thermal and flow analysis methodology was developed to predict the thermal/gas dynamic conditions in the field joint region of the Redesigned Solid Rocket Motor (RSRM). Transient temperature response, pressure history, and combustion gas flow rate (within the field joint region), were of principle interest, in the course of this study. The thermal environment in the field joint was modeled using SINDA, a finite difference based thermal network analyzer. The combustion gas flow boundary condition was generated using the FLAP code; this code performs a transient, lumped-parameter, control volume analysis to solve the mass, momentum, and energy conservation equations. The FLAP computer code was modified to account for erosion of the NBR insulation material, following ignition. An independent grid sensitivity study was conducted to determine an appropriate grid distribution near the wall. The predicted results, obtained using an optimum grid distribution and computer generated flow boundary condition, were compared with subscale test data.

Ghaffarian, Benny↗
DOE OSTI2026

Dataset_for_Conserved_macromolecular_architecture_of_Poplar_secondary_cell_walls_revealed_by_ssNMR_and_atomistic_modeling

This dataset contains solid-state 13C NMR data and atomistic molecular dynamics simulation files supporting the study of nanoscale secondary cell wall architecture across 13 genetically diverse Populus trichocarpa genotypes grown under uniform greenhouse conditions in 13C-enriched CO2 atmospheres (~89% 13C enrichment).The dataset contains two collections of solid-state 13C NMR data. (1) 200 MHz data (Bruker Avance III HD, 4 mm HX probe, 10 kHz MAS): raw Bruker TopSpin experiment folders and DMFIT-exported ascii spectra for selective and non-selective 1D 13C-13C spin diffusion experiments (3000 ms mixing) used to quantify inter-polymer spatial proximities, and short-mixing (1 ms) reference spectra used for polymeric abundance quantification by spectral deconvolution. (2) 600 MHz data (Bruker Avance III, 1.6 mm PhoenixNMR HXY probe, 30 kHz MAS): raw Bruker TopSpin experiment folders containing 2D CORD, 2D CP-INADEQUATE, and 13C/1H relaxation (T1, T1rho) experiments for all 13 genotypes, with processed Excel workbooks per experiment type. Molecular dynamics simulation code, coordinate files, and analysis scripts (NAMD/CHARMM/Python) for six atomistic cell wall models are included. Summarized ssNMR data are compiled into a single excel file and subjected to statistical analysis. Multivariate analysis code (PCA, Pearson correlation) and summary data are provided as excel worksheets and Jupyter notebooks (Python 3).

09 BIOMASS FUELS↗
DOE OSTI2024

Ultrafast high-harmonic spectroscopy of solids

High-harmonic spectroscopy, an ultrafast all-optical technique initially conceptualized in atomic and molecular systems, has now emerged as a powerful platform for studying the structure and dynamics of condensed matter. Unlike that in the gas phase, solid-state high-harmonic generation relies on the fundamental response from high atomic density and periodicity, leading to interband transitions and coherent driving of electrons and holes in their respective bands. These mechanisms make high-harmonic spectroscopy particularly sensitive to the electronic band structure, topological properties and many-body correlations in condensed media. An advantage of high-harmonic spectroscopy over other spectroscopic methods is its ability to probe ultrafast phenomena, capturing femto- to attosecond dynamics of multi-band and strongly correlated electron interactions in solids. Furthermore, in this Review, we discuss the latest experimental and theoretical advances in ultrafast high-harmonic spectroscopy of solids and provide perspectives for future research in this field.

High-harmonic generation↗
DOE OSTI2023

Atomistic Origins of Reversible Noncatalytic Gas–Solid Interfacial Reactions

Noncatalytic gas-solid reactions are a large group of heterogeneous reactions that are usually assumed to occur irreversibly because of the strong driving force to favor the forward direction toward the product formation. Using the example of Ni oxidation into NiO with CO 2 , herein we demonstrate the existence of the reverse element that results in the NiO reduction from the countering effect of the gaseous product of CO. Using in-situ electron microscopy observations and atomistic modeling, we show that the oxidation process occurs via preferential CO 2 adsorption along step edges that results in step-flow growth of NiO layers, and the presence of Ni atoms on the flat NiO surface promotes the nucleation of NiO layers. Simultaneously, the NiO reduction happens via preferential step-edge adsorption of CO that leads to the receding motion of atomic steps, and the presence of Ni vacancies in the NiO surface facilitates the CO adsorption induced surface pitting. Based on the in-situ TEM observations, temperature and CO 2 pressure effect maps are constructed to illustrate the spatiotemporal dynamics of the competing NiO redox reactions. We report, these results demonstrate the rich gas-solid surface reaction dynamics induced by the co-existing forward and reverse reaction elements and have practical applicability in manipulating gas-solid reactions via controlling the gas environment or atomic structure of the solid surface to steer the reaction toward the desired direction.

25 ENERGY STORAGE↗
DOE OSTI2021

Understanding Solid Electrolyte Interphase Nucleation and Growth on Lithium Metal Surfaces

Experiments and theory are needed to decode the exact structure and distribution of components of a passivation layer formed at the anode surface of Li metal batteries, known as the Solid Electrolyte Interphase (SEI). Due to the inherent dynamic behavior as well as the lithium reactivity, the SEI structure and its growth mechanisms are still unclear. This study uses molecular simulation and computational chemistry tools to investigate the initial nucleation and growth dynamics of LiOH and Li2O that provide us with thermodynamics and structural information about the nucleating clusters of each species. Following the most favorable pathways for the addition of each of the components to a given nascent SEI cluster reveals their preferential nucleation mechanisms and illustrates different degrees of crystallinity and electron density distribution that are useful to understand ionic transport through SEI blocks.

25 ENERGY STORAGE↗
DOE OSTI2020

Improved Single-Ion Conductivity of Polymer Electrolyte via Accelerated Segmental Dynamics

Single-ion conducting polymer electrolytes (SICPEs) have many advantages for solid-state battery applications, while low conductivities remain the bottleneck for their practical applications. In this paper, a strategy that can significantly improve the ionic conductivity of SICPEs based on the concept of accelerated segmental dynamics was demonstrated by the evaluation of the thermal property, rheological behavior, morphology analysis, molecular dynamics simulation, and conductivity performance. With both “soft” poly(dimethylsiloxane) (PDMS) backbone and poly(ethylene glycol) (PEG) side chains, the obtained SICPE possesses faster segmental relaxation and higher lithium-ion conductivity. With lithium “transference” number close to unity, the obtained SICPE also shows excellent electrochemical stability against lithium metal electrodes. The clear relationship established between the segmental dynamics of polymer electrolyte and its ionic conductivity should contribute to achieve a solid electrolyte with improved ionic conductivity toward the next-generation solid-state battery.

Segmental Dynamics↗
DOE OSTI2026

Elucidating Lithium Transport Mechanisms in Disordered LiF from Machine-Learning Molecular Dynamics Simulations

Lithium fluoride (LiF) is a ubiquitous component of solid- and cathode–electrolyte interphases, yet its functional role remains unclear under the structural and chemical heterogeneity typical of cycling batteries. Here, we systematically quantify how structural disorder, off-stoichiometry, and strain govern Li-ion transport in LiF. Using a machine-learning potential to enable extensive molecular-dynamics sampling, we compare crystalline and amorphous LiF, Li 0.95 F, and LiF 0.95 , and evaluate the impact of small homogeneous deformations. Defect-free crystalline LiF is effectively ion-blocking at 300−500 K, whereas amorphization generates free-volume–assisted percolation pathways that facilitates Li-ion diffusion. At elevated temperatures, thermodynamically driven crystallization disrupts these pathways, leading to non-Arrhenius behavior. In crystalline phases, Li deficiency activates vacancy-mediated diffusion, while in amorphous LiF, transport is governed primarily by network connectivity. Strain is found to have only a marginal effect on Li mobility in both crystalline and amorphous structures.

Batteries↗
NASA NTRS1992

An analysis of the flow field in the region of the ASRM field joints

The flow field in the region of a solid rocket motor field joint is very important since fluid dynamic and mechanical propellant stresses can couple to cause a motor failure at a joint. Presented here is an examination of the flow field in the region of the Advanced Solid Rocket Motor (ASRM) field joints. The analyses were performed as a first step in assessing the design of the ASRM forward and aft field joints in order to assure the proper operation of the motor prior to further development of test firing. The analyses presented here were performed by employing a two-dimensional axisymmetric assumption. Fluent/BFC, a three dimensional full Navier-Stokes flow field code, was used to make the numerical calculations. This code utilizes a staggered grid formulation along with the SIMPLER numerical algorithm. Wall functions are used to determine the character of the laminar sublayer, and a standard kappa-epsilon turbulence model is used to close the fluid dynamic equations. The analyses performed to this date verify that the ASRM field joint design operates properly. The fluid dynamic stresses at the field joints are small due to the inherent design of the field joints. A problem observed in some other solid rocket motors is that large fluid dynamic stresses are generated at the motor joint on the downstream propellant grain due to forward facing step geometries. The design of the ASRM field joints are such that this is not a problem as shown by the analyses. Also, the analyses of the inhibitor stub left protruding into the port flow from normal propellant burn back show that more information is necessary to complete these analyses. These analyses were performed as parametric analyses in relation to the height of the inhibitor stub left protruding into the motor port. A better estimate of the amount of the inhibitor stub remaining at later burn times must be determined since the height which the inhibitor stub protrudes into the port flow drastically affects the fluid dynamics induced stresses on the propellant grain at the field joints.

Dill, Richard A.↗
DOE OSTI2026

Dynamic magneto-chiral instability in photoexcited tellurium

In systems of charged chiral fermions out of equilibrium, an electric current parallel to a magnetic field can generate a dynamic instability that amplifies electromagnetic waves. Whether this mechanism also operates in chiral solid-state systems has remained uncertain. Here we observe signatures of a dynamic magneto-chiral instability in elemental tellurium, a structurally chiral crystal, using time-domain terahertz emission spectroscopy. Under transient photoexcitation in a moderate magnetic field, we observe terahertz radiation with coherent modes that grow in amplitude over time. We present a theoretical model that describes this behaviour based on a dynamic instability of electromagnetic waves interacting with infrared-active oscillators of acceptor states in tellurium, giving rise to an amplifying polariton. These results demonstrate that magneto-chiral instabilities can emerge in solid-state systems and establish a mechanism for terahertz-wave amplification in chiral materials.

42 ENGINEERING↗
NASA NTRS2002

New Method Developed to Measure Contact Angles of a Sessile Drop

The spreading of an evaporating liquid on a solid surface occurs in many practical processes and is of importance in a number of practical situations such as painting, textile dyeing, coating, gluing, and thermal engineering. Typical processes involving heat transfer where the contact angle plays an important role are film cooling, boiling, and the heat transfer through heat pipes. The biological phenomenon of cell spreading also is analogous to a drop spreading (ref. 1). In the study of spreading, the dynamic contact angle describes the interfacial properties on solid substrates and, therefore, has been studied by physicists and fluid mechanics investigators. The dynamic contact angle of a spreading nonvolatile liquid drop provides a simple tool in the study of the free-boundary problem, but the study of the spreading of a volatile liquid drop is of more practical interest because the evaporation of common liquids is inevitable in practical processes. The most common method to measure the contact angle, the contact radius, and the height of a sessile drop on a solid surface is to view the drop from its edge through an optical microscope. However, this method gives only local information in the view direction. Zhang and Yang (ref. 2) developed a laser shadowgraphy method to investigate the evaporation of sessile drop on a glass plate. As described here, Zhang and Chao (refs. 3 and 4) improved the method and suggested a new optical arrangement to measure the dynamic contact angle and the instant evaporation rate of a sessile drop with much higher accuracy (less than 1 percent). With this method, any fluid motion in the evaporating drop can be visualized through shadowgraphy without using a tracer, which often affects the field under investigation.

Chao, David F.↗