Survey of simulation and diagnostic techniques for hypersonic nonequilibrium flows
Explore the source record for details and available documents.
Engineering topics
Publications and source records attributed to Park, Chul.
Explore the source record for details and available documents.
An aerobraking orbital transfer vehicle (AOTV) (80) has aerobrake (82) with a blunted raked-off circular-cone configuration. The other components of the AOTV, including command/control module (95), fuel tanks (86, 88, 89 and 91), rocket engines (94) and afterbody (84), are positioned substantially along resultant force axis (104) of the AOTV (80). The axis (104) coincides with the resultant (sum of lift and drag) force vector. Afterbody (84) is mounted behind the aerobrake (82) with its length extending rearwardly from the aerobrake. The base flow clearance angle .phi. of the aerobrake (80) is 25.degree., thus allowing the afterbody (84) to extend rearwardly from the aerobrake (82) to a much greater extent than possible with a raked-off elliptic-cone aerobraking shield configuration. Afterbody size limitation and other problems associated with the raked-off elliptic-cone aerobraking shield configuration are alleviated by the combination of the aerobrake shape and positioning of the fuel tanks (86, 88, 89 and 91), rocket engines (94) and afterbody (84).
The computational equations governing the shock-layer flowfield over a blunt body dominated by radiation and ablation are derived. The flowfield considered is that expected around the vehicles entering the atmosphere of the planet Mars or the earth on return from Mars on a sprint-type manned mission. To determine the appropriate method of approach, the radiative transfer calculations are made first using the existing computer code RASLE (Nicolet et al., 1978), which uses the assumption of equilibrium and a coarse spectral band model for radiative transfer, and the code NONEQ which accounts for thermochemical nonequilibrium and carries out line-by-line calculation of radiation. By comparing the results, it is concluded that a thermochemical nonequilibrium model and a line-by-line spectral model must be used.
The experimental program for validating real gas hypersonic flow codes at NASA Ames Rsearch Center is described. Ground-based test facilities used include ballistic ranges, shock tubes and shock tunnels, arc jet facilities and heated-air hypersonic wind tunnels. Also included are large-scale computer systems for kinetic theory simulations and benchmark code solutions. Flight tests consist of the Aeroassist Flight Experiment, the Space Shuttle, Project Fire 2, and planetary probes such as Galileo, Pioneer Venus, and PAET.
Explore the source record for details and available documents.
The theoretical basis for the Landau-Teller equation commonly used to model rotational nonequilibrium is reviewed. Several assumptions underlying this model are indicated to be unrealistic for rarefied hypersonic flow. A new rotational nonequilibrium model based on recent measurements up to 2000K of quantum state transition rates is presented. The new model is applied to the continuum study of normal shock wave structure using both the Burnett equations and a simplified nonlinear constitutive relation. Comparisons are made between shock wave temperature profiles generated from the new model and those generated using the Landau-Teller model. Comparisons of shock reciprocal thicknesses between experimental data for nitrogen and continuum solutions using both rotational models are made. The new rotational model agrees well with experiment up to Mach 6, and under predicts shock thickness at higher Mach numbers. The Landau-Teller model agrees well with experimental shock thickness up to Mach 3, and over predicts the thickness at higher Mach numbers. A modification to the rotational collision number in the Landau-Teller model is found to give results which agree with experimental shock thicknesses at all Mach numbers up to 11.
The existing experimental data on the rate coefficients for the chemical reactions in nonequilibrium high temperature air are reviewed and collated, and a selected set of such values is recommended for use in hypersonic flow calculations. For the reactions of neutral species, the recommended values are chosen from the experimental data that existed mostly prior to 1970, and are slightly different from those used previously. For the reactions involving ions, the recommended rate coefficients are newly chosen from the experimental data obtained more recently. The reacting environment is assumed to lack thermal equilibrium, and the rate coefficients are expressed as a function of the controlling temperature, incorporating the recent multitemperature reaction concept.
The Cassini Mission, a NASA/ESA cooperative project which includes a deployment of probe into the atmosphere of Titan, is described, with particular attention given to the shock radiometer experiment planned for the Titan probe for the analysis of Titan's atmosphere. Results from a shock layer analysis are presented, demonstrating that the mole fractions of the major species (N2, CH4, and, possibly Ar) in the Titan atmosphere can be successfully determined by the Titan-probe radiometer, by measuring the intensity of the CN(violet) radiation emitted in the shock layer during the high velocity portion of the probe entry between 200 and 400 km altitude. It is shown that the sensitivity of the CN(violet) radiation makes it possible to determine the mole fractions of N2, CH4, and Ar to about 0.015, 0.003, and 0.01, respectively, i.e., much better than the present uncertainties in the composition of Titan atmosphere.
The aerothermodynamic problems associated with the aerobraking of the spacecraft proposed for the manned Mars mission are studied. The propulsive Delta V necessary at departure from earth and Mars and the velocities of the atmospheric entries into the two planets are deduced. It is shown that the propulsive Delta V can be reduced by increasing the entry velocities and that entry velocities up to about 15 km/sec are appropriate at both earth and Mars. L/D values of 0.8 and 2.0 are found to be necessary at earth and Mars, respectively. Density, pressure, and stagnation-point convective-heat-transfer rates are calculated for the typical aerobraking flights. Assuming the shock layer flow to be in equilibrium, the stagnation-point radiative-heat-transfer rates are calculated to be larger than the convective-heat-transfer rates. The possible impact of ablation, turbulence, and nonequilibrium are discussed.
The effect of vibrational excitation and dissociation at high temperatures on the trim angle of attack of a blunt lifting body is calculated for a nonequilibrium flow regime in air using a CFD technique. The vibrational-electronic temperature and the species densities are calculated assuming the flow to be in a nonequilibrium state. The forebody flow of a two-dimensional blunt body of the shape of the Apollo Command Module at a finite angle of attack is calculated. The results show that the pitching moment around a reference point is larger and the trim angle of attack is smaller for a reacting gas than for a perfect gas. The calculated shift in the trim angle due to the real-gas effect is of the same order as that seen during the Apollo flights.
The importance of thermochemical nonequilibrium in hypersonic flow is studied as well as the breakdown of a one-temperature kinetic model. Consideration is given to vibration relaxation, reaction rates, and conservation equations for chemical variables. Experimental data on the pitching moment characteristics of the raked-off blunted elliptic cone geometry proposed for use with an aeroassisted orbit transfer vehicle are presented.
A CFD technique is described in which the finite-rate chemistry in thermal and chemical nonequilibrium air is fully and implicitly coupled with the fluid motion. Developed for use in the suborbital hypersonic flight speed range, the method accounts for nonequilibrium vibrational and electronic excitation and dissociation, but not ionization. The steady-state solution to the resulting system of equations is obtained by using a lower-upper factorization and symmetric Gauss-Seidel sweeping technique through Newton iteration. Inversion of the left-hand-side matrices is replaced by scalar multiplications through the use of the diagonal dominance algorithm. The code, named CENS2H (Compressible-Euler-Navier-Stokes Two-Dimensional Hypersonic), is fully vectorized and requires about 8.8 x 10 to the -5th sec per node point per iteration using a Cray X-MP computer. Converged solutions are obtained after about 2400 iterations. Sample calculations are made for a circular cylinder and a 10 percent airfoil at 5 deg angle of attack. The calculated cylinder flow field agrees with that obtained experimentally. The code predicts a 10 percent change in lift, drag, and pitching moment for the airfoil due to the thermochemical phenomena.
Nonequilibrium phenomena in hypersonic flows are examined on the basis of theoretical models and selected experimental data, in an introduction intended for second-year graduate students of aerospace engineering. Chapters are devoted to the physical nature of gas atoms and molecules, transitions of internal states, the formulation of the master equation of aerothermodynamics, the conservation equations, chemical reactions in CFD, the behavior of air flows in nonequilibrium, experimental aspects of nonequilibrium flow, a review of experimental results, and gas-solid interaction. Diagrams, graphs, and tables of numerical data are provided.
The experimental program for validating real gas hypersonic flow codes at NASA Ames is described. Ground based test facilities used include ballistic ranges, shock tubes and shock tunnels, arcjet facilities and heated air hypersonic wind tunnels. Also included are large scale computer systems for kinetic theory simulations and benchmark code solutions. Flight tests consist of the Aeroassist Flight Experiment, the Space Shuttle, Project Fire 2, and planetary probes such as Galileo, Pioneer Venus and PAET.
The results of the solution of the equations that describe a hypersonic ionized flow about an elliptically blunted cone are presented. The flow conditions correspond to those of the proposed Aeroassist Flight Experiment (AFE) vehicle at altitudes between the perigee at 78 km and the approximate limit of the continuum regime at 90 km. For the free-stream velocities of interest, about 9 km/sec, the flowfield is out of thermo-chemical equilibrium, electronically excited, ionized and radiating. The gas consists of eight-chemical species including free electrons. The thermal state of the gas is modeled with a translational-rotational temperature, four vibrational temperatures for the diatomic species and an electron-electronic temperature. The electronic excitation of molecules is included. The nonequilibrium air radiation from each fluid element is computed and the radiative heat flux at the body surface is determined. The stagnation point radiative heating result agrees with previous calculations.
A theoretical study of vibrational excitations and dissociations of nitrogen undergoing a nonequilibrium relaxation process upon heating and cooling is reported. The rate coefficients for collisional induced vibrational transitions and transitions from a bound vibrational state into a dissociative state have been calculated using an extension of the theory originally proposed by Schwarz (SSH) et al. (1952). High-lying vibrational states and dissociative states were explicitly included but rotational energy transfer was neglected. The transition probabilities calculated from the SSH theory were fed into the master equation, which was integrated numerically to determine the population distribution of the vibrational states as well as bulk thermodynamic properties. The results show that: (1) the transition rates have a minimum near the middle of the bound vibrational levels, causing a bottleneck in the vibrational relaxation and dissociation rates; (2) high vibrational states are always in equilibrium with the dissociative state; (3) for the heating case, only the low vibrational states relax according to the Landau-Teller theory; (4) for the cooling case, vibrational relaxation cannot be described by a rate equation; (5) Park's (1985, 1988) two-temperature model is approximately valid; and (6) the average vibrational energy removed in dissociation is about 30 percent of the dissociation energy.
Explore the source record for details and available documents.
Explore the source record for details and available documents.