三维拉氏理想磁流体数值模拟方法

计算物理 ›› 2020, Vol. 37 ›› Issue (4): 403-412.

所属专题：第十届全国青年计算物理会议优秀论文

三维拉氏理想磁流体数值模拟方法

徐骁^1,2, 高志明¹, 戴自换¹

1. 北京应用物理与计算数学研究所邓稼先创新研究中心, 北京 100088;
2. 中国工程物理研究院研究生院, 北京 100088

收稿日期:2019-08-14 修回日期:2019-11-19 出版日期:2020-07-25 发布日期:2020-07-25
作者简介:徐骁(1989-),研究方向为三维磁流体方程组的ALE算法研究,E-mail:xu_xiao@iapcm.ac.cn
基金资助:
国家自然科学基金（11771052，11871009，11975054）及中物院创新发展基金（CX2019028）资助项目

3D Lagrangian Methods for Ideal Magnetohydrodynamics on Unstructured Meshes

XU Xiao^1,2, GAO Zhiming¹, DAI Zihuan¹

1. Institute of Applied Physics and Computational Mathematics, Beijing 100088, China;
2. Graduate School of China Academy of Engineering Physics, Beijing 100088, China

Received:2019-08-14 Revised:2019-11-19 Online:2020-07-25 Published:2020-07-25

摘要/Abstract

摘要： 在Z箍缩驱动ICF过程中，磁场流体耦合作用是其整个物理过程中非常重要的部分.针对Z箍缩过程中多介质、大变形的特点，发展了三维相容拉氏理想磁流体交错型以及单元中心型格式，两种格式均具有一阶时间与空间精度.通过数值算例，验证其精度和强壮性，并比较分析两种格式的特点.

关键词: 三维拉氏理想磁流体, 交错格式, 单元中心格式

Abstract: Interaction between magnetic field and fluids is crucial in Z pinch process. Since Z pinch involves with multi materials and severe deformation, we develop 3D compatible Lagrangian staggered and cell-centered schemes for ideal MHD on unstructured meshes. Both of the schemes are of first order in space and time discretization. Accuracy and robustness of the schemes are validated with typical numerical tests.

Key words: 3D Lagrangian magnetohydrodynamics, staggered scheme, cell-centered scheme

中图分类号:

O241

徐骁, 高志明, 戴自换. 三维拉氏理想磁流体数值模拟方法[J]. 计算物理, 2020, 37(4): 403-412.

