Microstructure-Based Thermochemical Ablation Model of Carbon/Carbon-Fiber Composites

The microstructure of carbon fiber–reinforced carbon-matrix composites (carbon/carbon composites) has important effects on its ablation performance. However, the traditional macro-ablation methods have underestimated the ablation recession rate and ignored the influence of microstructure. To simulate the ablation of large-sized structures while accounting for the influence of microstructure, it is necessary to modify these methods. In this work, a thermochemical ablation model for carbon/carbon composites is proposed based on the evolution behavior of their microstructure. The ablation recession rate and surface temperature predicted by this model are in good agreement with the experimental results. Through numerical analysis, we found that the ablation recession rate of the material without carbon fibers is much greater than that of the material containing carbon fibers. The ablation recession rate is influenced by the fiber orientation due to the change in thermal conductivity. The anti-ablation efficiency of carbon/carbon composites can be improved by increasing their fiber radius, radiation coefficient, specific heat capacity, interphase density, and thermal conductivity coefficient. The thermochemical ablation model provides a guide for the design of better anti-ablation carbon/carbon composites.


Introduction
Carbon fiber-reinforced carbon-matrix composites (carbon/carbon composites), which possess high latent heat of sublimation, excellent strength, and chemical stability at high temperatures, have been widely used as thermal protection materials in atmospheric re-entry vehicle heat shields and rocket engine nozzles [1][2][3][4][5][6][7][8][9][10][11]. High-temperature, highpressure, and high-velocity airflow are present in the typical service environments of carbon/carbon composites. Ablation, which is a material consumption behavior of carbon/carbon composite surfaces in service environments, is an effective protection mechanism for aircraft components [12][13][14]. The thermochemical ablation behavior of carbon/carbon composites involves the heat transfer, chemical reactions, surface recession and mass loss, which seriously affects some performance parameters of aircraft. Among them, the ablation recession rate is a key parameter for the characterization of surface recession behavior and the thickness design of thermal protection materials. With the progress of weaving technologies, the microstructure of carbon/carbon composites is becoming more and more complex. It changes not only the volume fraction of carbon fiber but also the fiber orientation and void volume fraction. Due to the different characteristics of the fibers and matrix [15,16], these microstructural changes have a great impact on the thermochemical ablation behavior. To reveal the ablation mechanisms, it is important to investigate the thermochemical ablation behavior of carbon/carbon composites at the micro scale.
Numerous researchers have focused on the thermochemical ablation behavior of carbon/carbon composites under high-temperature airflow [17][18][19], especially on its analysis models. These studies can generally be divided into the study of: (1) heat transfer analysis [20][21][22], (2) detailed chemical reactions [23][24][25], (3) surface ablation morphology [26,27], and (4) the effects of the material microstructure on ablation characteristics [28][29][30]. In the aspect of heat transfer analysis, the thermal response was simulated from the perspective of the matrix and fiber components of various braided carbon/carbon composites [31,32]. The results showed that the temperature distributions of the matrix and fiber, which presented different characteristics during ablation, were not uniform along the thickness direction of the material structure. In the aspect of detailed chemical reactions, the thermochemical ablation products were commonly studied based on the principle of minimizing the free energy and microscopic analysis by X-ray photoelectron spectroscopy (XPS) and X-ray diffraction (XRD) [33,34]. Thermogravimetric analysis (TGA) was used to determine the oxidation kinetics of carbon/carbon composites in oxidizing conditions [35,36]. Based on the analysis of detailed chemical reactions, the combined numerical analysis models of the fluid, thermal, and physicochemical reactions during the ablation process were established and used to predict/analyze the thermal response and ablation properties of the material [37]. Among them, the traditional macro-ablation model of physicochemical reactions assumes that the material surface is a flat plate and only considers the diffusion and reaction of gas in the boundary layer. It underestimates the ablation recession rate and ignores the influence of microstructure, which is inconsistent with the actual rough material surface. In the aspect of surface ablation morphology, the ablation microstructure of carbon/carbon composites has been studied by scanning electron microscopy (SEM) and microtomography [38,39]. These studies revealed the existence of a thin reactive interphase between the fibers and the matrix. The ablation recession rate of the interphase was higher than that of the fibers and the matrix. Based on the above phenomenon, Lachaud et al. [40] established a steady-state ablation morphology model for carbon/carbon composites to study the effects of reactivity and gas diffusion on the ablation morphology. The results might be beneficial for predicting of the surface roughness, predicting of the effective reactivity of the composite components and the inverting of the extremely harsh environment of composites. However, the morphology made it difficult to obtain specific data and to use as an index to evaluate the ablation performance of the materials, which limited the development of the morphology research. In the aspect of material microstructure, the effect of density and fiber orientation on the ablation behavior of carbon/carbon composites was investigated. The results showed that the fiber architecture and the density had more effects on the ablation behavior than the open porosity of the composites [41]. Liu et al. [42] investigated the impact of temperature and fiber volume fraction on the ablation performance of the carbon/carbon composites using numerical simulations. However, only the influence of the temperature was considered by the element exclusion criteria in this method, and the effect of gas diffusion and other factors were neglected. In fact, the temperature differential is very small at the microscopic scale. Li et al. [43] analyzed the influence of the braided angle, porosity, and reactivity contrast between the matrix and yarns on the mesoscopic thermochemical ablative behaviors. These works only considered the mesoscopic structure in the geometric modeling of materials but not in the theoretical model. It meant the models could only be used for material ablation simulation but not for large structures.
In summary, the thermochemical ablation recession rates of carbon/carbon composites are closely related to the properties and arrangement orientation of the components, the morphological characteristics of the reinforcing phase, and the interfacial properties between the reinforcing phase and the matrix. The traditional macro-ablation model only considers mesoscopic structure in geometric modeling, so it is impossible to consider the influence of detailed microstructure on ablation performance in large-sized structures. The microstructure such as fiber orientation has a significant effect on the ablation performance, it is necessary to modify the thermochemical ablation model to take the microscopic characteristics into account. Therefore, an effective microstructure-based thermochemical ablation analysis model for carbon/carbon composites is proposed in this work to study thermochemical ablation behavior. This paper is structured as follows. In Section 2, we analyze the ablation mechanism and proposed a microstructure-based thermochemical ablation model of carbon/carbon composites. In Section 3, we introduce the calculation process of the model. The model was validated by comparison to the experimental data of Reference [25]. Section 4 discusses the effects of the microstructure, material properties and external environment on the surface ablation recession rate. Section 5 concludes the main results. The research law of microstructure may be applied to the prediction of the macro-ablation performance of ablative materials in this paper, and this work may provide a reference for the design and preparation of better high-ablation resistance carbon/carbon composites.

