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 37 records · Page 2

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.↗

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↗

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↗

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.↗

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↗

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.↗

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.↗

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.↗

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.↗

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↗

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↗

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.↗

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.↗

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.↗

The entropy rate admissibility criterion for solution of hyperbolic conservation laws

The entropy rate admissibility criterion for solutions of hyperbolic conservation laws is numerically analyzed. The following admissibility criterion for solutions of hyperbolic conservation laws is proposed: a weak solution is admissible if the total entropy decays with the highest possible rate. The equivalence of this criterion and viscosity criterion is established for the single equation and the system of equations of one dimensional nonlinear elasticity.

Dafermos, C. M.↗

Thermodynamic and transport combustion properties of hydrocarbons with air. Part 1: Properties in SI units

Thermodynamic and transport combustion properties were calculated for a wide range of conditions for the reaction of hydrocarbons with air. Three hydrogen-carbon atom ratios (H/C = 1.7, 2.0, 2.1) were selected to represent the range of aircraft fuels. For each of these H/C ratios, combustion properties were calculated for the following conditions: Equivalence ratio: 0, 0.25, 0.5, 0.75, 1.0, 1.25 Water - dry air mass ratio: 0, 0.03 Pressure, kPa: 1.01325, 10.1325, 101.325, 1013.25, 5066.25 (or in atm: 0.01, 0.1, 1, 10, 50) Temperature, K: every 10 degrees from 200 to 900 K; every 50 degrees from 900 to 3000 K Temperature, R: every 20 degrees from 360 to 1600 R; very 100 degrees from 1600 to 5400 R. The properties presented are composition, density, molecular weight, enthalphy, entropy, specific heat at constant pressure, volume derivatives, isentropic exponent, velocity of sound, viscosity, thermal conductivity, and Prandtl number. Property tables are based on composites that were calculated by assuming both: (1) chemical equilibrium (for both homogeneous and heterogeneous phases) and (2) constant compositions for all temperatures. Properties in SI units are presented in this report for the Kelvin temperature schedules.

Gordon, S.↗

On the convergence of difference approximations to scalar conservation laws

A unified treatment is given for time-explicit, two-level, second-order-resolution (SOR), total-variation-diminishing (TVD) approximations to scalar conservation laws. The schemes are assumed only to have conservation form and incremental form. A modified flux and a viscosity coefficient are introduced to obtain results in terms of the latter. The existence of a cell entropy inequality is discussed, and such an equality for all entropies is shown to imply that the scheme is an E scheme on monotone (actually more general) data, hence at most only first-order accurate in general. Convergence for TVD-SOR schemes approximating convex or concave conservation laws is shown by enforcing a single discrete entropy inequality.

Osher, Stanley↗

An artificial energy method for calculating flows with shocks

The artificial-viscosity method, first proposed by von Neumann and Richtmyer, introduces an artificial viscous pressure term in regions of compression such that an increase in entropy occurs in shock transition zones. The paper describes how dissipative flows can be induced by reducing the total energy available for adiabatic processes in shock zones. A class of inviscid fluid flows, called semiflows, is described in which the flows exhibit thermodynamic differences. Induced dissipative flows modify the pressure in regions of compression in a manner analogous to the artificial-viscosity method and for a gas, the effect is equivalent to suitably modifying the gas constant in the equation of state. By employing MacCormack's method and the usual non-adiabatic equations, numerical solutions of a Riemann problem are compared with the modified artificial energy method, showing that the dissipation effect predicted by the analytical formulation is reflected in the numerical method as well.

Rose, M. E.↗