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,009 records · Page 56

Shock Hugoniot measurements of single-crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) compressed to 83 GPa

We present laser-driven shock Hugoniot measurements of single-crystal (SC) 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) between 15 and 83 GPa, spanning pressures below and well above the Chapman–Jouguet pressure of ~28 GPa for TATB formulations (TATB grains mixed with plastic binders at 5–10 wt. %). The new SC data are generally ~3% more compressible than previously published data on neat and formulated TATB measured in gas-gun and explosive-driven experiments. An exception is at compressions in the density of ~1.5 (~30–40 GPa), where our new SC data exhibit significantly lower pressures than previous results on overdriven TATB formulations, suggesting that our SC samples remain largely unreacted below 35 GPa over the short nanosecond-time scales inherent to our laser-driven experiments. We report these novel equation-of-state measurements are a critical step toward understanding TATB in its most fundamental form and improving predictive modeling of TATB-based explosives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ethane and methane at high pressures: Structure and stability

We have performed a combined experimental and theoretical study of ethane and methane at high pressures of up to 120 GPa at 300 K using x-ray diffraction and Raman spectroscopies and the USPEX ab initio evolutionary structural search algorithm, respectively. For ethane, we have determined the crystallization point, for room temperature, at 2.7 GPa and also the low pressure crystal structure (phase A). This crystal structure is orientationally disordered (plastic phase) and deviates from the known crystal structures for ethane at low temperatures. Moreover, a pressure induced phase transition has been identified, for the first time, at 13.6 GPa to a monoclinic phase B, the structure of which is solved based on good agreement with the experimental results and theoretical predictions. For methane, our x-ray diffraction measurements are in agreement with the previously reported high-pressure structures and equation of state (EOS). Finally, we have determined the EOSs of ethane and methane, which provides a solid basis for the discussion of their relative stability at high pressures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The behavior of methane–water mixtures under elevated pressures from simulations using many-body potentials

Non-polarizable empirical potentials have been proven to be incapable of capturing the mixing of methane–water mixtures at elevated pressures. Although density functional theory-based ab initio simulations may circumvent this discrepancy, they are limited in terms of the relevant time and length scales associated with mixing phenomena. Here, we show that the many-body MB-nrg potential, designed to reproduce methane–water interactions with coupled cluster accuracy, successfully captures this phenomenon up to 3 GPa and 500 K with varying methane concentrations. Two-phase simulations and long time scales that are required to fully capture the mixing, affordable due to the speed and accuracy of the MBX software, are assessed. Constructing the methane–water equation of state across the phase diagram shows that the stable mixtures are denser than the sum of their parts at a given pressure and temperature. As a result, we find that many-body polarization plays a central role, enhancing the induced dipole moments of methane by 0.20 D during mixing under pressure. Overall, the mixed system adopts a denser state, which involves a significant enthalpic driving force as elucidated by a systematic many-body energy decomposition analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

How macroscopic limits on neutron-star baryon loss yield microscopic limits on non-standard-model baryon decay

We investigate how our baryon-loss limits from anomalous binary-pulsar period lengthening can be interpreted microscopically to yield specific constraints on the particle physics of baryon number violation within a neutron star. We focus on the possibility of anomalous baryon disappearance via dark baryon processes and on scenarios in which the produced dark-sector particles do not survive to influence the response of the star to baryon-number-violating effects. We flesh out the conditions for which this may occur, as well as other key assumptions. We then turn to the analysis of particle processes in the dense nuclear medium found at the core of a neutron star, employing the techniques of relativistic mean-field theory. Using our study of in-medium effects and limits on macroscopic baryon number violation, we extract limits on in-vacuum baryon-number-violating processes, and we determine them for various equations of state. We conclude by noting the implications of our results for models of dark-sector-enabled baryogenesis. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Mirror neutron stars

The fundamental nature of dark matter is entirely unknown. A compelling candidate is twin Higgs mirror matter, invisible hidden-sector cousins of the Standard Model particles and forces. This predicts mirror neutron stars made entirely of mirror nuclear matter. Here we find their structure using realistic equations of state, robustly modified based on first-principle quantum chromodynamic calculations, for the first time. This allows us to predict their gravitational wave signals, demonstrating an impressive discovery potential and ability to probe dark sectors connected to the hierarchy problem.

79 ASTRONOMY AND ASTROPHYSICS↗

Direct Monte Carlo Simulations of Chemical Equilibrium Composition of Molecular Fluid Mixtures Under Shock Conditions

