留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

蒸气压对激光辐照靶材烧蚀速率的影响

姜学东 陈纪然 王彧 王超

姜学东, 陈纪然, 王彧, 等. 蒸气压对激光辐照靶材烧蚀速率的影响[J]. 强激光与粒子束, 2018, 30: 021002. doi: 10.11884/HPLPB201830.170271
引用本文: 姜学东, 陈纪然, 王彧, 等. 蒸气压对激光辐照靶材烧蚀速率的影响[J]. 强激光与粒子束, 2018, 30: 021002. doi: 10.11884/HPLPB201830.170271
Jiang Xuedong, Chen Jiran, Wang Yu, et al. Impact of vapor pressure on ablation rate of laser-irradiated ablate target[J]. High Power Laser and Particle Beams, 2018, 30: 021002. doi: 10.11884/HPLPB201830.170271
Citation: Jiang Xuedong, Chen Jiran, Wang Yu, et al. Impact of vapor pressure on ablation rate of laser-irradiated ablate target[J]. High Power Laser and Particle Beams, 2018, 30: 021002. doi: 10.11884/HPLPB201830.170271

蒸气压对激光辐照靶材烧蚀速率的影响

doi: 10.11884/HPLPB201830.170271
基金项目: 

国家自然科学基金项目 11472037

国家自然科学基金项目 1272042

中央高校基本科研业务费专项资金 2015JBM085

详细信息
    作者简介:

    姜学东(1970-), 男,硕士,从事开关电源以及等离子体方面研究;xdjiang@bjtu.edu.cn

    通讯作者:

    陈纪然(1993-), 女,硕士研究生,从事开关电源以及等离子体方面研究;15121395@bjtu.edu.cn

  • 中图分类号: O536;O552.3

