Next Article in Journal
Synthesis, Characterization, and Catalytic Behaviors in Isoprene Polymerization of Pyridine–Oxazoline-Ligated Cobalt Complexes
Previous Article in Journal
Rheological and Mechanical Characterization of 3D-Printable Solid Propellant Slurry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of the Thermal Diffusivity of Carbon/Phenolic Composites (CPCs) through Oxy-Acetylene Torch (OAT) Test—Part 1: Experimental Characterization and Preliminary Validation

Civil and Environmental Engineering Department, University of Perugia, Strada di Pentima 4, 05100 Terni, Italy
*
Authors to whom correspondence should be addressed.
Polymers 2024, 16(5), 577; https://doi.org/10.3390/polym16050577
Submission received: 7 December 2023 / Revised: 30 January 2024 / Accepted: 16 February 2024 / Published: 20 February 2024
(This article belongs to the Section Polymer Composites and Nanocomposites)

Abstract

:
Carbon/Phenolic Composites (CPCs) are essential to manufacture many portions of the nozzle assembly of Solid Rocket Motors (SRMs) which are essential both to preserve the independent access to space as well as for the homeland security. In our research, a feasible approach aimed at preliminary retrieving the in-plane and out-plane thermal diffusivity of CPCs through the Oxy-Acetylene Torch (OAT) tests was validated. The proposed approach showed to be effective and able to bypass some limitations of common protocols, especially in terms of capability to determine the thermal diffusivity of CPCs at high heating rates. A comprehensive work of comparison of the obtained data with state-of-the-art CPCs such as MX-4926 and FM-5014 has also been carried out, evidencing the effectiveness of the proposed method.

1. Introduction

During last decades, many Thermal Protection System (TPS) materials have been introduced and optimized for different applications. These materials are used to manufacture propulsion systems such as liquid-fueled rocket engines, Solid Rocket Motors (SRMs)—which are essential to manufacture civilian and military launchers—and, more recently, Hybrid Rocket Motors (HRMs) [1]. They are also used to manufacture the heat shields which protect probes and space vehicles during the hypersonic flight (also known as re-entry) through a planetary atmosphere [1]. Refractory Metals (RMs) and Ultra High Temperature Ceramics (UHTCs) have been evaluated or used as TPS materials. Moreover, carbon-based materials such as graphite, which possesses high heat capacity and high energy of vaporization, are also widely used as TPS material. However, even though graphite has a substantially lower cost and density than RMs and UHTCs, it has the drawbacks of possessing relatively poor mechanical properties and exhibiting problems in terms of resistance to thermal shocks. As a result, carbon fibers [2] are typically introduced in a carbon matrix to produce Carbon/Carbon Composites (CCCs) [1]. CCCs are used in the throat region of SRMs which is exposed to the harshest conditions generated by modern aluminized propellants and where the heat fluxes easily exceed ∼1000 W/cm2 [1].
Among TPS materials, Polymeric Ablative Materials (PAMs) possess the widest flexibility in terms of properties tunability [1]. Ablation is a self-regulating heat and mass transfer process in which the incident thermal energy is dissipated through a sacrificial material which undergoes endothermic degradation reactions [1]. Carbon/Phenolic Composites (CPCs) have shown the ability to resist heat fluxes up to 30,000 W/cm2 [1] and are used to produce heat shields of many civilian and military re-entry objects and many parts of SRM nozzle assemblies while CCCs are preferred in the portion of the motor where the erosive phenomena are more severe such as in the throat region. In fact, using a series of carbonization (up to about 1500 °C), graphitization (up to about 2500 °C) and re-impregnations cycles, CPCs can also be converted into CCCs having the desired density [1]. For all these reasons [2,3], CPCs are virtually the most important class of TPS materials.
In addition to the manufacturing challenges related to CPCs and CCCs reported elsewhere [2,3], another critical phase of the development and qualification of a CPC is related to the evaluation of thermo-physical and thermo-mechanical characterization at the proper heating rates and temperatures. In fact, the strength of PAMs such as CPCs strongly depends on temperature, with mechanical properties rapidly decreasing once the temperature exceeds the glass transition temperature of the polymer [4]. Moreover, the decomposition, the thermo-mechanical stresses and the thermo-physical properties of CPCs are strongly dependent on the type of the hyperthermal environment [1] and intimately related to the heating rate [1,5]. Once tested through proper facilities able to recreate in certain extent the hyperthermal conditions, in which the TPS material must operate [1], analytical models have been introduced to predict the response of the given formulation.
The necessity for reliable thermo-physical and thermo-mechanical models was recognized in 1960s/1970s when PAMs became essential in most of the high-temperature aerospace applications [4,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26]. The protocols introduced in these researches identified the methods to measure the decomposition kinetic parameters and the thermal properties [8,9,27,28] required for any ablation modeling [1,29,30] through the use of testing techniques such as Differential Scanning Calorimetry (DSC), Differential Thermal Analysis (DTA), Thermo-Gravimetric Analysis (TGA), Thermo-Mechanical Analysis (TMA), Laser Flash Analyzer (LFA), etc.
For TPS materials such as CPCs it is extremely important to evaluate the thermal diffusivity α defined as follows:
α = k ρ C p
where k is the thermal conductivity, C p is the heat capacity and ρ is the density of the studied material.
For PAMs such as CPCs, the density ρ of the given formulation as a function of temperature is usually determined through TGA test up to temperature of 1500/2000 °C and at heating rate from 10 up to about 500 °C/min [1,8,9,27,28].
The possibility to quantify the irreversible heat capacity [9] of TPS materials such as CPCs in a wide range of temperature is essential to properly evaluate the heat of decomposition of the given material. In a previous paper [31], we summarized the efforts to identify a feasible approach able to assess the heat capacity C p as a function of temperature of TPS materials such as CPCs through the use of common DSC and TGA analysis, taking into account the merits and limitations of common DSC/TGA equipments. Accordingly, the heat capacity was calculated through to the following equation:
C P T = 1 m T H T = 1 m T H t t T = 1 m T P β
where H / T is the heat flow, m T is the mass of the sample rescaled by the TGA in nitrogen, P is the signal of the DSC (W) and β is the heating rate.
However, most of the papers on the determination of the heat capacity C p of PAMs [10,11,12,13,14,15,16,17,18,19,20,21] were based on very old testing facilities [32] and, in many cases, special testing devices were developed ad hoc [1,5,9]. As an example, Lattimer et al. [9] designed a thermal decomposition apparatus that was able to retrieve the heat capacity of composite materials at high temperature. In order to validate the testing protocols, the authors also carried out a comparison with the DSC results.
Finally, the thermal conductivity of a material can be evaluated through different protocols as reported in Table 1 [33].
However, most of these protocols are not suitable to fully assess the determination of the thermal conductivity of a PAM such as a CPC especially due to the fact that it is very difficult to impose high heating rates through common electrical heaters. In fact, CPCs used in SRMs commonly experience heating rates in excess of 1000 °C/min and can reach values of 50,000 °C/min [1]. As a result, the main limitation of testing techniques such as the Laser Flash Analyzer (LFA), which are currently used to determine the thermal conductivity of CPCs, is related to the fact that the tests are carried out under steady-state conditions. On the other hand, during the 1960s/1970s, agencies such as NASA used ad hoc designed transient hot-wire equipment to determine the thermal conductivity of CPCs [34,35]. Similarly, Shilav et al. [33,36] developed a test rig aimed at evaluating the thermo-physical properties of Ethylene Propylene Diene Monomer (EPDM) based Elastomeric Heat Shielding Materials (EHSMs) for SRMs. Even though this device is extremely well designed and accurate, its main limitation is related to the heating rate (of about 10 °C/min).
In general, the methods and experimental devices for the determination of the thermal conductivity and heat capacity of ablative composite materials such as PAMs and CPCs differ from those used for metals. Composites differ from metals in the following ways: (1) they have considerably lower thermal conductivity; (2) there is substantial gas generation in the composite; (3) their thermal properties change non-monotonically; (4) the thermal characteristics of composites depend on rate of heating (or, generally speaking, on ‘history’ of heating) [37].
Moreover, in general, the manufacturing of these apparatuses for the determination of the heat capacity and of the thermal conductivity is not trivial. Certainly, the reports freely available through the NASA Technical Reports Server (NTRS) can help to retrieve some data on International Traffic in Arms Regulations (ITAR) materials (such as the MX-4926 or the FM-5014) and can provide guidance to design special testing apparatus. Nonetheless, also in light of the limitations of most of these protocols, we wanted to validate a simple approach to preliminary evaluate the in-plane and out-plane thermal diffusivity α of PAMs such as CPCs through the Oxy-Acetylene Torch (OAT) tests. From the knowledge of α and of the C p [31], it is possible to retrieve the value of the thermal conductivity k . Using K-type thermocouples to retrieve the in-depth temperatures (as in our experiments), this protocol should be able to measure the thermal diffusivity values from room temperature up to ~1300 °C. This upper limit is also related to the maximum surface temperature that can be typically imposed by an OAT torch to a given ablative material [38,39,40] and (3). Moreover, through the OAT is also possible to reproduce heating rates with values in excess of ~4000 °C/min, enabling the possibility to provide unique insights on the charring phenomena of ablative materials exposed to the harshest conditions [1,41].
In terms of heat transfer analysis and exploitation of the in-depth temperature profiles, the knowledge of the surface thermal boundary conditions in applications such as re-entry flight vehicles or SRM nozzle assemblies is critical to properly design the related TPS as well as for the modeling of the ablation phenomena of charring materials such as CPCs. This is a concern both in the real operational environment of the TPS material as well as in the laboratory test which exploits torches such as the OAT. In fact, harsh ambient conditions related to the given hyperthermal environment [1] or other practical limitations may preclude the use of surface-mounted sensors for measuring the surface temperature and heat flux. To overcome this problem, the temperature profiles obtained from TCs located below the exposed surface can be used to determine the surface thermal boundary conditions. The exploitation of the in-depth temperature data to estimate the surface boundary conditions is known as the Inverse Heat Conduction Problem (IHCP). The diffusion of the thermal front as it penetrates the solid domain—also in the case of charring materials such as CPCs—results in a given temperature profile at the considered depth. Thus, one of the challenges of the IHCP is to exploit the temperature profile at the TC site to estimate the temperature change that occurred at the surface. This implies that small measurement errors in the in-depth TC data are amplified when projected to the surface determining significant errors in the surface prediction. The amplification of measurement errors in the data makes the IHCP prone to errors [1]. Numerous regularization techniques have been developed to obtain a stable prediction of the surface condition. Many methods have been proposed to solve the IHCP. These methods have been summarized in [42].
Some practical applications allow the use of surface-mounted heat-flux gauges. Considerable research has been spent on the development of high-temperature heat-flux gauges that do not require active cooling [42]. However, in many applications, the use of surface-mounted sensors is not practical. In those applications, researchers have proposed the use of plug-type sensors that are instrumented with embedded TCs. For example, plug-type gauges have been used in the Mars Science Laboratory Entry, Decent, and Landing Instrumentation (MEDLI) project [1,43,44]. The MEDLI project used MEDLI Integrated Sensor Plugs (MISPs) that were made of TPS material with four embedded TCs. The measured in-depth TC data were projected in an inverse analysis to obtain an estimation of the surface temperature and heat flux of the MISP. A logical improvement to the use of a plug-type sensor is to perform a calibration test where the sensor is subjected to known surface thermal conditions and the in-depth TC response is measured. This calibration test data, along with in-depth TC data from a new test, can be used to determine the surface thermal conditions that caused the measured in-depth temperature response of the considered test. The Non-Integer System Identification (NISI) method [42] is a system calibration where an impulse function is obtained based on the fractional derivative formulation of the heat equation. The calibration test data are used to determine the calibration coefficients using a least squares method. The surface heat flux of an unknown heating process can be determined using the calibration coefficients along with the in-depth temperature data resulting from the unknown heating process. The Calibration Integral Equation Method (CIEM) is an alternative approach to the concept of system calibration which does not require the calculation of calibration coefficients [42].
Myrick et al. [42] presented an experimental investigation of the accuracy of CIEM approach testing an AISI 304 stainless-steel plug-type test sample heated with a laser power source. The one-probe CIEM requires the knowledge of the surface thermal boundary condition of a calibration test along with temperature data from one in-depth probe in order to predict the surface thermal boundary condition of an unknown (reconstruction) heating process given the temperature data from the same probe. Two modified versions of the CIEM that account for temperature-dependent properties were investigated. This study demonstrated that a calibrated plug-type sensor can be used in high temperature applications to accurately reconstruct the surface heat flux and temperature of an unknown surface heating process. The stainless-steel test sample was a proof of concept for a future prototype of this type of calibrated dual measurement gauge. It is worth noting that in a forthcoming work, we will report the effectiveness of the method summarized in this paper to determine the thermal diffusivity of CPCs on AISI 304 stainless-steel cylindrical plugs.
For several decades, Alifanov spent substantial efforts [45,46,47] on the determination of the thermo-physical properties of TPS materials (mostly non ablative) through the design of ad hoc testing facilities and then exploiting IHCP-based algorithms. As an example, in [47] Alifanov presented a work aimed at the verification of the procedure of determining the thermo-physical properties of non-ablative lightweight ceramic TPS materials [1] under representative experimental conditions through the use of the in-depth temperature profiles. The computational algorithm for the exploitation of the experimental data was based on the solution of the inverse heat transfer problem by the method of iteration regularization.
Zhou et al. [48] adopted an IHCP-based approach in a one-dimensional composite slab with rate-dependent pyrolysis chemical reaction and outgassing flow effects using the iterative regularization approach. The thermal properties of the composites were considered to be temperature dependent (like in a CPC and as in this paper). A non-linear conjugate gradient method formulation was developed and applied to solve the IHCP in an organic composite slab whose front-surface was subjected to high intensity periodic laser heating. It was found that for the low-thermal-conductivity composite materials, in addition to the knowledge of the cold-face temperature of the sample, it was necessary to measure the temperatures at an in-depth location as additional information to recover the exposed surface periodic heating condition.
To conclude, the complexity related to the proper evaluation of the thermo-physical and thermo-mechanical properties of TPS materials, as well as the necessity to have good heat transfer (and ablation) models, have been recently recognized both by the European Space Agency (ESA) and NASA in two space programs. In the Vega C program, ESA announced 2 October 2023 the completion of an independent investigation into an anomaly that took place during a static-fire test of a Zefiro 40 motor 28 June 2023. That test was part of efforts to return the Vega C to flight after a December 2022 launch failure blamed on that motor. The test was intended to confirm the performance of the throat insert, which prime contractor Avio had replaced as part of the recommendations into the December 2022 launch failure. That investigation, released in March, concluded that the CCC from the original supplier (Ukrainian company Yuzhnoye), did not meet specifications. ArianeGroup now supplies the throat insert. “The failure of the test is related to the design of the nozzle that was not upgraded with the change in the CCC supplier for the throat insert”, Giovanni Colangelo, ESA’s inspector general said. The geometry of the new throat insert and its different thermo-mechanical properties contributed to the failure [49]. In the Artemis program—in particular in the Artemis 1 test (December 2022)—NASA managers evidenced that during reentry they noted “more erosion” of the PAM used in the heat shield of the vehicles than expected [50]. For this reason, during a media conference on 9 January 2024, the NASA Administrator announced that Artemis 2, the first crewed mission that will send four astronauts around the moon, had been pushed back from the end of 2024 to no earlier than September 2025. Artemis 3, the first crewed landing, was in turn delayed from late 2025 to no earlier than September 2026. The NASA deputy associate administrator for the Moon to Mars Program said the agency was making good progress on understanding what caused the erosion and expected to find a root cause by spring but needed more time to synthesize the data and update models before flying again.

