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

*(English)*Zbl 1349.76353Summary: 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.

##### MSC:

76M12 | Finite volume methods applied to problems in fluid mechanics |

65M08 | Finite volume methods for initial value and initial-boundary value problems involving PDEs |

65M55 | Multigrid methods; domain decomposition for initial value and initial-boundary value problems involving PDEs |

76G25 | General aerodynamics and subsonic flows |

##### Keywords:

Reynolds averaged Navier-Stokes equations; high Reynolds number; agglomeration multigrid; implicit smoother; aerodynamic applications; eigenvalue analysis
Full Text:
DOI

##### References:

[1] | Pierce, Niles A.; Giles, Michael B., Preconditioned multigrid methods for compressible flow calculations on stretched meshes, J. Comput. Phys., 136, 425-445, (1997) · Zbl 0893.76061 |

[2] | Rossow, Cord-Christian, Efficient computation of compressible and incompressible flows, J. Comput. Phys., 220, 879-899, (2007) · Zbl 1235.76096 |

[3] | Swanson, R. C.; Turkel, E.; Rossow, C.-C., Convergence acceleration of Runge-Kutta schemes for solving the Navier-Stokes equations, J. Comput. Phys., 224, 365-388, (2007) · Zbl 1261.76036 |

[4] | Swanson, R. C.; Rossow, C.-C., An efficient solver for the RANS equation and a one-equation turbulence model, Comput. Fluids, 42, 13-25, (2011) · Zbl 1271.76196 |

[5] | Baldwin, B. S.; Lomax, H., Thin layer approximation and algebraic model for separated turbulent flows, AIAA J., 78, 257, (1978) |

[6] | Cagnone, J. S.; Sermeus, Kurt; Nadarajah, Siva K.; Laurendeau, Eric, Implicit multigrid schemes for challenging aerodynamic simulations, Comput. Fluids, 44, 314-327, (2011) · Zbl 1271.76179 |

[7] | B. van Leer, Characteristic time-stepping or local preconditioning of the Euler equations, AIAA Paper-91-1552-CP, 1991. |

[8] | Langer, Stefan, Investigation and application of point implicit Runge-Kutta methods to inviscid flow problems, Int. J. Numer. Methods Fluids, 69, 2, 332-352, (2012) · Zbl 1245.76049 |

[9] | Langer, Stefan, Hierarchy of preconditioning techniques for the solution of the Navier-Stokes equations discretized by 2nd order unstructured finite volume methods, (ECCOMAS 2012 Conference, Conference Proceeding Series, (September 2012)) |

[10] | Pierce, Niles A.; Giles, Michael B.; Jameson, A.; Martinelli, L., Accelerating three-dimensional Navier-Stokes calculations, (13th AIAA Computational Fluid Dynamics Conference, Conference Proceeding Series, vol. 1997-1953, (1997), AIAA) |

[11] | Mavriplis, Dimitri J., Directional coarsening and smoothing for anisotropic Navier-Stokes problems, Electron. Trans. Numer. Anal., 6, 182-197, (1997) · Zbl 0901.76068 |

[12] | Mavriplis, Dimitri J., Directional agglomeration multigrid techniques for high-Reynolds number viscous flows, (1998), ICASE Report No. 98-7 |

[13] | Eliasson, Peter; Weinerfelt, Per; Nordström, Jan, Application of a line-implicit scheme on stretched unstructured grids, (Proc. 47th AIAA Aerospace Sciences Meeting, Conference Proceeding Series, vol. 2009-0163, (2009), AIAA) |

[14] | Langer, Stefan, Application of a line implicit method to fully coupled system of equations for turbulent flow problems, Int. J. Comput. Fluid Dyn., 27, 3, 131-150, (March 2013) |

[15] | A. Jameson, W. Schmidt, E. Turkel, Numerical solutions of the Euler equations by finite volume methods using Runge-Kutta time-stepping schemes, AIAA Paper, AIAA-81-1259, 1981. |

[16] | Moinier, Pierre, Algorithm developments for an unstructured viscous flow solver, (1999), University of Oxford, PhD thesis |

[17] | Moinier, Pierre; Müller, Jens-Dominik; Giles, Michael B., Edge-based multigrid and preconditioning for hybrid grids, AIAA J., 99, 3339, (1999) |

[18] | Nishikawa, Hiroaki; Diskin, Boris; Thomas, J. L., Development and application of agglomerated multigrid methods for complex geometries, (40th AIAA Fluid Dynamics Conference and Exhibit, (2010)), AIAA Paper 2010-4731, Chicago |

[19] | Spalart, P. R.; Allmaras, S. R., A one-equation turbulence model for aerodynamic flows, (AIAA Computational Fluid Dynamics Conference, Conference Proceeding Series, vol. 1992-439, (1992), AIAA) |

[20] | Allmaras, Steven R.; Johnson, Forrester T.; Spalart, Philippe R., Modifications and clarifications for the implementation of the Spalart-Allmaras turbulence model, (International Conference on Computational Fluid Dynamics 7, Hawaii, Conference Proceeding Series, vol. ICCFD7-1902, (2012)) |

[21] | Rumsey, C. L., Apparent transition behavior of widely-used turbulence models, Int. J. Heat Fluid Flow, 28, 1460-1471, (2007) |

[22] | Turkel, E., Improving the accuracy of central difference schemes, (Dwoyer, D. L.; Hussaini, M. Y., 11th International Conference on Numerical Methods in Fluid Dynamics, Lecture Notes in Physics, vol. 323, (1989), Springer-Verlag) |

