Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equations 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 1,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Fluctuations in heavy ion collisions and global conservation effects

Subensemble is a type of statistical ensemble which is the generalization of grand canonical and canonical ensembles. The subensemble acceptance method (SAM) provides general formulas to correct the cumulants of distributions in heavy-ion collisions for the global conservation of all QCD charges. The method is applicable for an arbitrary equation of state and sufficiently large systems, such as those created in central collisions of heavy ions. The new fluctuation measures insensitive to global conservation effects are presented. The main results are illustrated in the hadron resonance gas and van der Waals fluid frameworks.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Analysis of static Wilson line correlators from lattice QCD at finite temperature with T -matrix approach

The thermodynamic T-matrix approach is used to study Wilson line correlators (WLCs) for a static quark-antiquark pair in the quark-gluon plasma (QGP). Selfconsistent results that incorporate constraints from the QGP equation of state can approximately reproduce WLCs computed in 2+1-flavor lattice-QCD (lQCD), provided the input potential exhibits less screening than in previous studies. Utilizing the updated potential to calculate pertinent heavylight T-matrices we evaluate thermal relaxation rates of heavy quarks in the QGP. We find a more pronounced temperature dependence for low-momentum quarks than in our previous results (with larger screening), which turns into a weaker temperature dependence of the (temperature-scaled) spatial diffusion coefficient, in fair agreement with the most recent lQCD data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

QCD material parameters at zero and non-zero chemical potential from the lattice

Using an eighth-order Taylor expansion in baryon-chemical potential, we recently obtained the (2+1) flavor QCD equation of state (EoS) at non zero conserved charge chemical potentials from the lattice. We focused on strangeness-neutral, isospin-symmetric QCD matter, which closely resembles the situation encountered in heavy-ion collision experiments. Using this EoS, we present here results on various QCD material parameters; in particular we compute the specific heat, speed of sound, and compressibility along appropriate lines of constant physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Polysulfone shock compressed above the decomposition threshold: Velocimetry and modeling of two-wave structures

Polysulfone was shock loaded to pressures of 14.6–26.2 GPa in a series of gas gun-driven plate-impact experiments measuring material response with embedded electro-magnetic particle velocity gauges and optical velocimetry. The embedded electro-magnetic particle velocity gauges did not show a distinct two-wave structure but did show rounding that suggested a reaction but not a distinct separation of the reactants and product waves. In contrast, the transmission experiments fielded with optical velocimetry, with product pressures ranging from 21.0 to 26.2 GPa, showed well-defined two-wave structures due to shock-driven chemical decomposition of the polymer to products at a higher density. Distinct two-wave structures have not previously been published in literature during polymer compression; here, we observed these two-wave structures at both the polymer/lithium fluoride and polymer/polymethylmethacrylate interfaces. Hydrodynamic simulations of the experiments were performed using a newly constructed SESAME equation of state (EOS) for the polymer, a thermochemical EOS for the decomposition products, and an Arrhenius reaction rate model for the kinetics of the reaction. Finally, simulation results also demonstrated two-wave structures but were unable to quantitatively reproduce either the embedded gauge or optical velocimetry data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Pulsed power driven cylindrical wire array explosions in different media

Cylindrical copper wire array explosions were carried out in de-ionized water, sodium polytungstate solution, nitromethane, and polyester in order to obtain high energy density conditions in the vicinity of implosion using the generated converging shock waves. The utilization of different materials in which the array is immersed can contribute to this goal with higher density resulting in higher shock velocities and possible combustion. The generated shock waves were captured by a framing and a streak camera and shock velocities were calculated and compared. Pressure behind the shock front was calculated using known hydrodynamic relations (for water, polytungstate, and polyester) and compared to two-dimensional hydrodynamic simulations coupled with equations of state (for water and polyester). It was demonstrated that despite lower shock wave velocity in polytungstate solution than water, the pressures generated are similar in both materials. In polyester, both shock velocities and generated pressures are 2-4 times higher than in water. Furthermore, it was shown that it is possible to carry out these explosions in a solid which has several advantages compared to liquids, such as not relying on waterproof systems and easier transportation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Spatially varying colloidal phase behavior on multi-dimensional energy landscapes