2. Materials and Methods

2.1. Matrices and Composites

Bakelite PFS121ME, supplied by Hexion (Duisburg, Germany), was selected as a phenolic matrix. This is a liquid phenolic resol with a viscosity of 0.75–1 mPas at 25 °C having methanol as a solvent. The neat CPC formulation referred as C-P was manufactured using the neat PFS121ME. Moreover, a hybrid phenolic-epoxy formulation called P10E was composed by 100 parts of PFS121ME and 10 parts of Epon 828. Epon 828 is an undiluted clear difunctional bisphenol A/epichlorohydrin derived liquid epoxy resin with a viscosity of 100–150 mPas at 25 °C. At the P10E formulation, a stoichiometric amount of hardener TETA (evaluated on the Epon 828 mass) was also added. The resulting CPC was referred as C-P10E. In both cases, a 2 × 2 Twill weave T300-like Polyacrylonitrile (PAN) based on Carbon Fibers (PAN-CFs) fabric (Microtex GG200T1253K with an areal density of 200 g/cm2) was used to produce the laminates. The pre-impregnated fabric was cut in suitable shape for being stacked in a mold: 45 layers of prepregs were overlapped to produce each laminate. The application of pressure and temperature allowed to obtain a square shaped (~100 × 100 mm) composite laminate, approximately 9 mm thick. It is worth remarking that all the composite formulations were consolidated at a very low pressure of 1 bar. The curing cycle used was the following: heating from 30 °C to 60 °C @ 2 °C/min and isothermal at 60 °C for 1 h; heating from 60 °C to 80 °C @ 2 °C/min and isothermal at 80 °C for 1 h; heating from 80 °C to 150 °C @ 2 °C/min and isothermal for 30 min. All the details on the manufacturing, characteristics, and mechanical properties of these two composites are reported in [2].

2.2. Oxy-Acetylene Torch (OAT) Test

From an experimental point of view, the data (in-depth temperature profiles) that will be exploited to determine the thermal diffusivity of the CPCs were based on the temperature profiles acquired with an OAT [1]. Two CPC formulations were considered in the OAT tests and in our study: the C-P and the C-P10E formulations described in previous section.
In general, the OAT test allows to produce high temperatures (∼3000 °C) and heat fluxes of over 800 W/cm2 [51,52,53,54,55,56,57]. Heating rates in excess of 1500 °C/min can also be recreated. In our tests, cubic shaped samples with a size (L) of 10 mm were used. The distance between the sample surface and the flame was set at 15 mm and the heat flux produced by the torch was set at 500 W/cm2 (Figure 1a): the calibration of the power of the torch was carried out using a copper slug calorimeter assuming a diameter of the torch plume of ~4 mm. Once the heat flux is set, the heat load can be imposed selecting the exposure time. A typical heat impulse is showed in Figure 1b.
In our test, an acetylene on oxidizer ratio equal to 1.33 was chosen in order to produce a neutral flame: this condition minimizes the thermo-oxidative ablation of the tested material satisfying the ASTM E-285-80 [1]. The in-depth temperature profiles were acquired at 3 mm (i.e., about 1/3 L, referred as T1) and 6 mm (i.e., about 2/3 L, referred as T2) from the hot surface on the corresponding isothermal surfaces. Fully sheathed, ungrounded, K-type thermocouples were placed in two blind holes which were drilled parallel to the exposed surface of the samples. To ensure a consistent acquisition of the temperatures, the depth of the holes was kept constant at 5 mm, i.e., in correspondence of the centerline of the sample, exactly where the torch impinged. A dual-band pyrometer (DIAS Infrared Systems model PYROSPOT DSR 55N) able to operate up to 3000 °C was used to record the surface temperature of the samples tested with the OAT (Figure 1c). It is worth noting that at a sample-to-pyrometer distance equal to 800 mm (as in our tests), the temperature on the exposed surface of the sample is read on an area (circle) having a diameter of about 3 mm i.e., in line with the size of the torch spot (see Figure 1a).
Finally, the loss of mass due to the ablation process was also evaluated.
Figure 2 reports the appearance of the C-P (a) and C-P10E (b) OAT virgin samples and the testing coupons of the same formulations after the torch test (C-P (c) and C-P10E (d)). For the out-plane tests, an exposure time of 30 s was selected. On the other hand, for the in-plane tests, an exposure time equal to 10 s was chosen due to the fact that the in-depth charring was quite faster than the out-plane case, as a result of the high thermal diffusivity of CPC samples oriented along the fiber direction.
The most representative temperature profiles of the surface temperature T S ( t ) , of the T1(t) (acquired at 3 mm from the surface) and of the T2(t) (acquired at 6 mm from the surface) of the C-P and C-P10E both in the out plane and in plane configurations are reported in Figure 3 and Figure 4.
From a qualitatively point of view, for both the C-P as well as for the C-P10E, considering the samples tested out plane, the approximate heating rate of the rising part of the temperature profile acquired at 3 mm—T1(t)—was of about 1600 °C/min (Figure 3a,b) while the approximate heating rate of the rising part of the temperature profile acquired at 6 mm—T2(t)—was of about 500 °C/min (Figure 3a,b). On the other hand, from a qualitatively point of view, for both the C-P as well as for the C-P10E, considering the samples tested in-plane, the approximate heating rate of the rising part of the temperature profile acquired at 3 mm—T1(t)—was of about 3800 °C/min (Figure 4a,b), while the approximate heating rate of the rising part of the temperature profile acquired at 6 mm—T2(t)—was of about 1900 °C/min (Figure 4a,b). The grey curves reported in Figure 3 and Figure 4 correspond to the surface temperature profiles registered by the pyrometer.

