# Behavior of Carbon-Carbon Composite under Intense Heating

Authors: R. Palaninathan
Journal: International Journal of Aerospace Engineering
 (2010)

---

## Abstract

This paper is concerned with modeling of ablation behavior of carbon-carbon composites used in hot spot areas of reentry space and hypersonic vehicles. Of the three modes of ablation (thermal, chemical and mechanical), the chemical (oxidation) is considered to influence the performance of the material. Aerodynamic heat flux need to be computed separately and is the input for this. The thermal field is obtained by 3D finite element method. Nonlinear transient thermal analysis is carried out, as the material properties are dependent on temperatures. Oxidation rates are computed using the analytical relations available in literature. The oxidation is divided into two regimes: reaction rate and diffusion rate controlled. Mainly the surface temperature controls the regime. The oxidation protected materials are considered by using the parameter “activation energy.” The variations of ambient temperature, pressure and oxygen concentration with altitude are taken into consideration. As the
recession takes place, newer surfaces are exposed to aerodynamic heating. Numerical examples are presented to show the effects of:
heat flux, altitude and oxidation protection on the recession characteristics. Change of regime from reaction to diffusion rate
control depends on parameters such as flow velocity and altitude. The latter has significant influence on ablation rate.


---

## Body

## 1. Introduction and a Brief Literature Review

Hypersonic and reentry vehicles are subjected to intense aerothermal loads. The loads consist of external surface pressure,
skin friction, and aerodynamic heating. The latter is the predominant one. The heating rates vary with the type of mission, as low as 0.57 to as high as 57 W/mm2. The corresponding radiation equilibrium temperatures are 1800 and 5610 K, respectively [1]. Thermal protection systems (TPS) made of suitable material are required so as to keep the excessive heat away from destroying or damaging the vehicle at a minimum weight penalty. There are basically five types of TPS [2] and the ablative is the one that is more effective in expending thermal energy. Ablation is a process by which thermal energy is expended by sacrificing the material and thereby, absorbing a considerable amount of energy. There are two types of ablative materials, namely, melting and nonmelting types [3]. In melting type (thermoplastics), the heat is expended in transformation of material phase from solid to liquid. The liquid on the surface is removed immediately after formation by aerodynamic surface shear, and newer surface is exposed to heating. In non-melting type, there are two: low temperature ablators (LTAs) and high temperature ablators (HTAs). The LTAs are of fiber reinforced composites, with thermosetts as the binder. This material decomposes at low temperatures, about 600 K, into partly char and partly gases. The char remains on the surface and offers thermal protection. However, the char layers are removed by the aerodynamic surface shear very fast, as the reduction of strength of char at high temperature is substantial. This type of erosion of surface is termed as “mechanical ablation” [4]. These are the materials suitable for ballistic missiles, “high heating rate-shorter duration.” The examples of HTAs are carbon-carbon (C/C) and carbon-silicon carbide ceramic matrix composites. These materials possess double bond in their molecular structure, and it takes a lot of energy to break the bond when they decompose at high temperatures. This process is termed as “chemical ablation.” In order to improve the ablative performance, oxidation protection methods have been developed. In view of this, HTAs have been used in the recent technological applications: (i) manned missions, (ii) hypersonic vehicles, (iii) missiles, at the nose tips and leading edges. This paper concerns with ablation modeling of C/C materials, which is presented in three steps: (i) nonlinear transient thermal analysis to compute the temperature field for assumed heat flux, (ii) chemical ablation to predict the surface recession and (iii) remeshing by moving the surface nodes to the extent of recession at the current time step, before going to the next time step.A review of fifty years of hypersonic vehicles has been presented [5] describing the improvements in structures, materials, and manufacturing technologies which include the development of light weight high strength materials that can be effectively manufactured into structures, capable of withstanding the hypersonic environment. A great interest has been focused on the study of oxidation behavior of carbon, spread over about 70 years [6]. Their findings can be summarized as follows. The oxidation is divided into two regimes: (i) chemical reaction rate controlled, at low temperature, reaction rate limits oxidation process as the quantity of oxygen diffusing through the boundary is more than what is required to sustain the reaction at its own pace; CO2 is the only product of oxidation; (ii) diffusion rate controlled, the availability of oxygen diffusing through boundary layer limits the oxidation process and there is not enough oxygen for oxidation to grow at its own pace; both CO andCO2 are products of oxidation. It has been determined that bare C/C composites cannot be used at temperatures above 800 K for extended period of time, because of erosion due to oxidation. Many approaches have been considered to protect carbon from oxidation, which can be classified into: (i) use of protective coatings. In this, thin layers of oxidation resistance coatings are applied on the surface of bare material. The coatings are usually based on SiO2, Al2O3, Cl2O3. At high temperatures, these coatings decompose to form protective oxide scale. The major problem with the coatings is the cracking due to mismatch of thermal properties [7]. (ii) use of matrix modifiers. Materials like boron and phosphorous are added to the matrix during fabrication and they get embedded within the material system. On exposures to high temperature, these modifiers react with oxygen and form oxides (B2O3) which inhibits/retards the oxidation [8]. Comparative studies on protective coating and matrix modification have been carried out [9, 10]. There are several texts on thermal modeling of solid bodies using finite element method [11].


