Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equation-of-state sensitivity”

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.

Cross-Code verification and sensitivity analysis to effectively model the electrothermal instability

This manuscript presents verification cases that are developed to study the electrothermal instability (ETI). Specific verification cases are included to ensure that the unit physics components necessary to model the ETI are accurate, providing a path for fluid-based codes to effectively simulate ETI in the linear and nonlinear growth regimes. Two software frameworks with different algorithmic approaches are compared for accuracy in their ability to simulate diffusion of a magnetic field, linear growth of the ETI, and a fully nonlinear ETI evolution. The nonlinear ETI simulations show early time agreement, with some differences emerging, as noted in the wavenumber spectrum, late into the nonlinear development of ETI. Here, a sensitivity study explores the role of equation-of-state (EOS), vacuum density, and vacuum resistivity. EOS and vacuum resistivity are found to be the most critical factors in the modeling of nonlinear ETI development.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Code Development Supporting a Non-Thermal Source of High Fluence Warm X-Ray

A six-month research effort has advanced the hybrid kinetic-fluid modeling capability required for developing non-thermal warm x-ray sources on Z. The three particle treatments of quasi-neutral, multi-fluid, and kinetic are demonstrated in 1D simulations of an Ar gas puff. The simulations determine required resolutions for the advanced implicit solution techniques and debug hybrid particle treatments with equation-of-state and radiation transport. The kinetic treatment is used in preliminary analysis of the non-Maxwellian nature of a gas target. It is also demonstrates the sensitivity of the cyclotron and collision frequencies in determining the transition from thermal to non-thermal particle populations. Finally, a 2D Ar gas puff simulation of a Z shot demonstrates the readiness to proceed with realistic target configurations. The results put us on a very firm footing to proceed to a full LDRD which includes continued development transition criteria and x-ray yield calculation.

36 MATERIALS SCIENCE↗

Dependence of Cosmological Constraints on Gray Photometric Zero-point Uncertainties of Supernova Surveys

Abstract Type Ia supernova (SN) measurements of the Hubble constant, H 0 ; cosmic mass density, Ω M ; and dark energy equation-of-state parameter, w , rely on heterogeneous SN surveys across three decades of observation. These distinct surveys may have undiagnosed, relative photometric zero-point errors. We determine the sensitivities of the SH0ES+Pantheon+ cosmological constraints to unknown gray systematics in the photometric zero-point calibration between the 19 surveys that comprise the Pantheon+ SN compendium. Varying the surveys’ gray zero-points simultaneously with cosmological parameters, we determine that the SH0ES+Pantheon+ measurement of H 0 is robust against gray intersurvey photometric miscalibration. Specifically, uncalibrated intersurvey systematics could represent a source of uncertainty no larger than 0.39 km s −1 Mpc −1 for H 0 . This modest increase in H 0 uncertainty could not account for the 7 km s −1 Mpc −1 “Hubble tension” between the SH0ES measurement of H 0 and the Planck ΛCDM-based inference of H 0 . However, the SH0ES+Pantheon+ best-fit values of Ω M and w are not robust against gray zero-point error, slipping by up to 0.16 and −0.63. Because measurements of Ω M and w depend on intrasurvey cross-band calibration, the hypothetical gray miscalibration underestimates the sensitivity of these measurements to zero-point miscalibration. Because the Pantheon+ compendium contains many surveys that share low- z Hubble flow (HF) and Cepheid-paired SNe, intersurvey photometric calibration errors do not significantly impede the joint use of SH0ES and Pantheon+ to measure H 0 to 1% accuracy. However, H 0 constraints that rely on one HF survey but numerous galactic distance calibration surveys are susceptible to intersurvey photometric miscalibration.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantitative x ray phase contrast imaging of oblique shock wave–interface interactions

Oblique shock wave–interface interactions of gases and liquids have been extensively studied in shock tubes using optical methods to measure equation-of-state (EOS) parameters. However, this is difficult with solids due to their opaqueness to visible light. X ray phase contrast imaging (XPCI) has the penetrative strength to probe solids while still being sensitive to mass density and enhancing the visibility of material boundaries. We investigate the accuracy and repeatability of measuring the mean value of the average mass density (areal density divided by thickness) over region S (B S ) and flow deflection angle (θ) from XPCI images of a sample. To that end, a Hough transform-based method for measuring θ is developed. To measure B S , the XPCI image intensity probability density function (PDF) is modeled accounting for the spatial distribution of x ray energy, scintillator response, and pulse-to-pulse variation in the x ray intensity. In addition, a Monte Carlo-based algorithm for computing the B S PDF is developed. Both methods are validated on an impact-generated oblique shock wave interacting at a solid polymer-to-polymer interface. Furthermore this is accomplished through a comparison to hydrodynamic simulations using well-established EOS. Under the modeling framework for the XPCI image intensity, B S is computed with an accuracy of < 0.1% and precision of 3%–5%, while θ has an uncertainty of 0.2°, respectively. This shows that the XPCI-based model that is developed here could be an invaluable tool for high-fidelity testing of hydrodynamic models in shock polar configurations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Equation of state at neutron-star densities and beyond from perturbative QCD

