Impact of vapor pressure on ablation rate of laser-irradiated ablate target
-
摘要: 研究了强激光辐照碳/碳复合材料靶材引起的烧蚀现象及蒸气压对烧蚀速率的影响。基于傅里叶定律,建立了强激光辐照靶材的热传导模型,模拟了忽略蒸气压影响时烧蚀速率随功率的变化;通过Mott-smith近似方法描述了Knudsen层间断区域,分析了间断两侧表面粒子状态参数;结合质量连续方程和蒸气压与温度关系方程,并由气体状态方程描述蒸气流状态,对蒸气压条件下激光烧蚀碳/碳复合材料靶材的速率随功率变化的关系进行了数值模拟。结果表明,在高能激光对靶材的烧蚀过程中,蒸气压力变化会导致靶材的饱和蒸气温度发生变化,进而影响烧蚀速率且使其随功率呈非线性变化,与忽略蒸气压作用时的线性变化规律相差较大,从理论上解释了忽略蒸气压导致的实验数据与理论结果的差异。Abstract: The ablation phenomena induced by high energy laser irradiation of carbon/carbon composites and the effect of vapor pressure on the ablation rate are investigated.First, the heat conduction model of the target irradiated by high power laser is set up based on Fourier law. The discontinuous region of Knudsen layer is described and the state parameters of the particles on both sides of discontinuities are analyzed by the Mott-smith approximation method. Then, combined with the mass continuity equation and the relation equation between vapor pressure and temperature equations, and the gas state equation to describe the vapor flow state, the numerical simulation about the relationship of the laser ablation rate of C/C composite material under the condition of vapor pressure with the power change is carried out. The results show that vapor pressure change will lead to the change of target vapor temperature in the process of high energy laser ablation of target, which will affect the ablation rate and make the rate change with power nonlinearly. It is very different from the linear variation when ignoring the vapor pressure. And this theoretically explains the difference between experimental data and theoretical results caused by ignoring the vapor pressure.
-
Key words:
- laser /
- vapor pressure /
- ablation rate /
- Knudsen layer
-
激光自被发明以来,在材料加工[1-3]、国防军事和其他科技领域都具有巨大潜力,研究激光烧蚀作用的影响因素具有很强的理论与实际意义。近年来对激光的研究手段日益成熟[4],其中关于激光辐照时间[5]、入射强度、靶材材料等[6]影响烧蚀速率因素的理论和实验研究已有成果,然而对气压条件的研究均局限于环境气体[7-8]以及冲击压强等因素[9],蒸气压在其中的影响大多忽略处理[10],这有碍于对烧蚀过程的进一步研究。关于蒸气压对激光烧蚀作用的影响,已有的文献很少涉及。本文采用傅里叶热传导模型计算激光对靶材的烧蚀速率,然后分析烧蚀靶材产生的蒸气粒子运动情况并列出质量连续方程,计算了Knudsen层两侧的间断状态参数和考虑蒸气压时激光对靶材的烧蚀速率。
1. 强激光辐照靶材的热传导模型
百兆瓦级功率激光辐照靶材时,沿着靶材厚度方向形成了梯度非常大的温度场,靶面温度在数秒内上升到气化温度并开始烧蚀。烧蚀掉的靶面物质高速喷溅,物质粒子相互碰撞,在靶面前方数个分子平均自由程(μm级)内达到平衡,并发生温度超过标准状况下气化温度的现象,这个区域被称作Knudsen层[11](靶材上方阴影部分)。粒子经过Knudsen层后形成稳定的高温蒸气流(如图 1)。
图 2上表面为靶材表面,发生烧蚀时以一定速度后退,沿激光辐照方向引起的温度场梯度比平行于靶材表面方向的温度场梯度要大得多,同时激光与靶材作用时间很短暂,故采用一维热传导控制方程[12],物性参数为常数,并且不考虑内热源。材料的热弛豫时间在ps到ns量级,而通常纳秒激光脉宽不小于10 ns, 一般不会发生非Fourier效应,因此计算采用一维的Fourier热传导控制方程[13]
ρc∂T∂t=K∂2T∂x2 (1) 式中:x是与热量传递方向相平行的坐标;K为热导率;ρ为介质密度;c为比热容;T为温度。假定材料出现升华烧蚀,且暂不考虑蒸气压对烧蚀的影响。考虑烧蚀时入射热流、向靶材内部传导的能量以及材料升华带走的能量三者能量守恒,得到边界条件:
(1)未发生烧蚀时:q=K∂T(x,t)∂x=Q−εσT4w,(T<Ts,x=h,t⩾0) (2) (2)发生烧蚀时:q=K∂T(x,t)∂x=Q−εσT4w−ρ˙slv,(T⩾Ts,x=L(t),t⩾0) (3) 靶材背面采用绝缘假设,边界条件为∂T(x,t)∂x=0(x=0,t⩾0),初始条件为T(x, t)=300 K(t=0, 0≤x≤L)。其中:Q=φI为激光功率等效的靶材入射热流,φ为吸收系数;ε为辐射系数,σ为Boltzman常数;Tw为受激光辐照时表面的瞬时温度;Ts为升华温度;h为初始板厚;L(t)为发生烧蚀后t时刻板厚; ˙s为线升华,即烧蚀速率;lv为受热靶材升华潜热。
对数理模型离散后运用基于Gauss消元法的Thomas算法(TDMA)可求解靶材在任意时刻任意位置的温度[14]。
2. 基于蒸气压的烧蚀过程分析
2.1 Knudsen层的间断条件
激光辐照使靶材表面气化形成稳定蒸气流,过渡区域为Knudsen层, 其内侧是靶材蒸发表面,外侧是稳定的蒸气流,两侧粒子温度、密度和压力形成了间断[15-16]。由Maxwell分布规律分别计算靶表面和Knudsen层外表面粒子速率:
f+s=ρ(2πRT)−3/2exp[−m(v2x+v2y+v2z)2σT] (4) f+k=ρ(2πRT)−3/2exp{−m[(vx−u)2+v2y+v2z]2σT} (5) 式中: vx,vy,vz为蒸气粒子速度分量;垂直离开靶表面的运动速度vx=√2σT/πm;m为原子质量;u为蒸气粒子平均速度;R是气体常数。
对两侧应用质量、动量和能量守恒条件,并将速率分布方程(4)~(5)代入,得到两侧温度和密度间断条件:
TkTs={[1+π(γ−1γ+1α2)2]12−π12γ−1γ+1α2}2 (6) ρkρs=√TsTk[(α2+12)eα2erfc(α)−α√π]+12TsTk[1−√παeα2erfc(α)] (7) 式中:下标k,s分别表示层外边界和靶材相变表面;α=u/√2kTk/m;比热比γ=Cp/Cv取5/3,Cp为定压比热容,Cv为定容比热容;erfc(α)为余误差函数。
2.2 连续方程
Knudsen层外表面的马赫数为
Mak=ukck=α√2γ (8) 式中:声速ck=√γkTk/m。
蒸气粒子流会推动外层空气形成冲击波,当激光功率足够强时,就必须考虑超声速情形。考虑各区域粒子运动情况并参考超声速波点燃阈值[17-18]后,取Mak=1。
由于蒸气粒子不全是沿垂直方向离开靶面方向,还有部分从界面上返回的,因此用β(β < 0)表示Knudsen层向内返回粒子的因子,可以列出质量连续方程:
ρcuc=(1−β)ρkuk (9) β=ρsρk√TsTkexp[(2α2+1)−α√πTsTk][α√πeα2erfc(α)−1] (10) 式中:ρcuc为靶材质量损失速率。
Knudsen层界面上的物理量可以用Clasius-Clayperon方程描述:
pkp0=exp[MlvRTv(1−TvTk)] (11) 式中:p0,Tv分别为ρk/ρs标准状况下饱和蒸气压力和温度。由该方程可知压力会影响气化温度进而影响质量损失率,因此计算时要耦合进去。
蒸气流的状态可以用气体状态方程描述:
p=nkT (12) 式中:n为粒子的密度。
结合式(9)~(12),并应用间断条件可以得到Knudsen层界面处的温度和压力,将新的温度和压力值带回到热传导方程(1)和边界条件(2)和(3),可以得到受蒸气压影响的烧蚀速率。
3. 计算分析
由于强激光辐照下,表面温度在很短时间内就会达到升华温度,所以参数取为常数。采用C/C复合材料作为靶材,工况参数如下:密度ρ=1800 kg/m3,比热容c=1870 J/(kg·K),热导率K=70 W/(m·K),潜热lv=5.903 4×107 J/kg,升华温度Tv=3873 K,表面热辐射系数ε=0.78,Boltzman常数σ=5.67×10-8 W/(m2·K4),吸收系数φ=0.7,初始温度T0=300 K。
对上述热传导模型进行模拟,在实验得到烧蚀初始时间的基础上,从105~108 W/cm2量级间选取不同功率密度值,研究了忽略蒸气压时烧蚀速率随激光功率密度的变化规律,数据点拟合结果见图 3(纵坐标为对数)。可以看出随着功率密度的提升,烧蚀速率的增长几乎呈线性。
取直径为1 cm的靶材,得到表面功率不同时靶材烧蚀规律,图 4和图 5分别为Knudsen层蒸气压和表面温度随激光功率密度变化的结果,均显示为随激光功率密度的变大而增长。这是因为功率变大时表面激光转化成的等效热流值变大,烧蚀速率随之增加,溅射出来的粒子数变大,导致蒸气的饱和压力变大,由Clasius-Clayperon方程可知温度也会随之增加。
图 6显示了表面蒸气压力对烧蚀速率的影响,实线是未考虑蒸气压的计算结果,而虚线考虑了蒸气压的影响。蒸气压对烧蚀速率的影响比较大,并使其随功率变化表现出非线性规律,这是因为功率密度较高时饱和蒸气压非常大,影响了靶材的升华烧蚀,而烧蚀速率的变化反过来又会影响饱和蒸气压的值[15, 19-20]。
4. 结论
(1) 分析了离开靶材表面粒子形成的蒸气流,包括Knudsen层的形成及其两侧参数变化,计算发现靶材表面饱和蒸气温度高于标准气压下材料气化温度,随着激光功率的增大,靶材表面蒸气压力变大,表面温度也随之增加。
(2) 随着表面温度的增加,进一步分析对烧蚀过程的影响,与不考虑蒸气压条件的情况对比发现:蒸气压对烧蚀速率的影响比较大,使其随功率变化表现出非线性规律。
-
[1] Guizard S, Semerok A, Gaudin J, et al. Femtosecond laser ablation of transparent dielectrics: measurement and modelisation of crater profiles[J]. Appl Surf Sci, 2002, 186: 364-368. doi: 10.1016/S0169-4332(01)00681-X [2] 陈国夫. 飞秒激光与透明介质的相互作用[J]. 物理, 2005, 34(10): 725-730. doi: 10.3321/j.issn:0379-4148.2005.10.005Chen Guofu. Interaction of femto-second laser pulses with transparent materials. Physics, 2005, 34(10): 725-730 doi: 10.3321/j.issn:0379-4148.2005.10.005 [3] Sahin R, Simsek E, Akturk S. Nanoscale patterning of graphene through femtosecond laser ablation[J]. Appl Phys Lett, 2014, 104: 053118. [4] 冯培培, 吴寒, 张楠. 超短脉冲激光烧蚀石墨产生的喷射物的时间分辨发射光谱研究[J]. 物理学报, 2015, 64: 214201. doi: 10.7498/aps.64.214201Feng Peipei, Wu Han, Zhang Nan. Study of the time-resolved emission spectra of the ejected plume generated by ultrashort laser ablation of graphite. Acta Physica Sinica, 2015, 64: 214201 doi: 10.7498/aps.64.214201 [5] 郭亚林, 梁国正, 丘哲明, 等. 激光参数对碳纤维复合材料质量烧蚀率的影响[J]. 复合材料学报, 2006, 23(5): 84-88. https://www.cnki.com.cn/Article/CJFDTOTAL-FUHE200605015.htmGuo Yalin, Liang Guozheng, Qiu Zheming, et al. Effect of laser parameters on mass ablative rate of carbon fiber reinforced composite. Acta Mater Compos, 2006, 23(5): 84-88 https://www.cnki.com.cn/Article/CJFDTOTAL-FUHE200605015.htm [6] 姚黎为, 王新兵, 刘璐宁, 等. 脉冲CO2激光烧蚀锡靶等离子体的数值模拟[J]. 强激光与粒子束, 2016, 28: 112007. doi: 10.11884/HPLPB201628.160170Yao Liwei, Wang Xinbing, Liu Luning, et al. Numerical simulation of pulsed CO2 laser produced tin plasma. High Power Laser and Particle Beams, 2016, 28: 112007 doi: 10.11884/HPLPB201628.160170 [7] 傅广生, 丁学成, 郭瑞强, 等. 脉冲激光沉积纳米硅晶粒流体模型的推广[J]. 物理学报, 2011, 60: 018102. doi: 10.7498/aps.60.018102Fu Guangsheng, Ding Xuecheng, Guo Ruiqiang, et al. The extended inertia fluid model to interpret the size distribution of Si nanoparticles prepared by pulsed laser ablation. Acta Physica Sinica, 2011, 60: 018102 doi: 10.7498/aps.60.018102 [8] 丁学成, 傅广生, 褚立志, 等. 环境气体种类对激光烧蚀粒子速度劈裂的影响[J]. 物理学报, 2012, 61: 155207. doi: 10.7498/aps.61.155207Ding Xuecheng, Fu Guangsheng, Chu Lizhi, et al. Influence of gas type on velocity splitting of ablated particles. Acta Physica Sinica, 2012, 61: 155207 doi: 10.7498/aps.61.155207 [9] 王文亭, 张楠, 王明伟, 等. 飞秒激光烧蚀金属靶的冲击温度[J]. 物理学报, 2013, 62: 170601. https://www.cnki.com.cn/Article/CJFDTOTAL-WLXB201321011.htmWang Wenting, Zhang Nan, Wang Mingwei, et al. Shock temperature of femtosecond laser ablation of solid target. Acta Physica Sinica, 2013, 62: 170601 https://www.cnki.com.cn/Article/CJFDTOTAL-WLXB201321011.htm [10] 梁亦寒, 胡广月, 袁鹏, 等. 纳秒激光烧蚀固体靶材产生的等离子体在外加横向磁场中膨胀时的温度和密度参数演化[J]. 物理学报, 2015, 64: 125204. doi: 10.7498/aps.64.125204Liang Yihan, Hu Guangyue, Yuan Peng, et al. Temporal evolutions of the plasma density and temperature of laser-produced plasma expansion in an external transverse magnetic field. Acta Physica Sinica, 2015, 64: 125204 doi: 10.7498/aps.64.125204 [11] 陆建, 倪晓武, 贺安之. 激光与材料相互作用物理学[M]. 北京: 机械工业出版社, 1996.Lu Jian, Ni Xiaowu, He Anzhi. Physics of interaction between laser and materials. Beijing: China Machine Press, 1996 [12] 黄海明, 孙岳. 脉冲强激光辐照下材料响应的非傅里叶效应[J]. 强激光与粒子束, 2009, 21(6): 808-812. http://www.hplpb.com.cn/article/id/4044Huang Haiming, Sun Yue. Non-Fourier response of target irradiated by multi-pulse high power laser. High Power Laser and Particle Beams, 2009, 21(6): 808-812 http://www.hplpb.com.cn/article/id/4044 [13] 胡汉平. 热传导理论[M]. 1版. 合肥: 中国科学技术大学出版社, 2010.Hu Hanping. Thermal conduction theory. 1st ed. Hefei: University of Science and Technology of China Press, 2010 [14] 赵卫, 乔玲, 韩晓林, 等. C/C-SiC复合材料的表面烧蚀模型及数值模拟[J]. 东南大学学报(自然科学版), 2011, 41(2): 365-369.Zhao Wei, Qiao Ling, Han Xiaolin, et al. Surface ablation model and numerical simulation of C/C-SiC composites. J Southeast Univ (Nat Sci Ed), 2011, 41(2): 365-369 [15] Chan C L, Mazumder J. One-dimensional steady-state mode for damage by vaporization and liquid expulsion due to laser-material interaction[J]. J App Phys, 1987, 62(11): 4579-4586. [16] Leonid V Z, Barbara J G. Velocity distributions of molecules ejected in laser ablation[J]. App Phys Lett, 1997, 71(4): 551-553. [17] 袁钢, 周光泉. 用于等离子体及LSD波点燃阈值的判据[J]. 高压物理学报, 1988, 2(2): 182-191. https://www.cnki.com.cn/Article/CJFDTOTAL-GYWL198802014.htmYuan Gang, Zhou Guangquan. Criteria of thresholds for plasma and LSD wave ignition. Chin J High Pressure Phys, 1988, 2(2): 182-191 https://www.cnki.com.cn/Article/CJFDTOTAL-GYWL198802014.htm [18] 杨波, 朱金荣, 杨雁南, 等. 由激光靶冲量耦合实验结果判定激光支持爆轰波点燃阈值[J]. 中国激光, 2007, 34(1): 139-144. https://www.cnki.com.cn/Article/CJFDTOTAL-JJZZ200701027.htmYang Bo, Zhu Jinrong, Yang Yannan, et al. Determination of the laser supported detonation wave ignition threshold from impulse coupling between laser and target. Chin J Lasers, 2007, 34(1): 139-144 https://www.cnki.com.cn/Article/CJFDTOTAL-JJZZ200701027.htm [19] Batanov V A, Bunkin F V, Prokhorov A M, et al. Evaporation of metallic targets caused by intense optical radiation[J]. Sov Phys JETP, 1973, 36(2): 311-322. [20] Kelly R. Gas dynamics of the pulsed emission of a perfect gas with applications to laser sputtering and to nozzle expansion[J]. Phys Rev A, 1992, 46(2): 860-874. 期刊类型引用(1)
1. 耿莉,成溯,付前刚,李贺军. 碳/碳复合材料的激光烧蚀行为与机制. 复合材料学报. 2022(09): 4337-4343 . 百度学术
其他类型引用(0)
-