XU Xiao, GAO Zhiming, DAI Zihuan. 3D Lagrangian Methods for Ideal Magnetohydrodynamics on Unstructured Meshes[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2020, 37(4): 403-412.

参考文献

[1] 贾祖朋, 张树道, 蔚喜军. 多介质流体动力学计算方法[M]. 北京:科学出版社, 2014, 1-17.
[2] TOTH G. The Δ·B=0 constraint in shock-capturing magnetohydrodynamics codes[J]. Comput Phys, 2000, 161(2):605-652.
[3] WENSENBURG M. Efficient MHD Riemann solvers for simulations on unstructured triangular grids[J]. J Numer Math, 2002, 10(1):37-71.
[4] GALLICE G. Positive and entropy stable Godunov-type schemes for gas dynamicsand MHD equations in Lagrangian or Eulerian coordinates[J]. Numer Math, 2003, 94(4):673-713.
[5] Von NEUMANN J, RICHTMYER R D. A method for the numerical calculations of hydrodynamical shocks[J]. J Appl Phys, 1950, 21(3):232-238.
[6] CARAMANA E J, BURTON D E, SHASHKOV M J, et al. The construction of compatible hydrodynamics algorithms utilizing conservation of total energy[J]. J Comput Phys, 1998, 146(1):227-262.
[7] CARAMANA E J, SHASHKOV M J, WHALEN P P. Formulations of artificial viscosity for multi-dimensional shock wave computations[J]. J Comput Phys, 1998, 144(1):70-97.
[8] CAMPBELL J C, SHASHKOV M J. A tensor artificial viscosity using a mimetic finite difference algorithm[J]. J Comput Phys, 2001, 172(2):739-765.
[9] KOLEV T V, RIEBEN R N. A tensor artificial viscosity using a finite element approach[J]. J Comput Phys, 2009, 228:8336-8366.
[10] DAI Z, WU J, LIN Z, et al. Subzonal pressure methods in Lagrangian algorithm of two-dimensional three-temperature radiation hydrodynamics[J]. Chinese Journal of Computational Physics, 2010, 27(3):326-334.
[11] GODUNOV S K. A difference scheme for numerical computation of discontinue solution of hydrodynamic equations[J]. Mathematics Sbornik, 1959, 47:271-306.
[12] MAIRE P H, ABGRALL R, BREIL J, et al. A cell-centered Lagrangian scheme for two-dimensional compressible flows problems[J]. SIAM J Sci Comput, 2007, 29:1782-1824.
[13] GEORGES G, BREIL J, MAIRE P H. A 3D GCL compatible cell-centered Lagrangian scheme for solving gas dynamics equations[J]. J Comput Phys, 2016, 305:921-941.
[14] LOUBERE R, MAIRE P H, VACHAL P, 3D staggered Lagrangian hydrodynamics scheme with cell-centered Riemann solver-based artificial viscosity[J]. Int J Numer Meth Fluids, 2013, 72(1):22-42.
[15] SUN Y, YU M, REN Y, et al. A finite volume method for 2D Lagrangian hydrodynamics based on characteristics theory[J]. Chinese Journal of Computational Physics, 2011, 28(1):19-26.
[16] SUN Y, JIA A, YU M, et al. A second order Lagrangian scheme based on characteristics theory for two-dimensional compressible flows[J]. Chinese Journal of Computational Physics, 2012, 29(6):791-798.
[17] XU X, GAO Z, DAI Z. A 3D staggered Lagrangian scheme for ideal magnetohydrodynamics equations on unstructured meshes[J]. Int J Numer Meth Fluids, 2019, 90:584-602.
[18] XU X, DAI Z, GAO Z. A 3D cell centered Lagrangian scheme for the ideal magnetohydrodynamics equations on unstructured meshes[J]. Comput Methods Appl Mech Engrg, 2018, 342:490-508.
[19] BENSON D J. Computational methods in Lagrangian and Eulerian hydrocodes[J]. Comput Method Appl Mech Eng, 1992, 99(2-3):235-394.
[20] MAIRE P H, NKONGA B. Multi-scale Godunov-type method for cell-centered discrete Lagrangian hydrodynamics[J]. J Comput Phys, 2008, 228:799-821.
[21] BRIO M, WU C C. An upwind difference scheme for the equations of ideal magnetohydrodynamics[J]. J Comput Phys, 1988, 75(2):400-422.
[22] BALSARA D S, SPICER D. A staggered mesh algorithm using high order Godunov fluxs to ensure solenoidal magnetic fields in magnetohydrodynamic simulations[J]. J Comput Phys, 1999, 149(1):270-292.
[23] BALSARA D S. Self-adjusting, positivity preserving high order schemes for hydrodynamics and magnetohydrodynamics[J]. J Comput Phys, 2012, 231:7504-7517.

三维拉氏理想磁流体数值模拟方法

3D Lagrangian Methods for Ideal Magnetohydrodynamics on Unstructured Meshes

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

作者中心

审稿中心

期刊浏览

期刊介绍

[1]	魏雪丹, 戴厚平, 李梦军, 郑洲顺. 一维空间Riesz分数阶对流扩散方程的格子Boltzmann方法[J]. 计算物理, 2021, 38(6): 683-692.
[2]	房尧立, 王一. 一种求解二维辐射流体力学方程组的显隐式格式[J]. 计算物理, 2021, 38(4): 401-417.
[3]	徐潜龙, 李晔. 有限深水域中三维频域自由面格林函数及其导数的计算方法[J]. 计算物理, 2021, 38(1): 35-46.
[4]	赵菲, 盛志强, 袁光伟. 基于二阶格式的有限体积保正格式[J]. 计算物理, 2020, 37(4): 379-392.
[5]	张林, 葛永斌. 二维半线性扩散反应方程的高精度全隐格式及其多重网格方法[J]. 计算物理, 2020, 37(3): 307-319.
[6]	刘妍, 茅德康. 相容拉氏方法中的ADRS近似Riemann解算器[J]. 计算物理, 2020, 37(2): 140-152.
[7]	丁琪, 尚月强. 非定常Navier-Stokes方程基于两重网格离散的有限元并行算法[J]. 计算物理, 2020, 37(1): 10-18.
[8]	刘全, 倪国喜, 牛霄, 胡军. Vlasov-Poisson方程半拉格朗日守恒算法[J]. 计算物理, 2020, 37(1): 19-25.
[9]	赵国忠, 蔚喜军, 郭虹平, 董自明. 求解含有高阶导数偏微分方程的局部间断Petrov-Galerkin方法[J]. 计算物理, 2019, 36(5): 517-532.
[10]	孔令海, 孔令波, 许海波, 贾清刚. 含Laplace-Gauss型混合噪声图像二阶正则化重建方法[J]. 计算物理, 2019, 36(3): 280-290.
[11]	郭子滔, 冯仁忠. 一种高精度的修正Hermite-ENO格式[J]. 计算物理, 2019, 36(2): 141-152.
[12]	戴自换. ALE计算中基于声速分布的网格优化技术[J]. 计算物理, 2019, 36(1): 15-24.
[13]	徐丞君, 徐胜利. SPH二阶粒子近似光滑函数的充分条件及数值验证[J]. 计算物理, 2019, 36(1): 25-38.
[14]	符芳芳, 孔令华, 王兰, 徐远, 曾展宽. 一维Gross-Pitaevskii方程的高阶紧致分裂步多辛格式[J]. 计算物理, 2018, 35(6): 657-667.
[15]	张守慧, 梁栋. 二维抛物型问题的Strang型交替分段区域分裂格式[J]. 计算物理, 2018, 35(4): 413-428.