想象一下，利用强大的磁场在瞬间压缩燃料，使其发生核聚变，为我们提供几乎无限的清洁能源——这就是磁化靶向聚变（Magnetized Liner Inertial Fusion, MagLIF）的宏伟蓝图。在MagLIF方案中，一个脉冲功率源驱动巨大的轴向电流通过一个圆柱形的金属壳（即靶），产生环绕靶的强磁场，磁场产生的洛伦兹力（Lorentz force）会像一只无形的大手，迅速挤压靶及其内部的氘氚燃料，使其内爆、升温直至点燃聚变反应。近年来，MagLIF因其在实验中获得了较高的劳森三乘积（Lawson triple product）以及向点火和能量增益（ignition and high gain）扩展的前景而备受关注。
然而，通往聚变能源的道路上布满荆棘，其中一个关键障碍便是流体不稳定性 。在靶被加速内爆的过程中，其内外表面任何微小的粗糙度或扰动，在巨大的加速“压力”下都可能被急剧放大，导致靶材破碎、燃料混合，最终使内爆失败，聚变“哑火”。在MagLIF靶中，主要有两种不稳定性“作祟”：一种是香肠不稳定性 ，源于圆柱形扰动表面上的磁场变化；另一种则是磁瑞利-泰勒不稳定性 ，其物理本质与经典的瑞利-泰勒不稳定性相似，即当较轻的物质推动较重的物质时，在密度梯度界面处产生的不稳定现象，只不过这里的驱动力换成了磁场。这两种不稳定性一旦失控，将严重制约聚变性能和中子产额。为了提高靶设计的鲁棒性，科学家们一直在探寻抑制这些不稳定性的方法。
其中一个潜在的抑制机制便是粘度 。粘度是流体内部抵抗形变的一种属性，就像蜂蜜比水更难搅动一样。在经典的瑞利-泰勒不稳定性理论中，粘度已被证明能起到阻尼、减缓不稳定性增长的作用。那么，在MagLIF这种涉及高温、高速、强磁场的极端物理条件下，粘度能否扮演“稳定器”的角色，有效抑制香肠不稳定性和磁瑞利-泰勒不稳定性的增长呢？这正是本篇发表在《Journal of Plasma Physics》上的论文所要探究的核心问题。
为了回答上述问题，研究人员开展了一项深入的理论与数值模拟研究。他们首先扩展了Dai等人先前对粘性MagLIF靶内爆的线性增长率分析工作，推导出了一个描述包含粘性效应的不稳定性增长的解析表达式。这个表达式可以简化为一个仅由单一无量纲数——伽利略数 Ga 来参数化的形式，这个数衡量了重力和粘性力之间的相对强弱，为理解粘度如何影响不稳定性增长提供了清晰的物理图像。为了验证理论的正确性，研究者使用了辐射磁流体动力学模拟代码FLASH进行了数值模拟。通过对比，他们证实了线性增长率的理论预测与模拟结果吻合良好。基于验证后的理论模型，他们进一步计算了能够将不稳定性增长率阻尼5%所需的临界粘度阈值 ηc 。最后，他们结合MagLIF实验的实际参数范围，评估了粘度在实际内爆条件下的潜在影响。
2. 香肠不稳定性和磁瑞利-泰勒不稳定性增长率
研究人员对一个具有恒定密度ρ和动态粘度η的长圆柱形靶进行了分析。考虑了在靶外表面存在几何扰动的模型。通过求解线性化的流体力学方程组，他们得到了考虑粘度效应的不稳定性增长率γ的色散关系方程。在短波长极限下，这个复杂的方程可以简化为一个更具物理直观性的形式。这个形式揭示出，有效重力加速度geff （由香肠不稳定性分量和磁瑞利-泰勒不稳定性分量组成）是促进增长的机制，而粘度（体现为νk2 项，其中ν是运动粘度，k是波数）是阻尼机制。两者相互竞争，共同决定了最终的不稳定性增长率。
研究者进一步引入了伽利略数 Ga = geff / (4ν2 k3 ) ，将增长率表达为仅依赖于Ga 的函数。这个无量纲数成为了判断粘性力与重力谁占主导的关键指标。当Ga << 1 时，粘性力主导，不稳定性增长受到强烈抑制；当Ga >> 1 时，重力主导，不稳定性增长接近无粘（理想）情况。他们指出，先前Dai等人的渐近表达式仅在Ga << 1 的“粘性主导流”区域有效，而这在实际的MagLIF条件下往往不成立。
研究还推导出了能最大化粘性增长率的波数kmax = (geff / 32ν2 )1/3 ，以及对应的最大增长率γmax = (geff 2 / 16ν)1/3 ，并指出此时Gamax = 8 。这意味着，在Ga < 8 的区域，增加波数（即减小扰动波长）会抑制增长；而在Ga > 8 的区域，增加波数会促进增长，这与无粘情况下的趋势一致。这项分析对于理解由靶表面粗糙度或电热不稳定性（electrothermal instability）等因素引入的短波长扰动的演化至关重要。
3. 对MagLIF的启示：临界粘度与实际粘度
研究的最终落脚点，是评估粘度在真实MagLIF内爆中的实际重要性。为此，研究者计算了在不同驱动电流、靶厚度和扰动波长条件下，能将不稳定性增长率降低5%所需的临界动态粘度ηc 。
分析显示，对于典型的实验室MagLIF实验驱动电流（约数十兆安），除了极短的扰动波长区域外，这个临界粘度阈值（对于由表面粗糙度引起的扰动，约为100 g cm-1 s-1 ；对于由电热不稳定性引起的扰动，约为0.01 g cm-1 s-1 ）远高于现实靶材能达到的粘度值。高温高能量密度等离子体的无磁化粘度系数通常难以达到如此高的量级。
结论与讨论
本研究系统探讨了粘度对MagLIF内爆中香肠不稳定性和磁瑞利-泰勒不稳定性线性增长的影响，并得出了几个重要结论：
首先，研究者通过引入伽利略数Ga ，清晰地量化了重力（驱动不稳定性）与粘性力（阻尼不稳定性）之间的竞争关系，为理解和预测粘性效应提供了有力的理论工具。
其次，研究表明，在标准的MagLIF驱动电流和内爆条件下，除非扰动波长极短，否则实际等离子体粘度值远低于能显著抑制初始线性不稳定性增长的临界阈值。这意味着，在MagLIF内爆的加速阶段，粘度并不 是抑制香肠不稳定性和磁瑞利-泰勒不稳定性初始增长的主要因素。
这一结论具有重要的实际意义。它提示MagLIF的研究者和工程师，在设计靶和优化内爆性能时，不应寄希望于通过提高靶材的粘度来抑制不稳定性。相反，控制靶的表面粗糙度、优化电流波形以减缓加速过程、或利用其他物理机制（如磁场扩散、热传导效应等）可能是更有效的稳定化途径。
最后，研究者也指出，虽然粘度在初始线性阶段作用有限，但其在非线性阶段和高温、高压缩状态下的作用仍需进一步探索。在靶减速阶段，当有效重力方向反转时，内表面也可能变得不稳定，粘度的作用可能需要重新评估。这项工作为未来深入研究更复杂物理条件下的粘性效应指明了方向，并强调了在磁惯性约束聚变这一复杂物理体系中，全面、定量地评估各种物理效应的重要性。
打赏
