Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equation of State”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effect of lot microstructure variations on detonation performance of the triaminotrinitrobenzene (TATB)-Based insensitive high explosive PBX 9502

PBX 9502 is an important insensitive high explosive due to its combination of safety properties and detonation performance. It is a polymer-bonded formulation consisting of 95 wt.% 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) as the high explosive crystal, bound with Kel F-800 (FK-800), a co-polymer of chlorotrifluoroethylene and vinylidene-fluoride. Two different types are used, one known as virgin that uses only pristine manufactured TATB, and a second known as recycled that has 50 wt.% of its TATB reclaimed from machining scraps of previously pressed virgin PBX 9502. Recycled lots have a higher percentage of fine particles compared to virgin lots, due to the fracturing and damage sustained by TATB crystals during pressing. We examine the influence that this TATB microstructure difference between virgin and recycled PBX 9502 has on detonation performance properties. New rate-stick geometry diameter effect, detonation front shape and cylinder expansion test data are obtained for two previously uncharacterized virgin lots and one recycled lot of PBX 9502. This is combined with previously published data for one virgin and one recycled lot for evaluation purposes. Detonation shock dynamics model calibrations are conducted on each of the five lots to provide an assessment of the detonation timing characteristics of virgin versus recycled PBX 9502 lots. For two of the virgin lots, detonations propagate slower in the rate-stick geometry than those in the two recycled lots. However, the other virgin lot tested has propagation rates comparable to that of the recycled lots for larger diameter rate-sticks. The latter result, though, is shown to be geometry dependent and depends on the range of detonation curvatures accessed in different geometries. New copper-confined, cylinder expansion tests are conducted on each of the three virgin and two recycled lots to obtain detonation product Jones-Wilkins-Lee equations of state, enabling an assessment of the metal push capabilities for each of the five lots. We find that the metal push capabilities, characterized by the evolution of the heat of detonation with volume, are similar between the virgin and recycled lots. Thus, changes in microstructure between different PBX 9502 lots seemingly affect the rate of reaction, but not the overall energy content.

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

Development of slurry targets for high repetition-rate x-ray free electron laser experiments

Combining an x-ray free electron laser with a high-power laser driver enables the study of equations-of-state, high strain-rate deformation processes, structural phase transitions, and transformation pathways as a function of pressure to hundreds of GPa along different thermodynamic compression paths. Future high repetition-rate laser operation will enable data to be accumulated at >1 Hz, which poses a number of experimental challenges, including the need to rapidly replenish the target. Here, we present a combined shock compression and an x-ray diffraction study on epoxy (50% vol.)-crystalline grains (50% vol.) slurry targets, which can be fashioned into extruded ribbons for high repetition-rate operation. For shock-loaded NaCl-slurry samples, we observe pressure, density, and temperature states within the embedded NaCl grains consistent with observations from shock-compressed single-crystal NaCl.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Modeling atomically mixed graded density impactors

Graded density impactors (GDIs) are multi-material composite impactors used in gas gun experiments to tailor the drive conditions imparted to a sample test material. Previous graded density impactors generally rely on thin, but discrete, layers of different materials. The thinner and the greater number of layers will result in smoother compression. Taken to the limit of very thin layers would be pure material 1 at one surface, such as the front surface of an impactor, smoothly transitioning at the atomic scale to pure material 2 on the back surface. Such an impactor can initially shock, then smoothly compresses a material during a dynamic experiment. This type of experiment can serve to explore a larger region of thermodynamic space than a single or even multi-shock experiments. An overview of how graded density impactors are made is reviewed and sample results are given. Additionally, a strategy for modeling these kinds of impactors is presented. The length scales of constituent mixing are given from the experimental build through electrochemical-deposition. Equation of state models for pure constituents and their subsequent mixtures are presented. It is demonstrated that the time scales for pressure and temperature equilibration, for atomically mixed GDIs, are short enough to be a justifiable closure for the resulting multiphase flow. Furthermore, we present simulation results of dynamic shock followed by a ramp compression, utilizing a silver/gold graded density impactor, onto a tantalum sample.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

