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 703 records · Page 39

Preconditioned methods for solving the incompressible and low speed compressible equations

Acceleration methods are presented for solving the steady state incompressible equations. These systems are preconditioned by introducing artificial time derivatives which allow for a faster convergence to the steady state. The compressible equations in conservation form with slow flow are also considered. Two arbitrary functions, alpha and beta, are introduced in the general preconditioning. An analysis of this system is presented and an optimal value for beta is determined given a constant, alpha. It is further shown that the resultant incompressible equations form a symmetric hyperbolic system and so are well posed. Several generalizations to the compressible equations are presented which generalize previous results.

Turkel, E.↗

The thermodynamic properties of oxygen and nitrogen. Part 1: Thermodynamic properties of nitrogen from 115 R to 3500 R with pressures to 150000 psia

An equation of state is presented for liquid and gaseous nitrogen for temperatures from 115 R to 3500 R and pressures to 150,000 psia. All of the pressure-density-temperature data available from the published literature have been reviewed, and appropriate corrections have been identified and applied to bring experimental temperatures into accord with the International Practical Temperature Scale of 1968. Comparisons of property values calculated from the equation of state to measured values are included to illustrate the accuracy of the equation in representing the data. The coefficients of the equation of state were determined by a weighted least squares fit to selected published data and, simultaneously, to constant volume data determined by corresponding states analysis from oxygen data, and to data which define the phase equilibrium criteria for the saturated liquid and saturated vapor. The methods of weighting the various data for simultaneous fitting are presented and discussed. The equation of state is estimated to be accurate to within 0.5 percent in the liquid region, to within 0.1 percent for supercritical isotherms up to 15,000 psia, and to within 0.3 percent from 15,000 to 150,000 psia.

Stewart, R. B.↗

Using tabulated NLTE data for Hohlraum simulations

Non-local thermodyamic equilibrium (NLTE) atomic kinetics is necessary in inertial confinement fusion Hohlraum simulations for adequately modeling high-Z walls and dopants but is computationally very expensive. We present here an approach for tabulating NLTE material data in an economical manner. Material properties for an arbitrary radiation field are provided through tabulated data for a limited set of radiation fields plus derivatives with respect to components of those radiation fields. We have implemented this method in a radiation-hydrodynamics code and compare Hohlraum simulations done with inline NLTE calculations with those done with tables constructed from the same atomic data. The results demonstrate that NLTE tables can replace inline calculations in Hohlraum simulations without a significant loss of accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The role of vectors in reheating

We explore various aspects concerning the role of vector bosons during the reheating process. Generally, reheating occurs during the period of oscillations of the inflaton condensate and the evolution of the radiation bath depends on the inflaton equation of state. For oscillations about a quadratic minimum, the equation of state parameter, w = p/ρ = 0, and the evolution of the temperature, T(a) with respect to the scale factor is independent of the spin of the inflaton decay products. However, for cases when w > 0, there is a dependence on the spin, and here we consider the evolution when the inflaton decays or scatters to vector bosons. We also investigate the gravitational production of vector bosons as potential dark matter candidates. Gravitational production predominantly occurs through the longitudinal mode. We compare these results to the gravitational production of scalars.

79 ASTRONOMY AND ASTROPHYSICS↗

Three-dimensional finite element analysis for high velocity impact

A finite element algorithm for solving unsteady, three-dimensional high velocity impact problems is presented. A computer program was developed based on the Eulerian hydroelasto-viscoplastic formulation and the utilization of the theorem of weak solutions. The equations solved consist of conservation of mass, momentum, and energy, equation of state, and appropriate constitutive equations. The solution technique is a time-dependent finite element analysis utilizing three-dimensional isoparametric elements, in conjunction with a generalized two-step time integration scheme. The developed code was demonstrated by solving one-dimensional as well as three-dimensional impact problems for both the inviscid hydrodynamic model and the hydroelasto-viscoplastic model.

Chan, S. T. K.↗

Structure of neutron star envelopes

The envelopes of nonmagnetic neutron stars are studied using the best available opacities and equation of state. The general relativistic equations of the structure and evolution of nonmagnetic neutron stars are discussed, and it is shown that they can be reduced to a single equation for calculating the thermal structure of neutron star envelopes. The physical input needed to solve the thermal structure equation is investigated and the numerical results of envelope model calculations are presented. It is shown that the thermal structure of neutron star envelopes is determined by the single parameter T(s) to the 4th/g(s), where T(s) is the effective surface temperature and g(s) the surface gravity of the star. This result is used to derive a number of other scaling relations, and the effects of general relativity on the envelope thermal structure are examined in detail. The results of a sensitivity analysis of the regional opacity needed to obtain a reliable relationship between the temperatures of the inner and outer boundaries of the envelope is presented.

