Search NASA⌕ Search

SEARCH · Search NASA

Results for “entropy viscosity”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

GASP: A computer code for calculating the thermodynamic and transport properties for ten fluids: Parahydrogen, helium, neon, methane, nitrogen, carbon monoxide, oxygen, fluorine, argon, and carbon dioxide

A FORTRAN IV subprogram called GASP is discussed which calculates the thermodynamic and transport properties for 10 pure fluids: parahydrogen, helium, neon, methane, nitrogen, carbon monoxide, oxygen, fluorine, argon, and carbon dioxide. The pressure range is generally from 0.1 to 400 atmospheres (to 100 atm for helium and to 1000 atm for hydrogen). The temperature ranges are from the triple point to 300 K for neon; to 500 K for carbon monoxide, oxygen, and fluorine; to 600 K for methane and nitrogen; to 1000 K for argon and carbon dioxide; to 2000 K for hydrogen; and from 6 to 500 K for helium. GASP accepts any two of pressure, temperature and density as input conditions along with pressure, and either entropy or enthalpy. The properties available in any combination as output include temperature, density, pressure, entropy, enthalpy, specific heats, sonic velocity, viscosity, thermal conductivity, and surface tension. The subprogram design is modular so that the user can choose only those subroutines necessary to the calculations.

Hendricks, R. C.↗
NASA NTRS1989

Shock capturing by the spectral viscosity method

A main disadvantage of using spectral methods for nonlinear conservation laws lies in the formation of Gibbs phenomenon, once spontaneous shock discontinuities appear in the solution. The global nature of spectral methods than pollutes the unstable Gibbs oscillations overall the computational domain, and the lack of entropy dissipation prevents convergences in these cases. The Spectral Viscosity method, which is based on high frequency dependent vanishing viscosity regularization of the classical spectral methods is discussed. It is shown that this method enforces the convergence of nonlinear spectral approximations without sacrificing their overall spectral accuracy.

Tadmor, Eitan↗
DOE OSTI2024

Thermodynamically consistent Cahn–Hilliard–Navier–Stokes equations using the metriplectic dynamics formalism

Cahn–Hilliard–Navier–Stokes (CHNS) systems describe flows with two-phases, e.g., a liquid with bubbles. Obtaining constitutive relations for general dissipative processes for such systems, which are thermodynamically consistent, can be a challenge. We show how the metriplectic 4-bracket formalism (Morrison and Updike, 2024) achieves this in a straightforward, in fact algorithmic, manner. First, from the noncanonical Hamiltonian formulation for the ideal part of a CHNS system we obtain an appropriate Casimir to serve as the entropy in the metriplectic formalism that describes the dissipation (e.g. viscosity, heat conductivity and diffusion effects). General thermodynamics with the concentration variable and its thermodynamics conjugate, the chemical potential, are included. Having expressions for the Hamiltonian (energy), entropy, and Poisson bracket, we describe a procedure for obtaining a metriplectic 4-bracket that describes thermodynamically consistent dissipative effects. The 4-bracket formalism leads naturally to a general CHNS system that allows for anisotropic surface energy effects. Furthermore, this general CHNS system reduces to cases in the literature, to which we can compare.

Cahn–Hilliard↗
NASA NTRS1988

Improving the accuracy of central difference schemes

General difference approximations to the fluid dynamic equations require an artificial viscosity in order to converge to a steady state. This artificial viscosity serves two purposes. One is to suppress high frequency noise which is not damped by the central differences. The second purpose is to introduce an entropy-like condition so that shocks can be captured. These viscosities need a coefficient to measure the amount of viscosity to be added. In the standard scheme, a scalar coefficient is used based on the spectral radius of the Jacobian of the convective flux. However, this can add too much viscosity to the slower waves. Hence, it is suggested that a matrix viscosity be used. This gives an appropriate viscosity for each wave component. With this matrix valued coefficient, the central difference scheme becomes closer to upwind biased methods.

Turkel, Eli↗
NASA NTRS1983

Anomalous bulk viscosity of two-phase fluids and implications for planetary interiors

The irreversible entropy production is calculated for the imposition of a pressure perturbation on a two-phase medium composed of a dilute suspension of droplets (or snowflakes) and a liquid phase of other materials. An absence of metastability is assumed, allowing the relaxation to be dominated by the solute finite diffusivity. The fluid medium was found to display a behavior suggestive of a bulk viscosity near 10 trillion P, a finding that is significant for studies of dissipation in planetary cores for tidal or seismic disturbances. A minimum quality factor for acoustic or tidal pressure oscillations and the accompanying frequency are calculated. An example is provided in terms of helium rain clouds in the deep interiors of giant planets. Additionally, a tidal quality factor of 10 to the 15th is found necessary to account for Io volcanism and resurfacing on Enceladus.

Stevenson, D. J.↗
NASA NTRS1987

The convergence of spectral methods for nonlinear conservation laws

The convergence of the Fourier method for scalar nonlinear conservation laws which exhibit spontaneous shock discontinuities is discussed. Numerical tests indicate that the convergence may (and in fact in some cases must) fail, with or without post-processing of the numerical solution. Instead, a new kind of spectrally accurate vanishing viscosity is introduced to augment the Fourier approximation of such nonlinear conservation laws. Using compensated compactness arguments, it is shown that this spectral viscosity prevents oscillations, and convergence to the unique entropy solution follows.

Tadmor, Eitan↗
DOE OSTI2020

Locally adaptive artificial viscosity strategies for Lagrangian hydrodynamics

To accurately model inviscid flow with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity is introduced to convert kinetic energy into internal energy, thereby providing a mechanism to generate the required entropy increase across shocks. In this paper, we propose a new method for constructing an adaptive, artificial viscosity in the context of one-dimensional, staggered-grid Lagrangian hydrodynamics. Our adaptive, artificial viscosity is defined in terms of two parameters that depend on density locally in a neighborhood around a shock, and hence, vary cell-by-cell. Our methodology is based on building a reference set of pre-computed optimal, globally constant artificial viscosity coefficients for a family of isolated shock test problems. For arbitrary flows, the evaluation of the unknown coefficients is automated by first estimating shock intensity locally, and second computing the corresponding adaptive parameter value via interpolation over the data from the pre-computed reference set of optimal, globally constant coefficient values. To illustrate the performance of our new approach, we compare our results against two existing methods. The first method is a limiter-based approach, which relies on estimating velocity gradients of the flow, and the second method is utilized in several commercial codes. Finally, we demonstrate that our new adaptive methodology produces more accurate results for a variety of tests with propagating shock waves, as well as for the aforementioned family of isolated shock problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Effective viscosities in a hydrodynamically expanding boost-invariant QCD plasma

Background: The near-equilibrium properties of a QCD plasma can be encoded into transport coefficients such as bulk and shear viscosity. In QCD, the ratio of these transport coefficients to entropy density, ζ/$\textit{s}$ and η/$\textit{s}$, depends nontrivially on the plasma's temperature. This is unlike in conformal systems where they take constant values such as η/$\textit{s}$ = 1/(4π). Purpose: In this work, we show that in a (0 + 1)D boost-invariant fluid with no transverse expansion, a temperature-dependent ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ can be described by an equivalent effective viscosity $\langleζ/s\rangle_{\text{eff}}$ or $\langleη/s\rangle_{\text{eff}}$. This effective viscosity combines the actual temperature-dependent ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ with the temperature profile of the fluid. Further, we extend the concept of effective viscosity in systems with transverse expansion and discuss how effective viscosities can be used to identify families of ζ/$\textit{s(T)}$ and η /$\textit{s(T)}$ that lead to similar hydrodynamic evolution. Methods: The Navier-Stokes relativistic hydrodynamic equations are used to provide a first definition of effective viscosity, in (0 + 1)D and (1 + 1)D. In the (0 + 1)D case, the analysis is extended to Israel-Stewart-type second-order hydrodynamics to clarify the effect of higher-order hydrodynamics corrections on the effective viscosity. Results: In a boost-invariant fluid with no transverse expansion [(0 + 1)D], the effective viscosity is expressed as a simple integral of ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ over temperature, with a weight determined by the speed of sound of the fluid. The result is general for any equation of state with a moderate temperature dependence of the speed of sound, including the QCD equation of state. This definition of effective viscosity can be used to identify infinite families of ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ that produce essentially indistinguishable temperature profiles. In a boost-invariant cylindrical system [(1 + 1)D], a similar definition of effective viscosity is obtained in terms of characteristic trajectories in time and transverse direction. This leads to an infinite number of constraints on an infinite functional space for ζ/$\textit{s(T)}$ and η/$\textit{s(T)}$. Realistic examples are presented by using a finite number of constraints on a finite functional space. Conclusions: The definition of effective viscosity in a (0 + 1)D system clarifies how infinite families of ζ/$\textit{s(T)}$ and η/$\textit{s(T)}$ can result in nearly identical hydrodynamic temperature profiles. By extending the study to a boost-invariant cylindrical [(1 + 1)D] fluid, we identify an approximate but more general definition of effective viscosity that highlights the potential and limits of the concept of effective viscosity in fluids with limited symmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Empirical correlations between the function of entropy ( Z S ) and net artificial viscous work in a shock physics hydrocode

