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

Theoretical study of production of unique glasses in space

The potential of producing the glassy form of selected materials in the weightless, containerless nature of space processing is examined through the development of kinetic relationships describing nucleation and crystallization phenomena. Transformation kinetics are applied to a well-characterized system (SiO2), an excellent glass former (B2O3), and a poor glass former (Al2O3) by conventional earth processing methods. Viscosity and entropy of fusion are shown to be the primary materials parameters controlling the glass forming tendency. For multicomponent systems diffusion-controlled kinetics and heterogeneous nucleation effects are considered. An analytical empirical approach is used to analyze the mullite system. Results are consistent with experimentally observed data and indicate the promise of mullite as a future space processing candidate.

Larsen, D. C.↗

The development of a CFD potential method for the analysis of tilt-rotors

A new CFD potential code, FPX (eXtended Full-Potential), has been developed for application to both helicopters and tilt-rotors. The code solves the unsteady, three-dimensional full potential equation and is an extension of the rotor code, FPR. Both entropy and viscosity corrections are included to enhance the physical modeling capabilities. A number of efficiency related modifications have yielded a factor of two speed-up in the code. An axial flow capability has been added to treat tilt-rotor in forward flight (cruise mode). In order to employ streamwise periodicity and accurately solve for the propagation of acoustic signals in the tip region, an H-H topology has been added to the basic O-H grid system. Computations are performed for the XV-15 Standard and ATB blades at high-speed conditions. Comparisons are made for the blade aerodynamics and the induced fuselage cabin pressure for a range of Mach numbers. Grid generation, wake treatment, and far-field wall treatment are identified as problem areas with recommendations for future research.

Bridgeman, John O.↗

On the Coupling of CDISC Design Method with FPX Rotor Code

A rotor section aerodynamics design package is developed by coupling Constrained Direct Iterative Surface Curvature (CDISC) design method with the FPX rotor code. The coupling between the CDISC design and the FPX flow analysis is fully automated. The CDISC design method employs a predictor-corrector procedure iteratively to determine a surface geometry which produces a target pressure distribution, where the target pressure distributions is either pre-defined or automatically generated through flow and geometry constraints. The FPX code is an eXtended Full-Potential rotor Computational Fluid Dynamics (CFD) code, which solves the three-dimensional unsteady full-potential equation in a strong conservative form using an implicit approximate factorization finite-difference scheme with entropy and viscosity corrections. Application of the CDISC design method coupled with the FPX rotor code is made for rotor blades in hovering motions. Several design examples are presented to demonstrate the capability of the new package in rotor section design.

Hu, Hong↗

Classical and quantum computing of shear viscosity for ( 2 + 1 ) D SU(2) gauge theory

We perform a nonperturbative calculation of the shear viscosity for ( 2 + 1 )-dimensional SU(2) gauge theory by using the lattice Hamiltonian formulation. The retarded Green’s function of the stress-energy tensor is calculated from real time evolution via exact diagonalization of the lattice Hamiltonian with a local Hilbert space truncation, and the shear viscosity is obtained via the Kubo formula. When taking the continuum limit, we account for the renormalization group flow of the coupling but no additional operator renormalization. We find the ratio of the shear viscosity and the entropy density η s is consistent with a well-known holographic result 1 4 π at several temperatures on a 4 × 4 honeycomb lattice with the local electric representation truncated at j max = 1 2 . We also find the ratio of the spectral function and frequency ρ x y ( ω ) ω exhibits a peak structure when the frequency is small. Both the exact diagonalization method and simple matrix product state classical simulation method beyond j max = 1 2 on bigger lattices require exponentially growing resources. So we develop a quantum computing method to calculate the retarded Green’s function and analyze various systematics of the calculation including j max truncation and finite size effects, Trotter errors and the thermal state preparation efficiency. Our thermal state preparation method still requires resources that grow exponentially with the lattice size, but with a very small prefactor at high temperature. We test our quantum circuit on both the Quantinuum emulator and the IBM simulator for a small lattice and obtain results consistent with the classical computing ones. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Computer programs for thermodynamic and transport properties of hydrogen (tabcode-II)

