The conservation equations for simulating hypersonic flows in thermal and chemical nonequilibrium and details of the associated physical models are presented These details include the curve fits used for defining thermodynamic properties of the 11 species air model, curve fits for collision cross sections, expressions for transport properties, the chemical kinetics models, and the vibrational and electronic energy relaxation models The expressions are formulated in the context of either a two or three temperature model Greater emphasis is placed on the two temperature model in which it is assumed that the translational and rotational energy models are in equilibrium at the translational temperature, T, and the vibrational, electronic, and electron translational energy modes are in equilibrium at the vibrational temperature, T sub v The eigenvalues and eigenvectors associated with the Jacobian of the flux vector are also presented in order to accommodate the upwind based numerical solutions of the complete equation set

/pdf/conservation-equations-and-physical-models-for-hypersonic-194duqack2.pdf

Conservation equations and physical models for hypersonic air flows in thermal and chemical nonequilibrium

An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for three-dimensional, viscous, compressible, perfect-gas flows is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. Viscous terms are discretized using central differences. The relaxation strategy is well suited for computers employing either vector or parallel architectures. It is also well suited to the numerical solution of the governing equations on unstructured grids. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large, block tri-diagonal systems. Convergence rates and grid refinement studies are conducted for Mach 5 flow through an inlet with a 10 deg compression ramp and Mach 14 flow over a 15 deg ramp. Predictions for pressure distributions, surface heating, and aerodynamics coefficients compare well with experiment data for Mach 10 flow over a blunt body.

https://ntrs.nasa.gov/api/citations/19900007726/downloads/19900007726.pdf

An upwind-biased, point-implicit relaxation algorithm for viscous, compressible perfect-gas flows

Blunt-body flow simulations

An upwind-biased point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for three-dimensional viscous hypersonic flows in chemical and thermal nonequilibrium is described. Details of the algorithm development, in the context of an 11-species two-temperature reacting gas model, are emphasized. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large block-tridiagonal systems. Predictions for the hypersonic flow of air in chemical and thermal nonequilibrium (velocity = 8917 m/s, altitude = 78 km) over the Aeroassist Flight Experiment configuration, obtained on a multidomain grid, are discussed.

Upwind-biased, point-implicit relaxation strategies for viscous, hypersonic flows

The code calibration program for the Langley Aerothermodynamic Upwind Relaxation Algorithm to be used as support for the Aeroassist Flight Experiment (AFE) is discussed. Comparisons between experimental data and numerical simulations are made which focus on perfect-gas tests involving a scale model of the AFE. Aspects of the thermochemical nonequilibrium model are called into question by the results of ground tests performed in a ballistic range and in a shock tunnel.

Code calibration program in support of the Aeroassist Flight Experiment

Program LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is a robust, finite volume, single-level storage, implicit upwind differencing algorithm which has been documented and tested on several three-dimensional blunt-body flows. The algorithm can run at unlimited Courant numbers (relaxing the steady-state equations) but requires the inversion of only a 5 x 5 matrix per computational cell. An alternating directional sweep Gauss-Seidel substitution strategy is used to relax the governing equations. At present, the Euler and thin-layer Navier-Stokes equations using Sutherland's law for viscosity have been modeled for a perfect gas, equilibrium air, and nonequilibrium air chemistry neglecting diffusion. The equilibrium and nonequilibrium air chemistry options have been described in a companion paper. Good comparisons with experimental data and another calculation method for pressure distributions, aerodynamic coefficients, and heat-transfer distributions have been demonstrated for three-dimensional blunt-body flows.

P. A. Gnoffo

Papers

Application of program LAURA to three-dimensional AOTV flowfields