3.1. Fabrication, structure, and morphology of SQRP
A two-step impregnation-curing process was employed in the fabrication of the SQRPs (Fig. 1 (a)). For reinforcement, a quartz felt with high porosity was chosen. Initially, the reinforcement underwent a shallow impregnation with a pre-mixed CR solution through capillary adsorption, followed by curing to form the surface region matrix. This matrix is crucial for withstanding ablative recession and maintaining shape integrity. Subsequently, the reinforcement underwent a secondary vacuum impregnation using a meticulously formulated precursor solution, resulting in the formation of an internal PRa matrix through polymer-induced phase separation. This internal matrix is instrumental in providing thermal insulation. Macroscopically, the thickness of the surface region matrix is approximately 3 mm (Fig. 1 (b)), and the integrated reinforcement can ensure the absence of significant delamination or cracking at the dual-matrix junction. In contrast to conventional ATPMs, the layered dual-matrix of the SQRPs serves different thermal protection functions, specifically, ablation resistance and thermal insulation (Fig. 1 (c)), as will be elaborated in the following sections. Figure 1 (d-g) illustrate the microscopic structure of the SQRP. The enriched CR matrix effectively occupies the surface region, with tightly bonded micron-sized ceramic particles (Fig. 1 (d)). The interphase interface of the dual matrix integrated well and aerogel particles formed directly on the CR at the junction (Fig. 1 (e)). In sharp contrast to the surface, the interior of the composite contains pores resulting from PRa particles clustering together and forming a highly porous network (Fig. 1 (f)). Further magnified SEM images reveal the presence of nanosized particles.
3.2. Chemical structure and thermal stability of dual matrix
The methyl-phenyl-silicone resin in the CR matrix formed -Si-O-Si- bonds through dehydro-condensation, while vacant hydroxyl groups created weak hydrogen bonds with ceramic particles, facilitating the homogeneous dispersion. Within the internal aerogel matrix, a highly cross-linked phenolic network was formed through various chemical bonds, such as -CH2- and -Si-O-C- [30]. This polymerization mechanism had been validated through XPS and FTIR. In the XPS spectra (C 1s), the characteristic peak at 284.2 eV for CR corresponded to the -C-Si- bond connecting the methyl and benzene [31], and the peak at 102.5 eV (Si 2p) corresponded to the -Si-O-Si- bond along the molecular main chain [32, 33] (Fig. 2 (b-c)). For PRa, the characteristic peaks at 285.0 and 285.4 eV (C 1s) corresponded to the methylene and -C-O- bonds, and the peak at 103.9 eV (Si 2p) represented the -Si-O-C- bond formed by the reaction between silane and phenol [32, 34, 35] (Fig. 2 (d-e)). FTIR spectrum Fig. 2 (f) reveals the presence of antisymmetric stretching vibration peaks associated with -Si-O-Si- at 1027 and 1135 cm− 1 for CR. The double peaks provided evidence of the formation of long molecular chains through polycondensation [36]. The spectrum also shows additional phenolic chain linkages, including -CH2- stretching vibration peaks (2924 and 1466 cm− 1), -C-N-C deformation vibration peaks (2228 and 1272 cm− 1), and methylene ether (-CH2-O-CH2-) stretching vibration peaks (1109 and 968 cm− 1) [37]. The thermal stability of CR with ceramic particles under an argon atmosphere surpassed that of the original silicone resin, with a residual weight exceeding 90% at 1100°C. In contrast, the PRa matrix exhibited a 63.1% residual weight (see Fig. 2 (g-h)). Weight loss primarily occurred due to the release of solvents, small terminal molecules, and molecular chain breakdown at high-temperature, with significant losses evident around 550°C. In the case of SiR, the loss was mainly attributed to -Si-O- and -C = C- bond breakage, while for PRa, it resulted from methylene and -C = C- bond breakage. Figure 2 (i-k) illustrates the phases of CR after high-temperature static oxidation treatment. Below 1400°C, the silicone resin was cleaved into amorphous particles, while ZrB2 particles oxidized to ZrO2 and SiC particles oxidized to SiO2 [38, 39]. Importantly, the higher the temperature, the less oxide was generated (Fig. 2 (j)), as the rapid oxidation of ceramic particles at elevated temperatures formed a stable oxide protective layer on their surfaces, thereby reducing ablation recession.
3.3. Evaluation of mechanical strength and thermal insulation
The mechanical properties of the composites were evaluated utilizing a universal mechanical testing machine. Initially, the tensile and flexural response stresses in the XY direction of the SQRPs exhibited linear growth, followed by a rapid decline upon reaching critical points, a behavior consistent with brittle fracture [40] (Fig. 3 (a, c)). During fracture, stress concentrations were particularly pronounced on the densely structured CR matrix at the surface region. The matrix disintegrated and detached from the fibres, leading to failure of the composites. The variations in ceramic particle types had minimal influence on the mechanical strength. The tensile strengths of the composites exceeded 4.62 MPa, and the flexural strengths reached 19.75 MPa (Fig. 3 (b, d)). On the contrary, when subjected to Z-direction compression, the porous PRa matrix experienced compaction rather than fracturing. Consequently, the response stress exhibited a consistent increase (Fig. 3 (e)). The point of offset yield stress at 50% strain defined the maximum compressive strength of the composites, which exceeded 3.07 MPa.
The thermal insulation properties of the dual-matrix composites primarily rely on the internal porous structure. Both surface and internal region thermal conductivities were examined, alongside the high-temperature conductivities within the critical interior. The interior exhibited a low thermal conductivity of only 0.04 at room temperature, which increased to 0.104 at 500°C (Fig. 4 (a)). Heat transfer within the porous interior is governed by the necessity to traverse an intricate network of air-filled pores and an infinitely extended aerogel wall, considerably inhibiting solid-phase and gas-phase conduction [41]. Notably, the pore sizes of both the PRa matrix and the internal composites were predominantly in the range of 30 to 40 nm (Fig. 4 (b-c)). The N2 adsorption-desorption isotherms conformed to Type IV isotherms and H3 hysteresis lines in accordance with the IUPAC classification of particle stacking structures (detailed pore data refer to Table S3). The dynamic thermal insulation performance of the composites was evaluated using a butane flame ablation, and the cloud diagram illustrating the backside response is presented in Fig. 4 (f). SQRPs featuring various CR matrices exhibited consistent backside response temperatures. During the initial phase, the temperature remained relatively stable, signifying minimal heat transfer to the backside. Subsequently, the pyrolysis of the surface matrix led to a slight increase in response temperature, followed by a second slight increase caused by the pyrolysis of the aerogel matrix. SQRPs demonstrated effective dynamic thermal insulation, maintaining backside temperatures at only 116°C after 600 seconds (as shown in Fig. 4 (g)) and reaching an equilibrium temperature of approximately 450°C after extended exposure to the butane flame (detail in supporting information).
3.4. Evaluation of ablative thermal protection
The thermal protection performance of SQRPs was assessed using oxyacetylene flame ablation, generating heat flows of 1.5 MW (~ 1700°C) and 3.6 MW (~ 2500°C). As depicted in Fig. 5 (b), the ceramicization of the CR matrix resulted in robust surfaces without edge shrinkage. Additionally, the complete vaporization of the pyrolytic carbon exposed voids originating from vertically arranged fibre bundles. Following 90 seconds of exposure to 1.5 MW ablation and 60 seconds of exposure to 3.6 MW ablation respectively, the backside temperatures of the specimens measured approximately 48°C and 75°C (measured depths of 25 mm), as shown in Fig. 5 (c-d), which can be attributed to the exceptional dynamic ablative insulating performance of composites. Notably, the S2QRP displayed the lowest ablation recession, with a weight loss rate of 9.22 mg/s and a linear ablation setback rate of 0.11 µm/s during 1.5 MW heat flow ablation. At 3.6 MW, the corresponding figures were 18.4 mg/s and 6 µm/s, respectively (Fig. 5 (e-f)). Detailed abaltion data are shown in Table S3 and S4. This remarkable performance derived from the high concentration of nanoscale ceramic particles present in the S2QRP, which significantly enhanced the resistance to erosion caused by heat and oxygen under flame. Furthermore, the surface-enriched CR matrix resulted in a slight increase of density, effectively reducing the ablative recession of the composites. The pyrolysis products originating from the ceramic particles and silicone resin formed an anti-ablation shielding layer when exposed to the flame. This shielding layer served to protect the internal porous matrix and fibres, thus ensuring the stability of the aerodynamic profile. In addition, the internal phenolic matrix underwent partial pyrolysis, providing thermal protection through mass injection and pyrolysis gas emission. As materials designed for ablative thermal protection, SQRPs offer a competitive combination of low density and exceptionally low linear ablative setback, as shown in Fig. 5 (g) [19, 21–26, 37, 42–48]. Slicing the ablated specimen (1.5 MW ablated S2QRP, 3 mm thick per layer) in the thickness direction, defined as ablated surface region, charred region, pyrolyzed region, and virgin region [6, 49] (Fig. 5 (h)). It can be observed from the FTIR spectra, the silicone resin at the ablated surface underwent pyrolysis, with the double peak of the long-chain Si-O transforming into a single peak at 1388 cm− 1. Additionally, the characteristic peaks of the residual -CH3 were observed at 1388 and 838 cm− 1. The pyrolysis of the internal phenolic aerogel matrix exhibited regressive characteristics, with the layer (L2) closest to the ablated surface undergoing complete charring and the disappearance of the characteristic peak of the methylene group connecting the benzene ring bonds at 1466 cm− 1. The layers of the pyrolysis region (L2-L6) did not fully char, with the characteristic peak of the methylene ether -CH2-O-CH2- weakly associated with the benzene ring being observed at 968 cm− 1. The characteristic peaks of the layers in the virgin region (L7-L9) remained well preserved and essentially free from pyrolysis. Similarly, the density of the ablated layers exhibited a regressive recession corresponding to the degree of matrix pyrolysis. The ablated surface of the composite contained the expected oxides of the ceramic particles, along with melt crystallization of the quartz fibres.
The micromorphology of the dual-matrix in the ablated S2QRP was further examined, as depicted in Fig. 6. It is evident that the CR matrix formed an ablative shielding layer on the surface, while the silicone resin underwent pyrolysis and melt, forming a liquid film that encapsulated the ceramic particles and oxides. This ceramic liquid film effectively preventing the intrusion of heat and oxygen, thereby averting significant ablative degradation of the composite. Furthermore, the ceramic particles, resembling islands, were embedded within the liquid film of the silicon phase, which restricted its flow Fig. 6 (b)). Upon inspection, the high-magnification SEM image (Fig. 6 (e)) reveals that the oxides on the surface of the ceramic particles had fused together. In contrast, the oxides in proximity to the charred region exhibited lower response temperatures, causing to melt into small droplets, as seen in Fig. 6 (f). Figure 6 (i-k) provide the micro-morphology of the charred region, the pyrolysis region, and the virgin region, respectively. In the charred region, extensive pyrolysis of the phenolic aerogel resulted in enlarged pores and exposed fibres. Conversely, the aerogel in the pyrolysis and virgin regions remained intact. The high-magnification SEM image clearly shows unchanged aerogel particles.
3.5. Simulation analysis under 1.5MW heat flow ablation
To delve deeper into the ablation process within the composite, a simulation analysis was conducted, and the results are presented in Fig. 7. The simulation model encompasses a dual-matrix composite (S2QRP) and a graphite sleeve (Fig. 7 (a)), with the temperature collection point positioned at 25 mm depth. Figure 7 (b) demonstrates a consistency between experimental and simulated temperature curves, affirming the precision of model. With prolonged simulated ablation time, the organic resin matrix gradually pyrolyzed and permeated to the bottom (Fig. 7 (c)). Figure 7 (d) illustrates temperature curves at various locations inside the composite, showcasing a rapid rise in surface temperature (above 1700°C) and a gradual increase in internal response temperature due to heat conduction. Notably, the backside temperature at (0,0) surpasses that at (0,5), primarily attributed to the high thermal conductivity of the external graphite sleeve. Figure 7 (e) tracks the pyrolysis rate with time at different depths. At 30 s, the thickness of the fully pyrolyzed region reaches 5 mm, and even after heating ceases at 90 s, the inner phenolic matrix continues to pyrolyze due to accumulated heat. Figure 7 (f) depicts the evolution of pore pressure at various depths over time. One second into ablation witnesses extensive pyrolysis of the surface-enriched CR matrix, resulting in a rapid surge in gas pressure on the ablation surface, peaking at 10,500 Pascals (Pa). Throughout heating, pore pressure within the composite consistently remains above 2000 Pa. Conversely, during cooling, internal pressure drops to approximately 0 Pa, indicating rapid diffusion of pyrolysis gases through the porous structure into the surrounding air. Moving to Fig. 7 (g), a cloud diagram illustrates temperature distribution from left to right at six time points: 1, 45, 90, 91, 135, and 180 seconds. During the heating phase, the high-temperature region concentrates in the upper part of the composite. As the cooling phase progresses, a "heat kernel" region forms as heat continues to transfer inward, but the temperature peak in this region significantly decreases. By 180 seconds, the peak internal temperature of the material stands at 260°C. Figure 7 (h) offers a cloud diagram of pyrolysis ratio distribution at the same time intervals. During heating, the resin at the edge undergoes faster pyrolysis compared to the central region, primarily due to exposure to environmental heat flow. In the cooling phase, heat continues to conduct inward, causing further pyrolysis of the resin. These findings align with the results presented in Fig. 7 (e). Finally, Fig. 7 (i) displays a cloud diagram of pyrolysis gas flow field distribution at the same time intervals. Notably, in the first second of ablation, the flow velocity of pyrolysis gases is substantial, primarily confined to the vicinity of the ablation surface. This phenomenon results from the rapid pyrolysis of resin on the surface induced by environmental heat. As heat propagates inward, the rate of pyrolysis gradually decreases.
3.6. Evaluation of thermal protection under repeated ablation
The S2QRP underwent a series of repeated ablation tests consisting of 10 cycles, each lasting 300 seconds with a butane flame. Similar to the behavior observed during oxygen-acetylene ablation, the pyrolysis of the composite showed a regression. This behavior created distinct regions within the composite, including the ablated surface, charred region, pyrolyzed region, and virgin area. In the early stages of the ablation cycle, the pyrolysis of the composite was incomplete, leading to a gradual increase in weight loss. After four cycles, the pyrolysis was essentially complete, resulting in only marginal additional weight loss. As a result, the composite retained over 86% of its initial weight. Due to the ablation-resistant CR matrix on the surface, the composite showed slight linear and radial setbacks across repeated ablation cycles. Overall, the composite maintained effective insulation, with backside temperatures remaining below 80°C after 10 ablation cycles. Nevertheless, there was a slight increase in temperature as the number of ablation cycles increased. The ablative setback of ATPMs significantly limits reusability. To assess the surface contour after ten ablation cycles, laser microscopy was employed to observe the flame ablation center, as shown in Fig. 8 (c-f). Figure 8 (d) displays the surface contour at the blue markers on the ablated surface, revealing the presence of high peaks. A height distribution of the surface was also obtained through laser intensity, and 3D surface reconstruction was performed. The combination of Fig. 8 (e) and (f) highlights that the surface peaks consist of warped fibers and ceramicized resin. The arithmetic mean curvature (Spc) of the peak peaks in the collected region measures 3167.231 1/mm, with an arithmetic mean height (Sa) of 8.293 and a maximum height (Sz) of 166.381 µm, respectively. Additionally, observable concave valleys originate from the exposed needled fibre bundles. Overall, the ablation center area of the composite was flat and there was no visible ablation setback.