The thermodynamic and transport properties of para and equilibrium hydrogen have been programmed into a series of computer routines. Input variables are the pair's pressure-temperature and pressure-enthalpy. The programs cover the range from 1 to 5000 psia with temperatures from the triple point to 6000 R or enthalpies from minus 130 BTU/lb to 25,000 BTU/lb. Output variables are enthalpy or temperature, density, entropy, thermal conductivity, viscosity, at constant volume, the heat capacity ratio, and a heat transfer parameter. Property values on the liquid and vapor boundaries are conveniently obtained through two small routines. The programs achieve high speed by using linear interpolation in a grid of precomputed points which define the surface of the property returned.

Roder, H. M.↗

Normal- and oblique-shock flow parameters in equilibrium air including attached-shock solutions for surfaces at angles of attack, sweep, and dihedral

Normal- and oblique-shock flow parameters for air in thermochemical equilibrium are tabulated as a function of shock angle for altitudes ranging from 15.24 km to 91.44 km in increments of 7.62 km at selected hypersonic speeds. Post-shock parameters tabulated include flow-deflection angle, velocity, Mach number, compressibility factor, isentropic exponent, viscosity, Reynolds number, entropy difference, and static pressure, temperature, density, and enthalpy ratios across the shock. A procedure is presented for obtaining oblique-shock flow properties in equilibrium air on surfaces at various angles of attack, sweep, and dihedral by use of the two-dimensional tabulations. Plots of the flow parameters against flow-deflection angle are presented at altitudes of 30.48, 60.96, and 91.44 km for various stream velocities.

Hunt, J. L.↗

Computer program for thermal and transport properties of parahydrogen from 20 to 10,000 K

A computer program was recently developed to provide thermal and transport properties for parahydrogen across a wide temperature and pressure range. The program, NBS+/-pH2, matches the most recent parahydrogen property data from the National Bureau of Standards up to 3000 K and property data from the NASA Lewis Research Center's Chemical Equilibrium Computer Program up to 10,000 K. The pressure range of NBS+/-pH2 is from 1 x 10(exp 4) to 1.6 x 10(exp 7) Pa. The program was developed to meet the need for accurate parahydrogen properties from liquid to dissociated conditions as required by propulsion simulation programs being developed under the Space Exploration Initiative. NBS+/-pH2 is a machine-independent, standard Fortran 77 program which provides density, thermal conductivity, viscosity, Prandtl number, entropy, specific heats, and speed of sound given pressure and either temperature or enthalpy. This program is described and a comparison to programs previously available is provided.

Walton, James T.↗

Transport Coefficients in weakly compressible turbulence

A theory of transport coefficients in weakly compressible turbulence is derived by applying Yoshizawa's two-scale direct interaction approximation to the compressible equations of motion linearized about a state of incompressible turbulence. The result is a generalization of the eddy viscosity representation of incompressible turbulence. In addition to the usual incompressible eddy viscosity, the calculation generates eddy diffusivities for entropy and pressure, and an effective bulk viscosity acting on the mean flow. The compressible fluctuations also generate an effective turbulent mean pressure and corrections to the speed of sound. Finally, a prediction unique to Yoshizawa's two-scale approximation is that terms containing gradients of incompressible turbulence quantities also appear in the mean flow equations. The form these terms take is described.

Rubinstein, Robert↗

A discontinuous Galerkin spectral element method for compressible reacting flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large-eddy simulations because of their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reacting Navier-Stokes equations. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of the DG approach. The framework, implemented in the spectral element code Nek5000, is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. An entropy-residual based artificial viscosity is added to smooth shocked regions of flow, and a positivity-preserving limiter is implemented to suppress non-physical oscillations. These enhancements support the numerical stability of the hydrodynamic sub-step, which is decoupled from the chemistry integration through a second-order operator splitting method. Here, a series of smooth and discontinuous validation cases are presented in increasing physical and computational complexity for both inviscid and viscous flows. In particular, simulations of canonical one-dimensional and two-dimensional detonations are performed, and the high-order numerical results are validated against available literature data. Additional validation studies are carried out for classical three-dimensional numerical simulations of incompressible and compressible turbulent flows.

Compressible reacting flows↗

Azimuthal dependence of two-particle transverse momentum current correlations

Two-particle transverse momentum correlation functions are a powerful technique for understanding the dynamics of relativistic heavy-ion collisions. Among these, the transverse momentum correlator G 2 (Δη,Δφ) is of particular interest for its potential sensitivity to the shear viscosity per unit of entropy density η/s of the quark-gluon plasma formed in heavy-ion collisions. We use the UrQMD, AMPT, and EPOS models for Au–Au at √ s NN = 200 GeV and Pb–Pb at √ s NN = 2760 GeV to investigate the long range azimuthal dependence of G 2 (Δη,Δφ), and explore its utility to constrain η/s based on charged particle correlations. We find that the three models yield quantitatively distinct transverse momentum Fourier harmonics coefficients $a$ $^{pT}_{n}$. We also observe these coefficients exhibit a significant dependence on η/s in the context of the AMPT model. These observations suggest that exhaustive measurements of the dependence of G 2 (Δφ) with collision energy, system size, collision centrality, in particular, offer the potential to distinguish between different theoretical models and their underlying assumptions. Exhaustive analyses of G 2 (Δφ) obtained in large and small systems should also be instrumental in establishing new constraints for precise extraction of η/s.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Perturbation solution of the Navier-Stokes equations and its relation to the Lighthill-Curle solution of aerodynamic sound

The aerodynamic sound described by the Lighthill-Curle solution is reexamined using a method of matched asymptotic expansions. The governing Navier-Stokes equations written in nondimensional form are expanded for a small Mach number. First- and second-order solutions for the pressure field are obtained, and the singular nature of the expansion at large distances is indicated. The nearfield pressure is governed by the Poisson equation, whereas the farfield equations describe a linear wave system in a dissipative medium. The pseudosound is related to the incompressible Reynolds stresses associated with a solenoidal velocity field, the velocity, the pressure perturbation, and their derivatives on the boundaries. A uniformly valid first-order solution for the pressure is obtained. It is shown that viscosity, thermal conductivity, and entropy in the flow do not contribute to the first-order noise generation, while the viscous stress contributes to noise only from some boundaries. The application of the proposed perturbation method to a subsonically moving surface and a hot jet is discussed.

Pan, Y. S.↗

Entropy conservative finite element schemes

The question of entropy stability for discrete approximations to hyperbolic systems of conservation laws is studied. The amount of numerical viscosity present in such schemes is quantified and related to their entropy stability by means of comparison. To this end, two main ingredients are used: entropy variables and the construction of certain entropy conservative schemes in terms of piecewise-linear finite element approximations. It is then shown that conservative schemes are entropy stable, if and (for three-point schemes) only if, they contain more numerical viscosity than the abovementioned entropy conservation ones.

Tadmor, E.↗

Energy considerations in computational aeroacoustics

A finite-volume multistage time-stepping Euler code is used to investigate the use of CFD algorithms for the direct calculation of acoustics. The 2D compressible inviscid flow about an accelerating or decelerating circular cylinder is used as a model problem. The time evolution of the energy transfer from the cylinder to the fluid, as the cylinder is moved from rest to some nonnegligible velocity, is clearly seen. By examining the temporal and spatial characteristics of the numerical solution, a distinction can be made between the propagating acoustic energy, the convecting energy associated with the entropy change in the fluid, and the energy contained in the local aerodynamic field. Systematic variation of the cylinder acceleration shows that the radiated acoustic energy depends strongly upon the rate of acceleration or deceleration. The computational grid has a large effect on the ratio of acoustic energy to nonphysical entropy associated energy, while the role of the explicit artificial viscosity seems to be of second order. The entropy term was nearly negligible in all cases the cylinder was started slowly.