A new Monte Carlo simulation method is implemented which holds the number and kind of atoms constant rather than the number and kind of molecules for a molecular fluid mixture. Chemical reactions are allowed as a correlated move of atoms such as an interchange of atoms between molecules. Equilibrium chemical composition is determined as an average over the simulation along with equation of state quantities. Results are presented for N 2 + O 2 → 2NO at 30GPa and 3000K. The equilibrium composition is shown to be sensitive to the potentials between unlike species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spall strength of additively repaired 304L stainless steel

Additive manufacturing has the potential to repair damaged parts, but the performance of additive materials under high strain rate loading is still uncertain—especially with the added complexity of an interface with an existing wrought material. In this work, 304L stainless steel samples were intentionally damaged and then repaired with wire-fed laser additive manufacturing. The samples were subjected to shock loading to generate incipient spall. Velocimetry and post-mortem metallography results show that when the additive repair process parameters are optimized to reduce porosity and match the equation of state of the original material, the influence of the repair region on the shock propagation is negligible. In conclusion, the free-surface velocity profile and internal damage morphology of the repaired sample are shown to be practically identical to the pristine material.

36 MATERIALS SCIENCE↗

Nanolamellar phase transition in an additively manufactured eutectic high-entropy alloy under high pressures

Much is unknown about how phase transitions link to micro-/nano-structures in high-entropy systems, especially under extreme pressure and temperature conditions. This work studies the evolution of dual-phase nanolamellar eutectic high-entropy alloy phases of AlCoCrFeNi 2.1 generated by laser powder-bed fusion (L-PBF) for pressures up to 42 GPa. We compare quasi-hydrostatic high pressure synchrotron x-ray diffraction studies on L-PBF printed cylindrical samples up to 5.5 GPa (large-volume Paris–Edinburgh cell) to those carried out on an L-PBF printed foil in a diamond anvil cell where the pressure reached 42 GPa. Our results show that the initially alternating face-centered cubic (FCC) and body-centered cubic (BCC) nanolamellar structure of AlCoCrFeNi 2.1 transformed into single-phase FCC nanolamellae under high pressure with BCC–FCC phase transformation completion at 21 ± 3 GPa. Our results indicate a diffusionless BCC–FCC transformation in this additively manufactured far-from-equilibrium microstructure and demonstrate that the FCC phase is stable up to very high pressures. The measured equation of state for the FCC phase of AlCoCrFeNi 2.1 is presented up to 42 GPa and shows excellent agreement between the data obtained in large-volume press and diamond anvil cell experiments.

36 MATERIALS SCIENCE↗

Nonequilibrium information entropy approach to ternary fission of actinides

Ternary fission of actinides probes the state of the nucleus at scission. Light clusters are produced in space and time very close to the scission point. Within the nonequilibrium statistical operator method, a generalized Gibbs distribution is constructed from the information given by the observed yields of isotopes. Using this relevant statistical operator, yields are calculated taking excited states and continuum correlations into account, in accordance with the virial expansion of the equation of state. Furthermore, clusters with mass number A≤10 are well described using the nonequilibrium generalizations of temperature and chemical potentials. Improving the virial expansion, in-medium effects may become of importance in determining the contribution of weakly bound states and continuum correlations to the intrinsic partition function. Yields of larger clusters, which fail to reach this quasiequilibrium form of the relevant distribution, are described by nucleation kinetics, and a saddle-to-scission relaxation time of about 7000 fm/c is inferred. Light-charged particle emission, described by reaction kinetics and virial expansions, may therefore be regarded as a very important tool to probe the nonequilibrium time evolution of actinide nuclei during fission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Thermochemical nonequilibrium in atomic hydrogen at elevated temperatures

A numerical study of the nonequilibrium flow of atomic hydrogen in a cascade arc was performed to obtain insight into the physics of the hydrogen cascade arc. A rigorous mathematical model of the flow problem was formulated, incorporating the important nonequilibrium transport phenomena and atomic processes which occur in atomic hydrogen. Realistic boundary conditions, including consideration of the wall electrostatic sheath phenomenon, were included in the model. The governing equations of the asymptotic region of the cascade arc were obtained by writing conservation of mass and energy equations for the electron subgas, an energy conservation equation for heavy particles and an equation of state. Finite-difference operators for variable grid spacing were applied to the governing equations and the resulting system of strongly coupled, stiff equations were solved numerically by the Newton-Raphson method.

Scott, R. K.↗

The Scaling of Coronal Models from One Star to Another

The requirements that must be met in order that stationary numerical corona models can be scaled from one star to another are discussed. A corona model is a solution of the conservation equations for mass, momentum, and energy, subject to appropriate boundary conditions, and of the equation of state. In general, the mass M and radius R of the star enter these equations and boundary conditions as free parameters. A given solution can be scaled to other stars only if all equations can be rewritten in such a form that M and R do no longer appear explicitly as free parameters, but only implicitly as scaling factors of the variables. An adequate means to find these scaling factors is a homologous transformation: one multiplies all variables and parameters by separate constants (i.e., scaling factors) and requires that the equations and boundary conditions remain valid. This leads to a set of nonlinear relations between the transformation constants. Only if in this set the two constants associated with M and R can be chosen independently, can a given numerical corona model be scaled to arbitrary stars.

