SEARCH · Search NASA

Results for “Entropy Viscosity Method”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2015

Leading-order anisotropic hydrodynamics for central collisions

In this work, we use leading-order anisotropic hydrodynamics to study an azimuthally symmetric boost-invariant quark-gluon plasma. We impose a realistic lattice-based equation of state and perform self-consistent anisotropic freeze-out to hadronic degrees of freedom. We then compare our results for the full spatiotemporal evolution of the quark-gluon plasma and its subsequent freeze-out to results obtained using 1+1D Israel-Stewart second-order viscous hydrodynamics. We find that for small shear viscosities, 4πη/s ~1, the two methods agree well for nucleus-nucleus collisions; however, for large-shear-viscosity-to-entropy-density ratios or proton-nucleus collisions we find important corrections to the Israel-Stewart results for the final particle spectra and the total number of charged particles. Finally, we demonstrate that the total number of charged particles produced is a monotonically increasing function of 4πη/s in Israel-Stewart viscous hydrodynamics, whereas in anisotropic hydrodynamics it has a maximum at 4πη/s ~10. For all 4πη/s > 0, we find that for Pb-Pb collisions Israel-Stewart viscous hydrodynamics predicts more dissipative particle production than anisotropic hydrodynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1987

The numerical viscosity of entropy stable schemes for systems of conservation laws. I

Discrete approximations to hyperbolic systems of conservation laws are 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, conservative schemes which are also entropy-conservative are constructed. These entropy-conservative schemes enjoy second-order accuracy; moreover, they can be interpreted as piecewise-linear finite-element methods, and hence can be formulated on various mesh configurations. It is then shown that conservative schemes are entropy stable, if and (for three-point schemes) only they contain more viscosity than that present in the above-mentioned entropy-conservative ones.

Tadmor, Eitan↗

DOE OSTI2017

CosmosDG: An hp -adaptive Discontinuous Galerkin Code for Hyper-resolved Relativistic MHD

We have extended Cosmos++, a multidimensional unstructured adaptive mesh code for solving the covariant Newtonian and general relativistic radiation magnetohydrodynamic (MHD) equations, to accommodate both discrete finite volume and arbitrarily high-order finite element structures. The new finite element implementation, called CosmosDG, is based on a discontinuous Galerkin (DG) formulation, using both entropy-based artificial viscosity and slope limiting procedures for the regularization of shocks. High-order multistage forward Euler and strong-stability preserving Runge–Kutta time integration options complement high-order spatial discretization. We have also added flexibility in the code infrastructure allowing for both adaptive mesh and adaptive basis order refinement to be performed separately or simultaneously in a local (cell-by-cell) manner. We discuss in this report the DG formulation and present tests demonstrating the robustness, accuracy, and convergence of our numerical methods applied to special and general relativistic MHD, although we note that an equivalent capability currently also exists in CosmosDG for Newtonian systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

Viscous anisotropic hydrodynamics for the Gubser flow

In this work we describe the dynamics of a highly anisotropic system undergoing boost-invariant longitudinal and azimuthally symmetric radial expansion (Gubser flow) for arbitrary shear viscosity to entropy density ratio. We derive the equations of motion of dissipative anisotropic hydrodynamics by applying to this situation the moments method recently derived by Moln´ar et al. (MNR) [1, 2], based on an expansion around an arbitrary anisotropic one-particle distribution function. One requires an additional evolution equation in order to close the conservation laws. This is achieved by selecting the relaxation equation for the longitudinal pressure with a suitable Landau matching condition. As a result one obtains two coupled differential equations for the energy density and the longitudinal pressure which respect the SO (3)q Ⓧ SO (1, 1) Ⓧ Z 2 symmetry of the Gubser flow in the deSitter space. These equations are solved numerically and compared with the predictions of the recently found exact solution of the relaxation-time-approximation Boltzmann equation subject to the same flow. We also compare our numerical results with other fluid dynamical models. As a result we observe that the MNR description of anisotropic fluid dynamics reproduces the space-time evolution of the system than all other currently known hydrodynamical approaches.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1982

Design of supercritical swept wings

Computational fluid dynamics are used to discuss problems inherent to transonic three-dimensional flow past supercritical swept wings. The formulation for a boundary value problem for the flow past the wing is provided, including consideration of weak shock waves and the use of parabolic coordinates. A swept wing code is developed which requires a mesh of 152 x 10 x 12 points and 200 time cycles. A formula for wave drag is calculated, based on the idea that the conservation form of the momentum equation becomes an entropy inequality measuring the drag, expressible in terms of a small-disturbance equation for a potential function in two dimensions. The entropy inequality has been incorporated in a two-dimensional code for the analysis of transonic flow over airfoils. A method of artificial viscosity is explored for optimum pressure distributions with design, and involves a free boundary problem considering speed over only a portion of the wing.

