Implicit discontinuous Galerkin method for RANS simulation utilizing pointwise relaxation algorithm

An efficient implicit procedure for the Discontinuous Galerkin (DG) method is developed utilizing a pointwise relaxation algorithm. In the pointwise relaxation, those contributions from the degrees of freedom in own computational cell are accounted for in the implicit matrix inversion. The resulting scheme is shown to be stable with very large CFL numbers for both the Euler and the Navier-Stokes equations for typical test problems. In order to achieve a faster convergence, efforts are also made to reduce computing time of the present method by utilizing a p-multigrid scheme and also by solving a simplified matrix instead of a fully loaded dense matrix in the implicit matrix inversion. A superior performance of the present implicit DG method on the parallel computer using up to 128 PEs is shown in terms of readily achievable scalability and high parallel efficiency. The RANS simulation of turbulent flowfield over AGARD-B model is carried out to show the convergence property and numerical stability of the present implicit DG method for engineering applications. AMS subject classifications: 76M10, 65M60, 65D30, 65B99

[1]  Vipin Kumar,et al.  A Fast and High Quality Multilevel Scheme for Partitioning Irregular Graphs , 1998, SIAM J. Sci. Comput..

[2]  T. Barth.  Aspects of Unstructured Grids and Finite-Volume Solvers for the Euler and Navier-Stokes Equations , 1994 .

[3]  Ralf Hartmann,et al.  Adaptive discontinuous Galerkin methods with shock‐capturing for the compressible Navier–Stokes equations , 2006 .

[4]  V. Schmitt,et al.  Pressure distributions on the ONERA M6 wing at transonic Mach numbers , 1979 .

[5]  Ralf Hartmann,et al.  Symmetric Interior Penalty DG Methods for the Compressible Navier-Stokes Equations II: Goal--Oriented A Posteriori Error Estimation , 2005 .

[6]  Vít Dolejší,et al.  A semi-implicit discontinuous Galerkin finite element method for the numerical solution of inviscid compressible flow , 2004 .

[7]  M. Y. Hussaini,et al.  An efficient implicit discontinuous spectral Galerkin method , 2001 .

[8]  S. Rebay,et al.  Numerical evaluation of two discontinuous Galerkin methods for the compressible Navier–Stokes equations , 2002 .

[9]  Antony Jameson,et al.  Lower-upper implicit schemes with multiple grids for the Euler equations , 1987 .

[10]  P. Spalart.  A One-Equation Turbulence Model for Aerodynamic Flows , 1992 .

[11]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[12]  R. Hartmann,et al.  Symmetric Interior Penalty DG Methods for the CompressibleNavier-Stokes Equations I: Method Formulation , 2005 .

[13]  V. Dolejší,et al.  Semi-Implicit Interior Penalty Discontinuous Galerkin Methods for Viscous Compressible Flows , 2008 .

[14]  George Em Karniadakis,et al.  A NEW TRIANGULAR AND TETRAHEDRAL BASIS FOR HIGH-ORDER (HP) FINITE ELEMENT METHODS , 1995 .

[15]  Dimitri J. Mavriplis,et al.  High-order discontinuous Galerkin methods using an hp-multigrid approach , 2006, J. Comput. Phys..

[16]  Meng-Sing Liou,et al.  A Flux Splitting Scheme with High-Resolution and Robustness for Discontinuities(Proceedings of the 12th NAL Symposium on Aircraft Computational Aerodynamics) , 1994 .

[17]  Rainald Löhner,et al.  Fast p-Multigrid Discontinuous Galerkin Method for Compressible Flows at All Speeds , 2008 .

[18]  Moshe Dubiner.  Spectral methods on triangles and other domains , 1991 .

[19]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .