激光烧蚀掺杂金属聚合物推力产生过程数值模拟

推进技术 ›› 2017, Vol. 38 ›› Issue (2): 257-265.

激光烧蚀掺杂金属聚合物推力产生过程数值模拟

段兴跃¹,李小康¹,程谋森¹,李干²

国防科学技术大学航天科学与工程学院，湖南长沙 410000,国防科学技术大学航天科学与工程学院，湖南长沙 410000,国防科学技术大学航天科学与工程学院，湖南长沙 410000,解放军理工大学爆炸冲击防灾减灾国家重点实验室，江苏南京 210007

发布日期:2021-08-15
作者简介:段兴跃，男，硕士生，研究领域为烧蚀模式激光推进技术。
基金资助:
国家自然科学基金青年科学基金(51306203)；国防科学技术大学科研计划项目(JC14-01-02)。

Numerical Simulation of Thrust Generation Process of Laser Ablating Polymer Doped Metal

Department of Aerospace Science and Engineering，National University of Science and Technology，Changsha 410000，China,Department of Aerospace Science and Engineering，National University of Science and Technology，Changsha 410000，China,Department of Aerospace Science and Engineering，National University of Science and Technology，Changsha 410000，China and National Key Laboratory of Explosion Shock and Disaster Prevention and Reduction， PLA University of Science and Technology，Nanjing 210007，China

Published:2021-08-15

摘要/Abstract

摘要： 为了增进对激光烧蚀推进中推力产生过程的认识，建立了激光烧蚀掺杂聚合物推力产生过程的模型，发展了一套包括激光能量沉积、工质烧蚀、烧蚀羽流飞散等过程的数值仿真程序。数值计算了真空中3～40J/cm2激光烧蚀掺杂微米铝颗粒聚甲醛工质的推力、烧蚀轮廓及质量、压强分布和比冲特性，且比冲与实验数据较吻合。计算结果表明:高激光能量密度(30.0J/cm2)较低激光能量密度(5.0J/cm2)的金属颗粒剥蚀情况严重；低激光能量密度(5.0J/cm2)下推力时间变化规律较简单，总体呈现先增后减的趋势，且从整体上看只有一个压强峰值；而相对较高的激光能量密度(30.0J/cm2)下，由于存在“烧蚀-屏蔽-烧蚀被削弱”的制约关系，推力时间变化规律复杂；流场产生的高压区较多，且呈现交替产生、并存发展到衰减消失的规律；羽流与激光的相互作用更为剧烈，峰值压强也更大。

关键词: 激光烧蚀推进；掺杂聚合物；推力产生机理；数值模拟

Abstract: In order to improve the understanding of thrust generation process in laser ablation propulsion，a thrust generation model of laser ablating doped polymer was established，and a numerical simulation program considering the process of laser energy deposition，propellant ablation and expansion of ablation products was developed. Characteristics of thrust，ablation profile，ablation mass，pressure distribution and specific impulse of laser ablating polyoxymethylene (POM) doped micro-scale aluminum particles between the laser fluence of 3～40J/cm2 in vacuum were calculated numerically，and the obtained specific impulse shows a good agreement with the corresponding experimental values. The numerical results show that the erosion of metal particles under a higher laser fluence (30.0J/cm2) is severer than that under a lower laser fluence (5.0J/cm2). The curve of the temporal variation of thrust is simple，demonstrating a trend of first rise and then declining with only one pressure peak as a whole under the 5.0J/cm2 case. By contrast，the curve of the temporal variation of thrust under a higher laser fluence (30.0J/cm2) is more complex because of the relationship of‘ablation-shielding-weakened ablation’during the laser ablation process. Moreover，the number of high pressure zones is much more than that under a lower laser fluence of 5.0J/cm2. Under the 30.0J/cm2 case，it shows a trend of alternate generation，coexistence，development，and eventual dissipation. Thus，the interaction between the ablation plume and the laser is stronger，and the peak pressure is much bigger as well.

Key words: Laser ablation propulsion；Doped polymer；Thrust generation mechanism；Numerical simulation

段兴跃,李小康,程谋森,李干. 激光烧蚀掺杂金属聚合物推力产生过程数值模拟[J]. 推进技术, 2017, 38(2): 257-265.

DUAN Xing-yue¹,LI Xiao-kang¹,CHENG Mou-seng¹,LI Gan². Numerical Simulation of Thrust Generation Process of Laser Ablating Polymer Doped Metal[J]. Journal of Propulsion Technology, 2017, 38(2): 257-265.

参考文献

