高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

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

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

组合脉冲激光铝等离子体的数值模拟研究

赵佳伟 吴涛 廖青 潘子萌 马子琦

引用本文:
Citation:

组合脉冲激光铝等离子体的数值模拟研究

    作者简介: 赵佳伟(1996-), 男, 硕士研究生, 现主要从事激光诱导等离子体的研究.
    通讯作者: 廖青, liaoqing@wit.edu.cn
  • 基金项目:

    广东省重大基础与应用基础研究资助项目 2019B030302003

  • 中图分类号: O539

Numerical simulation of combined pulsed lasers aluminum plasma

    Corresponding author: LIAO Qing, liaoqing@wit.edu.cn ;
  • CLC number: O539

  • 摘要: 为了获得组合脉冲激光诱导等离子体光谱性能增强的物理机制, 基于FLASH程序模拟计算了预脉冲参数对组合脉冲激光诱导的铝等离子体参数时空分布的影响, 获得了在不同预脉冲波长和不同预主脉冲延时下产生的铝等离子体电子温度、电子密度和烧蚀质量的空间演化规律。结果表明, 在预主脉冲总能量相同的情况下, 随着预脉冲波长从0.266μm变化到1.064μm, 高温环境氦气等离子体羽辉的空间范围从0.7cm增大到3.0cm, 但组合脉冲对靶材的烧蚀效率严重下降, 而铝等离子体的最大电子温度保持稳定; 此外, 组合脉冲的时间延迟低于100ns。该研究可为组合脉冲激光诱导击穿等离子体光谱增强技术提供理论参考。
  • Figure 1.  Comparison of FLASH and MED103 numerical simulation results

    Figure 2.  Comparison between simulation data and experimental results of FLASH

    Figure 3.  Theoretical model of laser induction

    Figure 4.  Variation of electron temperature on target surface with time under different pre-main pulse delay

    Figure 5.  Spatial and temporal distribution of electron temperature and electron density at 1064nm

    a—electron temperature b—electron density

    Figure 6.  Spatial and temporal distribution of electron temperature and electron density at 530nm

    a—electron temperature b—electron density

    Figure 7.  Spatial and temporal distribution of electron temperature and electron density at 266nm

    a—electron temperature b—electron density

    Figure 8.  Mass density distribution of target material at different pre-pulse wavelengths

    a—1064nm b—530nm c—266nm

    Figure 9.  Spatial distribution of laser energy deposition at different pre-pulse wavelengths

  • [1]

    HAO X J, SUN Y K. Simultaneous detection of multi-elements in coal based on laser-induced breakdown spectroscopy[J]. Laser Technology, 2020, 44(1): 119-124(in Chinese).
    [2]

    EASON R W. Pulsed laser deposition of thin films: Applications-led growth of functional materials[M]. New York, USA: John Wiley and Sons Inc, 2006: 368-382.
    [3]

    ZHANG Y, ZHANG D, WU J, et al. A thermal model for nanosecond pulsed laser ablation of aluminum[J]. AIP Advances, 2017, 7(7): 075010. doi: 10.1063/1.4995972
    [4]

    ZHANG P B, QIN Y, ZHAO J J, et al. Two-dimensional numerical simulation of nanosecond laser-ablation of aluminum material by nanosecond laser pulse[J]. Acta Physica Sinica, 2010, 59(10): 7120-7128(in Chinese). doi: 10.7498/aps.59.7120
    [5]

    KNIGHT A K, SCHERBARTH N L, CREMERS D A, et al. Characterization of laser-induced breakdown spectroscopy (LIBS) for application to space exploration[J]. Applied Spectroscopy, 2000, 54(3): 331-340. doi: 10.1366/0003702001949591
    [6]

    PHIPPS C R, LUKE J R, MCDUFF G G, et al. Laser-ablation-powered mini-thruster[J]. Proceedings of the SPIE, 2002, 4760: 833-842. doi: 10.1117/12.482038
    [7]

    ZHANG J, SUGIOKA K, WADA S, et al. Dual-beam ablation of fused quartz using 266nm and VUV lasers with different delay-times[J]. Applied Physics, 1997, A64(5): 477-481.
    [8]

    BARTHÉLEMY O, MARGOT J, CHAKER M, et al. Influence of the laser parameters on the space and time characteristics of an aluminum laser-induced plasma[J]. Spectrochimica Acta, 2005, B60(7/8): 905-914.
    [9]

    SHAIKH N M, HAFEEZ S, RASHID B, et al. Spectroscopic studies of laser induced aluminum plasma using fundamental, second and third harmonics of a Nd: YAG laser[J]. The European Physical Journal, 2007, D44(2): 371-379.
    [10]

    DROGOFF B L, MARGOT J, VIDAL F, et al. Influence of the laser pulse duration on laser-produced plasma properties[J]. Plasma Sources Science and Technology, 2004, 13(2): 223-230. doi: 10.1088/0963-0252/13/2/005
    [11]

    DIWAKAR P K, HARILAL S S, FREEMAN J R, et al. Role of laser pre-pulse wavelength and inter-pulse delay on signal enhancement in collinear double-pulse laser-induced breakdown spectroscopy[J]. Spectrochimica Acta, 2013, B87: 65-73.
    [12]

    JAFARABADI M A, MAHDIEH M H. Investigation of phase explosion in aluminum induced by nanosecond double pulse technique[J]. Applied Surface Science, 2015, 346: 263-269. doi: 10.1016/j.apsusc.2015.03.158
    [13]

    SCAFFIDI J, PENDER J, PEARMAN W, et al. Dual-pulse laser-induced breakdown spectroscopy with combinations of femtosecond and nanosecond laser pulses[J]. Applied Optics, 2003, 42(30): 6099-6106. doi: 10.1364/AO.42.006099
    [14]

    de GIACOMO A, DELL'AGLIO M, BRUNO D, et al. Experimental and theoretical comparison of single-pulse and double-pulse laser induced breakdown spectroscopy on metallic samples[J]. Spectrochimica Acta, 2008, B63(7): 805-816.
    [15]

    BENEDETTI P A, CRISTOFORETTI G, LEGNAIOLI S, et al. Effect of laser pulse energies in laser induced breakdown spectroscopy in double-pulse configuration[J]. Spectrochimica Acta, 2005, B60(11): 1392-1401.
    [16]

    AMOUYE F A, FÖRSTER D J, GHORBANFEKR H, et al. Atomistic simulation of ultra-short pulsed laser ablation of metals with single and double pulses: An investigation of the re-deposition pheno-menon[J]. Applied Surface Science, 2021, 537: 147775. doi: 10.1016/j.apsusc.2020.147775
    [17]

    POVARNITSYN M E, FOKIN V B, LEVASHOV P R, et al. Molecular dynamics simulation of subpicosecond double-pulse laser ablation of metals[J]. Physical Review, 2015, B92(17): 174104.
    [18]

    ORBAN C, FATENEJAD M, CHAWLA S, et al. A radiation-hydrodynamics code comparison for laser-produced plasmas: FLASH versus HYDRA and the results of validation experiments[J/OL]. (2013-06-07)[2021-09-23]. https://arxiv.org/abs/1306.1584.
    [19]

    WHITE J B S. Opening the extreme ultraviolet lithography source bottleneck developing a 13.5-nm laser-produced plasma source for the semiconductor industry[D]. Dublin, Ireland: University College Dublin, 2006: 158-169.
    [20]

    BOGAERTS A, CHEN Z, AUTRIQUE D. Double pulse laser ablation and laser induced breakdown spectroscopy: A modeling investigation[J]. Spectrochimica Acta, 2008, B63(7): 746-754.
    [21]

    MAO X, ZENG X, WEN S B, et al. Time-resolved plasma properties for double pulsed laser-induced breakdown spectroscopy of silicon[J]. Spectrochimica Acta, 2005, B60(7): 960-967.
    [22]

    YIN P Q, WANG X B, WU Y X, et al. Experimental study on water droplet plasma induced by pulse Nd: YAG laser[J]. Laser Technology, 2020, 44(6): 726-731(in Chinese).
  • [1] 占剑杨明江 . YAG激光诱导组合脉冲放电对材料表面强化的研究. 激光技术, 2010, 34(5): 577-580. doi: 10.3969/j.issn.1001-3806.2010.O5.001
    [2] 张萌方英武张广鹏 . 纳秒脉冲激光辐照铝靶碎片动态响应的仿真研究. 激光技术, 2023, 47(4): 541-546. doi: 10.7510/jgjs.issn.1001-3806.2023.04.015
    [3] 吴鹏范云茹郭嘉伟刘晓旭周强邓光伟宋海智王浟 . 纳秒脉冲激光加工高反射率铝膜. 激光技术, 2019, 43(6): 779-783. doi: 10.7510/jgjs.issn.1001-3806.2019.06.009
    [4] 袁红童慧峰李牧孙承纬 . 强激光加载真空中铝靶冲量耦合的数值模拟. 激光技术, 2012, 36(4): 520-523. doi: 10.3969/j.issn.1001-806.2012.04.021
    [5] 李洪敏左军卫徐健徐一汀彭钦军许祖彦 . 脉冲激光的非相干合成技术研究. 激光技术, 2015, 39(2): 237-241. doi: 10.7510/jgjs.issn.1001-3806.2015.02.020
    [6] 吕雪明李泽文张检民倪晓武 . 不同延时的组合脉冲激光致硅表面损伤研究. 激光技术, 2020, 44(6): 695-699. doi: 10.7510/jgjs.issn.1001-3806.2020.06.008
    [7] 张翔龚心弦丁亚茜王小音金操帆 . 基于嵌入式技术的可控光纤激光脉冲光源. 激光技术, 2016, 40(5): 711-715. doi: 10.7510/jgjs.issn.1001-3806.2016.05.019
    [8] 徐健梁勖刘冬生潘冰冰方晓罗乐 . 脉冲激光系统充电电源功率因数校正技术研究. 激光技术, 2020, 44(6): 732-737. doi: 10.7510/jgjs.issn.1001-3806.2020.06.015
    [9] 蒋晓东于纪言朱立坤 . 基于位置敏感探测器的组合导航技术研究. 激光技术, 2019, 43(3): 335-340. doi: 10.7510/jgjs.issn.1001-3806.2019.03.009
    [10] 李秋实罗洪胡永明倪明 . 三角波形激光脉冲在液体中致声的特性分析. 激光技术, 2010, 34(3): 409-412. doi: 10.3969/j.issn.1001-3806.2010.03.035
    [11] 杨洪亮金湘中修腾飞费鑫江叶颖 . 钢/铝异种金属光纤激光焊接数值模拟. 激光技术, 2016, 40(4): 606-609. doi: 10.7510/jgjs.issn.1001-3806.2016.04.031
    [12] 王超康轶凡田进寿刘虎林 . 两类单阿秒脉冲产生技术的相位依赖性分析. 激光技术, 2012, 36(4): 516-519. doi: 10.3969/j.issn.1001-806.2012.04.020
    [13] 吴涛王新兵唐建饶志明王少义 . 缓冲气压对CO2激光Al靶等离子体参量的影响. 激光技术, 2011, 35(6): 800-803. doi: 10.3969/j.issn.1001-3806.2011.06.021
    [14] 邱婉静姜三平郝晓剑 . 激光光幕靶坐标测量系统的设计和实现. 激光技术, 2012, 36(6): 735-737. doi: 10.3969/j.issn.1001-3806.2012.06.006
    [15] 田伟周惦武乔小杰李升 . 镁/铝异种金属激光焊气孔形成原因研究. 激光技术, 2013, 37(6): 825-828. doi: 10.7510/jgjs.issn.1001-3806.2013.06.027
    [16] 陈根余吴克如廖生慧邓辉张明军 . V型坡口对钢/铝激光对接熔钎焊性能的影响. 激光技术, 2014, 38(1): 11-16. doi: 10.7510/jgjs.issn.1001-3806.2014.01.003
    [17] 昝少平焦俊科王强陶俊张文武 . 镍铁铝混合粉末的激光熔覆冶金研究. 激光技术, 2018, 42(1): 131-135. doi: 10.7510/jgjs.issn.1001-3806.2018.01.026
    [18] 宋新华金湘中袁江张明军宋斌冯斌华 . 激光深熔钎焊车用铝/钢异种金属试验研究. 激光技术, 2014, 38(3): 342-346. doi: 10.7510/jgjs.issn.1001-3806.2014.03.012
    [19] 麻丁龙牛锐锋赵宏刚张俊华朱维兴王鹏飞 . Zn中间层对镁/铝异种金属激光焊接的影响. 激光技术, 2017, 41(6): 858-861. doi: 10.7510/jgjs.issn.1001-3806.2017.06.018
    [20] 周健刘双宇张福隆 . 激光自蔓延连接CFRTP/铝接头微观形貌及形成机理. 激光技术, 2019, 43(2): 147-153. doi: 10.7510/jgjs.issn.1001-3806.2019.02.001
  • 加载中