MechElastic: A Python library for analysis of mechanical and elastic properties of bulk and 2D materials

We report the MechElastic Python package evaluates the mechanical and elastic properties of bulk and 2D materials using the elastic coefficient matrix ( C ij ) obtained from any ab-initio density-functional theory (DFT) code. The current version of this package reads the output of VASP, ABINIT, and Quantum Espresso codes (but it can be easily generalized to any other DFT code) and performs the appropriate post-processing of elastic constants as per the requirement of the user. This program can also detect the input structure's crystal symmetry and test the mechanical stability of all crystal classes using the Born-Huang criteria. Various useful material-specific properties such as elastic moduli, longitudinal and transverse elastic wave velocities, Debye temperature, elastic anisotropy, 2D layer modulus, hardness, Pugh's ratio, Cauchy's pressure, Kleinman's parameter, and Lame's coefficients, can be estimated using this program. Another existing feature of this program is to employ the ELATE package (2016) [29] and plot the spatial variation of several elastic properties such as Poisson's ratio, linear compressibility, shear modulus, and Young's modulus in three dimensions. Further, the MechElastic package can plot the equation of state (EOS) curves for energy and pressure for a variety of EOS models such as Murnaghan, Birch, Birch-Murnaghan, and Vinet, by reading the inputted energy/pressure versus volume data obtained via numerical calculations or experiments. This package is particularly useful for the high-throughput analysis of elastic and mechanical properties of materials.

2D materials↗
DOE OSTI2024

Emulating ab initio computations of infinite nucleonic matter

We construct efficient emulators for the computation of the infinite nuclear matter equation of state. These emulators are based on the subspace-projected coupled-cluster method for which we here develop a new algorithm called small-batch voting to eliminate spurious states that might appear when emulating quantum many-body methods based on a non-Hermitian Hamiltonian. The efficiency and accuracy of these emulators facilitate a rigorous statistical analysis within which we explore nuclear matter predictions for > 10 6 different parametrizations of a chiral interaction model with explicit Δ -isobars at next-to-next-to leading order. Constrained by nucleon-nucleon scattering phase shifts and bound-state observables of light nuclei up to He 4 , we use history matching to identify nonimplausible domains for the low-energy coupling constants of the chiral interaction. Within these domains we perform a Bayesian analysis using sampling and importance resampling with different likelihood calibrations and study correlations between interaction parameters, calibration observables in light nuclei, and nuclear matter saturation properties. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Large and massive neutron stars: Implications for the sound speed within QCD of dense matter

Here, the NASA telescope NICER has recently measured x-ray emissions from the heaviest of the precisely known two-solar mass neutron stars, PSR J0740 + 6620. Analysis of the data suggests that PSR J0740 + 6620 has a radius in the range of R 2.0 ≈ (11.4 – 16.1) km at the 68% credibility level. In this article, we study the implications of this analysis for the sound speed in the high-density inner cores by using recent chiral effective field theory (χEFT) calculations of the equation of state at next-to-next-to-next-to-leading order to describe outer regions of the star at modest density. We find that the lower bound on the maximum speed of sound in the inner core, min {c$^{2}_{s,max}$}, increases rapidly with the radius of massive neutron stars. If χ EFT remains an efficient expansion for nuclear interactions up to about twice the nuclear saturation density, R 2.0 ≥ 13 km requires min {c$^{2}_{s,max}$} ≥ 0.562 and 0.442 at the 68% and 95% credibility level, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Deducing the EOS of dense neutron star matter with machine learning

