Benjamin S. Kirk

Development and Verification of the Charring Ablating Thermal Protection Implicit System Solver

The development and verification of the Charring Ablating Thermal Protection Implicit System Solver is presented. This work concentrates on the derivation and verification of the stationary grid terms in the equations that govern three-dimensional heat and mass transfer for charring thermal protection systems including pyrolysis gas flow through the porous char layer. The governing equations are discretized according to the Galerkin finite element method with first and second order implicit time integrators. The governing equations are fully coupled and are solved in parallel via Newtons method, while the fully implicit linear system is solved with the Generalized Minimal Residual method. Verification results from exact solutions and the Method of Manufactured Solutions are presented to show spatial and temporal orders of accuracy as well as nonlinear convergence rates.

Boundary Layer Transition Results From STS-114

The tool for predicting the onset of boundary layer transition from damage to and/or repair of the thermal protection system developed in support of Shuttle Return to Flight is compared to the STS-114 flight results. The Boundary Layer Transition (BLT) Tool is part of a suite of tools that analyze the aerothermodynamic environment of the local thermal protection system to allow informed disposition of damage for making recommendations to fly as is or to repair. Using mission specific trajectory information and details of each damage site or repair, the expected time of transition onset is predicted to help determine the proper aerothermodynamic environment to use in the subsequent thermal and stress analysis of the local structure. The boundary layer transition criteria utilized for the tool was developed from ground-based measurements to account for the effect of both protuberances and cavities and has been calibrated against flight data. Computed local boundary layer edge conditions provided the means to correlate the experimental results and then to extrapolate to flight. During STS-114, the BLT Tool was utilized and was part of the decision making process to perform an extravehicular activity to remove the large gap fillers. The role of the BLT Tool during this mission, along with the supporting information that was acquired for the on-orbit analysis, is reviewed. Once the large gap fillers were removed, all remaining damage sites were cleared for reentry as is. Post-flight analysis of the transition onset time revealed excellent agreement with BLT Tool predictions.

Review of Orbiter Flight Boundary Layer Transition Data

In support of the Shuttle Return to Flight program, a tool was developed to predict when boundary layer transition would occur on the lower surface of the orbiter during reentry due to the presence of protuberances and cavities in the thermal protection system. This predictive tool was developed based on extensive wind tunnel tests conducted after the loss of the Space Shuttle Columbia. Recognizing that wind tunnels cannot simulate the exact conditions an orbiter encounters as it re-enters the atmosphere, a preliminary attempt was made to use the documented flight related damage and the orbiter transition times, as deduced from flight instrumentation, to calibrate the predictive tool. After flight STS-114, the Boundary Layer Transition Team decided that a more in-depth analysis of the historical flight data was needed to better determine the root causes of the occasional early transition times of some of the past shuttle flights. In this paper we discuss our methodology for the analysis, the various sources of shuttle damage information, the analysis of the flight thermocouple data, and how the results compare to the Boundary Layer Transition prediction tool designed for Return to Flight.

Aeroheating Testing and Predictions for Project Orion CEV at Turbulent Conditions

An investigation of the aeroheating environment of the Project Orion Crew Exploration Vehicle was performed in the Arnold Engineering Development Center Hypervelocity Wind Tunnel No. 9 Mach 8 and Mach 10 nozzles and in the NASA Langley Research Center 20 - Inch Mach 6 Air Tunnel. Heating data were obtained using a thermocouple-instrumented approx.0.035-scale model (0.1778-m/7-inch diameter) of the flight vehicle. Runs were performed in the Tunnel 9 Mach 10 nozzle at free stream unit Reynolds numbers of 1x10(exp 6)/ft to 20x10(exp 6)/ft, in the Tunnel 9 Mach 8 nozzle at free stream unit Reynolds numbers of 8 x 10(exp 6)/ft to 48x10(exp 6)/ft, and in the 20-Inch Mach 6 Air Tunnel at free stream unit Reynolds numbers of 1x10(exp 6)/ft to 7x10(exp 6)/ft. In both facilities, enthalpy levels were low and the test gas (N2 in Tunnel 9 and air in the 20-Inch Mach 6) behaved as a perfect-gas. These test conditions produced laminar, transitional and turbulent data in the Tunnel 9 Mach 10 nozzle, transitional and turbulent data in the Tunnel 9 Mach 8 nozzle, and laminar and transitional data in the 20- Inch Mach 6 Air Tunnel. Laminar and turbulent predictions were generated for all wind tunnel test conditions and comparisons were performed with the experimental data to help define the accuracy of computational method. In general, it was found that both laminar data and predictions, and turbulent data and predictions, agreed to within less than the estimated 12% experimental uncertainty estimate. Laminar heating distributions from all three data sets were shown to correlate well and demonstrated Reynolds numbers independence when expressed in terms of the Stanton number based on adiabatic wall-recovery enthalpy. Transition onset locations on the leeside centerline were determined from the data and correlated in terms of boundary-layer parameters. Finally turbulent heating augmentation ratios were determined for several body-point locations and correlated in terms of the boundary-layer momentum Reynolds number.

Crew Exploration Vehicle (CEV) Crew Module shape selection analysis and CEV Aeroscience Project Overview