## 2. Formulation

### 2.1. Finite Element Thermal Modeling

As mostlyC/C composites are solid bodies (nose tips and leading edges), 3-dimensional solid elements are used. The nature of variations of the field variables (temperature) and the geometry determine the selection of number of nodes in the element. Figure 1 is a variable node 3d solid element, suitable for both quadratic and linear variations of field variables. The temperature and the geometry are expressed using the same approximating shape functions, Ni  [11]:(1)T(x,y,z,t)=∑i=1nNi(ε,η,ς)Ti(t)(2)(x,y,z)=∑i=1nNi(ε,η,ς)(xi,yi,zi),
where ε,η,ς are spatial natural coordinates, Ti are the time-dependent nodal temperatures, x,y,z are the global coordinates of any point within an element, and xi, yi, zi are the nodal coordinates. Using Fourier’s law, the component of heat flow rate for an anisotropic medium can be written in matrix form(3){qxqyqz}=-[kxxkxykxzkyxkyykyzkzxkzykzz]{T,xT,yT,z}=-[kxyz][B]{Ti},
where [kxyz] is the material conductivity matrix in global directions expressed as(4)[kxyz]=[λ]T[k123][λ]
where [k123] is the conductivity matrix in principal material directions, see Figure 2,  and [λ] is the matrix of direction cosines. [B] is the product of matrices of Jacobian inverse and derivatives of shape functions w.r.t natural coordinates. Applying the method of weighted residuals and using shape functions as weighting function, the governing equation is obtained as [12]
(5)[C]{Ṫ}+[[Kc]+[Kr]]{T}={Fq}+{Fr}.
The associated boundary conditions are (i) specified heat flux, q, and (ii) surface radiation, σε(Ts-T∞)4·Ts: surface temperature, T∞: medium temperature, σ: Stephen-Boltzmann constant and ε: emissivity. [C] is the element capacitance matrix. [Ke] and [Kr] are conductance matrices corresponding to material conductivity and radiation, respectively. {Fq} and {Fr} are the load vectors due to applied heat flux and radiation, respectively. As high temperatures are involved, the material properties are functions of temperature and hence function of time. The solution to this nonlinear transient problem in time domain is obtained in small time steps. The variation of the matrices over the time interval, Δt is assumed to be linear. Using residual method and minimization of error over time interval, Δt, equation (5) is transformed to:
(6)[K]{ΔTi}+{R}={0},
where {ΔTi}: incremental nodal temperatures and {R}: residual load. Equation (6) is solved by iterative process till convergence is achieved. The converged values are added to the accumulated values up to the previous time step.Variable node 3D solid element.
(a)
8 nodes(b)
12 nodes(c)
21 nodesFigure 2
3DC/C composites.


### 2.2. Surface Recession

Figure3 is the schematic of surface oxidation and boundary movement. As mentioned earlier, the oxidation is divided into two regimes, reaction rate and diffusion rate regimesChemical ablation and boundary movement.
(a)(b)(c)


#### 2.2.1. Reaction Rate Controlled Regime

The molar oxidation rate is given by [13](7)r=30.7e(-Ea1.99Ts)CO2molcm2s,
where Ts: surface temperature, CO2: concentration of oxygen, which is a function of altitude, Ea: activation energy, the value of which changes depending on the level of oxidation protection, beingequal to about 35 kcal/mole for unprotected and to about 56 kcal/mole for protected materialThe Oxidation rate is given by(8)Kr=r×12gcm2s.


#### 2.2.2. Diffusion Rate Controlled Regime

