Ablative materials represent a low cost and reliable means to insulate rocket engine components from high-temperature, corrosive combustion product environments. Besides their diffuse application in solid rocket nozzles, their use also emerges as a valid alternative in liquid rocket engines. Together with the growing interest in oxygen/methane liquid rocket engines, these materials have gained attention as possible insulator for small upperstage engines or in-space thrusters. In this framework, a validated approach for the study of carbon-based pyrolyzing and non-pyrolyzing materials, together with a novel boundary condition developed to analyze the silica-based material behavior, has been used to numerically reproduce the material response in the highly oxidizing environment generated by the combustion of oxygen and methane. At first, the validation against experimental data of the silica-based material erosion model is presented. Subsequently, the behavior and the response of different ablators in a oxygen/methane environment is numerically investigated for a wide range of operating conditions. Commonly made assumptions in the simulation of the material response are thoroughly analyzed and a critical overview of the results is presented.
Nomenclaturebinary diffusion coefficient, m 2 /s h = enthalpy, J/kg k = thermal conductivity, W/m · K N c = number of specieṡ m = mass blowing rate per unit area, kg/m 2 · s p = pressure, N/m 2 q radout = wall radiative heat flux, W/m 2 s = erosion rate, m/s T = temperature, K t = time, s v = flow velocity, m/s v = velocity component normal to surface, m/s * Ph.D. Student, Dipartimento di Ingegneria Meccanica e Aerospaziale, Via Eudossiana 18. AIAA Student Member. alessandro.turchi@uniroma1.it. ẇ i = species source term in control volume, kg/m 3 · s x = mole fraction y = mass fraction y + = dimensionless wall distance Subscripts c = charred state f = failed material g = pyrolysis gas i = species s = solid state v = virgin state w = gas properties at gas/solid interface Symbols η = outward coordinate normal to surface Γ m = mass resin fraction Γ V = volume resin fraction µ = viscosity, kg/m · s µ t = turbulent viscosity, kg/m · s ρ = density, kg/m 3 ρ A = primary resin component density, kg/m 3 ρ B = secondary resin component density, kg/m 3 ρ R = reinforcing material density, kg/m 3