A multi-dimensional finite element solver for decomposing and non-decomposing ablating materials has recently been developed and is discussed in this paper. The underlying mathematical and material models are presented along with its discretization via the finite element method. The governing equations and solution algorithm is based on the one-dimensional control-volume finite element method (CVFEM) Chaleur code, a successful ablation code in use at Sandia National Labs, and this paper represents a multi-dimensional extension of Chaleur. The Equilibrium Surface Thermochemistry (EST) code, an equilibrium gas/surface thermochemistry code for decomposing and non-decomposing materials that was previously developed by the authors is used in conjunction with this new multi-dimensional ablation code to provide ablation thermochemistry information (i.e. B c and enthalpy tables). This new multi-dimensional ablation response code is first applied to solve two established code-to-code comparison problems with tabular aeroheating data. Another aspect of this work has been to develop the ability to couple CFD-based aeroheating data to the ablation code as a spatial and time variant boundary condition. Towards this end, we have established a one-way passing of aeroheating data from a hypersonic CFD code to the ablation code. We then examine the problem of simulating the ablation response of non-decomposing and decomposing materials in two arc-jet facilities.
NomenclatureRoman symbols B c dimensionless ablation rate B g dimensionless pyrolysis gas transfer rate C H dimensionless heat transfer coefficient C M dimensionless mass transfer coefficient ρ e u e C H dimensional heat transfer coefficient (kg/s − m 2 ) ρ e u e C M dimensional mass transfer coefficient (kg/s − m 2 ) (ρv) w boundary layer mass flux (kg/s − m 2 ) c v volumetric specific heat (J/kg − K) D elastic constitutive tensor E Arrhenius activation energy (J/kg) e internal energy (J/kg) h enthalpy (J/kg) k thermal conductivity (W/m − K) or Arrhenius pre-exponent factor (s −1 ) M pyrolysis gas molecular weighṫ m mass flux (kg/s − m 2 ) n unit normal p pressure (P a) Q heat source term (J/s − m 3 or W/m 3 ) This material is declared a work of the U.S. Government and is not subject to copyright protection in the United States. AIAA Aviation q heat flux (J/s − m 2 or W/m 2 ) R gas constant (J/kg − K) T temperature (K) t time (s) u displacement (m) v velocity (m/s) W finite element test function x spatial coordinate (m) Greek symbols dimensionless radiative emissivity or elastic strain Γ dimensionless volume fraction or surface of a domain κ permeability (m 2 ) µ dynamic viscosity (kg/s − m) Ω interior of a domaiṅ ω mass source term (kg/s − m 3 ) φ dimensionless porosity ρ density (kg/m 3 ) σ radiative Stefan-Boltzmann constant (W/m 2 − K 4 ) or elastic stress (N/m 2 ) Subscripts cchar quantity e boundary layer edge quantity g gas quantity i, j, k, l index notation identifiers ∞ far-field quantity v virgin quantity w wall quantity