Gudmundsson, E. H.↗

Entropy Splitting and Numerical Dissipation

A rigorous stability estimate for arbitrary order of accuracy of spatial central difference schemes for initial-boundary value problems of nonlinear symmetrizable systems of hyperbolic conservation laws was established recently by Olsson and Oliger (1994) and Olsson (1995) and was applied to the two-dimensional compressible Euler equations for a perfect gas by Gerritsen and Olsson (1996) and Gerritsen (1996). The basic building block in developing the stability estimate is a generalized energy approach based on a special splitting of the flux derivative via a convex entropy function and certain homogeneous properties. Due to some of the unique properties of the compressible Euler equations for a perfect gas, the splitting resulted in the sum of a conservative portion and a non-conservative portion of the flux derivative. hereafter referred to as the "Entropy Splitting." There are several potential desirable attributes and side benefits of the entropy splitting for the compressible Euler equations that were not fully explored in Gerritsen and Olsson. The paper has several objectives. The first is to investigate the choice of the arbitrary parameter that determines the amount of splitting and its dependence on the type of physics of current interest to computational fluid dynamics. The second is to investigate in what manner the splitting affects the nonlinear stability of the central schemes for long time integrations of unsteady flows such as in nonlinear aeroacoustics and turbulence dynamics. If numerical dissipation indeed is needed to stabilize the central scheme, can the splitting help minimize the numerical dissipation compared to its un-split cousin? Extensive numerical study on the vortex preservation capability of the splitting in conjunction with central schemes for long time integrations will be presented. The third is to study the effect of the non-conservative proportion of splitting in obtaining the correct shock location for high speed complex shock-turbulence interactions. The fourth is to determine if this method can be extended to other physical equations of state and other evolutionary equation sets. If numerical dissipation is needed, the Yee, Sandham, and Djomehri (1999) numerical dissipation is employed. The Yee et al. schemes fit in the Olsson and Oliger framework.

Yee, H. C.↗

Converging towards atomic and nuclear pressures (Final Report)

his report details the development of transformational tools and techniques, allowing the accurate characterization of matter at the most extreme conditions yet studied in the high energy density (HED) domain. Before this work, most experiments quantitatively mapping the nature of matter in the HED realm were performed using planar geometry, which allows for the isolation and measurement of variables in a single thermodynamic state. Such measurements include equation of state variables, optical conductivity, heat transport and more. However, due to a variety of plasma processes, such experiments are limited to below ~10 TPa (100 Mbar) pressures. To achieve higher pressures, convergent experiments are needed. Historically, convergent experiments have not been used for benchmark data because they are integral measurements. That is each part of a convergent target undergoes a time dependent wide range of states, making it difficult to accurately isolate any particular quantity for a given thermodynamic state. This effort developed innovative convergent HED platforms and techniques enabling the exploration of matter from atomic to nuclear scale pressures. This effort also performed pioneering experiments that yielded the first rigorous benchmark data at these extreme conditions. This funding award began with the overarching goal to understand the behavior of matter at extreme (atomic-to-nuclear scale) pressures. The motivation for these goals lie in the fact that every time scientists explore matter beyond the threshold of an atomic unit, there is a fundamental shift in science. The atomic unit for energy, mass, charge, length, and time have each been explored, and each time such a threshold was crossed, a new sequence of discoveries was made resulting in significant awards such as the Nobel Prize. The only unexplored atomic unit is pressure, and this effort set the course for exploring matter at and beyond such pressures. That most of the recently discovered extrasolar planets and stars, as well as matter in the late implosion stages of inertial fusion targets, have deep internal pressures at and beyond atomic pressures amplifies the importance of this effort. Much of the initial work in developing techniques to create these atomic-to-nuclear pressures already existed within the HED community, predominantly at large scale laser facilities such as the National Ignition Facility (NIF) and the Omega60 laser the University of Rochester and although these experiments were routinely performed the ability to extract information about the underlying physical states and processes remained elusive due to the complexity of the experiments, extreme scales in both time and space, and the integrated nature of all the measurements. This work built a rigorous framework so such measurements can routinely be made. This was done in part through the introduction of Bayesian inference techniques into the field of HED science. These techniques allow for the self-consistent extraction of the relevant variables and their explicit and implicit correlations, so as to make use of integrated experimental data to constrain physical models and provide accurate uncertainty bounds for the data. The techniques developed here are fully transparent and proved immediately useful providing quantitative rigorous benchmarks for physical states and processes at some of the most extreme conditions yet explored on Earth. This funding is directly or partly responsible for 7 publications in major peer-reviewed scientific journals, a doctoral thesis, 3 invited talks at major conferences, and the training of a post-doc, 2 graduate students, and 1 undergraduate student. Beyond the effort supported by this award launched the introduction of Bayesian inference into the HED physics community and helped push the usage if modern data-science techniques within the physics community at large.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theoretical priors in scalar-tensor cosmologies: Thawing quintessence