Problem Statement
While flying at ultra-high speed in the atmosphere, an aircraft interacts violently with the surrounding air, creating a strong shock wave ahead, as shown in Figure 1a. The airflow is sharply compressed after the shock wave. The kinetic energy of the airflow is converted into internal energy of the air molecules, which heats the airflow to an ultra-high temperature. The flow of heat in the air acts on the material surface through convective and radiative heat transfer, and a part of the heat is transferred to the material. This creates a temperature gradient inside the material in which the surface temperature is higher and the internal temperature is lower, as shown in Figure 1b. In a high-temperature airflow environment, the surface material reacts with the air in a complex manner and becomes very rough. This process is called thermochemical ablation, which results in quality loss and surface recession.  Figure 2 shows the typical morphology of carbon/carbon composites before and after oxidation in air. Based on the main characteristics of the fiber composites shown in Figure  2, we create an abstract analysis model shown in Figure 1c,d for the numerical simulation. The morphology of the channels in the SEM image indicates that ablation occurs at the  The carbon/carbon composite is a gasification-type ablation material that exhibits complex thermochemical ablation characteristics under high-temperature airflow. The thermochemical ablation of the carbon/carbon composites involves complex physicochemical reactions and gas diffusion, such as oxidation reactions, carbon-nitrogen reactions, sublimation, and oxidation diffusion, as shown in Equations (1)-(5): CO(g) is the major ablation product of the carbon/carbon composites at temperatures below 3300 K because of the insufficient oxygen on the surface of the material [33]. Equation (1) describes a spontaneous chemical reaction at high temperatures; for example, because ∆G • (CO) = −654.726 kJ·mol −1 at 2500 K, this reaction can be considered as irreversible because the equilibrium partial pressure ratio between the reactant and product is P 2 CO /P O 2 = 4.796 × 10 13 Pa. Consequently, the thermochemical ablation model in this study considers only the CO(g) formation reaction. Figure 2 shows the typical morphology of carbon/carbon composites before and after oxidation in air. Based on the main characteristics of the fiber composites shown in Figure 2, we create an abstract analysis model shown in Figure 1c,d for the numerical simulation. The morphology of the channels in the SEM image indicates that ablation occurs at the fiber-matrix interface. In other words, the fiber-matrix interfaces are much easier to oxidize than the fiber and matrix. Under high-temperature airflow, the oxygen in the boundary layer impinges and reacts with the carbon on the material surface. Then, the CO gas product escapes from the material surface into the boundary layer. As the interphase continuously retreats, gas diffusion channels are formed. The channels expose the fibers and matrix to the heat airflow. As shown in Figure 1c, oxygen is transported from the boundary layer to the ablation surfaces through the gas diffusion channels. With the growth of the gas diffusion channels, the oxidation rate of the interphase is limited by gas diffusion. As a result, the thermochemical ablation mechanisms are diffusion-limited and/or reaction-limited. In addition, as shown in Figure 1d, the fiber orientation affects the fiber morphology and influences the gas diffusion.  Figure 2 shows the typical morphology of carbon/carbon composites before and after oxidation in air. Based on the main characteristics of the fiber composites shown in Figure  2, we create an abstract analysis model shown in Figure 1c,d for the numerical simulation. The morphology of the channels in the SEM image indicates that ablation occurs at the fiber-matrix interface. In other words, the fiber-matrix interfaces are much easier to oxidize than the fiber and matrix. Under high-temperature airflow, the oxygen in the boundary layer impinges and reacts with the carbon on the material surface. Then, the CO gas product escapes from the material surface into the boundary layer. As the interphase continuously retreats, gas diffusion channels are formed. The channels expose the fibers and matrix to the heat airflow. As shown in Figure 1c, oxygen is transported from the boundary layer to the ablation surfaces through the gas diffusion channels. With the growth of the gas diffusion channels, the oxidation rate of the interphase is limited by gas diffusion. As a result, the thermochemical ablation mechanisms are diffusion-limited and/or reaction-limited. In addition, as shown in Figure 1d, the fiber orientation affects the fiber morphology and influences the gas diffusion.

