NUMERICAL SIMULATION OF MICRO-ABLATION BEHAVIOR FOR CARBON FIBER REINFORCED COMPOSITES1)

doi:10.6052/0459-1879-18-351

Abstract

Abstract:

Carbon fiber reinforced composites are widely used in ablative thermal protection systems (TPS). The microstructures of the composites are greatly related with the ablation behavior. Thus, the study of micro-ablation mechanisms of the composites is significant for the material design and manufacturing. A mathematic micro-ablation model is developed by using the finite-volume-method (FVM) combing with the piecewise linear interface calculation (PLIC) method. Comparing with the analytical results, the numerical method is validated by calculating the ablation morphologies for single fiber embedded in carbon matrix. The effect of carbon fiber inclination on the microscopic ablation behavior is studied for the unidirectional carbon reinforced composites with different carbon fiber inclination. The results are found that the ablation morphology for carbon fiber is bamboo shoot shape if the oxygen diffusion rate is far greater than that of carbon oxygen reaction when the oxidation resistance of carbon fibers is stronger than that of matrix. If the carbon oxygen reaction rate is far greater than that of oxygen diffusion, carbon fiber and matrix will be ablated at the same rate. When the oxidation resistance of carbon fiber is stronger than that of the matrix, the inclination angle of the carbon fiber has a great influence on the micro-ablation behavior of the material. On the contrary, the effect is not significant.

Key words: micro-ablation, carbon fibre reinforced composites, numerical simulation, fiber inclination angle

CLC Number:

O346.1

Wei Li, Guodong Fang, Weijie Li, Bing Wang, Jun Liang. NUMERICAL SIMULATION OF MICRO-ABLATION BEHAVIOR FOR CARBON FIBER REINFORCED COMPOSITES¹⁾[J]. Chinese Journal of Theoretical and Applied Mechanics, 2019, 51(3): 835-844.

Figures/Tables 19

Fig. 1 Illustration of phase contact area

Fig. 2 Illustration of loosely coupled algorithm

Fig. 3 Sketch map of the validation of loose coupling algorithm

Table 1 Interface position

Fig. 4 Oxygen concentration distribution

Fig. 5 Schematic diagram of mode

Fig. 6 Comparison of the morphologies from simulation results and analytical results

Fig. 7 Ablation morphologies of single fiber at different Sh numbers ($A=5$)

Fig. 8 Typical micro-ablation morphologies

Fig. 9 Schematic diagram of material with different fiber inclination angle

Fig. 10 Micro-ablation morphologies at different times ($A=5,Sh=0.01$)

Fig. 11 Mass loss rate of materials at different fiber inclination angles ($A=5, Sh=0.01$)

Fig. 12 Micro-ablation morphologies at different times ($A=5, Sh=1$)

Fig. 13 Mass loss rate of materials at different fiber inclination angles ($A=5$, $Sh=1$)

Fig. 14 Micro-ablation morphologies at different times ($A=5$, $Sh=100$)

Fig. 15 Micro-ablation morphologies at different times ($A=0.2, Sh=0.01$)

Fig. 16 Mass loss rate of materials at different fiber inclination angles ($A=0.2$, $Sh =0.01$)

Fig. 17 Micro-ablation morphologies at different times ($A=0.2$, $Sh =1$)

Fig. 18 Mass loss rate of materials at different fiber inclination angles ($A=0.2$, $Sh =1$)

References 33

[1]	叶友达, 张涵信, 蒋勤学等. 近空间高超声速飞行器气动特性研究的若干关键问题. 力学学报, 2018, 50(6): 1292-1310 DOI URL
	(Ye Youda, Zhang Hanxin, Jiang Qinxue, et al.Some key problems in the study of aerodynamic characteristics of near-space hypersonic vehicles. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(6): 1292-1310 (in Chinese)) DOI URL
[2]	赵金山, 张志刚, 石义雷等. 高超声速飞行器气动热关联换算方法研究. 力学学报, 2018, 50(5): 1235-1243
	( Zhao Jinshan, Zhang Zhigang, Shi Yilei, et al.Research on the conversion method of aeroheating environment of hypersonic vehicle. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(5): 1235-1243 (in Chinese))
[3]	费庆国, 姜东, 陈素芳等. 高温下编织复合材料热相关参数识别方法研究. 力学学报, 2018, 50(3): 497-507
	(Fei Qingguo, Jiang Dong, Chen Sufang, et al.Thermal-related parameter identification of braided composites at high temperature. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(3): 497-507 (in Chinese))
