A NEW DISCRETELY DIVERGENCE-FREE POSITIVITY-PRESERVING HIGH-ORDER FINITE VOLUME METHOD FOR IDEAL MHD EQUATIONS

被引:6
作者
Ding, Shengrong [1 ,2 ]
Wu, Kailiang [1 ,2 ]
机构
[1] Southern Univ Sci & Technol, SUSTech Int Ctr Math, Shenzhen 518055, Peoples R China
[2] Southern Univ Sci & Technol, Dept Math, Shenzhen 518055, Peoples R China
基金
中国国家自然科学基金;
关键词
compressible MHD; positivity-preserving; divergence-free; finite volume method; high-order accuracy; hyperbolic conservation laws; DISCONTINUOUS GALERKIN METHODS; GEOMETRIC QUASI-LINEARIZATION; CONSTRAINED-TRANSPORT METHOD; DIFFERENCE SCHEME; MAGNETOHYDRODYNAMICS; RELAXATION; LIMITERS;
D O I
10.1137/23M1562081
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
This paper proposes and analyzes a novel efficient high -order finite volume method for the ideal magnetohydrodynamics (MHD). As a distinctive feature, the method simultaneously preserves two critical physical constraints: a discretely divergence -free (DDF) constraint on the magnetic field and the positivity -preserving (PP) property, which ensures the positivity of density, pressure, and internal energy. To enforce the DDF condition in each cell, we design a new discrete projection approach that projects the reconstructed point values at the cell interface into a DDF space, without using any approximation polynomials. This projection method is highly efficient, easy to implement, and particularly suitable for the high -order finite volume methods that return only the point values (no explicit approximation polynomials) in the reconstruction. Moreover, we also develop a new finite volume framework for constructing provably PP schemes for the ideal MHD system. The framework comprises the discrete projection technique, a suitable approximation to the Godunov-Powell source terms, and a simple PP limiter. We provide rigorous analysis of the PP property of the proposed finite volume method, demonstrating that the DDF condition and the proper approximation to the source terms eliminate the impact of magnetic divergence terms on the PP property. The analysis is challenging due to the internal energy function's nonlinearity and the intricate relationship between the DDF and PP properties. To address these challenges, we adopt the recently developed geometric quasilinearization approach [K. Wu and C. -W. Shu, SIAM Rev., 65 (2023), pp. 1031-1073], which transforms a nonlinear constraint into a family of linear constraints. Finally, we validate the effectiveness of the proposed method through several benchmark and demanding numerical examples. The results demonstrate that the proposed method is robust, accurate, and highly effective, confirming the significance of the proposed DDF projection and PP techniques.
引用
收藏
页码:A50 / A79
页数:30
相关论文
共 59 条
[41]   The del•B=0 constraint in shock-capturing magnetohydrodynamics codes [J].
Tóth, G .
JOURNAL OF COMPUTATIONAL PHYSICS, 2000, 161 (02) :605-652
[42]   POSITIVITY PRESERVING LIMITERS FOR TIME-IMPLICIT HIGHER ORDER ACCURATE DISCONTINUOUS GALERKIN DISCRETIZATIONS [J].
Van der Vegt, J. J. W. ;
Xia, Yinhua ;
Xu, Yan .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2019, 41 (03) :A2037-A2063
[43]   A robust numerical scheme for highly compressible magnetohydrodynamics: Nonlinear stability, implementation and tests [J].
Waagan, K. ;
Federrath, C. ;
Klingenberg, C. .
JOURNAL OF COMPUTATIONAL PHYSICS, 2011, 230 (09) :3331-3351
[44]   A positive MUSCL-Hancock scheme for ideal magnetohydrodynamics [J].
Waagan, K. .
JOURNAL OF COMPUTATIONAL PHYSICS, 2009, 228 (23) :8609-8626
[45]   Geometric Quasilinearization Framework for Analysis and Design of Bound-Preserving Schemes [J].
Wu, Kailiang ;
Shu, Chi-Wang .
SIAM REVIEW, 2023, 65 (04) :1031-1073
[46]   PROVABLY POSITIVE CENTRAL DISCONTINUOUS GALERKIN SCHEMES VIA GEOMETRIC QUASILINEARIZATION FOR IDEAL MHD EQUATIONS [J].
Wu, Kailiang ;
Jiang, Haili ;
Shu, Chi-Wang .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2023, 61 (01) :250-285
[47]   Provably physical-constraint-preserving discontinuous Galerkin methods for multidimensional relativistic MHD equations [J].
Wu, Kailiang ;
Shu, Chi-Wang .
NUMERISCHE MATHEMATIK, 2021, 148 (03) :699-741
[48]   Provably positive high-order schemes for ideal magnetohydrodynamics: analysis on general meshes [J].
Wu, Kailiang ;
Shu, Chi-Wang .
NUMERISCHE MATHEMATIK, 2019, 142 (04) :995-1047
[49]   A PROVABLY POSITIVE DISCONTINUOUS GALERKIN METHOD FOR MULTIDIMENSIONAL IDEAL MAGNETOHYDRODYNAMICS [J].
Wu, Kailiang ;
Shu, Chi-Wang .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2018, 40 (05) :B1302-B1329
[50]   POSITIVITY-PRESERVING ANALYSIS OF NUMERICAL SCHEMES FOR IDEAL MAGNETOHYDRODYNAMICS [J].
Wu, Kailiang .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2018, 56 (04) :2124-2147