We explore the consequences of imposing robust thermodynamic constraints arising from perturbative quantum chromodynamics (QCD) when inferring the dense-matter equation-of-state (EOS). We find that the termination density, up to which the EOS modeling is performed in an inference setup, strongly affects the constraining power of the QCD input. This sensitivity in the constraining power arises from EOSs that have a specific form, with drastic softening immediately above the termination density followed by a strong stiffening. Here we also perform explicit modeling of the EOS down from perturbative-QCD densities to construct a new QCD likelihood function that incorporates additional perturbative-QCD calculations of the sound speed and is insensitive to the termination density, which we make publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

In situ synchrotron diffraction of pressure-induced phase transition in DyPO 4 under variable hydrostaticity

In situ synchrotron x-ray diffraction was conducted on polycrystalline DyPO 4 to elucidate the details of the pressure-induced transition from the xenotime polymorph to the monazite polymorph. We used three different pressure-transmitting media (neon, a 16:3:1 methanol-ethanol-water mixture, and potassium chloride) to investigate the effect of hydrostaticity on the phase behavior. Specifically, our data clearly show a hydrostatic onset pressure of the xenotime-monazite transition of 9.1 GPa, considerably lower than the 15.3 GPa previously determined by Raman spectroscopy. Based on (quasi)hydrostatic data taken in a neon environment, third-order Birch-Murnaghan equation-of-state fits give a xenotime bulk modulus of 144 GPa and a monazite bulk modulus of 180 GPa (both with pressure derivatives of 4.0). Structural data and axial compressibilities show that DyPO 4 is sensitive to shear and has an anisotropic response to pressure. More highly deviatoric conditions cause the onset of the transition to shift to pressures at least as low as 7.0 GPa. We attribute early transition to shear-induced distortion of the PO4 tetrahedra. Finally, our characterization of the high-pressure behavior of DyPO 4 under variable hydrostaticity is critical for advancing rare earth orthophosphate fiber coating applications in ceramic matrix composites and may inform future tailoring of phase composition for controlled shear and pressure applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Constraints on the Maximum Densities of Neutron Stars from Postmerger Gravitational Waves with Third-Generation Observations

Using data from 289 numerical relativity simulations of binary neutron star mergers, we identify, for the first time, a robust quasiuniversal relation connecting the postmerger peak gravitational-wave frequency and the value of the density at the center of the maximum mass nonrotating neutron star. This relation offers a new possibility for precision equation-of-state constraints with next-generation ground-based gravitational-wave interferometers. Mock Einstein Telescope observations of fiducial events indicate that Bayesian inferences can constrain the maximum density to ~15% (90% credibility level) for a single signal at the minimum sensitivity threshold for a detection. If the postmerger signal is included in a full-spectrum (inspiral-merger-postmerger) analysis of such a signal, the pressure-density function can be tightly constrained up to the maximum density, and the maximum neutron star mass can be measured with an accuracy better than 12% (90% credibility level).

79 ASTRONOMY AND ASTROPHYSICS↗

Quantifying uncertainty in analysis of shockless dynamic compression experiments on platinum. II. Bayesian model calibration

Dynamic shockless compression experiments provide the ability to explore material behavior at extreme pressures but relatively low temperatures. Typically, the data from these types of experiments are interpreted through an analytic method called Lagrangian analysis. Here, in this work, alternative analysis methods are explored using modern statistical methods. Specifically, Bayesian model calibration is applied to a new set of platinum data shocklessly compressed to 570 GPa. Several platinum equation-of-state models are evaluated, including traditional parametric forms as well as a novel non-parametric model concept. The results are compared to those in Paper I obtained by inverse Lagrangian analysis. The comparisons suggest that Bayesian calibration is not only a viable framework for precise quantification of the compression path, but also reveals insights pertaining to trade-offs surrounding model form selection, sensitivities of the relevant experimental uncertainties, and assumptions and limitations within Lagrangian analysis. The non-parametric model method, in particular, is found to give precise unbiased results and is expected to be useful over a wide range of applications. The calibration results in estimates of the platinum principal isentrope over the full range of experimental pressures to a standard error of 1.6%, which extends the results from Paper I while maintaining the high precision required for the platinum pressure standard.

Brown, Justin Lee↗

A platform to measure isentropes from proton-heated warm dense matter on short pulse laser facilities

We describe the development of an experimental platform that measures the release isentrope of materials heated isochorically to temperatures of a few electron volts, using short-pulse laser-produced protons to heat the sample and long-pulse laser-produced x rays to perform streaked x-ray radiography. The density profiles derived from the radiography data are integrated to generate pressure–density isentropes, independent of prior knowledge of the equation of state of the sample material. In order to understand the sensitivities of isentrope extraction from radiography data, we analyze synthetic radiographs generated by a radiation hydrodynamics code. Noise reduction and high spatial resolution are critical for isentrope reconstruction, as demonstrated by the analysis of a proof-of-principle shot day on the OMEGA-EP facility. In conclusion, the data demonstrate the feasibility of the platform for characterizing isentropes, and we discuss the necessary improvements to enhance precision in differentiating between equation-of-state models.