Entropy is a state variable that may be obtained from any thermodynamically complete equation of state (EOS). However, hydrocode calculations that output the entropy often contain numerical errors; this is not because of the EOS, but rather the solution techniques that are used in hydrocodes (especially Eulerian) such as convection, remapping, and artificial viscosity. Here, in this work, empirical correlations are investigated to reduce the errors in entropy without altering the solution techniques for the conservation of mass, momentum, and energy. Specifically, these correlations are developed for the function of entropy Z S , and they depend upon the net artificial viscous work, as determined via Sandia National Laboratories’ shock physics hydrocode CTH. These results are a continuation of a prior effort to implement the entropy-based CREST reactive burn model in CTH, and they are presented here to stimulate further interest from the shock physics community. Future work is planned to study higher-dimensional shock waves, shock wave interactions, and possible ties between the empirical correlations and a physical law.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1989

Thermophysical properties of Helium-4 from 0.8 to 1500 K with pressures to 2000 MPa

Tabular summary data of the thermophysical properties of fluid helium are given for temperatures from 0.8 to 1500 K, with pressures to 2000 MPa between 75 and 300 K, or to 100 MPa outside of this temperature band. Properties include density, specific heats, enthalpy, entropy, internal energy, sound velocity, expansivity, compressibility, thermal conductivity, and viscosity. The data are calculated from a computer program which is available from the National Institute of Standards and Technology. The computer program is based on carefully fitted state equations for both normal and superfluid helium.

Arp, Vincent D.↗
NASA NTRS2006

The Curious Events Leading to the Theory of Shock Waves

We review the history of the development of the modern theory of shock waves. Several attempts at an early-theory quickly collapsed for lack of foundations in mathematics and thermodynamics. It is not until the works of Rankine and later Hugoniot that a full theory is established. Rankine is the first to show that within the shock a non-adiabatic process must occur. Hugoniot showed that in the absence of viscosity and heat conduction conservation of energy implies conservation of entropy in smooth regions and a jump in entropy across a shock. Even after the theory is fully developed, old notions continue to pervade the literature well into the early part of the 20th Century.

Salas, Manuel D.↗
NASA NTRS1971

Extrapolation of thermophysical properties data for oxygen to high pressures (5000 to 10,000 psia) at low temperatures (100-600 R)

Thermophysical properties data for oxygen at pressures below 5000 psia have been extrapolated to higher pressures (5,000-10,000 psia) in the temperature range 100-600 R. The tables include density, entropy, enthalpy, internal energy, speed of sound, specific heat, thermal conductivity, viscosity, thermal diffusivity, Prandtl number, and dielectric constant.

Weber, L. A.↗
NASA NTRS1979

Isentropic decompression of fluids from crustal and mantle pressures

Criteria are derived according to which the flow of single-phase magmatic fluids and the rarefaction expansion of low-viscosity liquids and gases may be considered approximately isentropic. Graphs of entropy vs. density with contours of constant pressure and mass fraction are used to examine the possible thermodynamic histories of H2O and CO2 decompressing isentropically from crustal and upper mantle pressures; these graphs offer a simple visual representation of a number of thermodynamic variables involved in isentropic processes. It is shown how the graphs can be used to examine the behavior of volatiles that (1) ascend in volcanic systems originating at different depths within the earth, and (2) decompress from a shock Hugoniot state. Entropy-density graphs are presented separately for H2O and CO2.

Kieffer, S. W.↗
NASA NTRS1994

Energy and maximum norm estimates for nonlinear conservation laws

