Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complete 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing

A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.

tubulent multicompenent-species↗
NASA NTRS1972

Effect of support flexibilty and damping on the dynamic response of a single mass flexible rotor in elastic bearings

A steady state analysis of the shaft and the bearing housing motion was made by assuming synchronous precession of the system. The conditions under which the support system would act as a dynamic vibration absorber at the rotor critical speed were studied; plots of the rotor and support amplitudes, phase angles, and forces transmitted were evaluated by the computer, and the performance curves were automatically plotted by a CalComp plotter unit. Curves are presented on the optimization of the support housing characteristics to attenuate the rotor unbalance response over the entire rotor speed range. The complete transient motion including rotor unbalance was examined by integrating the equations of motion numerically using a modified fourth order Runge-Kutta procedure, and the resulting whirl orbits were plotted by the CalComp plotter unit. The results of the transient analysis are discussed with regards to the design optimization procedure derived from the steady-state analysis.

Kirk, R. G.↗
DOE OSTI2021

Characterizing Topological Order with Matrix Product Operators

One of the most striking features of gapped quantum phases that exhibit topological order is the presence of long-range entanglement that cannot be detected by any local order parameter. The formalism of projected entangled-pair states is a natural framework for the parameterization of gapped ground state wavefunctions which allows one to characterize topological order in terms of the virtual symmetries of the local tensors that encode the wavefunction. In their most general form, these symmetries are represented by matrix product operators acting on the virtual level, which leads to a set of algebraic rules characterizing states with topological quantum order. This construction generalizes the concepts of $\mathsf G$- and twisted injectivity; the corresponding matrix product operators encode all topological features of the theory and provide a complete picture of the ground state manifold on the torus. We show how the string-net models of Levin and Wen fit within this formalism and in doing so provide a particularly intuitive interpretation of the pentagon equation for F-symbols as the pulling of matrix product operators through the string-net tensor network. Our approach paves the way to finding novel topological phases beyond string nets and elucidates the description of topological phases in terms of entanglement Hamiltonians and edge theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Accurate effective fluid approximation for ultralight axions

Ultralight axions are theoretically interesting and phenomenologically rich dark sector candidates, but they are difficult to track across cosmological timescales because of their fast oscillations. We resolve this problem by developing a novel method to evolve them efficiently and accurately. We first construct an exact effective fluid which at late times matches the axion but which evolves in a simple way. We then approximate this evolution with a carefully chosen equation of state and sound speed. With our scheme we find that we can obtain subpercent accuracy for the linear theory suppression of axion density fluctuations relative to that of cold dark matter without tracking even a single complete oscillation of the axion field. Furthermore, we use our technique to test other approximation schemes and to provide a fitting formula for the transfer function for the matter power spectrum in linear theory in axion models. Implementing our approach in existing cosmological axion codes is straightforward and will help unleash the potential of high-precision next-generation experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Block smoothers and generalized ideal interpolation in AMG (Final Report)

The Pennsylvania State University (“Subcontractor”) worked on developing new parallel algebraic multilevel methods suitable for solving PDEs. Specifically, work on the design of multigrid solvers for coupled systems of partial differential equations arising in numerical modeling of various applications was completed. A main emphasis was on the design of new ideal algebraic multigrid interpolation for problems such as Maxwell’s equations where block smoothers are needed and the standard form of ideal interpolation is not an effective choice.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1975

Numerical computation of viscous flows on the lee side of blunt shapes flying at supersonic speeds

A numerical method for solving the parabolic approximation to the steady-state compressible Navier-Stokes equations is examined. The approximation neglects only the streamwise gradients of shear stress. An implicit finite difference method is used which advances the solution downstream from an initial data surface and determines the complete viscous-inviscid flow between the body and bow shock wave. It is necessary that the inviscid portion of the flow field be supersonic. Crossflow separation is determined as part of the solution. The method is applied to a 15 deg sphere-cone at 15 deg angle of attack, and the results are compared with an inviscid method-of-characteristics calculation.

Rakich, J. V.↗
DOE OSTI2021

Near-optimal combination of high performance and insensitivity in a shock compressed high explosive single crystal

