Search NASASearch

Engineering topics

Dinesh K Prabhu

Publications and source records attributed to Dinesh K Prabhu.

Assessment of the Fluid Dynamics Boundary Condition in Ablating or Blowing Flows

Improved models of ablative thermal protection systems have enabled the treatment of materials and fluid behavior in a coupled manner. This paper reports a new approach to modeling the interface between fluid and material, with attention to the conservation of species mass flux and energy on the fluid side of the interface. The general equation is presented and is shown to recover the traditional uncoupled fluid/materials response interface. Including the chemical reaction terms on the CFD side of the interface makes the heat flux exchange independent of the thermodynamic reference state and, therefore, a measurable quantity. Doing so allows the material response solver to take as input the surface heat flux rather than a film coefficient. Removing the film coefficient approximation enables more direct solution of vehicle thermal response but requires consistency in the wall state. The mixing of the shock layer and pyrolysis gas is then computed with finite rate chemistry within the fluid solver. The boundary conditions described have been implemented in the DPLR v4.05.1 code. Char removal is captured using finite rate chemistry in DPLR’s gas surface interaction module. Aspects of coupling these solutions to material response are discussed.

Ablation

Analysis of MSL/MEDLI Entry Data with Coupled CFD and Material Response

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield using NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performances for protecting probes during planetary entry. The Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP) [2]. The objective of this work is to showcase and analyze the coupling between the material response and the aerothermal environment. As shown in Figure 1, the workflow is divided into the following steps. First, the aerothermal properties are computed in the Data Parallel Line Relaxation (DPLR) code [3] and used with the Nonequilibrium air radiation (NEQAIR) program [8] to compute radiative heating. Second, the thermal response inside the material is computed in the Porous material Analysis Toolbox based on OpenFOAM (PATO) [4,5,6] using a fixed blowing correction parameter. Third, the pyrolysis gases computed in PATO are used as inputs to a blowing boundary condition within DPLR. Fourth, the new environment properties from DPLR are used in NEQAIR to provide an updated solution, then both the updated aerothermal environment and radiative heating are used in PATO without blowing correction. The third and fourth steps are then repeated until convergence in surface temperature is obtained. Convergence in the radiative heating is generally achieved before surface temperature, at which point the radiative heating is no longer updated. Char mass loss rates are forced to zero to produce a non-receding surface condition. For early time points in the trajectory, where flow around the MSL aeroshell is rarefied, the Direct Simulation Monte Carlo (DSMC) code, SPARTA [7], is used to compute the aerothermal environment. Iteration between PATO and SPARTA is not performed due to the computational cost of DSMC simulations. Preliminary results of the coupling between PATO and DPLR for the MSL heatshield atmospheric entry model are presented in Figures 2-4 at 65 seconds after entry interface. Figure 2 shows the surface temperature results from an uncoupled simulation in PATO with the blowing correction parameter applied (left) along with the coupled surface temperature after iteration (right). Figure 3 shows the surface temperature along the centerline from windward to leeward for easier comparison. Figure 4 shows the coupled and uncoupled pyrolysis gas blowing rate. Mars 2020 used a similar heatshield consisting of PICA for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the predictive material response capability is benchmarked against flight data from MEDLI. This work represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes.

Thermal Protection Systems

MEDLI2 Material Response Model Development and Validation

During the entry of Mars Science Laboratory (MSL) the heatshield was equipped with the instrumentation suite Mars Entry, Descent, and Landing Instruments (MEDLI). In-depth thermocouple (TC) data was used to reconstruct the surface heating and temperatures. Discrepancies between MEDLI’s recorded in-depth heating data and the predicted thermal response and recession led to a campaign to characterize and test Mars Entry, Descent, and Landing Instrument 2 (MEDLI2) flight lot thermal protection materials (TPS) at relevant temperatures and pressures for Martian entry conditions. This paper covers the development and validation of MEDLI2 flight lot material response models for the heatshield material Phenolic Impregnated Carbon Ablator (PICA) and the backshell material Super-Lightweight Ablator (SLA-561V). Virgin and char thermal conductivities were updated as a function of temperature and pressure in the MEDLI2 specific material response models. Other material properties such as virgin and char density, emissivity, absorptivity, and specific heat capacitance was also characterized and compared to the “Heritage” models used during the TPS design phase of Mars 2020. Fully Implicit Ablation and Thermal response (FIAT) simulations were completed using the Heritage and MEDLI2 material response models to provide evidence of the increased accuracy of the MEDLI2 model. FIAT predicted in-depth temperatures were compared to flight lot certification ground-test arc jet test data and MEDLI2 flight thermocouple data.

Material response

Simulations of Hayabusa2 Atmospheric Entry and Comparisons with Data from the Imaging Campaign

Results from flow and radiation simulations for the Hayabusa2 Sample Return Capsule (SRC) along its best-estimated trajectory (BET) are presented. Flow fields are first computed at several points along the best-estimated trajectory for a freestream consisting of N2, O2, Ar, CO2, and H2O; the last three are trace species in Earth’s atmosphere. The thermal response of an assumed ablative material is estimated and flow fields are recomputed at each time point with modified surface temperature distributions. The spectral radiance, from the combination of an aeroheated capsule and the high-temperature shock layer around it, is computed on an imaging plane for a given view angle and slant range at each time point. The radiance is transported through the intervening atmosphere between the Hayabusa2 SRC and the observing aircraft and upon integration over solid angle representing the field of view of the instruments converted to a spectral irradiance. The spectral irradiance is convolved with the instrument function to obtain results that can be directly compared with measurement. Fair to good comparisons with measurements required shifting the time coordinate of measurements by 7.5 s. The disagreement between prediction and measurement is substantial after about 72.5 s along the BET. Additional computations seem to support the idea of flow transition over the capsule.

Atmospheric entry

Simulations of Hayabusa2 Atmospheric Entry and Comparisons with Data from the Imaging Campaign

Results from flow and radiation simulations for the Hayabusa2 Sample Return Capsule (SRC) along its best-estimated trajectory (BET) are presented. Flow fields are first computed at several points along the best-estimated trajectory for a freestream consisting of N2, O2, Ar, CO2, and H2O; the last three are trace species in Earth’s atmosphere. The thermal response of an assumed ablative material is estimated and flow fields are recomputed at each time point with modified surface temperature distributions. The spectral radiance, from the combination of an aeroheated capsule and the high-temperature shock layer around it, is computed on an imaging plane for a given view angle and slant range at each time point. The radiance is transported through the intervening atmosphere between the Hayabusa2 SRC and the observing aircraft and upon integration over solid angle representing the field of view of the instruments converted to a spectral irradiance. The spectral irradiance is convolved with the instrument function to obtain results that can be directly compared with measurement. Fair to good comparisons with measurements required shifting the time coordinate of measurements by 7.5 s. The disagreement between prediction and measurement is substantial after about 72.5 s along the BET. Additional computations seem to support the idea of flow transition over the capsule.

Atmospheric entry

MARS 2020 Backshell Radiative Heating Measurement and Shock Tube Verification

In February 2021, the Perseverance rover was brought to the surface of Mars by the Mars 2020 mission. A feature of the Mars 2020 capsule was instrumentation to measure its entry, descent and landing (EDL) with the so-called Mars EDL Instrumentation 2 (MEDLI2) [1]. The MEDLI2 introduced, among other things, backshell instrumentation, including a broadband radiometer. The radiometer was mounted on the leeside of the vehicle next to thermocouple plugs and a heat flux gauge. The data returned by the leeside MEDLI2 heat flux gauge is largely analogous to that measured by the COMARS gauge flown on the ExoMars Schiaparelli [2] entry in that it is measuring heat flux in an area that is entirely dominated by radiative heating. The COMARS measurement provided excellent validation of backshell radiative heating models for Mars entry, albeit at a limited number of points [3]. The MEDLI2 heat flux gauge measurement effectively confirmed the quality of the prediction, extended over the full trajectory, as will be presented in this paper. The MEDLI2 radiometer,however, was blocked by ablation products and suffered a loss of half of its signal. The second backshell heat flux gauge installed on the windside of the vehicle was also well predicted, although the heating had both radiative and convective contributions. It was desired to reproduce the conditions of the Mars 2020 entry via ground testing in the Electric Arc Shock Tube (EAST) at NASA Ames. Tests to verify stagnation line heating were previously reported in EAST, confirming the presence of shock layer radiation as the major discrepancy in heatshield temperature modelling [4]. Therefore, tests for stagnation line heating were not repeated. Instead, the shock tube informed bias method [5, 6] was used to identify test conditions that may produce similarity to streamlines that pass around the backside of the vehicle and are responsible for the radiation observed at the two heat flux gauge locations. This method was used to identify a range of velocities and densities in the shock tube that are relevant for confirming the radiative environment encountered. This paper reports the results obtained in the 10 cm diameter EAST shock tube, corresponding to later trajectory points at ambient pressures of 1.1-2.0 Torr and velocities from 1.2-3.5 km/s. The test series employed two primary diagnostics: emission spectroscopy and tunable diode laser absorption spectroscopy (TDLAS). The emission spectroscopy performed broadband measurements of the radiative emission of the 4.3 and 2.7 m bands of CO2 at flight similar conditions, obtaining both spectral and spatial data corresponding to the relaxation behind the shock front. The TDLAS measured the absorption of several lines of CO and CO2 and obtains species number densities and temperatures as a function of time behind the shock front. This paper will review highlights of this test series and analyses of the emission and absorption data. While the datasets generally show good agreement with predictions, a few discrepancies and items for additional investigation are identified and will be discussed

Brett A Cruden

Coupling CFD and Material Response for Analysis of Mars Entry

In computing the response of an ablating thermal protection system during atmospheric entry, the aerothermal environment and material response are generally computed separately with a blowing correction term in the material response model to account for the blowing of char and pyrolysis gases. In this work, we apply a coupled approach in which pyrolysis blowing gases, computed in the PATO material response code, are used with a blowing boundary condition in the DPLR hypersonic CFD code. This leads to an iterative method in which blowing products from PATO are input into DPLR to update surface heating estimates. The full iterative method, with the addition of radiative heating estimates using the NEQAIR radiation solver, is shown in Fig. 1. The method is demonstrated on a sphere case with the environment and material properties based on the Mars Science Laboratory entry. Future work includes utilizing this method in computing full 3D material response during the Mars 2020 entry.

Ablation

Parametric Analysis of Entry Vehicles for Giant Planet Missions

Summarize concept designs for entry vehicles for giant planet destinations identified as high priority in the latest decadal survey. Identify viable entry states, aerothermal environments, and heatshield thermal protection systems (TPS) for various vehicle dimensions. Viable TPS options and sizing focus on low-g (< 50g) entry conditions to avoid costly instrumentation requalification. TPS manufacturing capabilities are also taken into consideration while adhering to the manufacturing capabilities of the TPS options. Margined mass estimates for heatshield TPS options are provided for Uranus and Saturn test cases.

Giant Planet