Thermochemical Ablation Model
The thermochemical ablation mechanisms are diffusion-limited and/or reactionlimited. The ablation recession rate of the interphase v i can be obtained from the reactivity k i (m·s −1 ) and the molar concentration of oxygen C i (mol·m −3 ) at the surface of the interphase [40]: where Ω i is the solid molar volume of interphase (m 3 ·mol −1 ), J C is the molar rate of ablation (mol·m −2 ·s −1 ). Similarly, the ablation recession rate of the fiber v f (z) can be written as where the coordinates in the thermochemical ablation model are micro scale, while the coordinates in the one-dimensional thermal response model of Section 2.3 are macro scale. In this work, the macro-ablation recession rate of the material v macro is defined as the average ablation rate, and can be written as: where δ(t) and δ 0 are the sample thickness at time t and the initial thickness of specimen (m), respectively, and t is the ablation time (s). The micro-ablation recession rate v micro can be written as: , when the top of the fiber morphology presents a platform v i , when the top of the fiber morphology is needle-shaped . (11) The reaction to produce carbon monoxide can be viewed as a one-step reaction. Thus, the temperature-dependent reactivity k can be described by the first-order Arrhenius equation [44]: where A is the pre-exponential factor (m·s −1 ), Ea is the activation energy of the reaction between carbon and oxygen generating carbon monoxide (J·mol −1 ), R is the universal gas constant (J·mol −1 ·K −1 ), and T is the temperature (K). We assume that the oxygen transport can be described by pure diffusion. Thus, the mass conservation for the gas can be written as [45]: where D is the diffusion coefficient (m 2 ·s −1 ). The mutual diffusion of oxygen and carbon monoxide is given by [46]: where n is the number density (m −3 ), d 1 and d 2 are the molecular diameters of O 2 and CO (m), respectively, and m 1 and m 2 are the molar masses of O 2 and CO (kg·mol −1 ), respectively. It is assumed that the material quickly reaches a steady state at the microscopic scale, and that the oxygen concentration has a gradient only along the Z-direction. Hence, according to Equation (13), the oxygen concentration distribution C(z) can be written as: where δ c is the boundary layer thickness (m), C 0 is the molar concentration of the boundary layer (mol·m −3 ), and C i = C 0 /(1 + Da i ), where Da i = k i (h f + δ c )/D is the Damköhler number [47]. The top of the fiber is flat at first. Therefore, the fiber height h f can be written as: When the platform at the top of the fiber disappears, Equation (16) no longer works because the ablation recession velocity of the fiber is not parallel to the Z-axis. As shown in Figure 3, in the steady-state ablation process, the relationship between the ablation velocity of the inclined fiber and the interphase can be rewritten as: where α is the tilt angle (°), zα and lα are the coordinates of the fiber surface (m), χ = D/ki is the characteristic length of diffusion (m) (which is related to the interphase reactivity and the gas diffusion), and η = (kiΩi)/(kfΩf) is a dimensionless number that represents the relative reactivity. When the tilt angle α = 90°, hf is given by the steady-state characteristic height: where lf is the characteristic radius of the carbon fiber on the cross-section (m).

