A p-multigrid spectral difference method with explicit and implicit smoothers on unstructured triangular grids

The convergence of high-order methods, such as recently developed Spectral Difference (SD) method, can be accelerated using both implicit temporal advancement and a p-Multigrid approach. A good combination of these two can significantly improve the efficiency of the SD method for steady flow problems. In this paper, we demonstrate a p-Multigrid approach developed for a 2D Euler solver using the SD method. A fast preconditioned Lower-Upper Symmetric Gauss Seidel (LU-SGS) implicit method is developed and tested for both linear scalar, and nonlinear Euler equations. Meanwhile, a five-stage Runge-Kutta explicit method is employed for comparison. We are able to achieve significant speedup (up to two orders) using both p-Multigrid method and the implicit smoother while maintaining very stable convergence property and nearly ideal order of accuracy. The numerical results are very promising, and indicate that the approach has great potential for 3D flow problems.

[1]  S. Ta'asan.  Multigrid one-shot methods and design strategy , 1997 .

[2]  Zhi J. Wang,et al.  Fast, Block Lower-Upper Symmetric Gauss-Seidel Scheme for Arbitrary Grids , 2000 .

[3]  Chunlei Liang,et al.  Extension of the SD Method to Viscous Flow on Unstructured Grids , 2009 .

[4]  Thomas A. Manteuffel,et al.  Algebraic multigrid for higher-order finite elements , 2005 .

[5]  Anthony T. Patera,et al.  Spectral element multigrid. I. Formulation and numerical results , 1987 .

[6]  Rainald Löhner,et al.  A p-multigrid discontinuous Galerkin method for the Euler equations on unstructured grids , 2006 .

[7]  D. Gaitonde,et al.  Behavior of linear reconstruction techniques on unstructured meshes , 1995 .

[8]  Antony Jameson,et al.  Spectral Difference Method for Unstructured Grids II: Extension to the Euler Equations , 2007, J. Sci. Comput..

[9]  B. T. Helenbrook,et al.  Application of “ p ”-multigrid to discontinuous Galerkin formulations of the Poisson equation , 2008 .

[10]  David L. Darmofal,et al.  p-Multigrid solution of high-order discontinuous Galerkin discretizations of the compressible Navier-Stokes equations , 2005 .

[11]  D. Kopriva.  A Conservative Staggered-Grid Chebyshev Multidomain Method for Compressible Flows. II. A Semi-Structured Method , 1996 .

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

[13]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[14]  N. Ron-Ho,et al.  A Multiple-Grid Scheme for Solving the Euler Equations , 1982 .

[15]  Zhi J. Wang,et al.  High-Order Multidomain Spectral Difference Method for the Navier-Stokes Equations , 2006 .

[16]  Dimitri J. Mavriplis,et al.  Implicit Solution of the Unsteady Euler Equations for High-Order Accurate Discontinuous Galerkin Discretizations , 2006 .

[17]  A. Jameson,et al.  Lower-upper Symmetric-Gauss-Seidel method for the Euler and Navier-Stokes equations , 1988 .

[18]  George Em Karniadakis,et al.  Galerkin and discontinuous Galerkin spectral/hp methods , 1999 .

[19]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[20]  Francesco Bassi,et al.  Numerical Solution of the Euler Equations with a Multiorder Discontinuous Finite Element Method , 2003 .

[21]  Chris Lacor,et al.  Dispersion and dissipation properties of the 1D spectral volume method and application to a p-multigrid algorithm , 2007, J. Comput. Phys..

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

[23]  E. Turkel,et al.  Preconditioned methods for solving the incompressible low speed compressible equations , 1987 .

[24]  Chi-Wang Shu,et al.  The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV. The multidimensional case , 1990 .

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

[26]  Zhi J. Wang,et al.  Extension of the spectral volume method to high-order boundary representation , 2006 .

[27]  Michael A. Leschziner,et al.  Average-State Jacobians and Implicit Methods for Compressible Viscous and Turbulent Flows , 1997 .

[28]  Brian T. Helenbrook,et al.  Analysis of ``p''-Multigrid for Continuous and Discontinuous Finite Element Discretizations , 2003 .

[29]  Rainald Löhner,et al.  A Class of Matrix-Free Implicit Methods for Compressible Flows on Unstructured Grids , 2001 .

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

[31]  David A. Kopriva,et al.  A Staggered-Grid Multidomain Spectral Method for the Compressible Navier-Stokes Equations , 1998 .

[32]  Y. Maday,et al.  Spectral element multigrid. Part 2: Theoretical justification , 1988 .

[33]  Antony Jameson,et al.  How Many Steps are Required to Solve the Euler Equations of Steady, Compressible Flow: In Search of a Fast Solution Algorithm , 2001 .

[34]  Zhi J. Wang,et al.  The Spectral Difference Method for the 2D Euler Equations on Unstructured Grids , 2005 .

[35]  Marcel Vinokur,et al.  Spectral difference method for unstructured grids I: Basic formulation , 2006, J. Comput. Phys..

[36]  Georg May,et al.  A Spectral Dierence Method for the Euler and Navier-Stokes Equations on Unstructured Meshes , 2006 .