Convergence of defect-correction and multigrid iterations for inviscid flows

Convergence of multigrid and defect-correction iterations is comprehensively studied within different incompressible and compressible inviscid regimes on medium to high-density grids. Good smoothing properties of the defect-correction relaxation have been shown using both a modified Fourier analysis and a more general idealized-coarse-grid analysis. Single-grid defect correction alone has some slowly converging iterations on grids of medium density. The convergence is especially slow for near-sonic flows and for low Mach numbers. Additionally, the fast asymptotic convergence seen on medium density grids deteriorates on high-density grids. Certain downstream-boundary modes are slowly damped on high-density grids. Multigrid accelerates convergence of the slow defect-correction iterations to the extent determined by the coarse-grid correction. The two-level asymptotic convergence rates are stable and significantly below one in most of the regions but slow convergence is noted for near-sonic and low-Mach compressible flows. The multigrid solver has been applied to the NACA 0012 airfoil and to different flow regimes, such as near-tangency and stagnation. Certain convergence difficulties have been encountered within stagnat ion regions. Nonetheless, for the airfoil flow, with a sharp trailing-edge, residuals were fast converging for a subcritical flow on a sequence of grids. For supercritical flow, residuals converged slower on some intermedi ate grids than on the finest grid or the two coarsest grids. At either conditions, convergence of drag below the level of discretization errors occurs in a single cycle.

[1]  John C. Vassberg,et al.  In Pursuit of Grid Convergence for Two-Dimensional Euler Solutions , 2009 .

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

[3]  Boris Diskin,et al.  Towards Verification of Unstructured-Grid Solvers , 2008 .

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

[5]  R. C. Nelson,et al.  Experimental Investigation of the Confluent Boundary Layer of a High-Lift System , 2000 .

[6]  J. Edwards A low-diffusion flux-splitting scheme for Navier-Stokes calculations , 1997 .

[7]  W. K. Anderson,et al.  Comparison of Finite Volume Flux Vector Splittings for the Euler Equations , 1985 .

[8]  W. Mulder A high-resolution Euler solver based on multigrid, semi-coarsening, and defective correction , 1992 .

[9]  Barry Koren,et al.  Multigrid solution method for the steady RANS equations , 2007, J. Comput. Phys..

[10]  Boris Diskin,et al.  Towards Verification of Unstructured-Grid Solvers , 2008 .

[11]  Achi Brandt,et al.  Multigrid Techniques: 1984 Guide with Applications to Fluid Dynamics, Revised Edition , 2011 .

[12]  James L. Thomas,et al.  Half-Space Analysis of the Defect-Correction Method for Fromm Discretization of Convection , 2000, SIAM J. Sci. Comput..

[13]  J. Edwards,et al.  Low-Diffusion Flux-Splitting Methods for Flows at All Speeds , 1997 .

[14]  Jack R. Edwards,et al.  Development of O(Nm2) Preconditioned Multigrid Solvers for Euler and Navier-Stokes Equations , 2000 .

[15]  Achi Brandt,et al.  Textbook Multigrid Efficiency for Computational Fluid Dynamics Simulations , 2001 .

[16]  Wim A. Mulder,et al.  A Note on the Use of Symmetric Line Gauss-Seidel for the Steady Upwind Differenced Euler Equations , 1987, SIAM J. Sci. Comput..

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

[18]  W. K. Anderson,et al.  An implicit upwind algorithm for computing turbulent flows on unstructured grids , 1994 .

[19]  BORIS DISKIN,et al.  On Quantitative Analysis Methods for Multigrid Solutions , 2005, SIAM J. Sci. Comput..

[20]  W. K. Anderson,et al.  An O(Nm{2}) Plane Solver for the Compressible Navier-Stokes Equations , 1999 .