Search NASA⌕ Search

Engineering topics

Haas, Brian L.

Publications and source records attributed to Haas, Brian L..

Particle kinetic simulation of high altitude hypervelocity flight

Rarefied flows about hypersonic vehicles entering the upper atmosphere or through nozzles expanding into a near vacuum may only be simulated accurately with a direct simulation Monte Carlo (DSMC) method. Under this grant, researchers enhanced the models employed in the DSMC method and performed simulations in support of existing NASA projects or missions. DSMC models were developed and validated for simulating rotational, vibrational, and chemical relaxation in high-temperature flows, including effects of quantized anharmonic oscillators and temperature-dependent relaxation rates. State-of-the-art advancements were made in simulating coupled vibration-dissociation recombination for post-shock flows. Models were also developed to compute vehicle surface temperatures directly in the code rather than requiring isothermal estimates. These codes were instrumental in simulating aerobraking of NASA's Magellan spacecraft during orbital maneuvers to assess heat transfer and aerodynamic properties of the delicate satellite. NASA also depended upon simulations of entry of the Galileo probe into the atmosphere of Jupiter to provide drag and flow field information essential for accurate interpretation of an onboard experiment. Finally, the codes have been used extensively to simulate expanding nozzle flows in low-power thrusters in support of propulsion activities at NASA-Lewis. Detailed comparisons between continuum calculations and DSMC results helped to quantify the limitations of continuum CFD codes in rarefied applications.

Boyd, Iain↗

Particle kinetic simulation of high altitude hypervelocity flight

Rarefied flows about hypersonic vehicles entering the upper atmosphere or through nozzles expanding into a near vacuum may only be simulated accurately with a direct simulation Monte Carlo (DSMC) method. Under this grant, researchers enhanced the models employed in the DSMC method and performed simulations in support of existing NASA projects or missions. DSMC models were developed and validated for simulating rotational, vibrational, and chemical relaxation in high-temperature flows, including effects of quantized anharmonic oscillators and temperature-dependent relaxation rates. State-of-the-art advancements were made in simulating coupled vibration - dissociation - recombination for post-shock flows. Models were also developed to compute vehicle surface temperatures directly in the code rather than requiring isothermal estimates. These codes were instrumental in simulating aerobraking of NASA's Magellan spacecraft during orbital maneuvers to assess heat transfer and aerodynamic properties of the delicate satellite. NASA also depended upon simulations of entry of the Galileo probe into the atmosphere of Jupiter to provide drag and flow field information essential for accurate interpretation of an onboard experiment. Finally, the codes have been used extensively to simulate expanding nozzle flows in low-power thrusters in support of propulsion activities at NASA-Lewis. Detailed comparisons between continuum calculations and DSMC results helped to quantify the limitations of continuum CFD codes in rarefied applications.

Boyd, Iain↗

Aerothermodynamics of Pyrolizing Surfaces in Hypersonic Rarefied Flows

Direct simulation Monte Carlo (DSMC) calculations of rarefied flows about entry bodies typically employ a fixed surface temperature or a radiative-equilibrium energy balance to compute that temperature. Such boundary conditions neglect any effects of heat capacitance and heat conduction in the spacecraft heat shield and, therefore, provide an upper bound for the surface temperature. Such calculations also neglect pyrolysis from the heat shield which can be significant for a high-energy incident flow at very low densities. Accurate prediction of both heating and aerodynamic forces requires including pyrolysis and surface heat transfer in the models for gas-surface interaction employed in DSMC methods. Although these physical models have long appeared in various continuum flow calculation codes, they have only recently appeared in DSMC codes which are required to simulate rarefied flows during entry at high altitudes. In the current implementation, routines from the widely distributed Charring Material Thermal Response and Ablation (CMA) program are coupled into a DSMC code to calculate the one-dimensional heat transfer into the carbon phenolic heat shield at each point on a vehicle surface. Temperature-dependent material properties, surface re-radiation, and in-depth pyrolysis were included in the calculation, but surface ablation was neglected. Sample calculations for entry of the Galileo probe into the atmosphere of Jupiter demonstrate that including pyrolysis in the model leads to significant differences in predicted aerodynamics. Granted, the drag coefficient does not depend strongly on the surface temperature which can itself be significantly below the radiative equilibrium value during entry. However, the surface mass flux due to pyrolysis of the material is significant once the probe drops to altitudes characterized by transition flow. This leads to a noticeable increase in drag and a decrease in heating compared to a body without pyrolysis.