Achieving the desired combination of superior detonation performance and insensitivity to shock initiation has been a long-standing goal in high explosive (HE) science and technology. Having previously established the shock insensitivity of 1,1-diamino-2,2-dinitroethene (also known as DADNE or FOX-7) single crystals to 20 GPa (extended to 25 GPa in this work), the FOX-7 detonation response was determined through wave profile measurements in ~250 μm thick single crystals shock compressed to 64 GPa. Quite unexpectedly, FOX-7 demonstrated the classic Chapman–Jouguet (C–J) detonation response—reaction completion in the detonation front (<0.7 ns) at pressures of 44 GPa and higher—not observed in other insensitive high explosives. Here, the experimentally determined C–J pressure (35 GPa), detonation wave velocities and the detonation products equation of state—together with shock insensitivity to 25 GPa—demonstrate that FOX-7 single crystals display a near-optimal combination of high performance and shock insensitivity, not observed in another HE crystal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

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↗
DOE OSTI2022

PUSHing Core-collapse Supernovae to Explosions in Spherical Symmetry. V. Equation of State Dependency of Explosion Properties, Nucleosynthesis Yields, and Compact Remnants

Abstract In this fifth paper of the series, we use the parameterized, spherically symmetric explosion method PUSH to investigate the impact of eight different nuclear equations of state (EOS). We present and discuss the explosion properties and the detailed nucleosynthesis yields, and predict the remnant (neutron star or black hole) for all our simulations. For this, we perform two sets of simulations. First, a complete study of nonrotating stars from 11 to 40 M ⊙ at three different metallicities using the SFHo EOS; and, second, a suite of simulations for four progenitors (16 M ⊙ at three metallicities and 25 M ⊙ at solar metallicity) for eight different nuclear EOS. We compare our predicted explosion energies and yields to observed supernovae and to the metal-poor star HD 84937. We find EOS-dependent differences in the explosion properties and the nucleosynthesis yields. However, when comparing to observations, these differences are not large enough to rule out any EOS considered in this work.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2004

On the Inversion for Mass (Re)Distribution from Global (Time-Variable) Gravity Field

The well-known non-uniqueness of the gravitational inverse problem states the following: The external gravity field, even if completely and exactly known, cannot Uniquely determine the density distribution of the body that produces the gravity field. This is an intrinsic property of a field that obeys the Laplace equation, as already treated in mathematical as well as geophysical literature. In this paper we provide conceptual insight by examining the problem in terms of spherical harmonic expansion of the global gravity field. By comparing the multipoles and the moments of the density function, we show that in 3-S the degree of knowledge deficiency in trying to inversely recover the density distribution from external gravity field is (n+l)(n+2)/2 - (2n+l) = n(n-1)/2 for each harmonic degree n. On the other hand, on a 2-D spherical shell we show via a simple relationship that the inverse solution of the surface density distribution is unique. The latter applies quite readily in the inversion of time-variable gravity signals (such as those observed by the GRACE space mission) where the sources over a wide range of the scales largely come from the Earth's Surface.

Chao, Benjamin F.↗
NASA NTRS2000

Resin Film Infusion (RFI) Process Modeling for Large Transport Aircraft Wing Structures

This investigation completed the verification of a three-dimensional resin transfer molding/resin film infusion (RTM/RFI) process simulation model. The model incorporates resin flow through an anisotropic carbon fiber preform, cure kinetics of the resin, and heat transfer within the preform/tool assembly. The computer model can predict the flow front location, resin pressure distribution, and thermal profiles in the modeled part. The formulation for the flow model is given using the finite element/control volume (FE/CV) technique based on Darcy's Law of creeping flow through a porous media. The FE/CV technique is a numerically efficient method for finding the flow front location and the fluid pressure. The heat transfer model is based on the three-dimensional, transient heat conduction equation, including heat generation. Boundary conditions include specified temperature and convection. The code was designed with a modular approach so the flow and/or the thermal module may be turned on or off as desired. Both models are solved sequentially in a quasi-steady state fashion. A mesh refinement study was completed on a one-element thick model to determine the recommended size of elements that would result in a converged model for a typical RFI analysis. Guidelines are established for checking the convergence of a model, and the recommended element sizes are listed. Several experiments were conducted and computer simulations of the experiments were run to verify the simulation model. Isothermal, non-reacting flow in a T-stiffened section was simulated to verify the flow module. Predicted infiltration times were within 12-20% of measured times. The predicted pressures were approximately 50% of the measured pressures. A study was performed to attempt to explain the difference in pressures. Non-isothermal experiments with a reactive resin were modeled to verify the thermal module and the resin model. Two panels were manufactured using the RFI process. One was a stepped panel and the other was a panel with two 'T' stiffeners. The difference between the predicted infiltration times and the experimental times was 4% to 23%.

Loos, Alfred C.↗
DOE OSTI2022

