Search

A methodology to correlate flight aeroheating environments to the arc-jet environment is presented. For a desired hot-wall flight heating rate, the methodology provides the arcjet bulk enthalpy for the corresponding cold-wall heating rate. A series of analyses were conducted to examine the effects of the test sample model holder geometry to the overall performance of the test sample. The analyses were compared with arc-jet test samples and challenges and issues are presented. The transient flight environment was calculated for the Hypersonic Inflatable Aerodynamic Decelerator (HIAD) Earth Atmospheric Reentry Test (HEART) vehicle, which is a planned demonstration vehicle using a large inflatable, flexible thermal protection system to reenter the Earth's atmosphere from the International Space Station. A series of correlations were developed to define the relevant arc-jet test environment to properly approximate the HEART flight environment. The computed arcjet environments were compared with the measured arc-jet values to define the uncertainty of the correlated environment. The results show that for a given flight surface heat flux and a fully-catalytic TPS, the flight relevant arc-jet heat flux increases with the arc-jet bulk enthalpy while for a non-catalytic TPS the arc-jet heat flux decreases with the bulk enthalpy.

A model for simulating nonequilibrium radiation from Mars entry shock layers is presented. A new chemical kinetic rate model is developed that provides good agreement with recent EAST and X2 shock tube radiation measurements. This model includes a CO dissociation rate that is a factor of 13 larger than the rate used widely in previous models. Uncertainties in the proposed rates are assessed along with uncertainties in translational-vibrational relaxation modeling parameters. The stagnation point radiative flux uncertainty due to these flowfield modeling parameter uncertainties is computed to vary from 50 to 200% for a range of free-stream conditions, with densities ranging from 5e-5 to 5e-4 kg/m3 and velocities ranging from of 6.3 to 7.7 km/s. These conditions cover the range of anticipated peak radiative heating conditions for proposed hypersonic inflatable aerodynamic decelerators (HIADs). Modeling parameters for the radiative spectrum are compiled along with a non-Boltzmann rate model for the dominant radiating molecules, CO, CN, and C2. A method for treating non-local absorption in the non-Boltzmann model is developed, which is shown to result in up to a 50% increase in the radiative flux through absorption by the CO 4th Positive band. The sensitivity of the radiative flux to the radiation modeling parameters is presented and the uncertainty for each parameter is assessed. The stagnation point radiative flux uncertainty due to these radiation modeling parameter uncertainties is computed to vary from 18 to 167% for the considered range of free-stream conditions. The total radiative flux uncertainty is computed as the root sum square of the flowfield and radiation parametric uncertainties, which results in total uncertainties ranging from 50 to 260%. The main contributors to these significant uncertainties are the CO dissociation rate and the CO heavy-particle excitation rates. Applying the baseline flowfield and radiation models developed in this work, the radiative heating for the Mars Pathfinder probe is predicted to be nearly 20 W/cm2. In contrast to previous studies, this value is shown to be significant relative to the convective heating.

