Communications on Applied Mathematics and Computation >

2024 , Vol. 6 >Issue 2: 811 - 836

DOI: https://doi.org/10.1007/s42967-023-00254-0

ORIGINAL PAPERS

Iterative Subregion Correction Preconditioners with Adaptive Tolerance for Problems with Geometrically Localized Stiffness

Expand
  • 1. Department of Mathematics, University of California, Berkeley, Berkeley, CA 94720-3840, USA;
    2. Fariborz Maseeh Department of Mathematics and Statistics, Portland State University, Portland, OR 97201, USA

Received date: 2022-06-23

  Revised date: 2022-11-21

  Accepted date: 2023-01-07

  Online published: 2023-03-21

Fold

Abstract

We present a class of preconditioners for the linear systems resulting from a finite element or discontinuous Galerkin discretizations of advection-dominated problems. These preconditioners are designed to treat the case of geometrically localized stiffness, where the convergence rates of iterative methods are degraded in a localized subregion of the mesh. Slower convergence may be caused by a number of factors, including the mesh size, anisotropy, highly variable coefficients, and more challenging physics. The approach taken in this work is to correct well-known preconditioners such as the block Jacobi and the block incomplete LU (ILU) with an adaptive inner subregion iteration. The goal of these preconditioners is to reduce the number of costly global iterations by accelerating the convergence in the stiff region by iterating on the less expensive reduced problem. The tolerance for the inner iteration is adaptively chosen to minimize subregion-local work while guaranteeing global convergence rates. We present analysis showing that the convergence of these preconditioners, even when combined with an adaptively selected tolerance, is independent of discretization parameters (e.g., the mesh size and diffusion coefficient) in the subregion. We demonstrate significant performance improvements over black-box preconditioners when applied to several model convection-diffusion problems. Finally, we present performance results of several variations of iterative subregion correction preconditioners applied to the Reynolds number 2.25×106 fluid flow over the NACA 0012 airfoil, as well as massively separated flow at 30° angle of attack.

Key words： Subregion correction; Nested Krylov; Geometrically localized stiffness

Cite this article

Michael Franco, Per-Olof Persson, Will Pazner . Iterative Subregion Correction Preconditioners with Adaptive Tolerance for Problems with Geometrically Localized Stiffness[J]. Communications on Applied Mathematics and Computation, 2024 , 6(2) : 811 -836 . DOI: 10.1007/s42967-023-00254-0

References

[1] Alexander, R.: Diagonally implicit Runge-Kutta methods for stiff ODE’s. SIAM J. Numer. Anal. 14(6), 1006-1021 (1977). https://doi.org/10.1137/0714068
[2] Anderson, R., et al.: MFEM: a modular finite element methods library. Comput. Math. Appl. 81, 42-74 (2021). https://doi.org/10.1016/j.camwa.2020.06.009
[3] Antonietti, P.F., Houston, P.: A class of domain decomposition preconditioners for hp-discontinuous Galerkin finite element methods. J. Sci. Comput. 46(1), 124-149 (2011). https://doi.org/10.1007/s10915-010-9390-1
[4] Berger, M.J., Colella, P.: Local adaptive mesh refinement for shock hydrodynamics. J. Comput. Phys. 82(1), 64-84 (1989). https://doi.org/10.1016/0021-9991(89)90035-1
[5] Briggs, W.L., Henson, V.E., McCormick, S.F.: A Multigrid Tutorial. SIAM, Philadelphia (2000). https://doi.org/10.1137/1.9780898719505
[6] Cai, X.-C.: Some domain decomposition algorithms for nonselfadjoint elliptic and parabolic partial differential equations. PhD thesis, New York University, New York (1989)
[7] Cockburn, B., Shu, C.-W.: The local discontinuous Galerkin method for time-dependent convection-diffusion systems. SIAM J. Numer. Anal. 35(6), 2440-2463 (1998). https://doi.org/10.1137/s0036142997316712
[8] Constantinescu, E.M., Sandu, A.: Multirate timestepping methods for hyperbolic conservation laws. J. Sci. Comput. 33(3), 239-278 (2007). https://doi.org/10.1007/s10915-007-9151-y
[9] D’Azevedo, E.F., Forsyth, P.A., Tang, W.-P.: Ordering methods for preconditioned conjugate gradient methods applied to unstructured grid problems. SIAM J. Matrix Anal. Appl. 13(3), 944-961 (1992). https://doi.org/10.1137/0613057
[10] Desjardins, B., Grenier, E., Lions, P.-L., Masmoudi, N.: Incompressible limit for solutions of the isentropic Navier-Stokes equations with Dirichlet boundary conditions. Journal de Mathématiques Pures et Appliquées. 78(5), 461-471 (1999). https://doi.org/10.1016/s0021-7824(99)00032-x
[11] Falgout, R.D., Yang, U.M.: hypre: A library of high performance preconditioners. In: International Conference on Computational Science, pp. 632-641. Springer (2002). https://doi.org/10.1007/3-540-47789-6_66.
[12] Fernandez, P., Nguyen, N.-C., Peraire, J.: On the ability of discontinuous Galerkin methods to simulate under-resolved turbulent flows. Preprint at https://doi.org/10.48550/ARXIV.1810.09435 (2018)
[13] Golub, G.H., Ye, Q.: Inexact preconditioned conjugate gradient method with inner-outer iteration. SIAM J. Sci. Comput. 21(4), 1305-1320 (1999). https://doi.org/10.1137/s1064827597323415
[14] Griebel, M., Oswald, P.: On the abstract theory of additive and multiplicative Schwarz algorithms. Numerische Mathematik 70(2), 163-180 (1995). https://doi.org/10.1007/s002110050115
[15] Heil, M., Hazel, A.L., Boyle, J.: Solvers for large-displacement fluid-structure interaction problems: segregated versus monolithic approaches. Comput. Mech. 43(1), 91-101 (2008). https://doi.org/10.1007/s00466-008-0270-6
[16] Huang, D.Z., Persson, P.-O., Zahr, M.J.: High-order, linearly stable, partitioned solvers for general multiphysics problems based on implicit-explicit Runge-Kutta schemes. Comput. Methods Appl. Mech. Eng. 346, 674-706 (2019). https://doi.org/10.1016/j.cma.2018.09.015
[17] Jourdan de Araujo Jorge Filho, E., Wang, Z.J.: A study of p-multigrid approach for the high order FR/CPR method. In: AIAA Aviation 2019 Forum (2019). https://doi.org/10.2514/6.2019-3711
[18] Kanevsky, A., Carpenter, M.H., Gottlieb, D., Hesthaven, J.S.: Application of implicit-explicit high order Runge-Kutta methods to discontinuous-Galerkin schemes. J. Comput. Phys. 225(2), 1753-1781 (2007). https://doi.org/10.1016/j.jcp.2007.02.021
[19] Karabelas, S.J.: Large eddy simulation of high-Reynolds number flow past a rotating cylinder. Int. J. Heat Fluid Flow 31(4), 518-527 (2010). https://doi.org/10.1016/j.ijheatfluidflow.2010.02.010
[20] Knupp, P.M.: Algebraic mesh quality metrics. SIAM J. Sci. Comput. 23(1), 193-218 (2001). https://doi.org/10.1137/S1064827500371499
[21] Krivodonova, L.: An efficient local time-stepping scheme for solution of nonlinear conservation laws. J. Comput. Phys. 229(22), 8537-8551 (2010). https://doi.org/10.1016/j.jcp.2010.07.037
[22] Nguyen, C., Terrana, S., Peraire, J.: Wall-resolved implicit large eddy simulation of transonic buffet over the OAT15A airfoil using a discontinuous Galerkin method. In: AIAA Scitech 2020 Forum (2020). https://doi.org/10.2514/6.2020-2062
[23] Pazner, W., Franco, M., Persson, P.-O.: High-order wall-resolved large eddy simulation of transonic buffet on the OAT15A airfoil. In: AIAA Scitech 2019 Forum (2019). https://doi.org/10.2514/6.2019-1152
[24] Pazner, W., Persson, P.-O.: Stage-parallel fully implicit Runge-Kutta solvers for discontinuous Galerkin fluid simulations. J. Comput. Phys. 335, 700-717 (2017). https://doi.org/10.1016/j.jcp.2017.01.050
[25] Pazner, W., Persson, P.-O.: Approximate tensor-product preconditioners for very high order discontinuous Galerkin methods. J. Comput. Phys. 354, 344-369 (2018). https://doi.org/10.1016/j.jcp.2017.10.030
[26] Peraire, J., Persson, P.-O.: The compact discontinuous Galerkin (CDG) method for elliptic problems. SIAM J. Sci. Comput. 30(4), 1806-1824 (2008). https://doi.org/10.1137/070685518
[27] Persson, P.-O.: High-order LES simulations using implicit-explicit Runge-Kutta schemes. In: 49th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition (2011). https://doi.org/10.2514/6.2011-684
[28] Persson, P.-O., Peraire, J.: Newton-GMRES preconditioning for discontinuous Galerkin discretizations of the Navier-Stokes equations. SIAM J. Sci. Comput. 30(6), 2709-2733 (2008). https://doi.org/10.1137/070692108
[29] Roe, P.L.: Approximate Riemann solvers, parameter vectors, and difference schemes. J. Comput. Phys. 43(2), 357-372 (1981). https://doi.org/10.1016/0021-9991(81)90128-5
[30] Saad, Y.: A flexible inner-outer preconditioned GMRES algorithm. SIAM J. Sci. Comput. 14(2), 461-469 (1993). https://doi.org/10.1137/0914028
[31] Saad, Y.: Iterative Methods for Sparse Linear Systems. Society for Industrial and Applied Mathematics, Philadelphia (2003). https://doi.org/10.1137/1.9780898718003
[32] Saad, Y., Zhang, J.: BILUM: block versions of multielimination and multilevel ILU preconditioner for general sparse linear systems. SIAM J. Sci. Comput. 20(6), 2103-2121 (1999). https://doi.org/10.1137/S106482759732753X
[33] Sato, M., Okada, K., Nonomura, T., Aono, H., Yakeno, A., Asada, K., Abe, Y., Fujii, K.: Massive parametric study by LES on separated-flow control around airfoil using DBD plasma actuator at Reynolds number 63,000. In: 43rd AIAA Fluid Dynamics Conference (2013). https://doi.org/10.2514/6.2013-2750
[34] Sundar, H., Biros, G., Burstedde, C., Rudi, J., Ghattas, O., Stadler, G.: Parallel geometric-algebraic multigrid on unstructured forests of octrees. In: SC ’12: Proceedings of the International Conference on High Performance Computing, Networking, Storage and Analysis, pp. 1-11. IEEE (2012). https://doi.org/10.1109/SC.2012.91
[35] Toselli, A., Widlund, O.: Domain Decomposition Methods—Algorithms and Theory, vol. 34. Springer, Heidelberg (2005). https://doi.org/10.1007/b137868
[36] Uranga, A., Persson, P.-O., Drela, M., Peraire, J.: Implicit large eddy simulation of transition to turbulence at low Reynolds numbers using a discontinuous Galerkin method. Int. J. Numer. Methods Eng. 87(1/2/3/4/5), 232-261 (2011). https://doi.org/10.1002/nme.3036
[37] Van Den Eshof, J., Sleijpen, G.L., van Gijzen, M.B.: Relaxation strategies for nested Krylov methods. J. Comput. Appl. Math. 177(2), 347-365 (2005). https://doi.org/10.1016/j.cam.2004.09.024
[38] Wanner, G., Hairer, E.: Solving Ordinary Differential Equations II, vol. 375. Springer, Heidelberg (1996). https://doi.org/10.1007/978-3-642-05221-7
[39] Yang, U.M., et al.: BoomerAMG: a parallel algebraic multigrid solver and preconditioner. Appl. Numer. Math. 41(1), 155-177 (2002). https://doi.org/10.1016/S0168-9274(01)00115-5
Options
Outlines

/

〈 〉

Copyright © Editorial By Communications on Applied Mathematics and Computation
沪ICP备09014157