[23] | Swanson, R. C.; Turkel, E., On central difference and upwind schemes, J. Comput. Phys., 101, 292-306, (1992) · Zbl 0757.76044 |

[24] | Svärd, Magnus; Nordström, Jan, Stability of finite volume approximations for the Laplacian operator om quadrilateral and triangular grids, Appl. Numer. Math., 51, 101-125, (2004) · Zbl 1065.65111 |

[25] | Magnus Svärd; Gong, Jing; Nordström, Jan, Stable artificial dissipation operators for finite volume schemes on unstructured grids, Appl. Numer. Math., 56, 1481-1490, (2006) · Zbl 1103.65096 |

[26] | Eliasson, Peter; Weinerfelt, Per; Nordström, Jan, The influence of weak and strong solid wall boundary conditions on the convergence to steady-state of the Navier-Stokes equations, (Proc. 19th AIAA Computational Fluid Dynamics, Conference Proceeding Series, vol. 2009-3551, (2009), AIAA) |

[27] | Hairer, Ernst; Norsett, Syvert P.; Wanner, Gerhard, Solving ordinary differential equations II. stiff and differential-algebraic problems, (1992), Springer New York · Zbl 0638.65058 |

[28] | Langer, Stefan; Li, Dian, Application of point implicit Runge-Kutta methods to inviscid and laminar flow problems using AUSM and AUSM^+ upwinding, Int. J. Comput. Fluid Dyn., 25, 5, (2011) · Zbl 1271.76187 |

[29] | Mavriplis, Dimitri J., Multigrid techniques for unstructured meshes, (1995), ICASE Report 95-27 |

[30] | Trottenberg, Ulrich; Oosterlee, C. W.; Schuller, Anton, Multigrid, (2001), Academic Press, Inc. Orlando, FL, USA · Zbl 0976.65106 |

[31] | Moulitsas, Irene; Karypis, George, Multilevel algorithms for generating coarse grids for multigrid methods, (Johnson, Greg, SC, (2001), ACM), 45 |

[32] | Lassaline, J. V.; Zingg, D. W., Development of an agglomeration multigrid algorithm with directional coarsening, (14th AIAA Computational Fluid Dynamics Conference, Conference Proceeding Series, vol. 99-3338, (1999)) |

[33] | T. Gerhold, M. Galle, O. Friedrich, J. Evans, Calculation of complex three-dimensional configurations employing the DLR-Tau-code, AIAA Paper 1997-167, 1997. |

[34] | Nishikawa, H.; Diskin, B.; Thomas, J. L., Critical study of agglomerated multigrid methods for diffusion, AIAA J., 48, 4, 839-847, (2010) |

[35] | Thomas, James L.; Diskin, Boris; Nishikawa, Hiroaki, A critical study of agglomerated multigrid methods for diffusion on highly-stretched grids, Special Issue Dedicated to Drs. W.R. Briley and H. McDonald, Comput. Fluids, 41, 1, 82-93, (2011) · Zbl 1271.76198 |

[36] | Cook, P. H.; McDonald, M. A.; Firmin, M. C.P., Aerofoil RAE 2822 pressure distributions and boundary layer and wake measurements, (1979), AGARD-AR, 138 |

[37] | Klausmeyer, S. M.; Lin, J. C., Comparative results from a CFD challenge over a 2D three-element high-lift airfoil, (1997), NASA Technical Memorandum 112858 |

[38] | Martineau, D. G.; Stokes, S.; Munday, S. J.; Jackson, A. P., Anisotropic hybrid mesh generation for industrial RANS applications, (44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, Conference Proceeding Series, vol. 2006-534, (2006)) |

[39] | Vassberg, John C., A unified baseline grid about the common research model wing-body for the fifth AIAA CFD drag prediction workshop, (Proc. 29th AIAA Applied Aerodynamics Conference, Honolulu, HI, Conference Proceeding Series, vol. 2011-3508, (June 2011), AIAA) |

[40] | Levy, David W.; Laflin, Kelly R.; Tinoco, Edward N.; Vassberg, John C.; Mani, Mori; Rider, Ben; Rumsey, Christopher L.; Wahls, Richard A.; Morrison, Joseph H.; Brodersen, Olaf P.; Crippa, Simone; Mavriplis, Dimitri J.; Murayama, Mitsuhiro, Summary of data from the fifth AIAA drag prediction workshop, (Proc. 51st AIAA Aerospace Sciences Meeting, Gravepine (Dallas/Ft. Worth Region), Texas, January 2013, Conference Proceeding Series, vol. 2013-0046, (2013), AIAA) |

[41] | Nishikawa, Hiroaki; Diskin, Boris; Thomas, James L.; Hammond, D. H., Recent advances in agglomerated multigrid, (51st AIAA Aerospace Sciences Meeting, (2013)), AIAA Paper 2013-863, Texas |

[43] | Eriksson, Lars E.; Rizzi, Arthur, Computer-aided analysis of the convergence to steady state of discrete approximations to the Euler equations, J. Comput. Phys., 57, 1, 90-128, (1985) · Zbl 0624.76078 |

[44] | Saad, Y., Variations on Arnoldi’s method for computing eigenelements of large unsymmetric matrices, Linear Algebra Appl., 34, 0, 269-295, (1980) · Zbl 0456.65017 |

[45] | Saad, Y., Iterative methods for sparse linear systems, (1996), International Thomson Publishing Boston · Zbl 1002.65042 |

[46] | Kuijlaars, A. B., Which eigenvalues are found by the Lanczos method?, SIAM J. Matrix Anal. Appl., 22, 306-321, (2000) · Zbl 0969.65028 |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.