磁场对激光驱动Rayleigh-Taylor不稳定性影响的数值研究*

孙伟 吕冲 雷柱 仲佳勇

1) (中国原子能科学研究院核物理所,北京 102413)

2) (北京大学物理学院,HEDPS 应用物理与技术中心,核物理与技术国家重点实验室,北京 100871)

3) (北京师范大学天文系,北京 100875)

1 引言

经典的Rayleigh-Taylor 不稳定性(RTI)最初由Rayleigh[1]和Taylor[2]分别在19 世纪和20 世纪发现.RTI 主要发生在密度不均匀分布的流体界面处,当压力梯度的方向和密度梯度的方向相反时,界面扰动增长,逐渐失稳最终导致剧烈的湍流混合现象,其典型特征是气泡和尖钉结构[3].RTI是流体和等离子体中涉及多尺度、强非线性的基本物理过程.RTI 广泛存在于自然现象中,如卷状云以及地下矿床的形成等[4].在天体物理现象中,利用X 射线天文台对超新星遗迹RX J1713.7-3946的观测发现,其鞘层附近存在毫高斯(1 mG10−7T)量级强磁场,比周围星际环境磁场高3 个数量级[5].以往的研究工作表明,RTI 是导致超新星遗迹局域磁场增强的一种可能机制[6,7].此外,RTI 在超新星爆发[8]、太阳上丝状结构的形成[9]、太阳风等离子体与地球磁层等离子体之间的相互作用过程中扮演着重要角色[10].在工程领域,如惯性约束核聚变(ICF)的研究中发现,由于RTI 和速度剪切诱导的Kelvin-Helmholtz不稳定性以及冲击波驱动的Richtmyer-Meshkov不稳定性的存在,会严重影响靶丸内爆过程,造成不同燃料层间的物质混合和能量的耗散,是导致点火失败的关键因素之一[11−18].研究人员急需寻找抑制RTI 等流体不稳定性的方法.

理论方面,Chandrasekhar[19]基于理想磁流体给出了RTI 发生的条件.大量研究表明,磁场、黏度和表面张力对RTI 具有致稳作用[20].Wang 等[21−23]提出了预热烧蚀RTI 的弱非线性模型,很好地描述了气泡和尖钉的形成过程.Reckinger 等[24]和Wieland 等[25]研究了二维可压缩流体中的RTI,分析了不同初始状态下黏度对RTI 的致稳效应.李源和罗喜胜[26]基于Layzer 的势流模型,研究了理想磁流体中横向磁场对RTI 的影响,发现RTI的演化受磁压和动压的相对方向限制.近年高功率激光系统迅速发展,为研究磁流体动力学不稳定性和相关天体物理现象提供了新的平台.Rigon 等[27]在LULI2000 激光装置上开展了Ia 型超新星中的RTI 定标模拟实验,研究了年轻超新星遗迹尺度下RTI 的演化.Sauppe 等[28]利用OMEGA 激光装置和美国国家点火装置(NIF)开展圆柱型内爆靶平台实验,首次测量了减速期间RTI 的增长因子.Gao 等[29]首次利用OMEGA-EP激光装置结合质子背光成像诊断,实验测量了非线性RTI 中的Biermann 自生磁场强度分布.Khiar等[30]首次开展了激光驱动的等离子体羽流在横向20 T 磁场中的三维动力学,发现远离靶表面的膨胀等离子体羽流不受横向磁场的影响,在真空-等离子体界面可发展出磁化RTI.Sun 等[31]对黏性和垂直磁场作用下的RTI 进行了线性分析,结果发现RTI 是否被抑制取决于黏度和外加磁场强度的比值S,当S0.1时RTI 被完全抑制;当S >10时RTI 反而被增强.Barbeau 等[32]基于美国国家点火装置提出了一个实验方案,通过模拟分析了不同水平流向外加磁场对不同密度的调制层产生的RTI,结果表明强外加磁场(80 T)对低密度的材料(20 mg/cm3)具有更强的致稳效果.尽管目前已有很多关于高能量密度状态下RTI 的理论和实验研究进展,但没有从受力的角度来做详细分析,因此进一步探究磁场对RTI 的演化影响仍十分必要.