The late time acceleration of the Universe can be characterized in terms of an extra, time-dependent, component of the Universe - dark energy. The simplest proposal for dark energy is a scalar-tensor theory - quintessence - which consists of a scalar field, φ, whose dynamics is solely dictated by its potential, V(φ). Such a theory can be uniquely characterized by the equation of state of the scalar field energy momentum-tensor. We find the time dependence of the equation of state for a broad family of potentials and, using this information, we propose an analytic prior distribution for the most commonly used parametrization. Here, we show that this analytic prior can be used to accurately predict the distribution of observables for the next generation of cosmological surveys. Including the theoretical priors in the comparison with observations considerably improves the constraints on the equation of state.

79 ASTRONOMY AND ASTROPHYSICS↗

Neutron Star Radii, Deformabilities, and Moments of Inertia from Experimental and Ab Initio Theory Constraints of the 208Pb Neutron Skin Thickness

Recent experimental and ab initio theory investigations of the 208Pb neutron skin thickness have the potential to inform the neutron star equation of state. In particular, the strong correlation between the 208Pb neutron skin thickness and the pressure of neutron matter at normal nuclear densities leads to modified predictions for the radii, tidal deformabilities, and moments of inertia of typical 1.4M⊙ neutron stars. In the present work, we study the relative impact of these recent analyses of the 208Pb neutron skin thickness on bulk properties of neutron stars within a Bayesian statistical analysis. Two models for the equation of state prior are employed in order to highlight the role of the highly uncertain high-density equation of state. From our combined Bayesian analysis of nuclear theory, nuclear experiment, and observational constraints on the dense matter equation of state, we find at the 90% credibility level R1.4=12.36−0.73+0.38 km for the radius of a 1.4M⊙ neutron star, R2.0=11.96−0.71+0.94 km for the radius of a 2.0M⊙ neutron star, Λ1.4=440−144+103 for the tidal deformability of a 1.4M⊙ neutron star, and I1.338=1.425−0.146+0.074×1045gcm2 for the moment of inertia of PSR J0737-3039A whose mass is 1.338M⊙.

Lim, Yeunhwan↗

A comparison of the Method of Lines to finite difference techniques in solving time-dependent partial differential equations

Steady state solutions to two time dependent partial differential systems have been obtained by the Method of Lines (MOL) and compared to those obtained by efficient standard finite difference methods: (1) Burger's equation over a finite space domain by a forward time central space explicit method, and (2) the stream function - vorticity form of viscous incompressible fluid flow in a square cavity by an alternating direction implicit (ADI) method. The standard techniques were far more computationally efficient when applicable. In the second example, converged solutions at very high Reynolds numbers were obtained by MOL, whereas solution by ADI was either unattainable or impractical. With regard to 'set up' time, solution by MOL is an attractive alternative to techniques with complicated algorithms, as much of the programming difficulty is eliminated.

Kurtz, L. A.↗

The structure of the planets Jupiter and Saturn

Planetary models for Jupiter and Saturn are computed using a fourth-order theory and an equation of state for the molecular hydrogen and helium planetary envelopes taken from the Monte Carlo calculations of Slattery and Hubbard (1976). Models for Jupiter are found that have a small amount of heavy elements either mixed with hydrogen and helium throughout the interior of the planet or concentrated in a small dense core. Saturn is modeled with a solar-composition hydrogen and helium envelope and a small dense core. It is concluded that the molecular equation of state linked with suitable interior equations of state can produce Jovian models which satisfy the observational data. The planetary models show that the enrichment of heavy elements (relative to solar composition) is approximately 3 times for Jupiter and 10 times for Saturn.