图(9)
计量
  • 文章访问数:  3952
  • HTML全文浏览量:  2766
  • PDF下载量:  22
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-09-23
  • 录用日期:  2021-11-09
  • 刊出日期:  2022-11-25

组合脉冲激光铝等离子体的数值模拟研究

    通讯作者: 廖青, liaoqing@wit.edu.cn
    作者简介: 赵佳伟(1996-), 男, 硕士研究生, 现主要从事激光诱导等离子体的研究
  • 武汉工程大学 光学信息与模式识别湖北省重点实验室, 武汉 430205
基金项目:  广东省重大基础与应用基础研究资助项目 2019B030302003

摘要: 为了获得组合脉冲激光诱导等离子体光谱性能增强的物理机制, 基于FLASH程序模拟计算了预脉冲参数对组合脉冲激光诱导的铝等离子体参数时空分布的影响, 获得了在不同预脉冲波长和不同预主脉冲延时下产生的铝等离子体电子温度、电子密度和烧蚀质量的空间演化规律。结果表明, 在预主脉冲总能量相同的情况下, 随着预脉冲波长从0.266μm变化到1.064μm, 高温环境氦气等离子体羽辉的空间范围从0.7cm增大到3.0cm, 但组合脉冲对靶材的烧蚀效率严重下降, 而铝等离子体的最大电子温度保持稳定; 此外, 组合脉冲的时间延迟低于100ns。该研究可为组合脉冲激光诱导击穿等离子体光谱增强技术提供理论参考。