Abstract The interior of a neutron star is a unique astrophysical laboratory for studying matter at extreme densities and pressures beyond what is replicable in terrestrial experiments. While there is no direct way to simulate the interior of these stars, one promising avenue to learning more about the equation of state (EOS) of such matter is through X‐rays emitted from the star's surface. The current state‐of‐the‐art method for inference of EOS from a star's X‐ray spectra uses piece‐wise, simulation‐based likelihoods that rely on theoretical assumptions complicated by systematic uncertainties. To reduce the dimensionality of the problem, this method infers macroscopic properties of the star (mass and radius) from emitted X‐ray spectra, and from those quantities infers the EOS. This work approaches the same problem using machine learning techniques, demonstrating a series of enhancements to the current state‐of‐the‐art by realistic uncertainty quantification and reducing the need for theoretical assumptions. We also demonstrate novel inference of the EOS directly from high‐dimensional simulated X‐ray spectra from neutron stars that negate the need for a piece‐wise approach. This inference allows for a natural propagation of uncertainties from the X‐ray spectra by conditioning the discussed networks on realistic sources of uncertainty for each star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Combining Electromagnetic and Gravitational-Wave Constraints on Neutron-Star Masses and Radii

In this work, we perform a joint Bayesian inference of neutron-star mass and radius constraints based on GW170817, observations of quiescent low-mass x-ray binaries (QLMXBs), photospheric radius expansion x-ray bursting sources, and x-ray timing observations of J0030+0451. With this dataset, the form of the prior distribution still has an impact on the posterior mass-radius curves and equation of state (EOS), but this impact is smaller than recently obtained when considering QLMXBs alone. We analyze the consistency of the electromagnetic data by including an “intrinsic scattering” contribution to the uncertainties, and find only a slight broadening of the posteriors. This suggests that the gravitational-wave and electromagnetic observations of neutron-star structure are providing a consistent picture of the neutron-star mass-radius curve and the EOS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

T -matrix analysis of static Wilson line correlators from lattice QCD at finite temperature

Here, we utilize a previously constructed thermodynamic T -matrix approach to the quark-gluon plasma (QGP) to derive constaints on its input by using results on Wilson line correlators (WLCs) of a static quark-antiquark pair from 2 + 1-flavor lattice-QCD (lQCD) computations with realistic pion mass. The self-consistent T -matrix results, which include previous constraints from the lQCD equation of state in the light-parton sector, can describe the lQCD data for WLCs fairly well once refinements of its driving kernel are applied. In particular, the input potential requires less screening than what has been inferred from previous T-matrix analyses. Pertinent predictions for the spectral and transport properties of the QGP are discussed, including the spatial diffusion coefficient for heavy quarks; the latter turns out to have a rather weak temperature dependence, in approximate agreement with recent unquenched lQCD results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A reactive flow model for the 3,3'-diamino-4,4'-azoxyfurazan based plastic bonded explosive (PBX 9701)

Here, this paper describes the calibration of the Arrhenius Wescott–Stewart–Davis (AWSD) reactive flow model for the recently developed high explosive PBX 9701, which consists of 97% 3,3'-diamino-4,4'-azoxyfurazan (DAAF) and 3% FK-800 binder by weight. DAAF-based explosive formulations have several desirable qualities as they are relatively insensitive to non-shock insults but have higher performance than triaminotrinitrobenzene based formulations. Equations of state for the explosive reactants and products are calibrated using a combination of existing and new experimental data and theoretical calculations. The AWSD rate law calibration utilizes both one-dimensional shock-to-detonation and multidimensional rate stick experiments to capture the shock initiation and propagating detonation regimes. Validation of the calibrated model is demonstrated through comparison with recent gas-gun experiments.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2021

Hot and dense quark-gluon plasma thermodynamics from holographic black holes

We present new results on the equation of state and transition line of hot and dense strongly interacting QCD matter, obtained from a bottom-up Einstein-Maxwell-Dilaton holographic model. We considerably expand the previous coverage in baryon densities in this model by implementing new numerical methods to map the holographic black hole solutions onto the QCD phase diagram. We are also able to obtain, for the first time, the first-order phase transition line in a wide region of the phase diagram. Comparisons with the most recent lattice results for the QCD thermodynamics are also presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Verification of a specialized hydrodynamic simulation code for modeling deflagration and detonation of high explosives

A specialized hydrodynamic simulation code has been developed and verified for the simulation of one-dimensional unsteady problems involving the detonation and deflagration of high explosives. To model all the relevant physical processes in these problems, a code is required to simulate compressible hydrodynamics, unsteady thermal conduction, and chemical reactions with complex rate laws. Several verification exercises are presented which test the implementation of these capabilities. The code also requires models for physics processes such as equations of state and conductivity for pure materials and mixtures as well as rate laws for chemical reactions. Additional verification tests are required to ensure that these models are implemented correctly. Though this code is limited in the types of problems it can simulate, its computationally efficient formulation allows it to be used in calibration studies for reactive burn models for high explosives. Furthermore, this study demonstrates how a series of verification tests can be used to ensure that the various physics processes needed to simulate complex phenomenon can be tested to ensure that they are correctly implemented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A modified many-body dissipative particle dynamics model for mesoscopic fluid simulation: methodology, calibration, and application for hydrocarbon and water

The many-body dissipative particle dynamics (mDPD) is a prominent mesoscopic multiphase model for fluid transport in mesoconfinement. However, it has been a long-standing challenge for mDPD (and other multiphase-enabled DPD models) to accurately predict real-fluid static and dynamic properties simultaneously. We have developed a modified mDPD model to overcome the issue and a rigorous calibration approach that uses reference data, including experimental and/or molecular dynamics (MD) simulations to parameterise the modified mDPD for real fluids. We choose heptane as a representative example of hydrocarbon in source rocks to demonstrate the model's capability to accurately predict the equation of state (EOS), free surface tension, diffusivity, and viscosity. Our timing test shows that the modified mDPD is 400–500 times faster than its MD counterpart for simulating bulk heptane in equivalent volumes. To further demonstrate the robustness of the model, we revisited the benchmark problem of mesoscopic modelling of water, in which all the previous works on DPD reported only a limited portion of the water properties. Here, we show that the modified mDPD can provide accurate modelling of water static and dynamic properties and an EOS that matches the experimental data to a large range of confinement pressure.

42 ENGINEERING↗
DOE OSTI2021

Impact of NICER’s Radius Measurement of PSR J0740+6620 on Nuclear Symmetry Energy at Suprasaturation Densities

By directly inverting several neutron star observables in the three-dimensional parameter space for the Equation of State of super-dense neutron-rich nuclear matter, we show that the lower radius limit for PSR J0740 + 6620 of mass 2.08 ± 0.07 M ⊙ from Neutron Star Interior Composition Explorer (NICER)'s very recent observation sets a much tighter lower boundary than previously known for nuclear symmetry energy in the density range of (1.0 ~ 3.0) times the saturation density ρ 0 of nuclear matter. Furthermore, the super-soft symmetry energy leading to the formation of proton polarons in this density region of neutron stars is clearly disfavoured by the first radius measurement for the most massive neutron star observed reliably so far.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Avoiding pitfalls in molecular simulation of vapor sorption: Example of propane and isobutane in metal–organic frameworks for adsorption cooling applications

This study introduces recommendations for conducting molecular simulations of vapor adsorption, with an emphasis on enhancing the accuracy, reproducibility, and comparability of results. The first aspect we address is consistency in the implementation of some details of typical molecular models, including tail corrections and cutoff distances, due to their significant influence on generated data. We highlight the importance of explicitly calculating the saturation pressures at relevant temperatures using methods such as Gibbs ensemble Monte Carlo simulations and illustrate some pitfalls in extrapolating saturation pressures using this method. For grand canonical Monte Carlo (GCMC) simulations, the input fugacity is usually calculated using an equation of state, which often requires the critical parameters of the fluid. We show the importance of using critical parameters derived from the simulation with the same model to ensure internal consistency between the simulated explicit adsorbate phase and the implicit bulk phase in GCMC. We show the advantages of presenting isotherms on a relative pressure scale to facilitate easier comparison among models and with experiment. Extending these guidelines to a practical case study, we evaluate the performance of various isoreticular metal–organic frameworks (MOFs) in adsorption cooling applications. This includes examining the advantages of using propane and isobutane as working fluids and identifying MOFs with a superior performance.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Density of NaAlSi 2 O 6 Melt at High Pressure and Temperature Measured by In-Situ X-ray Microtomography

