Agglomeration multigrid methods with implicit Runge-Kutta smoothers applied to aerodynamic simulations on unstructured grids

For unstructured finite volume methods an agglomeration multigrid with an implicit multistage Runge-Kutta method as a smoother is developed for solving the compressible Reynolds averaged Navier-Stokes (RANS) equations. The implicit Runge-Kutta method is interpreted as a preconditioned explicit Runge-Kutta method. The construction of the preconditioner is based on an approximate derivative. The linear systems are solved approximately with a symmetric Gauss-Seidel method. To significantly improve this solution method grid anisotropy is treated within the Gauss-Seidel iteration in such a way that the strong couplings in the linear system are resolved by tridiagonal systems constructed along these directions of strong coupling. The agglomeration strategy is adapted to this procedure by taking into account exactly these anisotropies in such a way that a directional coarsening is applied along these directions of strong coupling. Turbulence effects are included by a Spalart-Allmaras model, and the additional transport-type equation is approximately solved in a loosely coupled manner with the same method. For two-dimensional and three-dimensional numerical examples and a variety of differently generated meshes we show the wide range of applicability of the solution method. Finally, we exploit the GMRES method to determine approximate spectral information of the linearized RANS equations. This approximate spectral information is used to discuss and compare characteristics of multistage Runge-Kutta methods.

[1]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[2]  H. Lomax,et al.  Thin-layer approximation and algebraic model for separated turbulent flows , 1978 .

[3]  Stefan Langer,et al.  Application of a line implicit method to fully coupled system of equations for turbulent flow problems , 2013 .

[4]  Magnus Svärd,et al.  Stable artificial dissipation operators for finite volume schemes on unstructured grids , 2006 .

[5]  Arno B. J. Kuijlaars,et al.  Which Eigenvalues Are Found by the Lanczos Method? , 2000, SIAM J. Matrix Anal. Appl..

[6]  Jan Nordström,et al.  The Influence of Weak and Strong Solid Wall Boundary Conditions on the Convergence to Steady-State of the Navier-Stokes Equations , 2009 .

[7]  Boris Diskin,et al.  A critical study of agglomerated multigrid methods for diffusion on highly-stretched grids , 2011 .

[8]  Antony Jameson,et al.  Accelerating three-dimensional navier-stokes calculations , 1997 .

[9]  David W. Zingg,et al.  Development of an agglomeration multigrid algorithm with directional coarsening , 1999 .

[10]  C. Rumsey Apparent Transition Behavior of Widely-Used Turbulence Models , 2006 .

[11]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[12]  R. C. Swanson,et al.  An efficient solver for the RANS equations and a one-equation turbulence model , 2011 .

[13]  I. Moulitsas,et al.  Multilevel Algorithms for Generating Coarse Grids for Multigrid Methods , 2001, ACM/IEEE SC 2001 Conference (SC'01).

[14]  Eli Turkel,et al.  On Central-Difference and Upwind Schemes , 1992 .

[15]  Magnus Svärd,et al.  Stability of finite volume approximations for the Laplacian operator on quadrilateral and triangular grids , 2004 .

[16]  T. Gerhold,et al.  Calculation of Complex Three-Dimensional Configurations Employing the DLR-tau-Code , 1997 .

[17]  John C. Vassberg,et al.  A Unified Baseline Grid about the Common Research Model Wing-Body for the Fifth AIAA CFD Drag Prediction Workshop , 2011 .

[18]  Cord-Christian Rossow,et al.  Efficient computation of compressible and incompressible flows , 2007, J. Comput. Phys..

[19]  Pierre Moinier,et al.  Edge-Based Multigrid and Preconditioning for Hybrid Grids , 1999 .

[20]  Philip L. Roe,et al.  Characteristic time-stepping or local preconditioning of the Euler equations , 1991 .

[21]  Jan Nordström,et al.  Application of a line-implicit scheme on stretched unstructured grids , 2009 .

[22]  D. Mavriplis Directional coarsening and smoothing for anisotropic Navier-Stokes problems. , 1997 .

[23]  Michael B. Giles,et al.  Preconditioned Multigrid Methods for Compressible Flow Calculations on Stretched Meshes , 1997 .

[24]  James L. Thomas,et al.  Recent Advances in Agglomerated Multigrid , 2013 .

[25]  Y. Saad Variations on Arnoldi's method for computing eigenelements of large unsymmetric matrices , 1980 .

[26]  A. Jameson,et al.  Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes , 1981 .

[27]  Stefan Langer,et al.  Investigation and application of point implicit Runge–Kutta methods to inviscid flow problems , 2012 .

[28]  Stefan Langer,et al.  Application of point implicit Runge–Kutta methods to inviscid and laminar flow problems using AUSM and AUSM+ upwinding , 2011 .

[29]  Dimitri J. Mavriplis,et al.  MULTIGRID TECHNIQUES FOR UNSTRUCTURED MESHES , 1995 .

[30]  E. Hairer,et al.  Solving Ordinary Differential Equations II , 2010 .

[31]  Steven M. Klausmeyer,et al.  Comparative Results From a CFD Challenge Over a 2D Three-Element High-Lift Airfoil , 1997 .

[32]  Andrew Jackson,et al.  Anisotropic Hybrid Mesh Generation for Industrial RANS Applications , 2006 .

[33]  Eli Turkel,et al.  Convergence acceleration of Runge-Kutta schemes for solving the Navier-Stokes equations , 2007, J. Comput. Phys..

[34]  Arthur Rizzi,et al.  Computer-aided analysis of the convergence to steady state of discrete approximations to the euler equations , 1985 .

[35]  Boris Diskin,et al.  Development and Application of Agglomerated Multigrid Methods for Complex Geometries , 2010 .

[36]  Eric Laurendeau,et al.  Implicit multigrid schemes for challenging aerodynamic simulations on block-structured grids , 2011 .

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

[38]  Eli Turkel Improving the accuracy of central difference schemes , 1989 .

[39]  Dimitri J. Mavriplis,et al.  Directional Agglomeration Multigrid Techniques for High-Reynolds Number Viscous Flows , 1998 .

[40]  Boris Diskin,et al.  A Critical Study of Agglomerated Multigrid Methods for Diffusion , 2009 .

[41]  Jens-Dominik Mueller,et al.  Edge-based multigrid and preconditioning for hybrid grids , 1999 .