The oxidation rate is given by [13]
(9)Kd=12(2φ+1φ+1)Ngcm2s,
where φ, ratio of CO and CO2 produced, is given by
(10)φ=(1+0.3215e(868/Ts)-1)/0.7656e(868/Ts)0.21-(1+0.3215e(868/Ts)-1)/0.7656e(868/Ts)
and N, the diffusion rate of oxygen through boundary layer, is given by
(11)N=2.07×10-6PpiPogRTs0.8d(2422.77(ρdu)0.7Ts0.55+3.95),
where P: total pressure, Pog: pressure of oxygen in the bulk air, pi: logarithmic pressure of inert, R: universal gas constant, d: characteristic dimension, u: flow velocity of air around the body, and ρ: density of air. Several of the above parameters are functions of altitude.Figure4 is the plot of oxidation rates (as per the two regimes, (8) and (9)) as function of temperatures [13]. At temperatures <1150K, the reaction rate regime controls and at higher temperatures,the diffusion rate regime controls the oxidation. The recession rate is obtained from the lower of the two as
(12)dRdt=Kρc,
where ρc: density of material, gm/cm3. If carbon fiber and carbon matrix do not oxidize simultaneously, the value of ρc needs to be corrected accordingly.Figure 4
Oxidation behavior of C/C composite.


### 2.3. Boundary Movement Scheme

With reference to Figures3 and 4, the recession depth at a surface node (ΔRc) is obtained by the product of recession rate and the time increment Δt. The boundary nodes are moved inward along the normal to the surface, by an amount equal to ΔRc, The new position of each boundary node is given by
(13)newx=x-ΔRc×lx,newy=y-ΔRc×ly,newz=y-ΔRc×lz,
where lx, ly, lz are the direction cosines. Also the mid-side nodes of the elements lying on the surface are moved by half the distance, ΔRc/2. The nodal positions of all other inner element are not changed. The temperatures of the   nodes at their new positions, (ξT=1-2/h*ΔRC), are obtained by interpolation, assuming quadratic variations; see Figure 5. During subsequent time steps, if the recession depth is sufficiently large, the surface elements and the nodes are removed from FE idealization by making them dummy.Figure 5
Boundary movement Scheme.


## 2.1. Finite Element Thermal Modeling

As mostlyC/C composites are solid bodies (nose tips and leading edges), 3-dimensional solid elements are used. The nature of variations of the field variables (temperature) and the geometry determine the selection of number of nodes in the element. Figure 1 is a variable node 3d solid element, suitable for both quadratic and linear variations of field variables. The temperature and the geometry are expressed using the same approximating shape functions, Ni  [11]:(1)T(x,y,z,t)=∑i=1nNi(ε,η,ς)Ti(t)(2)(x,y,z)=∑i=1nNi(ε,η,ς)(xi,yi,zi),
where ε,η,ς are spatial natural coordinates, Ti are the time-dependent nodal temperatures, x,y,z are the global coordinates of any point within an element, and xi, yi, zi are the nodal coordinates. Using Fourier’s law, the component of heat flow rate for an anisotropic medium can be written in matrix form(3){qxqyqz}=-[kxxkxykxzkyxkyykyzkzxkzykzz]{T,xT,yT,z}=-[kxyz][B]{Ti},
where [kxyz] is the material conductivity matrix in global directions expressed as(4)[kxyz]=[λ]T[k123][λ]
where [k123] is the conductivity matrix in principal material directions, see Figure 2,  and [λ] is the matrix of direction cosines. [B] is the product of matrices of Jacobian inverse and derivatives of shape functions w.r.t natural coordinates. Applying the method of weighted residuals and using shape functions as weighting function, the governing equation is obtained as [12]
(5)[C]{Ṫ}+[[Kc]+[Kr]]{T}={Fq}+{Fr}.
The associated boundary conditions are (i) specified heat flux, q, and (ii) surface radiation, σε(Ts-T∞)4·Ts: surface temperature, T∞: medium temperature, σ: Stephen-Boltzmann constant and ε: emissivity. [C] is the element capacitance matrix. [Ke] and [Kr] are conductance matrices corresponding to material conductivity and radiation, respectively. {Fq} and {Fr} are the load vectors due to applied heat flux and radiation, respectively. As high temperatures are involved, the material properties are functions of temperature and hence function of time. The solution to this nonlinear transient problem in time domain is obtained in small time steps. The variation of the matrices over the time interval, Δt is assumed to be linear. Using residual method and minimization of error over time interval, Δt, equation (5) is transformed to:
(6)[K]{ΔTi}+{R}={0},
where {ΔTi}: incremental nodal temperatures and {R}: residual load. Equation (6) is solved by iterative process till convergence is achieved. The converged values are added to the accumulated values up to the previous time step.Variable node 3D solid element.
(a)
8 nodes(b)
12 nodes(c)
21 nodesFigure 2
3DC/C composites.