In this study, the volumetric compression of jadeite (NaAlSi 2 O 6 ) melt at high pressures was determined by three-dimensional volume imaging using the synchrotron-based X-ray microtomography technique in a rotation-anvil device. Combined with the sample mass, measured using a high-precision analytical balance prior to the high-pressure experiment, the density of jadeite melt was obtained at high pressures and high temperatures up to 4.8 GPa and 1955 K. The density data were fitted to a third-order Birch-Murnaghan equation of state, resulting in a best-fit isothermal bulk modulus K T0 of $10.8$ $^{+1.9}_{–5.3}$ GPa and its pressure derivative K' T0 of $3.4$ $^{+6.6}_{–0.4}$. Comparison with data for silicate melts of various compositions from the literature shows that alkali-rich, polymerized melts are generally more compressible than alkali-poor, depolymerized ones. The high compressibility of jadeite melt at high pressures implies that polymerized sodium aluminosilicate melts, if generated by low-degree partial melting of mantle peridotite at ~250–400 km depth in the deep upper mantle, are likely denser than surrounding mantle materials, and thus gravitationally stable.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Shear strength measurements and hydrostatic compression of rhenium diboride under high pressures

Shear strength measurements have been carried out on rhenium diboride, ReB 2 , to a pressure of 74 GPa using a Radial X-ray Diffraction (R-XRD) technique in a diamond anvil cell using platinum as an internal x-ray pressure standard. The R-XRD result has provided a unique insight into the deformation of hexagonal crystal lattice under non-hydrostatic compression and variation of shear strength with increasing pressure. From R-XRD data, we have estimated hydrostatic component of compression to determine an equation of state of rhenium diboride yielding a bulk modulus of K 0 = 366 ± 25 GPa with a pressure derivative $K^{'}_{0}$ = 4.3 ± 0.5 in good agreement with hydrostatic density functional theory calculations. Here, the average lower bound of shear strength (τ) from various diffraction planes was then calculated using the measured interplanar d-spacing (d m ) and hydrostatic component of d-spacing (d p ) to be shown to approach 6.7 ± 0.4 GPa at 70 GPa. Our results show that the anisotropic compression effects observed in ReB 2 under hydrostatic compression are correlated to electronic structure changes under compression as predicted by theoretical calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Photoluminescence of Cr 3+ in β-Ga 2 O 3 and (Al 0.1 Ga 0.9 ) 2 O 3 under pressure

The effects of pressure on single crystals of Cr-doped gallium oxide (β-Ga 2 O 3 :Cr 3+ ) and aluminum gallium oxide [(Al 0.1 Ga 0.9 ) 2 O 3 ] were examined by measuring the wavelength shift in the spectral R lines. Photoluminescence (PL) spectra of these materials were collected from samples in diamond anvil cells at pressures up to 9 GPa. The β-Ga 2 O 3 :Cr 3+ R lines were found to shift linearly under hydrostatic pressure. The (Al 0.1 Ga 0.9 ) 2 O 3 R lines also show a linear shift but the R 1 line shifted less than for β-Ga 2 O 3 :Cr 3+ . The ratio of R 2 to R 1 peak areas vs pressure is dominated by nonradiative recombination. X-ray diffraction measurements of (Al 0.1 Ga 0.9 ) 2 O 3 indicate that its equation of state is similar to that of β-Ga 2 O 3 . β-Ga 2 O 3 :Cr 3+ was examined under non-hydrostatic conditions by using mineral oil as a pressure transmitting medium. Similar to the case in ruby, the R 1 line is much more sensitive to non-hydrostatic stress than R 2 . Spatially resolved PL of a sample at 8 GPa in mineral oil showed significant variations in the R 1 emission wavelength. Furthermore, these results suggest that the R 1 line can serve as a sensitive probe of alloy composition and non-hydrostatic stress, while the R 2 line is insensitive to these perturbations.

36 MATERIALS SCIENCE↗