DOCSLIB.ORG

Multigrid Method for Solving 2D-Poisson Equation with Sixth

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Multigrid Method for Solving 2D-Poisson Equation with Sixth International Mathematical Forum, Vol. 6, 2011, no. 55, 2729 - 2736 Multigrid Method for Solving 2D-Poisson Equation with Sixth Order Finite Difference Method Bouthina S. Ahmed1 and S. J. Monaquel2 1) Mathematics Department, Rabigh Faculty of Science & Arts King Abdul Aziz University,P.O. Box, 344, Rabigh 21911, Saudi Arabia [email protected] 2) Mathematics Department, Faculty of Science King Abdul Aziz University, Jeddah, Saudi Arabia [email protected] Abstract In this paper we develop a sixth order finite difference discretization strategy to solve the two dimensional Poisson equation. We use multigrid V-cycle procedure to built multiscale multigrid method which is similar to the full multigrid method. Numerical results are given to illustrate this method. Mathematics subject classification: 65Mxx Keywords: Poisson equation, six order finite difference method, multigrid method. 1-Introduction Poisson equation is a partial differential equation (PDF) with broad application s in mechanical engineering, theoretical physics and other fields. The two dimensional (2D) Poisson equation can be written in the form: uxx (x, y) + u yy (x, y) = f (x, y), (x, y) ∈ Ω (1) Where Ω is a rectangular domain, or union of rectangular domains, with suitable boundary conditions defined on ∂Ω . The solution u(x,y) and forcing function f(x,y) are 2730 B. S. Ahmed and S. J. Monaquel assumed to be sufficiently smooth and have the necessary derivatives up to certain orders. A second order accurate solution can be computed by applying standard second order centeral difference operators uxx (x, y) and u yy (x, y) in Eq (1). Higher order (more than two) accurate discertization methods need more complex procedure than the second order accurate discretization method to compute the coefficient matrix, but they usually generate linear systems of much smaller size, compared with that from the lower order accurate discetization method [1,4] . There has been growing interest in developing higher order accurate discretization methods, especially the high order compact difference scheme to solve partial differential equations (PDEs) [6, 8, 9, 10]. 2-Sixth order compact approximation This method is similar to a sixth-order accurate approximation to the derivatives calculated from Helmholtz equation [5]. We the scheme for the two dimensional uniform Cartesian grids with grid spacing Δx = Δy = h. The mesh points are (xi , y j ) with xi = ih a y j = jh, 0 ≤ i ≤ N, 0 ≤ j ≤ N where N is the number of uniform intervals in the x and y directions. We write the second order central difference operators as u − 2u + u u − 2u + u 2 i+1, j i, j i−1, j 2 i, j+1 i, j i, j−1 δ u = , δ u = , x i, j 2 y i, j 2 h h (2) Using Taylor series expansions at the grid point (xi , y j ) , we have 2 2 4 6 6 2 ∂ u h ∂ u h ∂ u 6 δ u = + + + O(h ) x i, j 2 4 6 ∂x 12 ∂x 360 ∂x (3) 2 2 4 6 6 2 ∂ u h ∂ u h ∂ u 6 δ u = + + + O(h ) y i, j 2 4 6 ∂y 12 ∂y 360 ∂y (4) The central difference difference for Eq. (1) can be written as 2 2 δ x ui, j + δ y ui, j + Ti, j = fi, j (5) where ui, j = u (xi , y j ) , fi, j = f (xi , y j ) and Multigrid method for solving 2D-Poisson equation 2731 2 4 4 4 6 6 h ⎛ ∂ u ∂ u ⎞ h ⎛ ∂ u ∂ u ⎞ 6 T = − ⎜ + ⎟ − ⎜ + ⎟ + O(h ) i, j ⎜ 4 4 ⎟ ⎜ 6 6 ⎟ 12 ⎝ ∂x ∂y ⎠ 360 ⎝ ∂x ∂y ⎠ (6) Using the following appropriate derivatives of Eq. (1) 4 2 2 4 4 2 2 4 ∂ u ∂ f 2 ∂ u ∂ u ∂ u ∂ f 2 ∂ u ∂ u = − k − , = − k − , 4 2 2 2 2 4 2 2 2 2 ∂x ∂x ∂x ∂x ∂ y ∂y ∂y ∂y ∂y ∂x (7) In Eq. (6) we get 2 ⎛ ⎡ 4 ⎤ ⎞ 4 ⎛ 6 6 ⎞ h ⎜ 2 ∂ u ⎟ h ⎜ ∂ u ∂ u ⎟ 6 Ti, j = − ∇ fi, j − 2⎢ ⎥ − + + O(h ) 12 ⎜ 2 2 ⎟ 360 ⎜ 6 6 ⎟ ⎝ ⎣⎢∂x ∂y ⎦⎥i, j ⎠ ⎝ ∂x ∂y ⎠i, j (8) ∂4u The fourth order approximation of in Eq. (8) can be written as ∂x2∂2 y 4 2 6 6 ⎡ ∂ u ⎤ 2 2 h ⎡ ∂ u ∂ u ⎤ 4 ⎢ ⎥ = δ x δ yui, j − ⎢ + ⎥ + O(h ) 2 2 12 4 2 2 4 ⎣⎢∂x ∂ y ⎦⎥i, j ⎣⎢∂x ∂ y ∂x ∂ y ⎦⎥i, j (9) Substituting Eq. (9) into Eq. (8), we get 2 4 6 6 6 6 h 2 2 2 h ∂ u ∂ u ∂ u ∂ u 6 Ti, j = (∇ fi, j + 2δ x δ yui, j ) 6 + 5 4 2 + 5 2 4 + 6 + O(h ) 12 360 x x y x y y ∂ ∂ ∂ ∂ ∂ ∂ i, j (10) Clearly, getting a compact sixth-order approximation requires compact expressions of the four derivatives of order six in Eq. (1) which can be done by further differentiating Eq. (10) that is 4 6 6 ∂ f ∂ u ∂ u = + 2 2 4 2 2 4 ∂x ∂y ∂x ∂y ∂x ∂y (11) 2732 B. S. Ahmed and S. J. Monaquel 6 6 6 6 ∂ u ∂ u 4 ⎛ ∂ u ∂ u ⎞ + = ∇ f − ⎜ + ⎟ 6 6 ⎜ 4 2 2 4 ⎟ ∂x ∂y ⎝ ∂x ∂y ∂x ∂y ⎠ (12) Substituting Eq. (7), (11) into Eq. (13) gives 6 6 4 ∂ u ∂ u 4 ∂ f + = ∇ f − 6 6 2 2 ∂x ∂y ∂x ∂y (13) Using Eq. (11), (13) we can eliminate all derivatives of u in Eq. (10) that is 2 4 ⎛ 4 ⎞ h 2 2 2 h ⎜ 4 ⎡ ∂ u ⎤ ⎟ 6 T = − ∇ f + 2δ δ u − ∇ f + 4 + O(h ) i, j ()i, j x y i, j ⎜ i, j ⎢ 2 2 ⎥ ⎟ 12 360 ⎜ ⎢∂x ∂ y ⎥ ⎟ ⎝ ⎣ ⎦i, j ⎠i, j (14) The compact sixth-order approximation of two dimensional can be written as 2 2 4 4 ⎡ 4 ⎤ h 2 2 h 2 h 4 h ⎢ ∂ f ⎥ δ δ u + ⎜⎛δ 2 + δ 2 ⎟⎞u = f + ∇ f + ∇ f + x y i, j x y i, j i, j i, j i, j ⎢ ⎥ 6 ⎝ ⎠ 12 360 90 ∂2∂ 2 ⎣⎢ x y ⎦⎥ (15) We can express Eq. (16) in the form 10 2 − u + ⎜⎛u + u + u + u ⎟⎞ 3 i, j 3 ⎝ i + 1, j i − 1, j i, j + 1 i, j − 1⎠ 1 + ⎜⎛u + u + u + u ⎟⎞ 6 ⎝ i + 1, j + 1 i + 1, j − 1 i − 1, j + 1 i − 1, j − 1⎠ 4 6 6 ⎡ 4 ⎤ 2 h 2 h 4 h ∂ f = h f + ∇ f + ∇ f + ⎢ ⎥ (16) i, j i, j i, j 12 360 90 ⎢∂2∂2 ⎥ ⎣⎢ x y ⎦⎥ 3- Multigrid method Algorithm 1 u 1- Let 2h be the solution on the coarse gridΩ 2h Multigrid method for solving 2D-Poisson equation 2733 2- Use some high order interpolation schemes here we use Newton difference h interpolation, to interpolate Ω2h , uh = I2h to the coarse grid (we interpolate even, even, odd, even and even, odd grids points Fig. 1). 3- Update every (odd-odd) grid point on From Eq. (16) for each point the updated solution is ⎡ 2 ⎤ F + ⎜⎛u + u + u + u ⎟⎞ ⎢ i, j i + 1, j i − 1, j i, j + 1 i, j − 1 ⎥ k + 1 3 ⎝ ⎠ ⎛10 ⎞ u = ⎢ ⎥ /⎜ ⎟ i, j 1 ⎛ ⎞ ⎝ 3 ⎠ ⎢+ ⎜u k + u k + u k + u k ⎟⎥ ⎣⎢ 6 ⎝ i + 1, j + 1 i − 1, j + 1 i + 1, j − 1 i − 1, j − 1 ⎠⎦⎥ Here, Fi, j represents the right-hand side of Eq. (16) 4- Update every (odd, even) grid point on Ωh from Eq. (16) 5- Update every (even, odd) grid point on Ωh from Eq. (16) h,k+1 h,k 6- Compute the l norm Ω R = u − u if not converged go back to 2 2h 2 step3. 7- Compute residual on the fine grid Ωh from Lh − fh = dh and use full weighting scheme to project residual from fine grid to the coarse grid. 8- Use interpolation to transfer corrections from the coarse grid to the fine grid. L u = f 9- Relax υ1 times on h h h . 10- Use from the previous step as the initial guess to run the multigrid algorithm until it Converges. 2734 B. S. Ahmed and S. J. Monaquel Fig. 1. Illustration of the operator based interpolation scheme for a 5x5 fine grid 4- Modified multiscale multigrid method The convergence rate of multigrid method is independent of the grid size [2, 3, 7]. Here we use a geometric multiscale multigrid method [3] which is similar to full multigrid method. Algorithm 2- Multiscale multigrid method h 1- Use the Newton interpolation difference method to interpolateuh = I 2hu2h . 2- Relax υ1 times on Lhuh = f h . 2h R r 3- Restrict the residual r to the coarse grid Ω2h by h h (we use full weighting scheme to restrict the residual from the fine grid to the coarse grid) L u = f 4- Relax υ 2 times on h h h . 5- Use uh from the previous step as the initial guess to run the multigrid V-cycle algorithm MG(uh , fh ) on the fine grid until it converges. Multigrid method for solving 2D-Poisson equation 2735 5- Numerical Results Problem A Consider Poisson equation 2 2 ∂ u ∂ u 2 + = −2π sin(πx)sin(πy) 0 ≤ x ≤ 1 and 0 ≤ y ≤ 1 ∂x2 ∂y2 with Dirichlet boundary conditions on all sides f unit square, that is u(0, y) = u(1, y) = u(x,0) = u(x,1) = 0 The exact solution is u(x, y) = sin πx sin πy Problem B Consider Poisson equation ∂2u ∂2u + = −sin π y, 0 ≤ x ≤ 1 and 0 ≤ y ≤ 1 ∂x2 ∂y2 with the Dirichlet boundary conditions u(x,0) = u(x,1) = 0 1 u((0, y) = 0, u(1, y) = − sin(πy)(eπ −1) π 2 1 The exact solution is u(x, y) = − sin(π y)(eπx −1) π In Table 1, 2 we compute the l2 - norm, maximum norm of the error and the order of accuracy, for different values of N for problem A and B respectively.
Load more

Recommended publications
  • Multigrid Methods
    Multigrid Methods Volker John Winter Semester 2013/14 Contents 1 Literature 2 2 Model Problems 3 3 Detailed Investigation of Classical Iterative Schemes 8 3.1 General Aspects of Classical Iterative Schemes . .8 3.2 The Jacobi and Damped Jacobi Method . 10 3.3 The Gauss{Seidel Method and the SOR Method . 15 3.4 Summary . 18 4 Grid Transfer 19 4.1 Algorithms with Coarse Grid Systems, the Residual Equation . 19 4.2 Prolongation or Interpolation . 21 4.3 Restriction . 23 5 The Two-Level Method 27 5.1 The Coarse Grid Problem . 27 5.2 General Approach for Proving the Convergence of the Two-Level Method . 30 5.3 The Smoothing Property of the Damped Jacobi Iteration . 32 5.4 The Approximation Property . 33 5.5 Summary . 36 6 The Multigrid Method 37 6.1 Multigrid Cycles . 37 6.2 Convergence of the W-cycle . 39 6.3 Computational Work of the Multigrid γ-Cycle . 45 7 Algebraic Multigrid Methods 49 7.1 Components of an AMG and Definitions . 49 7.2 Algebraic Smoothness . 52 7.3 Coarsening . 56 7.4 Prolongation . 58 7.5 Concluding Remarks . 60 8 Outlook 61 1 Chapter 1 Literature Remark 1.1 Literature. There are several text books about multigrid methods, e.g., • Briggs et al. (2000), easy to read introduction, • Hackbusch (1985), the classical book, sometimes rather hard to read, • Shaidurov (1995), • Wesseling (1992), an introductionary book, • Trottenberg et al. (2001). 2 2 Chapter 2 Model Problems Remark 2.1 Motivation. The basic ideas and properties of multigrid methods will be explained in this course on two model problems.
    [Show full text]
  • Hybrid Particle-Grid Water Simulation Using Multigrid Pressure Solver Per Karlsson
    LiU-ITN-TEK-G--14/006-SE Hybrid Particle-Grid Water Simulation using Multigrid Pressure Solver Per Karlsson 2014-03-13 Department of Science and Technology Institutionen för teknik och naturvetenskap Linköping University Linköpings universitet nedewS ,gnipökrroN 47 106-ES 47 ,gnipökrroN nedewS 106 47 gnipökrroN LiU-ITN-TEK-G--14/006-SE Hybrid Particle-Grid Water Simulation using Multigrid Pressure Solver Examensarbete utfört i Medieteknik vid Tekniska högskolan vid Linköpings universitet Per Karlsson Handledare George Baravdish Examinator Camilla Forsell Norrköping 2014-03-13 Upphovsrätt Detta dokument hålls tillgängligt på Internet – eller dess framtida ersättare – under en längre tid från publiceringsdatum under förutsättning att inga extra- ordinära omständigheter uppstår. Tillgång till dokumentet innebär tillstånd för var och en att läsa, ladda ner, skriva ut enstaka kopior för enskilt bruk och att använda det oförändrat för ickekommersiell forskning och för undervisning. Överföring av upphovsrätten vid en senare tidpunkt kan inte upphäva detta tillstånd. All annan användning av dokumentet kräver upphovsmannens medgivande. För att garantera äktheten, säkerheten och tillgängligheten finns det lösningar av teknisk och administrativ art. Upphovsmannens ideella rätt innefattar rätt att bli nämnd som upphovsman i den omfattning som god sed kräver vid användning av dokumentet på ovan beskrivna sätt samt skydd mot att dokumentet ändras eller presenteras i sådan form eller i sådant sammanhang som är kränkande för upphovsmannens litterära eller konstnärliga anseende eller egenart. För ytterligare information om Linköping University Electronic Press se förlagets hemsida http://www.ep.liu.se/ Copyright The publishers will keep this document online on the Internet - or its possible replacement - for a considerable time from the date of publication barring exceptional circumstances.
    [Show full text]
  • Multigrid Algorithms for High Order Discontinuous Galerkin Methods
    Multigrid algorithms for high order discontinuous Galerkin methods Paola F. Antonietti1, Marco Sarti2, and Marco Verani3 Abstract I n this paper we study the performance of h- and p-multigrid algorithms for high order Discontinuous Galerkin discretizations of elliptic problems. We test the performance of the multigrid schemes employing a wide class of smoothers and considering both two- and three-dimensional test cases. 1 Introduction In the framework of multigrid solvers for Discontinuous Galerkin (DG) schemes, the first contributions are due to Gopalakrishnan and Kanschat [16] and Brenner and Zhao [10]. In [16] a V-cycle preconditioner for a Sym- metric Interior Penalty (SIP) discretization of an elliptic problem is analyzed. They prove that the condition number of the preconditioned system is uni- formly bounded with respect to the mesh size and the number of levels. In [10] V-cycle, F-cycle and W-cycle multigrid schemes for SIP discretizations are presented and analyzed, employing the additive theory developed in [8, 9]. A uniform bound for the error propagation operator is shown provided the number of smoothing steps is large enough. All the previously cited works focus on low order, i.e., linear, DG approximations. With regard to high order DG discretizations, h- and p-multigrid schemes are successfully em- ployed for the numerical solution of many different kinds of problems, see e.g. [14, 20, 22, 21, 24, 6], even if only few theoretical results are available that show the convergence properties of the underlying algorithms. In the MOX, Dipartimento di Matematica, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milano ITALY.
    [Show full text]
  • Hybrid Multigrid/Schwarz Algorithms for the Spectral Element Method
    Hybrid Multigrid/Schwarz Algorithms for the Spectral Element Method James W. Lottes¤ and Paul F. Fischery February 4, 2004 Abstract We study the performance of the multigrid method applied to spectral element (SE) discretizations of the Poisson and Helmholtz equations. Smoothers based on finite element (FE) discretizations, overlapping Schwarz methods, and point-Jacobi are con- sidered in conjunction with conjugate gradient and GMRES acceleration techniques. It is found that Schwarz methods based on restrictions of the originating SE matrices converge faster than FE-based methods and that weighting the Schwarz matrices by the inverse of the diagonal counting matrix is essential to effective Schwarz smoothing. Sev- eral of the methods considered achieve convergence rates comparable to those attained by classic multigrid on regular grids. 1 Introduction The availability of fast elliptic solvers is essential to many areas of scientific computing. For unstructured discretizations in three dimensions, iterative solvers are generally optimal from both work and storage standpoints. Ideally, one would like to have computational complexity that scales as O(n) for an n-point grid problem in lRd, implying that the it- eration count should be bounded as the mesh is refined. Modern iterative methods such as multigrid and Schwarz-based domain decomposition achieve bounded iteration counts through the introduction of multiple representations of the solution (or the residual) that allow efficient elimination of the error at each scale. The theory for these methods is well established for classical finite difference (FD) and finite element (FE) discretizations, and order-independent convergence rates are often attained in practice. For spectral element (SE) methods, there has been significant work on the development of Schwarz-based methods that employ a combination of local subdomain solves and sparse global solves to precondition conjugate gradient iteration.
    [Show full text]
  • Multigrid and Saddle-Point Preconditioners for Unfitted Finite
    Multigrid and saddle-point preconditioners for unfitted finite element modelling of inclusions Hardik Kothari∗and Rolf Krause† Institute of Computational Science, Universit`adella Svizzera italiana, Lugano, Switzerland February 17, 2021 Abstract In this work, we consider the modeling of inclusions in the material using an unfitted finite element method. In the unfitted methods, structured background meshes are used and only the underlying finite element space is modified to incorporate the discontinuities, such as inclusions. Hence, the unfitted methods provide a more flexible framework for modeling the materials with multiple inclusions. We employ the method of Lagrange multipliers for enforcing the interface conditions between the inclusions and matrix, this gives rise to the linear system of equations of saddle point type. We utilize the Uzawa method for solving the saddle point system and propose preconditioning strategies for primal and dual systems. For the dual systems, we review and compare the preconditioning strategies that are developed for FETI and SIMPLE methods. While for the primal system, we employ a tailored multigrid method specifically developed for the unfitted meshes. Lastly, the comparison between the proposed preconditioners is made through several numerical experiments. Keywords: Unfitted finite element method, multigrid method, saddle-point problem 1 Introduction In the modeling of many real-world engineering problems, we encounter material discontinuities, such as inclusions. The inclusions are found naturally in the materials, or can be artificially introduced to produce desired mechanical behavior. The finite element (FE) modeling of such inclusions requires to generate meshes that can resolve the interface between the matrix and inclusions, which can be a computationally cumbersome and expensive task.
    [Show full text]
  • Multigrid Method Based on a Space-Time Approach with Standard Coarsening for Parabolic Problems
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by CWI's Institutional Repository Applied Mathematics and Computation 317 (2018) 25–34 Contents lists available at ScienceDirect Applied Mathematics and Computation journal homepage: www.elsevier.com/locate/amc Multigrid method based on a space-time approach with standard coarsening for parabolic problems Sebastião Romero Franco a,b, Francisco José Gaspar c, Marcio Augusto Villela ∗ Pinto d, Carmen Rodrigo e, a Department of Mathematics, State University of Centro-Oeste, Irati, PR, Brazil b Graduate Program in Numerical Methods in Engineering, Federal University of Paraná, Curitiba, PR, Brazil c Centrum Wiskunde & Informatica (CWI), Science Park 123, P.O. Box 94079, Amsterdam 1090, The Netherlands d Department of Mechanical Engineering, Federal University of Paraná, Curitiba, PR, Brazil e IUMA and Applied Mathematics Department, University of Zaragoza, María de Luna, 3, Zaragoza 50018, Spain a r t i c l e i n f o a b s t r a c t MSC: In this work, a space-time multigrid method which uses standard coarsening in both tem- 00-01 poral and spatial domains and combines the use of different smoothers is proposed for the 99-00 solution of the heat equation in one and two space dimensions. In particular, an adaptive Keywords: smoothing strategy, based on the degree of anisotropy of the discrete operator on each Space-time multigrid grid-level, is the basis of the proposed multigrid algorithm. Local Fourier analysis is used Local Fourier analysis for the selection of the crucial parameter defining such an adaptive smoothing approach.
    [Show full text]
  • Electrostatic PIC with Adaptive Cartesian Mesh
    Electrostatic PIC with adaptive Cartesian mesh Vladimir Kolobov1,2 and Robert Arslanbekov1 1 CFD Research Corporation, Huntsville, AL, USA 2 The University of Alabama in Huntsville, Huntsville, AL 35899, USA E-mail: [email protected] Abstract. We describe an initial implementation of an electrostatic Particle-in-Cell (ES-PIC) module with adaptive Cartesian mesh in our Unified Flow Solver framework. Challenges of PIC method with cell-based adaptive mesh refinement (AMR) are related to a decrease of the particle-per-cell number in the refined cells with a corresponding increase of the numerical noise. The developed ES-PIC solver is validated for capacitively coupled plasma, its AMR capabilities are demonstrated for simulations of streamer development during high-pressure gas breakdown. It is shown that cell-based AMR provides a convenient particle management algorithm for exponential multiplications of electrons and ions in the ionization events. 1. Introduction Particle-in-Cell (PIC) is a mature technology widely used for plasma simulations. There are three versions of PIC: an electrostatic (EC) PIC uses Poisson equation for calculating electric field, an electromagnetic (EM) PIC employs full-wave Maxwell solver for calculations of electric and magnetic fields, and a quasi-static (QS) PIC uses a quasi-static version of the Maxwell equations (such as Darwin model) for calculation of the EM fields. Surprisingly, PIC solvers with adaptive mesh refinement (AMR) are relatively rare. It is worth mentioning the WARP code for accelerator applications [1], and an explicit EM-PIC code of Fujimoto [2] for astrophysics applications, both using block-structured AMR. In this paper, we describe an initial implementation of an ES-PIC module in our Unified Flow Solver (UFS) framework [3].
    [Show full text]
  • Hp-Multigrid As Smoother Algorithm for Higher Order Discontinuous Galerkin Discretizations of Advection Dominated Flows Part I
    HP-MULTIGRID AS SMOOTHER ALGORITHM FOR HIGHER ORDER DISCONTINUOUS GALERKIN DISCRETIZATIONS OF ADVECTION DOMINATED FLOWS PART I. MULTILEVEL ANALYSIS J.J.W. VAN DER VEGT∗ AND S. RHEBERGENy AMS subject classifications. Primary 65M55; Secondary 65M60, 76M10. Abstract. The hp-Multigrid as Smoother algorithm (hp-MGS) for the solution of higher order accurate space-(time) discontinuous Galerkin discretizations of advection dominated flows is pre- sented. This algorithm combines p-multigrid with h-multigrid at all p-levels, where the h-multigrid acts as smoother in the p-multigrid. The performance of the hp-MGS algorithm is further improved using semi-coarsening in combination with a new semi-implicit Runge-Kutta method as smoother. A detailed multilevel analysis of the hp-MGS algorithm is presented to obtain more insight into the theoretical performance of the algorithm. As model problem a fourth order accurate space-time dis- continuous Galerkin discretization of the advection-diffusion equation is considered. The multilevel analysis shows that the hp-MGS algorithm has excellent convergence rates, both for low and high cell Reynolds numbers and on highly stretched meshes. Key words. multigrid algorithms, discontinuous Galerkin methods, higher order accurate dis- cretizations, space-time methods, Runge-Kutta methods, Fourier analysis, multilevel analysis 1. Introduction. Discontinuous Galerkin finite element methods are well suited to obtain higher order accurate discretizations on unstructured meshes. The use of basis functions which are only weakly coupled to neighboring elements results in a local discretization which allows, in combination with hp-mesh adaptation, the efficient capturing of detailed structures in the solution, and is also beneficial for parallel computing.
    [Show full text]
  • CONVERGENCE of ALGEBRAIC MULTIGRID BASED on SMOOTHED AGGREGATION∗ ˇ PETR VANEK† , MARIAN BREZINA‡ , and JAN MANDEL§ February 1998 Revised February 2000 Abstract
    CONVERGENCE OF ALGEBRAIC MULTIGRID BASED ON SMOOTHED AGGREGATION∗ ˇ PETR VANEKy , MARIAN BREZINAz , AND JAN MANDELx February 1998 Revised February 2000 Abstract. We prove an abstract convergence estimate for the Algebraic Multigrid Method with prolongator defined by a disaggregation followed by a smoothing. The method input is the problem matrix and a matrix of the zero energy modes of the same problem but with natural boundary conditions. The construction is described in the case of a general elliptic system. The condition number bound increases only as a polynomial of the number of levels, and requires only a uniform weak approximation property for the aggregation operators. This property can be a-priori verified computationally once the aggregates are known. For illustration, it is also verified here for a uniformly elliptic diffusion equations discretized by linear conforming quasiuniform finite elements. Only very weak and natural assumptions on the hierarchy of aggregates are needed. Key words. Algebraic multigrid, zero energy modes, convergence theory, computational mechanics, finite elements, iterative solvers. 1. Introduction. This paper establishes abstract convergence bounds for an Algebraic Multigrid Method (AMG) based on smoothed aggregation. The bounds are obtained by invoking the general convergence theory of [6]. Our main result is a bound on the condition number that grows only as a power of the number of levels, and requires only a weak approximation property for the aggregates, similar to the weak approximation condition in classical AMG investigations [9, 22, 28, 37]. Our weak approximation condition can be easily verified computationally, and we show that it holds for general unstructured meshes and under natural assumptions on aggregates used to construct the coarse levels.
    [Show full text]
  • Multigrid Solvers for Immersed Finite Element Methods and Immersed Isogeometric Analysis
    Computational Mechanics https://doi.org/10.1007/s00466-019-01796-y ORIGINAL PAPER Multigrid solvers for immersed finite element methods and immersed isogeometric analysis F. de Prenter1,3 · C. V. Verhoosel1 · E. H. van Brummelen1 · J. A. Evans2 · C. Messe2,4 · J. Benzaken2,5 · K. Maute2 Received: 26 March 2019 / Accepted: 10 November 2019 © The Author(s) 2019 Abstract Ill-conditioning of the system matrix is a well-known complication in immersed finite element methods and trimmed isogeo- metric analysis. Elements with small intersections with the physical domain yield problematic eigenvalues in the system matrix, which generally degrades efficiency and robustness of iterative solvers. In this contribution we investigate the spectral prop- erties of immersed finite element systems treated by Schwarz-type methods, to establish the suitability of these as smoothers in a multigrid method. Based on this investigation we develop a geometric multigrid preconditioner for immersed finite element methods, which provides mesh-independent and cut-element-independent convergence rates. This preconditioning technique is applicable to higher-order discretizations, and enables solving large-scale immersed systems at a computational cost that scales linearly with the number of degrees of freedom. The performance of the preconditioner is demonstrated for conventional Lagrange basis functions and for isogeometric discretizations with both uniform B-splines and locally refined approximations based on truncated hierarchical B-splines. Keywords Immersed finite element method · Fictitious domain method · Iterative solver · Preconditioner · Multigrid 1 Introduction [14–21], scan based analysis [22–27] and topology optimiza- tion, e.g., [28–34]. Immersed methods are useful tools to avoid laborious and An essential aspect of finite element methods and iso- computationally expensive procedures for the generation of geometric analysis is the computation of the solution to a body-fitted finite element discretizations or analysis-suitable system of equations.
    [Show full text]
  • HMG--Homogeneous Multigrid For
    HMG — HOMOGENEOUS MULTIGRID FOR HDG PEIPEI LU, ANDREAS RUPP, AND GUIDO KANSCHAT Abstract. We introduce a homogeneous multigrid method in the sense that it uses the same HDG discretization scheme for Poisson’s equation on all levels. In particular, we construct a stable injection operator and prove optimal convergence of the method under the assumption of elliptic regularity. Numerical experiments underline our analytical findings. 1. Introduction While hybridizable discontinuous Galerkin (HDG) methods have be- come popular in mathematics and applications over the last years, lit- erature on efficient solution of the resulting discrete systems is still scarce. In this article, we propose to our knowledge the first multigrid preconditioner for such methods which is homogeneous in the sense that it uses the same discretization scheme on all levels. Such methods are important, since they have the same mathematical properties on all levels. They can also be advantageous from a computational point of view, since their data structures and execution patterns are more regular. HDG methods have been gaining popularity in the last decade. Orig- inally, they were analyzed for the Laplacian, see for instance [12] for an overview. Meanwhile, they have been applied to stationary [11, 14, 15, 17, 28, 29] and instationary [25] Stokes systems, to the locking-free discretization of problems in elasticity [19], as well as to plates [9, 24] and beams [7]. However, only few results have been achieved for solv- ing the large systems of equations arising from this type of numerical method, while multigrid [21] and domain decomposition [18] solvers arXiv:2011.14018v1 [math.NA] 27 Nov 2020 for earlier discontinuous Galerkin (DG) schemes have been available for many years.
    [Show full text]
  • An Efficient Finite Element Based Multigrid Method for Simulations Of
    An efficient finite element based multigrid method for simulations of the mechanical behavior of heterogeneous materials using CT images Xiaodong Liu, Julien Réthoré, Marie-Christine Baietto, Philippe Sainsot, Antonius Lubrecht To cite this version: Xiaodong Liu, Julien Réthoré, Marie-Christine Baietto, Philippe Sainsot, Antonius Lubrecht. An effi- cient finite element based multigrid method for simulations of the mechanical behavior of heterogeneous materials using CT images. Computational Mechanics, Springer Verlag, 2020, 10.1007/s00466-020- 01909-y. hal-02919095 HAL Id: hal-02919095 https://hal.archives-ouvertes.fr/hal-02919095 Submitted on 21 Aug 2020 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Distributed under a Creative Commons Attribution| 4.0 International License An efficient finite element based multigrid method for simulations of the mechanical behavior of heterogeneous materials using CT images Xiaodong Liu1, Julien R´ethor´e∗1, Marie-Christine Baietto2, Philippe Sainsot2, and Antonius Adrianus Lubrecht2 1GeM, CNRS UMR 6183, Centrale Nantes, F-44321, France 2Univ Lyon, INSA-Lyon, CNRS UMR5259, LaMCoS, F-69621, France 19 August, 2020 Abstract X-ray tomography techniques give researchers the full access to material inner structures. With such ample information, employing numerical simulation on real material images becomes more and more common.
    [Show full text]