Hammer, R.↗

Convergence to steady state of solutions of Burgers' equation

Consider the initial boundary value problem for Burgers' equation. It is shown that its solutions converge, in time, to a unique steady state. The speed of the convergence depends on the boundary conditions and can be exponentially slow. Methods to speed up the rate of convergence are also discussed.

Kreiss, G.↗

The thermal catastrophe model of substorms

Equations are presented for the physical mechanisms involved in the resonant absorption of Alfven waves in the plasma sheet boundary layer. It is shown that energy absorbed by the plasma-sheet particles is a function of the central plasma sheet temperature. The heating curve, when coupled with convective transport, yielded an equation of state for the steady state plasma sheet, whose solution has the form of a mathematical catastrophe. The master equation includes dynamic terms describing the transition across the thermal catastrophe, making it possible to evaluate the time scale for the catastrophe to occur.

Goertz, C. K.↗

Calibration and validation of the foundation for a multiphase strength model for tin

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 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. 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 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 a fully liquid state, never fully lost strength, stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.

36 MATERIALS SCIENCE↗

Influence of heavy resonances in SMASH

Recent lattice QCD results, comparing to a hadron resonance gas model, have shown the need for hundreds of particles in hadronic models. These extra particles influence both the equation of state and hadronic interactions within hadron transport models. Here, we introduce the PDG21+ particle list, which contains the most up-to-date database of particles and their properties. We then convert all particles decays into 2 body decays so that they are compatible with SMASH in order to produce a more consistent description of a heavy-ion collision.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A robust three-parameter reference curve for condensed phase materials

Many isothermal equations of state lead to numerical problems under tension conditions due to negative bulk modulus. We summarize the properties of four existing models and devise a continuous non-piecewise reference curve to address this issue. The final functional form contains three fitting parameters that are directly related to the bulk modulus and its pressure derivative at the reference state and at infinite pressures. We also include bounds that guarantee the monotonicity of its pressure derivative in all volume space. To demonstrate the qualitative features of each isotherm, we conduct calibrations against experimental data for cyclotetramethylene tetranitramine, magnesium oxide, gold, and platinum. Furthermore, pressure residuals and root mean square deviations show that the proposed functional form is capable of adequately fitting the data compared to other three-parameter models.

36 MATERIALS SCIENCE↗

Comparing temperature convergence of shocked thin films of tin and iron to a bulk temperature source

An outstanding challenge in developing a complete equation of state for materials at elevated pressure and temperature is a robust method of determining the bulk temperature state under dynamic conditions. In metals, the determination of bulk temperature states by optical pyrometry is complicated by the small optical depth and thermal conduction effects. These effects lead to observed temperatures differing by 20% or more from the bulk temperature state. In this work, we show the presence of thermal conduction effects in temperature measurements of tin and iron coatings during dynamic compression experiments. We demonstrate that tin, in contrast to iron, coatings can fail to converge to a bulk temperature source over the time scale of the experiment, requiring the experimenter to modify assumptions, design, or analysis. This work bounds thermal transport at shocked conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Direct imaging of shock wave splitting in diamond at Mbar pressure

Understanding the behavior of matter at extreme pressures of the order of a megabar (Mbar) is essential to gain insight into various physical phenomena at macroscales—the formation of planets, young stars, and the cores of super-Earths, and at microscales—damage to ceramic materials and high-pressure plastic transformation and phase transitions in solids. Under dynamic compression of solids up to Mbar pressures, even a solid with high strength exhibits plastic properties, causing the induced shock wave to split in two: an elastic precursor and a plastic shock wave. This phenomenon is described by theoretical models based on indirect measurements of material response. The advent of x-ray free-electron lasers (XFELs) has made it possible to use their ultrashort pulses for direct observations of the propagation of shock waves in solid materials by the method of phase-contrast radiography. However, there is still a lack of comprehensive data for verification of theoretical models of different solids. Here, we present the results of an experiment in which the evolution of the coupled elastic–plastic wave structure in diamond was directly observed and studied with submicrometer spatial resolution, using the unique capabilities of the x-ray free-electron laser (XFEL). The direct measurements allowed, for the first time, the fitting and validation of the 2D failure model for diamond in the range of several Mbar. Our experimental approach opens new possibilities for the direct verification and construction of equations of state of matter in the ultra-high-stress range, which are relevant to solving a variety of problems in high-energy-density physics.

36 MATERIALS SCIENCE↗