Abstract
A mathematical model based on minimal thermal resistance and equal law of specific equivalent thermal conductivity is developed to discuss the heat transfer characteristics of ablative thermal insulating material from the mesoscopic scale. Based on the statistical results of mesoscopic parameters, the microstructure unit cell model was established to analyze the influence rule of mesoscopic parameterization which includes the size, distribution, and positional relation of microsphere and fiber. The results show that the equivalent thermal conductivity decreases with the density, size, distribution area, and distance of microsphere and the space distance and volume fraction of fiber decreasing. Besides, the equivalent thermal conductivity will become larger when more quality of heat transfers along the fiber direction. Exploring the relationship between the macroscopic heat transfer process and the microstructure is meaningful for exploring the heat transfer behavior of thermal insulating material and improvement of the processing technology.
1. Introduction
With the acceleration of the flight speed of an aircraft, and the strong aerodynamic heating of a reentry aircraft, to ensure that the internal instruments are not burnt out, it is necessary to take effective measures of heat insulation. Ablative thermal protection is a widely used thermal protection method. The ablative thermal insulating material will change a lot with the processing of ablation; it includes the pyrolysis, carbonization, chemical reaction, surface ablation, and second reaction of pyrolysis gas. The process is very complicated, and the material property will change. The microstructure of the ablative thermal insulating material will also change in the process of ablation. For example, the porosity increases, matrix deposits, fiber orientation changes, and distribution of microsphere and fiber will also change. The heat transfer process is very complicated. Solid and gas conduction, gas convection, and radiation exist in the process of ablation. However, according to Skochdopole’s research [1], when the diameter of a hole is less than four millimeters, the gas convection can be ignored. As the radiation is very complicated, here, we just consider the influence of microstructure on the solid and gas conduction mechanism at room temperature; this can help us better understand the influence of the microstructure parameter on the heat transfer mechanism. Combining with the microstructure changes during the process of ablation, we can better understand the relationship between the microstructure parameter and the heat transfer process, and it is very meaningful for the material process improvement.
Theoretical models of evaluating the thermal conductivity of composites mainly include effective medium approximation (EMA) and micromechanics method. The EMA method comes from the classical Maxwell model [2–4]; it is used in the electromagnetic field at the earliest. Then, several researches used it to get the thermal conductivity [5]. Fricke [6] extended the Maxwell model to ellipsoidal particles. Hasselman and Johnson [7] considered the effect of interface and particle size. Bruggeman [8–10] considered the thermal interaction to extend the Maxwell model to make it suitable for a variety of materials and high filler volume fractions. Every extended Bruggeman model considers the interfacial thermal resistance between particles and matrix [11, 12]. The micromechanics methods are based on the variational principle and the mean field approximation, for example, the Hashin-Shtrikman bounding model [13], Mori-Tanaka model [14], and Benveniste’s model [15]. Besides, there are also some semiempirical models [16–18]. The Maxwell model and Fricke model can only be used in the lower volume fraction of particles which guarantees that the particles are apart from each other. The Hasselman-Johnson model, Bruggeman model, and every model considered the thermal resistance between fillers, but it is difficult to determine thermal resistance parameters, and sometimes, they should be dependent on experiment. The Hashin-Shtrikman bounding model is based on the variational principle; the Mori-Tanaka model and Benveniste’s model are based on mean field approximation. They focus on one factor and are suitable for some specific material. Semiempirical models are dependent on experiment results, but they are often hard to realize.
The main theoretical models mostly focus on one factor about some kinds of material, and they have some limitation more or less. There is another theoretical model rooting from the electrical field. It is based on the law of minimal thermal resistance and the equal law of the specific equivalent thermal conductivity. Its main idea is to treat each section of microstructure as the thermal resistance. By analyzing the heat transfer mechanism, the series-parallel thermal resistance model is established; then, the equivalent thermal conductivity can be gotten. The challenge of this method is how to establish the unit cell model and the series-parallel thermal resistance model. But this method is convenient to consider the influence of microstructure on heat transfer characteristics. Literature reported so far gives some application of this model. For example, Zeng et al. [19], Cheng et al. [20], Li [21], Liang and Liu [22], Agrawal and Satapathy [23], and Devireddy and Biswas [24] used this model to predict the equivalent thermal conductivity.
The above theoretical models can all forecast the thermal conductivity of the composite and conform to the experiment values, but for the influence of mesoscopic parameters on the thermal conductivity of composites, there is less discussion. In the previous article of the authors [25], the mathematical and finite element models were used to analyze mesoscale heat transfer characteristics, and some of the mesoscopic parameterization was analyzed; however, it is not very detailed and full-scale. This paper discusses the heat transfer characteristics of thermal insulation materials with the establishment of a theoretical model, especially to explore the influence law of mesoscopic parameters from different scales. The purpose is to explore the relationship between material macroscopic heat transfer mechanism and material mesoscopic parameters and to provide reference for process design of ablative thermal insulating material.
2. Mesoscopic Observation and Statistical Analysis
Figure 1 shows the microstructure image of a kind of ablative thermal insulating material which is composed of hollow microspheres, resin matrix, fibers, and irregular pores. Statistics of the sizes of microspheres and fibers obtained through the Image-Pro Plus software are shown in Figure 2. Sizes of microspheres and fibers conform to the Gaussian distribution. The size range of microspheres and fibers can be obtained through the statistics results, and it is useful for the parameterization in the following passage that it determines the range of variables.
(a)
(b)
(c)
(a) Microsphere radius distribution
(b) Fiber length distribution
(c) Fiber diameter distribution
3. Mathematical Model
3.1. Mathematical Model of Microsphere Scale
The heat transfer model at the microsphere scale is established based on the law of minimal thermal resistance and the equal law of the specific equivalent thermal conductivity, as shown in Figure 3 [21]. In this article, this model is extended to be used for the heat transfer characteristics of thermal insulating material with complicated microstructure and the parameterization influence analysis.
Here, the unit cell model is the cube whose side length is . The outer radius and inner radius of the microsphere are separately and .
We consider the heat transfer of the unit cell in the direction. There are three main ways of heat transfer: (1) heat conduction through solid and gas, (2) natural convection through the gas in the microsphere and the pores between the microspheres, and (3) thermal radiation passing through the inner and outer surfaces of the microsphere. Since the thermal conductivity of materials under normal temperature is researched, the amount of thermal radiation is very small, so it is ignored here. At the same time, the size of microspheres and the pores between microspheres is also small, so convection heat transfer is not considered in the modeling process. The following basic assumptions are established: (1) the microsphere and the matrix are in good contact; that is, the role of additional pores is not considered; (2) the microspheres are assumed to be periodically distributed; and (3) the microsphere scale is the micron scale, and the influence of the thin transition zone on the microsphere heat transfer is not considered. In this way, the equivalent thermal conductivity of the unit cell of a single microsphere can be predicted according to the statistical law [21].
The macroscopic equivalent thermal conductivity of materials is derived based on the heat resistance series and parallel law. As shown in Figure 3, the unit cell is divided into two parts. One part is the part without microspheres, and the height is ; the other part is the part containing microspheres, and the height is .
Considering the integral along the direction, the unit cell is divided into thin slices whose height is . For the first part,
For the second part, it is obtained by the Fourier law: where is the area of the cross-section perpendicular to the heat transfer direction, is the quantity of conduction heat of the matrix in the second part, is the quantity of conduction heat of the air inside the microsphere, and is the quantity of conduction heat of the microsphere wall. The expressions of each quantity of conduction heat are as follows:
Then,
The average thermal conductivity of each part is as follows:
The thermal resistance of the two parts are
Since the heat is transferred along the direction and it is in a series structure, the total thermal resistance and the equivalent thermal conductivity of the composite are separately
Then,
The volume fraction of microsphere is
The relationship between the radius inside and outside of microspheres is
Besides,
Then,
The density of air is relatively small compared with the other components, so it can be ignored. Therefore, equation (12) can be simplified as
The above is the process deriving the equivalent thermal coefficient of thermal conductivity of the microsphere scale based on the law of minimal thermal resistance and the equal law of the specific equivalent thermal conductivity. At the fiber scale, we use the same method to get the macroscale thermal coefficient of thermal conductivity. The results of this theoretical method are verified by experimental results in literature [25]; thus, the correctness of this theory is certified. So, the method is used for the following discussion.
3.2. Heat Transfer Model of Matrix and Porosity
It can be seen from Figure 1 that there are four mesoscopic structures, namely, microspheres, matrix, fibers, and pores. At the microsphere scale, the heat transfer between matrix and pores is rather complicated. Based on the observation results of mesoscopic structure images and the actual possible heat transfer situation, several heat transfer models of matrix and pores are established.
3.2.1. Parallel Model of Microsphere and Matrix
It can be seen from Figure 4 that in the process of heat transfer, most of the quantity of heat is transferred along the microsphere wall and the matrix outside the microsphere, so it is approximate to a parallel heat transfer structure. Here, a parallel heat transfer model between the matrix and the pores can be adopted [26], and the expression is
In equation (14), is the thermal conductivity of the new equivalent matrix obtained in the parallel model between the pores and the matrix, is thermal conductivity of the resin matrix, is the thermal conductivity of air, is the value of porosity, and is the volume fraction of the resin matrix.
3.2.2. Spherical Pore Inclusion
The actual pores mostly exist among the microspheres, and they are irregular-shaped; here, they are assumed to be spherical pore inclusions. The established unit cell model is shown in Figure 5. The side length of the unit cell is , and the radius of the spherical pore is ; the equivalent thermal conductivity is derived by the series-parallel method.
The equivalent thermal conductivity of spherical pore inclusion is
Then,
3.2.3. Cubic Pore Inclusion
It is assumed that the pores are cubes, as shown in Figure 6. The side length of the unit cell is , and the side length of the cubic pore is .
The expression of thermal conductivity is the same as formula (15), where
Substituting into formula (15), we can get
The results of the above three models are listed in Table 1. A, B, and C are three kinds of ablative thermal insulating material.
It can be seen from Table 1 that the results of the parallel model and the spherical pore inclusion model are relatively close, while the results of the cubic pore inclusion model are somewhat different from the former two models. It is because the actual pore structure is irregular, and the angular shape of the cubic pore is relatively clear, so it is not quite consistent with the actual situation. Among the three models, the parallel model is relatively simple and more consistent with the real heat transfer mechanism. Therefore, the parallel model is used to represent the heat transfer mechanism between the pores and the matrix.
4. Parametric Analysis
It can be seen from the microstructure images of ablative thermal insulating materials in Figure 1 and the statistical results in Figure 2 that the length scale of fibers is much larger than that of microspheres and pores, so the analysis method of two scales can be adopted here. From the microsphere scale, the size and distribution of microspheres are random, but from the fiber scale, the influence of microsphere randomness is weak as the fiber size is much larger compared with the microsphere size. Besides, investigating the influence of microspheres and fiber parameters separately can analyze the influence of different parameters more accurately.
4.1. Parameterization Influence Rule of Microsphere Scale
Due to the different position relationships between microspheres, the size of microspheres, and the changes in distribution, the transfer of heat flow will be affected, which will affect the transfer process of heat. In this way, the equivalent thermal conductivity of the calculated unit cell will be different.
4.1.1. The Density of Microspheres
In the actual material production process, the density of microspheres cannot be obtained accurately, and due to the influence of material technology, the wall thickness of each microsphere will change, so its density is a variable quantity, generally ranging from 200 to 400 kg/m3. The following discussion is about the change of the equivalent thermal conductivity when the density of microsphere changes. When the density of microspheres changes, the volume content of each component also changes, but the mass content of each component and the density of other components remain unchanged. Taking these factors into consideration, the curve of the change of equivalent thermal conductivity with the density of microspheres is shown in Figure 7.
It can be seen from Figure 7 that the equivalent coefficient of thermal conductivity increases with the increasing density of microsphere. It is because the microsphere is a hollow structure; larger density means larger volume of the microsphere wall and less gas inside the microspheres. So, the quantity of heat can be transferred more through the microsphere wall, and then, it can get larger equivalent thermal conductivity. So, in order to increase heat insulating performance, the density of microspheres should be smaller. However, when the density of microspheres is smaller, the microsphere wall is thinner and the mechanical strength is lower. So, the thermal insulating performance should be improved under the condition of ensuring the mechanical strength.
4.1.2. Microsphere Size
The influence of microsphere size on its heat transfer performance is worth studying. Here we assume that the size of microspheres is the same in the same model. As shown in Figure 8, microspheres are distributed on the plane whose coordinate is zero. The change of the equivalent thermal conductivity along the direction is discussed. The equivalent thermal conductivity is derived below. And is the number of microspheres of each model and the total volume fraction of microspheres remains the same.
The equivalent thermal conductivity of direction is as follows:
The value of changes from 1 to 10, that is, under the condition that the volume fraction of the microspheres remains unchanged, the number of microspheres increases, and the volume of each microsphere decreases. The equivalent thermal conductivity of the unit cell varies with the number of microspheres, as shown in Figure 9.
It can be seen in Figure 9 that the microspheres distribute on the same plane, when the number of microspheres increases on the condition that the total volume of microspheres remains unchanged and the equivalent thermal conductivity of direction decreases. The reason is that when the number of microspheres increased, the size of microsphere diminished. The matrix between microspheres decreases, the hollow microspheres will obstruct heat transfer path, and the thermal conductivity of microspheres is very low, so the thermal conductivity of direction decreases when the size of microspheres decreases.
4.1.3. Microsphere Distribution Range
It can be seen from Figure 8 that the center of the microspheres is randomly distributed in a plane. The equivalent thermal conductivity which is perpendicular to the plane will decrease when the size of the microsphere decreases. If the size and number of the microspheres remain unchanged, the distribution area keeps expanding. In other words, when the distribution range of the microspheres increases from a certain height value (greater than the diameter of the spheres) to the whole height of the unit cell model which is shown in Figure 8, the change of the equivalent thermal conductivity of the unit cell is where is the proportion of the distribution area; equals to 1 means the microspheres distribute in all of the unit cell. The variation curve of the equivalent thermal conductivity of the unit cell model changing with is shown in Figure 10.
It can be seen in Figure 10 that when the number and the size of microspheres remain unchanged, the thermal conductivity along the direction increases with the distribution area increase of microspheres. This is because when the microsphere distribution area is small, its configuration will be relatively dense which blocks the transfer of heat. When the distribution area of microsphere increases, the distribution of the microspheres is relatively dispersed and the heat can be transferred through the matrix between the microspheres, so the equivalent thermal conductivity of the unit cell is relatively large. This effect is particularly pronounced when the distribution increase is fast at first and then increases slowly.
4.1.4. The Distance between Microspheres
The microstructure of ablative thermal insulating material is that the microspheres distribute in the matrix, so the distance between the microspheres varies greatly. The influence of the distance between the microspheres on the thermal conductivity is discussed below. The established model is shown in Figure 11. The influence of the distance in direction on the thermal conductivity is considered. It is assumed that the distance between the outer wall of microspheres along the direction is . Obviously, from the derivation process of thermal conductivity, it can be seen that when the distance of direction changes, the thermal conductivity along the and directions remains unchanged. The thermal conductivity only depends on the total volume fraction of the two microspheres. The influence of microsphere distance on thermal conductivity is discussed below.
It is supposed that the outer radius of the microsphere is , the inner radius of the microsphere is , the volume fraction of the microsphere is , the unit cell is the cube whose side length is , and the equivalent thermal conductivity along the direction is
The distance value range of two microspheres is , and the length value of the unit cell is μm. When the microsphere is large enough, the two microspheres could not appear as shown in Figure 12, so the microsphere volume fraction is taken as 5% here. So, the value range of is , and thermal conductivity of direction with respect to the distance between microspheres is shown in Figure 12.
It can be seen from Figure 12 that when the distance between two microspheres along direction increases, the equivalent thermal conductivity of the unit cell along direction increases. This result is in accord with what is shown in Figure 10. Here, we take the distance of two microspheres as variables, so we can study the influence between microspheres on the equivalent thermal conductivity more targeted.
4.2. Analysis of Heat Transfer Characteristics at Fiber Scale
The heat transfer characteristics at the microsphere scale have been discussed above. When the unit cell model is established at the fiber scale, the problem is equivalent to the case of fiber inclusion in a kind of homogeneous matrix.
From Figure 1, we can see that most of the fibers are straight and only a few of the fibers are bent regarding this kind of ablative thermal insulating material. This research is based on the microstructure of this kind of ablative thermal insulating material. So, we take the fibers as cylinder, and it is reasonable. For other materials, if most of fibers are bent, the model should be improved.
4.2.1. The Influence of Fiber Size
The shape of the fiber is generally cylindrical, but the fiber with the same volume content can be made into thin and long type or short and thick type.
The influence of the crudeness or fineness of the fiber on its thermal conductivity is discussed through the inclusion of a single fiber; the unit cell model considering the fiber size is shown in Figure 13. (1)When the heat transfer direction is along the direction of the fiber, its thermal conductivity is expressed as (2)When the heat transfer direction is perpendicular to the direction of fiber, its thermal conductivity is expressed as
In the process that the fiber changes from slender to short and thick, the bottom radius changes from small to large, the length of the fiber can be determined based on the fiber volume invariability, and its variation law is shown in Figure 14. The longer the heat transfer path along the direction of the fiber, the greater the thermal conductivity.
4.2.2. Discussion of the Position Relationship between Fibers
(1) Several Cases of Positional Relationships between Fibers. The orientation, distribution, and position relationship of the fibers have great influence on the heat transfer characteristics. The influence of the relationship between two fibers on heat transfer characteristics is discussed below. Figure 15 is the unit cell model of the relationship between two fibers.
The position relationship between two fibers is represented by the distance between the fiber centers in three directions and the angle, that is, . The average distance between two fibers is ; the angle range between two fibers is . The number of the two fibers is 1 and 2. The middle point of fiber 1 is , the middle point of fiber 2 is , and the angle of the two fibers is denoted as .
The distance between the two fibers is expressed as the distance between the midpoints:
The projection of fiber distance on three coordinate axes, that is, the distances in three directions, is, respectively,
The angle between the two fibers is complicated. As shown in Figure 16, the orientation of the two fibers is variable. When the midpoint of the two fibers is determined and the fiber rotates at any angle around the midpoint, the trajectories of the two endpoints of the fibers are spherical, denoted as spherical surface 1 and spherical surface 2.
(a)
(b)
(c)
When the distance between the two midpoints of the two fibers is greater than the length of the fiber, the sphere formed by the two ends of the fibers does not intersect, no matter what the angle is between the two fibers.
But in fact, the unit cell model is small and the fibers are always long enough which can intersect with each other. When the two spherical surfaces locate at the two corner points of the unit cell, the distance between the midpoints of the two fibers reaches the maximum.
Then,
When two fibers have one end in the volume formed by the intersection of two spheres, the two fibers can intersect; otherwise, the two fibers cannot intersect.
As shown in Figure 16(b), two fibers meet at the endpoints. And as shown in Figure 16(c), two fibers intersect with each other. The expression of the space angle between the two fibers is derived as follows.
It is assumed that the coordinate of any point on spherical surface 1 is ; the coordinate of any point on spherical surface 2 is . Then, the spherical equation is
The direction of fiber 1 is
The direction of fiber 2 is
Therefore, the inclined angle between the two fibers is
When the two fiber ends meet, then the intersection of spherical surface 1 and spherical surface 2 forms a circle of space. The equation is derived as follows:
Then,
This equation is the circular ring equation which is obtained by the two spherical surfaces intersecting. When two fibers are connected to the circular ring, the inclined angle is
The parameters of , , and in equation (39) satisfy equation (38). When two fibers translate to the intersect state, its inclined angle is the space angle of two fibers.
In order to intuitively describe influence of fiber distance and angle on thermal conductivity, two aspects were considered. First is the two fiber spatial disjunct, as shown in Figure 16(a), and the second is the two fiber intersect, including the endpoint connect, as shown in Figure 16(b). Figure 16(c) is the case that the rest part of the fiber is connected.
(2) The Relationship between Thermal Conductivity and Distance in Three Directions of Fiber. For two straight lines in space, they can determine a plane when they are shifted to the intersection state. So, two straight lines can only determine a plane which is parallel to the two straight lines. We can assume the two fibers parallel to the plane to establish a space rectangular coordinate system, as shown in Figure 17.
(1) The Distance between the Fibers Varies in the z Direction. First, the influence of the change of distance in the direction on the thermal conductivity along the direction was studied when and remained unchanged. The distance in the and directions was taken to be the average distance, that is, . Figure 18 is the front view of the spatial distribution model of two fibers.
The distance between the fibers in the direction is denoted as . We assume that the fiber length is the same, denoting as . The fiber is represented as cylinder, and the bottom radius is denoted as ; then, the equivalent thermal conductivity along the direction is
(2) The Distance between the Fibers Varies in the x Direction. The top view of the position relationship between the two fibers is shown in Figure 19.
In Figure 19, the inclined angle of two fibers is ; and are the inclined angles of fiber 1 and fiber 2 in the direction. The following is the discussion about the change of the thermal conductivity of the unit cell along the direction with the distance .
Here, the distance in the direction is taken as the fixed value, that is, . It can be obtained from Figure 19 that the distance range of the two fibers in the direction is
So, the thermal conductivity of the unit cell along direction is
(3) The Distance between the Fibers Varies in the y Direction. As shown in Figure 20, the distance between fibers in the direction is kept constant, which is equal to the average distance between fibers. The change in the equivalent thermal conductivity of the unit cell in the direction is discussed.
The value range of is
So, the equivalent thermal conductivity of the unit cell along direction is
In order to study the relationship between the thermal conductivity of the unit cell in the direction and the distance in the direction, we assume and substitute it into equation (44).
The curves of thermal conductivity changes in , , and directions are shown in Figure 21 [25].
As can be seen from Figure 21, the equivalent thermal conductivity in , , and directions increases with the distance, while the thermal conductivity in and directions is larger than that in direction. When the distance increases to a certain extent, the thermal conductivity in the three directions tends to be the same.
(3) Effect of Fiber Space Angle on Thermal Conductivity. When the distance in the , , and directions of the fibers remains the same, we take , the inclined angle changes. Since the plane of the coordinate system is parallel to the two fibers, the inclined angle of the two fibers is the inclined angle in the plane. When the inclined angle changes, the equivalent thermal conductivity in the direction of the unit cell remains unchanged, while the thermal conductivity in the and directions changes with the inclined angle of the two fibers.
As the equivalent thermal conductivity along the and directions is also related to the angle between the two fibers and the axis, without loss of generality, we can take this .
(1) Thermal Conductivity along x Direction. As shown in Figure 19, the equivalent thermal conductivity along direction is
Since the fiber orientation is stochastic, the fiber inclined angle obeys uniform distribution, and its probability density function is as follows:
Therefore, the equivalent thermal conductivity along direction with probability and statistical significance is
(2) Thermal Conductivity along y Direction. The derivation process of the thermal conductivity in the direction is the same as that in the direction, except that the trig function of the angle is different. The coefficient of thermal conductivity is as follows:
Since the fiber orientation is stochastic, the fiber inclined angle obeys uniform distribution, and its probability density function is as follows:
Therefore, the equivalent thermal conductivity along direction with probability and statistical significance is
It can be seen from Figure 22 that the thermal conductivity is large when the heat transfer direction is consistent with the orientation of the fiber, while the thermal conductivity is small when the orientation of the fiber is perpendicular to the orientation of the heat transfer. This conclusion has been applied in thermal protection materials such as PICA. It is a new type of thermal protection material. The fiber distributes in the plane which is perpendicular to the direction of heat transfer. Fiber orientation with this plane is random, the space distance among fibers is very small, and although there is matrix filling, there are still a lot of pores. The distribution form of the fiber can greatly increase the complexity of the heat transfer path, so the heat transfer is blocked in the process of transmission; then, this microstructure can improve the thermal insulating capability significantly.
4.2.3. Effect of Fiber Thermal Conductivity and Volume Fraction
The effects of the fiber orientation angle, distance between fibers, and angle between fibers on the equivalent thermal conductivity of the unit cell were discussed from the perspective of fiber position relationship. As can be seen from the results, when the above parameters change, the equivalent thermal conductivity of the unit cell changes, but the amount of change is small. The first reason is that the thermal conductivity of the fiber is not particularly large compared with the matrix, and the second reason is that the volume content of the fiber is very small. The following two aspects are discussed to study the change of the equivalent thermal conductivity of the unit cell when the thermal conductivity of the fiber and the volume fraction change. Here, the orientation of a single fiber and the intersection of fibers are considered.
(1) The Orientation of a Single Fiber. Considering that the thermal conductivity of the fiber increases from 0.6 W/(m·K) to 80 W/(m·K), the equivalent thermal conductivity changes with the fiber orientation angle as shown in Figure 23.
As can be seen from Figure 23, when the thermal conductivity of the fiber is small, its orientation angle has little effect on the thermal conductivity of the unit cell. When it increases to more than 20 W/(m·K), its different orientation can significantly affect the equivalent thermal conductivity of the unit cell.
Next, the increase of fiber volume fraction will be discussed. As shown in Figure 24, when the thickness of the unit cell decreases along the direction until it becomes the diameter of the fiber, the change of the equivalent thermal conductivity of the fiber along the direction with the angle between the fiber and the direction is shown in Figure 25. The thickness along the direction is denoted as . The equivalent thermal conductivity of the unit cell along the direction is expressed as follows:
It can be seen from Figure 25 that, when the unit cell is thin along the direction, the fiber volume fraction is large, so the fiber orientation angle has a great influence on the equivalent thermal conductivity of the unit cell. When is 1, it is a cubic unit cell, and the volume fraction of the fiber is small, so the orientation of the fiber has very little effect on the equivalent thermal conductivity of the unit cell.
(2) Fiber Intersection. The fiber intersection model is shown in Figure 26, different thermal conductivities of fibers are considered here, and the equivalent thermal conductivity is shown as equation (52). The equivalent thermal conductivity along the direction varies with the inclined angle, and the thermal conductivity of the fiber is shown in Figure 27:
(a)
(b)
As can be seen from Figure 27, when the thermal conductivity of the fibers is small, the fiber inclined angle has little effect on thermal conductivity of the unit cell. When it increases to more than 20 W/(m·K), the different inclined angles between the fibers can significantly affect the equivalent thermal conductivity.
When the thickness of the unit cell in the direction which is shown in Figure 26 becomes smaller, that is, when the volume fraction of the fiber increases, its equivalent thermal conductivity along the direction can be expressed as
Different proportional coefficients were taken; that is, the fiber volume fraction in Figure 26 was 2.2849%, 1.3405%, 0.8043%, 0.4022%, 0.1341%, and 0.08043%, respectively. The unit cell gradually thickens along the direction, so the fiber volume fraction decreases, and the change of the equivalent thermal conductivity of the unit cell along the direction with the inclined angle is shown in Figure 28.
When the value is small, that is, when the unit cell is thin, the fiber volume fraction is large, so the angle between the fibers has a great influence on the equivalent thermal conductivity. When is 1, it is a cubic unit cell, and the volume fraction of the fiber is small. Therefore, the angle between the fibers has very little effect on the equivalent thermal conductivity.
According to the analysis of the parameterization influence rule, the equivalent thermal conductivity is influenced greatly by the microstructure, and it is meaningful for the actual manufacturing process of ablative thermal insulating material. For example, we can use the hollow microspheres with low density and more and smaller microspheres to let them close to each other to make the equivalent thermal conductivity smaller. For fibers, they should be vertical to the heat transfer direction. In the actual manufacturing, when most fibers distribute randomly in the plane which is perpendicular to the heat transfer direction, the equivalent thermal conductivity will become smaller than other cases. So, the analysis of microstructural parameterization influence rule is meaningful for the manufacturing process improvement.
5. Conclusions
In this paper, the heat transfer characteristics of ablative thermal insulating materials at room temperature are analyzed, and the influence of various parameters on the thermal conductivity is discussed from two scales by establishing a kind of theoretical model. A fast method for predicting the heat transfer characteristics of ablative thermal insulating material was established. The conclusions are as follows: (1)The mathematical model established in this article which is based on the law of minimal thermal resistance and the equal law of the specific equivalent thermal conductivity can get relatively accurate thermal conductivity prediction results. And the mathematical model can also obtain intuitive results of the parameterization influence rule of mesostructure(2)At the microsphere scale, the equivalent thermal conductivity will decrease with the density, size, distribution area, and distance between microspheres decreasing(3)At the fiber scale, the longer the heat transfer path along the direction of the fiber, the greater the equivalent thermal conductivity. As the thermal conductivity of fibers is larger than the matrix and microsphere, the equivalent thermal conductivity will increase with the increasing of fiber volume fraction and fiber thermal conductivity(4)In the actual process of material technology design, the composite material with better performance can be designed considering the influence law of material microstructure parameters. Smaller equivalent thermal conductivity will be obtained if lower density and smaller size, distance, and distribution area of microspheres and fibers which distribute vertical to the heat transfer direction were adopted
Data Availability
(1) The data (table, figure, and curve graph) used to support the findings of this study are included within the article. (2) Previously reported data (table, curve graph, and SEM image) were used to support this study and are available at 10.1155/2019/8142532. These prior studies (and datasets) are cited at relevant places within the text as references [6, 7].
Conflicts of Interest
As authors, we certify that we have NO affiliations with or involvement in any organization or entity with any financial interest (such as honoraria; educational grants; participation in speakers’ bureaus; membership, employment, consultancies, stock ownership, or other equity interest; and expert testimony or patent-licensing arrangements), or non-financial interest (such as personal or professional relationships, affiliations, knowledge or beliefs) in the subject matter or materials discussed in this manuscript.
Acknowledgments
The authors acknowledge the fund of the State Key Basic Research and Development Plan: 2015CB655201.