本文提出了一个实验方案,用于产生磁化RTI.使用辐射磁流体力学(FLASH)程序[33],对激光驱动聚苯乙烯(CH)调制靶产生的RTI 进行了二维(2D)数值模拟.相较于前人的工作,本文针对磁场对RTI 演化影响的物理机理进行了更为全面的考察.系统地研究了Biermann 自生磁场、不同方向和不同强度的外加磁场对RTI 演化影响的物理机理,并做了平面靶和弱激光源注入的对比模拟.直观地展示了RTI 演化区域的电子密度分布、RTI气泡、尖钉、RTI 混合区域和冲击波界面高度随时间的变化以及RTI 演化区域的涡度分布和磁场强度分布等.研究结果可为下一步在强激光装置上开展强磁环境下的RTI 实验提供理论指导,为后续研究激光驱动的湍流磁场放大效应研究提供参考.这对深入研究ICF 具有重要意义,也为天体物理学中RTI 现象的定标模拟提供了必要的依据.本文的结构如下:第2 部分给出了理论模拟模型和实验方案;第3 部分给出了模拟结果分析;第4 部分是结论.

2 理论模拟模型和实验方案

这里使用的FLASH 程序是一个开源的、模块化、并行的多物理模拟代码,具有自适应网格细化功能[33].其中的三温辐射磁流体力学解算器(3T)可以处理激光能量沉积、多组辐射扩散、平衡电子、离子和辐射温度,适合模拟激光驱动的高能量密度物理实验[17,34].靶的状态方程(EOS)和不透明度表基于IONMIX4 数据库[35].为了便于分析,这里只考虑可压缩磁流体动力学(MHD)过程,而忽略热传导和霍尔效应的影响.涉及的主要控制方程包括:

(1)式—(3)式分别为连续性方程、动量方程和能量方程,其中ρ是密度,u为速度,j是电流密度,B是磁场强度,g是重力加速度,E是总能量,ϵ是特征内能.在FLASH 中,压力p包括辐射压、电子压力和离子压力,p*表征总压力.(4)式—(6)式为麦克斯韦方程组.

在本文提出的实验方案中,将四束总能量El为1000 J、波长λl为 0.351µm、脉宽τl为 2ns、上升沿和下降沿宽度 Δτ为 0.1 ns 的梯形平顶激光束从计算域的下方朝上正入射,如图1 所示.激光焦斑半径为Rl=100µm,功率密度P约为1.6×1015W·cm−2.2D 模拟区域为x-y平面,具体尺寸设置为Lx×Ly=2400µm×1000µm,最大分辨率为 4µm×4µm .此外,模拟中采用了等离子体自由流出的边界条件.图1 给出了驱动靶由下至上各层的结构特征,包括Au 屏蔽层、CH 调制层和低密度泡沫层.Au 屏蔽层的长度×宽度为2200µm×20µm.CH 调制层长度×宽度为 2000µm×120µm,其中基底厚度为 100µm,密度为1.05 g/cm3.基底上方引入初始单模正弦扰动,波长λl=60µm,峰谷振幅A20µm(y20 sin(2πx/60)).上方的泡沫层和CH 调制层相互嵌合,密度为 0.1 g/cm3.模拟区域的其余部分填充低密度(1×10−6g/cm3)氦气.本文中的模拟共包括7 种情况,如表1 所列.1)为无磁场、调制靶;2),3)考虑Biermann 自生磁场,其中3)为平面靶情况;4),5)有外加磁场、调制靶,分别对应x方向(垂直流向)外加 104G 和105G磁场;6)y方向(平行流向)外加 105G 磁场;7)激光能量降为原来的10%,x方向外加 104G 磁场.其中情况3)和7)为对比模拟情况,情况4)—7) 关闭了Biermann 自生磁场开关,即忽略自生磁场效应.

图1 实验方案初始设置图Fig.1.Experimental scheme initial setup.

表1 不同初始条件下模拟参数比较.Table 1.Comparison of simulation parameters under different initial conditions.

3 模拟结果分析

图2 为无磁场情况下不同时刻的电子密度分布图,这里延时时间从激光源注入时刻开始计时.当激光束从Au 屏蔽层的一侧注入(入射位置为(1000µm,0)),激光穿过Au 屏蔽层烧蚀CH 产生冲击波,通过CH 向低密度泡沫传播,在冲击波不断加速靶后扰动界面过程中引发Richtmyer-Meshkov不稳定性.激光脉冲结束后,冲击波演变为爆炸波,导致系统减速,从而在界面的参考系中引发RTI.这里将RTI 的演化过程分为3 个阶段,依次是线性增长过程、非线性增长过程、湍流混合过程.在线性增长阶段,不同密度流体间相互渗透,发生涡量沉积,初始扰动随时间以指数型增长;当扰动振幅达到和初始扰动波长相当时认为已经进入非线性增长过程,轻流体进入重流体形成气泡,重流体进入轻流体形成蘑菇状尖钉结构,后期尖钉尾部出现卷吸现象产生次级不稳定性,即Kelvin-Helmholtz 不稳定性(KHI);在湍流混合阶段,扰动界面剧烈变形,尖钉将破裂,逐渐向湍流转化.图2(a)和图2(b)对应RTI 的线性增长阶段,而图2(c)和图2(d)对应RTI 的非线性增长阶段,对于湍流混合过程将在更晚时刻出现.

图2 无磁场情况下,不同时刻电子密度分布图 (a) 2.5 ns时刻;(b) 5ns时刻;(c) 7.5 ns时刻;(d) 10 ns 时刻Fig.2.In the absence of a magnetic field,the electron density distribution at different times:(a) 2.5 ns;(b) 5ns;(c) 7.5 ns;(d) 10 ns .

图310 ns 时刻的电子密度分布图 (a)无磁场情况;(b)考虑Biermann 自生磁场情况;(c)考虑Biermann 自生磁场、平面靶情况Fig.3.Snapshots of the electron density distribution at the time of 10 ns :(a) Without magnetic field;(b) considering Biermann self-generated magnetic field;(c) Biermann self-generated magnetic field and plane target.

图3 所示为10 ns 时刻的电子密度分布图,这里延时时间从激光源注入时刻开始计时,图3(a)—(c)对应无磁场调制靶情况,仅考虑Biermann 自生磁场、调制靶情况和仅考虑Biermann 自生磁场、平面靶情况.其中Biermann 自生磁场是指在激光与固体靶相互作用中,由于等离子体的温度梯度和密度梯度方向非共线所产生的磁场,具体表示为,其中kB为玻尔兹曼常数,在二维模拟中表现为垂直于纸面的环形分量.图3(a)和图3(b)所示为典型的非线性阶段的RTI,当无磁场作用时,RTI 可充分地自由发展,其中RTI 的界面高度为750 µm,平均运动速度uint为 75 km/s .等离子体声速(Γ为绝热指数,P为压强,ρ为密度)约为64 km/s,最外侧等离子体出流界面高度为800 µm,其平均速度uk为 80 km/s,表明等离子体出流外边界处的弓形结构为冲击波结构,且RTI 在以超音速生长.对比图3(a)和图3(b)发现,两种情况下RTI 的尖钉结构几乎一致,表明Biermann 自生磁场对RTI 的影响可以忽略.此外发现,由于激光强度在空间分布不均匀,使得驱动后产生的冲击波亦存在空间上的不均匀性,后面导致在RTI 演化过程中,偏离激光中心注入位置的尖钉-气泡结构呈不对称排列,向两侧倾斜.图3(c)为对比模拟,当CH 烧蚀层不施加初始扰动时(基底厚度为120 µm),在 10 ns 时刻,并未出现RTI,表明初始扰动对RTI 的产生至关重要.

图4 为考虑不同外加磁场情况下,t10 ns 时刻RTI 演化区域的密度分布图,图4(a)—(d)分别对 应x方向外加 104G、x方向外加 105G磁场、y方向外加 105G 磁场、以及激光能量降为原来的10%时x方向外加 105G 磁场情况,这里涉及的4 种情况均忽略了Biermann 磁场.通过直观的对比可以发现,y方向外加磁场对RTI 和KHI的发展无影响,随着x方向外加磁场强度的增加,位于尖钉尾部的KHI 逐渐被抑制,而RTI 的发展不受影响.如图4(d)所示,当激光能量降低为原来的10%时,激光强度相应地下降为 1.6×1014W·cm−2,此时RTI 的界面高度明显降低(P430µm),界面速度低于声速,尖钉尾部附近的KHI 完全被抑制.这表明不稳定性的发展主要依赖于入射激光的强度,如果激光驱动产生的冲击波足够强,RTI 在强外加磁场环境下仍可以产生,反之亦然,这一判断和Sano 等[36]的相关结论也相符.

图4 不同初始外加磁场条件下,10ns时刻的电子密度分布图(a)x方向外加104 G磁场;(b) x 方向外加 105 G 磁场;(c) y 方向外加105 G磁场;(d)激光能量降为原来的10%时,x 方向外加105G磁场情况Fig.4.Snapshots of electron density distribution at 10 nsunder different initial applied magnetic field conditions:(a) With 104 G in the x-direction;(b) with 105 Gin the x-direction;(c) with 105 G in the y-direction;(d) when the laser energy is reduced to 10%,with 105 G in the x-direction.

对于RTI 的演化,混合区域(mixing zone,MZ)高度和界面高度是两个重要参数.图5(a)和图5(b)分别给出了6 种情况下RTI MZ 高度和界面高度P随时间的变化.其中黑线对应无磁场情况,红线对应只考虑Biermann 自生磁场情况,蓝线对应初始时刻x方向外加 104G 磁场,紫线对应初始时刻x方向外加 105G 磁场,绿线对应对应初始时刻y方向外加 105G 磁场,橙线对应激光能量下降为原来的 10%条件下x方向初始外加 105G 磁场情况.图5(a)中MZ 高度的变化分别对应前面描述的前两个演化阶段,0—2 ns 内,激光驱动的冲击波刚刚产生并向靶后传输,CH 等离子体没有渗透到低密度泡沫材料 中,RTI尚未形成;2—4 ns 对应RTI 线性增长阶段;4—10 ns 对应RTI 非线性增长过程,此时形成典型的泡沫-尖钉结构.图5 结果表明,当激光条件相同时,RTI 的界面动力学相同,但MZ 生长情况不同,受到x方向外加磁场的影响,MZ 高度明显低于其他情况,反映了磁场对RTI 的致稳作用;当激光能量下降时,RTI 的界面高度P和MZ 高度大幅降低.

图5 不同初始条件下,(a) RTI 的MZ 高度和(b)界面高度P 随时间的发展情况Fig.5.Development of (a) the RTI MZ height and (b) the RTI interface height with time under different initial conditions.

图6 所示为排除平面靶其余6 种情况下RTI尖钉高度和RTI 气泡高度随时间的变化情况.其中RTI 气泡为RTI 混合区中位于顶部的蘑菇状结构,RTI 尖钉则对应气泡包裹下的尖锐部分.通过对比图6(a)和图6(b)可发现,RTI 气泡和尖钉具有相近的增长规律.同样地,延时 2—4 ns 对应RTI的线性增长阶段,4—10 ns 对应RTI 的非线性增长阶段.对于无磁场情况,10 ns 内RTI 气泡的平均增长速率约 7km/s,而RTI 尖钉平均增长速率约9 km/s.随着x方向初始外加磁场强度的增加,对RTI 气泡和尖钉的抑制效果越加显著,如x方向外加 105G 磁场时,RTI 气泡的平均增长速率仅为4.2 km/s .图7 所示为排除平面靶其余6 种情况下冲击波界面高度和KHI 涡旋高度H随时间的变化.当激光条件相同时,冲击波的界面动力学相同,KHI 涡旋高度存在差异,随着x方向外加磁场强度的增加,KHI 涡旋高度逐渐下降,尖钉附近的小尺度结构被抑制.这表明x方向的外加磁场对RTI MZ 区域和尖钉尾部的KHI 有致稳作用,且KHI 涡旋更容易被抑制,从定量的角度佐证了之前的论述.

图6 不同初始条件下,(a) RTI 尖钉高度和(b) RTI 气泡高度随时间的变化情况Fig.6.Changes of (a) RTI spike height and (b) RTI bubble height with time under different initial conditions.

图7 不同初始条件下,(a)外围冲击波界面高度和(b) KHI 涡旋高度H 随时间的发展情况Fig.7.Development of (a) the peripheral shock wave interface height and (b) KHI vortex height H with time under different initial conditions.

由于在RTI 演化过程中总是伴随KHI 涡旋的形成,为了更好地理解RTI 和KHI 涡旋在不同磁场环境下的演化,图8 给出了 10 ns 时刻不同初始条件下的涡度场分布,分别为无磁场、x方向初始外加 105G 磁场、激光能量降为原来的10%条件下x方向初始外加 105G 磁场情况.涡度场可以反映流体之间的混合程度,定义为ω∇×u.在二维模拟x-y平面中仅有ωz分量,即ωz单位为 ns−1,和RTI 以及KHI 的线性增长率正相关.根据MHD 运动方程即(2)式可知,无磁场时,界面两侧流体之间的混合作用仅受负压力梯度力的影响.图8(a)中KHI 涡旋呈周期性对称排布,以黑线框中的KHI 涡旋为例,发现左右两侧KHI涡旋的涡度相反,左侧为正涡度即逆时针旋转,右侧为负涡度即顺时针旋转,由内向外涡度逐渐减小,极值分布在涡旋中心,表现为猫眼状,这和RTI尖钉演化彼此对应.当x方向加外 105G 磁场时,以图8(b)中红线框标注的RTI 尖钉为例,结果表明x方向外加磁场可以显著减弱流体之间的混合过程,此时KHI 涡旋的涡度极值分布在流体界面处,且涡度分布变得复杂.推断这是由于外加磁场和外侧流体相互压缩,破环了KHI 涡旋原有的演化进程,致使KHI 涡旋发生断裂,形成涡流片结构.此外通过对比图8(b)和图8(c)发现,二者的涡度场分布情况类似,表明降低激光能量可以减缓RTI 的演化进程.可以推断,如果需要完全抑制RTI,需要更强的水平流向外加磁场.

图8 有 无磁场 时,10 ns时刻RTI 演化区域涡度分布情况 (a)无磁场;(b)在x 方向施加 105 G 磁场;(c)当相应激光的入射能量降低到10%时,在x 方向外加 105 G 磁场Fig.8.Vorticity distribution in the RTI evolution region at 10 ns with or without magnetic field:(a) Without the magnetic field;(b) a 105 G magnetic field is applied in the xdirection;(c) when the incident energy of the corresponding laser is reduced to 10%,a 105 G is applied in the xdirection.

通过以上分析,可以直观地发现x方向的外加磁场对RTI 和尖钉尾部的KHI 有致稳作用,且KHI 涡旋更容易被抑制.为了更详细地探讨磁场的作用机制,首先分析了磁场强度分布.图9 分别给出了5 和 10 ns 两个时刻的磁场强度分布情况.图9(a)和图9(b)对应激光入射能量为 1000 J 时,x方向初始外加 105G磁场情况,在 10 ns 时刻,磁场强度被放大20 倍,约为 2×106G,这主要源于RTI 尖钉界面附近等离子体流的动能部分转化为磁场能量,从而导致磁场的有效放大,如图9(b)和图9(d) 所示,蘑菇状尖钉界面处的磁场剪切最为剧烈.对于图9(b),在t=10 ns 时,电子拉莫尔半径为rLeume/(eB)2.18µm,相应的离子拉莫尔半径rLi11.38µm,这里以10 ns 时RTI 界面高度作为系统尺度(L75µm),电子拉莫尔半径与系统尺寸的比值rLe/L约为0.03,离子拉莫尔半径与系统尺寸的比值rLi/L约为0.15.因此RTI 和KHI 演化区域的等离子体被完全磁化,表明磁场已穿透靶后的等离子体流区域,其中外部磁场被不断拉伸、挤压,最终被放大.另一方面,不稳定性界面将受到由磁场产生的洛伦兹力作用,并贯穿整个RTI 和KHI 的演化过程.

图9 x 方向初始外加 105 G磁场时,不同激光条件、不同时刻磁场强度分布图 (a),(b)分别对应入射激光能量为 1000 J时,5和 10 ns时刻磁场强度分布图;(c),(d)分别对应入射激光 100 J时 5和 10 ns 时刻磁场强度分布图Fig.9.Magnetic field intensity distribution at different times under different laser conditions,when 105 G is initially applied in xdirection:(a),(b) The magnetic field intensity distribution at 5and 10 nscorresponding to incident laser 1000 J,respectively;(c),(d) the magnetic field intensity distribution at 5and 10 nscorresponding to incident laser 100 J,respectively.