Slattery, W. L.↗

Conservative-variable average states for equilibrium gas multi-dimensional fluxes

Modern split component evaluations of the flux vector Jacobians are thoroughly analyzed for equilibrium-gas average-state determinations. It is shown that all such derivations satisfy a fundamental eigenvalue consistency theorem. A conservative-variable average state is then developed for arbitrary equilibrium-gas equations of state and curvilinear-coordinate fluxes. Original expressions for eigenvalues, sound speed, Mach number, and eigenvectors are then determined for a general average Jacobian, and it is shown that the average eigenvalues, Mach number, and eigenvectors may not coincide with their classical pointwise counterparts. A general equilibrium-gas equation of state is then discussed for conservative-variable computational fluid dynamics (CFD) Euler formulations. The associated derivations lead to unique compatibility relations that constrain the pressure Jacobian derivatives. Thereafter, alternative forms for the pressure variation and average sound speed are developed in terms of two average pressure Jacobian derivatives. Significantly, no additional degree of freedom exists in the determination of these two average partial derivatives of pressure. Therefore, they are simultaneously computed exactly without any auxiliary relation, hence without any geometric solution projection or arbitrary scale factors. Several alternative formulations are then compared and key differences highlighted with emphasis on the determination of the pressure variation and average sound speed. The relevant underlying assumptions are identified, including some subtle approximations that are inherently employed in published average-state procedures. Finally, a representative test case is discussed for which an intrinsically exact average state is determined. This exact state is then compared with the predictions of recent methods, and their inherent approximations are appropriately quantified.

Iannelli, G. S.↗

Novel scalings of neutron star properties from analyzing dimensionless Tolman–Oppenheimer–Volkoff equations

The Tolman–Oppenheimer–Volkoff (TOV) equations govern the radial evolution of pressure and energy density in static neutron stars (NSs) in hydrodynamical equilibrium. Using the reduced pressure and energy density with respect to the NS central energy density, the original TOV equations can be recast into dimensionless forms. While the traditionally used integral approach for solving the original TOV equations require an input nuclear Equation of State (EOS), the dimensionless TOV equations can be anatomized by using the reduced pressure and energy density as polynomials of the reduced radial coordinate without using any input nuclear EOS. It has been shown in several of our recent works that interesting and novel perspectives about NS core EOS can be extracted directly from NS observables by using the latter approach. Our approach is based on intrinsic and perturbative analyses of the dimensionless (IPAD) TOV equations (IPAD-TOV). In this review article, we first discuss the length and energy density scales of NSs as well as the dimensionless TOV equations for scaled variables and their perturbative solutions near NS cores. We then review several new insights into NS physics gained from solving perturbatively the scaled TOV equations. Whenever appropriate, comparisons with the traditional approach from solving the original TOV equations will be made. In particular, we first show that the nonlinearity of the TOV equations basically excludes a linear EOS for dense matter in NS cores. We then show that perturbative analyses of the scaled TOV equations enable us to reveal novel scalings of the NS mass, radius and the compactness with certain combinations of the NS central pressure and energy density. Thus, observational data on either mass, radius or compactness can be used to constrain directly the core EOS of NS matter independent of the still very uncertain nuclear EOS models. As examples, the EOS of the densest visible matter in our Universe before the most massive neutron stars collapse into black holes (BHs) as well as the central EOS of a canonical or a 2.1 solar mass NS are extracted without using any nuclear EOS model. In addition, we show that causality in NSs sets an upper bound of about 0.374 for the ratio of pressure over energy density and correspondingly a lower limit for trace anomaly in supra-dense matter. We also demonstrate that the strong-field gravity plays a fundamental role in extruding a peak in the density/radius profile of the speed of sound squared (SSS) in massive NS cores independent of the nuclear EOS. Finally, some future perspectives of NS research using the new approach reviewed here by solving perturbatively the dimensionless TOV equations are outlined.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Theoretical models of X-ray emission from rich clusters of galaxies

After assessing the theory of cluster gas with polytropic equations of state, and discussing the prediction of local polytropic indices in the gas distribution, attention is given to the consequences of modifying the equation of state by thermal cooling in the central regions of the cluster. It is suggested that, if the equation of state can be understood locally, an attempt may be made at (1) predicting the X-ray luminosities of the clusters at distant cosmological epochs and (2) measuring the evolution of the cluster's potential, through observation of the evolution of X-ray luminosity and surface brightness.

Cowie, L. L.↗

Techniques for studying materials under extreme states of high energy density compression

The properties of materials under extreme conditions of pressure and density are of key interest to a number of fields, including planetary geophysics, materials science, and inertial confinement fusion. In geophysics, the equations of state of planetary materials, such as hydrogen and iron, under ultrahigh pressure and density provide a better understanding of their formation and interior structure [Celliers et al., “Insulator-metal transition in dense fluid deuterium,” Science 361, 677–682 (2018) and Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)]. The processes of interest in these fields occur under conditions of high pressure (100 GPa–100 TPa), high temperature (>3000 K), and sometimes at high strain rates (>103 s−1) depending on the process. With the advent of high energy density (HED) facilities, such as the National Ignition Facility (NIF), Linear Coherent Light Source, Omega Laser Facility, and Z, these conditions are reachable and numerous experimental platforms have been developed. To measure compression under ultrahigh pressure, stepped targets are ramp-compressed and the sound velocity, measured by the velocity interferometer system for any reflector diagnostic technique, from which the stress-density of relevant materials is deduced at pulsed power [M. D. Knudson and M. P. Desjarlais, “High-precision shock wave measurements of deuterium: Evaluation of exchange-correlation functionals at the molecular-to-atomic transition,” Phys. Rev. Lett. 118, 035501 (2017)] and laser [Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)] facilities. To measure strength under high pressure and strain rates, experimenters measure the growth of Rayleigh–Taylor instabilities using face-on radiography [Park et al., “Grain-size-independent plastic flow at ultrahigh pressures and strain rates,” Phys. Rev. Lett. 114, 065502 (2015)]. The crystal structure of materials under high compression is measured by dynamic x-ray diffraction [Rygg et al., “X-ray diffraction at the national ignition facility,” Rev. Sci. Instrum. 91, 043902 (2020) and McBride et al., “Phase transition lowering in dynamically compressed silicon,” Nat. Phys. 15, 89–94 (2019)]. Medium range material temperatures (a few thousand degrees) can be measured by extended x-ray absorption fine structure techniques, Yaakobi et al., “Extended x-ray absorption fine structure measurements of laser-shocked V and Ti and crystal phase transformation in Ti,” Phys. Rev. Lett. 92, 095504 (2004) and Ping et al., “Solid iron compressed up to 560 GPa,” Phys. Rev. Lett. 111, 065501 (2013), whereas more extreme temperatures are measured using x-ray Thomson scattering or pyrometry. This manuscript will review the scientific motivations, experimental techniques, and the regimes that can be probed for the study of materials under extreme HED conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Newton's method applied to finite-difference approximations for the steady-state compressible Navier-Stokes equations

Finite-difference approximations for steady-state compressible Navier-Stokes equations, whose two spatial dimensions are written in generalized curvilinear coordinates and strong conservation-law form, are presently solved by means of Newton's method in order to obtain a lifting-airfoil flow field under subsonic and transonnic conditions. In addition to ascertaining the computational requirements of an initial guess ensuring convergence and the degree of computational efficiency obtainable via the approximate Newton method's freezing of the Jacobian matrices, attention is given to the need for auxiliary methods assessing the temporal stability of steady-state solutions. It is demonstrated that nonunique solutions of the finite-difference equations are obtainable by Newton's method in conjunction with a continuation method.

Bailey, Harry E.↗

Impact of moment-based, energy integrated neutrino transport on microphysics and ejecta in binary neutron star mergers

We present an extensive study of the effects of neutrino transport in three-dimensional general relativistic radiation hydrodynamics (GRHD) simulations of binary neutron star (BNS) mergers using our moment-based, energy-integrated neutrino radiation transport (M1) scheme. Here, we consider a total of eight BNS configurations, while varying equation of state models, mass ratios, and grid resolutions, for a total of 16 simulations. We find that M1 neutrino transport is crucial in modeling the local absorption of neutrinos and the deposition of lepton number throughout the medium. We provide an in-depth look at the effects of neutrinos on the fluid dynamics and luminosity during the late inspiral and postmerger phases, the properties of ejecta and outflow, and the postmerger nucleosynthesis. The simulations presented in this work comprise an extensive study of the combined effect of the equation of state and M1 neutrino transport in GRHD simulations of BNS mergers, and establish that the solution provided by our M1 scheme is robust across system properties.

150 ≤ A ≤ 18959 ≤ A ≤ 8990 ≤ A ≤ 149↗