A method is reported to determine equilibrium concentration profiles and local phase behavior for colloids on multi-dimensional energy landscapes. A general expression is derived based on local particle concentration and osmotic pressure differences that are balanced by forces on colloids due to energy landscape gradients. This analysis is applied to colloidal particles in high frequency AC electric fields within octupolar electrodes, where the energy landscape can be shaped in two dimensions. These results are also directly applicable to any particles having induced dipoles in spatially non-uniform electromagnetic fields. Predictions based on modelling colloids with an effective hard disk equation of state indicate inhomogeneous solid and fluid states coexisting on different shaped energy landscapes including multiple minima. Model predictions show excellent agreement with time-averaged Brownian Dynamic simulations at equilibrium. Findings demonstrate a general approach to understand colloidal phase behavior on energy landscapes due to external fields, which could enable control of colloidal microstructure on morphing energy landscapes and the inverse design of fields to assemble hierarchically structured colloidal materials.

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

Single Crystal Plasticity Model with Deformation Twinning for the High Rate Deformation of β-HMX

Deformation twinning is an active deformation mechanism during loading along certain crystallographic directions of β-cyclotetramethylene tetranitramine (β-HMX). In this work, a finite strain thermomechanical model is extended to include twinning as a deformation mechanism in addition to plastic slip. The stress is derived from the free energy expression including a term representing the equation of state. The crystal plasticity framework is used to divide the total strain into inelastic and elastic component, where the elastic part is used in the expression for the free energy. The twin systems are treated as pseudo slip systems and the shear rate on the twin systems is evaluated in terms of the twin resistance and the projection of the stress tensor. The model parameters were calibrated against a set of plate impact experiments performed on β-HMX by Dick et al. The remaining plate impact experiments are used to evaluate the predictive capability of the model. The quality of the model fits and predictions is discussed from a physical and modeling perspective. Particularly, the role of twin modeling on results is highlighted. Finally, the model is used to explore the relationship between the propensity for twinning and crystal orientation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Pressure-induced transformation of Nb 2 O 5 under shock compression from first principles

Ab initio molecular dynamics (AIMD) simulations were carried out to investigate the equation of state of Nb 2 O 5 and its pressure-density relationship under shock conditions. Here, the focus of this study is on the monoclinic B-Nb 2 O 5 (C2/c) polymorph. Enthalpy calculations from AIMD trajectories at 300 K show that the pressure-induced transformation between the thermody-namically most stable crystalline monoclinic parent phase H-Nb 2 O 5 (P2/m) and B-Nb 2 O 5 occurs at ~1.9 GPa. This H→B transition is energetically more favorable than the H→L(Pmm2) pressure-induced transition recently observed at ~5.9-9.0 GPa. The predicted shock properties of Nb 2 O 5 polymorphs are also compared to their Nb and NbO 2 counterparts to assess the impact of niobium oxidation on shock response.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Structural disorder, sublattice melting, and thermo-elastic properties of anti-perovskite Li 3 OBr under high pressure and temperature

Here, in situ synchrotron x-ray diffraction experiments were conducted in the system Li-O-Br at pressures up to 6.5 GPa and temperatures up to the melting points. The thermal equation of state determined from the P–V–T measurements indicates that crystalline anti-perovskite Li 3 OBr exhibits similar thermo-elastic properties to common salts. At temperatures that are 50–75 °C below the crystal melting, a pre-melting behavior was inferred based on weakened Bragg intensity and a substantial increase in the diffuse scattering. Concurrent with the onset of pre-melting, an anomalous increase in the lattice volume was observed and is interpreted as anion sublattice disorder. These findings support previous ab initio molecular dynamics simulations based on a simple anion interchange mechanism, indicating that structural disorder below crystal melting temperature can be a main driving force for Li-sublattice melting and superionic transition in the Li 3 OBr anti-perovskite.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Vibrational spectroscopy of shock compressed condensed phase nitrous oxide: Frequency shifts in the ν 1 mode