Equations of state↗

Diagnosing the DARHT Electron Beam-Target Interaction and Hydrodynamic Expansion [Dissertation]

Quantitative electron beam-target interaction studies on the Dual-Axis Radiographic Hydrodynamic Test (DARHT) electron linear induction accelerators have only been performed recently. This includes characterization of the temperature, density, pressure, extent, and expansion velocity of the plasma plume. The results presented in this dissertation present a detailed and unified overview of the accelerator systems, target heating physics, beam transport, and diagnostic tools. Additional information includes calibration sources, radiation hydrodynamics, spectroscopic-quality radiation transport modeling, experimental measurements, and analyses of electron beam driven aluminum experiments. The first set of spatially and temporally resolved spectroscopic measurements of electron beam driven aluminum are presented. Contamination quantification analyses are used to understand the origin of the strong Na-I 3p-3s lines that are observed in absorption within the aluminum plasma continuum. These results inform the creation of the first spectroscopic-quality radiation transport model that links several atomic physics codes to interpret the conditions from which the Na-I lines originate. A good agreement is found between the surface analysis results and the model which confirms the concentration of the sodium present within the aluminum alloy foil material. It also demonstrates, for the first time on electron beam-driven target experiments, the ability to interpret plasma conditions from measured absorption lines. In a second experimental campaign which focuses on pure aluminum, the Al-I 3p-4s and 3p3d doublets are both measured in emission. A detailed analysis of the Al-I 3p-4s doublet reveals that the lines undergo moderate self-absorption. A simple model of the self-absorption effect is successfully used to match measured spectra at various temperature/density/plasma scale length combinations. These measurements led to the realization of the minimum density that can be resolved by the spectrometers for the Al-I 3p-4s lines due to the large slit width required to observe a signal on aluminum. These measurements also demonstrate the sensitivity of visible and long wave UV spectroscopy to minor changes in both temperature and density. The simple self-absorption model will be useful for analysis of other beam-target interaction experiments with spectra exhibiting either self-absorption or full self-reversal. Substantial headway has been made on the modeling front by linking together several codes needed to model both the energy deposition, hydrodynamic motion, and atomic kinetics to produce synthetic spectral calculations that are compared with experimental measurements. There exists ample space for improvement, especially with benchmarking the hydrodynamics codes and equation-of-state tables with experimental measurements. The X-ray diagnostics required to make these new measurements along with the simulation capabilities required to interpret the results are under development and will be the subject of future studies.

43 PARTICLE ACCELERATORS↗

Laser-driven ramp-compression experiments on the national ignition facility

This report details the analyses and related uncertainties in measuring longitudinal-stress–density paths in indirect laser-driven ramp equation-of-state (EOS) experiments. Experiments were conducted at the National Ignition Facility (NIF) located at the Lawrence Livermore National Laboratory. The NIF can deliver up to 2 MJ of laser energy over 30 ns and provide the necessary laser power and control to ramp compress materials to TPa pressures (1 TPa = 10 × 10 6 atmospheres). These data provide low-temperature solid-state EOS data relevant to the extreme conditions found in the deep interiors of giant planets. In these experiments, multi-stepped samples with thicknesses in the range of 40–120 µm experience an initial shock compression followed by a time-dependent ramp compression to peak pressure. Interface velocity measurements from each thickness combine to place a constraint on the Lagrangian sound speed as a function of particle velocity, which in turn allows for the determination of a continuous stress–density path to high levels of compressibility. In this report, we present a detailed description of the experimental techniques and measurement uncertainties and describe how these uncertainties combine to place a final uncertainty in both stress and density. Here, we address the effects of time-dependent deformation and the sensitivity of ramp EOS techniques to the onset of phase transformations.

47 OTHER INSTRUMENTATION↗

Sensitivity of an integrated experiment to uncertainty in the high explosive equations of state

Traditionally, hydrodynamics simulations are performed with a single equation-of-state (EOS) to describe each material. These EOSs typically have a physics-informed functional form with adjustable parameters that are calibrated in order to replicate small-scale data. However, because the calibration data have uncertainty and there are typically inherent degeneracies in fitting the EOS, there are actually multiple EOSs that might be consistent with calibration data. In this work, we perform uncertainty quantification (UQ) for the reactant and product equations of state for the high explosive PBX 9501 to yield an ensemble of EOSs that match the uncertain small-scale calibration data. We then simulate an experiment of an explosively formed penetrator repeatedly with different EOSs to both validate the UQ analysis and determine the effects of EOS uncertainty on the prediction of quantities of interest in the experiment. In general, we find good agreement between the simulation predictions and the experimental measurements, and we identify an EOS variable that contributes most directly to the spread in the predictions as the EOSs are varied.

36 MATERIALS SCIENCE↗