搜索

文章查询

x

留言板

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

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

束匀滑光束偏折现象的模拟

李斌 刘占军 郝亮 郑春阳 蔡洪波 何民卿

束匀滑光束偏折现象的模拟

李斌, 刘占军, 郝亮, 郑春阳, 蔡洪波, 何民卿
PDF
HTML
导出引用
导出核心图
  • 达到高能量密度物理状态后, 光束在介质中的传播行为与经典光学研究范畴相比, 会出现一些新现象. 比如在各向同性介质内可出现光束传播方向改变的现象. 另一方面, 高能量密度物理实验中由高功率激光器产生的束匀滑光束较为常见. 本文分析了空间和时间束匀滑光束在各向同性等离子体传播中出现束偏折现象的机制和条件, 并利用三维激光等离子体相互作用程序LAP3D进行了验证. 模拟表明只有当同时发生成丝不稳定性和存在离子声速量级的横向流时束匀滑激光才会产生显著的束偏折现象.
      通信作者: 刘占军, liuzj@iapcm.ac.cn ; 郝亮, hao_liang@iapcm.ac.cn
    • 基金项目: 国家级-激光等离子体次级不稳定性及其耦合过程研究(11875093)
    [1]

    Lindl J D, Amendt P, Berger R L, Glendinning S G, Glenzer S H, Haan S W, Kauffman R L, Landen O L, Suter L J 2004 Phys. Plasmas 11 339

    [2]

    德雷克 著 (孙承纬 译) 2013 高能量密度物理基础、惯性约束聚变和实验天体物理学(北京: 国防工业出版社) 第245—286页

    Drake R P (translated by Sun C W) 2013 High-Energy-Density Physics Fundamentals, Inertial Fusion, and Experimental Astrophysics (Beijing: National Defense Industry Press) pp245−286 (in Chinese)

    [3]

    Montgomery D S 2016 Phys. Plasmas 23 055601

    [4]

    Kirkwood R K, Moody J D, Kline J, Dewald E, Glenzer S, Divol L, Michel P, Hinkel D, Berger R, Williams E, Milovich J, Yin L, Rose H, MacGowan B, Landen O, Rosen M, Lindl J 2013 Plasma Phys. Control. Fusion 55 103001

    [5]

    张家泰 1999 激光等离子体相互作用物理与模拟(郑州: 河南科学技术出版社) pp269−316

    Zhang J T 1999 Physics and Simulations of Laser Plasma Interactions (Zhengzhou: Henan Science and Technology Press) 第269—316页 (in Chinese)

    [6]

    Moody J D, MacGowan B J, Rothenberg J E, Berger R L, Young P E 2001 Phys. Rev. Lett. 86 2810

    [7]

    Craxton R S, Anderson K S, Boehly T R, Goncharov V N, Harding D R, Knauer J P, McCrory R L, McKenty P W, Meyerhofer D D, Myatt J F, Schmitt A J, Sethian J D, Short R W, Skupsky S, Theobald W, Kruer W L, Tanaka K, Betti R, Collins T J B, Delettrez J A, Hu S X, Marozas J A, Maximov A V, Michel D T, Radha P B, Regan S P, Sangster T C, Seka W, Solodov A A, Soures J M, Stoeckl C, Zuegel J D 2015 Phys. Plasmas 22 110501

    [8]

    周煜梁, 隋展, 刘兰琴, 粟敬钦, 李平, 张锐, 许立新, 王文义, 莫磊 2011 激光与光电子学进展 48 101407

    Zhou Y L, Sui Z, Liu L Q, Su J Q, Li P, Zhang R, Xu L X, Wang W Y, Mo L 2011 Laser & Optoelectronics Progress 48 101407

    [9]

    Young P E, Still C H, Hinkel D E, Kruer W L, Estabrook K G 1998 Phys. Rev. Lett. 81 1425

    [10]

    Moody J D, MacGowan B J, Hinkel D E, Kruer W L, Williams E A, Estabrook K, Berger R L, Kirkwood R K, Montgomery D S, Shepard T D 1996 Phys. Rev. Lett. 77 1294

    [11]

    李斌, 刘占军, 郑春阳, 胡晓燕 2014 强激光与离子束 26 122005

    Li B, Liu Z J, Zheng C Y, Hu X Y 2014 High Power Laser and Particle Beams 26 122005

    [12]

    Hu X Y, Hao L, Liu Z J, Zheng C Y, Li B, Guo H 2015 AIP Advances 5 087174

    [13]

    Liu Z J, Li B, Hu X Y, Xiang J, Zheng C Y, Cao L H, Hao L 2016 Phys. Plasmas 23 022705

    [14]

    李平, 马驰, 粟敬钦, 程文雍, 刘兰琴, 王文义, 莫磊, 周丽丹 2008 强激光与粒子束 20 1114

    Li P, Ma C, Su J Q, Cheng W Y, Liu L Q, Wang W Y, Mo L, Zhou L D 2008 High Power Laser and Particle Beams 20 1114

    [15]

    李斌, 胡晓燕, 郑春阳, 刘占军 2016 强激光与离子束 28 112004

    Li B, Hu X Y, Zheng C Y, Liu Z J 2016 High Power Laser and Particle Beams 28 112004

    [16]

    李斌, 刘占军, 郝亮, 胡晓燕, 郑春阳, 项江 2017 中国激光 44 1201004

    Li B, Liu Z J, Hao L, Hu X Y, Zheng C Y, Xiang J 2017 Chinese J. Lasers 44 1201004

    [17]

    Kruer W L 2003 The Physics of Laser Plasma Interactions (Colorado: Westview Press) pp70–71

    [18]

    Tomson J J, Kruer W L, Bodner S E, DeGroot J S 1974 Phys. Fluids 17 849

    [19]

    Valeo E J, Oberman C 1973 Phys. Rev. Lett. 30 1035

    [20]

    Thomson J J, Karush J I 1974 Phys. Fluids 17 1608

    [21]

    Williams E A, Albitton J R, Rosenbluth M N 1979 Phys. Fluids 22 139

  • 图 1  不同激光强度下空间束匀滑光束对应的束偏折模拟结果 (a) Φ200模型2加横向流; (b) Φ200模型5加横向流. 图中横纵坐标对应模拟空间坐标zx, 其量纲为激光波长. 横向流速等于离子声速

    Fig. 1.  Beam deflection simulation results at different incident intensity: (a) Transverse flow and average intensity lower than filamentation threshold; (b) transverse flow and average intensity higher than filamentation threshold. x and yaxes of two figures corresponding to xand z axes of simulation coordinates, respectively. The spatial scale is in unit of laser wave length. The transverse flow speed equals ion sound speed.

    图 2  对比空间束匀滑光束发生束偏折时Φ200模型5加横向流的入射面和出射面内光斑电场幅值分布 (a)入射面; (b)出射面. 图中横纵坐标对应模拟空间坐标xy, 其量纲为激光波长. 横向流速等于离子声速

    Fig. 2.  Comparison of spatial distribution of laser electric field between laser entrance and exit planes as beam deflection presents: (a) Laser entrance plane; (b) laser exit plane. x and y axes of two figures corresponding to x and y axes of simulation coordinates, respectively. The spatial scale is in unit of laser wave length. The transverse flow speed equals ion sound speed.

    图 3  调制频率为${10^{ - 3}}{\omega _0}$的时间束匀滑光束传播行为 (a)对应11000激光周期; (b)对应13750激光周期. 图中横纵坐标对应模拟空间坐标zy, 其量纲为激光波长

    Fig. 3.  Propagation of SSD beam at modulation frequency of 10–3ω0: (a) Corresponding simulation result at 11000 th laser periods; (b) corresponding simulation result at 13750 th laser periods. x and y axes of two figures corresponding to z and y axes of simulation coordinates, respectively. The spatial scale is in unit of laser wave length.

    图 4  有横向离子声速量级等离子体流时调制频率为${10^{ - 3}}{\omega _0}$的时间束匀滑光束的传播行为 (a)对应11000激光周期; (b)对应13750激光周期. 图中横纵坐标对应模拟空间坐标zy, 其量纲为激光波长. 横向流速等于离子声速

    Fig. 4.  Propagation of SSD beam with transverse flow at modulation frequency of 10–3 ω0: (a) Corresponding simulation result at 11000 th laser periods; (b) corresponding simulation result at 13750th laser periods. x and y axes of two figures corresponding to z and y axes of simulation coordinates, respectively. The spatial scale is in unit of laser wave length. The transverse flow speed equals ion sound speed.

    图 5  对比时间束匀滑光束在调制频率为${10^{ - 4}}{\rm{ }}{\omega _0}$时的光束传播行为 (a)等离子体横向流速为零; (b)等离子体横向流速等于离子声速. 图中横纵坐标对应模拟空间坐标yz, 其量纲为激光波长

    Fig. 5.  Propagation of SSD beam at modulation frequency of 10–4ω0: (a) No transverse flow; (b) the transverse flow speed equals ion sound speed. x and y axes of two figures corresponding to y and z axes of simulation coordinates, respectively. The spatial scale is in unit of laser wave length.

    表 1  空间束匀滑光束Φ200成丝和束偏折现象模拟结果

    Table 1.  Simulation results for filamentaion and beam deflection in the case of CPP smoothed beam Φ200.

    模型光斑平均强度/${\rm{W} } \cdot {\rm{c} }{ {\rm{m} }^{ {\rm{ - 2} } } }$成丝现象束偏折现象
    14.30 × 1013
    23.86 × 1014
    31.07 × 1015
    42.11 × 1015
    53.49 × 1015
    下载: 导出CSV
  • [1]

    Lindl J D, Amendt P, Berger R L, Glendinning S G, Glenzer S H, Haan S W, Kauffman R L, Landen O L, Suter L J 2004 Phys. Plasmas 11 339

    [2]

    德雷克 著 (孙承纬 译) 2013 高能量密度物理基础、惯性约束聚变和实验天体物理学(北京: 国防工业出版社) 第245—286页

    Drake R P (translated by Sun C W) 2013 High-Energy-Density Physics Fundamentals, Inertial Fusion, and Experimental Astrophysics (Beijing: National Defense Industry Press) pp245−286 (in Chinese)

    [3]

    Montgomery D S 2016 Phys. Plasmas 23 055601

    [4]

    Kirkwood R K, Moody J D, Kline J, Dewald E, Glenzer S, Divol L, Michel P, Hinkel D, Berger R, Williams E, Milovich J, Yin L, Rose H, MacGowan B, Landen O, Rosen M, Lindl J 2013 Plasma Phys. Control. Fusion 55 103001

    [5]

    张家泰 1999 激光等离子体相互作用物理与模拟(郑州: 河南科学技术出版社) pp269−316

    Zhang J T 1999 Physics and Simulations of Laser Plasma Interactions (Zhengzhou: Henan Science and Technology Press) 第269—316页 (in Chinese)

    [6]

    Moody J D, MacGowan B J, Rothenberg J E, Berger R L, Young P E 2001 Phys. Rev. Lett. 86 2810

    [7]

    Craxton R S, Anderson K S, Boehly T R, Goncharov V N, Harding D R, Knauer J P, McCrory R L, McKenty P W, Meyerhofer D D, Myatt J F, Schmitt A J, Sethian J D, Short R W, Skupsky S, Theobald W, Kruer W L, Tanaka K, Betti R, Collins T J B, Delettrez J A, Hu S X, Marozas J A, Maximov A V, Michel D T, Radha P B, Regan S P, Sangster T C, Seka W, Solodov A A, Soures J M, Stoeckl C, Zuegel J D 2015 Phys. Plasmas 22 110501

    [8]

    周煜梁, 隋展, 刘兰琴, 粟敬钦, 李平, 张锐, 许立新, 王文义, 莫磊 2011 激光与光电子学进展 48 101407

    Zhou Y L, Sui Z, Liu L Q, Su J Q, Li P, Zhang R, Xu L X, Wang W Y, Mo L 2011 Laser & Optoelectronics Progress 48 101407

    [9]

    Young P E, Still C H, Hinkel D E, Kruer W L, Estabrook K G 1998 Phys. Rev. Lett. 81 1425

    [10]

    Moody J D, MacGowan B J, Hinkel D E, Kruer W L, Williams E A, Estabrook K, Berger R L, Kirkwood R K, Montgomery D S, Shepard T D 1996 Phys. Rev. Lett. 77 1294

    [11]

    李斌, 刘占军, 郑春阳, 胡晓燕 2014 强激光与离子束 26 122005

    Li B, Liu Z J, Zheng C Y, Hu X Y 2014 High Power Laser and Particle Beams 26 122005

    [12]

    Hu X Y, Hao L, Liu Z J, Zheng C Y, Li B, Guo H 2015 AIP Advances 5 087174

    [13]

    Liu Z J, Li B, Hu X Y, Xiang J, Zheng C Y, Cao L H, Hao L 2016 Phys. Plasmas 23 022705

    [14]

    李平, 马驰, 粟敬钦, 程文雍, 刘兰琴, 王文义, 莫磊, 周丽丹 2008 强激光与粒子束 20 1114

    Li P, Ma C, Su J Q, Cheng W Y, Liu L Q, Wang W Y, Mo L, Zhou L D 2008 High Power Laser and Particle Beams 20 1114

    [15]

    李斌, 胡晓燕, 郑春阳, 刘占军 2016 强激光与离子束 28 112004

    Li B, Hu X Y, Zheng C Y, Liu Z J 2016 High Power Laser and Particle Beams 28 112004

    [16]

    李斌, 刘占军, 郝亮, 胡晓燕, 郑春阳, 项江 2017 中国激光 44 1201004

    Li B, Liu Z J, Hao L, Hu X Y, Zheng C Y, Xiang J 2017 Chinese J. Lasers 44 1201004

    [17]

    Kruer W L 2003 The Physics of Laser Plasma Interactions (Colorado: Westview Press) pp70–71

    [18]

    Tomson J J, Kruer W L, Bodner S E, DeGroot J S 1974 Phys. Fluids 17 849

    [19]

    Valeo E J, Oberman C 1973 Phys. Rev. Lett. 30 1035

    [20]

    Thomson J J, Karush J I 1974 Phys. Fluids 17 1608

    [21]

    Williams E A, Albitton J R, Rosenbluth M N 1979 Phys. Fluids 22 139

  • [1] 晏骥, 郑建华, 陈黎, 林稚伟, 江少恩. X射线相衬成像技术应用于高能量密度物理条件下内爆靶丸诊断. 物理学报, 2012, 61(14): 148701. doi: 10.7498/aps.61.148701
    [2] 裴晓星, 仲佳勇, 张凯, 郑无敌, 梁贵云, 王菲鹿, 李玉同, 赵刚. 实验室天体物理的验证特例:W43A磁喷流. 物理学报, 2014, 63(14): 145201. doi: 10.7498/aps.63.145201
    [3] 周冰洁, 钟哲强, 张彬. 光束运动特性对焦斑束匀滑效果的影响. 物理学报, 2012, 61(21): 214202. doi: 10.7498/aps.61.214202
    [4] 晏骥, 韦敏习, 蒲昱东, 刘慎业, 詹夏雨, 林稚伟, 郑建华, 江少恩. 新型针孔点背光发光模型与实验研究. 物理学报, 2013, 62(1): 015204. doi: 10.7498/aps.62.015204
    [5] 晏骥, 郑建华, 陈黎, 涂绍勇, 韦敏习, 余波, 刘慎业, 江少恩. 基于神光Ⅲ原型装置的新型针孔点背光实验. 物理学报, 2013, 62(4): 045203. doi: 10.7498/aps.62.045203
    [6] 张扬, 戴自换, 孙奇志, 章征伟, 孙海权, 王裴, 丁宁, 薛创, 王冠琼, 沈智军, 李肖, 王建国. FP-1装置铝套筒内爆动力学过程的一维磁流体力学模拟. 物理学报, 2018, 67(8): 080701. doi: 10.7498/aps.67.20172300
    [7] 尹传磊, 王伟民, 廖国前, 李梦超, 李玉同, 张杰. 超强圆偏振激光直接加速产生超高能量电子束. 物理学报, 2015, 64(14): 144102. doi: 10.7498/aps.64.144102
    [8] 刘兰琴, 张颖, 耿远超, 王文义, 朱启华, 景峰, 魏晓峰, 黄晚晴. 小宽带光谱色散匀滑光束传输特性研究. 物理学报, 2014, 63(16): 164201. doi: 10.7498/aps.63.164201
    [9] 钟哲强, 侯鹏程, 张彬. 基于光克尔效应的径向光束匀滑新方案. 物理学报, 2016, 65(9): 094207. doi: 10.7498/aps.65.094207
    [10] 田博宇, 钟哲强, 隋展, 张彬, 袁孝. 基于涡旋光束的超快速角向集束匀滑方案. 物理学报, 2019, 68(2): 024207. doi: 10.7498/aps.68.20181361
    [11] 钟哲强, 周冰洁, 叶荣, 张彬. 多频多色光谱角色散束匀滑新方案. 物理学报, 2014, 63(3): 035201. doi: 10.7498/aps.63.035201
    [12] 李福建, 高妍琦, 赵晓晖, 季来林, 王伟, 黄秀光, 马伟新, 隋展, 裴文兵. 诱导空间非相干束匀滑技术的近区特性及改善技术. 物理学报, 2018, 67(17): 175201. doi: 10.7498/aps.67.20180533
    [13] 李腾飞, 钟哲强, 张彬. 用于超快束匀滑的动态波前调控新方案. 物理学报, 2018, 67(17): 174206. doi: 10.7498/aps.67.20172527
    [14] 熊皓, 钟哲强, 张彬, 隋展, 张小民. 基于束间动态干涉的快速匀滑新方法. 物理学报, 2020, 69(6): 064206. doi: 10.7498/aps.69.20190962
    [15] 高妍琦, 赵晓晖, 贾果, 李福建, 崔勇, 饶大幸, 季来林, 刘栋, 冯伟, 黄秀光, 马伟新, 隋展. 基于低相干光的阵列透镜束匀滑技术研究. 物理学报, 2019, 68(7): 075201. doi: 10.7498/aps.68.20182138
    [16] 温圣林, 姚欣, 高福华, 张怡霄, 郭永康, 林祥棣. 激光惯性约束聚变驱动器终端光学系统中束匀滑器件前置的条件研究. 物理学报, 2009, 58(5): 3130-3134. doi: 10.7498/aps.58.3130
    [17] 姚欣, 高福华, 高博, 张怡霄, 黄利新, 郭永康, 林祥棣. 惯性约束聚变驱动器终端束匀滑器件前置时频率转换系统优化研究. 物理学报, 2009, 58(7): 4598-4604. doi: 10.7498/aps.58.4598
    [18] 邹其徽, 吕百达. 等束宽超短脉冲光束的远场特性. 物理学报, 2005, 54(12): 5642-5647. doi: 10.7498/aps.54.5642
    [19] 沈晓鹏, 韩 奎, 李海鹏, 沈义峰, 王子煜. 光子晶体自准直光束偏振分束器. 物理学报, 2008, 57(3): 1737-1741. doi: 10.7498/aps.57.1737
    [20] 肖玲, 程小劲, 徐剑秋. 分数自成像平面波导的光束组束. 物理学报, 2009, 58(6): 3870-3876. doi: 10.7498/aps.58.3870
  • 引用本文:
    Citation:
计量
  • 文章访问数:  288
  • PDF下载量:  16
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-10-28
  • 修回日期:  2020-02-01
  • 刊出日期:  2020-04-01

束匀滑光束偏折现象的模拟

    基金项目: 国家级-激光等离子体次级不稳定性及其耦合过程研究(11875093)

摘要: 达到高能量密度物理状态后, 光束在介质中的传播行为与经典光学研究范畴相比, 会出现一些新现象. 比如在各向同性介质内可出现光束传播方向改变的现象. 另一方面, 高能量密度物理实验中由高功率激光器产生的束匀滑光束较为常见. 本文分析了空间和时间束匀滑光束在各向同性等离子体传播中出现束偏折现象的机制和条件, 并利用三维激光等离子体相互作用程序LAP3D进行了验证. 模拟表明只有当同时发生成丝不稳定性和存在离子声速量级的横向流时束匀滑激光才会产生显著的束偏折现象.

English Abstract

    • 在经典光学范围内, 光束在各向同性介质内的传播方向不会发生改变, 传播方向改变的现象通常发生在不同介质的交界面处. 但当能量驱动源如激光, 使介质能量增加并达到高能量密度状态后, 激光强度足以改变介质折射率, 在这种介质内也会发生光束传播方向改变的现象(束偏折)[1]. 另一方面在高能量密度物理研究中, 随着近年来高功率激光器技术的发展, 采用激光束辐照产生高能量密度状态成为一种常用的方式. 许多采用大型高功率激光器的实验, 通常输出束匀滑光束[2-8]. 针对高能量密度状态下, 束匀滑光束发生束偏折的机制和条件, 本文给出了相关研究结果.

      束匀滑光束按光束匀滑方式可分为空间束匀滑和时间束匀滑两种. 它们的共同特点是光束包络较为平滑, 但内部结构较为复杂, 内部包含数千个小光斑, 纵向长度为$8{F^2}{\lambda _0}$[1,2], 横向长度为$2 F{\lambda _0}$[1,2], 这里F为透镜光学F数, ${\lambda _0}$为激光波长. 分析束匀滑光束发生束偏折的条件, 要从光束内部多个小光束相互作用角度考虑, 涉及到多体相互作用, 仅应用理论分析较为困难. 如果不考虑光束内部结构, 仅从整个光束角度考虑, 其发生束偏折的条件在一定程度上可类比单个光束, 如高斯型光束发生束偏折的机制和条件, 本文从这个角度展开了对束匀滑光束发生偏折现象所需条件的分析和研究.

      实验研究[9,10]和近期利用LAP3D程序对高斯型光束的模拟研究[11]都表明, 高斯型激光束发生显著偏折需满足两个条件: 发生成丝不稳定性和存在声速量级的等离子体横向流[11]. 前者在聚焦光斑处形成密度凹坑并改变此区域折射率, 后者使这种密度空间结构能在横向整体移动. 据此认为研究束匀滑光束发生束偏折的条件可进一步细化为研究以下两点: 束匀滑光束发生成丝不稳定性的条件以及声速量级的等离子体横向流对成丝不稳定性的影响.

      本文根据以上两个条件对束匀滑光束偏折现象的产生机制和条件进行了研究, 给出了用于模拟三维束匀滑光束偏折现象的数值模拟程序LAP3D的物理建模, 分析了横向流大小对不稳定性发展的影响, 给出了研究空间束匀滑光束和时间束匀滑光束发生束偏折条件的数值模拟结果和分析.

    • 为描述光束在等离子体中的传输和等离子体响应, 三维激光等离子体相互作用程序LAP3D分别采用了激光包络方程和双温流体方程组[12,13]的物理模型. 由于束偏折现象发生与否跟光束成丝现象紧密相关, 因此在激光传播模型中只考虑衍射项和折射项, 忽略激光受激散射项, 对应的包络化激光传播方程为

      $\left( {\frac{\partial }{{\partial t}} + {v_{\rm{g}}}\frac{\partial }{{\partial z}} - \frac{{{\rm{i}}{c^2}}}{{2{\omega _0}}}D_ \bot ^2 + \nu } \right)E = \frac{{ - {\rm{i}}2{\text{π}}{{\rm{e}}^2}}}{{{m_{\rm e}}{\omega _0}}}\delta {n_{\rm e{\rm{f}}}}E,$

      其中, E为激光电场强度; ${v_{\rm{g}}}$ 为激光波包传播群速度; c为真空中激光光速; ${\omega _0}$为激光频率; $\nu $为逆轫致吸收率; ${m_{\rm e}}$为电子质量; $\delta {n_{\rm e{\rm{f}}}}$为电子扰动数密度; $D_ \bot ^2$为广义衍射算子[12], 在下文推导中为简化用$\nabla _ \bot ^2$来代替.

    • 从前期对单个高斯光束偏折的模拟研究[11]发现, 产生束偏折的空间位置与离子声速大小有以下关系: 当横向流为离子声速时, 发生束偏折时的光束传播距离最短; 当横向流速逐渐小于离子声速时, 发生束偏折的位置也逐渐“靠后”, 即光束传播距离逐渐增加. 当横向流速远大于离子声速后, 束偏折现象消失. 通过下面分析可知, 这实际上是横向流速大小对成丝不稳定性发展影响的一种体现.

      设沿激光传播方向为z方向, 存在横向(沿x方向)流速为${u_ \bot }$的等离子体流场, 考虑流体方程组中的质量连续方程:

      ${\partial \rho }/\partial t + \nabla \cdot ({ u}) = 0, $

      和动量方程:

      $ \dfrac{\partial { S}}{\partial t} + \nabla \cdot ({ {S u}}) + \nabla Q = - \nabla (p + p_{\rm e}) - \rho \nabla \psi , $

      其中, ${{S}}$为电子动量, ${{u}}$为流速, ${{Q}}$为人为粘性应力张量, $\psi $为有质动力势, p${p_e}$为离子压力和电子压力.

      对上面两个方程线性化, 设$E = {E_0} + \delta E$, $n = {n_0} + \delta {{{n}}_{\rm ef}}$, 则质量连续方程和动量方程经化简后, 可得$\delta {n_{e{\rm{f}}}}$满足的密度涨落方程为

      $\begin{split}&\left( {\frac{{{\partial ^2}}}{{\partial {t^2}}} + u_ \bot ^2\frac{\partial }{{\partial {x^2}}} - c_{\rm{s}}^2{\nabla ^2}} \right)\delta {n_{\rm e{\rm{f}}}}\\ =& \;\frac{{Z{n_{\rm e}}{e^2}}}{{mM{\omega _0}}}{\nabla ^2}\left( {{E_0} \cdot \delta E} \right),\end{split}$

      其中, ${c_{\rm{s}}}$为离子声速, M为离子质量, 对(1)式和(4)式做傅里叶分析可得色散关系为

      $\begin{split} &{\omega ^2} - k_ \bot ^2\left( {c_{\rm{s}}^2 - u_ \bot ^2} \right) = \frac{{k_ \bot ^2v_{{\rm{os}}}^2\omega _{{\rm{pi}}}^2}}{4}\\\; & \times\left[ {\frac{1}{{D\left( {\omega \!-\! {\omega _0},k \!-\! {k_0}} \right)}} \!+\! \frac{1}{{D\left( {\omega \!+\! {\omega _0},k \!+\! {k_0}} \right)}}} \right].\end{split}$

      考虑$\omega = {\rm{i}}\gamma \ll {\omega _0}$并且${k} \cdot { k} _0 = 0$, 化简得到光束传播过程中存在横向小扰动及有横向流场发生成丝不稳定性对应的空间增长率:

      $\frac{{{\gamma ^2}}}{{k_ \bot ^2}} = {{\frac{1}{8}\frac{{{n_{\rm e}}}}{{{n_{{\rm{cr}}}}}}{{\left( {\frac{{{v_{{\rm{os}}}}}}{{{v_{\rm e}}}}} \right)}^2}k_ \bot ^2}\bigg/{\left( {1 - \frac{{u_ \bot ^2}}{{c_{\rm{s}}^2}}} \right)}} - \frac{{k_ \bot ^4}}{{4k_0^2}}, $

      式中, $\gamma $为成丝增长率, ${k_ \bot }$为横向扰动波数, ${k_0}$为入射光束波数, ${v_{{\rm{os}}}}$${v_{\rm e}}$分别为电子在光场中的抖动速度和热速度.

      对于在各向同性等离子体中传播的高斯光束, 可认为在给定参数下(6)式中物理量${v_{\rm e}}$${v_{{\rm{os}}}}$${n_{\rm e}}$${n_{{\rm{cr}}}}$${k_0}$${c_{\rm{s}}}$均为定值. 在光束的扰动波数${k_ \bot }$范围内, 以${k_{{\rm{opt}}}}$为例考察(6)式中横向流大小对不稳定性发展的影响. 这里${k_{{\rm{opt}}}}$为没有横向流时, 在此条件下使成丝不稳定性增长率取最大值时对应的扰动波数. 进一步分析(6)式可知, 当横向流速为离子声速时, 对应不稳定增长率为极大值; 当流速逐渐小于离子声速时, 不稳定性增长率逐渐减小, 发生成丝的空间位置也逐渐“靠后”; 当流速大于离子声速时, 成丝发展受抑制.

      结合本节第一段所描述的模拟结果[11], 上面分析表明, 对高斯型光束, 横向流对其束偏折的影响本质上是横向流对成丝不稳定性发展的影响. 考虑到束匀滑光束包络内部包含多个小焦斑, 其光强分布可近似等同于高斯光束, 因此可推断束匀滑光束发生显著的束偏折现象也需满足发生成丝不稳定性和存在声速量级横向流的条件. 本文下面两节将采用数值模拟研究的方法来进一步分析和判断.

    • 激光器产生的入射光束通过位于聚焦透镜前的相位板后, 被分割成多个子光束, 在远场叠加形成空间束匀滑光束. 本文模拟中应用连续相位板(continuous phase plate, CPP), 它使子光束间相位差在0到${\text{π}}$之间连续变化, 对比随机相位板, 不仅在远场可产生更为平滑的包络分布, 也使焦斑中能量提高到入射光能量的95%[14].

      焦平面内经CPP后生成的空间束匀滑光束的电场幅值分布可表示为下式:

      $\begin{split} E(x',y') =\; & \frac{1}{{{\lambda _0}f}}\int E(x,y)t(x,y)\\ & \times\exp \left[ { - {\rm{i}}\frac{{2{\text{π}}}}{{\lambda f}}\left( {xx' + yy'} \right)} \right] {\rm{d}}x{\rm{d}}y,\end{split}$

      其中, f为透镜焦距, $E(x, y)$为入射光复振幅, $t\left( {x, y} \right)$为CPP的相位信息.

      LAP3D程序采用笛卡尔坐标系, 模拟三维大尺度上发生的空间束匀滑光束偏折现象. 模拟中激光从(x, y)面入射, 沿z方向传播. (x, y, z)方向上对应的模拟尺度分别为$\left( {2048{\lambda _0}, 2048{\lambda _0}, 2048{\lambda _0}} \right)$, 单位为激光波长${\lambda _0}$, ${\lambda _0} = $ 0.351 μm. 模拟空间内等离子体均匀分布, 初始等离子体密度为$0.{\rm{1 }}{n_{\rm{c}}}$, 电子温度${T_e}$$3\;{\rm{ keV}}$, 离子温度${T_{\rm{i}}}$$1\;{\rm{ keV}}$, 离子电离度$Z = 1$.

      模拟中采用光学F数为11, 长轴约为200 μm, 短轴约为100 μm的椭圆焦斑Φ200, 椭圆焦斑内的电场幅值分布如文献[15]中图1所示. 根据空间束匀滑光束入射面平均光强大小, 设计了5个模型. 首先模拟了在没有横向等离子体流的条件下, 空间束匀滑光束在均匀等离子体中的传播行为, 重点考察是否发生了光束成丝现象. 然后引入声速量级的横向流, 再次计算上述5个模型, 重点考察是否发生了光束偏折现象. 等离子体横向流由x方向左边界流入, 从x方向右边界流出, 垂直于激光传播方向. 这些模拟结果如表1所示.

      图  1  不同激光强度下空间束匀滑光束对应的束偏折模拟结果 (a) Φ200模型2加横向流; (b) Φ200模型5加横向流. 图中横纵坐标对应模拟空间坐标zx, 其量纲为激光波长. 横向流速等于离子声速

      Figure 1.  Beam deflection simulation results at different incident intensity: (a) Transverse flow and average intensity lower than filamentation threshold; (b) transverse flow and average intensity higher than filamentation threshold. x and yaxes of two figures corresponding to xand z axes of simulation coordinates, respectively. The spatial scale is in unit of laser wave length. The transverse flow speed equals ion sound speed.

      模型光斑平均强度/${\rm{W} } \cdot {\rm{c} }{ {\rm{m} }^{ {\rm{ - 2} } } }$成丝现象束偏折现象
      14.30 × 1013
      23.86 × 1014
      31.07 × 1015
      42.11 × 1015
      53.49 × 1015

      表 1  空间束匀滑光束Φ200成丝和束偏折现象模拟结果

      Table 1.  Simulation results for filamentaion and beam deflection in the case of CPP smoothed beam Φ200.

      对于空间束匀滑激光在均匀等离子体传播过程中能否发生成丝现象, 可根据下面公式[15]来判断,

      $ {I_{{\rm{av}}}} \leqslant {\left( {0.8{\text{π}}\frac{{{n_{\rm e}}}}{{{n_{\rm{c}}}}}\frac{{{F^2}\lambda _0^2}}{{{T_{\rm e}}({\rm{keV}})}}} \right)^{ - 1}}, $

      其中, ${I_{{\rm{av}}}}$, ${\lambda _0}$, ${n_{\rm e}}$, ${n_{\rm{c}}}$, ${T_{\rm e}}$F分别为光斑平均光强、激光波长、电子密度、等离子体临界密度、电子温度和光学F数. 当空间束匀滑入射光束平均光强大于${I_{{\rm{av}}}}$时, 空间束匀滑光束在传播过程中会发生成丝现象; 反之, 小于此值时, 空间束匀滑光束不成丝. 根据上文给出的模拟参数, 可知表1中模型3, 4, 5可发生成丝不稳定性情况, LAP3D模拟结果也证实了这一点.

      对于束偏折现象的模拟表明, 如果空间束匀滑光束发生了成丝现象, 则引入声速量级的横向流后, 光束会发生显著的束偏折现象; 反之, 如果光束没有发生成丝现象, 则不会出现束偏折现象. 如图1所示, 模型2中没有发生激光成丝现象, 加横向流时对应的模拟结果中光束传播方向没有变化; 而模型5中发生了激光成丝现象, 加横向流后对应模拟结果显示光束传播方向随传播距离增加明显偏向等离子体流方向, 即发生了束偏折现象, 这些模拟结果与物理预期相同.

      此外从空间束匀滑光束入射面和出射面内光斑电场幅值空间分布情况也可判断光束是否发生了明显的束偏折现象. 如图2所示, 对比Φ200模型5加横向流时模拟结果给出的电场幅值在入射面和出射面内的分布, 可以发现在出射面内光斑整体偏向等离子体流的方向, 表明发生了束偏折现象.

      图  2  对比空间束匀滑光束发生束偏折时Φ200模型5加横向流的入射面和出射面内光斑电场幅值分布 (a)入射面; (b)出射面. 图中横纵坐标对应模拟空间坐标xy, 其量纲为激光波长. 横向流速等于离子声速

      Figure 2.  Comparison of spatial distribution of laser electric field between laser entrance and exit planes as beam deflection presents: (a) Laser entrance plane; (b) laser exit plane. x and y axes of two figures corresponding to x and y axes of simulation coordinates, respectively. The spatial scale is in unit of laser wave length. The transverse flow speed equals ion sound speed.

      Φ200对应的模拟结果表明, 影响空间束匀滑光束明显发生束偏折的物理条件与高斯光束时的主要物理因素相一致, 即空间束匀滑光束满足成丝条件和存在流速在离子声速量级的等离子体横向流.

    • 采用时间束匀滑的典型应用为谱色散匀滑化技术, 它利用光栅、三倍频晶体, 展宽入射光束频率, 然后使光束通过相位板和聚焦透镜, 形成远场的谱色散匀滑(smoothing by spectral dispersion, SSD)光束, 其传播主要特征表现为光束在横向上的“扫描”行为[16]. SSD光束在焦平面上的光强分布为

      $\begin{split} & E(x,y)= \\ \; &{E_0}\exp \left( {{\rm{i}}3\omega t} \right)\sum\limits_n {{J_n}\left( {3\delta } \right)} \frac{{\sin \left( {n\gamma + q} \right)}}{{n\gamma + q}}\frac{{\sin p}}{p}\sum\limits_{KL}\\ &\times{\exp \left[ {{\rm{i}}\left( {n{\Omega _m}t - 2n\gamma L - 2Lq - 2Kp + {\phi _{KL}}} \right)} \right]} , \end{split} $

      式中, $3\delta $为调制振幅; ${\varOmega _m}$为调制频率; $(p, q) = $$(x, y){{k\varDelta }/({2 f}})$, $x, y$为焦平面上的空间坐标, k是三倍频后的波数, $\varDelta $为相邻相位板面元之间的距离; ${\phi _{KL}}$$\left( {K, L} \right)$通过相位板面元造成的相移; $\gamma = {{\alpha \varDelta }/2}$, 是由不同频率以不同角度入射造成的位移. 修改(9)式可得LAP3D程序中采用的入射面上光强分布的数值计算模型为[16]

      $E(x,y) = \sum\limits_n {{\rm{ETEMP}}\left( {x,y,n} \right)} {J_n}\left( {3\delta } \right)\exp \left( {{\rm{i}}n{\varOmega _m}t} \right), $

      其中

      $\begin{split} &{\rm{ETEMP}}(x,y,n) \\= \; &{E_0}{\rm{sinc}}[0.5{k_{0y}}(n \cdot {N_{\rm{c}}}F\lambda + y)] \cdot {\rm{sinc}}[0.5{k_{0x}}x] \\ & \times \sum\limits_{KL}\exp[ - {\rm{i}}L{k_{0y}}(n \cdot {N_{\rm{c}}}F\lambda + y) \\ & - {\rm{i}}K{k_{0x}}x + {\rm{i}}{\varphi _{KL}}] , \end{split}$

      这里, ${k_{0 y}} \equiv \dfrac{{2{\text{π}}}}{{F{N_y}\lambda }}$, ${k_{0 x}} \equiv \dfrac{{2{\text{π}}}}{{F{N_x}\lambda }}$; ${N_x}$, ${N_y}$分别为沿x, y方向上相位板的面元个数. 相较于(9)式, LAP3D程序计算每一时间步采用(11)式, 可大幅节约SSD光束传播过程的计算时间并提高程序并行计算能力[16].

      仿照上节研究空间束匀滑光束出现束偏折现象条件的思路, 可认为对于SSD匀滑光束, 发生束偏折的前提是发生成丝不稳定性. 因此这里首先考察SSD光束发生成丝的条件. 美国学者Kruer等[17-21]研究分析了有一定带宽的激光抑制不稳定性发展的条件, 指出当不稳定性增长率为${\gamma _0}$的波引入带宽$\Delta \omega $时, 不稳定性引发的振幅增长可被削弱为没有带宽时的${{{\gamma _0}}/{\Delta \omega }}$倍, 即

      $\big\langle {\tilde f} \big\rangle = \beta \exp \left( {\frac{{\gamma _0^2t}}{{\Delta \omega }}} \right),$

      这里, $\tilde f$对应不稳定性增长的幅值, $\beta $由初始条件确定. 此式表明当激光带宽远大于不稳定性增长率${\gamma _0}$时, 可有效抑制不稳定性增长. 由此可认为当SSD光束带宽大于成丝不稳定性增长率时, 成丝不稳定性的发展会受到抑制.

      为考察SSD带宽对其在等离子体传播过程中发生束偏折的影响, 用LAP3D模拟了三种不同调制频率的SSD光束在均匀等离子体中的传播, 这里取调制频率分别为${10^{ - 3}}{\omega _0}$, ${10^{ - 4}}{\omega _0}$${10^{ - 5}}{\omega _0}$. 模拟中取入射SSD光束光斑平均光强为$2.11 \times $$ {10^{15}}\;{{\rm{W}}/{{\rm{c}}{{\rm{m}}^2}}}$, 其余模拟参数与上节中算例相同.

      根据模拟参数, 可分析哪种调制频率会出现束偏折现象. 成丝不稳定性最大增长率${\gamma _{\max }}$

      ${\gamma _{{\rm{max}}}} = \frac{1}{8}{\left( {\frac{{{v_{{\rm{os}}}}}}{{{v_{\rm e}}}}} \right)^2}\frac{{{\omega _{{\rm{p}}\rm e}}}}{{{k_0}{c^2}}},$

      其中, ${v_{{\rm{os}}}}$, ${v_{\rm e}}$, ${\omega _{{\rm{p}}e}}$, ${k_0}$, c分别为电子在光场中的抖动速度、电子热速度、等离子体频率、光束波矢和真空光速. 根据上文模拟参数可知成丝不稳定性最大增长率${\gamma _{\max }} \approx 2.56 \times {10^{{\rm{ - }}4}}{\rm{ }}{\omega _0}$. SSD光束带宽$\Delta \omega $

      $\Delta \omega = 2 \cdot \delta \cdot {f_{{\rm{SSD}}}},$

      其中, $\delta $为调制深度, 模拟中取$\delta = 12$, ${f_{{\rm{SSD}}}}$为调制频率. 由本文算例所取的SSD光束调制频率可知, 当SSD光束调制频率${f_{{\rm{SSD}}}}$${10^{ - 3}}{\rm{ }}{\omega _0}$${10^{ - 4}}{\rm{ }}{\omega _0}$时, SSD光束带宽均远大于算例对应的成丝最大增长率${\gamma _{\max }} \approx 2.56 \times {10^{ - 4}}{\rm{ }}{\omega _0}$. 但当SSD光束调制频率${f_{{\rm{SSD}}}}$${10^{ - 5}}{\rm{ }}{\omega _0}$时, SSD光束带宽小于成丝最大增长率.

      根据上面的理论分析, 可以预期在取${10^{ - 3}}{\rm{ }}{\omega _0}$${10^{ - 4}}{\rm{ }}{\omega _0}$这两种调制频率时, 时间束匀滑光束可以抑制成丝不稳定性; 相应束偏折现象也不会出现. 当SSD光束调制频率为${10^{ - 5}}{\rm{ }}{\omega _0}$时, 可发生成丝不稳定性; 当同时存在声速量级横向流时会出现显著的束偏折现象. 下面给出的模拟结果也证实了这些判断.

      图3给出了当SSD光束调制频率为${10^{ - 3}}{\rm{ }}{\omega _0}$且没有横向流时的传播过程的模拟结果. 模拟发现, 在给定的传播距离内(约2000个激光波长), 与空间束匀滑光束成直线传播相比, SSD光束在空间传播中表现出明显的“扫描”行为, 这和已有模拟结果[16]和理论预期相符.

      图  3  调制频率为${10^{ - 3}}{\omega _0}$的时间束匀滑光束传播行为 (a)对应11000激光周期; (b)对应13750激光周期. 图中横纵坐标对应模拟空间坐标zy, 其量纲为激光波长

      Figure 3.  Propagation of SSD beam at modulation frequency of 10–3ω0: (a) Corresponding simulation result at 11000 th laser periods; (b) corresponding simulation result at 13750 th laser periods. x and y axes of two figures corresponding to z and y axes of simulation coordinates, respectively. The spatial scale is in unit of laser wave length.

      由上一节给出的模拟结果可知在给定的等离子体条件下, 本节算例中的SSD光束强度远大于成丝阈值, 但图3中显示SSD光束并没有出现束发散现象. 这表明当SSD光束取调制频率为${10^{ - 3}} {\omega _0}$时, 抑制了成丝不稳定性. 由此可以预期, 当加入等离子体横向流时, 也不会明显发生束偏折现象, 模拟结果证实了这一点(图4).

      图  4  有横向离子声速量级等离子体流时调制频率为${10^{ - 3}}{\omega _0}$的时间束匀滑光束的传播行为 (a)对应11000激光周期; (b)对应13750激光周期. 图中横纵坐标对应模拟空间坐标zy, 其量纲为激光波长. 横向流速等于离子声速

      Figure 4.  Propagation of SSD beam with transverse flow at modulation frequency of 10–3 ω0: (a) Corresponding simulation result at 11000 th laser periods; (b) corresponding simulation result at 13750th laser periods. x and y axes of two figures corresponding to z and y axes of simulation coordinates, respectively. The spatial scale is in unit of laser wave length. The transverse flow speed equals ion sound speed.

      除调制频率为${10^{ - 3}} {\omega _0}$的情况, 还模拟了SSD光束调制频率为${10^{ - 4}}{\omega _0}$时, 无横向流和有横向流时的激光传播情况. 模拟结果表明, 无横向流时算例没有出现激光成丝现象, 相应模拟流速在离子声速量级的横向流时的算例如图5所示, 也没有出现束偏折现象. 当调制频率为${10^{ - 5}} {\omega _0}$时, 模拟结果表明: 无横向流时, SSD光束会出现成丝现象; 有横向流时, 光束会出现束偏折现象. 由于其调制频率较低, 所以在本文所选的模拟空间尺度内, SSD光束传播行为类似于前一节的空间束匀滑激光, 如束发散时光场空间分布的模拟结果与图1(b)类似, 故这里略去.

      图  5  对比时间束匀滑光束在调制频率为${10^{ - 4}}{\rm{ }}{\omega _0}$时的光束传播行为 (a)等离子体横向流速为零; (b)等离子体横向流速等于离子声速. 图中横纵坐标对应模拟空间坐标yz, 其量纲为激光波长

      Figure 5.  Propagation of SSD beam at modulation frequency of 10–4ω0: (a) No transverse flow; (b) the transverse flow speed equals ion sound speed. x and y axes of two figures corresponding to y and z axes of simulation coordinates, respectively. The spatial scale is in unit of laser wave length.

      综上理论分析和模拟结果都表明, 当SSD光束带宽远大于成丝稳定性最大增长率时, 通过抑制成丝不稳定性发展, 可抑制束偏折现象出现.

    • 通过对高斯型光束发生束偏折所需条件进行引申, 提出了束匀滑光束发生束偏折的条件, 并应用数值模拟程序进行了验证和研究. 模拟证实, 对于空间束匀滑激光和时间束匀滑激光来说, 只有同时满足发生成丝不稳定性和存在离子声速量级的横向流时, 才会产生显著的束偏折现象. 这表明束匀滑光束发生束偏折的物理成因与高斯型光束类似, 即成丝不稳定性产生的密度空间结构变化——“凹坑”, 以及在等离子体横向流作用下产生横向移动. 结合初步理论分析表明, 横向流速为离子声速时束偏折现象最显著.

      研究表明, 在均匀等离子体及存在离子声速量级横向流时, 如果空间束匀滑光束入射光斑平均光强大于成丝阈值, 可引发成丝不稳定性从而发生束偏折现象; 对于时间束匀滑光束, 当其带宽远大于成丝不稳定性最大增长率时, 可以抑制成丝不稳定性从而避免束偏折现象产生, 反之则会发生束偏折现象.

参考文献 (21)

目录

    /

    返回文章
    返回