百万发平台注册地址

搜索

文章查询

x
中国物理学会期刊

留言板

ZUNJINGDEDUZHE、ZUOZHE、SHENGAOREN, GUANYUBENKANDETOUGAO、SHENGAO、BIANJIHECHUBANDERENHEWENTI, NINKEYIBENYETIANJIALIUYAN。WOMENJIANGJINKUAIJININDAFU。XIEXIENINDEZHICHI!

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

JIGUANGJUBIANHEIQIANGZHONGDENGLIZITIDERELIUYANJIU

张恩浩 蔡洪波 杜报 田建民 张文帅 康洞国 朱少平

激光聚变黑腔中等离子体的热流研究

张恩浩, 蔡洪波, 杜报, 田建民, 张文帅, 康洞国, 朱少平
PDF
HTML
导出引用
导出核心图
  • 辐射流体采用限流的局域Spitzer-Härm (S-H)电子热流近似, 在预估等离子体状态时可能与实验观察存在偏差. 利用一维(1D3V)含碰撞的粒子模拟程序, 研究了激光聚变黑腔中金等离子体的电子分布函数和电子热流. 分析表明, 在等离子体的冕区, $\alpha=Z(v_{\rm{os}}/v_{\rm{te}})^2 > 1$, 电子分布函数表现为超高斯分布(m = 3.34), 克努森数$\lambda_{\rm{e}}/L_{\rm{e}}=0.011$大于局域S-H理论的临界值$2\times10^{-3}$. 这导致了局域S-H电子热流远大于实际热流. 这种实际热流受限现象将导致辐射流体模拟给出的冕区电子温度高于神光实验测量值. 而在等离子体的高密度区域, 电子分布函数仍表现为超高斯分布(m = 2.93), 克努森数$\lambda_{\rm{e}}/L_{\rm{e}}=7.58\times10^{-4}$小于局域S-H理论的临界值, 限流的局域S-H电子热流具有一定的适用性. 但电子热流严重依赖于限流因子$f_{\rm{e}}$, 辐射流体模拟需要根据不同位置的光强和电子温度调整$f_{\rm e}$的大小.
      通信作者: 蔡洪波, cai_hongbo@iapcm.z28sp.com ; 朱少平, zhu_shaoping@iapcm.z28sp.com
    • 基金项目: 科学挑战专题(资助号: TZ2016005), 国家重点研发计划(资助号: 2016YFA0401100), 国家自然科学基金联合基金(批准号: U1730449), 国家自然科学基金(批准号: 11975055和11575030)
    [1]

    Lindl J 1995 Phy. Plasma 2 3933

    [2]

    Hinkel D E, Edwards M J, Li Y T, et al. 2013 Plasma Phys. Controlled Fusion 55 124015

    [3]

    Batani D, Baton S, Casner A, et al. 2014 Nucl. Fusion 54 054009

    [4]

    Kirkwood R K, Moody J D, Kline J, et al. 2013 Plasma Phys. Controlled Fusion 55 103001

    [5]

    Spitzer L, Härm R 1953 Phys.Rev. 89 977

    [6]

    Hurricane O A, Callahan D A, Springer P T, et al. 2019 Plasma Phys. Controlled Fusion 61 014033

    [7]

    Hall G N, Jones O S, Strozzi D J, et al. 2017 Phys. Plasmas 24 052706

    [8]

    Decker C D, Mori W B, Dawson J M, et al. 1994 Phys. Plasmas 1 4043

    [9]

    Porshnev P, Bivona S, Ferrante G 1994 Phys. Rev. E 50 3943

    [10]

    Fourkal E, Bychenkov V Y, Rozmus W, et al. 2001 Phys. Plasmas 8 550

    [11]

    Matte J P, Lamoureux M, Moller C, et al. 1988 Plasma Phys. and Controlled Fusion 30 1665

    [12]

    Liu J M, Groot J S, Matte J P, et al. 1994 Phys. Rev. Lett. 72 2717

    [13]

    Langdon A B 1980 Phys. Rev. Lett. 44 575

    [14]

    Bendib A 2017 Phys. Plasmas 24 070702

    [15]

    Weng S M, Sheng Z M, Zhang J 2009 Phys. Rev. E 80 056406

    [16]

    Mora P, Yahi H 1982 Phys. Rev. A 26 2259

    [17]

    Bibi F A, Matte J P 2002 Phys. Rev. E 66 066414

    [18]

    Holstein P A, Delettrez J, Skupsky S, Matte J P 1986 J. Appl. Phys. 60 2296

    [19]

    Epperlein E M, Short R W 1991 Phys. Fluids B 3 3092

    [20]

    Zhu S P, Gu P J 1999 Chin. Phys. Lett. 16 520

    [21]

    赵斌, 郑坚 2012 中国科学: 物理学 力学 天文学 42 1196

    Zhao B, Zheng J 2012 Sci. Sin. Phys. Mech. Astron. 42 1196

    [22]

    朱少平, 古培俊 1999 强激光与粒子束 11 687

    Zhu S P, Gu P J 1999 High Power Laser Part. Beams 11 687

    [23]

    Klimo O, Tikhonchuk V T, Ribeyre X, et al. 2011 Phys. Plasmas 18 082709

  • 百万发平台注册地址

    图 1  电子的能量分布函数(PIC统计的电子分布函数严重偏离了同一个温度对应的麦克斯韦平衡分布, 而在低能段与相应的超高斯分布非常符合, 其中$0.02\ n_{\rm c}\;(100\lambda_0)$Te = 1.53 keV, α = 1.49, m = 3.34; $0.23\;n_{\rm c}\;(300\lambda_0)$Te = 2.52 keV, α = 0.66, m = 2.93; 在高能段, 电子分布函数介于麦克斯韦分布和超高斯分布之间, 这部分电子对热流贡献最大)

    Fig. 1.  Electron distribution function (dotted lines) from PIC simulation in comparison with a Maxwell equilibrium distribution (dashed and dotted lines) and a super-Gaussian distribution (solid lines). The low-energy electron distribution from PIC simulation is well fitted to the super-Gaussian distribution. PIC simulation shows that the number of energetic electrons ($3 v_{\rm {te}} - 4 v_{\rm {te}} $百万发平台注册地址) is more than that from Maxwell equilibrium distribution, but less than that from super-Gaussian distribution. These energetic electrons play an important role in electron heat flux.

    图 2  等离子体中的电子热流${{F}}_{\rm e}$ (热流以$(1\;n_{\rm c}, 1\;{\rm{keV}})$对应的电子自由流$Q_{\rm{fs0}}=1.936\times10^{19}\;{\rm J}/({\rm{cm}}^2\cdot {\rm{s}})$为单位; 限流的S-H电子热流无法解释等离子体冕区存在的热流受限现象, 而与高密度区域的电子热流比较符合; 但电子热流严重依赖于限流因子$f_{\rm{e}}$, 需要根据不同位置的光强和电子温度调整$f_{\rm{e}}$的大小)

    Fig. 2.  Electron heat flux (black line) from PIC simulation in comparison with that from flux-limited S-H model for $f_{\rm{e}}=0.05, 0.07, 0.09$. The unit of the electron heat flux ${{F}}_{\rm{e}}$ is the electron free stream Qfs0 (1nc, 1 keV) = 1.936 × 1019 J/(cm2·s). In the high density region, the electron heat flux from the limited S-H agrees well with PIC simulation result. But in the coronal region, the electron heat flux from the flux-limited S-H is much larger than that from PIC simulation.

    图 3  电子温度的空间分布(上方三条粗线对应$5000T_0$, 下方三条细线对应$1300T_0$, 蓝色点线Theory表示无热流的光路追踪理论模型; 等离子体冕区温度梯度大, 克努森数超过临界值, 限流的局域S-H电子热流无法准确描述电子热流, 造成了RDMG对冕区电子热流的高估, 从而高估了电子温度; 高密度部分区域由于实际IBA系数低于理论值, 导致RDMG和无热流光路追踪模型估计的电子温度高于PIC结果)

    Fig. 3.  百万发平台注册地址 Electron temperature from PIC simulation (solid lines) in comparison with that from radiation hydrodynamic simulation RDMG (dashed lines) and that from the optical path tracking model without heat flux (dashed and dotted lines). In the coronal region, the limited-flux S-H predicts too large electron heat flux and results in too high electron temperature from RDMG.

    图 4  等离子体中的光强衰减图(上方三条粗线对应$5000T_0$, 下方三条细线对应$1300T_0$百万发平台注册地址, 蓝色点线Theory表示无热流的光路追踪理论模型; 由于高密度区域实际IBA系数低于理论值, 导致RDMG和无热流光路追踪模型高估了能量沉积, 激光衰减更快而无法到达临界面)

    Fig. 4.  Decay curves of laser intensity from PIC simulation (solid lines) in comparison with that from radiation hydrodynamics simulation RDMG (dashed lines) and that from the optical path tracking model without heat flux (dashed and dotted lines). In the high density region, the inverse Bremsstraws absorption coefficient is smaller than the classical expression $\kappa_{\rm{s}}=\dfrac{\nu_{\rm{ei}}}{c}\dfrac{n_{\rm{e}}/n_{\rm{c}}}{\sqrt{1-n_{\rm{e}}/n_{\rm{c}}}}$百万发平台注册地址, so the laser intensity from PIC simulation is higher than that from RDMG and the optical path tracking model.

                      <dfn id='v4cwt'><optgroup id='v4cwt'></optgroup></dfn><tfoot id='v4cwt'><bdo id='v4cwt'><div id='v4cwt'></div><i id='v4cwt'><dt id='v4cwt'></dt></i></bdo></tfoot>

                      <ul id='v4cwt'></ul>

                      • [1]

                        Lindl J 1995 Phy. Plasma 2 3933

                        [2]

                        Hinkel D E, Edwards M J, Li Y T, et al. 2013 Plasma Phys. Controlled Fusion 55 124015

                        [3]

                        Batani D, Baton S, Casner A, et al. 2014 Nucl. Fusion 54 054009

                        [4]

                        Kirkwood R K, Moody J D, Kline J, et al. 2013 Plasma Phys. Controlled Fusion 55 103001

                        [5]

                        Spitzer L, Härm R 1953 Phys.Rev. 89 977

                        [6]

                        Hurricane O A, Callahan D A, Springer P T, et al. 2019 Plasma Phys. Controlled Fusion 61 014033

                        [7]

                        Hall G N, Jones O S, Strozzi D J, et al. 2017 Phys. Plasmas 24 052706

                        [8]

                        Decker C D, Mori W B, Dawson J M, et al. 1994 Phys. Plasmas 1 4043

                        [9]

                        Porshnev P, Bivona S, Ferrante G 1994 Phys. Rev. E 50 3943

                        [10]

                        Fourkal E, Bychenkov V Y, Rozmus W, et al. 2001 Phys. Plasmas 8 550

                        [11]

                        Matte J P, Lamoureux M, Moller C, et al. 1988 Plasma Phys. and Controlled Fusion 30 1665

                        [12]

                        Liu J M, Groot J S, Matte J P, et al. 1994 Phys. Rev. Lett. 72 2717

                        [13]

                        Langdon A B 1980 Phys. Rev. Lett. 44 575

                        [14]

                        Bendib A 2017 Phys. Plasmas 24 070702

                        [15]

                        Weng S M, Sheng Z M, Zhang J 2009 Phys. Rev. E 80 056406

                        [16]

                        Mora P, Yahi H 1982 Phys. Rev. A 26 2259

                        [17]

                        Bibi F A, Matte J P 2002 Phys. Rev. E 66 066414

                        [18]

                        Holstein P A, Delettrez J, Skupsky S, Matte J P 1986 J. Appl. Phys. 60 2296

                        [19]

                        Epperlein E M, Short R W 1991 Phys. Fluids B 3 3092

                        [20]

                        Zhu S P, Gu P J 1999 Chin. Phys. Lett. 16 520

                        [21]

                        赵斌, 郑坚 2012 中国科学: 物理学 力学 天文学 42 1196

                        Zhao B, Zheng J 2012 Sci. Sin. Phys. Mech. Astron. 42 1196

                        [22]

                        朱少平, 古培俊 1999 强激光与粒子束 11 687

                        Zhu S P, Gu P J 1999 High Power Laser Part. Beams 11 687

                        [23]

                        Klimo O, Tikhonchuk V T, Ribeyre X, et al. 2011 Phys. Plasmas 18 082709

                      • [1] 翁苏明, 盛政明, 张杰. 等离子体在任意强度的直流电场中产生电流的过程. 物理学报, 2009, 58(12): 8454-8460. doi: 10.7498/aps.58.8454
                        [2] 刘宏, 朱京平, 王凯. 基于随机表面微面元理论的二向反射分布函数几何衰减因子修正. 物理学报, 2015, 64(18): 184213. doi: 10.7498/aps.64.184213
                        [3] 杨 涓, 刘文一, 朱国强, 毛根旺. 真空中微波等离子体喷流电子数密度分布规律实验研究. 物理学报, 2007, 56(1): 366-370. doi: 10.7498/aps.56.366
                        [4] 江建生, 李方华. 电子的原子散射因子解析表达式的拟合. 物理学报, 1984, 33(6): 845-849. doi: 10.7498/aps.33.845
                        [5] 孙 恺, 辛 煜, 黄晓江, 袁强华, 宁兆元. 60MHz电容耦合等离子体中电子能量分布函数特性研究. 物理学报, 2008, 57(10): 6465-6470. doi: 10.7498/aps.57.6465
                        [6] 陈建文, 傅淑芬. 放电激励KrF,ArF激光器中电子能量分布函数和输运系数的计算. 物理学报, 1981, 30(9): 1165-1173. doi: 10.7498/aps.30.1165
                        [7] 朱文浩, 朱南强, 陈跃山. 射频低压等离子体电子能量分布函数的探针诊断. 物理学报, 1989, 38(2): 236-246. doi: 10.7498/aps.38.236
                        [8] 周前红, 董志伟. 弱电离大气等离子体电子能量分布函数的理论研究. 物理学报, 2013, 62(1): 015201. doi: 10.7498/aps.62.015201
                        [9] 卿绍伟, 李梅, 李梦杰, 周芮, 王磊. 二次电子分布函数对绝缘壁面稳态鞘层特性的影响. 物理学报, 2016, 65(3): 035202. doi: 10.7498/aps.65.035202
                        [10] 范海福, 韩福森. 起始相位的受限排列. 物理学报, 1981, 30(7): 921-927. doi: 10.7498/aps.30.921
                        [11] 曹昌祺. 关于层子模型中的电磁形式因子与结构波函数的积分表示. 物理学报, 1976, 151(5): 423-432. doi: 10.7498/aps.25.423
                        [12] 谭毅, 李新阳. 光束相干合成中填充因子对远场光强分布的影响. 物理学报, 2014, 63(9): 094202. doi: 10.7498/aps.63.094202
                        [13] 吴全德. 在发射式电子光学系统中,初速度按麥克斯韋分布时像的品质. 物理学报, 1957, 34(1): 90-100. doi: 10.7498/aps.13.90
                        [14] 韩韬, 刘香莲, 李璞, 郭晓敏, 郭龑强, 王云才. 线宽增强因子对光反馈半导体激光器混沌信号生成随机数性能的影响. 物理学报, 2017, 66(12): 124203. doi: 10.7498/aps.66.124203
                        [15] 赵朋程, 廖成, 杨丹, 钟选明, 林文斌. 基于流体模型和非平衡态电子能量分布函数的高功率微波气体击穿研究. 物理学报, 2013, 62(5): 055101. doi: 10.7498/aps.62.055101
                        [16] 王俊, 王涛, 唐成双, 辛煜. 甚高频激发容性耦合氩等离子体的电子能量分布函数的演变. 物理学报, 2016, 65(5): 055203. doi: 10.7498/aps.65.055203
                        [17] 罗 莹, 王若桢, 马本堃. 介电受限对CdSe量子点中受限激子的影响. 物理学报, 1999, 48(4): 729-734. doi: 10.7498/aps.48.729
                        [18] 吴自勤, 段建中. 俄歇电子能谱中背散射修正因子的计算. 物理学报, 1984, 33(3): 419-424. doi: 10.7498/aps.33.419
                        [19] 苏娜娜, 韩庆邦, 蒋謇. 无限流体中孔隙介质圆柱周向导波的传播特性. 物理学报, 2019, 68(8): 084301. doi: 10.7498/aps.68.20182300
                        [20] 陈济舟, 王俊桥, 韩甫田. 取向对径向分布函数的影响. 物理学报, 1994, 43(7): 1105-1110. doi: 10.7498/aps.43.1105
                      • 引用本文:
                        Citation:
                        计量
                        • 文章访问数:  148
                        • PDF下载量:  6
                        • 被引次数: 0
                        出版历程
                        • 收稿日期:  2018-09-18
                        • 修回日期:  2018-11-25
                        • 刊出日期:  2019-02-01

                        激光聚变黑腔中等离子体的热流研究

                        • 1. 中国工程物理研究院研究生院, 北京 100088
                        • 2. 北京应用物理与计算数学研究所, 北京 100094
                        • 3. 北京大学应用物理与技术研究中心, 北京 100871
                        • 4. 中国工程物理研究院激光聚变研究中心, 等离子体物理重点实验室, 绵阳 621900
                        • 通信作者: 蔡洪波, cai_hongbo@iapcm.z28sp.com ; 朱少平, zhu_shaoping@iapcm.z28sp.com
                          基金项目: 科学挑战专题(资助号: TZ2016005), 国家重点研发计划(资助号: 2016YFA0401100), 国家自然科学基金联合基金(批准号: U1730449), 国家自然科学基金(批准号: 11975055和11575030)

                        摘要: 辐射流体采用限流的局域Spitzer-Härm (S-H)电子热流近似, 在预估等离子体状态时可能与实验观察存在偏差. 利用一维(1D3V)含碰撞的粒子模拟程序, 研究了激光聚变黑腔中金等离子体的电子分布函数和电子热流. 分析表明, 在等离子体的冕区, $\alpha=Z(v_{\rm{os}}/v_{\rm{te}})^2 > 1$, 电子分布函数表现为超高斯分布(m = 3.34), 克努森数$\lambda_{\rm{e}}/L_{\rm{e}}=0.011$大于局域S-H理论的临界值$2\times10^{-3}$. 这导致了局域S-H电子热流远大于实际热流. 这种实际热流受限现象将导致辐射流体模拟给出的冕区电子温度高于神光实验测量值. 而在等离子体的高密度区域, 电子分布函数仍表现为超高斯分布(m = 2.93), 克努森数$\lambda_{\rm{e}}/L_{\rm{e}}=7.58\times10^{-4}$小于局域S-H理论的临界值, 限流的局域S-H电子热流具有一定的适用性. 但电子热流严重依赖于限流因子$f_{\rm{e}}$, 辐射流体模拟需要根据不同位置的光强和电子温度调整$f_{\rm e}$的大小.

                        English Abstract

                          • 国家点火装置(NIF)未能如期实现惯性约束聚变点火[1,2], 主要原因之一就是黑腔内激光等离子体相互作用(LPI)过程过于复杂[3,4]. 这些复杂的LPI过程大大增加了点火的困难, 也导致实验结果更加难以用数值模拟进行预测. 在点火靶物理设计中, 主要采用辐射流体力学程序来开展模拟研究. 为了保障物理设计的可靠性, 辐射流体力学程序需要与聚变实验相互验证, 这需要对大多数物理过程具有非常完善的物理建模, 主要物理过程包括: 激光沉积、辐射输运、靶丸内爆、热核燃烧等. 但是受限于计算能力, 激光能量沉积往往只能采用较为粗糙的光路追踪的方法, 其中反常吸收是人为给定的份额; 等离子体中的电子热流则是采用限流的局域Spitzer-Härm (S-H)电子热流[5], 其中限流因子$ f_{\rm e} $则是人为给定的常值. 辐射流体模拟中的这些近似处理会造成对等离子体状态的预估出现偏差, 从而影响点火靶物理设计的可靠性[6,7], 因此很有必要对辐射流体程序的逆韧致吸收过程与电子热传导物理建模进行细致研究和评估.

                            在强激光的辐照下, 激光光路上的电子分布函数(EDF)表现为非麦克斯韦分布, 这将导致麦克斯韦平衡分布假设下导出的理论结果与实验数据出现明显偏差[8-10]. Matte等[11]发现当激光较强时, EDF表现为超高斯分布$ f_m({{v}})\propto {\rm e}^{-(v/v_m)^m} $, $ \alpha = Z(v_{\rm{os}}/v_{\rm{te}})^2 $越大则$ m = 2+3/(1+1.66\alpha^{-0.724}) $越大, EDF偏离平衡分布的程度越高. 其中$ v_{\rm{te}} = $$ \sqrt{(T_{\rm{e}}/m_{\rm{e}})} $为电子热速度, $ v_{\rm{os}} = eE/m_{\rm{e}}\omega_0 $为电子在激光场中的振荡速度. Liu等[12]采用等离子体探针发现电子分布表现为超高斯分布. Langdon[13]和Bendib[14]认为电子-电子碰撞项$ R_{\rm{eq}} = \dfrac{1}{2}m_{\rm{e}}v_{\rm{te}} ^2\nu_{\rm{ee}} $使电子趋于平衡分布, 而与逆韧致吸收有关的电子-离子碰撞项$ R_{\rm{heat}} = \dfrac{1}{2}m_{\rm{e}}v_{\rm{os}} ^2\nu_{\rm{ei}} $则使得电子偏离平衡分布, 当$ R_{\rm{eq}}<R_{\rm{heat}} $时, 电子分布函数由IBA主导, 能谱中电子主要集中于中间能段, 使得电子数目在中间能段变得平坦. Weng等[15]发现激光强度越大, EDF越偏离平衡分布, 并造成电子-离子碰撞频率和逆韧致吸收系数的下降, 在$ v_{\rm{os}}\gg v_{\rm{te}} $条件下EDF不满足超高斯分布, 需要在$ {{v}}_{\rm{os}} $的振荡坐标系下考虑. Mora和Yahi[16]以及Bibi和Matter[17]发现超高斯分布函数的$ 3 v_{\rm{te}}—4 v_{\rm{te}} $ 这部分电子数目明显少于麦克斯韦分布, 而这部分电子(尤其是v = $ 3.7 v_{\rm te} $) 是热流的主要携带者, 从而实际电子超高斯分布给出的热流小于电子麦克斯韦分布假设下给出的局域S-H热流.

                            在电子分布表现为超高斯特征的前提下, 研究发现, 局域S-H电子热流给出的热流${{Q}}_{\rm{SH}} = $$ -\kappa_{\rm{SH}}\nabla T_{\rm{e}} $在温度梯度较大时可能超过电子自由流$ Q_{\rm{fs}} = n_{\rm{e}}T_{\rm{e}}\sqrt{T_{\rm{e}}/m_{\rm{e}}} $, 不符合物理预期[16]. 为了解决这个难题, 辐射流体程序通常采用局域SH电子热流与电子自由流的并联限流$ {{Q}} = {{Q}}_{\rm{SH}}\cdot f_{\rm{e}}Q_{\rm{fs}}/(\mid{{Q}}_{\rm{SH}}\mid+ $$f_{\rm{e}}Q_{\rm{fs}}) $, 其中$ f_{\rm{e}} $是人为取定的限流因子, 这个热流也被称为限流的局域S-H电子热流. Holstein等[18]在对比限流的局域S-H电子热传导模型、非局域电子热传导模型和Fokker-Planck模拟时发现局域的限流S-H电子热传导模型不能准确描述低密度区域电子温度, 部分区域甚至存在热传导系数为负的现象, 而非局域热传导模型给出的温度分布更符合Fokker-Planck结果. Epperlein和Short[19]也研究发现了冕区存在热传导系数为负的现象. Mora和Yahi[16]以及Zhu和Gu[20]从局域S-H热传导理论推导出发, 研究发现当克努森数$ \epsilon = $$\lambda_{\rm{e}}/L_{\rm{e}} > 2\times10^{-3} $时, 局域S-H电子热流完全不适用, 其中, $ 2\times10^{-3} $为局域S-H理论的临界值, $ \lambda_{\rm{e}} = $ $ T_{\rm{e}} ^2/[4{\text{π}}(Z+1)n_{\rm{e}}e^4\ln\varLambda],\; L_{\rm{e}} = T_{\rm{e}}/({\rm{d}}T_{\rm{e}}/{\rm{d}}l) $分别为电子平均自由程和电子温度梯度标长. 赵斌和郑坚[21]用Fokker-Planck程序模拟了电子的非局域热输运问题, 发现等离子体冕区存在热流受限的现象, 即Fokker-Planck给出的电子热流几乎为0, 而限流的局域S-H电子热传导模型给出了很大的热流, 非局域电子热传导模型给出的热流更符合Fokker-Planck模拟结果. 在激光聚变实验中, Hall等[7]在研究NIF不同密度气体填充实验的受激拉曼散射(SRS) 散射光光谱时发现, 实验上$ 525\ {\rm nm }$散射光对应于Deplete程序的$540\ {\rm nm }$, 散射光存在红移现象, 主要原因可能是流体模拟给出的等离子体状态温度过高.

                            本文将针对黑腔等离子体的特征, 利用1D3V的粒子模拟程序(particle-in-cell, PIC) Ascent、辐射流体程序RDMG和无热流的简单光路追踪理论模型来研究强度为$ I_0 = 10^{15} \;{\rm W}/{\rm{cm}}^2$百万发平台注册地址的三倍频激光在Au等离子体中的能量沉积、电子分布函数和电子热流. 通过对比分析, 考察辐射流体程序在预估等离子体状态时的可靠性.

                          • QIANGJIGUANGZAIDENGLIZITIDECHUANBOGUOCHENGZHONG, ZHUYAOTONGGUONIRENZHIGUOCHENGJIANGNENGLIANGJIAOJIDENGLIZITIZHONGDEDIANZI, RANHOUTONGGUODIANZI-LIZIPENGZHUANGJIANGYIBUFENNENGLIANGJIAOJILIZI, CONGERZHUANHUAWEIDENGLIZITIDEDONGNENG. TONGSHI, JIGUANGHAIKEYITONGGUOJIFADENGLIZITIZHONGDEBOLIXIANGHUZUOYONGFASHENGNENGLIANGJIAOHUAN, TONGCHANGBEICHENGWEIFANCHANGXISHOU. FUSHELIUTIMOXINGDEJIBENFANGCHENGZUWEI

                            $ \left\{ {\begin{aligned}& {\frac{{\rm{d}}}{{\rm{d}}t}{{u}} = -\frac{1}{\rho}\nabla (P_{\rm{e}}+P_{\rm{i}}+P_{\rm{r}}+q)},\\ & \frac{{\rm{d}}}{{\rm{d}}t}(C_{\rm{Ve}}T_{\rm{e}}+V) = -\frac{1}{\rho}\nabla\cdot{{F}}_{\rm{e}}-P_{\rm{e}}\frac{{\rm{d}}}{{\rm{d}}t}\left(\frac{1}{\rho}\right)\\ & \qquad \qquad \qquad\qquad +W_{\rm{l}}-W_{\rm{ei}}-W_{\rm{r}},\\ & {\frac{{\rm{d}}}{{\rm{d}}t}(C_{\rm{Vi}}T_{\rm{e}}) = -\frac{1}{\rho}\nabla\cdot{{F}}_{\rm{i}}-P_{\rm{i}}\frac{{\rm{d}}}{{\rm{d}}t}\left(\frac{1}{\rho}\right)+W_{\rm{ei}}}. \end{aligned}} \right. $

                            辐射流体模型假设电子、离子各自处于麦克斯韦平衡分布, $ W_{\rm{l}} $表示激光通过等离子体的逆韧致吸收过程沉积的能量, $ W_{\rm{ei}} = W_{{\rm{ei}}0}(T_{\rm{e}}-T_{\rm{i}}) $表示电子-离子交换的能量. $ {{F}}_{\rm{e}} $表示电子热流, 采用S-H电子热流$ {{Q}}_{\rm{SH}} = -K_{\rm{e}}\nabla T_{\rm{e}} $和电子自由流$Q_{\rm{fs}} = $$ n_{\rm{e}}T_{\rm{e}}\sqrt{T_{\rm{e}}/m_{\rm{e}}} $的并联限流. 其中$ f_{\rm{e}} $百万发平台注册地址为限流因子, 是一个人为可调的常数参量.

                            $ {{F}}_{\rm{e}} = \frac{f_{\rm{e}}Q_{\rm{fs}}{{Q}}_{\rm{SH}}}{f_{\rm{e}}Q_{\rm{fs}}+|{{Q}}_{\rm{SH}}|}. $

                            为了评估电子热流对激光能量沉积和电子温度分布的影响, 本文另外给出了一个无流动无热流的简单光路追踪理论模型. 激光沿着光路通过IBA过程在等离子体中沉积能量, 等离子体中电子温度上升, 通过能量弛豫过程将一部分能量交给离子. 这里电子、离子仍假定各自处于麦克斯韦平衡分布, $ \tau_{\epsilon} ^{\rm{ei}} $为电子-离子能量弛豫时间.

                            $\left\{ {\begin{aligned} & {{\frac{{\rm{\partial}} T_{\rm{e}}}{{\rm{\partial}} t} = -(T_{\rm{e}}-T_{\rm{i}})/\tau_\epsilon ^{\rm{ei}}}},\\ & {{\frac{\partial T_{\rm{i}}}{\partial t} = Z(T_{\rm{e}}-T_{\rm{i}})/\tau_\epsilon ^{\rm{ei}}}},\\ & {\frac{3}{2}n_{\rm{e}}\frac{\partial T_{\rm{e}}}{\partial t} = \kappa_{\rm{s}} I},\\ & {\frac{{\rm{d}}I}{{\rm{d}}x} = -\kappa_{\rm{s}} I}. \end{aligned}} \right.$

                            为了避免普遍意义下的逆韧致吸收系数$ \kappa_{\rm{s}} $在临界面出现奇点, 这里仍采用WKB近似, 而通过电磁波在等离子体中的色散关系求解波数k的实部与虚部时, 不去假设$ k_{\rm{R}}\gg k_{\rm{I}} $而取$ k_{\rm{R}} = b k_{\rm{I}} $,

                            $ \left\{ {\begin{aligned} & {b = \left(\frac{n_{\rm{c}}}{n_{\rm{e}}}-1\right)\left(\frac{\omega_0}{\nu_{\rm{ei}}}\right)+\sqrt{\left(\frac{n_{\rm{c}}}{n_{\rm{e}}}-1\right)^2\left(\frac{\omega_0}{\nu_{\rm{ei}}}\right)^2+1}},\\ & {\kappa_{\rm{s}} = \sqrt{\frac{2n_{\rm{e}}\nu_{\rm{ei}}}{bn_{\rm{c}}\omega_0}}\left(\frac{\omega_0}{c}\right)}. \end{aligned}} \right. $

                            $ n_{\rm{e}}\ll n_{\rm{c}} $的低密度区域, (4)式回归普遍意义下的$ \kappa_{\rm{s}} = \dfrac{\nu_{\rm{ei}}}{c}\dfrac{n_{\rm{e}}/n_{\rm{c}}}{\sqrt{1-n_{\rm{e}}/n_{\rm{c}}}} $; 在临界面附近, (4)式与朱少平和古培俊[22]$ \kappa_{\rm{s}} $的处理一致. 考虑到临界面附近的共振吸收, 本文近似认为$ \eta = 10\% $份额的能量通过隧穿效应沉积在临界面.

                          • 在间接驱动方式激光聚变实验中, 激光从两端注入黑腔, 穿过黑腔充气区域烧蚀金壁材料并膨胀形成指数分布密度特征的高Z等离子体. 黑腔中的填充气体, 通常是惰性气体He或者辛戊烷气体, 会在X射线辐照或者激光辐照下迅速离化, 形成低密度等离子体抑制高Z等离子体的迅速膨胀[2]. 本文针对黑腔等离子体的特征, 利用含Monte-Carlo两体碰撞建模的PIC程序Ascent1D、辐射流体程序RDMG和无热流的理论光路追踪模型对比研究了强度为I0 = 1015 W/cm2的三倍频激光(λ0 = 351 nm)从左侧($ x = x_{\rm{l}} = 0 $)辐射密度为指数分布$ n_{\rm{e}} = n_{\rm{c}}{\rm e}^{(x-x_{\rm{c}})/L} $的Au等离子体的激光能量沉积和电子热流. 模拟空间尺度为(0—500λ0), 初始等离子体分布在50λ0—450λ0, 临界面位置为$ x_{\rm{c}} = 418\lambda_0 $, 等离子体标长为$ L = 80 \lambda_0 $, 模拟时间尺度为$ 5000T_0 $, 初始温度分别为Te = 1 keV, Ti = 0.25 keV. 在这个条件下, 参量不稳定过程受到抑制没有增长起来. 等离子体的电子分布函数具有良好的各向同性$ (v_{\rm{te}}\approx5 v_{\rm{os}}) $, 且临界面的等离子体也没有强耦合效应.

                            在强激光的辐照下, 等离子体电子分布函数逐渐偏离麦克斯韦分布, 表现为低能段电子数目减少. 原因是低能电子更容易从高频激光场中得到能量, 而电子-电子碰撞的能量弛豫时间为$\tau_{\epsilon} ^{\rm{ee}} = $$ 3/4\sqrt{2{\text{π}}}\cdot \sqrt{m_{\rm{e}}T_{\rm{e}} ^3}/(n_{\rm{e}}e^4\ln\varLambda) \!\approx\! 0.476 \; {\rm ps} \!= \!407 T_0\!\gg \!T_0, $ 故激光加载过程中电子-电子碰撞不足以使光路上的电子回到麦克斯韦平衡分布$\left( R_{\rm{heat}} = \dfrac{1}{2}m_{\rm{e}}v_{\rm{os}} ^2\nu_{\rm{ei}}> \right.$$\left. R_{\rm{eq}} = \dfrac{1}{2}m_{\rm{e}}v_{\rm{te}} ^2\nu_{\rm{ee}}\right)$[14]. 在过去的电子分布函数研究工作中, 通常用特征量$ \alpha = Z(v_{\rm{os}}/v_{\rm{te}})^2 $来描述电子分布函数的非平衡性[8,11]. 其中, $ v_{\rm{os}} = eE/(m_{\rm{e}}\omega_0 )$表示电子在激光场中的振荡速度, $ v_{\rm{te}} = \sqrt{{k_{\rm{B}}} {T_{\rm{e}}}/m_{\rm{e}}} $表示电子热运动的速度. 特征量$ \alpha = Z(v_{\rm{os}}/v_{\rm{te}})^2 = $$ R_{\rm{heat}}/R_{\rm{eq}} $可以很好地表征电子分布函数的非平衡性. 在强逆韧致吸收的条件下, 电子分布函数仍保持着各向同性, 并表现为超高斯分布[11]

                            $ \left\{ {\begin{aligned} & { f({{v}}) = C_m{\rm e}^{-(v/v_m)^m} },\\ & {m = 2+\frac{3}{1+1.66/\alpha^{0.724}}}, \end{aligned}} \right.$

                            $ \alpha\rightarrow0 $(无激光)时, $ m = 2 $, 电子分布函数回归麦克斯韦分布; 当$ \alpha\rightarrow\infty $(等离子体热速度趋于0) 时, $ m = 5 $百万发平台注册地址, 电子分布函数表现为极端超高斯分布.

                            图1所示, PIC模拟给出的电子分布函数与超高斯分布非常符合, 在低密度区域($ x = 100\lambda_0 $), Te= 1.53 keV, α = 1.49, m = 3.34, 在高密度区域($ x = 300\lambda_0 $), Te = 2.52 keV, α = 0.66, m = 2.93. 这与理论上低密度区域光强更大, 电子分布函数更偏离麦克斯韦分布的结果一致. 从图1可以发现, PIC统计的高能段电子数目高于超高斯分布, 而低于麦克斯韦分布. 实际EDF尾端处于$ 3 v_{\rm{te}}—4 v_{\rm{te}} $的电子数目低于麦克斯韦分布近似下的电子数目, 而这部分电子对电子热流贡献最大[16,17], 因此实际热流小于局域S-H热流$ {{Q}}_{\rm {SH}} $. 辐射流体通过简单的并联限流来减小$ {{Q}}_{\rm {SH}} $, 从而与实际电子热流更为符合.

                            图  1  电子的能量分布函数(PIC统计的电子分布函数严重偏离了同一个温度对应的麦克斯韦平衡分布, 而在低能段与相应的超高斯分布非常符合, 其中$0.02\ n_{\rm c}\;(100\lambda_0)$Te = 1.53 keV, α = 1.49, m = 3.34; $0.23\;n_{\rm c}\;(300\lambda_0)$Te = 2.52 keV, α = 0.66, m = 2.93; 在高能段, 电子分布函数介于麦克斯韦分布和超高斯分布之间, 这部分电子对热流贡献最大)

                            Figure 1.  Electron distribution function (dotted lines) from PIC simulation in comparison with a Maxwell equilibrium distribution (dashed and dotted lines) and a super-Gaussian distribution (solid lines). The low-energy electron distribution from PIC simulation is well fitted to the super-Gaussian distribution. PIC simulation shows that the number of energetic electrons ($3 v_{\rm {te}} - 4 v_{\rm {te}} $) is more than that from Maxwell equilibrium distribution, but less than that from super-Gaussian distribution. These energetic electrons play an important role in electron heat flux.

                            为了研究辐射流体模拟中限流的局域S-H电子热流与实际电子热流是否一致, 本文统计了PIC模拟各个区域的热流, 一种按热流的定义式$ {{F}}_{\rm{e}} = n_{\rm{e}}\left\langle \epsilon_{\rm{e}}{{v}}_{\rm{e}}\right\rangle $进行统计, 另一种按限流的局域S-H电子热流((2)式)进行计算. 如图2所示, 在高密度区域, 限流的局域S-H电子热流具有一定的适用性, 但热流$ {{F}}_{\rm e} $严重依赖于限流因子$ f_{\rm{e}} $. 0.23nc—0.43nc (300λ0—350λ0)区域的限流因子取0.05较为合适, 而0.43nc—0.80nc (350λ0—400λ0)区域的限流因子则取0.09较为合适. 这与理论上的限流因子$ f_{\rm e} $随着驱动激光的增强而减小, 随着热阻增强而减小相一致. 辐射流体对空间尺度较大的等离子体中的限流因子取一个固定值会引入明显的误差. 在低密度区域, PIC模拟出现了热流受限现象: 按热流定义式$ {{F}}_{\rm{e}} = n_{\rm{e}}\left\langle \epsilon_{\rm{e}}{{v}}_{\rm{e}}\right\rangle $统计的热流远小于限流的局域S-H电子热流. 这与赵斌和郑坚[21]的Fokker-Planck模拟结果一致. 在0.12nc—0.23nc (250λ0—300λ0)之间甚至出现了热传导系数为负的现象, 需要用非局域热传导模型做出解释[18,19].

                            图  2  等离子体中的电子热流${{F}}_{\rm e}$ (热流以$(1\;n_{\rm c}, 1\;{\rm{keV}})$对应的电子自由流$Q_{\rm{fs0}}=1.936\times10^{19}\;{\rm J}/({\rm{cm}}^2\cdot {\rm{s}})$为单位; 限流的S-H电子热流无法解释等离子体冕区存在的热流受限现象, 而与高密度区域的电子热流比较符合; 但电子热流严重依赖于限流因子$f_{\rm{e}}$, 需要根据不同位置的光强和电子温度调整$f_{\rm{e}}$的大小)

                            Figure 2.  Electron heat flux (black line) from PIC simulation in comparison with that from flux-limited S-H model for $f_{\rm{e}}=0.05, 0.07, 0.09$. The unit of the electron heat flux ${{F}}_{\rm{e}}$ is the electron free stream Qfs0 (1nc, 1 keV) = 1.936 × 1019 J/(cm2·s). In the high density region, the electron heat flux from the limited S-H agrees well with PIC simulation result. But in the coronal region, the electron heat flux from the flux-limited S-H is much larger than that from PIC simulation.

                            本文进一步统计了等离子体各个位置的电子温度和激光到达各个位置的相对光强, 如图3图4所示. 一方面, 在低密度区域, 光强衰减速度一致, 无热流的简单光路追踪理论模型的电子温度与PIC较为符合, 而采用限流的局域S-H电子热流的RDMG程序给出的电子温度则明显高于PIC模拟结果. 这是因为RDMG程序的限流的局域S-H电子热流加热了低密度区域等离子体, 而PIC模拟发生了热流受限现象. 另一方面, 低密度区域的温度梯度比高密度区域更大, 温度梯度标长更小; 低密度区域的密度远低于高密度区域, $ \nu_{\rm{ee}}, \nu_{\rm{ei}} $更小, 电子平均自由程$ \lambda_{\rm{e}} = v_{\rm{te}}/(\nu_{\rm{ee}}+\nu_{\rm{ei}}) $更大, 因此低密度区域克努森数$ \epsilon = \lambda_{\rm{e}}/L_{\rm{e}} $更大. 当克努森数$ \epsilon $大于S-H理论的临界值$ 2\times10^{-3} $时, 局域S-H电子热流完全不适用[16,20]. $ \epsilon\,(x \!=\! 100 \lambda_0) \!=\! 0.011$ $ > 2 \times 10^{-3}, \;\epsilon\,(x = 300 \lambda_0) = 7.58 \times 10^{-4} < 2\times10^{-3}, $ 因此局域S-H电子热流在低密度区域完全不适用, 即使采用并联限流的处理方式也不能使理论热流与实际热流相符合; 而高密度区域限流的局域S-H电子热流具有一定的适用性. 热流受限导致辐射流体高估了等离子体冕区的热流, 从而高估了冕区的电子温度, 造成SRS散射光光谱出现红移[7]. 而高密度区域RDMG和无热流光路追踪理论模型的电子温度高于PIC模拟结果是因为理论IBA系数$ \kappa_{{\rm s}} = \dfrac{\nu_{{\rm {ei}}}}{\rm c}\dfrac{n_{{\rm e}}/n_{{\rm c}}}{\sqrt{1-n_{{\rm e}}/n_{{\rm c}}}} $随密度急剧增大, 而实际IBA系数随电子-离子碰撞频率$ \nu_{{\rm {ei}}} $的变化相对平缓[23].

                            图  3  电子温度的空间分布(上方三条粗线对应$5000T_0$, 下方三条细线对应$1300T_0$百万发平台注册地址, 蓝色点线Theory表示无热流的光路追踪理论模型; 等离子体冕区温度梯度大, 克努森数超过临界值, 限流的局域S-H电子热流无法准确描述电子热流, 造成了RDMG对冕区电子热流的高估, 从而高估了电子温度; 高密度部分区域由于实际IBA系数低于理论值, 导致RDMG和无热流光路追踪模型估计的电子温度高于PIC结果)

                            Figure 3.  Electron temperature from PIC simulation (solid lines) in comparison with that from radiation hydrodynamic simulation RDMG (dashed lines) and that from the optical path tracking model without heat flux (dashed and dotted lines). In the coronal region, the limited-flux S-H predicts too large electron heat flux and results in too high electron temperature from RDMG.

                            图  4  等离子体中的光强衰减图(上方三条粗线对应$5000T_0$, 下方三条细线对应$1300T_0$百万发平台注册地址, 蓝色点线Theory表示无热流的光路追踪理论模型; 由于高密度区域实际IBA系数低于理论值, 导致RDMG和无热流光路追踪模型高估了能量沉积, 激光衰减更快而无法到达临界面)

                            Figure 4.  Decay curves of laser intensity from PIC simulation (solid lines) in comparison with that from radiation hydrodynamics simulation RDMG (dashed lines) and that from the optical path tracking model without heat flux (dashed and dotted lines). In the high density region, the inverse Bremsstraws absorption coefficient is smaller than the classical expression $\kappa_{\rm{s}}=\dfrac{\nu_{\rm{ei}}}{c}\dfrac{n_{\rm{e}}/n_{\rm{c}}}{\sqrt{1-n_{\rm{e}}/n_{\rm{c}}}}$百万发平台注册地址, so the laser intensity from PIC simulation is higher than that from RDMG and the optical path tracking model.

                          • 本文利用含Monte-Carlo两体碰撞的一维1D3V的PIC程序Ascent、辐射流体程序RDMG以及无热流的光路追踪理论模型模拟了激光聚变黑腔中激光在金等离子体中的能量沉积以及激光光路上金等离子体的电子热流. 对比发现, 等离子体冕区存在热流受限现象, 限流的局域S-H电子热流不适用的原因是温度梯度标长较小而电子平均自由程较大, 克努森数大于局域S-H理论的临界值. 在高密度区域, RDMG采用限流的局域S-H电子热流基本上是成立的, 但不同区域的限流因子应有不同, 需要根据当地光强和电子温度估计限流因子. 实际IBA系数随$ \nu_{\rm{ei}} $的变化比理论$ \kappa_{\rm{s}} = \dfrac{\nu_{\rm{ei}}}{c}\times$$ \dfrac{n_{\rm{e}}/n_{\rm{c}}}{\sqrt{1-n_{\rm{e}}/n_{\rm{c}}}} $百万发平台注册地址更为平缓, 辐射流体模拟的IBA系数在高密度区域需要乘一个小于1的修正因子. 本文可促进理解在等离子体状态研究中辐射流体模拟和实验观察之间的差别.

                        参考文献 (23)

                        目录

                          /

                          返回文章
                          返回