2.3. More Insights on the Heat Capacity of CPCs

In general, the heat capacity of a CPCs can be considered composed by a reversible component (Debye or Einstein like [8]) with a superimposed irreversible part (degradation, i.e., charring of the polymeric matrix). To note the differences between the reversible and irreversible part of the heat capacity, and in order to qualitatively evaluate the degree of conversion of the matrix from an organic material to a carbonaceous medium, a cycle of four DSC tests (through a TA Instruments, model Q200 [31]) was carried out on the P10E matrix used to produce the C-P10E. For the first scan, the mass of the sample was considered temperature-dependent in Equation (2) (and obtained through a TGA test) while for all the following four scans, the mass of the material was considered constant and with a value equal to the original mass of the sample reduced by a factor given by the value of the TGA of the P10E matrix tested in nitrogen at 550 °C (all the technical details on this analysis are reported in [31]). Figure 5 shows that after a first scan, the original shape of the P10E is completely altered and the irreversible part of the heat capacity is basically no longer present. However, even after four additional scans operated at the maximum temperature allowed by our DSC (550 °C), the capability to convert the powdered P10E into a charred medium was limited. This conclusion can be drawn considering the values of the heat capacity of a fully carbon material [58] reported in the same figure (black dots) which were very different from the 5th scan on the P10E (blue curve). Since the maximum temperature allowed by the DSC used in this analysis was not sufficient to operate a change of the P10E in a fully charred medium, i.e., the organic-to-char conversion rate was below 1 [8], we carried out a TGA dynamic scan in nitrogen up to 1000 °C on the powders obtained at the end of the 5th scan and the resulting powders were studied again into the DSC. As reported in Figure 5, the TGA scan reduced the difference between the heat capacity of a fully carbon-based material (black dots) [58] and the charred powders of the P10E (black line), i.e., after the TGA test, the resulting degree of conversion was closer to 1. Consequently, at the end of all these tests, the original organic matrix is fully charred and the irreversible part of the C p has been completely removed.
Figure 6a reports the C p of the C-P10E as a function of temperature as compared to state-of-the-art CPC formulations. We also exploited the data of the heat capacity of a phenolic matrix studied by NASA in [32] and, through the rule of mixture [31], the C p of a CPC made of 30 wt% of carbon fiber and 70 wt% of matrix (this method was validated in [31,59]), was extrapolated. Figure 6b also reports the ρ of the C-P10E as a function of temperature as compared to state-of-the-art CPC formulations. In this figure, we considered an upper temperature limit of 800 °C, i.e., a value in which all the charring processes can be considered completed.
Moreover, Figure 7 reports the behavior of the C p and of the ρ of the materials considered in Figure 6 in the entire range of temperature typically considered to properly model the ablation response of CPCs, i.e., up to 3000 °C (in this case we were forced to use a log scale in the temperature axis). From these graphs, it is possible to evidence that in terms of heat capacity above about 800 °C, the value of the C p is considered to possess a gradually increasing, monotonic behavior with no presence of irreversible phenomena [8,9]. For the density of CPCs, it is possible to conclude that above about 800 °C, the value of the ρ of the CPC formulations tend to be considered constant as reported in the literature [60,61,62,63]. In general, in a cubic CPC sample exposed to the torch (Figure 7c), the thermal degradation depends on the depth of the sample we consider with the surface exposed to the harshest conditions (see in-depth temperature profiles in Figure 3 or Figure 4). As a result, assuming our OAT testing conditions, the degradation is highest on the surface of the sample while in the cold part of the test coupon, the material is still virgin (Figure 7c,d) and at room temperature (RT), i.e., the thermo-physical properties are not uniform through the sample. More in detail, if we consider the OAT testing condition reported in Section 2.2 for an out-plane CPC sample, it is possible to correlate the maximum temperatures recorded at 0 mm (Ts), 3 mm (T1), 6 mm (T2), and 10 mm (RT), at the beginning (0 s, Figure 7c) and at the end of the test (30 s, Figure 7d) with the values of the C p and of the ρ at these given temperatures (Figure 7a,b) evaluating the values of these properties at the given depths (Figure 7c,d). The green arrows are used to graphically indicate the values of the C p and of the ρ (in Figure 7a,b, respectively) at the beginning of the tests (Figure 7c), while the red arrows are used to graphically indicate the heat capacity and density of the same testing coupon (in Figure 7a,b, respectively) at the end of the OAT test (Figure 7d). At 0 s, independently from the depth, the values of the C p and of the ρ are constant ( C p ( R T ) and ρ R T ,   depth), while at the end of the test (30 s), the values of these thermo-physical properties possess a dependence on the depth ( C p ( T s m a x ) and ρ T s m a x at 0 mm, C p ( T 1 m a x ) and ρ T 1 m a x at 3 mm, C p ( T 2 m a x ) and ρ T 2 m a x at 6 mm, and C p ( R T ) and ρ R T at 10 mm).
Since the values of the C p and of the ρ are not uniform through a CPC sample subjected to an OAT as explained in Section 2.2., if we want to evaluate the mean values of these thermo-physical properties of the entire tested coupon as a function of the exposure time (from 0 s to 30 s), once we define a mesh of the sample, we have to consider the contribution of each slice in terms of value of heat capacity and of density as a function of the temperature (or of the charring degree). As a result, if as an example we consider an out-plane sample of C-P10E exposed to the OAT for 30 s, if we evaluate the evolution of the mean C p and ρ as function of the exposure time, then is possible to obtain the graphs reported in Figure 8.
Figure 8 demonstrates that even considering the dependence of C p and ρ with the temperature, the denominator of Equation (1) can be considered nearly constant, i.e., during the OAT test, the thermal diffusivity depends on the temperature manly through the k .

2.4. Theoretical Background on the Adopted Heat Conduction Models

2.4.1. Ideal Semi-Infinite Solid

The proposed method to study the thermal diffusivity of CPCs is based on the assumption of the semi-infinite solid. Since in principle it extends infinitely in all directions except one, this solid is defined by a singular identifiable surface. If a sudden change in the conditions is imposed at this surface, transient, one-dimensional conduction will occur within the solid. The semi-infinite solid provides a useful idealization for many practical problems. Among the uniqueness of the protocol proposed in this paper it is worth mentioning the relative simplicity to run the experiments while allowing to carry out test that cannot be reproduced by any other common protocol designed to evaluate the thermal diffusivity such as an LFA (especially in terms of heating rates).
Two scenarios will be considered in this approach: constant surface temperature or constant heat flux [64]. The first scenario of constant surface temperature is summarized in Figure 9a and responds to the following equation:
T x , t T S T i T S = e r f x 2 α t
where e r f is the error function, T i is the temperature of the slab at t = 0 , T S is the surface temperature imposed by the torch plume, T x , t is the temperature recorded at x = 3 mm or T1(t) and at x = 6 mm or T2(t) and α is the thermal diffusivity of the tested expressed as in Equation (1). In most of the cases, Equation (3) is used to derive the in-depth temperature profiles T x , t of a given material exposed to a constant surface temperature once the thermal diffusivity of the given material is known:
T x , t = e r f x 2 α t T i T S + T S .
The surface net heat flux q n e t can be obtained by applying Fourier’s law at x = 0 , in which case:
q n e t = k T x x = 0 q n e t = k T S T i α t 1 2 1 π 1 2 = α ρ C p T s T i α t 1 / 2 1 π 1 2
where k is the mean thermal conductivity of the sample. From the knowledge of the T S ( t ) (acquired with a pyrometer as described in Section 2.2), of the T1(t) and of the T2(t) (acquired with thermocouples as described in Section 2.2) and exploiting only the heating phase of the test (Figure 3 and Figure 4), Equation (3) can be solved to retrieve α .
The second scenario of constant heat flux is summarized in Figure 9b and responds to the following equation:
T x , t T i = 2 q 0 α t π k exp x 2 4 α t q 0 x k e r f c x 2 α t
where T i is the temperature of the slab at t = 0 . It is worth remembering that the complementary error function, e r f c w , is defined as e r f c w 1 e r f w , and q 0 is the surface heat flux which can be obtained by applying Fourier’s law at x = 0, in which case:
q n e t = q 0 .
As for the first scenario, from the knowledge of the experimental values of the T S ( t ) (acquired with a pyrometer as described in Section 2.2), of the T1(t) and of the T2(t) (acquired with thermocouples as described in Section 2.2) and exploiting only the heating phase of the test (Figure 1), Equation (6) can be numerically solved to evaluate the α   ( T ) . In Equation (6), the thermal conductivity can be expressed as k = α   ρ   C p and, as a result, from the knowledge of the T x , t , this equation can be solved in α as a function of the temperature directly from the knowledge of ρ and of the C p :
T x , t = q 0 k 2 α t π exp x 2 4 α t x e r f c x 2 α t + T i .
Considering that k = α   ρ   C p , Equation (7) can be directly expressed in terms of α   ρ   C p :
T x , t = q 0 α ρ C p 2 α t π exp x 2 4 α t x e r f c x 2 α t + T i .
Even though Equation (7) has been found considering a constant surface heat flux, we wanted to relate q n e t on the T S t . An explicit expression for the q n e t as a function of the T S t can be obtained from Equation (6) evaluating the T x , t on the surface:
T 0 , t T i = q n e t k 2 α t π
but since T 0 , t is the T S t , Equation (8) can be rewritten as:
k T S t T i = q n e t 4 α t π 1 2 .
Then, the q n e t can be expressed in terms of T S t as follows:
q n e t = k T S t T i α t 1 2 π 4 1 2 .
As a result, if we compare Equations (5) to (9), we can see that in the expression of the   q n e t there is a common term:
k T S t T i α t 1 2
times the term ( 1 / π ) 1 / 2 in the constant surface temperature scenario, or times ( π / 4 ) 1 / 2 in the case of the constant heat flux framework.
Ideally, in a material in which the thermo-physical properties are constant, with a step change in the surface temperature (case 1), the temperatures within the medium monotonically approach T s with increasing t, while the magnitude of the surface temperature gradient, and hence the surface heat flux, decreases as t 1 / 2 . On the other hand, for a constant surface heat flux (case 2), Equation (6) reveals that T 0 , t = T s t increases monotonically as t 1 / 2 .