One-Dimensional Thermal Response Model
In the thermal response model, we have made appropriate simplifications. There are the following assumptions: The fiber surface morphology can then be described as: where α is the tilt angle ( • ), z α and l α are the coordinates of the fiber surface (m), χ = D/k i is the characteristic length of diffusion (m) (which is related to the interphase reactivity and the gas diffusion), and η = (k i Ω i )/(k f Ω f ) is a dimensionless number that represents the relative reactivity. When the tilt angle α = 90 • , h f is given by the steady-state characteristic height: where l f is the characteristic radius of the carbon fiber on the cross-section (m).

One-Dimensional Thermal Response Model
In the thermal response model, we have made appropriate simplifications. There are the following assumptions: (a) The material surface is subjected to uniform heat flow, and the other surfaces satisfy the adiabatic boundary conditions.
(b) The temperature changes only in the sample thickness direction, that is, it satisfies the one-dimensional thermal response model.
(c) The residual resin pyrolysis in carbon/carbon composites can be neglected.
In the one-dimensional heat transfer problem, the energy conservation equation in the ablation material control body can be written as [48,49]: where, ρ is the density, C p is the specific heat capacity, and λ is thermal conductivity coefficient.
Equation (20) describes the change of internal energy due to heat conduction. The sample is spatially divided into M + 1 elements, as shown in Figure 4a. Element 1 is the upper surface element of the sample, which includes the endothermic heat capacity Q p , energy of heat-flow heating Q w , thermal radiation Q r , chemical reaction heat Q c , and heat conduction Q cond . As shown in Figure 4b, the energy conservation equation for element 1 is given by [50]: where where q w , q cold , h w , and h re are the actual heat flux considering the thermal blockade effect, the cold wall heat flux, the wall enthalpy, and the recovery enthalpy of the boundary layer, respectively; ψ, M CO , . m CO , M C , ∆H C , ε, σ, T w and S are the mass ejector factor, molar mass of CO, mass flux of CO, molar mass of C, reaction heat, radiation coefficient, Stefan-Boltzmann constant, wall temperature, and the area of the sample section, respectively. Thus, the energy conservation equation for element 1 can also be written as: There are only endothermic heat capacity and the heat conduction in the internal elements of the material. As shown in Figure 4c, the energy conservation equation for elements 2 to M is given by Because of the adiabatic boundary condition, the heat transfer between element M + 1 and the outside environment is ignored. As shown in Figure 4d, the energy conservation equation for element M + 1 can be expressed as: 1 and the outside environment is ignored. As shown in Figure 4d, the energy conservation equation for element M + 1 can be expressed as:

Moving Boundary
As the surface temperature of the carbon/carbon composite increases continuously during the ablation process, the oxidation of carbon on the surface leads to a reduction of the ablation material thickness. The reduction in material thickness causes the change in the heat transfer. Consequently, it is necessary to consider the boundary movement caused by ablation in the one-dimensional thermal response model. The thickness of the composite and elements at time t can be written as: The temperature of each element considering the boundary movement is corrected by the interpolation method.

Moving Boundary
As the surface temperature of the carbon/carbon composite increases continuously during the ablation process, the oxidation of carbon on the surface leads to a reduction of the ablation material thickness. The reduction in material thickness causes the change in the heat transfer. Consequently, it is necessary to consider the boundary movement caused by ablation in the one-dimensional thermal response model. The thickness of the composite and elements at time t can be written as: The temperature of each element considering the boundary movement is corrected by the interpolation method.

Calculation Process
The heat transfer and surface recession should be considered when solving for the temperature distribution in the carbon/carbon composite. Figure 5 shows the flow chart for the solution of the carbon/carbon composite ablation model using MATLAB. The first step is to discretize the spatial and time domains, as described in Section 2.3. Then, based on the initial boundary conditions, the actual heat flux on the material surface q w is calculated. In the third step, the temperature distribution at time t is solved using the energy conservation equations in the one-dimensional thermal response model. Because the energy equations are written as a block tridiagonal matrix, the chasing method is used to solve this problem. Subsequently, the surface temperature can be obtained. The fourth step is to compute the micro-ablation recession rate v micro , the ablation recession thickness of the material δ ac , and the macro-ablation recession rate of the material v macro . The parameters of each element are then updated in the fifth step. If the cumulative time has not reached the set requirements, steps 2-5 are repeated and t is set to t + ∆t; otherwise, the results are output.
on the initial boundary conditions, the actual heat flux on the material surface qw is calcu-lated. In the third step, the temperature distribution at time t is solved using the energy conservation equations in the one-dimensional thermal response model. Because the energy equations are written as a block tridiagonal matrix, the chasing method is used to solve this problem. Subsequently, the surface temperature can be obtained. The fourth step is to compute the micro-ablation recession rate vmicro, the ablation recession thickness of the material δac, and the macro-ablation recession rate of the material vmacro. The parameters of each element are then updated in the fifth step. If the cumulative time has not reached the set requirements, steps 2-5 are repeated and t is set to t + Δt; otherwise, the results are output.  Table 1 displays the relevant parameters of the carbon/carbon composite. We additionally assume that the boundary layer thickness at the stagnation point δc is 0 m [40]. Using the above solution procedure, the material surface and back surface temperatures, the interphase micro-ablation recession rate, and the macro-ablation average recession of carbon/carbon composites (Φ 24 mm × 100 mm) ablated for 30 s were calculated. The heat flux qcold, gas recovery enthalpy at the boundary layer hre, and stagnation point pressure Ps were respectively set to 18,338 kW·m −2 , 8051 kJ·kg −1 , and 0.598 MPa for the calculation. The results and discussion are shown in Section 4.  Table 1 displays the relevant parameters of the carbon/carbon composite. We additionally assume that the boundary layer thickness at the stagnation point δ c is 0 m [40]. Using the above solution procedure, the material surface and back surface temperatures, the interphase micro-ablation recession rate, and the macro-ablation average recession of carbon/carbon composites (Φ 24 mm × 100 mm) ablated for 30 s were calculated. The heat flux q cold , gas recovery enthalpy at the boundary layer h re , and stagnation point pressure P s were respectively set to 18,338 kW·m −2 , 8051 kJ·kg −1 , and 0.598 MPa for the calculation. The results and discussion are shown in Section 4.

Parameters Values
Pre-exponential factor of fiber, A f (m·s −1 ) [51] 2002.2 × T Pre-exponential factor of interphase, A i (m·s −1 ) 241,000 Activation energy of fiber, Ea f (J·mol −1 ) [51] 353,365 Activation energy of interphase, Ea i (J·mol −1 ) [44] 245,000 Density of fiber and interphase, ρ f and ρ i (kg·m −3 ) [14] 2000 Thermal conductivity coefficient, λ (W·m −1 ·K −1 ) [52] 151.55 Specific heat capacity, C p (kJ·kg −1 ·K −1 ) [25] 0.997 + 2.93 × 10 −5 × T Radiation coefficient, ε (-) [53] 0.9 Table 2 and Figure 6 show a comparison between the calculated and experimental results for the ablation performance under an oxyacetylene flame. The macro-ablation recession rate and surface temperature predicted by the current model are consistent with the experimental results in Reference [25]. As shown in Figure 6a, the surface temperature of the material increased rapidly to approximately 2000 K under the action of the heat flux. The surface temperature and back surface temperature of the material after ablation for 30 s were 2967 K and 976 K, respectively. As shown in Figure 6b, the micro-and macro-ablation recession rate of the material increased slowly during the first 6 s of ablation because of the low surface temperature and transient-state ablation. After ablation for 6 s, the micro-and macro-ablation recession rates increased significantly along with the surface temperature. This is because the surface temperature is still low during the first 6 s of ablation, and the ablation morphology of the fiber has not yet reached its steady state during the first 6 s of ablation and shows a plateau. The micro-and macro-ablation behaviors are thus described by the recession rate of the fiber, which is relatively small. The final macro-ablation recession rate of the material reached after ablation for 30 s was 0.2003 mm/s. The prediction ablation recession rate of this work is greater than that of the traditional model, which is approximately 6.5%. The reason is that the ablation velocity direction of the fiber is not parallel to the macroscopic recession direction, and this factor is considered in our model. The rationality of this result can be explained in Figure 3. Therefore, it is necessary to consider the characteristics of fiber morphology to modify the thermochemical ablation model.

Factors Affecting Micro-Ablation Recession Rate
The surface morphology reaches a steady state on the microscopic scale when the ablation surface temperature is given. When the carbon fibers are not platform-shaped, it is difficult to accurately find a recession velocity parallel to the object surface on the fiber surface. Thus, the micro-ablation recession rate of the carbon fibers cannot be directly used to describe macro-ablation recession. However, due to the interphase in the model assumed to be parallel to the object plane, the micro-ablation recession rate of the interphase is a crucial parameter that directly reflects the magnitude of the macro-ablation average recession rate. Thus, the influence of the fiber orientation, fiber radius, interphase density, and ablation surface temperature on the micro-ablation recession rate of the interphase was studied, and the results are shown in Figure 7. The results indicate that the interphase

Factors Affecting Micro-Ablation Recession Rate
The surface morphology reaches a steady state on the microscopic scale when the ablation surface temperature is given. When the carbon fibers are not platform-shaped, it is difficult to accurately find a recession velocity parallel to the object surface on the fiber surface. Thus, the micro-ablation recession rate of the carbon fibers cannot be directly used to describe macro-ablation recession. However, due to the interphase in the model assumed to be parallel to the object plane, the micro-ablation recession rate of the interphase is a crucial parameter that directly reflects the magnitude of the macro-ablation average recession rate. Thus, the influence of the fiber orientation, fiber radius, interphase density, and ablation surface temperature on the micro-ablation recession rate of the interphase was studied, and the results are shown in Figure 7. The results indicate that the interphase micro-ablation recession rate changes slowly with an increase in the fiber orientation, radius, and interphase density when the temperature is lower than 2000 K. The material reactivity at low temperatures is lower than that at high temperatures, and the change in the micro-ablation recession rate induced by the change in the influencing factors is not significant. An increase in the interphase density indicates a greater number of carbon atoms per unit volume. If the reactivity on the interphase surface is constant, it takes longer for the dense interphase to decay to the same thickness. Thus, the micro-ablation recession rate decreases with an increase in the interphase density, as shown in Figure 7d. However, the decrease in the micro-ablation recession rate of the interphase causes the steady-state morphology of the fiber to decrease. The oxygen concentration on the surface of the interphase increases, resulting in an increase in the micro-ablation recession rate. This reduces the decrease in the micro-ablation rate.
In summary, the anti-ablation efficiency of the material can be enhanced by increasing the fiber radius and the interphase density.

Factors Affecting Macro-Ablation Recession Rate
The macro-ablation average recession rate is the most important indicator used in engineering to evaluate the properties of ablation materials, which are affected by the heat flux, pressure, radiation coefficient, specific heat capacity, fiber radius, fiber orientation, As shown in Figure 7a, the fiber orientation does not have a significant influence on the micro-ablation recession rate at different temperatures. The micro-ablation recession rate was reduced by only approximately 4% over the fiber orientation range studied, as shown in Figure 7b. This is because the fiber height difference induced by the fiber orientation is not significant, as has been reported in Reference [40]. However, 4% optimization is also very important and cannot be ignored in the aerospace field.
As shown in Figure 7c, the micro-ablation recession rate is very sensitive to changes in the fiber radius, especially when the fiber radius is small and the temperature is high. The fiber height in the steady state increases with the fiber radius. This causes the oxygen concentration on the interphase surface to decrease. Thus, the micro-ablation recession rate decreases significantly as the fiber radius increases. At 3000 K, the micro-ablation recession rate of the ablation material without carbon fibers is much greater than of the material with carbon fibers. The results also indicate that the interphase micro-ablation recession rate does not change with an increase in the fiber radius, when the temperature is 3000 K and the fiber radius is larger than 22.5 µm. The reason is that the steady-state ablation morphology of the fiber presents a platform shape, and the fiber height no longer increases with an increase in the fiber radius in this case.
An increase in the interphase density indicates a greater number of carbon atoms per unit volume. If the reactivity on the interphase surface is constant, it takes longer for the dense interphase to decay to the same thickness. Thus, the micro-ablation recession rate decreases with an increase in the interphase density, as shown in Figure 7d. However, the decrease in the micro-ablation recession rate of the interphase causes the steady-state morphology of the fiber to decrease. The oxygen concentration on the surface of the interphase increases, resulting in an increase in the micro-ablation recession rate. This reduces the decrease in the micro-ablation rate.
In summary, the anti-ablation efficiency of the material can be enhanced by increasing the fiber radius and the interphase density.

