基于有限体积法数值模拟双撕裂模非线性演化-核聚变与等离子体物理.PDFVIP

第 37 卷 第 3 期 核 聚 变 与 等 离 子 体 物 理 Vol.37, No.3 2 0 1 7 年 9 月 Nuclear Fusion and Plasma Physics Sep. 2017 文章编号：0254−6086(2017)03−0257−05 DOI: 10.16568/j.0254-6086.201703002 基于有限体积法数值模拟双撕裂模非线性演化 1 *1 2 1 雷晓晨 ，李新霞 ，马 骏 ，龚学余 (1. 南华大学核科学技术学院，衡阳 421001 ；2. 中国科学院等离子体物理研究所，合肥 230026) 摘 要：以磁流体理论为基础，采用基于有限体积法的通量差分分裂格式数值求解具有双曲保守律形式的电 阻磁流体方程组。编写 C++程序对平板几何位形下的等离子体双撕裂模进行了长时间数值模拟，得到双撕裂模不 稳定性的演化图景，捕捉到了双撕裂模非线性发展过程中磁场重联的几个典型阶段，讨论了等离子体电阻和两个 有理面之间的距离对双撕裂模不稳定性非线性发展的影响。为研究磁流体动力学提供了一种可行的高精度数值算 法。 关键词：有限体积法；双撕裂模不稳定性；数值模拟 中图分类号：TL61+2.2 文献标志码 ：A 1 引言 高保真度，已被成功的应用在航空设计、天气预报、 [6, 7] 未来的托卡马克聚变反应堆要想实现稳态运 水力学等研究领域 。 行，磁流体力学(MHD)中的非线性效应将变得更加 本文基于有限体积法的通量差分分裂格式 重要。然而，磁流体力学中非线性效应常常伴随着 (FVM-FDS) ，联合高阶多态黎曼求解器(HLL 和 异常复杂的物理现象，至今仍缺乏充分的理论解释 HLLC) 、高精度重构格式(2 阶 MUSCL 和 5 阶 [1, 2] 。在目前高温等离子体物理研究领域，对磁流 WENO)和总变差递减龙格-库塔格式(TVD-RK) ，针 体动力学 (MHD)行为数值计算与模拟已然成为最 对 2 维电阻 MHD 问题(双撕裂模不稳定性)编写 重要的研究手段之一。近年来，研究人员开发如 C++程序数值求解相应的磁流体力学方程组，得到 NIMROD 、M3D 等磁流体动力学代码[3, 4] ，在解释 平板几何位形下电阻引起的磁流体双撕裂模不稳 等离子体约束、输运等方面取得了极大的进展。同 定性的发展演化物理图景，捕捉得到其非线性发展 时，这些代码普遍基于传统的算符离散算法(有限差 过程中磁场重联的 4 个典型阶段。此外，本文还讨 分法)来编写。该算法虽然数学表达直观，但是在较 论了等离子体电阻和两个有理面之间的距离对双 复杂的几何位形和非线性物理现象中难以保证离 撕裂模不稳定性非线性发展的影响。 散化的方程仍然满足守恒规律，并且计算精度有待 进一步提高。基于有限体积法的通量差分分裂格式 2 数值算法及物理模型 (FDS)是在各种计算流体力学(CFD)算法中一个非 2.1 数值算法 常重要的分支[5] ，因其理想精确度及高精度黎曼近 基于有限体积法的基本思路，采用 FDS(通量 似解，可以保证长时间尺度下模拟结果的精确度和 差分分裂)方法来数值求解具有双曲型守恒律形式 收稿日期：2016−07−13；修订日期：2017−06−07 基金项目：国