丝阵 Z-Pinch 消融及内爆动力学过程的二维数值模拟研究

1. 丝阵Z-Pinch消融及内爆动力学过程的二维数值模拟研究丝阵Z-Pinch消融及内爆动力学过程的二维数值模拟研究 黄俊 丁宁 宁成 孙顺凯 北京应用物理与计算数学研究所 磁驱动聚变研究中心 2012. 4. 19

2. 报告内容 研究背景 二维(r, q)丝阵内爆动力学模型 模拟结果分析 小结与讨论

3. 研究背景 Z-pinch ：当强电流通过轻质量的气柱、箔壳、单丝或丝阵时，产生等离子体，其在角向磁场的作用下发生径向内爆压缩，最终在轴线处碰撞滞止，释放出大量X射线的过程。 应用前景：高效强X光辐射源 惯性约束聚变、高能密度物理、实验室天体物理等领域。

4. 研究背景 • 实验表明，丝阵Z-pinch过程可分为以下几个阶段: I.金属丝相变、电离，形成“芯-晕”结构及先驱等离子体； II.主体等离子体内爆； III.内爆到心、滞止、辐射。

5. 研究背景 • 丝间隙的减小大大提高了X光功率 丝阵的角向特征（丝间隙、丝等离子体直径等）如何影响Z-Pinch过程？ 需要开展(r, q)平面内的二维数值模拟研究 [ T. W. L. Sanford et. al，Phys. Rev. Lett. 77, 5063 (1996)]

6. 研究背景 • 国内外研究现状： • 英国Imperial College：Gorgon程序*( RMHD欧拉+双温近似+辐射损失) 丝间隙的减小降低了先驱等离子体份额及丝阵的角向不均匀性，内爆呈现出壳层特性。 *J. P. Chittenden et al., Phys. Rev. Lett. 83, 100, 1999

7. 研究背景 • 美国Sandia实验室：ALEGRA-HEDP程序*(3D、RMHD+ALE+辐射输运) • 丝阵消融-内爆动力学过程及其与 X光辐射功率的关系 • 拖尾质量形成机制 * C. J. Garasi et al., Phys. Plasmas 11, 2729(2004)

8. 研究背景 • 国内：刚刚起步 • 北京应用物理与计算数学研究所Z箍缩组: 采用二维理想MHD模型，研究了丝阵Z-Pinch早期消融阶段的等离子体动力学过程，给出了合理的电流分布及磁场位形*。 随着中物院10MA水平PTS装置的建成和即将投入使用，必须开展丝阵Z－Pinch内爆过程的二维(r, q)数值模拟研究，以对未来的实验进行优化设计。 * Jun Huanget al., Phys. Plasmas 18, 042704(2011)

9. 二维(r,q)丝阵内爆动力学模型 • 控制方程 • 质量守恒方程： • 动量守恒： • 能量守恒： • 磁感应方程： • 状态方程： • 坐标系为二维极坐标(r, q)，矢量有三个分量。 • 电阻率模型：h=h0(1+kr0/r) * 其中h0：Spitzer电阻率； h0k：“真空区”电阻率， k取100； r0： “真空区”密度，取1.0×10-6。电离度的求解采用Thomas-Fermi模型。 * M. R. Martin et al., Phys. Plasmas 17, 052706(2010)

