The direct discontinuous Galerkin (DDG) viscous flux scheme for the high order spectral volume method

被引:34
作者
Kannan, R. [1 ]
Wang, Z. J. [2 ]
机构
[1] CFD Res Corp, Huntsville, AL 35805 USA
[2] Dept Aerosp Engn, Ames, IA 50011 USA
关键词
Spectral volume; Penalty; DDG; DDG2; DDG3; High order; FINITE-ELEMENT-METHOD; ESSENTIALLY NONOSCILLATORY SCHEMES; CONSERVATIVE DIFFERENCE SCHEME; UNSTRUCTURED GRIDS; LAWS; EXTENSION; SYSTEMS; EQUATIONS; FLOW;
D O I
10.1016/j.compfluid.2010.07.006
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
The direct discontinuous Galerkin (DDG) method was developed by Liu and Yan to discretize the diffusion flux. It was implemented for the discontinuous Galerkin (DG) formulation. In this paper, we perform four tasks: (i) implement the direct discontinuous Galerkin (DDG) scheme for the spectral volume method (SV) method, (ii) design and implement two variants of DDG (called DDG2 and DDG3) for the SV method, (iii) perform a Fourier type analysis on both methods when solving the 1D diffusion equation and combine the above with a non-linear global optimizer, to obtain modified constants that give significantly smaller errors (in 1D), (iv) use the above coefficients as starting points in 2D. The dissipation properties of the above schemes were then compared with existing flux formulations (local discontinuous Galerkin, Penalty and BR2). The DDG, DDG2 and DDG3 formulations were found to be much more accurate than the above three existing flux formulations. The accuracy of the DDG scheme is heavily dependent on the penalizing coefficient for the odd ordered schemes. Hence a loss of accuracy was observed even for mildly non-uniform grids for odd ordered schemes. On the other hand, the DDG2 and DDG3 schemes were mildly dependent on the penalizing coefficient for both odd and even orders and retain their accuracy even on highly irregular grids. Temporal analysis was also performed and this yielded some interesting results. The DDG and its variants were implemented in 20 (on triangular meshes) for Navier-Stokes equations. Even the non-optimized versions of the DDG displayed lower errors than the existing schemes (in 2D). In general, the DDG and its variants show promising properties and it indicates that these approaches have a great potential for higher dimension flow problems. (C) 2010 Elsevier Ltd. All rights reserved.
引用
收藏
页码:2007 / 2021
页数:15
相关论文
共 34 条
[1]   ON ESSENTIALLY NONOSCILLATORY SCHEMES ON UNSTRUCTURED MESHES - ANALYSIS AND IMPLEMENTATION [J].
ABGRALL, R .
JOURNAL OF COMPUTATIONAL PHYSICS, 1994, 114 (01) :45-58
[2]   AN INTERIOR PENALTY FINITE-ELEMENT METHOD WITH DISCONTINUOUS ELEMENTS [J].
ARNOLD, DN .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 1982, 19 (04) :742-760
[3]  
Barth T., 1990, 28 AER SCI M AM I AE
[4]   A high-order accurate discontinuous finite element method for the numerical solution of the compressible Navier-Stokes equations [J].
Bassi, F ;
Rebay, S .
JOURNAL OF COMPUTATIONAL PHYSICS, 1997, 131 (02) :267-279
[5]   High-order accurate discontinuous finite element solution of the 2D Euler equations [J].
Bassi, F ;
Rebay, S .
JOURNAL OF COMPUTATIONAL PHYSICS, 1997, 138 (02) :251-285
[6]  
BAUMANN CE, 1997, THESIS U TEXAS AUSTI
[7]  
Brezzi F, 2000, NUMER METH PART D E, V16, P365, DOI 10.1002/1098-2426(200007)16:4<365::AID-NUM2>3.0.CO
[8]  
2-Y
[10]   The local discontinuous Galerkin method for time-dependent convection-diffusion systems [J].
Cockburn, B ;
Shu, CW .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 1998, 35 (06) :2440-2463