Over many years, we have published vibrational spectra of molecular materials at extreme pressure and temperature obtained using shock compression and coherent Raman spectroscopy methods. For diatomic molecules, we were able to extract vibrational temperatures from the intensities of spectrally resolved vibrational hot bands. Larger molecules in the condensed phase suffer band broadening effects that obscure the vibrational hot bands as the anharmonicities are typically smaller than the widths of the bands. This inability to resolve the hot bands inhibits the extraction of vibrational temperature and the measurement of the vibrational frequency of the fundamental. Here, we use a hot band model based on gas-phase anharmonic coupling coefficients to fit coherent anti-Stokes Raman spectra of the ν 1 vibrational mode of shock compressed condensed phase N 2 O with shock pressures and temperatures estimated from literature equations of state and compare to fits from a model using a single Gaussian peak. We report the resulting vibrational frequency shifts with shock pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Near-field models and simulations of pellet ablation in tokamaks

Here, detailed numerical studies of the ablation of a single neon pellet in the plasma disruption mitigation parameter space have been performed. Simulations were carried out using FronTier, a hydrodynamic/low magnetic Reynolds number MHD code with explicit tracking of material interfaces. FronTier's physics models resolve the pellet surface ablation and the formation of a dense, cold cloud of ablated material, the deposition of energy from hot plasma electrons, the expansion of the ablation cloud along magnetic field lines, and the radiation losses. A local thermodynamic equilibrium model based on Saha equations has been used to resolve atomic processes in the cloud and Redlich–Kwong corrections to the ideal gas equation of state for cold and dense gases, which have been used near the pellet surface. The FronTier pellet code is a next generation of the code described. It has been validated against the semi-analytic improved neutral gas shielding model in the 1D spherically symmetric approximation. Main results include quantification of the influence of atomic processes and Redlich–Kwong corrections on the pellet ablation in spherically symmetric approximation and verification of analytic scaling laws in a broad range of pellet and plasma parameters. Using axially symmetric MHD simulations, the properties of ablation channels and the reduction of pellet ablation rates in magnetic fields of increasing strength have been studied. While the main emphasis has been given to neon pellets for the plasma disruption mitigation, selected results on deuterium fueling pellets have also been presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The ultrahigh pressure stability of silver: An experimental and theoretical study

We measured the atomic volume of Ag in a toroidal diamond anvil cell to a maximum pressure of 416 GPa and calculated the atomic volume and elastic constants of Ag up to 750 and 460 GPa, respectively. Our density functional theory calculations at 0 K utilize an all-electron fully relativistic method and agree well with our volume measurements, particularly at pressures above ~75 GPa. In this work, we corrected our experimental results for non-hydrostaticity using a line shift analysis, and the resulting Vinet equation of state (EOS) parameters are reported. We find that the uniaxial stress sustained by Ag increases linearly up to 4.5 GPa at a pressure of 416 GPa. Our experimental results indicate that the fcc structure of Ag remains stable to at least 416 GPa at room temperature. Our theoretical results show that C44 increases as pressure increases and reaches a maximum at ~100 GPa above which it begins to decrease, a sign that the fcc structure of Ag is becoming unstable, and at V/V 0 = 0.30, the bcc structure is lower in energy than fcc.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Shock to detonation transition of pentaerythritol tetranitrate (PETN) initially pressed to 1.65 g/cm 3

A novel set of experiments and reactive flow modeling of pentaerythritol tetranitrate (PETN) is presented. In this work, the specific phenomenon of shock to detonation transition is examined, where an initial, relatively weak shock is propagated into pressed PETN powder at 1.65 g/cm 3 and the subsequent buildup to detonation is observed experimentally. These experiments, in conjunction with reactant and products’ equations of state, are utilized for building reactive flow models.

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

Simultaneous inference of the compressibility and inelastic response of tantalum under extreme loading

We study the deformation of tantalum under extreme loading conditions. Experimental velocity data are drawn from both ramp loading experiments on Sandia’s Z-machine and gas gun compression experiments. The drive conditions enable the study of materials under pressures greater than 100 GPa. We provide a detailed forward model of the experiments including a model of the magnetic drive for the Z-machine. Utilizing these experiments, we simultaneously infer several different types of physically motivated parameters describing equation of state, plasticity, and anelasticity via the computational device of Bayesian model calibration. Characteristics of the resulting calculated posterior distributions illustrate relationships among the parameters of interest via the degree of cross correlation. The calibrated velocity traces display good agreement with the experiments up to experimental uncertainty as well as improvement over previous calibrations. Examining the Z-shots and gun-shots together and separately reveals a trade-off between accuracy and transferability across different experimental conditions. Implications for model calibration, limitations from model form, and suggestions for improvements are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Code-to-code comparison and validation of the radiation-hydrodynamics capabilities of the FLASH code using a laboratory astrophysical jet

The potential for laser-produced plasmas to yield fundamental insights into high energy density physics (HEDP) and deliver other useful applications can sometimes be frustrated by uncertainties in modeling the properties and behavior of these plasmas using radiation-hydrodynamics codes. In an effort to overcome this and to corroborate the accuracy of the HEDP capabilities in the publicly available FLASH radiation-hydrodynamics code, we present detailed code-to-code comparisons between FLASH and the HYDRA code developed at Lawrence Livermore National Laboratory using previously published HYDRA simulations from Grava et al. [Phys. Rev. E 78, 016403 (2008)]. That study describes a laser experiment that produced a jet-like feature that the authors compare to astrophysical jets. Importantly, the Grava et al. [Phys. Rev. E 78, 016403 (2008)] experiment included detailed x-ray interferometric measurements of electron number densities and a time-integrated measurement of the soft x-ray spectrum. Despite markedly different methods for treating the computational mesh, and different equations of state and opacity models, the FLASH results resemble the results from HYDRA and, most importantly, the experimental measurements of electron density. Having validated the FLASH code in this way, we use the code to further investigate and understand the formation of the jet seen in the Grava et al. [Phys. Rev. E 78, 016403 (2008)] experiment and discuss its relation to the Wan et al. [Phys. Rev. E 55, 6293 (1997)] experiment at the NOVA laser.

Orban, Chris (ORCID:0000000312004538)↗
DOE OSTI2022

The role of detonation condensates on the performance of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) detonation

Thermochemical models of detonation are widely used to estimate energy delivery, but they are based on the assumption that the carbon-rich condensates (soot) formed during detonation are very similar to bulk carbon. Here, we present an analytic equation of state (EOS) based on experimental detonation data for soot formed during the detonation of triaminotrinitrobenzene (TATB)-based high explosives. X-ray photoelectron spectra of several detonation soots are used to determine the elemental nitrogen abundance, with surprisingly high values for TATB. The proposed TATB soot EOS is highly compressible at low pressures and shares some features of glassy carbon, exhibiting graphite- and diamond-like behavior as a function of pressure. We demonstrate the influence of formed soot on detonation performance, including a lowering of the detonation velocity at typical charge densities, and a more compressive product Hugoniot at overdriven conditions. The soot model improves the accuracy of thermochemical calculations for TATB-based explosives across a wide range of states. Detonation velocity predictions for HMX (cyclotetramethylene-tetranitramine)-TATB blends with 80% or more TATB content, as well predictions for 1,3-diamino-2,4,6-trinitrobenzene (DATB) and 3-nitro-1,2,4-triazol-5-one (NTO), which share some features with TATB, are also improved.

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

Explosive dispersal of particles in high speed environments

In this paper, we present the results of the explosive dispersal of particles in high-speed environments. We carry out Euler–Lagrange numerical simulations of a source at quiescent ambient conditions as well as moving at Mach numbers of 3 and 6. Particle volume fractions of 0%, 1%, and 4.5% are presented. The detonation profile is computed with the Jones–Wilkins–Lee equation of state using a reactive burn model. Non-static cases provide a framework to consider the effect of a bow shock and pre-existing high-speed flow conditions on the dispersal process. We also compute averages of both static and dynamic pressures, as well as impulse density histories on virtual probe planes to characterize the momentum of the flow and particles that would deposit on a target. Results suggest that the presence of the particles can have a substantial effect on the pressure average of the virtual target planes.

Durant, B. A.↗