10. B2i,j B1i-1,j B1i,j ε3i,j hJ3i,j B2i,j-1 二维(r,q)丝阵内爆动力学模型 • 流体力学方程和磁场对流方程：Zeus2D*程序求解。 • 磁扩散方程的求解：(参考Constrained Transport方法） Faraday定律 积分形式 将hJ3定义在网格角点位置(同e3)，则上式有差分形式（仅对扩散项） ： • 边界的hJ3在每个方程中只出现一次且方向相反，在对某一封闭曲面磁通量求和时可以自动消去，磁场的无源性条件自动得到满足； * J. M. Stone & M. L. Norman, ApJS 80, 791 (1992).

11. 二维(r,q)丝阵内爆动力学模型 • 磁扩散程序模块的编制及检验 人为构造初始磁场分布，比较模拟和解析结果 • 对两个方向的磁扩散进行检验，模拟计算结果与解析解符合的很好。

12. 二维(r,q)丝阵内爆动力学模型 • 模拟对象： • Saturn装置上的单层Al丝阵 • 电流波形 峰值：约7MA， 上升时间(10%～90％)：35ns • 模拟区域 0.2mm≤r<10.25mm, 0 ≤q< q0

13. 二维(r,q)丝阵内爆动力学模型 • 边界条件 • r=0.2mm, 对称轴 • r=R0=10.25mm，固壁 • q=0 and q0，周期边界

14. 模拟结果分析 • 丝阵等离子体内爆动力学 以16根铝丝阵为例，假定初始等离子体直径d0=0.4mm, 温度T0=104K • 内爆过程： • 丝等离子体膨胀，径向位置几乎不变； • 丝周围的低密度等离子体首先向轴运动，形成先驱等离子体； • 主等离子体内爆； • 内爆到心，反弹。

15. 模拟结果分析 • 早期动力学特征：t=10ns 白色箭头代表速度矢量 • 早期热压远大于磁压，等离子体向周围缓慢扩散； • 边界出现密度梯度。

16. 模拟结果分析 • 消融阶段：以t=65ns为例 • 丝芯周围低密度等离子体在洛伦兹力的作用下率先向轴运动，部分磁场通过对流进入丝阵内部。 • 电流密度主要集中在丝的外侧低密度区，先驱等离子体中携带少量电流； • 磁力线在高密度区呈向外凸的形状 • Gorgon程序计算得到的等离子体分布及磁力线的形状[J. P. Chittenden etal., Phys. Rev. Lett., 101, 055005, 2008]

17. “强光一号”实验使用磁探针对r=0.5R0内部总电流的测量结果， • 模拟得到的电流增长趋势与实验基本一致！ 模拟结果分析 • 丝阵内部电流分布随时间的变化 • 某径向位置r内部总电流增长过程： 1. 阶段一，主体等离子体尚未到达r，先驱携带电流进入该位置，其内部总电流快速上升； 2. 阶段二，主体等离子体经过测点，其内部总电流继续上升，且增长速度加快。

18. 模拟结果分析 • 丝根数（丝间隙）的变化对内爆过程的影响 • 32根铝丝阵内爆过程： • 64根铝丝阵内爆过程： • 丝等离子体在初始位置缓慢膨胀； • 丝阵形成先驱等离子体,并开始内爆； • 在内爆过程中，主体等离子体产生明显的角向融合； • 内爆到心，反弹。 • 丝等离子体在初始位置缓慢膨胀，在内爆开始之前即发生明显的角向融合； • 丝阵内爆，未观察到先驱等离子体的形成，内爆过程类似等离子体壳的内爆。 • 内爆到心，反弹。

19. 模拟结果分析 • 平均径向密度分布(N=16 VS. N=64)

20. 模拟结果分析 • 先驱等离子体（定义为70ns时刻0.6cm内部的等离子体） • 份额 先驱质量份额[cm] 丝间隙[cm] • 先驱质量占丝阵总质量的份额随着丝间隙的增大而增大，其规 • 律与丝阵消融模型的模拟结果基本一致*。 * Jun Huanget al., Phys. Plasmas 18, 042704(2011)

21. 模拟结果分析 • N=64丝阵的内爆轨迹 0D • 实验中根据0维计算得到的最大内爆速度约为5.0×107cm/s* • 由此得到最大内爆速度约为4.6×107cm/s *C. Deeney et al., Phys. Rev. E., 56(5), 5945(1997)

22. 模拟结果分析 • 磁场位形（t=65ns） N=16 N=32 N=64 • 磁场分布 当丝数N=64，丝间隙δ~0.9mm, 丝阵融合成等离子体壳，磁场具有较好的角向均匀性。此时，采用等离子体壳及一维磁场近似对丝阵进行描述是合理的。

23. 模拟结果分析 • 初始丝等离子体直径d0对等离子体角向融合过程的影响（Gi为丝阵融合成壳的丝间隙阈值） d0 =0.1mm, Gi=0.5mm d0 =0.4mm, Gi=1.0mm d0 =0.2mm, Gi=0.7mm • 丝阵融合成壳的丝间隙阈值随d0的增大而增大 • d0的确定需要对早期金属丝的电爆炸过程进行深入研究 [ T. W. L. Sanford et. al，Phys. Rev. Lett. 77, 5063 (1996)]

24. 小结与讨论 • 建立了二维(r,q)丝阵内爆动力学模型，研制了相应的电阻MHD程序，对丝阵的消融－内爆动力学过程进行了一系列研究。模拟得到了丝阵消融－内爆过程的基本物理图象，其中磁场位形和电流分布与国外的模拟结果以及实验结果符合的较好； • 对不同丝数情况下的Al丝阵Z箍缩过程进行了试算，得出以下主要结论：随着丝根数的增加： • 等离子体在内爆过程中发生明显的角向融合； • 先驱等离子占丝阵总质量的份额减少; • 磁场沿角向的不均匀性逐渐减弱； • 存在一个依赖于初始等离子体直径的丝间隙阈值，当丝间隙小于该阈值，丝阵的内爆过程呈现出等离子体壳层的特性。

25. 谢 谢！