## 2.2. Surface Recession

Figure3 is the schematic of surface oxidation and boundary movement. As mentioned earlier, the oxidation is divided into two regimes, reaction rate and diffusion rate regimesChemical ablation and boundary movement.
(a)(b)(c)


### 2.2.1. Reaction Rate Controlled Regime

The molar oxidation rate is given by [13](7)r=30.7e(-Ea1.99Ts)CO2molcm2s,
where Ts: surface temperature, CO2: concentration of oxygen, which is a function of altitude, Ea: activation energy, the value of which changes depending on the level of oxidation protection, beingequal to about 35 kcal/mole for unprotected and to about 56 kcal/mole for protected materialThe Oxidation rate is given by(8)Kr=r×12gcm2s.


### 2.2.2. Diffusion Rate Controlled Regime

The oxidation rate is given by [13]
(9)Kd=12(2φ+1φ+1)Ngcm2s,
where φ, ratio of CO and CO2 produced, is given by
(10)φ=(1+0.3215e(868/Ts)-1)/0.7656e(868/Ts)0.21-(1+0.3215e(868/Ts)-1)/0.7656e(868/Ts)
and N, the diffusion rate of oxygen through boundary layer, is given by
(11)N=2.07×10-6PpiPogRTs0.8d(2422.77(ρdu)0.7Ts0.55+3.95),
where P: total pressure, Pog: pressure of oxygen in the bulk air, pi: logarithmic pressure of inert, R: universal gas constant, d: characteristic dimension, u: flow velocity of air around the body, and ρ: density of air. Several of the above parameters are functions of altitude.Figure4 is the plot of oxidation rates (as per the two regimes, (8) and (9)) as function of temperatures [13]. At temperatures <1150K, the reaction rate regime controls and at higher temperatures,the diffusion rate regime controls the oxidation. The recession rate is obtained from the lower of the two as
(12)dRdt=Kρc,
where ρc: density of material, gm/cm3. If carbon fiber and carbon matrix do not oxidize simultaneously, the value of ρc needs to be corrected accordingly.Figure 4
Oxidation behavior of C/C composite.


## 2.2.1. Reaction Rate Controlled Regime

The molar oxidation rate is given by [13](7)r=30.7e(-Ea1.99Ts)CO2molcm2s,
where Ts: surface temperature, CO2: concentration of oxygen, which is a function of altitude, Ea: activation energy, the value of which changes depending on the level of oxidation protection, beingequal to about 35 kcal/mole for unprotected and to about 56 kcal/mole for protected materialThe Oxidation rate is given by(8)Kr=r×12gcm2s.


## 2.2.2. Diffusion Rate Controlled Regime

The oxidation rate is given by [13]
(9)Kd=12(2φ+1φ+1)Ngcm2s,
where φ, ratio of CO and CO2 produced, is given by
(10)φ=(1+0.3215e(868/Ts)-1)/0.7656e(868/Ts)0.21-(1+0.3215e(868/Ts)-1)/0.7656e(868/Ts)
and N, the diffusion rate of oxygen through boundary layer, is given by
(11)N=2.07×10-6PpiPogRTs0.8d(2422.77(ρdu)0.7Ts0.55+3.95),
where P: total pressure, Pog: pressure of oxygen in the bulk air, pi: logarithmic pressure of inert, R: universal gas constant, d: characteristic dimension, u: flow velocity of air around the body, and ρ: density of air. Several of the above parameters are functions of altitude.Figure4 is the plot of oxidation rates (as per the two regimes, (8) and (9)) as function of temperatures [13]. At temperatures <1150K, the reaction rate regime controls and at higher temperatures,the diffusion rate regime controls the oxidation. The recession rate is obtained from the lower of the two as
(12)dRdt=Kρc,
where ρc: density of material, gm/cm3. If carbon fiber and carbon matrix do not oxidize simultaneously, the value of ρc needs to be corrected accordingly.Figure 4
Oxidation behavior of C/C composite.


## 2.3. Boundary Movement Scheme

