Construction of low dissipative high-order well-balanced filter schemes for non-equilibrium flows

被引:16
作者
Wang, Wei [2 ]
Yee, H. C. [3 ]
Sjoegreen, Bjoern [4 ]
Magin, Thierry [2 ]
Shu, Chi-Wang [1 ]
机构
[1] Brown Univ, Div Appl Math, Providence, RI 02912 USA
[2] Stanford Univ, Ctr Turbulence Res, Stanford, CA 94305 USA
[3] NASA, Ames Res Ctr, Moffett Field, CA 94035 USA
[4] Lawrence Livermore Natl Lab, Livermore, CA 94551 USA
关键词
High-order filter methods; WENO schemes; Well-balanced schemes; Non-equilibrium flow; Chemical reactions; 1D turbulence/shock interactions; DIFFERENCE WENO SCHEMES; THERMAL RATE CONSTANTS; EFFICIENT IMPLEMENTATION; HYPERBOLIC SYSTEMS; CONSERVATION-LAWS; SOURCE TERMS;
D O I
10.1016/j.jcp.2010.04.033
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
The goal of this paper is to generalize the well-balanced approach for non-equilibrium flow studied by Wang et al. (2009)129] to a class of low dissipative high-order shock-capturing filter schemes and to explore more advantages of well-balanced schemes in reacting flows. More general 1D and 2D reacting flow models and new examples of shock turbulence interactions are provided to demonstrate the advantage of well-balanced schemes. The class of filter schemes developed by Yee et al. (1999) 1331, Sjogreen and Yee (2004) [27] and Yee and Sjogreen (2007) [381 consist of two steps, a full time step of spatially high-order non-dissipative base scheme and an adaptive non-linear filter containing shock-capturing dissipation. A good property of the filter scheme is that the base scheme and the filter are stand-alone modules in designing. Therefore, the idea of designing a well-balanced filter scheme is straightforward, i.e. choosing a well-balanced base scheme with a well-balanced filter (both with high-order accuracy). A typical class of these schemes shown in this paper is the high-order central difference schemes/predictor-corrector (PC) schemes with a high-order well-balanced WENO filter. The new filter scheme with the well-balanced property will gather the features of both filter methods and well-balanced properties: it can preserve certain steady-state solutions exactly; it is able to capture small perturbations, e.g. turbulence fluctuations; and it adaptively controls numerical dissipation. Thus it shows high accuracy, efficiency and stability in shock/turbulence interactions. Numerical examples containing 1D and 2D smooth problems, 1D stationary contact discontinuity problem and 1D turbulence/shock interactions are included to verify the improved accuracy, in addition to the well-balanced behavior. (C) 2010 Elsevier Inc. All rights reserved.
引用
收藏
页码:4316 / 4335
页数:20
相关论文
共 40 条
[1]  
BAYLISS A, 1985, AIAA FLUID DYN PLASM
[2]   Thermal rate constants of the O-2+N->NO+O reaction based on the (2)A' and (4)A' potential-energy surfaces [J].
Bose, D ;
Candler, GV .
JOURNAL OF CHEMICAL PHYSICS, 1997, 107 (16) :6136-6145
[3]   Thermal rate constants of the N-2+O->NO+N reaction using ab initio (3)A'' and (3)A' potential energy surfaces [J].
Bose, D ;
Candler, GV .
JOURNAL OF CHEMICAL PHYSICS, 1996, 104 (08) :2825-2833
[4]  
CIMENT M, 1975, MATH COMPUT, V29, P985, DOI 10.1090/S0025-5718-1975-0416049-2
[5]   A well-balanced flux-vector splitting scheme designed for hyperbolic systems of conservation laws with source terms [J].
Gosse, L .
COMPUTERS & MATHEMATICS WITH APPLICATIONS, 2000, 39 (9-10) :135-159
[6]  
GOTTLIEB D, 1976, MATH COMPUT, V30, P703, DOI 10.1090/S0025-5718-1976-0443362-6
[7]   A well-balanced scheme for the numerical processing of source terms in hyperbolic equations [J].
Greenberg, JM ;
Leroux, AY .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 1996, 33 (01) :1-16
[8]  
GURVICH LV, 1989, THERMODYNAMIC PROP 2, V1
[9]  
GURVICH LV, 1989, THERMODYNAMIC PROP 1, V1
[10]  
HARTEN A, 1978, MATH COMPUT, V32, P363, DOI 10.2307/2006149