2.4.2. Heat Transfer Models in Presence of Non-Constant Thermo-Physical Properties

In a polymeric ablative material such a CPC the dependence of α , ρ and of the   C p from the temperature (they are not constant) alter the shape of the temperature profiles. As an example, for a CPC, the extent of the irreversible part of the heat capacity profile of a phenolic matrix which is associated with the decomposition (or charring) of the matrix and which works as a heat sink, influences the in-depth temperature profile of the composite virgin/charred material. Equation (4) is analytically derived from the assumption of constant surface temperature and constant thermo-physical properties. Moreover, the above introduced models also imply uniform α through the thickness of the sample. Nonetheless, it has been showed that this analytical model can be used also in the case of non-constant α , i.e., α ( T ) = k ( T ) / ρ T C p T [65]. We will show that with the proper assumptions, the same conclusion can be drawn for Equation (7).
It should be pointed out that the value of α determined from this type of calculation is not the classical thermal diffusivity that is a conduction property of a material as summarized in Equation (1). Rather, it is a pseudothermal diffusivity that includes the effects of conduction, chemical reactions within the char, and, most important, radiation within the char. The value of α determined from Equations (4) or (7) starting from the in-depth temperature profile is the average value between the exposed surface ( T S ) and the internal isotherm ( T x , t ) for which the thermal diffusivity is being calculated [65]. As a result, from a mathematical point of view, keeping in mind the considerations introduced in Section 2.2, the values of the thermal diffusivity derived from Equations (4) or (7) can be considered as mean values ( α ) of the region comprised between the surface and the given ( x ). However, in terms of heat transfer analysis the mean values of the thermal diffusivity determined from Equations (4) or (7) keep into account the effect of the entire sample.

3. Results and Discussion

3.1. Post-OAT Appearance of the Tested CPCs

Even though this paper is focused on the evaluation of the thermal diffusivity of CPCs, some comments on the post-OAT appearance of the tested formulations will be made. In fact, in a previous paper [2], we assessed the dependence of the Inter Laminar Shear Strength (ILSS) of the C-P and C-P10E formulations as a function of the temperature. In [2], it was reported the positive effect of the epoxy phase in the phenolic matrix in terms of capability to boost the ILSS of the resulting composites. However, in that study, we did not focus the attention on the capability of the C-P and C-P10E formulations tested in-plane with the OAT. According to Figure 2, it is possible to evidence that in the in-plane OAT tests the neat C-P tended to behave better than the C-P10E since the sample of the hybrid epoxy–phenolic-based CPC tended to be more delaminated in correspondence of the exposed surface. This behavior which can be evidently attributed to the polymeric phase, i.e., once fully charred, the hybrid epoxy–phenolic matrix P10E tended to degrade more than the neat phenolic matrix and, in other terms, even the addition of 10 wt% of epoxy on the phenolic matrix tended to reduce the char yield of the resulting P10E. This effect could not be directly evidenced in the out-plane with hybrid virgin/charred coupons studied in the out-plane configuration in [2,3]. To better understand these results, it is worth noting that during the firing of the motor [66] the inner layers of the nozzle wall start to char (Figure 10a, in red) dramatically reducing the local ILSS (over 90%). As a result, the structural integrity of the nozzle wall is enabled by the virgin portion (Figure 10a, in green) of the overlapping fiber layers. The charred part does not contribute to the structural integrity of the nozzle wall which is subjected to ply-lift. Depending on the combustion time, every design and layout of the fibers have to keep into account the fact that a certain portion of the CPC wall has to be left virgin to avoid the failure of the nozzle. That is the reason why it is imperative to have CPCs with the highest ILSS and the lower thermal conductivity. The part of throat evidenced in the circular red area of Figure 10b is completely charred as compared to other portions of the same part which still possess a certain level of virgin material—as resumed in Figure 10a. In Figure 10b, the triangular part evidenced in red is the portion that does not possess any virgin material and thus is completely removed.
Consequently, although in the in-plane OAT tests, the fully charred formulation C-P10E appears worse than the completely carbonized C-P neat material (in terms of delamination), it is worth remarking that the erosion resistance of the given part of the SRM nozzle assembly has to be evaluated in the hybrid virgin/charred configuration (as in Figure 10a). In other words, the overall delamination (especially in the in-plane configuration as reported in Figure 10c) and erosion resistance of the given item depend on the ratio between the ILSS of the charred CPC as compared to the ILSS of virgin backup material (Figure 10b). With this respect, since this paper is not focused on the study of the relative merit of the C-P as compared to the C-P10E in terms of ILSS, a discussion on the loss of the mechanical properties of fully charred samples of CPCs (under isothermal condition, i.e., similarly to the in-plane coupons studied in this research) as compared to hybrid virgin/charred coupons (out-plane) has been addressed in [2].

3.2. Evaluation of the Thermal Diffusivity: Constant Surface Temperature Scenario

The most representative temperature profiles of each formulation C-P and C-P10E, in the in-plane and out-plane configurations (Figure 3 and Figure 4) were used to evaluate the in-plane and out-plane thermal diffusivity through Equation (3). The numerical solution of this equation with MATLAB (R2023a) did not imply any particular difficulty. Figure 11 reports the out-plane thermal diffusivity α   ( T ) obtained with these OAT test exploiting the T1(t) and T2(t) temperature profiles, as compared to the corresponding data referred to the MX-4926 and FM-5014 [60,62,63,67,68]. The out-plane thermal diffusivity of the C-P formulation evaluated with an LFA test is also reported (unfortunately, the LFA results in-plane are not available). In the range of temperatures considered in our LFA tests, the C-P and the C-P10E have an out-plane thermal diffusivity qualitatively in line with the state-of-the-art CPC formulations. However, from a quantitative point of view, when adopted to retrieve the thermal diffusivity, Equation (3) produced numerical results lower than the values of state-of-the-art CPCs (such as the FM-5014) manufactured with low-fired stretch-broken PAN-CFs and even of those made with Rayon-CFs (such as the MX-4926). Moreover, the obtained curves of α   ( T ) tended to exhibit a dependence from the temperature less marked than the in the case of the MX-4926 and FM-5014, also considering that, in principle, the CPCs tested in this research should exhibit a thermal diffusivity higher than the values referred to the FM-5014. These results encouraged further efforts to increase the effectiveness of the protocol aimed at retrieving the thermal diffusivity using the in-depth temperature profiles and using the modified constant heat flux scenario.

3.3. Evaluation of the Thermal Diffusivity: Modified Constant Heat Flux Scenario

As for the constant surface temperature scenario, the most representative temperature profiles of each formulation C-P and C-P10E, in plane and out-plane (Figure 3 and Figure 4) were used again to evaluate the out-plane and in-plane thermal diffusivity through Equation (7). The numerical solution of this equation with MATLAB resulted in being more difficult than the previous case. In fact, analytically speaking, while for the determination of α , Equation (3) basically depends on the experimental values of x, T S ( t ) , T 1 ( t ) and T 2 t , on the other hand, Equation (7) explicitly also keeps into account the density ρ ( T ) , the heat capacity C p ( T ) , as well as the values of the net surface heat flux q n e t on the surface of the sample, increasing the difficulties to determine the values of α . However, as reported in Figure 8, since the term ρ ( T ) C p ( T ) tends to be nearly constant during the entire exposure time of the OAT test, most of the criticalities on the determination of α though Equation (7) is related to the values of the net surface heat flux q n e t . Due to the difficulties in solving Equation (7) using the correct values of q n e t and the reading of the heat flux from the slug calorimeter, we addressed the possibility of using Equation (9) for the determination of the net surface heat flux. In other words, since we did not explicitly know the exact values of the q n e t , we decided to correlate the values of the net surface heat flux to the surface temperature T S ( t ) through Equation (9) also in the constant heat flux scenario. As a result, on the surface of the exposed sample, if we introduce Equation (9) into Equation (7), we analytically obtain Equation (10):
T x , t = T S t T i exp x 2 4 α t π 4 α t 1 2 x e r f c x 2 α t + T i
i.e., a new equation which correlates the q n e t at the T S ( t ) as in the first scenario. Solving this equation, we obtain the results reported in Figure 12.
Figure 12 shows that the values of thermal diffusivity obtained in the second scenario are higher than the corresponding obtained in the constant surface temperature framework, i.e., they align more closely with the values derived from the LFA, highlighting the relatively superior predictive abilities of the modified constant heat flux model. However, especially for the out-plane samples, in both scenarios it is possible to see that above about 600 °C the results from our protocol tend to be lower than the values obtained from LFA or from the literature. This difference can be explained as follows.

4. Concluding Considerations

Starting from the equations for the in-depth temperatures obtained in the framework of semi-infinite solid and constant surface temperature and constant heat flux, two protocols aimed at retrieving the thermal diffusivity as a function of temperature of carbon/phenolic composites tested with an oxy-acetylene torch were validated. Both models showed to be feasible and exploitable to determine the thermal diffusivity of CPCs starting from the sole experimental in-depth temperature profiles. Up to a temperature of about 600 °C, the values of thermal diffusivity of the studied CPCs were found to be consistent with state-of-the-art formulations (MX-4926 and FM-5014) both in the in-plane as well as in the out-plane configurations.
The first reason why in general the experimental values of the thermal diffusivity tend to show a limited dependency from the temperature above about 600 °C is due to the fact that when Equations (4) and (10) are used to numerically retrieve α , these relationships are not able to directly find a function that strongly depends on the temperature. In other words, once Equations (4) or (10) is solved, i.e., reversed to determine α starting from the given temperature profiles, both models (constant surface temperature or constant heat flux) tend to produce (at the given time t ) a thermal diffusivity value that fits for the whole sample independently from the degradation degree. At the give time, and depth (temperature), these equations work determining a single value of thermal diffusivity that, on the other hand, is the same for the entire sample. To make it very short, Equations (4) and (10) derive from the assumption of constant thermal diffusivity. This point will be further addressed in another (second forthcoming) paper.
The second main limitation related to our experiments is that considering samples like the ones reported in Figure 2—which are significantly different from being regarded as a semi-infinite solid—the experimental in-depth temperature profiles tend to be lower than in the ideal case. In fact, Figure 13a reports the boundary conditions in terms of input and output heat fluxes (or better power) for the exposed sample during the OAT test, i.e., when the torch is turned on [69]:
q n e t A i = q s A i q r a d A o q c o n v A o
where q n e t is the net heat flux that enters in the CPC sample at the surface x = 0 , q s is the heat flux of the torch, and q r a d and q c o n v are the heat fluxes losses by radiation and convection, A i is the surface in which we consider entering the heat flux ( A i = L 2 , i.e., the surface of a face of the cubic sample) and A o is the surface that we consider responsible for the output of the heat flux since we are not in an ideal condition of semi-infinite solid.
Equation (11) defines the net heat flux q n e t entering into the CPC sample tested with an OAT having an exposed surface at the temperature T 0 , t = T s ( t ) . q n e t is then defined as a difference between q s , i.e., the heat flux entering in the surface from the energy source (OAT) and the heat loss due to radiation ( q r a d ) and convection ( q c o n v ).
As an example, if we set the entire output surface A o as equal to A i + 6 × 4 × 0.1 × A i (see Figure 13b) and each slice of the cubic sample lateral surface is put at a temperature as reported in the same figure—trying to mimic the real in-depth (and lateral, supposing to have perfectly plane iso-thermal surfaces for the tested CPC sample (Figure 13c)—the contribution of the different heat fluxes in Equation (11) could be defined as follows:
q n e t t A i = q s t A i + ε σ T s 4 t T s u r r 4 A i + ε σ T s t 4 T s u r r 4 A L + ε σ T s t + T 1 t 2 4 T s u r r 4 A L + ε σ T 1 t 4 T s u r r 4 A L + ε σ T 1 t + T 2 t 2 4 T s u r r 4 A L + ε σ T 2 t 4 T s u r r 4 A L + ε σ T 2 t + R T 2 4 T s u r r 4 A L + h T s t T s u r r A i + h T s t T s u r r A L + h T s t + T 1 t 2 T s u r r A L + h T 1 t T s u r r A L + h T 1 t + T 2 t 2 T s u r r A L + h T 2 t T s u r r A L + h T 2 t + R T 2 T s u r r A L
where A i is the surface in which we consider entering the heat flux, while each isothermal lateral slice defined as A L = 4 0.1 L L = 0.4 A i —in which we consider exiting the heat flux—is set at a specific temperature as reported in Equation (12), ε is the emissivity of the tested material, σ is the Stefan–Boltzmann constant, h is the average surface convective heat transfer (which was determined experimentally from the cooling part of the temperature profiles), T s u r r , the surrounding temperature of the sample, corresponds to the room temperature ( T s u r r = R T ). In other words, in our example of possible loss of energy of the real sample (as compared to a semi-infinite solid), to properly assess the heat loss due to the lateral faces of the tested coupon, the contribution of Equation (12) between square brackets is the entire loss by radiation while the contribution between brace brackets is the entire loss by convection.
Starting from Equations (11) and (12), we preliminary determined that the sum of the heat loss (power loss) term between square brackets ( q r a d   A L ) and of the contribution between brace brackets ( q c o n v   A L ) are extremely relevant and can also be substantially higher that the net entering heat term ( q n e t A i ). This evidence would support the considerations on the lower temperature profiles related to the real samples as compared to the ideal semi-infinite slab.
In fact, starting from the assumptions at the beginning of this section on the heat losses summarized in Figure 13, if we define the term ratio as
R ( t ) = q n e t ( t ) A i + ( q r a d ( t ) + q c o n v ( t ) ) A L q n e t ( t ) A i
and we plot it as a function of the time for the C-P10E out-plane and in-plane using the constant heat flux scenario data, we obtained the graph presented in Figure 14. For the C-P10E out-plane, the loss of heat on the lateral surfaces of the sample is higher than in the case of the in-plane test. This result can be justified considering that the in-plane thermal diffusivity of the C-P10E is higher than the out-plane counterpart. In other words, as compared to the C-P10E in-plane samples, for the out-plane coupons a higher amount of heat is spread towards the lateral faces of the samples and then dissipated by radiation and convection on the four lateral faces of the coupon. This behavior would also imply the general better capability of the protocol to determine the thermal diffusivity in the case of the in-plane samples as compared to the out-plane coupons.
As a conclusion, it is worth highlighting that in the scenario of an ideal semi-infinite slab, all the temperatures introduced in this work ( T s t , T 1 t , T 2 t ) would be rescaled and increased thus modifying the values of thermal diffusivity. We would expect that a higher value of T 1 t / T s t would increase the value of α . In a forthcoming additional paper on the same topic, we will try to evaluate the effect of the loss of heat related to the lateral faces of the cubic sample on the values of the thermal diffusivity.
In further research, we will also try to consider the effect of the pyrolysis gases on the in-depth temperature profiles. In fact, phenolics possess a substantial amount of carbon-hydrogen bonds, which supply the evolution gases (such as H2 and CH4 and others) that contribute to the reduction in the local temperature and also provide the transpiration cooling [70,71,72,73,74,75,76,77]. At the same time, once the gas reached the surface, this gaseous film also works as a reaction barrier; however, when this gaseous layer is removed, due to the interaction of the hyper-thermal environment gases, chemical erosion tends to increase.

5. Conclusions

A protocol aimed at determining the thermal diffusivity of Carbon/Phenolic Composites (CPCs) through the Oxy-Acetylene Torch (OAT) tests was validated and presented in this paper. The proposed approach was demonstrated to bypass some limitations of common or commercial protocols, especially in terms of capability in determining the thermal diffusivity of CPCs at high heating rates and through experiments which reconstruct the real hyperthermal conditions of the TPS material. The proposed method to retrieve the thermal diffusivity of CPCs is based on the assumption of the semi-infinite solid. Two scenarios were considered in this approach: constant surface temperature or constant heat flux. The values of in-plane and out-plane thermal diffusivity measured for some experimental CPC laminates based on Polyacrylonitrile (PAN) based on Carbon Fibers (PAN-CFs) fabric resulted to be in line with the values of a state-of-the-art CPC such as the FM-5014 (manufactured with low-fired stretch-broken PAN-CFs). Finally, a thorough comparison of the comparison of the obtained data with state-of-the-art CPCs such as MX-4926 and FM-5014 has also been carried out.

Author Contributions

Conceptualization, M.N. and L.T.; Methodology, M.N.; Software, M.R.; Validation, M.N. and M.R.; Formal analysis, M.N.; Investigation, M.N. and M.R.; Data curation, M.N. and M.R.; Writing—original draft, M.N.; Writing—review & editing, M.N., L.T. and M.R.; Supervision, M.N. and L.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Correction Statement

This article has been republished with a minor correction to the image quality of all figures and equation (12). This change does not affect the scientific content of the article.

References

  1. Natali, M.; Kenny, J.M.; Torre, L. Science and Technology of Polymeric Ablative Materials for Thermal Protection Systems and Propulsion Devices: A Review. Prog. Mater. Sci. 2016, 84, 192–275. [Google Scholar] [CrossRef]
  2. Natali, M.; Torre, L.; Rallini, M. Interlaminar Properties, Ablation Resistance and Thermal Diffusivity of Carbon/Phenolic Composites (CPCs) Derived from Optimized High Char Yield Matrices: A Comparison with MX-4926 and FM-5014. Compos. Part. A Appl. Sci. Manuf. 2023, 175, 107801. [Google Scholar] [CrossRef]
  3. Natali, M.; Rallini, M.; Torre, L.; Puglia, D. High Temperature Composites from Renewable Resources: A Perspective on Current Technological Challenges for the Manufacturing of Non-Oil Based High Char Yield Matrices and Carbon Fibers. Front. Mater. 2022, 9, 805131. [Google Scholar] [CrossRef]
  4. Mouritz, A.; Gibson, A. Structural Properties of Composites in Fire. In Fire Properties of Polymer Composite Materials; Springer: Dordrecht, The Netherlands, 2006; pp. 163–213. [Google Scholar]
  5. Clayton, W.A.; Kennedy, P.B.; Evans, R.J.; Cotton, J.E.; Francisco, A.C.; Fabish, T.J.; Eldridge, E.A.; Lagedrost, J.F. Thermal properties of ablative chars. Report by the aerospace group of the boeing company under USAF, Contract Number AF 33(615)-3804. 1967. [Google Scholar]
  6. Mouritz, A.; Gibson, A. Modelling the Thermal Response of Composites in Fire. In Fire Properties of Polymer Composite Materials; Springer: Dordrecht, The Netherlands, 2006; pp. 133–161. [Google Scholar]
  7. Mouritz, A.P.; Feih, S.; Kandare, E.; Mathys, Z.; Gibson, A.G.; Des Jardin, P.E.; Case, S.W.; Lattimer, B.Y. Review of Fire Structural Modelling of Polymer Composites. Compos. Part A Appl. Sci. Manuf. 2009, 40, 1800–1814. [Google Scholar] [CrossRef]
  8. Yu, B.; Till, V.; Thomas, K. Modeling of Thermo-Physical Properties for FRP Composites under Elevated and High Temperature. Compos. Sci. Technol. 2007, 67, 3098–3109. [Google Scholar] [CrossRef]
  9. Lattimer, B.Y.; Ouellette, J.; Trelles, J. Thermal Response of Composite Materials to Elevated Temperatures. Fire Technol. 2011, 47, 823–850. [Google Scholar] [CrossRef]
  10. Henderson, J.B.; Wiebelt, J.A.; Tant, M.R. A Model for the Thermal Response of Polymer Composite Materials with Experimental Verification. J. Compos. Mater. 1985, 19, 579–595. [Google Scholar] [CrossRef]
  11. Henderson, J.B.; Wiecek, T.E. A Mathematical Model to Predict the Thermal Response of Decomposing, Expanding Polymer Composites. J. Compos. Mater. 1987, 21, 373–393. [Google Scholar] [CrossRef]
  12. Florio, J., Jr.; Henderson, J.B.; Test, F.L. Measurement of the Thermochemical Expansion of Porous Composite Materials. High Temp. High Press. 1989, 21, 157–165. [Google Scholar]
  13. Tant, M.R.; Henderson, J.B.; Boyer, C.T. Measurement and Modelling of the Thermochemical Expansion of Polymer Composites. Composites 1985, 16, 121–126. [Google Scholar] [CrossRef]
  14. Henderson, J.B.; Dohert, M.P. Measurement of Selected Properties of a Glass-Filled Polymer Composite. High. Temp. High. Press. 1987, 19, 95–102. [Google Scholar]
  15. Florio, J., Jr.; Henderson, J.B.; Test, F.L.; Hariharan, R. A Study of the Effects of the Assumption of Local-Thermal Equilibrium on the Overall Thermally-Induced Response of a Decomposing, Glass-Filled Polymer Composite. Int. J. Heat. Mass. Transf. 1991, 34, 135–147. [Google Scholar] [CrossRef]
  16. Sullivan, R.M. A Finite Element Model for Thermochemically Decomposing Polymers; The Pennsylvania State University: State College, PA, USA, 1990. [Google Scholar]
  17. Sullivan, R.M.; Salamon, N.J. A Finite Element Method for the Thermochemical Decomposition of Polymeric Materials—I. Theory. Int. J. Eng. Sci. 1992, 30, 431–441. [Google Scholar] [CrossRef]
  18. Sullivan, R.M.; Salamon, N.J. A Finite Element Method for the Thermochemical Decomposition of Polymeric Materials—II. Carbon Phenolic Composites. Int. J. Eng. Sci. 1992, 30, 939–951. [Google Scholar] [CrossRef]
  19. Sullivan, R.M. A Coupled Solution Method for Predicting the Thermostructural Response of Decomposing, Expanding Polymeric Composites. J. Compos. Mater. 1993, 27, 408–434. [Google Scholar] [CrossRef]
  20. Mcmanus, H.L.N.; Springer, G.S. High Temperature Thermomechanical Behavior of Carbon-Phenolic and Carbon-Carbon Composites, I. Analysis. J. Compos. Mater. 1992, 26, 206–229. [Google Scholar] [CrossRef]
  21. Mcmanus, H.L.N.; Springer, G.S. High Temperature Thermomechanical Behavior of Carbon-Phenolic and Carbon-Carbon Composites, II. Results. J. Compos. Mater. 1992, 26, 230–255. [Google Scholar] [CrossRef]
  22. Dimitrienko, Y.I. Thermal Stresses and Heat-Mass Transfer in Ablating Composite Materials. Int. J. Heat. Mass. Transf. 1995, 38, 139–146. [Google Scholar] [CrossRef]
  23. Dimitrienko, Y.I. Thermomechanical Behaviour of Composite Materials and Structures under High Temperatures: 1. Materials. Compos. Part A Appl. Sci. Manuf. 1997, 28, 453–461. [Google Scholar] [CrossRef]
  24. Dimitrienko, Y.I. Thermomechanical Behaviour of Composite Materials and Structures under High Temperatures: 2. Structures. Compos. Part A Appl. Sci. Manuf. 1997, 28, 463–471. [Google Scholar] [CrossRef]
  25. Gibson, A.G.; Wu, Y.-S.; Chandler, H.W.; Wilcox, J.A.D.; Bettess, P. Model for the Thermal Performance of Thick Composite Laminates in Hydrocarbon Fires. Rev. Institute Fr. Pet. 1995, 50, 69–74. [Google Scholar] [CrossRef]
  26. Dodds, N.; Gibson, A.G.; Dewhurst, D.; Davies, J.M. Fire Behaviour of Composite Laminates. Compos. Part A Appl. Sci. Manuf. 2000, 31, 689–702. [Google Scholar] [CrossRef]
  27. Zheng, F.; Ren, Z.; Xu, B.; Wan, K.; Cai, J.; Yang, J.; Zhang, T.; Wang, P.; Niu, B.; Zhang, Y.; et al. Elucidating Multiple-Scale Reaction Behaviors of Phenolic Resin Pyrolysis via TG-FTIR and ReaxFF Molecular Dynamics Simulations. J. Anal. Appl. Pyrolysis 2021, 157, 105222. [Google Scholar] [CrossRef]
  28. Jiang, H.; Wang, J.; Wu, S.; Wang, B.; Wang, Z. Pyrolysis Kinetics of Phenol–Formaldehyde Resin by Non-Isothermal Thermogravimetry. Carbon. N.Y. 2010, 48, 352–358. [Google Scholar] [CrossRef]
  29. Torre, L.; Kenny, J.M.; Maffezzoli, A.M. Degradation Behaviour of a Composite Material for Thermal Protection Systems Part I—Experimental Characterization. J. Mater. Sci. 1998, 33, 3137–3143. [Google Scholar] [CrossRef]
  30. Torre, L.; Kenny, J.M.; Maffezzoli, A.M. Degradation Behaviour of a Composite Material for Thermal Protection Systems. Part II—Process Simulation. J. Mater. Sci. 1998, 33, 3145–3149. [Google Scholar] [CrossRef]
  31. Natali, M.; Torre, L.; Puri, I.; Rallini, M. Thermal Degradation of Phenolics and Their Carbon Fiber Derived Composites: A Feasible Protocol to Assess the Heat Capacity as a Function of Temperature through the Use of Common DSC and TGA Analysis. Polym. Degrad. Stab. 2022, 195, 109793. [Google Scholar] [CrossRef]
  32. Sykes, G.F. Decomposition Characteristics of a Char-Forming Phenolic Polymer Used for Ablative Composites; National Aeronautics and Space Administration: Washington, DC, USA, 1967.
  33. Shilav, R.; Leitner, A.; Gany, A. Development of Thermal Conductivity Apparatus for Ablative Composite Materials. In Proceedings of the 57th Israel Annual Conference on Aerospace Sciences, Tel-Aviv & Haifa, Israel, 15–16 March 2017. [Google Scholar]
  34. ASTM C1113/C1113M; Standard Test Method for Thermal Conductivity of Refractories by Hot Wire (Platinum Resistance Thermometer Technique). iTeh Standards: Toronto, Italy, 2013.
  35. Roder, H.M.; Perkins, R.A.; Laesecke, A.; Nieto de Castro, C.A. Absolute Steady-State Thermal Conductivity Measurements by Use of a Transient Hot Wire System. J. Res. Natl. Inst. Stand. Technol. 2000, 105, 221–253. [Google Scholar] [CrossRef] [PubMed]
  36. Shilav, R.; Leitner, A.; Gany, A. New Measurement Technique for Thermal Conductivity of Ablative Materials. Int. J. Energetic Mater. Chem. Propuls. 2018, 17, 171–182. [Google Scholar] [CrossRef]
  37. Dimitrienko, Y.I. Thermomechanics of Composites under High Temperatures; Springer: Dordrecht, The Netherlands, 1999; Volume 65, ISBN 978-90-481-5122-6. [Google Scholar]
  38. Bahramian, A.R.; Kokabi, M. Ablation Mechanism of Polymer Layered Silicate Nanocomposite Heat Shield. J. Hazard. Mater. 2009, 166, 445–454. [Google Scholar] [CrossRef] [PubMed]
  39. Paul, A.; Binner, J.G.P.; Vaidhyanathan, B.; Heaton, A.C.J.; Brown, P.M. Heat Flux Mapping of Oxyacetylene Flames and Their Use to Characterise Cf-HfB2 Composites. Adv. Appl. Ceram. 2016, 115, 158–165. [Google Scholar] [CrossRef]
  40. Zibitsker, A.; Berreby, M.; Michaels, D.; Shilav, R.; Frisman, I. Ultrasonic Temperature Compensating Method for Tracking Decomposition Front in Silica-Phenolic Thermal Protection Material. J. Thermophys. Heat. Trans. 2021, 35, 770–787. [Google Scholar] [CrossRef]
  41. Minges, M.L. Thermophysical Characteristics of High-Performance Ablative Composites. In Ablative Plastics; D’Aelio, G.F., Parker, J.A., Eds.; Marcel Dekker Inc.: New York, NY, USA, 1971; pp. 287–312. [Google Scholar]
  42. Myrick, J.A.; Keyhani, M.; Frankel, J.I. Determination of Surface Heat Flux and Temperature Using in-Depth Temperature Data—Experimental Verification. Int. J. Heat. Mass. Transf. 2017, 111, 982–998. [Google Scholar] [CrossRef]
  43. Gazarik, M.J.; Wright, M.J.; Little, A.; Cheatwood, F.M.; Herath, J.A.; Munk, M.M.; Novak, F.J.; Martinez, E.R. Overview of the MEDLI Project. In Proceedings of the 2008 IEEE Aerospace Conference, Big Sky, MT, USA, 1–8 March 2008; IEEE: Piscataway, NJ, USA, 2008; pp. 1–12. [Google Scholar]
  44. Bose, D.; Santos, J.A.; Rodriguez, E.; White, T.R.; Mahzari, M. Mars Science Laboratory Heat Shield Instrumentation and Arc Jet Characterization. In Proceedings of the 44th AIAA Thermophysics Conference, San Diego, CA, USA, 24–27 June 2013; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 2013. [Google Scholar]
  45. Morzhukhina, A.V.; Alifanov, O.M.; Budnik, S.A.; Nenarokomov, A.V.; Titov, D.M.; Delfini, A.; Pastore, R.; Santoni, F.; Albano, M.; Marchetti, M. Estimation of the Thermophysical Properties of C/C Plates with Ceramic Nanocoating at Different SiO2 Filling for Aerospace Applications. Aerotec. Missili Spaz. 2023, 103, 3–15. [Google Scholar] [CrossRef]
  46. Alifanov, O.M.; Budnik, S.A.; Nenarokomov, A.V.; Netelev, A.V.; Okhapkin, A.S. Thermophysical Characteristics of Fibrous Thermoprotective Materials at High Temperatures. J. Eng. Phys. Thermophys. 2021, 94, 1052–1062. [Google Scholar] [CrossRef]
  47. Alifanov, O.M.; Budnik, S.A.; Mikhaylov, V.V.; Nenarokomov, A.V.; Titov, D.M.; Yudin, V.M. An Experimental–Computational System for Materials Thermal Properties Determination and Its Application for Spacecraft Structures Testing. Acta Astronaut. 2007, 61, 341–351. [Google Scholar] [CrossRef]
  48. Zhou, J.; Zhang, Y.; Chen, J.K.; Feng, Z.C. Inverse Heat Conduction in a Composite Slab with Pyrolysis Effect and Temperature-Dependent Thermophysical Properties. J. Heat Transf. 2010, 13, 034502. [Google Scholar] [CrossRef]
  49. The European Space Agency. N° 46–2023: Vega-C Zefiro40 Test: Independent Enquiry Commission Announces Conclusions. 2023. Available online: https://www.esa.int/Newsroom/Press_Releases/Vega-C_Zefiro40_Test_Independent_Enquiry_Commission_announces_conclusions (accessed on 19 January 2024).
  50. NASA. NASA Shares Progress Toward Early Artemis Moon Missions with Crew. 2024. Available online: https://www.nasa.gov/news-release/nasa-shares-progress-toward-early-artemis-moon-missions-with-crew (accessed on 19 January 2024).
  51. Natali, M.; Monti, M.; Kenny, J.; Torre, L. Synthesis and Thermal Characterization of Phenolic Resin/Silica Nanocomposites Prepared with High Shear Rate-Mixing Technique. J. Appl. Polym. Sci. 2011, 120, 2632–2640. [Google Scholar] [CrossRef]
  52. Natali, M.; Monti, M.; Kenny, J.M.; Torre, L. A Nanostructured Ablative Bulk Molding Compound: Development and Characterization. Compos. Part. A Appl. Sci. Manuf. 2011, 42, 1197–1204. [Google Scholar] [CrossRef]
  53. Natali, M.; Koo, J.H.; Allcorn, E.; Ezekoye, O.A. An In-Situ Ablation Recession Sensor for Carbon/Carbon Ablatives Based on Commercial Ultra-Miniature Thermocouples. Sens. Actuators B Chem. 2014, 196, 46–56. [Google Scholar] [CrossRef]
  54. Vojvodich, N.S. Hypervelocity Heat Protection—A Review of Laboratory Experiments. In Ablative Plastics; Marcel Dekker: New York, NY, USA, 1971; pp. 41–68. ISBN 13: 9780824711207. [Google Scholar]
  55. Koubek, F.J. A Review of Ablative Studies of Interest to Naval Applications. In Ablative Plastics; Marcel Dekker: New York, NY, USA, 1971; pp. 69–83. [Google Scholar]
  56. Natali, M.; Puri, I.; Kenny, J.M.; Torre, L.; Rallini, M. Microstructure and Ablation Behavior of an Affordable and Reliable Nanostructured Phenolic Impregnated Carbon Ablator (PICA). Polym. Degrad. Stab. 2017, 141, 84–96. [Google Scholar] [CrossRef]
  57. Natali, M.; Rallini, M.; Puglia, D.; Kenny, J.; Torre, L. EPDM Based Heat Shielding Materials for Solid Rocket Motors: A Comparative Study of Different Fibrous Reinforcements. Polym. Degrad. Stab. 2013, 98, 2131–2139. [Google Scholar] [CrossRef]
  58. Toulouklan, Y.S.; Saxens, S.C.; Klemena, P.G. Thermophysical Properties of Matter—The TPRC Data 1976, 6s.
  59. Natali, M.; Puri, I.; Rallini, M.; Kenny, J.; Torre, L. Ablation Modeling of State of the Art EPDM Based Elastomeric Heat Shielding Materials for Solid Rocket Motors. Comput. Mater. Sci. 2016, 111, 460–480. [Google Scholar] [CrossRef]
  60. Shieds, S. Scout Nozzle Data Book (No. NASA-CR-145136), 1976; NASA: Washington, DC, USA, 1976.
  61. Schaefer, J.W.; Dahm, T.J. Studies of Nozzle Ablative Material Performance for Large Solid Boosters. NASA CR-72080 Aerotherm Report No. 66-2; NASA: Washington, DC, USA, 1966.
  62. Clayton, W.A.; Kennedy, P.B.; Evans, R.J.; Cotton, J.E.; Francisco, A.C.; Fabish, T.J.; Eldridge, E.A.; Lagedrost, J.F. AFML-TR-67-413; Battelle Memorial Institute: Columbus, OH, USA, 1968.
  63. Pears, C.D.; Engelke, W.T.; Thornburgh, J.D. The Thermal and Mechanical Properties of Five Ablative Reinforced Plastics from Room Temperature to 750°F. U.S. Air Force Report; AFML-TR-65-133; Southern Research Institute: Birmingham, AL, USA, 1965. [Google Scholar]
  64. Incropera, F.P.; Dewitt, D.P.; Bergman, T.L.; Lavine, A.S. Transient Conduction. In Fundamentals of Heat and Mass Transfer; John Wiley & Sons: Hoboken, NJ, USA, 2007; pp. 255–319. [Google Scholar]
  65. Liston, E.M. Arc-Image Testing of Ablation Materials. In Ablative Plastics; D’Aelio, G.F., Parker, J.A., Eds.; Marcel Dekker Inc.: New York, NY, USA, 1971; pp. 379–407. [Google Scholar]
  66. Mathis, J.R.; Laramee, R.C. Development of Low Cost Ablative Nozzles for Solid Propellant Rocket Motors; Thiokol Chemical Corporation: Brigham City, UT, USA, 1970; Volume II. [Google Scholar]
  67. Wang, W.-X.; Dgany, O.; Wolf, S.G.; Levy, I.; Algom, R.; Pouny, Y.; Wolf, A.; Marton, I.; Altman, A.; Shoseyov, O. Aspen SP1, an Exceptional Thermal, Protease and Detergent-Resistant Self-Assembled Nano-Particle. Biotechnol. Bioeng. 2006, 95, 161–168. [Google Scholar] [CrossRef]
  68. Lachaud, J.; Cozmuta, I.; Mansour, N.N. Multiscale Approach to Ablation Modeling of Phenolic Impregnated Carbon Ablators. J. Spacecr. Rockets 2010, 47, 910–921. [Google Scholar] [CrossRef]
  69. Lyon, R.E. Heat Release Kinetics. Fire Mater. 2000, 24, 179–186. [Google Scholar] [CrossRef]
  70. D’alelio, G.F. Structural Design and Thermal Properties of Polymers. J. Macromol. Sci. Part A Chem. 1969, 3, 411–452. [Google Scholar] [CrossRef]
  71. Fyfe, C.A.; McKinnon, M.S.; Rudin, A.; Tchir, W.J. Investigation of the Mechanism of the Thermal Decomposition of Cured Phenolic Resins by High-Resolution 13C CP/MAS Solid-State NMR Spectroscopy. Macromolecules 1983, 16, 1216–1219. [Google Scholar] [CrossRef]
  72. Rastogi, R.P.; Deepak, D. Char Formation in Ablatives. AIAA J. 1974, 12, 114–116. [Google Scholar] [CrossRef]
  73. Morterra, C.; Low, M.J.D.I.R. Studies of Carbons—VII. The Pyrolysis of a Phenol-Formaldehyde Resin. Carbon. N. Y. 1985, 23, 525–530. [Google Scholar] [CrossRef]
  74. Lum, R.; Wilkins, C.W.; Robbins, M.; Lyons, A.M.; Jones, R.P. Thermal Analysis of Graphite and Carbon-Phenolic Composites by Pyrolysis-Mass Spectrometry. Carbon. N. Y. 1983, 21, 111–116. [Google Scholar] [CrossRef]
  75. Blazso, M.; Tóth, T. Sequence of Phenolic Units in Phenol—Formaldehyde Polycondensates Studied by Pyrolysis—Gas Chromatography/Mass Spectrometry. J. Anal. Appl. Pyrolysis 1991, 19, 251–263. [Google Scholar] [CrossRef]
  76. Schmidt, D. Hypersonic Atmospheric Flight. In Environmental Effects on Polymeric Materials, Environments; Rosato, D.V., Schwartz, R.T., Eds.; Interscience Publishers: Geneva, Swizterland, 1968; pp. 488–578. [Google Scholar]
  77. Schmidt, D. Ablative Polymers in Aerospace Technology. In Ablative Plastics; D’Aelio, G.F., Parker, J.A., Eds.; Marcel Dekker Inc.: New York, NY, USA, 1971; pp. 1–35. [Google Scholar]
Figure 1. CPC sample and torch tip (side-view) (a). The inner core of the torch has an estimated diameter of 4 mm. A typical heat impulse produced by the torch (b). Laser pointer of the pyrometer indicating the exact point in which the surface temperature was read (c).
Figure 1. CPC sample and torch tip (side-view) (a). The inner core of the torch has an estimated diameter of 4 mm. A typical heat impulse produced by the torch (b). Laser pointer of the pyrometer indicating the exact point in which the surface temperature was read (c).
Polymers 16 00577 g001
Figure 2. Virgin in-plane and out-plane OAT samples of the C-P (a) and C-P10E (b). Post-OAT testing appearance of the in-plane and out-plane samples of the C-P (c) and C-P10E (d).
Figure 2. Virgin in-plane and out-plane OAT samples of the C-P (a) and C-P10E (b). Post-OAT testing appearance of the in-plane and out-plane samples of the C-P (c) and C-P10E (d).
Polymers 16 00577 g002
Figure 3. Out-plane temperature profiles of the C-P (a) and C-P10E (b). The approximate heating rates of the rising parts of the curves are also reported.
Figure 3. Out-plane temperature profiles of the C-P (a) and C-P10E (b). The approximate heating rates of the rising parts of the curves are also reported.
Polymers 16 00577 g003
Figure 4. In-plane temperature profiles of the C-P (a) and C-P10E (b). The approximate heating rates of the rising parts of the curves are also reported.
Figure 4. In-plane temperature profiles of the C-P (a) and C-P10E (b). The approximate heating rates of the rising parts of the curves are also reported.
Polymers 16 00577 g004
Figure 5. Heat capacity of P10E after four dynamic scans compared to the C p of carbon (Black ■ from [58]) and of the material obtained after a sixth scan in TGA test up to 1000 °C (black continuous line).
Figure 5. Heat capacity of P10E after four dynamic scans compared to the C p of carbon (Black ■ from [58]) and of the material obtained after a sixth scan in TGA test up to 1000 °C (black continuous line).
Polymers 16 00577 g005
Figure 6. Heat capacity (a) and density (b) as a function of temperature for the developed CPCs as compared to state of art formulations. (Green ● from [32]; Yellow ▲ from [31]; Red and Blue x from [60]; Red and Blue ▲ from [61]; Red and Blue ⯁ from [62]; Red ■ from [63]).
Figure 6. Heat capacity (a) and density (b) as a function of temperature for the developed CPCs as compared to state of art formulations. (Green ● from [32]; Yellow ▲ from [31]; Red and Blue x from [60]; Red and Blue ▲ from [61]; Red and Blue ⯁ from [62]; Red ■ from [63]).
Polymers 16 00577 g006
Figure 7. Heat capacity (a) and density (b) as a function of temperature for the developed CPCs as compared to state-of-the-art formulations up to 3000 °C. Sample of a CPC in the virgin state (c) and at the end of the heat impulse imposed by the OAT test (d). The different depths at which the temperatures are measured are reported and then correlated with the heat capacity and density profiles. (Green ● from [32]; Yellow ▲ from [31]; Red and Blue x from [60]; Red and Blue ▲ from [61]; Red and Blue ⯁ from [62]; Red ■ from [63]).
Figure 7. Heat capacity (a) and density (b) as a function of temperature for the developed CPCs as compared to state-of-the-art formulations up to 3000 °C. Sample of a CPC in the virgin state (c) and at the end of the heat impulse imposed by the OAT test (d). The different depths at which the temperatures are measured are reported and then correlated with the heat capacity and density profiles. (Green ● from [32]; Yellow ▲ from [31]; Red and Blue x from [60]; Red and Blue ▲ from [61]; Red and Blue ⯁ from [62]; Red ■ from [63]).
Polymers 16 00577 g007
Figure 8. Mean values of heat capacity and density as a function of time of a C-P10E coupon tested with the OAT.
Figure 8. Mean values of heat capacity and density as a function of time of a C-P10E coupon tested with the OAT.
Polymers 16 00577 g008
Figure 9. Transient temperature distributions in a semi-infinite solid for constant surface temperature (a) and constant heat flux (b).
Figure 9. Transient temperature distributions in a semi-infinite solid for constant surface temperature (a) and constant heat flux (b).
Polymers 16 00577 g009
Figure 10. Charring phenomena within the CPC walls of a nozzle of a solid rocket motor (a). Post-burning appearance of a composite segmented submerged nozzle assembly made of CPCs and related phenomena of erosion and delamination (b). Appearance of an in-plane CPC sample tested with the OAT reproducing the delamination processes evidenced in (c), i.e., in the upper portion of the throat made of MX-4926 (image produced starting from an image taken from [66], credit NASA).
Figure 10. Charring phenomena within the CPC walls of a nozzle of a solid rocket motor (a). Post-burning appearance of a composite segmented submerged nozzle assembly made of CPCs and related phenomena of erosion and delamination (b). Appearance of an in-plane CPC sample tested with the OAT reproducing the delamination processes evidenced in (c), i.e., in the upper portion of the throat made of MX-4926 (image produced starting from an image taken from [66], credit NASA).
Polymers 16 00577 g010
Figure 11. Out-plane thermal diffusivity of the C-P and C-P10E formulations obtained under the assumption of constant surface temperature as compared to the corresponding data of the state-of-the-art counterparts (MX-4926 and FM-5014) (a). In-plane thermal diffusivity of the C-P and C-P10E formulations obtained under the assumption of constant surface temperature as compared to the corresponding data of the state-of-the-art counterparts (MX-4926 and FM-5014) (b). (Red and Blue x from [60]; Red and Blue ▲ from [61]; Red and Blue ⯁ from [62]; Red ■ from [63]).
Figure 11. Out-plane thermal diffusivity of the C-P and C-P10E formulations obtained under the assumption of constant surface temperature as compared to the corresponding data of the state-of-the-art counterparts (MX-4926 and FM-5014) (a). In-plane thermal diffusivity of the C-P and C-P10E formulations obtained under the assumption of constant surface temperature as compared to the corresponding data of the state-of-the-art counterparts (MX-4926 and FM-5014) (b). (Red and Blue x from [60]; Red and Blue ▲ from [61]; Red and Blue ⯁ from [62]; Red ■ from [63]).
Polymers 16 00577 g011
Figure 12. Out-plane thermal diffusivity of the C-P and C-P10E formulations obtained under the assumption of constant surface heat flux as compared to the corresponding data of the state-of-the-art counterparts (MX-4926 and FM-5014) (a). In-plane thermal diffusivity of the C-P and C-P10E formulations obtained under the assumption of constant surface heat flux as compared to the corresponding data of the state-of-the-art counterparts (MX-4926 and FM-5014) (b). (Red and Blue x from [60]; Red and Blue ▲ from [61]; Red and Blue ⯁ from [62]; Red ■ from [63]).
Figure 12. Out-plane thermal diffusivity of the C-P and C-P10E formulations obtained under the assumption of constant surface heat flux as compared to the corresponding data of the state-of-the-art counterparts (MX-4926 and FM-5014) (a). In-plane thermal diffusivity of the C-P and C-P10E formulations obtained under the assumption of constant surface heat flux as compared to the corresponding data of the state-of-the-art counterparts (MX-4926 and FM-5014) (b). (Red and Blue x from [60]; Red and Blue ▲ from [61]; Red and Blue ⯁ from [62]; Red ■ from [63]).
Polymers 16 00577 g012
Figure 13. Control surface for the transient heat transfer problem (a). Possible schematization of the heat losses by radiation and convection on the surface and lateral faces of the square shaped sample (b). Real temperature distribution on the lateral surfaces of a tested cubic CPC sample (c).
Figure 13. Control surface for the transient heat transfer problem (a). Possible schematization of the heat losses by radiation and convection on the surface and lateral faces of the square shaped sample (b). Real temperature distribution on the lateral surfaces of a tested cubic CPC sample (c).
Polymers 16 00577 g013
Figure 14. Behavior of the R(t) for the C-P10E out-plane and in-plane using the constant heat flux scenario.
Figure 14. Behavior of the R(t) for the C-P10E out-plane and in-plane using the constant heat flux scenario.
Polymers 16 00577 g014
Table 1. Summary of devices for determination of thermal conductivity (from [34]).
Table 1. Summary of devices for determination of thermal conductivity (from [34]).
MethodTemperature Range
(K)
Thermal Conductivity Range
(W/mK)
MaterialReference
Guarded Hot Plate80–8000.001–0.5PolymersASTM C177
Flynn et al. (2002) [35]
Heat Flow Meter80–8005–400SolidsASTM C518 [36]
Radial Heat Flow4–1000>0.1Solids, PowdersMaglic et al. (1992) [37]
Panel Test300–12500.05–15SolidsDaryabeigi (1999) [38]
Laser Flash80–2800(*)SolidsASTM E1461 [39]
Hot Wire300–1800<15Solids, LiquidsASTM C1113
Roder et al. (2000) [40]
(*) LFA measures thermal diffusivity, but it is able to retrieve the thermal conductivity through temperature-dependent heat capacity and density data.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Natali, M.; Torre, L.; Rallini, M. Evaluation of the Thermal Diffusivity of Carbon/Phenolic Composites (CPCs) through Oxy-Acetylene Torch (OAT) Test—Part 1: Experimental Characterization and Preliminary Validation. Polymers 2024, 16, 577. https://doi.org/10.3390/polym16050577

AMA Style

Natali M, Torre L, Rallini M. Evaluation of the Thermal Diffusivity of Carbon/Phenolic Composites (CPCs) through Oxy-Acetylene Torch (OAT) Test—Part 1: Experimental Characterization and Preliminary Validation. Polymers. 2024; 16(5):577. https://doi.org/10.3390/polym16050577

Chicago/Turabian Style

Natali, Maurizio, Luigi Torre, and Marco Rallini. 2024. "Evaluation of the Thermal Diffusivity of Carbon/Phenolic Composites (CPCs) through Oxy-Acetylene Torch (OAT) Test—Part 1: Experimental Characterization and Preliminary Validation" Polymers 16, no. 5: 577. https://doi.org/10.3390/polym16050577

APA Style

Natali, M., Torre, L., & Rallini, M. (2024). Evaluation of the Thermal Diffusivity of Carbon/Phenolic Composites (CPCs) through Oxy-Acetylene Torch (OAT) Test—Part 1: Experimental Characterization and Preliminary Validation. Polymers, 16(5), 577. https://doi.org/10.3390/polym16050577

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop