Processing math: 100%

留言板

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

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

通过理论定标律研究激光间接驱动的参数空间

杨宸 范征锋 刘杰 张存波 于承新 任国利

杨宸, 范征锋, 刘杰, 等. 通过理论定标律研究激光间接驱动的参数空间[J]. 强激光与粒子束, 2019, 31: 062001. doi: 10.11884/HPLPB201931.180339
引用本文: 杨宸, 范征锋, 刘杰, 等. 通过理论定标律研究激光间接驱动的参数空间[J]. 强激光与粒子束, 2019, 31: 062001. doi: 10.11884/HPLPB201931.180339
Yang Chen, Fan Zhengfeng, Liu Jie, et al. Parameter spaces of laser indirect-driven target by theoretic scaling laws[J]. High Power Laser and Particle Beams, 2019, 31: 062001. doi: 10.11884/HPLPB201931.180339
Citation: Yang Chen, Fan Zhengfeng, Liu Jie, et al. Parameter spaces of laser indirect-driven target by theoretic scaling laws[J]. High Power Laser and Particle Beams, 2019, 31: 062001. doi: 10.11884/HPLPB201931.180339

通过理论定标律研究激光间接驱动的参数空间

doi: 10.11884/HPLPB201931.180339
基金项目: 

国家自然科学基金项目 11475032

国家自然科学基金项目 11775030

国家自然科学基金项目 11871114

中国工程物理研究院院长基金项目 201402037

详细信息
    作者简介:

    杨宸(1993—), 男,博士研究生,从事内爆物理的研究;yangchen18@gscaep.ac.cn

    通讯作者:

    刘杰(1964—),男,博士,研究员,从事激光聚变物理的研究;liu_jie@iapcm.ac.cn

  • 中图分类号: O532.13