外加磁场的引入给RTI 以及KHI 的演化带来了两个后果,即外加种子磁场的放大和对不稳定性的抑制作用.其中的致稳作用主要来自于洛伦兹力的两个分量,如(9)式所示,等式右侧依次对应磁压力梯度和磁张力梯度,主要表现为和负压力梯度力之间的抗衡作用.图10 分别给出了10 ns 时刻的磁压力梯度与磁张力梯度分布情况,展示的是矢量值.其中图10(a)和图10(c)对应激光入射能量为 1000 J时,x方向初始外加 105G 磁场条件下,磁压力梯度与磁张力梯度的分布情况;图10(b)和图10(d)对应激光能量降为原来的1/10 时x方向外加 105G 磁场条件下,磁压力梯度与磁张力梯度分布情况.结合上述分析可知,两种情况下磁压力梯度及磁张力梯度分布规律类似.以图10(a)和图10(b)为例,其中磁压集中在整个扰动界面上,磁能梯度越大,磁压力梯度越大,极值分布在蘑菇尖钉附近,可达 1×1015dyn/cm3(1 dyn=10–5N).磁张力主要分布在弯曲界面处,呈现不均匀分布状态,界面曲率越大,磁张力梯度越大,最大值约4×1014dyn/cm3,其中磁压力梯度明显大于磁张力梯度,在磁压力和磁张力的共同作用下对KHI和RTI 起到致稳作用.图11 给出了10 ns 时刻的负压力梯度力的分布情况,可以发现在RTI 尖钉附近,负压力梯度力∇p极值在 2×1015dyn/cm3量级,洛伦兹力不足以抗衡负压力梯度力,使得RTI未完全被抑制,如需完全抑制其发展,初始时刻需要引入更强的水平背景磁场.如图10(a)和图10(b)中黑线框标注的RTI 尖钉尾部的KHI 涡旋附近,磁压力梯度在 1×1014dyn/cm3量级,和压力梯度力相当,使得KHI 涡旋等小尺度结构被完全抑制,其中磁压力在抑制KHI 中起主导作用:

图10 x 方向初始外加 105 G磁场时,不同激光条件下,10 ns 时刻磁压力与磁张力分布图 (a),(b)分别对应入射激光能量为1000和 100 J时磁压力分布情况;(c),(d)分别对应入射激光能量为 1000和 100 J 时磁张力分布情况Fig.10.Distribution of magnetic pressure and magnetic tension at 10 nsunder different laser conditions when 105 G is initially applied in x-direction:(a),(b) The magnetic pressure distribution corresponding to the incident laser energy of 1000and 100 J respectively;(c),(d) the magnetic tension distribution corresponding to the incident laser energy of 1000and 100 J respectively.

图11 x 方向初始外加 105 G磁场时,不同激光条件下,10 ns时刻负压力梯度力分布 (a)入射激光能量为 1000 J ;(b)入射激光能量为100 JFig.11.When a 105 Gis initially applied in the x-direction,the negative pressure gradient force distribution at 10 ns under different laser conditions:(a) when the incident laser energy is 1000 J;(b) when the incident laser energy is 100 J .

4 结论

本文提出了一种高能量密度物理中产生RTI的实验设计方案.主要用FLASH 程序对激光驱动的CH 调制靶产生的RTI 进行二维数值模拟.系统地研究了Biermann 自生磁场、外加磁场对RTI的演化影响的物理机理.模拟结果表明,Biermann自生磁场和平行流向即y方向的外加磁场在RTI演化过程中基本不会改变其界面动力学.x方向外加 105G 磁场情况下,在RTI 演化区域附近的等离子体完全被磁化.外加磁场和冲击波后的流体之间发生显著的相互作用,在0—10 ns 期间,被放大20 倍(约 2×106G),放大过程中的磁场对RTI MZ区域演化和RTI 尖钉尾部的KHI 具有明显的抑制作用.抑制效果主要来自于磁压力、磁张力和负压力梯度力之间的平衡,其中磁压力占主导作用.抑制程度和x方向初始外加磁场强度呈正相关.另外,不稳定性的发展严重依赖于入射激光的强度,当激光驱动产生的冲击波足够强,RTI 可以在相对更强的外加磁场环境下产生.后续我们将尝试开展相关的实验研究,一方面给出磁场对RTI 以及KHI演化影响的实验证据,另一方面也可以重点分析不稳定性发展过程中的磁场放大现象,磁场起源及其放大机制一直是天体物理学、宇宙学和高能量密度物理等领域的重要挑战.研究结果将为后续开展和ICF 相关的靶物理研究以及流体混合过程提供借鉴,有助于加深对星际湍流和磁场演化的理解.