The development of natural gas and hydrogen pipeline capital cost estimating equations

Natural gas pipeline cost data collected by the Oil and Gas Journal (O&GJ) for interstate pipelines constructed from 1980 through 2017 were used to develop capital cost estimating equations that are a function of pipeline diameter, length, and U.S. region. Equations were developed for material, labor, miscellaneous, and right-of-way costs, the four cost components in the O&GJ data, for six different regions of the United States (U.S.). Each equation is a function of pipeline diameter and length. Adjustment mechanisms were then developed for converting the natural gas pipeline equations into equations for estimating the costs of hydrogen pipelines. These adjustments were based in part on an analysis completed by the National Institute for Standards and Technology (NIST). Here the results of this work were used to update cost models in the Hydrogen Delivery Scenario Analysis Model (HDSAM), developed by Argonne National Laboratory for the U.S. Department of Energy’s Hydrogen Program. Our analysis shows a wide range of pipeline cost across different U.S. regions, especially with respect to labor and right-of-way costs. The developed cost formulas for hydrogen pipelines are both important and timely as hydrogen is being considered as a zero-carbon energy carrier with the potential to decarbonize all energy sectors, and the cost of hydrogen transportation is essential for techno-economic analysis of its potential use in these sectors.

03 NATURAL GAS↗
NASA NTRS1984

The Stability of the Low Degree 5 Minute Solar Oscillations

The decay rate for many of the low degree p modes observed as 5 minute oscillations of the Sun is discussed and decay rates for modes ranging from radial (l=0) to the nonradial ones with l=5 for overtones 10 through 28. Parameters needed for he solar model are given and the hydrogen mass fraction composition structure is shown as well as the structure given by Christensen-Dalsgaard (1982) for an evolved solar model. The special equation of state and opacity table with X=0.74 for the hydrogen mass fraction in the outer 0.40 of the mass needs slightly more hydrogen in the central regions than obtained by Christensen-Dalsgaard in order to give a complete and consistent model. The difference in helium production between these two models is about 10%, meaning that the total energy radiated by the Sun during its lifetime thus far agrees satisfactorily with accurately calculated evolution sequences.

Kidman, R. B.↗
DOE OSTI2024

Toward UV models of kinetic mixing and portal matter. VI. A more complex dark matter sector?

Portal matter (PM), having both Standard Model (SM) and dark sector charges, can induce kinetic mixing between the 𝑈⁢(1) 𝐷 dark photon and the SM gauge fields at the 1-loop level offering an attractive mechanism by which light (≲1 GeV) thermal dark matter (DM) can interact with visible matter and obtain its observed relic density. In doing so, if the DM is fermionic, the CMB and other astrophysical observations inform us that it must be Majorana/pseudo-Dirac in nature to avoid velocity/temperature-independent 𝑠-wave annihilation to SM final states. How does this idea fit into a more UV-complete picture also including the SM interactions? There are some reasons to believe that at least a first step along this path may not lie too far away in energy due to the renormalization group equations running of the dark gauge coupling, which for a significant range of parameters, becomes nonperturbative at/before the ∼10’s of TeV energy range. This implies that 𝑈⁢(1) 𝐷 must become embedded in an asymptotically free, non-Abelian group, 𝐺 𝐷 , before this can occur. The breaking of this larger group then produces the masses for the PM and the additional gauge fields associated with 𝐺 𝐷 then can lead to new interactions between the SM and the dark sector. Following several bottom-up approaches, we have examined a set of distinctive and testable phenomenological features associated with this general setup, based upon a number of simplifying assumptions. Clearly, it behooves us to explore the impact of these specific assumptions on these predictions for the array of possible experimental tests of this class of models. In most past analyses it has been assumed that DM is a vectorlike, complex singlet under the group 𝐺 𝐷 . If this assumption is relaxed, the dark sector must be augmented by additional fermion(s) and the associated scalar fields needed to break the gauge symmetries while generating the needed Majorana-like mass terms for the DM. In this paper, we analyze the simplest extension of this kind wherein the DM lies in a vectorlike doublet of 𝐺 𝐷 , which we take to have the structure 𝑆⁢𝑈⁢(2) 𝐼 ×𝑈⁢(1) 𝑌 𝐼 as in earlier work, leading to new phenomenological implications. We find, for example, that given the current LHC search constraints on the masses of heavy gauge bosons, the production of these new dark states with large rates is unlikely to occur at colliders unless they are produced singly in 𝑔⁡𝑔 fusion or their pair production cross sections are resonantly enhanced. Here, we also find that an additional mechanism arises to generate hierarchal neutrino masses in such a setup.