With reference to Figures3 and 4, the recession depth at a surface node (ΔRc) is obtained by the product of recession rate and the time increment Δt. The boundary nodes are moved inward along the normal to the surface, by an amount equal to ΔRc, The new position of each boundary node is given by
(13)newx=x-ΔRc×lx,newy=y-ΔRc×ly,newz=y-ΔRc×lz,
where lx, ly, lz are the direction cosines. Also the mid-side nodes of the elements lying on the surface are moved by half the distance, ΔRc/2. The nodal positions of all other inner element are not changed. The temperatures of the   nodes at their new positions, (ξT=1-2/h*ΔRC), are obtained by interpolation, assuming quadratic variations; see Figure 5. During subsequent time steps, if the recession depth is sufficiently large, the surface elements and the nodes are removed from FE idealization by making them dummy.Figure 5
Boundary movement Scheme.


## 3. Numerical Examples

Numerical examples are presented to bring out the effects of: reradiation, heating rate, altitude, and oxidation protection. The problem considered is that of rectangular slab of thickness equal to 150 mm, assumed to be made ofC/C composite. It is subjected to specified heat flux on one surface, in thickness direction. All other (five) surfaces are insulated; that is, no boundary condition is specified. This is essentially one-dimensional problem. However it is solved using 3d solid elements, 12 nodes per element with 4 midside nodes in thickness direction to account for nonlinear temperature profiles. The initial temperature is taken as 300 K. On the heated surface, also radiation boundary condition is applied to consider the heat loss from the hot surface due to reradiation. The properties of 3D C/C composites as functions of temperature are obtained from [14].


### 3.1. Effect of Surface Radiation

The heat transfer from the hot body by radiation is governed by fourth power law and is usually neglected in many cases as seen from the literature. However, in the case of aerospace vehicles subjected to aerodynamic heating of higher intensity, the radiation from the hot body needs to be considered. In order to show this effect, problem is analysed with and without radiation boundary condition on the heated surface. Figure6 shows variation of surface temperature with time for an assumed heat flux, q=0.45W/mm2. Oxidation and surface recession are not considered. The surface temperature goes up to 2300 K at t=500seconds for the case without radiation boundary condition, when this condition is included the temperature goes up to 1482 K only. This difference has significant influence on the oxidation and recession rates. These values are higher, when radiation not included. Subsequent examples are presented with radiation included.Figure 6
Surface temperature versus time(q=0.45W/mm2,u=19000cm/s).


### 3.2. Effect of Incident Heat Intensity

In aerospace vehicles, the intensity of heat flux varies as function of surface coordinates and time with maximum value at the stagnation point, as the vehicle moves in its flight path. Here, heat fluxes,q=0.25, 0.50, 0.75 and 1.00 W/mm2 assumed uniformly distributed on the surface and constant with time, are considered with oxidations at sea level and at high altitudes. Figure 7 shows the recession depth as a function of time at sea level. For the highest heat flux considered (q=1.0 W/mm2), the recession depth at t=500 seconds is 70 mm, that is, about half the thickness of slab has been eroded due to oxidation; whereas for q=0.25 W/mm2, the recession is negligibly small. Figure 8 shows the recession depth at altitude equal to 32.5 km. The recession reduces considerably to 3.5 mm at t=500 seconds for heat intensity q=1.00 W/mm2 against 70 mm at sea level, with other conditions remaining same.Figure 7
Surface recession with time at sea level,h=0.Figure 8
Surface recession with time at altitude,h=32.5 km.


### 3.3. Effect of Altitude

The altitude has greater influence on the oxidation. The variations of ambient temperature, pressure and density of air with altitude are the influencing factors. A study has been carried out for three altitudes,h=0, 10 and 30 km (body stationary at these altitudes). Heat flux applied is 0.45 W/mm2. Figure 9 shows the variation of recession depth with time. At t=500 seconds the depth is 9.4 mm at sea level which reduces to 0.2 mm at altitude equal to 30 km.Figure 9
Recession depth at different altitudes(q=0.45W/mm2).


### 3.4. Effect of Oxidation Protection

The results presented in the earlier three sections are for the cases of unprotectedC/C composites. In this section, the results are presented for protected material, obtained by matrix modification [15].The difference between the protected and unprotected materials comes through the values of activation energy, Ea in (7). The activation energy increase from 37500 to 56000 cal/mole is assumed for protected material. Figures 10 and 11 show the effects of oxidation protection at sea level and at altitude, h=32.5 km. At sea level the recession depth reduces from 9.4 mm for unprotected material to 2.6 mm for protected material. At higher altitude, the influence is much more and the reduction is from 0.13 to 0.03 mm.Figure 10
Effect of Oxidation protection at sea level,h=0.Figure 11
Effect of Oxidation protection at higher altitude,h=32.5km.


