Positivity-preserving high order finite difference WENO schemes for compressible Navier-Stokes equations

被引:26
作者
Fan, Chuan [1 ]
Zhang, Xiangxiong [2 ]
Qiu, Jianxian [3 ,4 ]
机构
[1] Xiamen Univ, Sch Math Sci, Xiamen 361005, Fujian, Peoples R China
[2] Purdue Univ, Dept Math, W Lafayette, IN 47907 USA
[3] Xiamen Univ, Sch Math Sci, Xiamen 361005, Fujian, Peoples R China
[4] Xiamen Univ, Fujian Prov Key Lab Math Modeling & High Performa, Xiamen 361005, Fujian, Peoples R China
关键词
Positivity-preserving; Compressible; Navier-Stokes equations; WENO; Finite difference; High order accuracy; NUMERICAL-SIMULATION; WEIGHTS;
D O I
10.1016/j.jcp.2022.111446
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
In this paper, we construct a high order weighted essentially non-oscillatory (WENO) finite difference discretization for compressible Navier-Stokes (NS) equations, which is rendered positivity-preserving of density and internal energy by a positivity-preserving flux splitting and a scaling positivity-preserving limiter. The novelty of this paper is WENO reconstruction performed on variables from a positivity-preserving convection diffusion flux splitting, which is different from conventional WENO schemes solving compressible NS equations. The core advantages of our proposed method are robustness and efficiency, which especially are suitable for solving tough demanding problems of both compressible Euler and NS equation including low density and low pressure flow regime. Moreover, in terms of computational cost, it is more efficient and easier to implement and extend to multi-dimensional problems than the positivity-preserving high order discontinuous Galerkin schemes and finite volume WENO scheme for solving compressible NS equations on rectangle domain. Benchmark tests demonstrate that the proposed positivity-preserving WENO schemes are high order accuracy, efficient and robust without excessive artificial viscosity for demanding problems involving with low density, low pressure, and fine structure. (c) 2022 Elsevier Inc. All rights reserved.
引用
收藏
页数:21
相关论文
共 36 条
[1]   ANALYSIS OF WENO SCHEMES FOR FULL AND GLOBAL ACCURACY [J].
Arandiga, F. ;
Baeza, A. ;
Belda, A. M. ;
Mulet, P. .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2011, 49 (02) :893-915
[2]   On the choice of wavespeeds for the HLLC Riemann solver [J].
Batten, P ;
Clarke, N ;
Lambert, C ;
Causon, DM .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 1997, 18 (06) :1553-1570
[3]   The Runge-Kutta discontinuous Galerkin method for conservation laws V - Multidimensional systems [J].
Cockburn, B ;
Shu, CW .
JOURNAL OF COMPUTATIONAL PHYSICS, 1998, 141 (02) :199-224
[4]   ON GODUNOV-TYPE METHODS NEAR LOW-DENSITIES [J].
EINFELDT, B ;
MUNZ, CD ;
ROE, PL ;
SJOGREEN, B .
JOURNAL OF COMPUTATIONAL PHYSICS, 1991, 92 (02) :273-295
[5]   Positivity-preserving high order finite volume hybrid Hermite WENO schemes for compressible Navier-Stokes equations [J].
Fan, Chuan ;
Zhang, Xiangxiong ;
Qiu, Jianxian .
JOURNAL OF COMPUTATIONAL PHYSICS, 2021, 445
[6]  
Fedkiw RP, 1999, J COMPUT PHYS, V152, P457, DOI 10.1006/jcph.1999.6136
[7]   NUMERICAL SIMULATION OF THE XZ TAURI SUPERSONIC ASTROPHYSICAL JET [J].
Gardner, Carl L. ;
Dwyer, Steven J. .
ACTA MATHEMATICA SCIENTIA, 2009, 29 (06) :1677-1683
[8]  
Grapsas D., 2016, SMAI J. Comput. Math., V2, P51
[9]   Positivity of flux vector splitting schemes [J].
Gressier, J ;
Villedieu, P ;
Moschetta, JM .
JOURNAL OF COMPUTATIONAL PHYSICS, 1999, 155 (01) :199-220
[10]   Second-order invariant domain preserving approximation of the compressible Navier-Stokes equations [J].
Guermond, Jean-Luc ;
Maier, Matthias ;
Popav, Bojan ;
Tomas, Ignacio .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2021, 375 (375)