Extensions of Higgs sector↗
NASA NTRS1987

Finite Element Flux-Corrected Transport (FEM-FCT) for the Euler and Navier-Stokes equations

A high resolution finite element method for the solution of problems involving high speed compressible flows is presented. The method uses the concepts of flux-corrected transport and is presented in a form which is suitable for implementation on completely unstructured triangular or tetrahedral meshes. Transient and steady state examples are solved to illustrate the performance of the algorithm.

Loehner, Rainald↗
DOE OSTI2022

Solving quasi-free and quadratic Lindblad master equations for open fermionic and bosonic systems

The dynamics of Markovian open quantum systems are described by Lindblad master equations. For fermionic and bosonic systems that are quasi-free, i.e. with Hamiltonians that are quadratic in the ladder operators and Lindblad operators that are linear in the ladder operators, we derive the equation of motion for the covariance matrix. This determines the evolution of Gaussian initial states and the steady states, which are also Gaussian. Using ladder super-operators (a.k.a. third quantization), here we show how the Liouvillian can be transformed to a many-body Jordan normal form which also reveals the full many-body spectrum. Extending previous work by Prosen and Seligman, we treat fermionic and bosonic systems on equal footing with Majorana operators, shorten and complete some derivations, also address the odd-parity sector for fermions, give a criterion for the existence of bosonic steady states, cover non-diagonalizable Liouvillians also for bosons, and include quadratic systems. In extension of the quasi-free open systems, quadratic open systems comprise additional Hermitian Lindblad operators that are quadratic in the ladder operators. While Gaussian states may then evolve into non-Gaussian states, the Liouvillian can still be transformed to a useful block-triangular form, and the equations of motion for k-point Green’s functions form a closed hierarchy. Based on this formalism, results on criticality and dissipative phase transitions in such models are discussed in a companion paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Thermoelastic properties of tungsten at simultaneous high pressure and temperature

In this work, the compressional ( P ) and shear wave velocities ( S ) and unit cell volumes (densities) of polycrystalline tungsten (W) have been measured simultaneously up to 10.5 GPa and 1073 K using ultrasonic interferometry in conjunction with x-ray diffraction and x-radiography techniques. Thermoelastic properties of W were derived using different methods. We obtained the isothermal bulk modulus K T 0 = 310.3(1.5) GPa, its pressure derivative K ' T 0 = 4.4(3), its temperature derivative at constant pressure ( ∂ K T / ∂ T ) P = – 0.0138 ( 1 ) GPa K – 1 and at constant volume ( ∂ K T / ∂ T ) V = – 0.0050 GPa K – 1 , the thermal expansion α ( 0, T ) = 1.02(27) × 10 –5 + 7.39(3.2) × 10 –9 T (K –1 ), as well as the pressure derivative of thermal expansion ( ∂ α / ∂ P ) T = – 1.44 ( 1 ) × 10 – 7 K – 1 GP a – 1 based on the high-temperature Birch–Murnaghan equation of state (EOS), the Vinet EOS, and thermal pressure approach. Finite strain analysis allowed us to derive the elastic properties and their pressure/temperature derivatives independent of the choice of pressure scale. A least-squares fitting yielded K S 0 = 314.5(2.5) GPa, K S 0 ' = 4.45(9), (∂ K S /∂ T ) P = – 0.0076(6) GPa K –1 , G 0 = 162.4(9) GPa, G 0 ' = 1.8(1), (∂ G /∂ T ) P = – 0.0175(9) GPa K –1 , and α 298 K = 1.23 × 10 – 5 K – 1 . Fitting current data to the Mie–Grüneisen–Debye EOS with derived θ 0 = 383.4 K yielded γ 0 = 1.81 ( 6 ) and q = 0.3 . The thermoelastic parameters obtained from various approaches are consistent with one another and comparable with previous results within uncertainties. Our current study provides a complete and self-consistent dataset for the thermoelastic properties of tungsten at high P – T conditions, which is important to improve the theoretical modeling of these materials under dynamic conditions.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Finite element flux-corrected transport (FEM-FCT) for the Euler and Navier-Stokes equations

A high resolution finite element method for the solution of problems involving high speed compressible flows is presented. The method uses the concepts of flux-corrected transport and is presented in a form which is suitable for implementation on completely unstructured triangular or tetrahedral meshes. Transient and steady-state examples are solved to illustrate the performance of the algorithm.

Lohner, Rainald↗