Parameter spaces of laser indirect-driven target by theoretic scaling laws

  • 摘要: 系统地梳理了激光间接驱动点火靶内爆压缩的物理过程,使用理论方法和一维流体力学模拟给出了靶丸内爆过程中的关键定标律公式。通过这些定标律公式获得了在给定黑腔辐射温度、飞行熵增因子、整形速度和烧蚀材料的条件下,靶丸装量——半径参数空间的点火岛区域。研究了靶丸性能参数随辐射温度、飞行熵增因子等的变化规律:当靶丸所处黑腔辐射温度升高时,内爆的稳定性将变好;设计上在靶丸装量不变的条件下,靶丸半径需要减小。当靶丸的飞行熵增因子增大时,内爆增益略微减小,内爆稳定性变好;但是点火阈值因子减小导致点火岛的区域变窄。当靶丸的整形速度增大时,点火岛的区域略微变大,内爆稳定性变化不显著;设计上在靶丸装量不变的条件下,需要增大靶丸半径,这会导致靶丸壳层形状因子变大。当改变靶丸烧蚀材料,提高质量烧蚀速率与烧蚀压时,能量增益变大且稳定性增强;设计上在靶丸装量不变的条件下,需要减小靶丸半径。
  • 惯性约束聚变(ICF)[1-4]的目标是约束氘氚聚变燃料并点火,实现可控核聚变反应并产生可观的能量。激光间接驱动[5-6]是目前最流行的方式之一。将激光注入装有氘氚燃料靶丸的金或铀黑腔中转换为X射线辐射,辐射均匀地辐照到靶丸上,烧蚀驱动靶丸内爆实现氘氚燃料的高密度压缩,并在靶丸芯部产生高温热斑,达到热核点火条件,并且在燃料飞散前充分燃烧、产生足够的聚变能量。2018年美国劳伦斯利弗莫尔国家实验室(LLNL)在国家点火装置(NIF)[7]上实现了聚变能输出大于内爆动能的实验[8]。我国2016年实现了间接驱动内爆集成实验[9],2018年在黑腔中进行了内爆实验[10]。激光聚变间接驱动目前的研究难点是控制黑腔内的激光等离子体不稳定性(LPI)、控制内爆中的流体力学不稳定性以及解决内爆混合[11-14]的问题。黑腔中靶丸加持膜[15]和充气管引起的不稳定性[16]也限制了内爆压缩达到高增益点火条件。

    直接数值模拟是点火靶物理研究与设计的重要手段,为此美国发展了LASNEX[5], HYDRA[17]等辐射流体力学模拟程序,我国也发展了系列辐射流体力学程序[18]。这些程序物理建模复杂、代码量庞大、网络分辨率高,具有较高的模拟置信度。然而,对数值模拟的过多依赖导致ICF内爆过程中的基本物理机制无法简明扼要地获得,不利于系统地分析内爆过程中影响靶丸能量增益的基本物理机制和规律。

    本文的工作是用定标方法建立简单、易懂的内爆动力学模型,分析影响内爆和能量增益的基本物理机制,在物理规律上研究设计参数对靶丸性能的影响。内爆加速阶段使用火箭模型[1, 2, 19]描述,减速阶段依靠劳森判据和定标律[20]描述,热斑形成与点火采用一维流体力学开源程序MULTI-IFE[21]辅助定标。由此获得靶丸设计所需的参数空间,研究参数空间在不同初始条件下的变化规律。在一定的靶丸性能要求下,给出了靶丸基本的设计参数范围。当改变靶丸的驱动辐射温度、飞行熵增因子、整形速度和烧蚀体材料的条件下,靶丸的装量、半径以及壳层厚度应当作相应地调整。

    图 1(a)所示,实际的间接驱动点火靶丸由同心的球型壳层组成。壳层分为外层烧蚀体和内层主燃料两部分。外层烧蚀体由聚苯乙烯塑料(CH)或金属铍(Be)等低原子序数物质制成,为了防止预热、控制流体力学不稳定性等,烧蚀体一般由不同的掺杂层构成,掺杂的元素有硅、锗、钨等。内层主燃料为冷冻氘氚(DT-ice)固体。壳层内部的空腔中充满稀薄的氘氚蒸汽(DT-gas)。壳层厚度比靶丸半径小很多,经过整形驱动压缩后,壳层变得更薄,因此靶丸壳层可以假设成一个具有初始装量M0且位置在R0处的薄壳,如图 1(b)所示。在薄壳假设的基础上,1.1至1.3部分分阶段对靶丸内爆过程进行物理建模。

    图  1  实际靶丸构型与简化模型
    Figure  1.  The real target configuration and the simplified model

    当靶丸受到辐射烧蚀时,壳层的一部分烧蚀体向外喷射,反作用力推动壳层聚心内爆。辐射对靶丸壳层产生的质量烧蚀速率和烧蚀压为˙m=aTr3pa=bTr3.5,其中˙m的单位是g·cm-2·μs,pa的单位是Mbar,辐射温度Tr的单位是heV。烧蚀定标律系数a, b由烧蚀体的材料决定:(1)文献[5]式(6-2)(6-3)给出了CH的烧蚀定标律系数;(2)文献[22]给出了Ge-doped CH、HDC和Cu-doped Be的质量烧蚀速率与辐射温度之间的实验数据图并拟合其数据给出这三种材料的系数a;文献[22]也给出了HDC和Cu-doped Be的烧蚀压与辐射温度之间的数据图,拟合可得bHDC=5;bBe=6.31,但是未给出Ge-doped CH烧蚀压数据。为了方便计算和讨论,我们假设有b/a=10。由此获得表 1中的数据。

    表  1  不同材料的烧蚀定标律系数
    Table  1.  The parameters for different ablator materials
    materials a/(g·cm-2·μs·heV-3) b/(Mbar·heV-3.5)
    CH 0.30 3.0
    Ge-doped CH 0.35 3.5
    HDC 0.50 5.0
    Cu-doped Be 0.75 7.5
    下载: 导出CSV 
    | 显示表格

    驱动壳层内爆的辐射由整形脉冲和主脉冲两部分构成。本文从主脉冲驱动的加速阶段开始研究,假设整形脉冲结束后,壳层获得初始内爆速度为v0。壳层的初始半径位置设为R0,壳层的初始装量(包括烧蚀层和DT主燃料)设为M0。随后,辐射主脉冲烧蚀靶丸壳层并驱动其加速内爆,根据质量守恒、牛顿第二定律可得[19]

    {dM(t)dt=4πR2(t)˙mM(t)dv(t)dt=4πR2(t)padR(t)dt=v(t) (1)

    由式(1)有dv(t)=-(pa/˙m)(dM(t)/M(t)),对其积分可得加速到t时刻的内爆速度v=v0+(pa/˙m)ln(1/δt)。物理上,随着壳层向内聚心,靶丸体积逐渐变小,内部气体压力逐步上升,当内部气体压力等于靶丸驱动压力时,靶丸壳层获得最大内爆速度,之后壳层开始减速。最大内爆速度时刻,记壳层剩余质量份额为δ=ms/M0,其中ms为剩余壳层质量,则最大的内爆速度为vimp=v0+(pa/˙m)ln(1/δ)。

    仍由式(1)得

    C0(1δ)+δlnδ=˙m2paM04π3R30=D0;C0=v0˙mpa+1 (2)

    用初始参量表达最大内爆速度有vimp=v0-10lnδTr0.5=v0δ-1+[10Tr0.5(δ-1-1)-(˙m/M0)(4π/3)R03δ-1]。进一步简化计算,将此加速过程近似为均加速的,使用v≈(v0+vimp)/2可得壳层的加速时间τ=R0/v=R0/(v0-5lnδTr0.5)和平均加速度g=2(R0-v0τ)/τ2=10(5lnδTr0.5-v0)lnδTr0.5/R0

    假设加速过程中密度为常数,壳层压力等于烧蚀压,加速过程飞行熵增因子为αif,可得壳层密度为ρ=[pa/(αifAdeg)]3/5=Adeg-0.6b0.6αif-0.6Tr2.1,其中Adeg=2.17×1012 (erg/cm3)/(g/cm3)5/3是氘氚固体的特征常数[2]。壳层的初始厚度为ΔR0=M0/(ρR02)=(4π)-1Adeg0.6b-0.6αif0.6Tr-2.1M0R0-2,靶丸的初始形状因子为Ar0=R0R0=(4π)1Adeg-0.6b0.6αif-0.6Tr2.1M0-1R03

    在靶丸加速过程中,烧蚀面上压力梯度和密度梯度反向,引发烧蚀瑞利—泰勒不稳定性[23](A-RTI)增长,此时界面附近的物理量可表示为f(x,y;t)=f0(x)+kfk(x)eikyeγtt,其中

    γ1=0.9k1g1+k1L1.4k1va (3)

    是A-RTI线性增长率[24-25]。不稳定性增长率由RT增长项0.9k1g/(1+k1L)和烧蚀致稳项-1.4klva组成,其中波数为kl=2π/λ=l/R=2lR0-1; R=R0/2。式(3)中L表示密度梯度标长[24],本文用辐射自由程[2]来近似,即L=ρ/(dρ/dr)=1/(κpρ),其中κp=0.43(Z3/A2)ρTr-7/2=0.43Z3A-2Adeg-0.6b0.6αif-0.6Tr-1.4表示辐射的平均光厚[2]。烧蚀致稳项中的烧蚀速度可定义为

    va=M0(1δ)4πρ(34R0)2τ=49πA0.6deg(v05lnδT0.5r)(1δ)b0.6α0.6ifT2.1rM0R40 (4)

    这里近似加速阶段壳层的表面积从4πR02收缩到0,故va中有R=3R0/4。A-RTI发展的时间就是壳层加速内爆的时间τ。则A-RTI线性增长的指数因子为

    γlτ=0.94l1+4.65lZ3A2A1.2degb1.2α1.2ifT0.7rR105lnδT0.5r(5lnδT0.5rv0)0.396lA0.6deg(1δ)b0.6α0.6ifT2.1rM0R40 (5)

    l模扰动线性增长的放大倍数为Gfl=eγlτ,记Gfl=1~100阶模放大倍数的最大值。设计上要求A-RTI的幅值不超过壳层厚度,则所容忍的最大初始扰动幅值可以近似为ηR0/Gf

    减速过程[19, 26]中,壳层对气体做功,内部气体的压力和能量升高,壳层的动能减小。当壳层动能消耗殆尽时,达到阻滞时刻,之后靶丸开始向外飞散。阻滞时刻,壳层的熵增因子记为αs(称为阻滞熵增因子),其与飞行熵增因子的关系为αs=αifM=α0.85ifv0.5impp0.1aγ0.25A0.15deg,其中M=vimp/cif=vimp/γpa/ρ是最大内爆速度时刻的马赫数[26, 27]cif是最大内爆速度时刻壳层的声速。又记阻滞时刻靶丸的半径为rs、壳层的密度和压力分别为ρsp。并且我们假设阻滞时刻:(1)壳层和热斑等压,热斑的质量可忽略;(2)壳层压力可用阻滞熵增因子与壳层密度表达;(3)用理想气体的内能描述壳层的内能。我们还假设最大内爆速度时刻壳层的动能可以全部转化为阻滞时刻壳层的内能。由这些假设可得

    {ms=4π3r3sρsp=αsAdegρ5/3s4π3r3spγ1=12msv2imp (6)

    化简后有

    {p=(γ12)5/2v5impα3/2sA3/2degρs=(γ12)3/2v3impα3/2sA3/2degrs=(4π3)1/3(γ12)1/2m1/3sv1impα1/2sA1/2deg (7)

    再将阻滞熵增αs换为飞行熵增αif,可得

    {p=(γ12)2.5γ0.375A1.275degv4.25impα1.275ifp0.15aρs=(γ12)1.5γ0.375A1.275degv2.25impα1.275ifp0.15ars=(4π3)1/3(γ12)0.5γ0.125A0.425degm1/3sv0.75impα0.425ifp0.05a (8)

    所以,我们获得了壳层内爆到阻滞时刻的物理量p, ρs, rs, αs对于最大内爆速度时刻的物理量ms, vimp, αif, pa和常数Adeg, γ的定标律。

    在理想情况下,阻滞或阻滞时刻之前,壳层内部(即靶丸中心)已形成点火热斑,引发聚变点火与燃烧。这样的燃烧从中心热斑开始,逐渐向外传播并引燃DT燃料壳层,进而引发所有的燃料燃烧。对于燃烧中的DT燃料等离子体,热斑和壳层等压且有均匀的压强p,此时约束热斑和壳层等离子体的是物质惯性。随后,这团等离子体将以稀疏波的形式向外膨胀,膨胀速度为阻滞时刻的声速cs,则惯性约束时间[4]可表示为

    τc=rscs=rsγp/ρs (9)

    用最大内爆速度时刻的物理量表示它,得到

    τc=(4π3)1/3(γ12)1γ0.625A0.425degm1/3sv1.75impα0.425ifp0.05a (10)

    我们需要进一步描述靶丸燃烧后的聚变反应产能和能量增益,这需要知道中心热斑和燃料壳层的总面密度。对于一个聚变点火并燃烧的靶丸,要求中心点火热斑的面密度(ρr)h>0.3 g/cm2、离子温度Th>4 keV且燃料壳层的面密度(ρr)s>1 g/cm2[28]。为了得到总面密度(包括热斑面密度和壳层面密度)的定标律[3],我们用MULTI-IFE程序模拟了从最大内爆速度时刻开始到靶丸最后飞散的过程中,热斑和壳层各个物理量的演化。我们在如下范围内,随机生成500组初始参数:壳层的初始质量ms=200~400 μg、密度ρs=10~30 g/cm3、温度Ts=1~31 eV、速度vimp=200~400 km/s,热斑的半径rh=50~300 μm、密度ρh=0.01~0.26 g/cm3、温度Th=0.1~1.6 keV。统计MULTI-IFE的输出结果后,得到了中间变量和总面密度ρr之间的数据,中间变量有壳层的质量ms,动能ks=msvimp2/2,烧蚀压pa,飞行熵增因子αif,热斑的体积Vh=4πrh3/3,质量mh=ρhVh,压力ph=ΓBρhTh,其中ΓB=7.66×1014 erg/(g·keV)是氘氚气体常数。拟合数据得到

    ρr=(4.5927×106)m0.476sk0.619sp0.0325aα0.369ifV0.603hm0.0975hp0.274h (11)

    图 2所示是拟合结果的准确性效果,横坐标是MULTI-IFE计算得到的总面密度,纵坐标是用拟合后的定标律给出的数值,可见点火面密度的定标律式(11)拟合的很好。再忽略式(11)中次方小的项,可以得到ρr=(8.2968×10-5)ms0.476ks0.624αif-0.342Vh-0.694pa-0.327,由于加速阶段的火箭模型无法确定最大内爆速度时刻热斑的参数,所以取Vh为500组随机Vh的平均值,最后得总面密度的定标律为

    ρr=(0.087325)m0.476sk0.624sα0.342ifp0.327a=(0.056663)(M0δ)1.100v1.248impα0.342ifp0.327a (12)
    图  2  500组随机初始参数拟合的定标面密度效果
    Figure  2.  The scaling MULTI-IFE data vs of areal density

    文献[3]从广义劳森判据推得了总面密度定标律,文献[25]从流体动力学相似律等价理论出发推导得到了壳层面密度和热斑面密度的定标律。

    由总面密度,我们可知氘氚燃料被点燃的比例即燃烧损耗[2, 25]

    Φ=ρrρr+HB (13)

    其中HB=6 g/cm2是燃烧参数,氘氚燃料反应后的能量产额

    Q=QDTms5mpΦ=(3.37184×1018erg/g)msΦ (14)

    其中QDT=17.6 MeV是一对氘氚原子发生核反应放出的能量,mp是单位原子质量。最后,靶丸内爆的能量增益为

    G=Q2.746F (14)

    其中,已将通过靶丸初始表面的辐射能F=4σTr4πR02τ=4πσTr4R03/(v0-5lnδTr0.5)耦合为间接驱动给与的激光能量Elaser=2.746F,即激光转换效率是36.42%,这个效率参考了文献[5, 30]的理论与数据。

    我们还需要判断靶丸内爆点火的裕量,即是否能顺利点火。由劳森判据[2-4],阻滞时刻壳层的压强p和约束时间τc乘积必须满足

    pτc=(4π3)1/3(γ12)1.5γ0.25A0.85degm1/3sv2.5impα0.85ifp0.1a>C (16)

    所以对于能顺利点火的靶丸,其壳层加速完成时的质量

    ms>mmin=C3(4π3)[(γ12)4.5γ0.75A2.55deg]v7.5impα2.55ifp0.3a (17)

    故最小点火能量

    Eigt=C32(4π3)[(γ12)4.5γ0.75A2.55deg]v5.5impα2.55ifp0.3a (18)

    定义出一维点火阈值因子

    ITF1D=12msv2impEigt=[C3(4π3)1(γ12)4.5γ0.75A2.55deg]msα2.55ifv7.5impp0.3a (19)

    ITF1D>1的靶丸是可以点火的靶丸。再使用ms=M0δ,并且在一个典型的点火点(ms=0.3 mg; vimp=360 km/s; αif=1.4; pa=131.6 Mbar; ITF1D=3.0)上确定常系数后有

    ITF1D=(2.9145×1057)M0δv7.5impα2.55ifp0.3a (20)

    ITF1D描述了靶丸能否顺利点火。文献[3]从广义劳森判据推得了点火阈值因子,文献[29]讨论了考虑热斑扰动和热斑纯度后的点火阈值因子。

    通过上述内爆动力学模型,我们获得了靶丸设计参量与性能参量之间的定标关系。关键的靶丸设计参量有:靶丸烧蚀体的材料系数a, b,靶丸的装量M0与初始半径R0及其关联的剩余质量份额δ与形状因子Ar0,靶丸所处环境参量即黑腔的辐射温度Tr,靶丸的内爆过程参量飞行熵增因子αif与初始内爆速度v0。关键的靶丸性能参量有:点火阈值因子ITF1D,核反应能量增益G,通过靶丸初始表面的辐射能F,靶丸的不稳定性增长因子Gf和允许的最大初始扰动η。通过分析这些设计参量与性能参量之间的变化关系,我们研究了靶丸设计参数空间的变化规律。

    靶丸的性能参量随靶丸装量——半径的变化规律是基本的关注点。在确定的靶丸烧蚀材料CH,飞行熵增因子αif=1.4,辐射温度Tr=3 heV和初始内爆速度v0=0 km/s的条件下,我们获得图 3。由图 3(a),对于确定的靶丸装量,靶丸的初始半径越大,被烧蚀的外表面越大,壳层加速完成时剩余质量份额越小。为了保留足够的燃料质量,需要δ>0.1。由图 3(b),点火阈值因子ITF1D>1才能顺利点火,考虑一定裕量,取ITF1D>2。δ>0.1且ITF1D>2的区域即为靶丸点火所考虑的范围,如图 3(c)以黄色所示,即点火岛区域。

    图  3  剩余质量比δ和点火阈值因子ITF1D在靶丸装量——半径参数空间中的变化及点火岛区域
    Figure  3.  The contour of residual mass fraction δ and ignition threshold factor ITF1D in the capsule mass-radius parameter space and the volcanic islands

    考虑点火岛中的性能参量。如图 4所示,初始装量为3 mg,初始半径为1 mm是主要参考点。该点靶丸需要消耗1.49 MJ激光能量转换为0.542 MJ通过靶丸初始表面的X射线能量,预计会产生8.46 MJ的聚变能得到5.69倍的能量增益。并且该靶丸在获得386.5 km/s的最大内爆速度时,壳层剩余质量份额为10.7%。靶丸的点火阈值因子ITF1D=5.59,形状因子是41.76。加速阶段产生的烧蚀瑞利—泰勒不稳定性增长因子为346.3,故最大初始扰动不能超过0.07 μm。

    图  4  点火岛中的参量
    Figure  4.  The parameters in the volcanic islands

    辐射温度Tr由3 heV增大到3.3 heV,质量烧蚀速率和烧蚀压增大,对装量相同的靶丸,为了保证烧蚀效果不变,必须减小靶丸受到烧蚀的面积,使得壳层所有的质量不被辐射烧蚀掉。如图 5所示,点火岛位置将向靶丸初始半径减小的方向移动,点火岛区域大小变化不显著。

    图  5  辐射温度对点火岛的影响
    Figure  5.  The effect of radiation temperature on the volcanic islands

    图 6所示,增大辐射温度后,虽然质量烧蚀速率和烧蚀压增大,但是在相同剩余质量份额的条件下,加速时间减小,使得通过靶丸初始表面的辐射能基本不变,需要的激光能量不变。靶丸产生的聚变能量与辐射温度无关,致使提高辐射温度后能量增益不变。增大辐射温度后,不稳定性的烧蚀致稳作用增强,不稳定性增长率和发展时间都变小,靶丸的稳定性变好,可以容忍更大的初始扰动。设计上,在靶丸装量不变的条件下,需要减小靶丸的半径,同时保证形状因子不变,所以靶丸的厚度也相应需要减小。

    图  6  不同辐射温度下性能参量的比较
    Figure  6.  The comparisons of parameters in different radiation temperatures

    飞行熵增因子αif由1.4增大到2.0,虽然不影响壳层的最大内爆动能,但是阻滞时刻可以获得的压力变小密度变小,最小点火能量的增大使得点火阈值因子变小,最终导致点火岛变窄。如图 7所示,点火岛的位置不变,但是范围变小。

    图  7  飞行熵增对点火区域的影响
    Figure  7.  The effect of in-flight adiabat on the volcanic islands

    图 8所示,增大飞行熵增因子后,虽然通过靶丸初始表面的辐射能不变,需要的激光能量不变,但是点火面密度略微减小,产能因此减小,所以能量增益略微减小。增大飞行熵增因子后,加速阶段烧蚀体被压缩的密度减小,其平均光厚减小,密度梯度标长增大,烧蚀速度增加,故此烧蚀致稳作用增强。不稳定性增长因子变小,允许的最大初始扰动变大,内爆稳定性变好。另外,由于减小了壳层的压缩密度,设计上需要增加靶丸的初始厚度,形状因子因此变小。

    图  8  不同飞行熵增下性能参量的比较
    Figure  8.  The comparisons of parameters in different in-flight adiabat

    整形速度v0由0 km/s增大到40 km/s,根据式(2)在靶丸装量不变的条件下,相同剩余质量份额对应的靶丸初始半径增大,又由于整形速度对最大内爆速度的贡献,点火阈值因子提高。如图 9所示,点火岛的位置向靶丸初始半径增大的方向移动,点火岛的区域范围变大。

    图  9  整形速度对点火区域的影响
    Figure  9.  The effect of shaping velocity on the volcanic islands

    图 10所示,增大初始整形速度后,虽然最大内爆速度增大,导致最大内爆动能增大。但是装量相同的靶丸初始半径增大,导致通过靶丸初始表面的辐射能增大,需要的激光能量增大。虽然内爆速度的增大对点火面密度有贡献,致使产能提高,但是权衡也增大的激光能量后,能量增益基本不变。增大初始整形速度后,内爆不稳定性增长因子的增加不超过vimp/(vimpv0)倍,对vimp=300~400 km/s, v0=40 km/s,不稳定性增长因子增加不会超过5.4%~7.4%,即内爆稳定性变化不显著。另外,由于相同装量下靶丸初始半径增大,设计上壳层初始厚度变小、形状因子变大。壳层厚度的变小使得允许的最大初始扰动变小。

    图  10  不同整形速度下性能参量的比较
    Figure  10.  The comparisons of parameters in different shaping velocity

    当改变靶丸的烧蚀材料,提高质量烧蚀速率与烧蚀压时,根据式(2)对于相同的δ参数D0是不变的,质量烧蚀速率˙m的提高使得R0减小,有R0~a-1/3b~R0-3。点火阈值因子不受影响,只是跟着δ移动位置。如图 11所示,点火岛向靶丸初始半径减小的方向移动,点火岛的区域范围大小不变。

    图  11  靶丸烧蚀材料对点火区域的影响
    Figure  11.  The effect of ablator material on the volcanic islands

    图 12所示,当改变靶丸的烧蚀材料,提高质量烧蚀速率与烧蚀压时,由于初始半径R0明显减小, 通过靶丸初始表面的辐射能明显减小,需要的激光能量明显减小。点火面密度不受影响产能不变,由于节约了驱动能量,所以能量增益变大。另外,虽然烧蚀材料烧蚀性能的改变,增大了烧蚀体的密度ρ~b0.6,使得烧蚀体的平均光厚增大κp~b0.6,密度梯度标长减小L~b-1.2,但是初始半径减小的效应更强。最终效果是增大了烧蚀速度va~b-0.6R0-4~b0.73,增加了烧蚀致稳~b1.06的效果,不稳定性增长因子明显减小,可以容忍的初始扰动增大,稳定性明显变好。设计上,对于装量相同的靶丸由于初始半径的减小,初始厚度需要增加,靶丸形状因子减小。

    图  12  不同靶丸烧蚀材料下性能参量的比较
    Figure  12.  The comparisons of parameters in different ablator material

    本文通过理论方法和一维流体力学模拟,给出了激光间接驱动点火靶丸内爆过程的定标律公式。通过这些定标律公式获得了在给定黑腔辐射温度、飞行熵增因子、整形速度和烧蚀材料的条件下,靶丸装量——半径参数空间的点火岛区域。研究了靶丸性能参数随辐射温度、飞行熵增因子等的变化规律:(1)当靶丸所处黑腔辐射温度升高时,内爆的稳定性将变好;设计上在靶丸装量不变的条件下,靶丸半径需要减小。(2)当靶丸的飞行熵增因子增大时,内爆增益略微减小,内爆稳定性变好;但是点火阈值因子减小导致点火岛的区域变窄。(3)当靶丸的整形速度增大时,点火岛的区域略微变大,内爆稳定性变化不显著;设计上在靶丸装量不变的条件下,需要增大靶丸半径,这会导致靶丸壳层形状因子变大。(4)当改变靶丸烧蚀材料,提高质量烧蚀速率与烧蚀压时,能量增益变大且稳定性增强;设计上在靶丸装量不变的条件下,需要减小靶丸半径。

    通过这些分析,我们获得了改进靶丸性能的几种方案:增大辐射温度、增大飞行熵增因子可增强A-RT不稳定性的烧蚀致稳作用,减弱内爆不稳定性;增大整形速度,可提高产能,但是并不能有效提高增益;改变靶丸的烧蚀材料提高质量烧蚀速率与烧蚀压,可以节约激光能量提高能量增益并控制不稳定性。在后续工作中,我们将进一步精细化分析影响靶丸增益和内爆不稳定性的物理机制,寻找更佳的优化设计方案。

  • 图  1  实际靶丸构型与简化模型

    Figure  1.  The real target configuration and the simplified model

    图  2  500组随机初始参数拟合的定标面密度效果

    Figure  2.  The scaling MULTI-IFE data vs of areal density

    图  3  剩余质量比δ和点火阈值因子ITF1D在靶丸装量——半径参数空间中的变化及点火岛区域

    Figure  3.  The contour of residual mass fraction δ and ignition threshold factor ITF1D in the capsule mass-radius parameter space and the volcanic islands

    图  4  点火岛中的参量

    Figure  4.  The parameters in the volcanic islands

    图  5  辐射温度对点火岛的影响

    Figure  5.  The effect of radiation temperature on the volcanic islands

    图  6  不同辐射温度下性能参量的比较

    Figure  6.  The comparisons of parameters in different radiation temperatures

    图  7  飞行熵增对点火区域的影响

    Figure  7.  The effect of in-flight adiabat on the volcanic islands

    图  8  不同飞行熵增下性能参量的比较

    Figure  8.  The comparisons of parameters in different in-flight adiabat

    图  9  整形速度对点火区域的影响

    Figure  9.  The effect of shaping velocity on the volcanic islands

    图  10  不同整形速度下性能参量的比较

    Figure  10.  The comparisons of parameters in different shaping velocity

    图  11  靶丸烧蚀材料对点火区域的影响

    Figure  11.  The effect of ablator material on the volcanic islands

    图  12  不同靶丸烧蚀材料下性能参量的比较

    Figure  12.  The comparisons of parameters in different ablator material

    表  1  不同材料的烧蚀定标律系数

    Table  1.   The parameters for different ablator materials

    materials a/(g·cm-2·μs·heV-3) b/(Mbar·heV-3.5)
    CH 0.30 3.0
    Ge-doped CH 0.35 3.5
    HDC 0.50 5.0
    Cu-doped Be 0.75 7.5
    下载: 导出CSV
  • [1] Pfalzner S. An introduction to inertial confinement fusion[M]. New York: Taylor & Francis Group LLC, 2006.
    [2] Atzeni S, Meyer-ter-Vehn J. The Physics of inertial Fusion[M]. New York: Oxford University Press, 2004.
    [3] Lindl J, Landen O, Edwards J, et al. Review of the national ignition campaign 2009-2012[J]. Phys Plasmas, 2014, 21: 020501. doi: 10.1063/1.4865400
    [4] Betti R, Hurricane O A. Inertial-confinement fusion with lasers[J]. Nat Phys, 2016, 12(5): 435-448. doi: 10.1038/nphys3736
    [5] Lindl J D, Amendt P, Berger R L, et al. The physics basis for ignition using indirect-drive targets on the National Ignition Facility[J]. Phys Plasmas, 2004, 11: 339. doi: 10.1063/1.1578638
    [6] Glenzer S H, Spears B K, Edwards M J, et al. First implosion experiments with cryogenic thermonuclear fuel on the National Ignition Facility[J]. Plasma Phys Control Fusion, 2012, 54: 045013. doi: 10.1088/0741-3335/54/4/045013
    [7] Moses E I, Boyd R N, Remington B A, et al. The National Ignition Facility: Ushering in a new age for high energy density science[J]. Phys Plasmas, 2009, 16: 041006. doi: 10.1063/1.3116505
    [8] Le Pape S, Berzak Hopkins L F, Divol L, et al. Fusion energy output greater than the kinetic energy of an imploding shell at the National Ignition Facility[J]. Phys Rev Lett, 2018, 120: 245003. doi: 10.1103/PhysRevLett.120.245003
    [9] 江少恩, 董云松, 黄天晅, 等. 激光装置首次间接驱动内爆集成实验[J]. 强激光与粒子束, 2016, 28: 080101. doi: 10.11884/HPLPB201628.160111

    Jiang Shaoen, Dong Yunsong, Huang Tianxuan, et al. Initial indirect-driven implosion integrated experiment on laser facility. High Power Laser and Particle Beams, 2016, 28: 080101 doi: 10.11884/HPLPB201628.160111
    [10] 李三伟, 杨冬, 李欣, 等. 我国激光间接驱动黑腔物理实验研究进展[J]. 中国科学: 物理学力学天文学, 2018, 48: 065202. https://www.cnki.com.cn/Article/CJFDTOTAL-JGXK201806002.htm

    Li Sanwei, Yang Dong, Li Xin, et al. Recent progress of hohlraum physics experiments in indirect driven ICF in China. Sci China Ser G—Phys Mech Astron, 2018, 48: 065202 https://www.cnki.com.cn/Article/CJFDTOTAL-JGXK201806002.htm
    [11] Regan S P, Epstein R, Hammel B A, et al. Hot-spot mix in ignition-scale implosions on the NIF[J]. Phys Plasmas, 2012, 19: 056307. doi: 10.1063/1.3694057
    [12] Ma T, Patel P K, Izumi N, et al. Onset of hydrodynamic mix in high-velocity, highly compressed inertial confinement fusion implosions[J]. Phys Rev Lett, 2013, 111: 085004. doi: 10.1103/PhysRevLett.111.085004
    [13] Remington B A, Atherton L J, Benedetti L R, et al. Hydrodynamic instabilities and mix studies on NIF: Predictions, observations, and a path forward[J]. J Phys Conf Ser, 2016, 688: 012090. doi: 10.1088/1742-6596/688/1/012090
    [14] 于承新, 范征锋, 刘杰, 等. 惯性约束聚变中内爆混合的模型构建[J]. 计算物理, 2017, 34(4): 379-386. doi: 10.3969/j.issn.1001-246X.2017.04.002

    Yu Chengxin, Fan Zhengfeng, Liu Jie, et al. Modeling of shell-mixing into central hotspot in inertial confinement fusion implosion. Chinese Journal of Computational Physics, 2017, 34(4): 379-386 doi: 10.3969/j.issn.1001-246X.2017.04.002
    [15] Martinez D A, Smalyuk V A, MacPhee A G, et al. Hydro-instability growth of perturbation seeds from alternate capsule-support strategies in indirect-drive implosions on National Ignition Facility[J]. Phys Plasmas, 2017, 24: 102707. doi: 10.1063/1.4995568
    [16] Clark D S, Kritcher A L, Yi S A, et al. Capsule physics comparison of National Ignition Facility implosion designs using plastic, high density carbon, and beryllium ablators[J]. Phys Plasmas, 2018, 25: 032703. doi: 10.1063/1.5016874
    [17] Marinak M M, Kerbel G D, Gentile N A, et al. Three-dimensional HYDRA simulations of National Ignition Facility targets[J]. Phys Plasmas, 2001, 8: 2275. doi: 10.1063/1.1356740
    [18] 宋鹏, 翟传磊, 李双贵, 等. 激光间接驱动惯性约束聚变二维总体程序—LARED集成程序[J]. 强激光与粒子束, 2015, 27: 032007. doi: 10.11884/HPLPB201527.032007

    Song Peng, Zhai Chuanlei, Li Shuanggui, et al. LARED—Integration code for numerical simulation of the whole process of the indirect-drive laser inertial confinement fusion. High Power Laser and Particle Beams, 2015, 27: 032007 doi: 10.11884/HPLPB201527.032007
    [19] Saillard Y. Acceleration and deceleration model of indirect drive ICF capsules[J]. Nucl Fusion, 2006, 46(12): 1017-1035. doi: 10.1088/0029-5515/46/12/005
    [20] Herrmann M C, Tabak M, Lindl J D. Ignition scaling laws and their application to capsule design[J]. Phys Plasmas, 2001, 8: 2296. doi: 10.1063/1.1364516
    [21] Ramis R, Meyer-Ter-Vehn J. MULTI-IFE - A one-dimensional computer code for Inertial Fusion Energy (IFE) target simulations[J]. Comput Phys Commun, 2016, 203: 226-237. doi: 10.1016/j.cpc.2016.02.014
    [22] Olson R E, Rochau G A, Landen O L, et al. X-ray ablation rates in inertial confinement fusion capsule materials[J]. Phys Plasmas, 2011, 18: 032706. doi: 10.1063/1.3566009
    [23] Lu Y, Fan Z, Lu X, et al. The analysis of harmonic generation coefficients in the ablative Rayleigh-Taylor instability[J]. Phys Plasmas, 2017, 24: 102705. doi: 10.1063/1.5007076
    [24] Betti R, Goncharov V N, McCrory R L, et al. Growth rates of the ablative Rayleigh-Taylor instability in inertial confinement fusion[J]. Phys Plasmas, 1998, 5: 1446. doi: 10.1063/1.872802
    [25] Nora R, Betti R, Anderson K S, et al. Theory of hydro-equivalent ignition for inertial fusion and its applications to OMEGA and the National Ignition Facility[J]. Phys Plasmas, 2014, 21: 056316. doi: 10.1063/1.4875331
    [26] Betti R, Anderson K, Goncharov V N, et al. Deceleration phase of inertial confinement fusion implosions[J]. Phys Plasmas, 2002, 9: 2277. doi: 10.1063/1.1459458
    [27] Kemp A, Meyer-ter-Vehn J, Atzeni S. Stagnation pressure of imploding shells and ignition energy scaling of inertial confinement fusion targets[J]. Phys Rev Lett, 2001, 86(15): 3336-3339. doi: 10.1103/PhysRevLett.86.3336
    [28] Hurricane O A, Callahan D A, Casey D T, et al. Fuel gain exceeding unity in an inertially confined fusion implosion[J]. Nature, 2014, 506: 7488.
    [29] Haan S W, Lindl J D, Callahan D A, et al. Point design targets, specifications, and requirements for the 2010 ignition campaign on the National Ignition Facility[J]. Phys Plasmas, 2011, 18: 051001. doi: 10.1063/1.3592169
  • 加载中
图(12) / 表(1)
计量
  • 文章访问数:  1444
  • HTML全文浏览量:  225
  • PDF下载量:  119
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-11-26
  • 修回日期:  2019-03-26
  • 刊出日期:  2019-07-15

目录

/

返回文章
返回