Garabedian, P.↗

NASA NTRS1988

Application of a nonisentropic full potential method to AGARD standard airfoils

An entropy-correction method for the unsteady full potential equation is presented. The unsteady potential equation is modified to model the entropy jumps across shock waves. The conservative form of the modified equation is solved in generalized coordinates using an implicit, approximate factorization method. A flux-biasing differencing method, which generates the proper amounts of artificial viscosity in supersonic regions, is used to discretize the flow equations in space. Calculated results are presented for the NLR 7301, NACA 0012, and NACA 64A010A airfoils. Comparisons of the present method and solutions of the Euler equations are presented for the NLR 7301 airfoil, and comparisons of the present method and experimental data are presented for all three airfoils. The comparisons show that the present method more accurately models solutions of the Euler equations and experiment than does the isentropic potential formulation. In addition, it is shown that modeling shock-generated entropy extends the range of validity of the full potential method.

Whitlow, Woodrow, Jr.↗

NASA NTRS1988

Application of a nonisentropic full potential method to AGARD standard airfoils

Whitlow, Woodrow, Jr.↗

DOE OSTI2018

Viscosity models for pure hydrocarbons at extreme conditions: A review and comparative study

Here, viscosity is a critical fundamental property required in many applications in the chemical and oil industries. In this review the performance of seven select viscosity models, representative of various predictive and correlative approaches, is discussed and evaluated by comparison to experimental data of 52 pure hydrocarbons including straight-chain alkanes, branched alkanes, cycloalkanes, and aromatics. This analysis considers viscosity data to extremely high-temperature, high-pressure conditions up to 573 K and 300 MPa. Unsatisfactory results are found, particularly at high pressures, with the Chung-Ajlan-Lee-Starling, Pedersen-Fredenslund, and Lohrenz-Bray-Clark models commonly used for oil reservoir simulation. If sufficient experimental viscosity data are readily available to determine model-specific parameters, the free volume theory and the expanded fluid theory models provide generally comparable results that are superior to those obtained with the friction theory, particularly at pressures higher than 100 MPa. Otherwise, the entropy scaling method by Lötgering-Lin and Gross is recommended as the best predictive model.

02 PETROLEUM↗

DOE OSTI2025

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↗

DOE OSTI2016

Stokes–Einstein relation and excess entropy in Al-rich Al-Cu melts

We investigate the conditions for the validity of the Stokes-Einstein relation that connects diffusivity to viscosity in melts using entropy-scaling relationships developed by Rosenfeld. Employing ab initio molecular dynamics simulations to determine transport and structural properties of liquid Al{sub 1−x}Cu{sub x} alloys (with composition x ≤ 0.4), we first show that reduced self-diffusion coefficients and viscosities, according to Rosenfeld's formulation, scale with the two-body approximation of the excess entropy except the reduced viscosity for x = 0.4. Then, we use our findings to evidence that the Stokes-Einstein relation using effective atomic radii is not valid in these alloys while its validity can be related to the temperature dependence of the partial pair-excess entropies of both components. Finally, we derive a relation between the ratio of the self-diffusivities of the components and the ratio of their pair excess entropies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1975

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

DOE OSTI2014

Traveling waves in Hall-magnetohydrodynamics and the ion-acoustic shock structure

Hall-magnetohydrodynamics (HMHD) is a mixed hyperbolic-parabolic partial differential equation that describes the dynamics of an ideal two fluid plasma with massless electrons. We study the only shock wave family that exists in this system (the other discontinuities being contact discontinuities and not shocks). We study planar traveling wave solutions and we find solutions with discontinuities in the hydrodynamic variables, which arise due to the presence of real characteristics in Hall-MHD. We introduce a small viscosity into the equations and use the method of matched asymptotic expansions to show that solutions with a discontinuity satisfying the Rankine-Hugoniot conditions and also an entropy condition have continuous shock structures. The lowest order inner equations reduce to the compressible Navier-Stokes equations, plus an equation which implies the constancy of the magnetic field inside the shock structure. We are able to show that the current is discontinuous across the shock, even as the magnetic field is continuous, and that the lowest order outer equations, which are the equations for traveling waves in inviscid Hall-MHD, are exactly integrable. We show that the inner and outer solutions match, which allows us to construct a family of uniformly valid continuous composite solutions that become discontinuous when the diffusivity vanishes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1994