[4]	Han JC, He XD, Du SY.Oxidation and ablation of 3D carbon-carbon composite at up to 3000 $^\circ$C. Carbon, 1995, 33(4): 473-478 DOI URL
[5]	Paglia L, Tirill J, Marra F, et al.Carbon-phenolic ablative materials for re-entry space vehicles: Plasma wind tunnel test and finite element modeling. Materials & Design, 2016, 90(11): 70-80 DOI URL
[6]	Tran H, Johnson C, Rasky AD, et al.Phenolic impregnated carbon ablators (PICA) for discovery class missions //Proceedings of the Thermophysics Conference, F, 2013 DOI URL
[7]	Dimitrienko YI.Thermomechanics of Composite Structures under High Temperatures. Netherlands: Springer, 2016 DOI URL
[8]	Edquist KT, Hollis BR, Johnston CO, et al.Mars science laboratory heatshield aerothermodynamics: Design and reconstruction // Proceedings of the AIAA Thermophysics Conference, F, 2014 DOI URL
[9]	Kumar S, Kushwaha J, Mondal S, et al.Fabrication and ablation testing of 4D C/C composite at 10 MW/m$^{2}$ heat flux under a plasma arc heater. Materials Science and Engineering: $A$, 2013, 566: 102-111
[10]	Moyer CB, Wool MR.Aerotherm Charring Material Thermal Response and Ablation Program, Version 3. Volume 1. Program Description and Sample Problems. Aerotherm Charring Material Thermal Response & Ablation Program Version, 1970
[11]	Chen YK, Milos F, Rasky DJ.Fully implicit ablation and thermal response program for spacecraft heatshield analysis. Journal of Spacecraft & Rockets, 1997, 36(3): 475-83 DOI URL
[12]	Lachaud J, Mansour NN.Porous-material analysis toolbox based on OpenFOAM and applications. Journal of Thermophysics & Heat Transfer, 2014, 28(2): 191-202. DOI URL
[13]	Li W, Huang H, Tian Y, et al.Nonlinear analysis on thermal behavior of charring materials with surface ablation. International Journal of Heat & Mass Transfer, 2015, 84(2): 45-52 DOI URL
[14]	中国人民解放军总装备部军事训练教材工作委员会. 高超声速气动热和热防护. 北京: 国防工业出版社, 2003
	(Working Committee on Military Training Textbooks of the General Armament Department of the Chinese People's Liberation Army. Hypersonic Aerodynamic Thermal and Thermal Protection. Beijing: National Defense Industry Press, 2003 (in Chinese))
[15]	姜贵庆, 刘连元. 高速气流传热与烧蚀热防护. 北京: 国防工业出版社, 2003
	(Jiang Guiqing, Liu Lianyuan.High Speed Airstream Heat Transfer and Ablation Thermal Protection. Beijing: National Defense Industry Press, 2003 (in Chinese))
[16]	姚承照, 胡宝刚, 冯志海等. 三维整体编织碳/酚醛复合材料烧蚀表面状态测试与分析. 宇航材料工艺, 2001, 31(6): 72-76 DOI URL
	(Yao Chengzhao, Hu Baogang, Feng Zhihai, et al.Test and analysis of a blation surface of 3D braided carbon/phenolic resin composites. Aerospace Materials Technology, 2001, 31(6): 72-76 (in Chinese)) DOI URL
[17]	姜贵庆, 李鸿权. 热化学烧蚀理论及其应用. 高超声速前沿问题研讨班, 2002
	(Jiang Guiqing, Li Hongquan.Thermochemical ablation theory and its application. Seminar on Advanced Hypersonic Problems, 2002 (in Chinese))
[18]	Duffa G, Vignoles GL, Goyhénèche JM, et al.Ablation of carbon-based materials: Investigation of roughness set-up from heterogeneous reactions. International Journal of Heat & Mass Transfer, 2005, 48(16): 3387-3401 DOI URL
[19]	Aspa Y, Quintard M, Plazanet F, et al.Ablation of Carbon/Carbon Composites: Direct Numerical Simulation and Effective Behavior. John & Sons, Inc, 2008 DOI URL
[20]	Vignoles GL, Aspa Y, Quintard M.Modelling of carbon--carbon composite ablation in rocket nozzles. Composites Science & Technology, 2010, 70(9): 1303-1311 DOI URL
[21]	Lachaud J, Bertrand N, Vignoles GL, et al.A theoretical/experimental approach to the intrinsic oxidation reactivities of C/C composites and of their components. Carbon, 2007, 45(14): 2768-2776 DOI URL
[22]	Lachaud J, Aspa Y, Vignoles GL.Analytical modeling of the steady state ablation of a 3D C/C composite. International Journal of Heat & Mass Transfer, 2008, 51(9): 2614-2627 DOI URL
[23]	Lachaud J, Vignoles GL.A Brownian motion technique to simulate gasification and its application to C/C composite ablation. Computational Materials Science, 2009, 44(4): 1034-1041 DOI URL
[24]	黄飞, 俞继军, 李秀涛等. 碳/碳材料细观尺度的氧扩散特性与烧蚀分析. 宇航学报, 2011, 32(8): 1848-1853 DOI URL
	(Huang Fei, Yu Jijun, Li Xiutao, et al.Oxygen diffusion characteristics and ablation analysis at mesoscale for carbon/carbon materials. Acta Astronautica Sinica, 2011, 32(8): 1848-1853 (in Chinese)) DOI URL
[25]	Deng D, Yu J, Yan X, et al.Engineering method for the thermal mechanical erosion of C/C composite with the mesoscale ablation model. Polymers & Polymer Composites, 2014, 22(2): 181-6
[26]	卞荫贵, 徐立功. 气动热力学. 北京: 中国科学技术大学出版社, 2011
	(Bian Yingui, Xu Li Gong.Aerodynamics. Beijing: University of Science & Technology China Press, 2011 (in Chinese))
[27]	易法军, 梁军, 孟松鹤等. 防热复合材料的烧蚀机理与模型研究. 固体火箭技术, 2000, 23(4): 48-56 DOI URL
	(Yi Fajun, Liang Jun, Meng Songhe, et al.Study on ablation mechanism and model of thermal-resistant composites. Solid Rocket Technology, 2000, 23(4): 48-56 (in Chinese)) DOI URL
[28]	Jamaluddin AS, Truelove JS, Wall TF.Intrinsic reactivity of carbons to oxygen. Fuel, 1983, 62(4): 484-486 DOI URL
[29]	Teigen KE, Li X, Lowengrub J, et al.A diffuse-interface approach for modeling transport, diffusion and adsorption/desorption of material quantities on a deformable interface. Communications in Mathematical Sciences, 2009, 4(7): 1009 DOI URL PMID
[30]	Gueyffier D, Li J, Nadim A, et al.Volume-of-fluid interface tracking with smoothed surface stress methods for three-dimensional flows. Journal of Computational Physics, 1999, 152(2): 423-456 DOI URL
[31]	杨云军, 崔尔杰, 周伟江. 细长三角翼摇滚运动数值研究. 空气动力学学报, 2007, 25(1): 34-44 DOI URL
	(Yang Yunjun, Cui Erjie, Zhou Weijiang.Numerical study of rock and roll movement of slender delta wing. Journal of Aerodynamics, 2007, 25(1): 34-44 (in Chinese)) DOI URL
[32]	Farhat C, Zee KGVD, Geuzaine P.Provably second-order time-accurate loosely-coupled solution algorithms for transient nonlinear computational aeroelasticity. Computer Methods in Applied Mechanics & Engineering, 2006, 195(17): 1973-2001 DOI URL
[33]	许阳阳. 低密度碳/酚醛复合材料烧蚀条件下的热$\!$-$\!$-$\!$化学耦合分析. [硕士论文]. 哈尔滨:哈尔滨工业大学, 2016
	(Xu Yangyang.Thermal-chemical coupling analysis of low density carbon/phenolic composites under ablation condition. [Master Thesis]. Harbin:Harbin Institute of Technology, 2016 (in Chinese))

NUMERICAL SIMULATION OF MICRO-ABLATION BEHAVIOR FOR CARBON FIBER REINFORCED COMPOSITES¹⁾

RichHTML

PDF (PC)

Knowledge

Abstract

Cite this article

share this article

Figures/Tables 19

References 33

Related Articles 15

Recommended Articles

Metrics

Comments

[1]	Wang Yang, Dong Gang. INTERFACE EVOLUTIONS AND GROWTH PREDICTIONS OF MIXING ZONE ON PREMIXED FLAME INTERFACE DURING RM INSTABILITY ¹⁾ [J]. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(6): 1655-1665.
[2]	Chen Yunmin,Ma Pengcheng,Tang Yao. CONSTITUTIVE MODELS AND HYPERGRAVITY PHYSICAL SIMULATION OF SOILS$^{\bf 1)}$ [J]. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(4): 901-915.
[3]	Ren Huilan, Chu Zhuxin, Li Jianqiao, Ma Tianbao. RESEARCH ON ELECTROMAGNETIC RADIATION DURING THE EXPLOSION PROGRESS OF COMPOSITION B EXPLOSIVES$^{\bf 1)}$ [J]. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(4): 1199-1210.
[4]	Wang Wenbo, Huang Ning, Dun Hongchao. ANALYSIS OF WIND-SAND MOVEMENT OVER SAND DUNE WITH DIFFERENT RAILWAY FORMS DOWNSTREAM¹⁾ [J]. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(3): 680-688.
[5]	Deng Bin,Wang Mengfei,Huang Zongwei,Wu Zhiyuan,Jiang Changbo. NUMERICAL SIMULATION OF THE HYDRODYNAMIC CHARACTERISTICS OF VIOLENT AERATED FLOWS NEAR VERTICAL STRUCTURE UNDER WAVE ACTION ¹⁾ [J]. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(2): 408-419.
[6]	Wang Yuerou, Wang Junfeng, Liu Hailong. NUMERICAL SIMULATION ON BUBBLE RINSING BEHAVIORS UNDER ELECTRIC FIELD ¹⁾ [J]. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(1): 31-39.
[7]	Zhan Jiemin, Li Yihua. A HYBRID TURBULENCE MODEL FOR WAVE BREAKING SIMULATION¹⁾ [J]. Chinese Journal of Theoretical and Applied Mechanics, 2019, 51(6): 1712-1719.
[8]	Du Xiaoqing, Qiu Tao, Zhao Yan. NUMERICAL INVESTIGATION OF MASS RATIO EFFECT ON FLOW-INDUCED VIBRATION OF TWO TANDEM SQUARE CYLINDERS AT LOW REYNOLDS NUMBER¹⁾ [J]. Chinese Journal of Theoretical and Applied Mechanics, 2019, 51(6): 1740-1751.
[9]	Xu Jianyu, Huang Shenghong. NUMERICAL SIMULATION OF CYLINDRICAL CONVERGING SHOCK INDUCED RICHTMYER-MESHKOV INSTABILITY WITH SPH¹⁾ [J]. Chinese Journal of Theoretical and Applied Mechanics, 2019, 51(4): 998-1011.
[10]	Hong Zheng, Ye Zhengyin. STUDY ON EVOLUTION CHARACTERISTICS OF ISOTROPIC TURBULENCE PASSING THROUGH A NORMAL SHOCK WAVE¹⁾ [J]. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(6): 1356-1367.
[11]	Wu Jian, Zhang Mengqi, Tian Fang-Bao. NUMERICAL ANALYSIS OF THREE-DIMENSIONAL ELECTRO-CONVECTION OF DIELECTRIC LIQUIDS IN A CUBICAL CAVITY¹⁾ [J]. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(6): 1458-1469.
[12]	Li Jianqiao, Ma Tianbao, Ning Jianguo. MECHANISM OF EXPLOSION-INDUCED DISTURBANCE IN NATURAL MAGNETIC FIELD¹⁾ [J]. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(5): 1206-1218.
[13]	Zhou Jianhong, Tong Baohong, Wang Wei, Su Jialei. DEFORMATION AND RUPTURE OF BUBBLE WHEN THE HOLLOW DROPLET IMPACTS ON THE OIL FILM [J]. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(2): 427-437.
[14]	Tong Fulin, Li Xin, Yu Changping, Li Xinliang. DIRECT NUMERICAL SIMULATION OF HYPERSONIC SHOCK WAVE AND TURBULENT BOUNDARY LAYER INTERACTIONS ¹⁾ [J]. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(2): 197-208.
[15]	Liu Zhaomiao, Wang Wenkai, Pang Yan. INFLUENCE OF EXPANDED CAVITY ON MIXING PERFORMANCE OF SQUARE-WAVE MICRO-MIXER [J]. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(2): 254-262.