We have devised a technique that makes it possible to obtain energy estimates for initial-boundary value problems for nonlinear conservation laws. The two major tools to achieve the energy estimates are a certain splitting of the flux vector derivative f(u)(sub x), and a structural hypothesis, referred to as a cone condition, on the flux vector f(u). These hypotheses are fulfilled for many equations that occur in practice, such as the Euler equations of gas dynamics. It should be noted that the energy estimates are obtained without any assumptions on the gradient of the solution u. The results extend to weak solutions that are obtained as point wise limits of vanishing viscosity solutions. As a byproduct we obtain explicit expressions for the entropy function and the entropy flux of symmetrizable systems of conservation laws. Under certain circumstances the proposed technique can be applied repeatedly so as to yield estimates in the maximum norm.

Olsson, Pelle↗
NASA NTRS1959

Approximations for the Thermodynamic and Transport Properties of High-Temperature Air

The thermodynamic and transport properties of high-temperature air are found in closed form starting from approximate partition functions for the major components in air and neglecting all minor components. The compressibility, enthalpy, entropy, the specific heats, the speed of sound, the coefficients of viscosity and of thermal conductivity, and the Prandtl numbers for air are tabulated from 500 degrees to 15,000 degrees Kelvin over a range of pressure from 0.0001 to 100 atmospheres. The energy of air and the mol. fractions of the major components of air can be found from the tabulated values for compressibility and enthalpy. It is predicted that the Prandtl number for fully ionized air, which is in complete equilibrium, will become small compared to unity, the order of 0.01, and this implies that boundary layers in such flow will be very transparent to heat flux.

Hansen, C. Frederick↗
NASA NTRS1991

Entropy-Based Approach To Nonlinear Stability

NASA technical memorandum suggests schemes for numerical solution of differential equations of flow made more accurate and robust by invoking second law of thermodynamics. Proposes instead of using artificial viscosity to suppress such unphysical solutions as spurious numerical oscillations and nonlinear instabilities, one should formulate equations so that rate of production of entropy within each cell of computational grid be nonnegative, as required by second law.

Merriam, Marshal L.↗
NASA NTRS2004

A Novel Liquid-Liquid Transition in Undercooled Ti-Zr-Ni Liquids

If crystallization can be avoided, liquids enter a metastable (undercooled) state below their equilibrium liquidus temperatures, T(sub l), finally 'freezing' into a glass below a characteristic temperature called the glass transition temperature, T(sub g). In rare cases, the undercooled liquid may undergo a liquid-liquid phase transition (liquid polymorphism) before entering the glassy state. This has been suggested from experimental studies of H2O and Si. Such phase transitions have been predicted in some stable liquids, ie. above T(sub l) at atmospheric pressure, for SiO2 and BeF2, but these have not been verified experimentally. They have been observed in liquids of P, Si and C, but only under high pressure. In this letter we present the first experimental evidence for a phase transition in a low viscosity metallic liquid that is driven by an approach to a constant entropy configuration state and correlated with a growing icosahedral order in the liquid. A maximum in the specific heat at constant pressure, similar to what is normally observed near T(sub g), is reported for undercooled liquids of quasicrystal-forming Ti-Zr-Ni alloys. A two-state excitation model that includes cooperativity by incorporating a temperature-dependent excitation energy, fits the specific heat data well, signaling a phase transition. An inflection in the liquid density with decreasing temperature instead of a discontinuity indicates that this is not a typical first order phase transition; it could be a weakly first order or higher order transition. While showing many similarities to a glass transition, this liquid-liquid phase transition occurs in a mobile liquid, making it novel.

Lee, G. W.↗
NASA NTRS1972

Thermophysical properties of parahydrogen from the freezing liquid line to 5000 R for pressures to 10000 psia

The tables include entropy, enthalpy, internal energy, density, volume, speed of sound, specific heat, thermal conductivity, viscosity, thermal diffusivity, Prandtl number, and the dielectric constant for 65 isobars. Quantities of special utility in heat transfer and thermodynamic calculations are also included in the isobaric tables. In addition to the isobaric tables, tables for the saturated vapor and liquid are given, which include all of the above properties, plus the surface tension. Tables for the P-T of the freezing liquid, index of refraction, and the derived Joule-Thomson inversion curve are also presented.

Mccarty, R. D.↗