Factors Affecting Macro-Ablation Recession Rate
The macro-ablation average recession rate is the most important indicator used in engineering to evaluate the properties of ablation materials, which are affected by the heat flux, pressure, radiation coefficient, specific heat capacity, fiber radius, fiber orientation, interphase density, and thermal conductivity of the material. Thus, the factors influencing the macro-ablation average recession rate are investigated in this section, and the results are shown in Figure 8.
The relationship between the macro-ablation average recession rate and the fiber orientation is shown in Figure 8a, in which the change in the thermal conductivity coefficient is not considered. The fiber orientation has little effect on the steady-state morphology of the fiber. Thus, the macro-ablation average recession rate is not significantly changed by the variation of the fiber orientation, similar to the micro-ablation recession rate. However, this result is inconsistent with the previous experimental results. Figure 8b shows the results of considering the change of thermal conductivity coefficients with fiber orientation. The results indicate that the fiber orientation has a significant influence on the macroablation average recession rate when the ratio of the radial and axial thermal conductivity coefficients λ 1 /λ 2 is considered. The results explain the reason why the fiber orientation affects the ablation performance. The key reason for this is because the fiber orientation affects the thermal conductivity of the material along the thickness direction, which in turn modifies the material surface temperature and the reactivity. The lowest macroablation average recession rate occurs in the direction of the maximum thermal conductivity coefficient. However, increasing the fiber orientation angle increases the temperature of the back surface. Hence, a proper arrangement of fibers can effectively improve the ablation resistance of the material, albeit at the cost of an increased back surface temperature.
The material and environment parameters in Figure 8c are the same as those in Section 3, except for the heat flux and fiber radius. The results indicate that the average macro-ablation recession rate increases rapidly with the heat flux. The reason is that a larger heat flux heats the surface of the material to a higher temperature faster and results in higher reactivity.  As shown in Figure 8d, the effect of pressure on the macro-ablation average recession rate is similar to that of the heat flux. The pressure affects the ablation rate by changing the oxygen concentration and the diffusion rate. A higher pressure increases the oxygen concentration in the environment but reduces the diffusion coefficient.
As shown in Figure 8e-g, the macro-ablation average recession rate decreases steadily with the increase of the radiation coefficient, specific heat capacity, interphase density, and thermal conductivity coefficient. Hence, the ablation performance of the material can be improved by increasing its radiation coefficient, specific heat capacity, interphase density, or thermal conductivity coefficient. The parameters that influence the macro-ablation average recession rate within the context of parameter changes discussed are, in deceasing order of significance, the specific heat capacity, thermal conductivity coefficient, radiation coefficient, and interphase density. The mechanism through which an increase in the thermal conductivity coefficient, specific heat capacity, or radiation coefficient reduces the ablation rate is mainly the reduction of the surface temperature. Therefore, the reactivity of each component is reduced. The effect of the interphase density on the macro-ablation average recession rate is similar to that of the micro-ablation recession rate, which has been discussed in Section 4.1.
Meanwhile, the influence of the fiber radius on the macro-ablation average recession rate is significant, as shown in Figure 8c-g. For example, the macro-ablation average recession rate of the ablation material without carbon fibers is approximately 1.57 times greater than that of the same material with 30 µm radius carbon fibers under high heat flux (25,000 kW·m −2 ), and it is approximately 4.71 times greater than that of the same material with 30 µm radius carbon fibers under low heat flux (10,000 kW·m −2 ). In addition, it is worth pointing out that when the fiber radius is larger than 20 µm, the impact of the fiber radius on the macro-ablation average recession rate is not significant. This is because although increasing the fiber radius still increases the fiber height and decreases the oxygen concentration on the interphase surface when the fiber radius is greater than 20 µm, the influence of the increased radius is greatly diminished. This implies that it is difficult to efficiently increase the ablation performance of a material by increasing the fiber radius when the fiber radius is greater than 20 µm.
In addition, the relationship between the macro-ablation recession rate of cyclic ablation and time is shown in Figure 9. The results indicate that the macro-ablation recession rate of ablation for 30 s increases rapidly with the increase of the number of cycles. The macro-ablation recession rates of cycle 1, cycle 2, and cycle 3 are 0.200, 0.224, and 0.233 mm/s respectively. The reason is that the surface of the sample after ablation is rougher than that of the original material, which increases the oxidation rate. Meanwhile, the thickness of the sample after ablation decreases, resulting in a higher surface temperature and an increased oxidation rate. rate of ablation for 30 s increases rapidly with the increase of the number of cycles. The macro-ablation recession rates of cycle 1, cycle 2, and cycle 3 are 0.200, 0.224, and 0.233 mm/s respectively. The reason is that the surface of the sample after ablation is rougher than that of the original material, which increases the oxidation rate. Meanwhile, the thickness of the sample after ablation decreases, resulting in a higher surface temperature and an increased oxidation rate.

Conclusions
A thermochemical ablation model of microstructure-based carbon/carbon composites was established to investigate the effects of the fiber radius, fiber orientation, inter-

Conclusions
A thermochemical ablation model of microstructure-based carbon/carbon composites was established to investigate the effects of the fiber radius, fiber orientation, interphase density, surface temperature, radiation coefficient, specific heat capacity, and thermal conductivity coefficient on the ablation recession rate. The predicted results from this model are consistent with the experimental results. The differences in the thermal conductivity coefficient caused by the fiber orientation can significantly affect the ablation performance of the material. Thus, a suitable arrangement of fibers can easily improve the anti-ablation ability of the material. The ablation recession rate of the ablation material without carbon fibers is much greater than that of the material containing carbon fibers. The anti-ablation efficiency of the material can be improved by increasing the radiation coefficient, specific heat capacity, interphase density, and thermal conductivity coefficient of the material. The parameters that influence the macro-ablation average recession rate are, in deceasing order of significance, the specific heat capacity, thermal conductivity coefficient, radiation coefficient, and interphase density. The influence of the fiber radius on the ablation recession rate is not significant when the fiber radius is sufficiently large (above 20 µm). The ablation resistance decreases with increases in the number of ablation cycles. The thermochemical ablation model in this work offers a guide for the design and preparation of better ablation materials.