辐射磁流体力学不稳定性的数值模拟方法研究

Radiation magnetohydrodynamic instability is an important fundamental problem in Z-Pinch research with enough challenge. Initiated by the need of numerically simulating the plasma implosion driven by the Z-pinch, this project will study on the high-accuracy and high-performance schemes for radiation magnetohydrodynamic instability on two-dimensional complex meshes, which include the high-accuracy Lagrangian hydrodynamic algorithm, the high order discrete schemes for the diffusion operator, the domain decomposition method and the rezoning and remapping method. This project will also develop an application program for simulating the Z-Pinch. Algorithms and the program developted by this project, can hopefully improve our numerical simulation ability for Z-pinch. At the same time, it can provide a beneficial reference for others radiation hydrodynamic problems such as the ablative instability in the laser-drived ICF.

辐射磁流体力学不稳定性是Z-Pinch 研究中的一个十分重要且富有挑战的基础性问题。本项目以Z-Pinch 等离子体内爆的数值模拟需求为牵引，研究二维复杂网格体系下辐射磁流体力学不稳定性的具有高效率和高置信度的数值模拟方法，包括复杂网格体系下的高精度拉氏流体动力学计算格式、高精度扩散算子离散格式、区域分解算法以及网格优化与物理量重映方法，并研制相关应用程序，为提升我国的Z-pinch 等离子体内爆数值模拟能力做出贡献，同时为其它辐射流体力学相关问题如激光驱动ICF 中的烧蚀不稳定性等的数值模拟提供有益的借鉴。

辐射磁流体力学不稳定性是Z-Pinch 研究中的一个十分重要且富有挑战的基础性问题。针对我国Z-Pinch等离子体内爆数值模拟的实际需求，项目组系统地开展了大量创新性研究工作，提出了二维复杂网格体系下辐射磁流体力学不稳定性问题的具有较高置信度的一整套数值模拟方法，包括精确保持磁场散度约束条件的单元中心型拉氏磁流体力学离散格式，扩散方程（组）的三种单元中心型有限体积格式，扩散方程（组）的一种收敛率不依赖于Robin 参数的区域分解算法，一种基于声速分布的拉氏多边形网格优化技术，研制了大型应用程序MARED-U并应用于我国Z-pinch等离子体内爆物理的数值模拟研究，取得了令人满意的研究成果。