This paper details how NASA selected the shape of the CEV (Crew Exploration Vehicle) Crew Module (CM) and describes the approach used to develop the associated aerodynamic and aerothermodynamic databases. The shape study discussion provides information on the analysis performed to derive the CM shape. Many classes of vehicles were assessed, including capsules, slender bodies, lifting bodies, and winged vehicles. It was determined that a capsule shape similar to Apollo provided the best balance of risk, cost, and performance. After selecting the shape, the CEV Aerosciences Project (CAP) was formed to be responsible for producing the aerodynamic and aerothermodynamic databases for all phases of flight (on-orbit aero, RCS plume environments, nominal entry, ascent aborts, etc). The CAP team includes equal representation from NASA’s Lyndon B. Johnson Space Center (JSC), Ames Research Center (ARC), and Langley Research Center (LaRC). This team leverages NASA’s expertise in both analytical and experimental techniques in the fields of aerodynamics and aerothermodynamics.

Aerothermal Testing for Project Orion Crew Exploration Vehicle

The Project Orion Crew Exploration Vehicle aerothermodynamic experimentation strategy, as it relates to flight database development, is reviewed. Experimental data has been obtained to both validate the computational predictions utilized as part of the database and support the development of engineering models for issues not adequately addressed with computations. An outline is provided of the working groups formed to address the key deficiencies in data and knowledge for blunt reentry vehicles. The facilities utilized to address these deficiencies are reviewed, along with some of the important results obtained thus far. For smooth wall comparisons of computational convective heating predictions against experimental data from several facilities, confidence was gained with the use of algebraic turbulence model solutions as part of the database. For cavities and protuberances, experimental data is being used for screening various designs, plus providing support to the development of engineering models. With the reaction-control system testing, experimental data were acquired on the surface in combination with off-body flow visualization of the jet plumes and interactions. These results are being compared against predictions for improved understanding of aftbody thermal environments and uncertainties.

Overview of the CHarring Ablator Response (CHAR) Code

An overview of the capabilities of the CHarring Ablator Response (CHAR) code is presented. CHAR is a one-, two-, and three-dimensional unstructured continuous Galerkin finite-element heat conduction and ablation solver with both direct and inverse modes. Additionally, CHAR includes a coupled linear thermoelastic solver for determination of internal stresses induced from the temperature field and surface loading. Background on the development process, governing equations, material models, discretization techniques, and numerical methods is provided. Special focus is put on the available boundary conditions including thermochemical ablation and contact interfaces, and example simulations are included. Finally, a discussion of ongoing development efforts is presented.

BLIMPK/Streamline Surface Catalytic Heating Predictions on the Space Shuttle Orbiter

This paper describes the results of an analysis of localized catalytic heating effects to the U.S. Space Shuttle Orbiter Thermal Protection System (TPS). The analysis applies to the High-temperature Reusable Surface Insulation (HRSI) on the lower fuselage and wing acreage, as well as the critical Reinforced Carbon-Carbon on the nose cap, chin panel and the wing leading edge. The object of the analysis was to use a modified two-layer approach to predict the catalytic heating effects on the Orbiter windward HRSI tile acreage, nose cap, and wing leading edge assuming localized highly catalytic or fully catalytic surfaces. The method incorporated the Boundary Layer Integral Matrix Procedure Kinetic (BLIMPK) code with streamline inputs from viscous Navier-Stokes solutions to produce heating rates for localized fully catalytic and highly catalytic surfaces as well as for nominal partially catalytic surfaces (either Reinforced Carbon-Carbon or Reaction Cured Glass) with temperature-dependent recombination coefficients. The highly catalytic heating results showed very good correlation with Orbiter Experiments STS-2, -3, and -5 centerline and STS-5 wing flight data for the HRSI tiles. Recommended catalytic heating factors were generated for use in future Shuttle missions in the event of quick-time analysis of damaged or repaired TPS areas during atmospheric reentry. The catalytic factors are presented along the streamlines as well as a function of stagnation enthalpy so they can be used for arbitrary trajectories.

Coupled CFD-Ablation Response Model Simulations using the libMesh Framework

The DPLR Navier-Stokes flow solver is coupled to two Ablation Response Model (ARM) codes, CHAR and TITAN, using a modular approach where a central handler code runs the analysis codes iteratively and passes the required boundary condition values back and forth between the codes. The handler code is based on the libMesh software libraries. The libMesh mesh-free interpolation routines allow for the coupled analysis codes to be written in different programming languages and for the boundary point data to be non-point-matched. The boundary data is interpolated using a K-D Tree mesh-free interpolation approach. The basic execution flow for the coupled DPLR-ARM code is presented, and the coupled DPLR-ARM code is applied to arc jet test cases. The libMesh mesh-free interpolation successfully transferred the required boundary condition data between the fluid dynamic and material response codes, and the coupled DPLRARM surface recession rates matched experimental measurements as well as previous computations.

Overset Grid Applications in Hypersonic Flow Using the DPLR Flow Solver

This paper describes the addition of an overset grid capability to the DPLR flow solver for hypersonic flow in thermochemical nonequilibrium. Modifications to the preexisting flow solver were simplified through the use of DiRTlib, a “solver neutral” library of overset utilities. The new capability is demonstrated on a series of examples, including the Orion Crew Module and other reentry vehicles. For the overset grids used in these examples, the hole cutting and interpolation stencils were determined using SUGGAR, a generalized grid assembly code that can naturally accommodate both the three-dimensional and true two-dimensional cell-centered discretization schemes in DPLR. First a series of building-block examples are presented which highlight aspects of the new capability and assess the technique with comparisons to baseline, block-structured discretizations. The new capability is then exercised for the specific case of a tension tie geometry protruding from the Orion heatshield at both wind tunnel and flight conditions. The addition of overset capability to the DPLR flow solver is seen to be an essential feature for analyzing increasingly complex geometries in thermochemical nonequilibrium.