## 3.1. Effect of Surface Radiation

The heat transfer from the hot body by radiation is governed by fourth power law and is usually neglected in many cases as seen from the literature. However, in the case of aerospace vehicles subjected to aerodynamic heating of higher intensity, the radiation from the hot body needs to be considered. In order to show this effect, problem is analysed with and without radiation boundary condition on the heated surface. Figure6 shows variation of surface temperature with time for an assumed heat flux, q=0.45W/mm2. Oxidation and surface recession are not considered. The surface temperature goes up to 2300 K at t=500seconds for the case without radiation boundary condition, when this condition is included the temperature goes up to 1482 K only. This difference has significant influence on the oxidation and recession rates. These values are higher, when radiation not included. Subsequent examples are presented with radiation included.Figure 6
Surface temperature versus time(q=0.45W/mm2,u=19000cm/s).


## 3.2. Effect of Incident Heat Intensity

In aerospace vehicles, the intensity of heat flux varies as function of surface coordinates and time with maximum value at the stagnation point, as the vehicle moves in its flight path. Here, heat fluxes,q=0.25, 0.50, 0.75 and 1.00 W/mm2 assumed uniformly distributed on the surface and constant with time, are considered with oxidations at sea level and at high altitudes. Figure 7 shows the recession depth as a function of time at sea level. For the highest heat flux considered (q=1.0 W/mm2), the recession depth at t=500 seconds is 70 mm, that is, about half the thickness of slab has been eroded due to oxidation; whereas for q=0.25 W/mm2, the recession is negligibly small. Figure 8 shows the recession depth at altitude equal to 32.5 km. The recession reduces considerably to 3.5 mm at t=500 seconds for heat intensity q=1.00 W/mm2 against 70 mm at sea level, with other conditions remaining same.Figure 7
Surface recession with time at sea level,h=0.Figure 8
Surface recession with time at altitude,h=32.5 km.


## 3.3. Effect of Altitude

The altitude has greater influence on the oxidation. The variations of ambient temperature, pressure and density of air with altitude are the influencing factors. A study has been carried out for three altitudes,h=0, 10 and 30 km (body stationary at these altitudes). Heat flux applied is 0.45 W/mm2. Figure 9 shows the variation of recession depth with time. At t=500 seconds the depth is 9.4 mm at sea level which reduces to 0.2 mm at altitude equal to 30 km.Figure 9
Recession depth at different altitudes(q=0.45W/mm2).


## 3.4. Effect of Oxidation Protection

The results presented in the earlier three sections are for the cases of unprotectedC/C composites. In this section, the results are presented for protected material, obtained by matrix modification [15].The difference between the protected and unprotected materials comes through the values of activation energy, Ea in (7). The activation energy increase from 37500 to 56000 cal/mole is assumed for protected material. Figures 10 and 11 show the effects of oxidation protection at sea level and at altitude, h=32.5 km. At sea level the recession depth reduces from 9.4 mm for unprotected material to 2.6 mm for protected material. At higher altitude, the influence is much more and the reduction is from 0.13 to 0.03 mm.Figure 10
Effect of Oxidation protection at sea level,h=0.Figure 11
Effect of Oxidation protection at higher altitude,h=32.5km.


## 4. Concluding Remarks

Finite element method-based software with boundary movement capability for modeling of oblation behavior of carbon-carbon composites has been developed. Analytical expressions for oxidation rates have been taken from literature. Based on the numerical studies on example problems, a few conclusions are drawn. The change of regime from reaction to-diffusion-rate controlled oxidation takes place at different temperatures depending on parameters such as velocity, altitude. The altitude has got significant influence on the recession. The protected material offers much better oxidation resistance compared to the unprotected material. It may be observed that in the present work, the oxidation protection has been modeled using the parameter, “activation energy.” This term appears only in the expression for the reaction rate controlled regime. However, a similar term for the protected material is not present in the expression for the diffusion controlled regime. This calls for improvement in the presently available model. With this software available, it is possible to model the ablation behavior of nose tips and leading edges used in the reentry space vehicles with aerodynamic loads computed separately. However, if the present software is integrated with a suitable aerothermodynamic package, it becomes possible to cocompute the aerodynamic loads along with the shape change due to ablation, taking in to account the varying local atmospheric data, as the vehicle moves in its trajectory.



---

*Source: 257957-2010-04-29.xml*