This paper investigates the shock-layer radiative heating uncertainty for hyperbolic Earth entry, with the main focus being a Mars return. In Part I of this work, a baseline simulation approach involving the LAURA Navier-Stokes code with coupled ablation and radiation is presented, with the HARA radiation code being used for the radiation predictions. Flight cases representative of peak-heating Mars or asteroid return are de ned and the strong influence of coupled ablation and radiation on their aerothermodynamic environments are shown. Structural uncertainties inherent in the baseline simulations are identified, with turbulence modeling, precursor absorption, grid convergence, and radiation transport uncertainties combining for a +34% and ..24% structural uncertainty on the radiative heating. A parametric uncertainty analysis, which assumes interval uncertainties, is presented. This analysis accounts for uncertainties in the radiation models as well as heat of formation uncertainties in the flow field model. Discussions and references are provided to support the uncertainty range chosen for each parameter. A parametric uncertainty of +47.3% and -28.3% is computed for the stagnation-point radiative heating for the 15 km/s Mars-return case. A breakdown of the largest individual uncertainty contributors is presented, which includes C3 Swings cross-section, photoionization edge shift, and Opacity Project atomic lines. Combining the structural and parametric uncertainty components results in a total uncertainty of +81.3% and ..52.3% for the Mars-return case. In Part II, the computational technique and uncertainty analysis presented in Part I are applied to 1960s era shock-tube and constricted-arc experimental cases. It is shown that experiments contain shock layer temperatures and radiative ux values relevant to the Mars-return cases of present interest. Comparisons between the predictions and measurements, accounting for the uncertainty in both, are made for a range of experiments. A measure of comparison quality is de ned, which consists of the percent overlap of the predicted uncertainty bar with the corresponding measurement uncertainty bar. For nearly all cases, this percent overlap is greater than zero, and for most of the higher temperature cases (T >13,000 K) it is greater than 50%. These favorable comparisons provide evidence that the baseline computational technique and uncertainty analysis presented in Part I are adequate for Mars-return simulations. In Part III, the computational technique and uncertainty analysis presented in Part I are applied to EAST shock-tube cases. These experimental cases contain wavelength dependent intensity measurements in a wavelength range that covers 60% of the radiative intensity for the 11 km/s, 5 m radius flight case studied in Part I. Comparisons between the predictions and EAST measurements are made for a range of experiments. The uncertainty analysis presented in Part I is applied to each prediction, and comparisons are made using the metrics defined in Part II. The agreement between predictions and measurements is excellent for velocities greater than 10.5 km/s. Both the wavelength dependent and wavelength integrated intensities agree within 30% for nearly all cases considered. This agreement provides confidence in the computational technique and uncertainty analysis presented in Part I, and provides further evidence that this approach is adequate for Mars-return simulations. Part IV of this paper reviews existing experimental data that include the influence of massive ablation on radiative heating. It is concluded that this existing data is not sufficient for the present uncertainty analysis. Experiments to capture the influence of massive ablation on radiation are suggested as future work, along with further studies of the radiative precursor and improvements in the radiation properties of ablation products.

Starting in January 2004, the NESC has received several communications from knowledgeable technical experts at NASA expressing shared concerns (mainly at the Langley Research Center (LaRC) and Ames Research Center (ARC)) about Huygens mission success. It was suggested that NASA become more technically involved directly in the analysis of Huygens' entry, descent and landing (EDL) focusing on the parachute deployment trigger performance and the resultant effects on the operation of the parachute system, and the determination of the radiative heating environment at Titan by ESA and the corresponding thermal protection system (TPS) response. A NESC Team was formed and tasked to provide an independent assessment of these concerns. The results of that assessment are documented in this report.

Using the coupled ablation and radiation capability recently included in the LAURA flowfield solver, this paper investigates the influence of ablation on the shock-layer radiative heating for Earth entry. The extension of the HARA radiation model, which provides the radiation predictions in LAURA, to treat a gas consisting of the elements C, H, O, and N is discussed. It is shown that the absorption coefficient of air is increased with the introduction of the C and H elements. A simplified shock layer model is studied to show the impact of temperature, as well as the abundance of C and H, on the net absorption or emission from an ablation contaminated boundary layer. It is found that the ablation species reduce the radiative flux in the vacuum ultraviolet, through increased absorption, for all temperatures. However, in the infrared region of the spectrum, the ablation species increase the radiative flux, through strong emission, for temperatures above 3,000 K. Thus, depending on the temperature and abundance of ablation species, the contaminated boundary layer may either provide a net increase or decrease in the radiative flux reaching the wall. To assess the validity of the coupled ablation and radiation LAURA analysis, a previously analyzed Mars-return case (15.24 km/s), which contains significant ablation and radiation coupling, is studied. Exceptional agreement with previous viscous shock-layer results is obtained. A 40% decrease in the radiative flux is predicted for ablation rates equal to 20% of the free-stream mass flux. The Apollo 4 peak-heating case (10.24 km/s) is also studied. For ablation rates up to 3.4% of the free-stream mass flux, the radiative heating is reduced by up to 19%, while the convective heating is reduced by up to 87%. Good agreement with the Apollo 4 radiometer data is obtained by considering absorption in the radiometer cavity. For both the Mars return and the Apollo 4 cases, coupled radiation alone is found to reduce the radiative heating by 30 60% and the convective heating by less than 5%.

This paper investigates the non-Boltzmann modeling of the radiating atomic and molecular electronic states present in lunar-return shock-layers. The Master Equation is derived for a general atom or molecule while accounting for a variety of excitation and de-excitation mechanisms. A new set of electronic-impact excitation rates is compiled for N, O, and N2+, which are the main radiating species for most lunar-return shock-layers. Based on these new rates, a novel approach of curve-fitting the non-Boltzmann populations of the radiating atomic and molecular states is developed. This new approach provides a simple and accurate method for calculating the atomic and molecular non-Boltzmann populations while avoiding the matrix inversion procedure required for the detailed solution of the Master Equation. The radiative flux values predicted by the present detailed non-Boltzmann model and the approximate curve-fitting approach are shown to agree within 5% for the Fire 1634 s case.

