Solution strategy for three-dimensional configurations at hypersonic speeds
Explore the source record for details and available documents.
Engineering topics
Publications and source records attributed to Gnoffo, Peter A..
Explore the source record for details and available documents.
The numerical simulation of hypersonic flow in chemical nonequilibrium over cylindrically blunted 6 degree wedge is described. The simulation was executed on a Cray C-90 with Program LAURA 92-vl. Code setup procedures and sample results, including grid refinement studies and variations of species number are discussed. This simulation relates to a study of wing leading edge heating on transatmospheric vehicles.
A multiblock, laminar heating analysis for the shuttle orbiter at three trajectory points ranging from Mach 24.3 to Mach 12.86 on re-entry is described. The analysis is performed using the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) with both a seven species chemical nonequilibrium model and an equilibrium model. A finite-catalytic-wall model appropriate for shuttle tiles at a radiative equilibrium wall temperature is applied. Computed heating levels are generally in good agreement with the flight data though a few rather large discrepancies remain unexplained. The multiblock relaxation strategy partitions the flowfield into manageable blocks requiring a fraction of the computational resources (time and memory) required by a full domain approach. In hot, the computational cost for a solution at even a single trajectory point would be prohibitively expensive at the given resolution without the multiblock approach. Converged blocks are reassembled to enable a fully coupled converged solution over the entire vehicle, starting from a nearly converged initial condition.
An analysis of the longitudinal aerodynamics of the Shuttle Orbiter in the hypersonic flight regime is made through the use of computational fluid dynamics (CFD). Particular attention is given to establishing the cause of the 'pitching moment anomaly' which occurred on the Orbiter's first flight and to computing the aerodynamics of a complete Orbiter configuration at flight conditions. Data from ground based facilities as well as Orbiter flight data are used to validate the computed results. Analysis shows that the 'pitching moment anomaly' is a real gas chemistry effect which cannot be simulated in ground-based facilities. Computed flight aerodynamics for the Orbiter are within 5 percent of the measured flight values and trim bodyflap deflections are predicted to within 10 percent.
A heating analysis for a 15.2 m diameter Lunar Transfer Vehicle (LTV) at 0 and 10.6 deg angle of attack for a nominal trajectory through the earth's atmosphere is described. The analysis utilizes the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) with thin-layer, Navier-Stokes, thermochemical nonequilibrium options. Radiative heating levels are calculated using the Langley Optimized RAdiative Nonequilibrium (LORAN) and the Non-EQuilibrium AIr Radiation (NEQAIR) codes. At peak heating, the shock layer is substantially in equilibrium. Comprehensive spatial and spectral grid convergence studies have been implemented to quantify grid effects on the convective and radiative heating levels. Axisymmetric tests including the coupled effects of radiative energy transfer show negligible change to the convective heating and a 20 percent reduction in the radiative heating.
Lunar and Mars return conditions are examined using the LAURA flow field code and the LORAN radiation code to assess the effect of radiative coupling on axisymmetric thermochemical nonequilibrium flows. Coupling of the two codes is achieved iteratively. Special treatment required to couple radiation in a shock-capturing method is discussed. Results indicate that while coupling effects are generally the same as occur in equilibrium flows, under certain conditions radiation can modify the chemical kinetics of a nonequilibrium flow and thus alter relaxation processes. Coupling effects are found to be small for all cases considered, except for a five meter diameter aerobrake returning from Mars at 13.6 kilometers per second.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Lunar and Mars return conditions are examined using the LAURA flowfield code and the LORAN radiation code to assess the effect of radiative coupling on axisymmetric thermochemical nonequilibrium flows. Coupling of the two codes is achieved iteratively. Special treatment required to couple radiation in a shock-capturing method is discussed. Results indicate that while coupling effects are generally the same as occur in equilibrium flows, under certain conditions radiation can modify the chemical kinetics of a nonequilibrium flow and thus alter relaxation processes. Coupling effects are found to be small for all cases considered, except for a five meter diameter aerobrake returning from Mars at 13.6 km/sec.
A procedure which reduces the memory requirements for computing the viscous flow over a modified Orbiter geometry at a hypersonic flight condition is presented. The Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) code which incorporates a thermochemical nonequilibrium chemistry model, a finite rate catalytic wall boundary condition and wall temperature distribution based on radiation equilibrium is used in this study. In addition, the effect of choice of 'min mod' function, eigenvalue limiter and grid density on surface heating is investigated. The surface heating from a flowfield calculation at Mach number 22, altitude of 230,000 ft and 40 deg angle of attack is compared with flight data from three Orbiter flights.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The basic features of a Microwave Reflectometer Ionization Sensor (MRIS) as designed for utilization on the Aeroassist Flight Experiment are described. A parametric study of the effects of trajectory and unknowns in the thermochemical nonequilibrium models for translational and vibrational-electronic energy exchange rates, reaction rates, and the average electronic excitation level of atoms is performed to illustrate how the MRIS data may be employed for code validation. This study, implemented with program LAURA, shows a particular sensitivity of the onset and severity of an electron avalanche phenomena associated with changes in these physical models.
Explore the source record for details and available documents.
Simulations of the near wake flowfield behind three aerobrakes have been implemented with Program LAURA, an algorithm for obtaining the numerical solution to the governing equations for three-dimensional, viscous, hypersonic flows in chemical and thermal nonequilibrium. Emphasis is placed on understanding the conditons which are likely to cause the shear layer to impinge on a payload positioned behind the aerobrake. A linear relationship between shear layer deflection angle and angle of attack (or lift-to-drag ratio) has been identified in several ground based tests. A similar relation appears in the numerical simulations, though there is some evidence that deflection angels may increase somewhat due to the effects of gas chemistry. Shear layer impingement can raise local heating levels a factor of 10 higher than levels present without impingement. Payload heating levels near impingement points are a larger percentage of stagnation point heating levels at higher altitudes and peak payload heating levels are likely to occur earlier in the trajectory than peak heating on the forebody stagnation point.
A nonequilibrium radiative heating prediction method has been used to evaluate several energy exchange models used in nonequilibrium computational fluid dynamics methods. The radiative heating measurements from the FIRE II flight experiment supply an experimental benchmark against which different formulations for these exchange models can be judged. The models which predict the lowest radiative heating are found to give the best agreement with the flight data. Examination of the spectral distribution of radiation indicates that despite close agreement of the total radiation, many of the models examined predict excessive molecular radiation. It is suggested that a study of the nonequilibrium chemical kinetics may lead to a correction for this problem.
A point-implicit, asynchronous macrotasked relaxation of the steady, thin-layer, Navier-Stokes equations is presented. The method employs multidirectional, single-level storage Gauss-Seidel relaxation sweeps, which effectively communicate perturbations across the entire domain in 2n sweeps, where n is the dimension of the domain. In order to enhance convergence the application of relaxation factors to specific components of the Jacobian is examined using a stability analysis of the advection and diffusion equations. Attention is also given to the complications associated with asynchronous multitasking. Solutions are generated for hypersonic flows over blunt bodies in two and three dimensions with chemical reactions, utilizing single-tasked and multitasked relaxation strategies.
Program LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is an upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3D viscous hypersonic flows in chemical and thermal nonequilibrium. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with a modified Roe's averaging and with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. The code has been applied to Problem 8.2 of this workshop for the case of thermochemical nonequilibrium flow through a nozzle. Chemical reaction rates are defined with the model of Park (1987). Thermal nonequilibrium is modeled using a two-temperature approximation in which the vibrational energies of all molecules are assumed to be in equilibrium at a single temperature which is generally different from the translational-rotational temperature. Two grids were used to define the flow for the original problem, with a stagnation temperature of 6500 K. A third case with a stagnation temperature of 10,000 K is also presented. The solution domain includes the converging nozzle, subsonic flow domain in which the gas is substantially in thermochemical equilibrium and the diverging nozzle, hypersonic flow domain in which the gas is substantially in thermochemical nonequilibrium.