An analysis of a charring ablation thermal protection system
Computer program for predicting charring ablation performance of heat shielding
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Computer program for predicting charring ablation performance of heat shielding
Char formation and surface recession of composite ablation material with silicone resin base
Ablation chars, carbon and graphite spectral emittance and reflectance as function of wavelength and temperature, noting applications to atmospheric entry heat shielding
Thermal conductivity determination for charring materials by dynamic technique, with application to carbon, graphite and silica cloth phenolics
Gas phase chemical reactions effect on heat transfer to charring ablator, deriving numerical solution for multicomponent stagnation point flow
Analysis of charring ablation with description of associated computing program
Ablative heat shield char layer, examining reacting nonequilibrium pyrolysis gas flow
Nonequilibrium modeling of pyrolysis gas products flow through char layer of ablative heat shield
Comparisons of ablation calculations with the charring ablation computer code and ablation test data are presented over a wide range of environmental conditions in air for three materials: low-density nylon phenolic, Avcoat 5026-39HC/G, and a filled silicon elastomer. Heat fluxes considered range from over 500 Btu/sq ft-sec to less than 50 Btu/sq ft-sec. Pressures range from 0.5 atm to .004 atm. Enthalpies range from about 2000 Btu/lb to 18000 Btu/lb. Predictions of recession, pyrolysis penetration, and thermocouple responses are considered. Recession predictions for nylon phenolic are good as steady state is approached, but strongly transient cases are underpredicted. Pyrolysis penetrations and thermocouple responses are very well predicted. Recession amounts for Avcoat and silicone elastomer are less well predicted, although high heat flux cases near steady state are fairly satisfactory. Pyrolysis penetrations and thermocouple responses are very well predicted.
The Hertz-Knudsen analysis is shown to accurately predict the sublimation rate from a charring ablator. Porosity is shown to have a significant effect on the surface temperature. The predominant carbon species found in the vapor is C3, which agrees well with the results of previous investigations.
General analyses implemented in FORTRAN IV for studying the energy transfer in the combined decomposition/char zone of a charring ablator are described.
Boron powder stabilizes char, so burned substances are shiny, smooth, and free of loose graphite fibers. Resin weight loss of laminates during burning in air is identical for the first three minutes for unfilled and boron-filled samples, then boron samples stabilize.
Class of imide-amine curing agents preserves structural integrity, prevents fiber release, and is fully compatible with conventional epoxy resins; agents do not detract from composite properties while greatly reducing char yield. Materials utilizing curing are used in aerospace, automotive, and other structural components where deterioration must be minimized and fiber release avoided in event of fire.
Analyzing the remains of Space Shuttle Columbia has proven technically beneficial years after the vehicle breakup. This investigation focused on charred deposits on fragments of Columbia overhead windowpanes. Results were unexpected relative to the engineering understanding of material performance in a reentry environment. The TEM analysis demonstrated that the oxides of aluminum and titanium mixed with silicon oxides to preserve a history of thermal conditions to which portions of the vehicle were exposed. The presence of Ti during the beginning of the deposition process, along with the thermodynamic phase precipitation upon cool down, indicate that temperatures well above the Ti melt point were experienced. The stratified observations implied that additional exothermic reaction, expectedly metal combustion of a Ti structure, had to be present for oxide formation. Results are significant for aerospace vehicles where thermal protection system (TPS) breaches cause substructures to be in direct path with the reentry plasma. 1
Reconstruction of flight aerothermal environments often requires the solution of an inverse heat transfer problem, which is an ill-posed problem of specifying boundary conditions from discrete measurements in the interior of the domain. This paper will present the algorithms implemented in the CHAR code for use in reconstruction of EFT-1 flight data and future testing activities. Implementation nuances will be discussed, and alternative hybrid-methods that are permitted by the implementation will be described. Results will be presented for a number of one-dimensional and multi-dimensional problems
Reconstruction of flight aerothermal environments often requires the solution of an inverse heat transfer problem, which is an ill-posed problem of determining boundary conditions from discrete measurements in the interior of the domain. This paper will present the algorithms implemented in the CHAR code for use in reconstruction of EFT-1 flight data and future testing activities. Implementation details will be discussed, and alternative hybrid-methods that are permitted by the implementation will be described. Results will be presented for a number of problems.
In this study, mathematical formulations to model charring ablation problems were numerically implemented using finite element analysis (FEA) with ABAQUS, which account for the material decomposition and progressive surface removal in the heat conduction and the surface energy balance equations. FEA was performed for a one-dimensional model to predict the temperature and ablation histories of a phenolic-impregnated carbon ablator sample (i.e., a common heat shield material for hypersonic vehicles and spacecraft) subjected to oxy-acetylene torch flame (i.e., 0.8 SLPM acetylene gas to 2.7 SLPM oxygen gas). The recovery enthalpy and convective heat transfer coefficient for the ablation model were calculated based on gas compositions and two assumed surface conditions (i.e., equilibrium and frozen). Simulations using the calculated recovery enthalpy and convective heat transfer coefficient resulted in a recession rate of 6.38 times (equilibrium) and 14.08 times (frozen) higher than the experimental data, despite fair agreement of the surface temperature. In addition, the effect of the heat transfer coefficient was investigated through a steady-state ablation analysis. The results of the analysis indicate that there is not one single value for the heat transfer coefficient that would allow the prediction to match both measured recession rate and surface temperature. Possible reasons for such an inconsistency are provided and discussed.
Computational modeling of Thermal Protection System (TPS) materials, used for aerospace applications, provides numerous advantages in preliminary selection and design of a heatshield material and shape for atmospheric entry vehicles. However, to serve as a reliable tool for prediction of material thermal and ablative behavior, the modeling approach needs to be validated against real experimental and flight data, preferably at a range of applied conditions. The validation study is typically very complex as it requires reliable measured data not only for the material thermal response and surface recession, but also well characterized environmental conditions. The validation problem becomes even more complex when the material thermal response is dictated by multi-physics effects such as solid conduction, in-depth thermal decomposition, pyrolysis gas flow and chemical reactions. The multi-physics effects complicate not only the modeling effort, but also the experimental measurement for validation of various aspects of the highly coupled problem. In this study, an attempt is made to identify suitable experimental data that could serve as a source for validation of material thermal response modeling tools. To reduce the computational complexity, this study focuses only on non-charring ablators, where the material thermal response could be modeled with a single governing equation for solid conduction and the ablation is limited only to the surface of the material. With a well characterized and publicly available experimental data being sparse, the study is limited in presenting test cases for only three materials: camphor, graphite and FiberForm® in the sequence of increased modeling complexity. Graphite is a commonly used TPS material for aerospace applications, both for leading edges of high-speed vehicles and internal insulation of solid rocket motors. FiberForm® is a porous carbon pre-form used in preparation of the well known PICA material Tran et al. [1996]. Inclusion of camphor into the list is conditioned with the relative simplicity in modeling the material thermal and chemical response and the low-enthalpy flow environment. In addition, camphor has been used as a simple test material for study of flow transition behavior by Stock and Ginoux [1973] and assessment of a heatshield shape change at flight relevant conditions by Rotondi et al. [2022]. In this work, the identified experimental data was extracted from the public literature and test cases that yet have been published. As it was found from the review, not a single test case contains an exhaustive set of data that would validate every aspect of the material physics. However, in the data collected, various aspects of the material behavior can be still validated, such as surface and in-depth temperature, amount of recession and a shape change. The identified experimental data for each case is accompanied with a characterized flow environment and simulated boundary conditions predicted by a Data-Parallel Line Relaxation (DPLR) code Wright et al. [1998]. In addition, material thermal response numerical simulations in each test case are performed with Kentucky Aerothermodynamics and Thermal Response System (KATS-MR) Zibitsker et al. [2022] providing a comparative study and a sanity check for the proposed validation data. Sample results from the performed numerical study are shown below. Figure 1 shows distribution of surface heat flux and pressure values on a hemi-cylinder model made of FiberForm® and tested in HyMETS arc-jet facility. The results are shown for the high pressure condition among the two tests. Flow simulation was performed with DPLR code on a quarter of original geometry. In the figure, the quarter shape was mirrored across zx and xy planes to show the complete distribution. Figure 2 shows the material response results for the high pressure case (7500 Pa), simulated with KATS-MR and a comparison to the experimental data for the surface temperature and shape shape. The simulation was performed on a 2-D slice, extracted in the xy plane at the middle of the sample. Figure 3 shows the material response simulation for the low pressure case (3500 Pa) and a comparison to the experimental data for surface temperature and shape change.