SEARCH · Engineering Papers
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.
In this poster we present our preliminary analysis on constraining cosmology parameters such as the Dark Energy Equation of State parameter using strong lensing. We use Simulations based inference with Machine Learning algorithms such as Neural Posterior density estimators.
Shortly after the Big Bang the entire universe was filled with a nearly perfect fluid known as the Quark Gluon Plasma. Relativistic heavy ion collisions can now reproduce this fluid in the laboratory where the Quark Gluon Plasma exhibits a rapid but smooth cross-over phase transition into hadrons at vanishing net-baryon densities. Recent experiments plan to explore finite baryon densities where a critical point is expected. If found, this would mark the first discovery of a critical point in a relativistic system described by a fundamental theory of nature, which would have far-reaching consequences for high-energy nuclear physics and nuclear astrophysics (such as in neutron star mergers). Characteristic temperatures of equilibrium (e.g. the inflection point of the entropy density) and transport coefficients (e.g. minimum of the shear viscosity over entropy density) vary widely at a cross-over phase transition but converge at a critical point, and extracting the behavior of these characteristic temperatures is a major focal point of this research. Specifically, the interplay between strange and light hadrons is exploited to study the flavor hierarchy in the cross-over region. To investigate this, a viscous relativistic hydrodynamics framework with two conserved charges is being developed into a new open-source code along with initial conditions that contain baryon number and strangeness. Flow observables sensitive to the equation of state and transport coefficients are calculated across beam energies. New techniques are being developed to study this flavor hierarchy from first principles and to extract the characteristic temperatures from experimental data. Through this new dynamical framework, this project provides essential guidance to the Beam Energy Scan II runs at Relativistic Heavy-Ion Collider and the future Facility for Antiproton and Ion Research facility in the search for the Quantum Chromodynamic critical point and subsequent investigation of the baryon-rich Quark Gluon Plasma.
This is the Final report of the grant Materials at Extreme Pressures and Temperatures. The proposal for this grant defined efforts to understand the role of inner-shell ionization on the equation of state of light elements under conditions of extreme pressure and temperature. The conditions of high energy and high density are described by several theories, and this work was aimed at comparing theories with experimental results from Discovery Science experiments on the National Ignition Facility (NIF) at the Lawrence Livermore National Laboratory.
Lithium Hydride (LiH) atomic properties make it an excellent candidate for thermal energy storage, hydrogen storage, and nuclear reactor technology. High Energy Density (HED) experiments on LiH at the National Ignition Facility (NIF) can also provide critical Equation of State (EOS) and Hugoniot data. Density requirements for LiH vary per application and therefore physical and structural requirements for LiH are non-uniform. Historical production of LiH relies on casting processesthat are limited by density and grain-size control, which leads to unfavorable machinable characteristics. LiH manufacturing is further hindered due to its chemical reactivity and propensity to readily react with moisture. This work aimed at exploring advanced fabrication techniques for LiH such as 1) uniaxial powder pressing methods and 2) laser powder bed fusion (L-PBF) additive manufacturing. Uniaxial powder compaction offers advantageous tailorable mechanical and physical properties via density control, and L-PBF can produce net-shaped complex parts with unique microstructure. Uniaxial pressing proved successful in readily producing LiH with density control.
Lithium Hydride (LiH) atomic properties make it an excellent candidate for thermal energy storage, hydrogen storage, and nuclear reactor technology. High Energy Density (HED) experiments on LiH at the National Ignition Facility (NIF) can also provide critical Equation of State (EOS) and Hugoniot data. Density requirements for LiH vary per application and therefore physical and structural requirements for LiH are non-uniform. Historical production of LiH relies on casting processes that are limited by density and grain-size control, which leads to unfavorable machinable characteristics. LiH manufacturing is further hindered due to its chemical reactivity and propensity to readily react with moisture. This work aimed at exploring advanced fabrication techniques for LiH such as 1) uniaxial powder pressing methods and 2) laser powder bed fusion (L-PBF) additive manufacturing. Uniaxial powder compaction offers advantageous tailorable mechanical and physical properties via density control, and L-PBF can produce net-shaped complex parts with unique microstructure. Uniaxial pressing proved successful in readily producing LiH with density control.
In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a 10 multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory (DFT). The shear modulus model for the ambient β phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the β phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress-strain response of tin was measured at strain rates of 10 -3 to 3 x 10 3 s -1 and temperatures ranging from 87 to 373 K. The Preston-Tonks-Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to 10 5 s -1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability 20 of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid-liquid but never fully liquid state, never fully lost strength, 25 stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.
A cylinder of the experimental HE is detonated inside of an aquarium and sequential fast-frame images are taken of the resulting shock wave. We wish to use the sequence of images to extract quantitative data regarding the detonation velocity and velocity of the shock wave which can then be used to tune the parameters of the equation of state for the experimental HE material.
Machine learning methods are increasingly deployed to construct surrogate models for complex physical systems at a reduced computational cost. However, the predictive capability of these surrogates degrades in the presence of noisy, sparse, or dynamic data. We introduce an online learning method empowered by optimizer-driven sampling that has two advantages over current approaches: it ensures that all local extrema (including endpoints) of the model response surface are included in the training data and it employs a continuous validation and update process in which surrogates undergo retraining when their performance falls below a validity threshold. We find, using benchmark functions, that optimizer-directed sampling generally outperforms traditional sampling methods in terms of accuracy around local extrema even when the scoring metric is biased towards assessing overall accuracy. Finally, the application to dense nuclear matter demonstrates that highly accurate surrogates for a nuclear equation of state model can be reliably auto-generated from expensive calculations using few model evaluations.
In this work, we present a scalable approach for inferring the dark energy equation-of-state parameter ($w$) from a population of strong gravitational lens images using Simulation-Based Inference (SBI). Strong gravitational lensing offers crucial insights into cosmology, but traditional Monte Carlo methods for cosmological inference are computationally prohibitive and inadequate for processing the thousands of lenses anticipated from future cosmic surveys. New tools for inference, such as SBI using Neural Ratio Estimation (NRE), address this challenge effectively. By training a machine learning model on simulated data of strong lenses, we can learn the likelihood-to-evidence ratio for robust inference. Our scalable approach enables more constrained population-level inference of $w$ compared to individual lens analysis, constraining $w$ to within $1\sigma$.
This report presents analyses of the AB5 and AB6 ABCOVE sodium spray fire experiments with the MELCOR code. This code simulates the progression of accident events for analysis and auditing purposes of nuclear facilities during accident conditions. Historically, the ABCOVE experiments have contributed to the validation of aerosol physics and related phenomena. Given advancements in sodium-cooled reactor designs, characterization of the sodium spray combustion may further the review and validation of newly incorporated sodium properties and physics packages, namely, the sodium equations of state (EOS) and the sodium combustion (NAC) package within MELCOR. By analyzing the AB5 and AB6 experiments with and without the NAC package, sodium specificity for spray combustion and aerosol formation as well as speciation of the combustion products are reviewed with the new packages. This effort provides code users with a demonstration of the current code capabilities. This report provides the current best practices for the NAC package as well as a discussion of any issues observed while performing the presented analyses.
A long-standing area of research for Eulerian shock wave physics codes has been the treatment of strength and damage for materials. Here we present a method that will aid in the analysis of strength and failure in shock physics applications where excessive diffusion of critical variables can occur and control the solution outcome. Eulerian methods excel for large deformation simulations in general but are inaccurate in capturing structural behavior. Lagrangian methods provide better structural response, but finite element meshes can become tangled. Therefore, a technique for merging Lagrangian and Eulerian treatments of material response, within a single numerical framework, was implemented in the Multiple Component computational shock physics hydrocode. The capability is a Lagrangian/Eulerian Particle Method (LEPM) that uses particles to interface a Lagrangian treatment of material strength with a more traditional Eulerian treatment of the Equation of State (EOS). Lagrangian numerical methods avoid the advection diffusion found in Eulerian methods, which typically strongly affects strength constitutive law internal variables, such as equivalent plastic strain, porosity and/or damage. The Lagrangian capability enhances existing capabilities and permits accurate predictions of high rate, large deformation and/or shock of mechanical structures.
Piezoelectric materials are used as a power source in stockpile components for safety and reliability assurances. The objective of this project is to gain insights into the fundamental pressure-induced behavior and electromechanical response of lead zirconate titanate (PZT). Specifically, to establish the basis for an accurate, physics-based material model for robust model-based design to rapidly optimize PZT-based materials and components for performance studies. Established ab-initio methods are used to interrogate and understand the dynamic behavior of PZT as a function of composition (50/50, 65/35, 80/20, 95/5) and dopants (La, Nb) to overcome the costly and time-consuming experimental methodologies. New cold curves for pure and doped single-crystal PZT are obtained as the reference equation-of-state (EOS). A negligible change in the pressure responses was observed for the systems and strains studied. Dielectric and piezoelectric responses of pure and doped single-crystal PZT were also calculated as a function of pressure. For undoped PZT, there is a clear and orderly pressure and composition dependence. For doped PZT, there is a significant increase in the responses, but the behavior is very disordered and inconclusive.
The shock-induced reaction of high explosives to gaseous and solid reaction products is a rapid, complex process. Understanding solid reaction product structure and formation kinetics is essential in determining the high-pressure equation of state of these multicomponent systems. We use the National Ignition Facility (NIF) to shock compress ~500-µm thick pressed powder high-explosive TATB samples to 70-135 GPa and collect in situ structural data on detonation byproducts using the TARDIS X-ray diffraction diagnostic. Velocimetry is used to record the transmitted compression wave profile, which is well described by Cheetah hydrocode simulations coupled with a reactive flow model, providing strong evidence of reaction in the TATB sample. While an unambiguous determination of the product phases was not possible owing to the low signal-to-noise quality of the diffraction signal, X-ray diffraction data of the product phases formed within the first 50 ns of the reaction is most consistent with a mixture of amorphous products and crystalline hexagonal diamond.
During this funding period I continued my work on far-from-equilibrium formulations of hydrodynamics and bottomonium suppression in the quark-gluon plasma. I published 23 peer-reviewed articles and 9 proceedings contributions. The publications produced during this period have thus far received 705 citations. I delivered 22 talks, seminars, and colloquia about my work, most of which were invited. In addition to the work on hydrodynamics and heavy quarkonia, I also worked on high-order perturbative calculations of N=4 supersymmetric Yang-Mills thermodynamics and contributed to the efforts of the MUSES collaboration whose focus is developing a modular equation-of-state framework for use in modeling neutron stars and neutron star mergers.
This report critically reviews the flow behavior of hydrogen-natural gas (H 2 -NG) mixtures in pipelines and examines the critical factors of hydrogen integration into existing natural gas infrastructure. It addresses the choking behavior characterized by velocity increase and pressure drop, as well as the effects of flow restrictions and pressure losses during hydrogen transport. Computational and analytical models are used to investigate these effects, and their effects on thermodynamic properties and system performance are evaluated. The study also reviews the energy efficiency and flow dynamics of hydrogen and methane-hydrogen mixtures and optimizes the hydrogen flow rate. In addition, the effects of these mixtures on the flow characteristics are discussed in detail, with special emphasis on the compressibility factor (z factor) and fluid properties based on equations of state for hydrogen-natural gas mixtures. The study also analyzes the mixture ratios and highlights the thermophysical properties, flow dynamics, and hydrogen-blended natural gas application potential. These investigations assess flow stability, material interactions, and operational feasibility of transporting hydrogen mixtures through natural gas pipelines, which contribute to developing sustainable and efficient energy systems.
This report presents analyses of the AB7 ABCOVE sodium spray fire experiment with the MELCOR code. This code simulates the progression of accident events for analysis and auditing purposes of nuclear facilities during accident conditions. Historically, the ABCOVE experiments have contributed to the validation of aerosol physics and related phenomena. Given advancements in sodium-cooled reactor designs, characterization of the sodium spray combustion may further the review and validation of newly incorporated sodium properties and physics packages, namely, the sodium equations of state (EOS) and the sodium combustion (NAC) package within MELCOR. Previously, the AB5 and AB6 experiments were analyzed with and without the NAC package. This work builds on the previous analyses with a demonstration of the current code capabilities of MELCOR with a more mild Na spray and pool fire scenario.
Brittle behavior of metal alloys is often critical to modeling ballistic impact and penetration. The ALEGRA multiphysics finite element software incorporates calibrated models for the equation of state, elasticity, yield stress, plasticity and fracture, but simulations do not always capture expected metal fracture. Here we report concerted efforts to do so for one important case where experiments clearly show shear fractures: a tungsten sphere impacting a steel plate at various angles. Our best simulations show fractures that are qualitatively similar to experiments, but there are significant differences in quantitative metrics. Specifically, velocities of tracers used to quantify simulated plug parameters consistently fall short of measured plug velocities. Also, simulated plugs break apart more than expected from experimental evidence. We attribute these shortfalls to the lack of an explicit shear fracture mechanism in the material models, leading to over-estimated resistance to plug formation and movement.
The dynamic strength and equation of state of an epoxy (Epon 828 cured with diethanolamine) is characterized through combined experimental and numerical analysis.