English Abstract

    • 激光诱导等离子体指的是以某种样品(固体、液体等)为目标靶体,用超过特定阈值的激光辐射以得到等离子体。激光诱导等离子体包含着大量的物理过程,如融化、汽化、电离等,可应用于元素痕量分析[1]、脉冲激光沉积薄膜[2]、激光切割[3-4]等行业,还包括太空科学中的火星化学传感器[5]以及用于卫星精准定位的微推进技术[6]。经过几十年的飞速发展,单脉冲激光诱导等离子体的发展陷入瓶颈,例如,基于激光诱导等离子体的元素痕量检测限无法进一步提高。使用多脉冲技术可以通过调控等离子体的电子温度和电子密度等参数的时空分布做一些改善。

      自ZHANG等人[7]提出多脉冲激光技术以后,通过增加脉冲链来改善激光诱导等离子体性能的实验探索和理论研究不断发展,其中,对于双脉冲激光诱导等离子体的研究较为普遍。双脉冲激光诱导等离子体在除波长、脉宽等激光参数[8-10]外,需要着重考虑两脉冲之间的时间延迟、脉冲波形等。双脉冲激光诱导等离子体主要有两种激光组合方式:激光均垂直于靶面的共线脉冲以及其中一束激光平行于靶面的正交脉冲。从实用角度来说,共线脉冲是最简单的方法。对于共线脉冲[11-12]而言,预脉冲烧蚀部分样品,形成预等离子体,在一段时间延迟后,主脉冲抵达预等离子体并与之相互作用。

      2003年,SCAFFIDI等人[13]通过组合飞秒和纳秒激光,研究了双脉冲激光下激光诱导击穿光谱增强的原因。2013年,DIWAKAR等人[11]较为系统地研究了共线双纳秒激光脉冲下,脉冲间延迟及预脉冲波长对等离子体性能的影响,发现短波长的预脉冲提供相对较大的烧蚀质量、较高的等离子体电子温度和电子密度。de GIACOMO等人[14]通过实验与理论相结合的方式对比研究了单脉冲与双脉冲下激光诱导击穿光谱,发现双脉冲比单脉冲性能好的原因有几点:(1)激光对靶材更有效的烧蚀;(2)预脉冲为主脉冲创造了相对低的低压环境;(3)等离子有更高的电子温度;(4)流体力学因素的影响,单脉冲需要对环境做功,属于开放系统;而双脉冲在单脉冲营造的环境中作用,在一定程度上属于封闭系统。还有部分学者[15]关注了预主脉冲激光能量占比对双脉冲激光诱导等离子体状态参数的影响。

      由于实验测试条件的限制,目前对双脉冲激光诱导等离子体的研究大多集中在等离子体演化的中后期,对等离子体早期演化的研究较少。本文中使用FLASH程序模拟得到了双脉冲激光诱导等离子体在300ns内的等离子体参数演化图像,有助于进一步了解双脉冲激光诱导等离子体光谱性能增强的物理机制。

    • 激光诱导等离子体属于高能量密度物理过程,激光与固体靶相互作用产生的早期等离子体可视为流体,其中脉冲激光在流体动力学的能量守恒方程中被视为一个能量来源项。FLASH程序可用于模拟高能量密度物理实验过程, 以进行分析和指导实验参数设置。因为等离子体中离子温度和电子温度不相等,通常需要对等离子体进行多温度处理,此外, 热辐射效应在许多高能量密度等离子体中也很重要。对于温度的处理,部分理论研究使用的是双温度[16-17]模型:(1)电子通过逆韧致辐射吸收激光能量;(2)电子通过和离子的碰撞将能量传递给离子。FLASH程序则采用三温度模型,除了考虑等离子体流体中电子和离子之间的能量转移外,还通过追踪固定数量光子能量群中的辐射能量密度来确定辐射的扩散。

      描述激光诱导等离子体的三温流体动力学方程组可以写为:

      $ \frac{\partial \rho}{\partial t}+\nabla \cdot(\rho v)=0 $

      (1)

      $ \frac{\partial}{\partial t}(\rho v)+\nabla \cdot\left(\rho v^2\right)+\nabla p_{\mathrm{t}}=0 $

      (2)

      $ \frac{\partial}{\partial t}\left(\rho e_{\mathrm{t}}\right)+\nabla \cdot\left[\left(\rho e_{\mathrm{t}}+p_{\mathrm{t}}\right) v\right]=Q_1-\nabla \cdot q_{\mathrm{t}} $

      (3)

      式中, ρ是总质量密度; v是平均流体速度; t是时间; $ \nabla$为算子,表示对后面的量求散度; Ql表示激光加热产生的能量源; pt是电子压力pe、离子压力pi和辐射压力pr的总和:

      $ p_{\mathrm{t}}=p_{\mathrm{i}}+p_{\mathrm{e}}+p_{\mathrm{r}} $

      (4)

      比动能为v2/2, 则总比能et为:

      $ e_{\mathrm{t}}=e_{\mathrm{i}}+e_{\mathrm{e}}+e_{\mathrm{r}}+v^2 / 2 $

      (5)

      式中,ee, ei, er分别是电子、离子和辐射场的比内能。qt是总热流,被定义为电子传导qe和辐射能流qr之和:

      $ q_{\mathrm{t}}=q_{\mathrm{e}}+q_{\mathrm{r}} $

      (6)

      为了验证FLASH程序的准确性,2013年,ORBAN等人[18]分别使用FLASH和HYDRA程序进行了激光诱导等离子体的2维数值模拟,发现使用FLASH与HYDRA得到的模拟结果吻合得很好,尤其是在低密度固体靶(铝)的情况下,FLASH可以更好地再现实验结果。为了进一步验证FLASH程序是否适用于1维数值模拟,作者选取了与JOHN[19]工作相同的数值模拟参数, 将其设置为初始条件,分别使用FLASH和MED103模拟了激光辐射铝靶材的演化过程。从图 1可以看出,两者的数值模拟结果在脉冲峰值(1.5ns)前偏差较小, 但FLASH程序所使用的环境气体为具有高热导率的氦气,因此激光功率下降后(2.3ns),等离子体羽辉仍可以保持较高的电子温度。考虑到FLASH是三温度处理,MED103是一个双温度模型,因此FLASH在处理1维问题时应具有更高的精度。同时,参考BARTHÉLEMY[8]的实验数据,模拟计算出一个大气压力下激光诱导铝靶材电子温度的空间分布。如图 2所示,对于微秒量级的演化过程,FLASH仍能很好地再现实验结果。延时2μs时刻,与实验结果相比,计算得到的电子温度空间大致均匀分布,考虑到实验测量误差,FLASH数值计算结果可能更符合真实情况。延时3μs时刻,计算结果与实验测量结果基本一致。

      Figure 1.  Comparison of FLASH and MED103 numerical simulation results

      Figure 2.  Comparison between simulation data and experimental results of FLASH

      本文中所用的1维模型如图 3所示。铝靶材的初始厚度设置为20μm,模拟实验在充满环境氦气的腔室中进行,压强为660Pa。两束脉冲激光以共线的方式先后辐射在铝靶材的左表面。预脉冲作为制造预等离子体的直接能量来源,它的激光参数不仅决定了预等离子的状态而且严重影响主脉冲的再加热效果。故此,本文中选取Nd ∶YAG激光器的3个谐波充当预脉冲,波长分别为1064nm、532nm、266nm,探究预脉冲在双脉冲激光诱导等离子体中发挥的作用。其中,预脉冲脉宽为5ns,激光峰值功率密度P1=5×109W/cm2,激光通量F1=26.6J/cm2。对于主脉冲,脉宽为30ns,峰值功率密度P2=4.2×109W/cm2,波长λ=10.6μm,激光通量F2=132.8J/cm2

      Figure 3.  Theoretical model of laser induction

    • 预主脉冲延时在双脉冲激光诱导等离子体中起着重要的作用,它决定了预等离子体以何种状态与主脉冲相互作用,从而影响靶材的烧蚀效率和等离子体的再加热过程。脉冲延时过短,预等离子体未充分膨胀,达不到最佳效果。相反,过长的脉冲延时造成预等离子体过度冷却。因此,一个最佳的延时时间随实验条件而定。预主脉冲延时分别为32ns,48ns,98ns,148ns和248ns时,靶表面温度随时间的变化曲线(预脉冲波长均为1064nm)如图 4所示。由于预脉冲的参数一致,在前30ns,靶表面温度相同,随着预脉冲激光强度的变化而变化。当脉冲延时小于50ns时,靶表面温度随时间的曲线基本相同,第2个电子温度峰值出现的时间仅有微小差异。值得注意的是,脉冲间延时为98ns时,靶表面温度在主脉冲作用之前已经降到了接近室温,此时可视为靶面重新冷却为固态。随着主脉冲的到来,靶面温度回升,但极值远远小于短延时下第2个电子温度峰值, 特别当脉冲延时大于150ns时,靶面在主脉冲作用前过度冷却,且主脉冲激光对靶材的再加热效果大幅下降。

      Figure 4.  Variation of electron temperature on target surface with time under different pre-main pulse delay

      BOGAERTS等人[20]的结果显示,在100ns的脉冲延时下,靶材几乎在第2个激光脉冲到来之前已经冷却凝固。de GIACOMO等人[14]也报道了类似的行为,指出两个激光脉冲的耦合效应仅对低于100ns的脉冲延时有效,因为对于微秒量级的脉冲延时,靶材表面有足够的时间达到冷却平衡。MAO等人[21]讨论了这一点,随着脉冲延时的增加,第1个等离子体的影响逐渐变得不那么重要。

    • 为了研究预脉冲波长对等离子体的影响,固定预主脉冲延时为48ns,图 5是预脉冲波长为1064nm时电子温度Te及电子密度ne(单位是cm-3)的时空分布图。时间零点对应于激光脉冲的初始时刻,坐标零点为铝靶的左表面位置。如图 5a所示,主脉冲的能量吸收区域分布在距离靶面0.3cm外,时间起点约为60ns,最高电子温度达到109.7eV。对于波长为10.6μm的CO2激光而言,临界电子密度为1019cm-3,通过观察电子密度轮廓图可知(见图 5b),能量吸收区域对应的电子密度在临界电子密度附近,因此主脉冲的激光能量主要在此处沉积。激光结束后,在电子热传导的作用下,形成了陡峭的电子温度分布梯度:等离子体羽辉的温度逐渐均匀化且远离靶面[22],但等离子体的电子温度渐渐降低。

      Figure 5.  Spatial and temporal distribution of electron temperature and electron density at 1064nm

      此外,在电子密度时空分布图(见图 5b)中可观察到两道电子密度较高的屏障。第1道处在距离靶面较近的位置(x=0.3cm), 对应于上文中的临界电子密度,形成于48ns并持续到300ns。第2道屏障位于距靶面约1.2cm处,初步显现于150ns并在后续的时间不断增强。通过追踪靶材粒子可知,此屏障对应铝等离子体的边缘位置。由于预脉冲所用的脉宽和波长较小,可有效烧蚀靶面。在相爆炸的作用下,大量的靶材飞溅,随着主脉冲的到来,飞溅的靶材被进一步电离,并在预等离子体创造的低压环境中迅速膨胀。由于氦气的热导率较高(146.2mW·m-1·K-1),脉冲激光结束后,在热传导的作用下,温度从靶材转移到周围的缓冲气体上,缓冲气体在高温状态下迅速膨胀得以反向挤压等离子体,形成了这道等离子体屏障。

      图 6图 7分别为预脉冲波长为530nm和266nm时,电子温度和电子密度的时空分布图。预脉冲波长为530nm时,烧蚀动力学发生了些许变化。对于电子温度而言(见图 6a),不仅最大电子温度降低到106.0eV,而且高温氦气等离子体的范围也有略微缩小,但总体看来,等离子体及羽辉的时空分布基本保持一致。预脉冲波长为266nm时,最大电子温度进一步降低到101.9eV(如图 7a所示),但更值得注意的是,高温羽辉的宽度大幅度降低。与预脉冲波长为1064nm的结果相比,200ns时等离子体的羽辉宽度从3.0cm缩减到了0.7cm,缩减幅度约为3倍,而核心区域的电子温度仅下降约10eV。此外必须指出,尽管预脉冲的波长不同,但300ns时刻铝等离子体的空间范围基本相同。

      Figure 6.  Spatial and temporal distribution of electron temperature and electron density at 530nm

      Figure 7.  Spatial and temporal distribution of electron temperature and electron density at 266nm

      对于不同预脉冲波长的等离子体参数时空分布而言,其中最为显著的是电子温度的变化。虽然改变的是预脉冲的波长,但在预脉冲阶段,电子温度并未发生明显差异,反而在主脉冲阶段有较为明显的改变。随着预脉冲波长的增加,等离子体的最高电子温度有轻微提升,且高温羽辉的范围也有所扩大。这些数值模拟结果可以从能量吸收机制的角度来解释,预脉冲主要靠多光子电离吸收能量,但对主脉冲激光,逆韧致辐射在能量吸收机制中起主要作用。在FLASH程序中,通过逆轫致辐射(inverse bremsstrahlung,IB)频率因子νIB计算逆韧致辐射过程:

      $ \nu_{\mathrm{IB}}=\frac{n_{\mathrm{e}}}{n_{\mathrm{c}}} \nu_{\mathrm{e}, \mathrm{i}} $

      (7)

      式中, ne是电子密度, nc表示临界电子密度,νe, i是电子离子碰撞频率:

      $ \nu_{\mathrm{e}, \mathrm{i}}=\frac{4}{3}\left(\frac{2 {\rm{ \mathsf{ π} }}}{m_{\mathrm{e}}}\right)^{1 / 2} \frac{n_{\mathrm{e}} Z e^4 \ln \mathit{\Lambda}}{\left(k_{\mathrm{B}} T_{\mathrm{e}}\right)^{1 / 2}} $

      (8)

      式中,me是电子质量,Z是等离子体的平均电离数,e是电子电荷,lnΛ是库仑对数,Λ为德拜数,kB是玻尔兹曼常数, Te是电子温度。因此,逆轫致辐射系数取决于电子密度和电子温度以及激光波长:

      $ \nu_{\mathrm{IB}}=\frac{4 \sqrt{2} Z e^6}{{\rm{ \mathsf{ π} }}^{1 / 2} m_{\mathrm{e}}{ }^{3 / 2} k_{\mathrm{B}}{ }^{3 / 2} c^2} \lambda^2 n_{\mathrm{e}}{ }^2 T_{\mathrm{e}}^{3 / 2} \ln \mathit{\Lambda} $

      (9)

      式中,c是真空中的波速,λ为激光的波长。

      从(9)式可以看出,逆轫致辐射系数随着激光波长、电子密度和电子温度的增大而增大,在本文中,就主脉冲峰值时刻而言,图 5所示的预脉冲波长、电子温度及电子密度均高于其它方案,因此νIB系数越大,预等离子体对主脉冲的能量吸收率也越大,解释了图 5~图 7中双脉冲激光等离子体参数随预脉冲波长的变化规律。

    • 共线双脉冲激光诱导等离子体必须考虑激光对靶材的烧蚀效果,尤其是主脉冲对靶材的再加热效果。本节中通过研究靶材及其附近质量密度的变化来探索不同的预脉冲波长下靶材的烧蚀情况。预脉冲波长为1064nm时,质量密度分布如图 8a所示。初期,预脉冲烧蚀铝靶材表面,少量靶材飞溅,同时产生的冲击波挤压靶材到更高的质量密度。在主脉冲峰值时刻,膨胀的等离子体产生压力将一部分烧蚀材料推向靶材,融化的靶材质量密度再次升高。激光辐照结束后,在热传导的作用下,高温的等离子体持续融化靶材。预脉冲波长变短后,波长为530nm时的靶材质量密度分布图如图 8b所示。当预脉冲波长为266nm时,情况发生了明显变化(见图 8c),预脉冲作用于铝靶材的融化效果大于冲击波的压缩能力,靶材的质量密度ρ(单位为g·cm-3)初步下降,当主脉冲抵达时,靶材质量密度进一步降低形成铝蒸气。

      Figure 8.  Mass density distribution of target material at different pre-pulse wavelengths

      尽管预脉冲波长影响烧蚀出的靶材质量,但预脉冲激光创造的等离子体屏障才是影响主脉冲与靶材能量耦合效率的重要因素。因此,计算出36ns时激光能量沉积的空间位置。如图 9所示,预脉冲波长为266nm时,激光能量沉积在靠近靶面的位置相对强烈,而预脉冲波长为1064nm时,激光能量沉积的峰值位置在远离靶面的同时, 强度降低到前者的约1/4。分析可知:长的预脉冲波长,压缩靶材的同时部分靶材产生飞溅,远离靶材表面的电子密度临界面对主脉冲产生屏蔽效应,则主脉冲主要与预烧蚀的等离子体相互作用,从而导致等离子体的温度提高;短的预脉冲波长,等离子体电子密度临界面靠近靶材表面,主脉冲可以通过热传导与靶材进行较高效的能量耦合,靶材吸收了部分能量进一步融化,因此远离靶面的等离子体的电子温度相对降低。

      Figure 9.  Spatial distribution of laser energy deposition at different pre-pulse wavelengths

    • 使用FLASH程序包模拟了一维纳秒组合双脉冲激光辐射铝靶材的演化过程。研究发现,组合脉冲的预主脉冲延应低于100ns,以保障主脉冲的再加热效应。通过电子密度、电子温度的时空变化确定了逆韧致辐射在能量吸收机制中的主导地位:预脉冲波长越长,等离子体有较高的电子温度和范围更大的羽辉。从电子密度时空分布图可看出,预等离子体的屏障作用对主脉冲的再加热效果具有重要影响:预脉冲波长越长,预等离子体电子临界密度越低,等离子体屏蔽效应越强且位置越远离靶面。通过计算烧蚀质量发现:长的预脉冲波长,因为较强的屏蔽效应减弱了主脉冲与靶材的耦合效率,故靶材的融化效果较差。

参考文献 (22)

目录

    /

    返回文章
    返回