Ablation behavior and multi-physical field numerical simulation of ultra-lightweight quartz/phenolic composite
-
摘要: 基于超轻质石英/酚醛复合材料的烧蚀防热机制,建立了包含二氧化硅相变的复合材料烧蚀多物理场模型,预测了超轻质石英/酚醛复合材料的表背面温度、热解度、不同层厚度及孔隙压力分布。数值模拟通过计算获得了表面二氧化硅液态层厚度变化规律,模型预测的温度结果与烧蚀实验中测量结果吻合。根据各项传热模式的热流分析结果可见,对防/隔热机制影响最主要的是热辐射、热阻塞和二氧化硅气化。针对超轻质石英/酚醛复合材料的典型服役工况,采用0.5~2.5 MW/m2之间热流密度作为环境输入参数,研究了不同热流密度条件下超轻质石英/酚醛复合材料的烧蚀行为,结果表明:超轻质石英/酚醛复合材料表面烧蚀后退量随热流密度的增加而增加;热流密度小于1.5 MW/m2时表面液态层厚度随热流密度的增加而增加,热流密度大于1.5 MW/m2时表面液态层厚度基本保持不变。该模型为深入研究超轻质石英/酚醛复合材料的烧蚀机理提供一定的指导。
-
关键词:
- 超轻质石英/酚醛复合材料 /
- 高温烧蚀试验 /
- 多物理场 /
- 有限元法 /
- 二氧化硅相变
Abstract: Based on the ablation and heat protection mechanism of ultra-lightweight quartz/phenolic composite, a multi-physics field model for the ablation of composite containing silica phase transition was established. The surface and backside temperature, pyrolysis degree, different layer thicknesses and pore pressure distribution of ultra-lightweight quartz/phenolic composite were predicted. The numerical simulation obtained the variation law of the thickness of the liquid layer of surface silica through calculation, and the temperature results predicted by the model are consistent with the measurement results in the ablation experiment. According to the heat flow analysis results of various heat transfer modes, it can be seen that the most important factors affecting the prevention/insulation mechanism are thermal radiation, thermal blockage and silica gasification. Aiming at the typical application conditions of ultra-lightweight quartz/phenolic composite, the heat flux densities between 0.5 MW/m2 and 2.5 MW/m2 were used as environmental input parameters to study the ablation behavior of ultra-lightweight quartz/phenolic composite. The results show that the surface ablation retreat of ultra-light quartz/phenolic composite increases with the increase of heat flux density; When the heat flux density is less than 1.5 MW/m2, the thickness of the surface liquid layer increases with the increase of heat flux density. When the heat flux density is greater than 1.5 MW/m2, the thickness of the surface liquid layer remains unchanged. This model provides certain guidance for in-depth research on the ablation mechanism of ultra-lightweight quartz/phenolic composite. -
图 5 超轻质石英纤维/酚醛复合材料烧蚀后试验结果分析:(a)烧蚀后材料截面图;(b)表面液态层微观形貌;(c)炭化层微观形貌;(d)热解层微观形貌;(e)原始层微观形貌;(f)炭化层硅元素分布;(g)炭化层氧元素分布;(h)炭化层碳元素分布;(i)炭化层不同元素含量占比
Figure 5. Analysis of experimental results after ablation of ultra-lightweight quartz/phenolic composite: (a) Cross section view of composite after ablation; (b) Microstructure of surface liquid layer; (c) Microstructure of charring layer; (d) Microstructure of pyrolysis layer; (e) Microstructure of virgin layer; (f) Distribution of silicon element in charring layer; (g) Distribution of oxygen element in charring layer; (h) Distribution of carbon element in charring layer; (i) Proportion of different element contents in charring layer
图 6 超轻质石英纤维/酚醛复合材料在1.5 MW/m2热流密度条件下不同物理量演化规律:(a)不同深度温度曲线;(b)不同层分布规律;(c)不同深度热解度分布;(d)不同深度和不同层孔隙压力分布
Figure 6. Evolution of different physical quantities in ultra-lightweight quartz/phenolic composite under 1.5 MW/m2: (a) Temperature curves at different depths; (b) Different layer distribution patterns; (c) Distribution of pyrolysis degree at different depths; (d) Pore pressure distribution at different depths
图 11 不同热流密度条件下90 s时刻超轻质石英/酚醛复合材料不同层分层规律(加热时间:90 s):(a)0.5 MW/m2;(b)1.0 MW/m2;(c)1.5 MW/m2;(d)2.0 MW/m2;(e)2.5 MW/m2
Figure 11. Pattern of different layers of ultra-lightweight quartz /phenolic composite at 90 s under different heat flux density conditions (heating time: 90 s): (a) 0.5 MW/m2; (b) 1.0 MW/m2; (c) 1.5 MW/m2; (d) 2.0 MW/m2; (e) 2.5 MW/m2
表 1 超轻质石英/酚醛复合材料烧蚀后各层厚度实验与模拟结果(单位:mm)
Table 1. Experimental and simulation results of the thickness of each layer of ultra-lightweight quartz/phenolic composite material after ablation (Unit: mm)
Different layer Experiment Simulation Difference Ablation retreat 1.85 1.10 0.75 Liquid layer 0.97 0.90 0.07 Charring layer 12.50 10.03 2.47 Pyrolysis layer 7.29 8.76 1.47 Virgin layer 7.39 9.22 2.00 表 1 超轻质石英/酚醛复合材料性能参数
Table 1. Performance parameters of ultra-lightweight quartz/phenolic composite
Performance parameter Value Unit The virgin material density of silica as a solid state, $ {\rho _{{\text{vs}}}} $ 312.56 kg/m3 The charring material density of silica as a solid state, $ {\rho _{{\text{cs}}}} $ 250.17 kg/m3 The charring material density of silica as a liquid state, $ {\rho _{{\text{cl}}}} $ 550.24 kg/m3 The porosity of virgin material of silica as a solid state, $ {\phi _{{\text{vs}}}} $ 0.78 1 The porosity of charring material of silica as solid state, $ {\phi _{{\text{cs}}}} $ 0.88 1 The porosity of charring material of silica as liquid state, $ {\phi _{{\text{cl}}}} $ 0.45 1 The specific heat capacity of virgin material of silica as a solid state, $ {c_{{\text{pvs}}}} $ 1628.7 +0.976×TJ/(kg·K) The specific heat capacity of charring material of silica as a solid state, $ {c_{{\text{pcs}}}} $ 954.80+0.764×T J/(kg·K) The specific heat capacity of charring material of silica as a liquid state, $ {c_{{\text{pcl}}}} $ 1250.45 +0.544×TJ/(kg·K) The latent heat of silica phase transition, $ L $ 333 kJ/kg The pyrolysis gas specific heat capacity of phenolic, $ {c_{{\text{pg}}}} $ 2103.35 +1.05×TJ/(kg·K) The gas specific heat rate, $ \gamma $ 1.401-8.63×105×T 1 The permeability of virgin material of silica as a solid state, $ {\lambda _{{\text{vs}}}} $ 1.33×10−11 m2 The permeability of charring material of silica as a solid state, $ {\lambda _{{\text{cs}}}} $ 3.81×10−11 m2 The permeability of charring material of silica as a liquid state, $ {\lambda _{{\text{cl}}}} $ 0.84×10−11 m2 The thermal conductivity of virgin material of silica as a solid state, $ {k_{{\text{vs}}}} $ 0.05+4.05×10−5 W/(m·K) The thermal conductivity of charring material of silica as a solid state, $ {k_{{\text{cs}}}} $ 0.316+9.76×10−5 W/(m·K) The thermal conductivity of virgin material of silica as a liquid state, $ {k_{{\text{cl}}}} $ 0.15+4.05×10−5 W/(m·K) The thermal conductivity of pyrolysis gas, $ {k_{\text{g}}} $ 0.0246 W/(m·K) The emissivity of virgin material of silica as a solid state, $ {\varepsilon _{{\text{vs}}}} $ 0.80 1 The emissivity of charring material of silica as a solid state, $ {\varepsilon _{{\text{cs}}}} $ 0.90 1 The emissivity of virgin material of silica as a liquid state, $ {\varepsilon _{{\text{cl}}}} $ 0.95 1 The phenolic pyrolysis pre exponential factor, $ A $ 3.3×105 1/s The activation energy of phenolic pyrolysis, $ E $ 74.8 kJ/mol The pyrolysis gas viscosity coefficient, $ {\bar \mu _{\text{g}}} $ 14.8×10−5+25×10−8×T Pa·s The pyrolysis gas average molar mass, $ {\bar M_{\text{g}}} $ 0.01835 kg/mol The latent heat of silica gasification, $ {H_{\text{v}}} $ 12868 kJ/kg The standard atmospheric pressure, $ {P_0} $ 1×105 Pa The gasification temperature of silica, $ {H_{\text{b}}} $ 2500 ℃ The density of silica, $ {\rho _{\text{f}}} $ 2200 kg/m3 The volume ablation rate of silica, $ {J_{\text{f}}} $ 6.6×103 kg/(m3·s−1) The thermal conductivity of silica, $ {k_{\text{f}}} $ 1.3 W/(m·K) The specific heat capacity of silica, $ {C_{\text{f}}} $ 650 J/(kg·K) The tensile strength of composite, $ {\sigma _{{\text{mT}}}} $ 220 MPa The mechanical erosion activation energy, $ {E_{\text{A}}} $ 9.8×103 J/mol The average molar mass of air, $ {M_{\text{a}}} $ 0.02896 kg/mol The boundary layer air recovery enthalpy, $ {h_{{\text{re}}}} $ 35 MJ/kg The gas constant, $ R $ 8.314 J/(mol·K) The Boltzmann constant, $ m $ 1.38×10−23 J/K The gasification coefficient of silica, $ \beta $ 0.84 1 The molecular weight of silica, $ m $ 60×1.674×10−27 kg The average specific heat of surface ablative materials, $ {c_{\text{p}}} $ 700 J/(kg·K) The content of phenolic, $ {f_{\text{p}}} $ 0.455 1 The carbon residue rate of phenolic, $ \xi $ 0.564 1 The relative molar mass of carbon atoms, $ {M_{\text{C}}} $ 0.012 kg/mol The relative molar mass of oxygen atoms, $ {M_{\text{O}}} $ 0.016 kg/mol -
[1] 姜贵庆, 刘连元. 高速气流传热与烧蚀热防护[M]. 北京: 国防工业出版社, 2003: 52-93.JIANG Guiqing, LIU Lianyuan. Heat transfer of hypersonic gas and ablation thermal protection[M]. Beijing: National Defense Industry Press, 2003: 52-93(in Chinese). [2] 张睿, 祝文祥, 张澳, 等. 烧蚀型热防护系统概率设计与可靠性评估方法研究[J]. 计算力学学报, 2023, 21(8): 4708.ZHANG Rui, ZHU Wenxiang, ZHANG Ao, et al. Probabilistic design and reliability assessment methods for ablative thermal protection systems[J]. Chinese Journal of Computational Mechanics, 2023, 21(8): 4708(in Chinese). [3] 韩杰才, 洪长青, 张幸红, 等. 新型超轻质热防护复合材料的研究进展[J]. 载人航天, 2015, 21(4): 315-321. doi: 10.3969/j.issn.1674-5825.2015.04.001Han Jiecai, Hong Changqing, Zhang Xinghong, et al. Research progress of novel lightweight thermal protection composites[J]. Manned Spaceflight, 2015, 21(4): 315-321(in Chinese). doi: 10.3969/j.issn.1674-5825.2015.04.001 [4] JIN X Y, LIU C, HUANG H, et al. Multiscale, elastic, and low-density carbon fibre/siliconoxycarbide-phenolic interpenetrating aerogel nanocomposite for ablative thermal protection[J]. Composites Part B: Engineering, 2022, 217(1): 109100. [5] Yan X J, Huang H, Fan Z L, et al. Assessment of a 3D ablation material response model for ultra-lightweight quartz/phenolic composite[J]. Polymer Composites, 2022, 43(4): 8341-8355. [6] Liang H R, Li W J, Wang T K, et al. Optimal design of three-dimensional thermal protection structure considering orthotropic properties of woven composites based on Micro-CT image[J]. International Journal of Thermal Sciences, 2023, 194(7): 108579. [7] 张军, 许阳阳, 张运法, 等. 石英灯辐射加热条件下低密度碳/酚醛复合材料高温响应及分析[J]. 装备环境工程, 2020, 17(1): 51-57.Zhang Jun, Xu Yangyang, Zhang Fayun, et al. High temperature response and analysis of low density carbon fiber/phenolic composites under quartz lamp radiation heating[J]. Equipment Environmental Engineering, 2020, 17(1): 51-57(in Chinese). [8] Pedro G S P, Homero d P e S, Cristian C P R, et al. The effect of environment in carbon-phenolic composite ablation[J]. Journal of the Brazilian Society of Mechanical Sciences and Engineering, 2021, 500(43): 1-16. [9] Pan Y W, Jin X Y, Wang H B, et al. Nano-TiO2 coated needle carbon fiber reinforced phenolic aerogel composite with low density, excellent heat-insulating and infrared radiation shielding performance[J]. Journal of Materials Science & Technology, 2023, 152(12): 181-189. [10] Wang H B, Quan X , Yin L , et al. Lightweight quartz fiber fabric reinforced phenolic aerogel with surface densified and graded structure for high temperature thermal protection[J]. Composites Part A: Applied Science and Manufacturing, 2022, 109798. [11] Cai H X, Niu B, Qian Zh, et al. Mechanical, thermal insulation, and ablation behaviors of needle-punched fabric reinforced nanoporous phenolic composites: The role of anisotropic microstructure[J]. Composites Science and Technology, 2024, 245(8): 110325. [12] 陈海龙, 方国东, 李林杰, 等. 高硅氧/酚醛复合材料热-力-化学多物理场耦合计算[J]. 复合材料学报, 2014, 81(8): 533-540.Chen Hailong, Fang Guodong, Li Linjie, et al. Multi-physical field coupling calculation of silica/phenolic composites under thermal-mechanical-chemical condition[J]. Acta Materiae Compositae Sinica, 2014, 81(8): 533-540(in Chinese). [13] 李林杰, 方国东, 易法军, 等. 高硅氧/酚醛复合材料热变形实验测试及表面烧蚀形貌分析[J]. 固体火箭技术, 2015, 38(3): 445-450.Li Linjie, Fang Guodong, Yi Fajun, et al. Thermal deformation test and surface ablative morphology analysis of silica / phenolic composites[J]. Journal of Solid Rocket Technology, 2015, 38(3): 445-450(in Chinese). [14] Shi S B, Li L J, Liang J, et al. Surface and volumetric ablation behaviors of SiFRP composites at high heating rates for thermal protection applications[J]. International Journal of Heat and Mass Transfer, 2016, 102(11): 1190-1198. [15] 宋若康, 张梦珊, 戴珍, 等. 烧蚀型防热/吸波多功能一体化复合材料的制备及性能[J]. 复合材料学报, 2024, 41(1): 252-261.Song Ruokang, Zhang Mengshan, Dai Zhen, et al. Preparation and properties of multi-functional composite integrated with heat-shielding and radar-absorbing[J]. Acta Materiae Compositae Sinica, 2024, 41(1): 252-261(in Chinese). [16] 白炳越, 刘丹, 杨进军, 等. 酚醛树脂浸渍石英纤维复合材料的热解行为[J]. 灭火剂与助燃材料, 2023, 42(9): 1275-1279.Bai Bingyue, Liu Dan, Yang Jinjun. Thermal decomposition behavior of phenolic resin impregnated quartz fiber composite materials[J]. Fire extinguishing agents and combustion materials, 2023, 42(9): 1275-1279(in Chinese). [17] Huang H, Yan X J, Jin X Y, et al. Flexible and interlocked quartz fibre reinforced dual polyimide network for high-temperature thermal protection[J]. Journal of Materials Chemistry A, 2023, 11(3): 9931-9941. [18] 王湘阳, 年永乐, 刘娜, 等. 考虑C-SiO2反应的新型硅基材料烧蚀分析模型[J]. 化工学报, 2021, 72(6): 3270-3277.Wang Xiangyun, Nian Yongle, Liu Na, et al. Novel ablation model of silica-reinforced composites considering C-SiO2 reaction[J]. CIESC Journal, 2021, 72(6): 3270-3277(in Chinese). [19] 许善成. 切向空气流下连续激光辐照玻璃纤维树脂基复合材料烧蚀特性的研究[D]. 南京: 南京理工大学, 2019.Xu Shancheng. Study on ablation characteristics of glass fiber reinforced resin matrix composites by CW laser irradiation under tangential air flow[D]. Nanjing: Nanjing University of Science & Technology, 2019(in Chinese). [20] 智伟. 柔性防热材料辐射热流下烧蚀行为的数值计算[D]. 天津: 天津工业大学, 2019.Zhi Wei. Numerical calculation of ablation behavior of flexible heat-resistant materials under radiation heat flux[D]. Tianjin: Tiangong University, 2019(in Chinese). [21] Humberto A M, Sonia F C Silva, Edison B. Simulation of ablation in composite via an interface tracking method[J]. Journal of the Brazilian Society of Mechanical Sciences and Engineering, 2017, 39: 1453-1467. [22] 黄鹏, 郑振荣, 毛科铸, 等. 热阻塞效应在有机硅树脂-碳纤织物复合材料烧蚀防热中的作用[J]. 复合材料学报, 2021, 38(9): 3045-3055.Huang Peng, Zhen Zhengrong, Mao Kezhu, et al. Effect of heat blockage on ablative thermal protection of silicone resin-carbon fiber fabrics[J]. Acta Materiae Compositae Sinica, 2021, 38(9): 3045-3055(in Chinese). [23] 万涛. 超轻质碳/酚醛复合材料的烧蚀行为和仿真预报研究[D]. 哈尔滨: 哈尔滨工业大学, 2017.Wang Tao. Research on ablation behavior and simulation prediction of ultra-lightweight carbon phenolic composite[D]. Harbin: Harbin institute of Technology, 2017(in Chinese). [24] 时圣波. 高硅氧/酚醛复合材料的烧蚀机理及热-力学性能研究[D]. 哈尔滨: 哈尔滨工业大学, 2013.Shi Shengbo. Ablation Mechanism and Thermo-Mechanical Behavior of Silica/Phenolic Composites[D]. Harbin: Harbin institute of Technology, 2013(in Chinese). -