The radiative heating environment for the Huygens probe near peak heating conditions for Titan entry is investigated in this paper. The task of calculating the radiation-coupled flowfield, accounting for non-Boltzmann and non-optically thin radiation, is simplified to a rapid yet accurate calculation. This is achieved by using the viscous-shock layer (VSL) technique for the stagnation-line flowfield calculation and a modified smeared rotational band (SRB) model for the radiation calculation. These two methods provide a computationally efficient alternative to a Navier-Stokes flowfield and line-by-line radiation calculation. The results of the VSL technique are shown to provide an excellent comparison with the Navier-Stokes results of previous studies. It is shown that a conventional SRB approach is inadequate for the partially optically-thick conditions present in the Huygens shock-layer around the peak heating trajectory points. A simple modification is proposed to the SRB model that improves its accuracy in these partially optically-thick conditions. This modified approach, labeled herein as SRBC, is compared throughout this study with a detailed line-by-line (LBL) calculation and is shown to compare within 5% in all cases. The SRBC method requires many orders-of-magnitude less computational time than the LBL method, which makes it ideal for coupling to the flowfield. The application of a collisional-radiative (CR) model for determining the population of the CN electronic states, which govern the radiation for Huygens entry, is discussed and applied. The non-local absorption term in the CR model is formulated in terms of an escape factor, which is then curve-fit with temperature. Although the curve-fit is an approximation, it is shown to compare well with the exact escape factor calculation, which requires a computationally intensive iteration procedure.

This paper presents a detailed analysis of the shock-layer radiative heating to the Fire II vehicle using a new air radiation model and a viscous shock-layer flowfield model. This new air radiation model contains the most up-to-date properties for modeling the atomic-line, atomic photoionization, molecular band, and non-Boltzmann processes. The applied viscous shock-layer flowfield analysis contains the same thermophysical properties and nonequilibrium models as the LAURA Navier-Stokes code. Radiation-flowfield coupling, or radiation cooling, is accounted for in detail in this study. It is shown to reduce the radiative heating by about 30% for the peak radiative heating points, while reducing the convective heating only slightly. A detailed review of past Fire II radiative heating studies is presented. It is observed that the scatter in the radiation predicted by these past studies is mostly a result of the different flowfield chemistry models and the treatment of the electronic state populations. The present predictions provide, on average throughout the trajectory, a better comparison with Fire II flight data than any previous study. The magnitude of the vacuum ultraviolet (VUV) contribution to the radiative flux is estimated from the calorimeter measurements. This is achieved using the radiometer measurements and the predicted convective heating. The VUV radiation predicted by the present model agrees well with the VUV contribution inferred from the Fire II calorimeter measurement, although only when radiation-flowfield coupling is accounted for. This agreement provides evidence that the present model accurately models the VUV radiation, which is shown to contribute significantly to the Fire II radiative heating.

An investigation of the aerothermodynamic environment of the Huygens entry probe has been conducted. A Monte Carlo simulation of the trajectory of the probe during entry into Titan's atmosphere was performed to identify a worst-case heating rate trajectory. Flowfield and radiation transport computations were performed at points along this trajectory to obtain convective and radiative heat-transfer distributions on the probe's heat shield. This investigation identified important physical and numerical factors, including atmospheric CH4 concentration, transition to turbulence, numerical diffusion modeling, and radiation modeling, which strongly influenced the aerothermodynamic environment.

Convective and radiative heating environments have been computed for a three-dimensional ellipsled configuration which would perform an aerocapture maneuver at Neptune. This work was performed as part of a one-year Neptune aerocapture spacecraft systems study that also included analyses of trajectories, atmospheric modeling, aerodynamics, structural design, and other disciplines. Complementary heating analyses were conducted by separate teams using independent sets of aerothermodynamic modeling tools (i.e. Navier-Stokes and radiation transport codes). Environments were generated for a large 5.50 m length ellipsled and a small 2.88 m length ellipsled. Radiative heating was found to contribute up to 80% of the total heating rate at the ellipsled nose depending on the trajectory point. Good agreement between convective heating predictions from the two Navier-Stokes solvers was obtained. However, the radiation analysis revealed several uncertainties in the computational models employed in both sets of codes, as well as large differences between the predicted radiative heating rates.

This report documents the results of a computational study done to compute the inviscid longitudinal aerodynamic characteristics of the Space Shuttle Orbiter for Mach numbers 10 and 15 at angles of attack of 40, 50, 55, and 60 degrees. These computations were done to provide limited aerodynamic data in support of the Orbiter contingency abort task. The Orbiter had all the control surfaces in the undeflected position. The unstructured grid software FELISA was used for these computations with the equilibrium air option. Normal and axial force coefficients and pitching moment coefficients were computed. The hinge moment coefficients of the body flap and the inboard and outboard elevons were also computed. These results were compared with Orbiter Air Data Book (OADB) data and those computed using GASP. The comparison with the GASP results showed very good agreement in Cm and Ca at all the points. The computed axial force coefficients were smaller than those computed by GASP. There were noticeable differences between the present results and those in the OADB at angles of attack greater than 50 degrees.

This report documents the results of a study conducted to compute the inviscid longitudinal aerodynamic characteristics of three aeroshell configurations of the proposed '07 Mars lander. This was done in support of the activity to design a smart lander for the proposed '07 Mars mission. In addition to the three configurations with tabs designated as the shelf, the canted, and the Ames, the baseline configuration (without tab) was also studied. The unstructured grid inviscid CFD software FELISA was used, and the longitudinal aerodynamic characteristics of the four configurations were computed for Mach number of 2.3, 2.7, 3.5, and 4.5, and for an angle of attack range of -4 to 20 degrees. Wind tunnel tests had been conducted on scale models of these four configurations in the Unitary Plan Wind Tunnel, NASA Langley Research Center. Present computational results are compared with the data from these tests. Some differences are noticed between the two results, particularly at the lower Mach numbers. These differences are attributed to the pressures acting on the aft body. Most of the present computations were done on the forebody only. Additional computations were done on the full body (forebody and afterbody) for the baseline and the Shelf configurations. Results of some computations done (to simulate flight conditions) with the Mars gas option and with an effective gamma are also included.

This report is a review of a team effort that focuses on advanced aerospace concepts of the 21st Century. The paper emphasis advanced technologies, rather than cataloging every unusual aircraft that has ever been attempted. To dispel the myth that "aerodynamics is a mature science" an extensive list of "What we cannot do, or do not know" was enumerated. A zeit geist, a feeling for the spirit of the times, was developed, based on existing research goals. Technological drivers and the constraints that might influence these technological developments in a future society were also examined. The present status of aeronautics, space exploration, and non-aerospace applications, both military and commercial, including enabling technologies are discussed. A discussion of non-technological issues affecting advanced concepts research is presented. The benefit of using the study of advanced vehicles as a tool to uncover new directions for technology development is often necessary. An appendix is provided containing examples of advanced vehicle configurations currently of interest.

The accuracy and complexity of solving multicomponent gaseous diffusion using the detailed multicomponent equations, the Stefan-Maxwell equations, and two commonly used approximate equations have been examined in a two part study. Part I examined the equations in a basic study with specified inputs in which the results are applicable for many applications. Part II addressed the application of the equations in the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) computational code for high-speed entries in Earth's atmosphere. The results showed that the presented iterative scheme for solving the Stefan-Maxwell equations is an accurate and effective method as compared with solutions of the detailed equations. In general, good accuracy with the approximate equations cannot be guaranteed for a species or all species in a multi-component mixture. 'Corrected' forms of the approximate equations that ensured the diffusion mass fluxes sum to zero, as required, were more accurate than the uncorrected forms. Good accuracy, as compared with the Stefan- Maxwell results, were obtained with the 'corrected' approximate equations in defining the heating rates for the three Earth entries considered in Part II.

Results are obtained for the forebody of a planetary exploration vehicle entering the earth's atmosphere. A viscous-shock-layer analysis is used assuming the flow to be laminar and in chemical equilibrium. Presented results include coupled radiation and ablation injection. This study further includes the effect of different transport and thermodynamic properties and radiation models. A Lewis number of 1.4 appears adequate for the radiation-dominated flows. Five velocities corresponding to different possible trajectory points at an altitude of 70 km have been further analyzed in detail. Sublimation and radiative equilibrium wall temperatures are employed for cases with and without coupled injection, respectively. For the cases analyzed here, the mass injection rates are small. However, the rates could become large if a lower altitude is used for aerobraking and/or the body size is increased. A comparison of the equilibrium results with finite-rate chemistry calculation shows the flowfield to be in chemical equilibrium.