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 条
[1]   Self-adjusting, positivity preserving high order schemes for hydrodynamics and magnetohydrodynamics [J].
Balsara, Dinshaw S. .
JOURNAL OF COMPUTATIONAL PHYSICS, 2012, 231 (22) :7504-7517
[2]   Divergence-free reconstruction of magnetic fields and WENO schemes for magnetohydrodynamics [J].
Balsara, Dinshaw S. .
JOURNAL OF COMPUTATIONAL PHYSICS, 2009, 228 (14) :5040-5056
[3]   A staggered mesh algorithm using high order Godunov fluxes to ensure solenoidal magnetic fields in magnetohydrodynamic simulations [J].
Balsara, DS ;
Spicer, DS .
JOURNAL OF COMPUTATIONAL PHYSICS, 1999, 149 (02) :270-292
[4]   Second-order-accurate schemes for magnetohydrodynamics with divergence-free reconstruction [J].
Balsara, DS .
ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES, 2004, 151 (01) :149-184
[5]   Maintaining pressure positivity in magnetohydrodynamic simulations [J].
Balsara, DS ;
Spicer, D .
JOURNAL OF COMPUTATIONAL PHYSICS, 1999, 148 (01) :133-148
[6]   An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws [J].
Borges, Rafael ;
Carmona, Monique ;
Costa, Bruno ;
Don, Wai Sun .
JOURNAL OF COMPUTATIONAL PHYSICS, 2008, 227 (06) :3191-3211
[7]  
Bouchut F, 2007, NUMER MATH, V108, P7, DOI [10.1007/s00211-007-0108-8, 10.1007/S00211-007-0108-8]
[8]   A multiwave approximate Riemann solver for ideal MHD based on relaxation II: numerical implementation with 3 and 5 waves [J].
Bouchut, Francois ;
Klingenberg, Christian ;
Waagan, Knut .
NUMERISCHE MATHEMATIK, 2010, 115 (04) :647-679
[9]   THE EFFECT OF NONZERO-DEL.B ON THE NUMERICAL-SOLUTION OF THE MAGNETO-HYDRODYNAMIC EQUATIONS [J].
BRACKBILL, JU ;
BARNES, DC .
JOURNAL OF COMPUTATIONAL PHYSICS, 1980, 35 (03) :426-430
[10]   ENTROPY STABLE FINITE VOLUME SCHEME FOR IDEAL COMPRESSIBLE MHD ON 2-D CARTESIAN MESHES [J].
Chandrashekar, Praveen ;
Klingenberg, Christian .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2016, 54 (02) :1313-1340