Haas, Brian L.↗

Rotational Relaxation Analysis of Nitrogen in Low-Density Freejet Expansions

The phenomena of rotational relaxation of nitrogen has been examined by numerous investigators over many years. One of the experiments which has been performed examines nonequilibrium flow in low-density free jet expansions. Data have been taken in such flows using a variety of techniques, including time-of-flight methods and electron beam fluorescence spectroscopy. The direct flow properties measured in these different investigations, such as density and translational, rotational and vibrational temperatures generally show reasonable agreement. However, this kind of correlation from experiment to experiment tends to be lost when these data are analyzed to obtain rotational relaxation time or collision number. The goal of such data analyses is to generate a succinct model for rotational relaxation in nitrogen which is essential for the computation of nonequilibrium rarefied flows. The objective of the present work is to process a large body of experimental data in a consistent manner to yield relaxation model parameters of the greatest utility for flow computations.

Gochberg, Lawrence A.↗

Particle kinetic simulation of high altitude hypervelocity flight

In this grant period, the focus has been on enhancement and application of the direct simulation Monte Carlo (DSMC) particle method for computing hypersonic flows of re-entry vehicles. Enhancement efforts dealt with modeling gas-gas interactions for thermal non-equilibrium relaxation processes and gas-surface interactions for prediction of vehicle surface temperatures. Both are important for application to problems of engineering interest. The code was employed in a parametric study to improve future applications, and in simulations of aeropass maneuvers in support of the Magellan mission. Detailed comparisons between continuum models for internal energy relaxation and DSMC models reveals that several discrepancies exist. These include definitions of relaxation parameters and the methodologies for implementing them in DSMC codes. These issues were clarified and all differences were rectified in a paper (Appendix A) submitted to Physics of Fluids A, featuring several key figures in the DSMC community as co-authors and B. Haas as first author. This material will be presented at the Fluid Dynamics meeting of the American Physical Society on November 21, 1993. The aerodynamics of space vehicles in highly rarefied flows are very sensitive to the vehicle surface temperatures. Rather than require prescribed temperature estimates for spacecraft as is typically done in DSMC methods, a new technique was developed which couples the dynamic surface heat transfer characteristics into the DSMC flow simulation code to compute surface temperatures directly. This model, when applied to thin planar bodies such as solar panels, was described in AIAA Paper No. 93-2765 (Appendix B) and was presented at the Thermophysics Conference in July 1993. The paper has been submitted to the Journal of Thermophysics and Heat Transfer. Application of the DSMC method to problems of practical interest requires a trade off between solution accuracy and computational expense and limitations. A parametric study was performed and reported in AIAA Paper No. 93-2806 (Appendix C) which assessed the accuracy penalties associated with simulations of varying grid resolution and flow domain size. The paper was also presented at the Thermophysics Conference and will be submitted to the journal shortly. Finally, the DSMC code was employed to assess the pitch, yaw, and roll aerodynamics of the Magellan spacecraft during entry into the Venus atmosphere at off-design attitudes. This work was in support of the Magellan aerobraking maneuver of May 25-Aug. 3, 1993. Furthermore, analysis of the roll characteristics of the configuration with canted solar panels was performed in support of the proposed 'Windmill' experiment. Results were reported in AIAA Paper No. 93-3676 (Appendix D) presented at the Atmospheric Flight Mechanics Conference in August 1993, and were submitted to Journal of Spacecraft and Rockets.

Haas, Brian L.↗

Models of thermal relaxation mechanics for particle simulation methods

A comparative study is conducted of the simulation of molecular vibrational relaxation in rarefied flows by means of existing models and a novel phenomenological microscopic model. This improved technique iterates between translation-rotation and rotation-vibration exchanges; equilibrium is in this way promoted if the model for the latter process is compatible with quantized oscillators. The iteration-equipartition model retains computational simplicity and promotes thermal equilibrium, even in the case of nondegenerate anharmonic quantized oscillators and multispecies gas mixtures.

Haas, Brian L.↗

Models for dynamic surface temperatures during rarefied aeropass maneuvers

New models for calculating dynamic surface temperatures in a direct simulation Monte Carlo particle method during atmospheric entry are described. The models take into account the radiative emissivity, thermal conductivity, and heat capacity of the material. They are applied to thin panels in a highly rarefied flows and yield the transient thermal response of individual surface facets. Simulation results for radiation panels with internal conductivity are verified by comparing them with free molecule theory.

Haas, Brian L.↗

Flow resolution and domain of influence in rarefied hypersonic blunt-body flows

The study assesses the effects of the upstream domain size and grid resolution upon flow properties and body aerodynamics computed for rarefied flows over cold blunt bodies with a direct simulation Monte Carlo (DSMC) particle method. Empirical correlations are suggested for aerodynamic coefficients for two-dimensional flows past a perpendicular flat plate. Free-stream parameters which were varied in the study include the Mach number, Knudsen number, surface temperature, and intermolecular potential. Insufficient grid resolution leads to overprediction of aerodynamic heating and forces in the DSMC method. Solution accuracy correlates well with the Reynolds number defined at the wall temperature and the stagnation mean free path relative to the cell dimension. Insufficient upstream domain size in the DSMC method leads to overprediction of heating and drag. Errors in aerodynamic coefficients correlate well with the distance ahead of the body where flow temperature reaches half of its peak value. Simulation of a hard-sphere gas is more sensitive to grid resolution, while simulation of a Maxwell gas is more sensitive to upstream domain size.

Haas, Brian L.↗

Models for direct Monte Carlo simulation of coupled vibration-dissociation

A new model for reactive collisions is developed within the framework of a particle method, which simulates coupled vibration-dissociation (CVD) behavior in high-temperature gases. The fundamental principles of particle simulation methods are introduced with particular attention given to the probability functions employed to select thermal and reactive collisions. Reaction probability functions are derived which favor vibrationally excited molecules as reaction candidates. The new models derived here are used to simulate CVD behavior during thermochemical relaxation of constant-volume O2 reservoirs, as well as the dissociation incubation behavior of postshock N2 flows for comparisons with previous models and experimental data.

Haas, Brian L.↗

Particle kinetic simulation of high altitude hypervelocity flight

In this grant period, the focus has been on the effects of thermo-chemical nonequilibrium in low-density gases, and on interactions between such gases and solid surfaces. Such conditions apply to hypersonic flows of re-entry vehicles, and to the expansion plumes of small rockets. Due to the nonequilibrium nature of these flows, a particle approach has been adopted. The method continues to undergo refinement and application to typical flows of interest. A number of studies have been performed for flows in thermo-chemical nonequilibrium. The effects of vibrational nonequilibrium on the rate of dissociation were studied for diatomic nitrogen. It was found that a new model reproduced the nonequilibrium behavior observed experimentally.

Heinemann, Klaus↗

Simulated rarefied aerodynamics of the Magellan spacecraft during aerobraking

Aerodynamic loads upon the Magellan spacecraft during aerobraking through the atmosphere of Venus are computed at off-design attitudes with a direct simulation Monte Carlo (DSMC) particle method. Simulated rarefied flows at nominal altitudes near 140 km and an entry speed of 8.6 km/s were compared to simulated and analytic free molecular results. Aerodynamic moments, forces, and heating for rarefied entry at all attitudes were 7-10 percent below free molecular results. All moments acted to restore the vehicle to its nominal zero-pitch, zero-yaw attitude. Suggested canting of the solar panels is an innovative configuration to assess gas-surface interaction during aerobraking. The resulting roll torques about the central body-axis as predicted in rarefied flow simulations were nearly twice that predicted for free molecular flow, although differences became less distinct for thermal accommodation coefficients well below unity. Roll torques increased dramatically with reduced accommodation coefficients employed in the simulation. In the DSMC code, periodic free-molecule boundary conditions and a coarse computational grid and body resolution served to minimize the simulation size and cost while retaining solution validity.

Haas, Brian L.↗

Particle simulation of rarefied aeropass maneuvers of the Magellan spacecraft

Rarefied and free molecular flowfields corresponding to proposed aeropass maneuvers through the atmosphere of Venus by the Magellan spacecraft are computed with a vectorized particle simulation method. The significance of rarefaction and reaction effects are assessed and surface heat flux, drag coefficients, and flowfield properties are computed. A simple surface heat transfer model is coupled directly to the simulation and assumes that each surface element is in radiative equilibrium with deep space. This allows direct computation of surface temperature distributions rather than requiring prescribed isothermal boundary conditions. Uncoupled dual-node heat transfer models, which account for heat capacity and thermal conductivity within each structural component, permit more accurate determination of surface temperatures. Such temperatures restrict the allowable altitudes and entry speeds of the aeropass maneuvers, and must therefore be estimated accurately. Simulations with an entry velocity of 8600 m/s at altitudes between 125 and 140 km reveal that allowable surface temperatures occur only at the highest altitudes where the flowfield is appropriately modeled as free-molecular. Excessive surface heat flux occurs at lower altitudes where molecular collisions are significant.

Haas, Brian L.↗

Resolution of differences between collision number definitions in particle and continuum simulations

The use of the same symbol, Z, representing a 'collision number' for thermal relaxation, has led to confusion regarding its definition in the context of both continuum and particle simulations. Examination of the relaxation mechanics employed in particle simulations demonstrates that these definitions differ by a numerical factor that depends upon the intermolecular potential. Particle and continuum simulations employing appropriate definitions of Z lead to identical results during isothermal and adiabatic stationary relaxation.

Lumpkin, Forrest E., III↗

Validation of chemistry models employed in a particle simulation method

The chemistry models employed in a statistical particle simulation method, as implemented in the Intel iPSC/860 multiprocessor computer, are validated and applied. Chemical relaxation of five-species air in these reservoirs involves 34 simultaneous dissociation, recombination, and atomic-exchange reactions. The reaction rates employed in the analytic solutions are obtained from Arrhenius experimental correlations as functions of temperature for adiabatic gas reservoirs in thermal equilibrium. Favorable agreement with the analytic solutions validates the simulation when applied to relaxation of O2 toward equilibrium in reservoirs dominated by dissociation and recombination, respectively, and when applied to relaxation of air in the temperature range 5000 to 30,000 K. A flow of O2 over a circular cylinder at high Mach number is simulated to demonstrate application of the method to multidimensional reactive flows.

Haas, Brian L.↗

Fundamentals of chemistry modeling applicable to a vectorized particle simulation

This paper describes the fundamentals of extending the vectorized particle simulation method derived by Baganoff and McDonald (1990), McDonald and Baganoff (1988), and McDonald (1989) for modeling chemically reacting flows. Details of reaction mechanics per reaction are presented, with particular attention given to the quantum nature of the vibrational mode. The models of reactive flows developed here were verified through a simulation of a superheated diatomic gas relaxing thermochemically to equilibrium in a reservoir.

Haas, Brian L.↗