[1] Kantrowitz A. Propulsion to Orbit by Ground Based Lasers[J]. Astronautics and Aeronautics,1972, 10(5):74-76.
[2] Phipps C, Birkan M, Bohn W, et al. Review: Laser-Ablation Propulsion[J]. Journal of Propulsion and Power, 2010, 26(4): 609-637.
[3] Mead Jr F B, Myrabo L N, Messitt D G . Flight Experiments and Evolutionary Development of a Laser-Propelled Transatmospheric Vehicle[C]. Santa Fe, NM, United States: High-Power Laser Ablation, 1998.
[4] Watanabe K, Sasoh A. Laser Impulse Generation Required for Space Debris Deorbiting[C]. Taos, NM, United States: High-Power Laser Ablation, 2004.
[5] Schall W O, Tegel J, Eckel H A. Ablation Performance Experiments with Metal Seeded Polymers[C]. Troy, NY, United States: 3rd International Symposium on Beamed Energy Propulsion, 2006.
[6] 常浩, 金星, 叶继飞, 等. 激光功率密度对纳秒激光烧蚀冲量耦合系数影响的数值模拟[J]. 推进技术, 2013, 34(10): 1426-1431. (CHANG Hao, JIN Xing, YE Ji-fei, et al. Numerical Simulation of Laser Power Density Effect on Nanosecond Laser Ablation Impulse Coupling [J]. Journal of Propulsion Technology, 2013, 34(10): 1426-1431.)
[7] 彭杰, 郑航, 胡晓军, 等. 激光烧蚀掺杂金属粉工质的推进性能[J]. 强激光与粒子束, 2009, 21(6): 821-825.
[8] 胡春波, 何国强, 程翔, 等. 不同辅助启动材料对激光推进器冲量的影响 [J]. 推进技术, 2007, 28(5): 475-477. (HU Chun-bo, HE Guo-qiang, CHENG Xiang, et al. Influence of Different Aid Initiation Materials on Laser Thruster Impulse[J]. Journal of Propulsion Technology, 2007, 28(5): 475-477.)
[9] Kundrapu M, Keidar M. Laser Ablation: Beam Wavelength Optimization and Communication Through Plasma[C]. San Antonio, Texas, United States: 40th AIAA Plasmadynamics and Lasers Conference, 2009.
[10] 童慧峰. 烧蚀模式激光推进的机理及实验研究[D]. 合肥:中国科学技术大学, 2006.
[11] Sakai T, Ichihashi K, Matsuda A, et al. Calculation of Pulsed Laser-Ablative Impulse on Polyacetal[C]. San Antonio, Texas, United States: 40th AIAA Plasmadynamics and Lasers Conference, 2009.
[12] Sinko J E, Phipps C R. Modeling CO₂ Laser Ablation Impulse of Polymers in Vapor and Plasma Regimes[J]. Applied Physics Letters, 2009, 95(13).
[13] Li Nanlei, Yanji Hong, Xiuqian Li. Numerical Study of Thrust Generation in the Process of Laser Ablated Doped Polymer[C]. Ludwigsburg, Germany: Seventh International Symposium on Beamed Energy Propulsion, 2011.
[14] B?uerle D. Laser Processing and Chemistry, Fourth Edition[M]. Berlin: Springer-Verlag, 2011.
[15] 李干. 脉冲激光辐照聚合物工质烧蚀推进机理与性能优化研究[D]. 长沙:国防科学技术大学, 2014.
[16] Bityurin N, Luk Yanchuk B S , Hong M H , et al. Models for Laser Ablation of Polymers[J]. Chemical Reviews, 2003, 103(2): 519-552.
[17] Sakai T, Sasoh A, Anju K, et al. Computational Simulation of Local Impulse Generation in Intense Laser-Solid Interaction[C]. Reno, Nevada, United States: 46th AIAA Aerospace Sciences Meeting and Exhibit, 2008.
[18] Pekker L, Keidar M, Cambier J L. Effect of Thermal Conductivity on the Knudsen Layer at Ablative Surfaces[J]. Journal of Applied Physics, 2008, 103(3).
[19] 孙承纬. 激光辐照效应[M]. 北京:国防工业出版社, 2002.
[20] Lippert T. Laser Application of Polymers[M]. Berlin: Polymers and Light, Springer, 2004.
[21] 常浩, 叶继飞, 周伟静. 纳秒激光烧蚀铝靶羽流流场能量转化过程分析[J]. 推进技术, 2015, 36(10): 1588-1594. (CHANG Hao, YE Ji-fei, ZHOU Wei-jing. Energy Balance of Plume Flow Field for Nanosecond Laser Ablation Al[J]. Journal of Propulsion Technology, 2015, 36(10): 1588-1594.)
[22] André P. Comparison and Thermodynamic Properties of Ablated Vapors of PMMA, PA6-6, PETP, POM and PE [J]. Journal of Physics D: Applied Physics, 1996, 29(7): 1963-1972.
[23] 李干, 程谋森, 李小康. 基于压力间断解析解对AUSM+-up格式的修正[J]. 国防科技大学学报, 2014, 36(3): 1-5.
[24] E K Donald. 计算机程序设计艺术. 第二卷, 半数值算法(第三版)[M]. 苏运霖, 译. 北京:国防工业出版社, 2003.
[25] 彭杰. 烧蚀模式激光推力器的实验研究[D]. 合肥:中国科学技术大学, 2009.
[26] Schall W O, Eckel H A, Tegel J, et al. Properties of Laser Ablation Products of Delrin with CO₂ Laser[R]. FA8655-03-1-3061.　　　　　　　　　　　　　　收稿日期:2016-09-28；修订日期:2016-10-21。基金项目:国家自然科学基金青年科学基金(51306203)；国防科学技术大学科研计划项目(JC14-01-02)。作者简介:段兴跃,男,硕士生,研究领域为烧蚀模式激光推进技术。E-mail: duanxingyue@163.com(编辑:张荣莉)