Brentner, Kenneth S.↗

WASP: A flexible FORTRAN 4 computer code for calculating water and steam properties

A FORTRAN 4 subprogram, WASP, was developed to calculate the thermodynamic and transport properties of water and steam. The temperature range is from the triple point to 1750 K, and the pressure range is from 0.1 to 100 MN/m2 (1 to 1000 bars) for the thermodynamic properties and to 50 MN/m2 (500 bars) for thermal conductivity and to 80 MN/m2 (800 bars) for viscosity. WASP accepts any two of pressure, temperature, and density as input conditions. In addition, pressure and either entropy or enthalpy are also allowable input variables. This flexibility is especially useful in cycle analysis. The properties available in any combination as output include temperature, density, pressure, entropy, enthalpy, specific heats, sonic velocity, viscosity, thermal conductivity, surface tension, and the Laplace constant. The subroutine structure is modular so that the user can choose only those subroutines necessary to his calculations. Metastable calculations can also be made by using WASP.

Hendricks, R. C.↗

Arrhenian to non-Arrhenian crossover in glass melt viscosity

The activation energy for glass melt viscosity is a function of temperature, i.e., non-Arrhenian, within a middle range of viscosity values (typically, between 10 2 and 10 10 Pa s) and constant, i.e., Arrhenian, outside this range. At the low-viscosity end, the Arrhenian to non-Arrhenian transition (termed the crossover) falls within the glass processing range: glass melting and fining occur within the Arrhenian range while glass forming occurs within the non-Arrhenian range. By the Adam-Gibbs equation, the configuration entropy is nearly constant at high temperatures (low viscosities), where the glass melt turns into a simple liquid, and decreases with decreasing temperature (increasing viscosity) as the glass is increasingly polymerized. In the present study, the power-law function and hyperbolic tangent function are used to express the configuration entropy as a function of temperature. The viscosity crossover is defined as the point at which onset of structural grouping is causing the configuration entropy to deviate from its high-temperature value. In conclusion, it appears that the limited range of the configuration entropy of the fully depolymerized melt imposes a restriction on the crossover viscosity; for well-defined glass families, such as float glasses or nuclear waste glasses, the crossover point can be defined using a fixed, i.e., composition-independent, viscosity value.

36 MATERIALS SCIENCE↗

Optimization of artificial viscosity in production codes based on Gaussian Regression surrogate models

To accurately model flows with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity has to be introduced to convert kinetic energy into internal energy, thereby increasing the entropy across shocks. Determining the appropriate strength of the artificial viscosity is an art and strongly depends on the particular problem and experience of the researcher. The objective of this study is to pose the problem of finding the appropriate strength of artificial viscosity as an optimization problem and solve this problem using machine learning (ML) tools, specifically using surrogate models based on Gaussian Process regression and Bayesian analysis. We describe the optimization method and discuss various practical details of its implementation. The shock-containing problems for which we apply this method all have been implemented in the LANL code FLAG. First, we apply ML to find optimal values to isolated shock problems of different strengths. Second, we apply ML to optimize viscosity for a 1D propagating detonation problem based on Zel’dovich-von Neumann-Doring (ZND) detonation theory using a reactive burn model. We compare results for default (currently used values in FLAG) and optimized values of artificial viscosity for these problems demonstrating the potential for significant improvement in the accuracy of computations.

42 ENGINEERING↗

Computer program for calculating water and steam properties

Computer subprogram calculates thermodynamic and transport properties of water and steam. Program accepts any two of pressure, temperature, and density as input conditions. Pressure and either entropy or enthalpy are also allowable input variables. Output includes any combination of temperature, density, pressure, entropy, enthalpy, specific heats, sonic velocity, viscosity, thermal conductivity, surface tension, and the Laplace constant.

Hendricks, R. C.↗