Fourth order difference methods for hyperbolic IBVP's

Fourth order difference approximations of initial-boundary value problems for hyperbolic partial differential equations are considered. We use the method of lines approach with both explicit and compact implicit difference operators in space. The explicit operator satisfies an energy estimate leading to strict stability. For the implicit operator we develop boundary conditions and give a complete proof of strong stability using the Laplace transform technique. We also present numerical experiments for the linear advection equation and Burgers' equation with discontinuities in the solution or in its derivative. The first equation is used for modeling contact discontinuities in fluid dynamics, the second one for modeling shocks and rarefaction waves. The time discretization is done with a third order Runge-Kutta TVD method. For solutions with discontinuities in the solution itself we add a filter based on second order viscosity. In case of the non-linear Burger's equation we use a flux splitting technique that results in an energy estimate for certain different approximations, in which case also an entropy condition is fulfilled. In particular we shall demonstrate that the unsplit conservative form produces a non-physical shock instead of the physically correct rarefaction wave. In the numerical experiments we compare our fourth order methods with a standard second order one and with a third order TVD-method. The results show that the fourth order methods are the only ones that give good results for all the considered test problems.

Gustafsson, Bertil↗

NASA NTRS1986

An entropy correction method for unsteady full potential flows with strong shocks

An entropy correction method for the unsteady full potential equation is presented. The unsteady potential equation is modified to account for entropy jumps across shock waves. The conservative form of the modified equation is solved in generalized coordinates using an implicit, approximate factorization method. A flux-biasing differencing method, which generates the proper amounts of artificial viscosity in supersonic regions, is used to discretize the flow equations in space. Comparisons between the present method and solutions of the Euler equations and between the present method and experimental data are presented. The comparisons show that the present method more accurately models solutions of the Euler equations and experiment than does the isentropic potential formulation.

Whitlow, W., Jr.↗

NASA NTRS1986

An entropy correction method for unsteady full potential flows with strong shocks

Whitlow, W., Jr.↗

DOE OSTI2023

Turbulence Modeling with Nek5000/RS, SOD2D and Alya

We present Validation and Verification (V & V) study of two high-order spectral element method (SEM) based Computational Fluid Dyanimcs (CFD) codes that we will utilize for work on turbulence modeling: Nek5000 and SOD2D. While the former solves the incompressible form of Navier-Stokes equation, the latter works with compressible set of equations and uses an entropy-viscosity formulation to account for the discontinuities for high Mach number flows. We demonstrate the accuracy of these codes for two benchmark sub-sonic turbulent flows: periodic channel and pipe flow, by carrying out first and second order statistical analysis, grid convergence and turbulent structure analysis using wall-resolved large eddy simulations (WRLES). Later, we report on the implementation and testing of various wall-modeling strategies for large eddy simulation of turbulent flows in Nek5000. These include both classical log-law based and decision-tree based machine-learning models. Accuracy of these closure strategies are analyzed and necessary future work is outlined.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1991

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

NASA NTRS2002

Entropy Splitting for High Order Numerical Simulation of Compressible Turbulence

A stable high-order numerical scheme for direct numerical simulation (DNS) of shock-free compressible turbulence is presented. The method is applicable to general geometries. It contains no upwinding, artificial dissipation, or filtering. Instead the method relies on the stabilizing mechanisms of an appropriate conditioning of the governing equations and the use of compatible spatial difference operators for the interior points (interior scheme) as well as the boundary points (boundary scheme). An entropy-splitting approach splits the inviscid flux derivatives into conservative and nonconservative portions. The spatial difference operators satisfy a summation-by-parts condition, leading to a stable scheme (combined interior and boundary schemes) for the initial boundary value problem using a generalized energy estimate. A Laplacian formulation of the viscous and heat conduction terms on the right hand side of the Navier–Stokes equations is used to ensure that any tendency to odd–even decoupling associated with central schemes can be countered by the fluid viscosity. The resulting methods are able to minimize the spurious high-frequency oscillations associated with pure central schemes, especially for long time integration applications such as DNS. For validation purposes, the methods are tested in a DNS of compressible turbulent plane channel flow at low values of friction Mach number, where reference turbulence data bases exist. It is demonstrated that the methods are robust in terms of grid resolution, and in good agreement with published channel data. Accurate turbulence statistics can be obtained with moderate grid sizes. Stability limits on the range of the splitting parameter are determined from numerical tests.

N.D. Sandham↗