Impact of vapor pressure on ablation rate of laser-irradiated ablate target

  • 摘要: 研究了强激光辐照碳/碳复合材料靶材引起的烧蚀现象及蒸气压对烧蚀速率的影响。基于傅里叶定律,建立了强激光辐照靶材的热传导模型,模拟了忽略蒸气压影响时烧蚀速率随功率的变化;通过Mott-smith近似方法描述了Knudsen层间断区域,分析了间断两侧表面粒子状态参数;结合质量连续方程和蒸气压与温度关系方程,并由气体状态方程描述蒸气流状态,对蒸气压条件下激光烧蚀碳/碳复合材料靶材的速率随功率变化的关系进行了数值模拟。结果表明,在高能激光对靶材的烧蚀过程中,蒸气压力变化会导致靶材的饱和蒸气温度发生变化,进而影响烧蚀速率且使其随功率呈非线性变化,与忽略蒸气压作用时的线性变化规律相差较大,从理论上解释了忽略蒸气压导致的实验数据与理论结果的差异。
  • 激光自被发明以来,在材料加工[1-3]、国防军事和其他科技领域都具有巨大潜力,研究激光烧蚀作用的影响因素具有很强的理论与实际意义。近年来对激光的研究手段日益成熟[4],其中关于激光辐照时间[5]、入射强度、靶材材料等[6]影响烧蚀速率因素的理论和实验研究已有成果,然而对气压条件的研究均局限于环境气体[7-8]以及冲击压强等因素[9],蒸气压在其中的影响大多忽略处理[10],这有碍于对烧蚀过程的进一步研究。关于蒸气压对激光烧蚀作用的影响,已有的文献很少涉及。本文采用傅里叶热传导模型计算激光对靶材的烧蚀速率,然后分析烧蚀靶材产生的蒸气粒子运动情况并列出质量连续方程,计算了Knudsen层两侧的间断状态参数和考虑蒸气压时激光对靶材的烧蚀速率。

    百兆瓦级功率激光辐照靶材时,沿着靶材厚度方向形成了梯度非常大的温度场,靶面温度在数秒内上升到气化温度并开始烧蚀。烧蚀掉的靶面物质高速喷溅,物质粒子相互碰撞,在靶面前方数个分子平均自由程(μm级)内达到平衡,并发生温度超过标准状况下气化温度的现象,这个区域被称作Knudsen层[11](靶材上方阴影部分)。粒子经过Knudsen层后形成稳定的高温蒸气流(如图 1)。

    图  1  激光致等离子体示意图
    Figure  1.  Schematic diagram of laser induced plasma

    图 2上表面为靶材表面,发生烧蚀时以一定速度后退,沿激光辐照方向引起的温度场梯度比平行于靶材表面方向的温度场梯度要大得多,同时激光与靶材作用时间很短暂,故采用一维热传导控制方程[12],物性参数为常数,并且不考虑内热源。材料的热弛豫时间在ps到ns量级,而通常纳秒激光脉宽不小于10 ns, 一般不会发生非Fourier效应,因此计算采用一维的Fourier热传导控制方程[13]

    ρcTt=K2Tx2
    (1)
    图  2  热传导模型示意图
    Figure  2.  Model of thermal conduction

    式中:x是与热量传递方向相平行的坐标;K为热导率;ρ为介质密度;c为比热容;T为温度。假定材料出现升华烧蚀,且暂不考虑蒸气压对烧蚀的影响。考虑烧蚀时入射热流、向靶材内部传导的能量以及材料升华带走的能量三者能量守恒,得到边界条件:

    1q=KT(x,t)x=QεσT4w,(T<Ts,x=h,t0)
    (2)
    (2)q=KT(x,t)x=QεσT4wρ˙slv,(TTs,x=L(t),t0)
    (3)

    靶材背面采用绝缘假设,边界条件为T(x,t)x=0(x=0,t0),初始条件为T(x, t)=300 K(t=0, 0≤xL)。其中:Q=φI为激光功率等效的靶材入射热流,φ为吸收系数;ε为辐射系数,σ为Boltzman常数;Tw为受激光辐照时表面的瞬时温度;Ts为升华温度;h为初始板厚;L(t)为发生烧蚀后t时刻板厚; ˙s为线升华,即烧蚀速率;lv为受热靶材升华潜热。

    对数理模型离散后运用基于Gauss消元法的Thomas算法(TDMA)可求解靶材在任意时刻任意位置的温度[14]

    激光辐照使靶材表面气化形成稳定蒸气流,过渡区域为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[(vxu)2+v2y+v2z]2σT}
    (5)

    式中: vxvyvz为蒸气粒子速度分量;垂直离开靶表面的运动速度vx=2σT/πmm为原子质量;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(α)为余误差函数。

    Knudsen层外表面的马赫数为

    Mak=ukck=α2γ
    (8)

    式中:声速ck=γkTk/m

    蒸气粒子流会推动外层空气形成冲击波,当激光功率足够强时,就必须考虑超声速情形。考虑各区域粒子运动情况并参考超声速波点燃阈值[17-18]后,取Mak=1。

    由于蒸气粒子不全是沿垂直方向离开靶面方向,还有部分从界面上返回的,因此用β(β < 0)表示Knudsen层向内返回粒子的因子,可以列出质量连续方程:

    ρcuc=(1β)ρkuk
    (9)
    β=ρsρkTsTkexp[(2α2+1)απTsTk][απeα2erfc(α)1]
    (10)

    式中:ρcuc为靶材质量损失速率。

    Knudsen层界面上的物理量可以用Clasius-Clayperon方程描述:

    pkp0=exp[MlvRTv(1TvTk)]
    (11)

    式中:p0Tv分别为ρk/ρs标准状况下饱和蒸气压力和温度。由该方程可知压力会影响气化温度进而影响质量损失率,因此计算时要耦合进去。

    蒸气流的状态可以用气体状态方程描述:

    p=nkT
    (12)

    式中:n为粒子的密度。

    结合式(9)~(12),并应用间断条件可以得到Knudsen层界面处的温度和压力,将新的温度和压力值带回到热传导方程(1)和边界条件(2)和(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(纵坐标为对数)。可以看出随着功率密度的提升,烧蚀速率的增长几乎呈线性。

    图  3  烧蚀速率随激光功率密度变化规律
    Figure  3.  Ablation velocity varies with along laser power density

    取直径为1 cm的靶材,得到表面功率不同时靶材烧蚀规律,图 4图 5分别为Knudsen层蒸气压和表面温度随激光功率密度变化的结果,均显示为随激光功率密度的变大而增长。这是因为功率变大时表面激光转化成的等效热流值变大,烧蚀速率随之增加,溅射出来的粒子数变大,导致蒸气的饱和压力变大,由Clasius-Clayperon方程可知温度也会随之增加。

    图  4  Knudsen层压力随激光功率密度的变化
    Figure  4.  Pressure of Knudsen layer vs laser power density
    图  5  表面温度随激光功率密度的变化
    Figure  5.  Surface temperature vs laser power density

    图 6显示了表面蒸气压力对烧蚀速率的影响,实线是未考虑蒸气压的计算结果,而虚线考虑了蒸气压的影响。蒸气压对烧蚀速率的影响比较大,并使其随功率变化表现出非线性规律,这是因为功率密度较高时饱和蒸气压非常大,影响了靶材的升华烧蚀,而烧蚀速率的变化反过来又会影响饱和蒸气压的值[15, 19-20]

    图  6  蒸气压对烧蚀速率的影响
    Figure  6.  Effect of vapor pressure on ablation velocity

    (1) 分析了离开靶材表面粒子形成的蒸气流,包括Knudsen层的形成及其两侧参数变化,计算发现靶材表面饱和蒸气温度高于标准气压下材料气化温度,随着激光功率的增大,靶材表面蒸气压力变大,表面温度也随之增加。

    (2) 随着表面温度的增加,进一步分析对烧蚀过程的影响,与不考虑蒸气压条件的情况对比发现:蒸气压对烧蚀速率的影响比较大,使其随功率变化表现出非线性规律。

  • 图  1  激光致等离子体示意图

    Figure  1.  Schematic diagram of laser induced plasma

    图  2  热传导模型示意图

    Figure  2.  Model of thermal conduction

    图  3  烧蚀速率随激光功率密度变化规律

    Figure  3.  Ablation velocity varies with along laser power density

    图  4  Knudsen层压力随激光功率密度的变化

    Figure  4.  Pressure of Knudsen layer vs laser power density

    图  5  表面温度随激光功率密度的变化

    Figure  5.  Surface temperature vs laser power density

    图  6  蒸气压对烧蚀速率的影响

    Figure  6.  Effect of vapor pressure on ablation velocity

  • [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.005

    Chen 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.214201

    Feng 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.htm

    Guo 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.160170

    Yao 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.018102

    Fu 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.155207

    Ding 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.htm

    Wang 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.125204

    Liang 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/4044

    Huang 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.htm

    Yuan 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.htm

    Yang 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)

  • 加载中
图(6)
计量
  • 文章访问数:  1281
  • HTML全文浏览量:  328
  • PDF下载量:  192
  • 被引次数: 1
出版历程
  • 收稿日期:  2017-06-30
  • 修回日期:  2017-09-26
  • 刊出日期:  2018-02-15

目录

/

返回文章
返回