Rocky Mountain Journal of Mathematics

A geometric full multigrid method and fourth-order compact scheme for the 3D Helmholtz equation on nonuniform grid discretization

Zhenwei Yang, Xiaobin Li, Lei Feng, and Xinxin Hou

Full-text: Access denied (no subscription detected)

We're sorry, but we are unable to provide you with the full text of this article because we are not able to identify you as a subscriber. If you have a personal subscription to this journal, then please login. If you are already logged in, then you may need to update your profile to register your subscription. Read more about accessing full-text


A full multigrid method and fourth-order compact difference scheme is designed to solve the 3D Helmholtz equation on unequal mesh size. Three dimensional restriction and prolongation operators of the multigrid method on unequal grids could be constructed based on volume law. Two numerical experiments are implemented, and the results show the computational efficiency and accuracy of the full multigrid method. The study also illustrates that the full multigrid method with fourth-order compact difference scheme has great advantages in computation, which has been time consuming, and in iterative convergence efficiency.

Article information

Rocky Mountain J. Math., Volume 49, Number 3 (2019), 1029-1047.

First available in Project Euclid: 23 July 2019

Permanent link to this document

Digital Object Identifier

Mathematical Reviews number (MathSciNet)

Zentralblatt MATH identifier

Primary: 65F10: Iterative methods for linear systems [See also 65N22]
Secondary: 65N06: Finite difference methods 65N55: Multigrid methods; domain decomposition

Full multigrid method Helmholtz equation restriction operator nonuniform grid


Yang, Zhenwei; Li, Xiaobin; Feng, Lei; Hou, Xinxin. A geometric full multigrid method and fourth-order compact scheme for the 3D Helmholtz equation on nonuniform grid discretization. Rocky Mountain J. Math. 49 (2019), no. 3, 1029--1047. doi:10.1216/RMJ-2019-49-3-1029.

Export citation


  • W.L. Briggs, V.E. Henson and S.F. McCormick, A multigrid tutorial, 2nd ed., SIAM, Philadelphia, 2000.
  • W. Cecot and M. Oleksy, High order FEM for multigrid homogenization, Comp. Math. Appl. 70 (2015), 1391–1400.
  • Z.Y. Chen, D.S. Cheng and W. Feng, T.T. Wu and H.Q. Yang, A multigrid-based preconditioned Krylov subspace method for the Helmholtz equation with PML, J. Math. Anal. Appl. 383 (2011), 522–540.
  • M. Commer and G.A. Newman, An accelerated time domain finite difference simulation scheme for three-dimensional transient electromagnetic modeling using geometricmultigrid concepts, Radio Sci. 41 (2016), 349–353.
  • M. Engel and M. Griebel, A multigrid method for constrained optimal control problems, J. Comp. Appl. Math. 235 (2011), 4368–4388.
  • Y.B. Ge, F.J. Cao and J. Zhang, A transformation-free HOC scheme and multigrid method for solving the 3D Poisson equation on nonuniform grids, J. Comp. Phys. 234 (2013), 199–216.
  • F. Ghaffar, N. Badshah, S. Islam and M.A. Khan, Multigrid method based on transformation-free high-order scheme for solving 2D Helmholtz equation on nonuniform grids, Adv. Diff. Eqs. 2016 (2016), 1–16.
  • B. Gmeiner, H. Köstler, M. Stürmer and U. Rüde, Parallel multigrid on hierarchical hybrid grids: a performance study on current high performance computing clusters, Concur. Comp. Pract. Exp. 26 (2013), 217–240.
  • X. Jia, J.G. Sun and H.H. Xie, A multigrid method for Helmholtz transmission Eigenvalue problems, J. Sci. Comp. 60 (2014), 276–294.
  • H. Knibbe, C.W. Oosterlee and C. Vuik, GPU implementation of a Helmholtz-Krylov solver preconditioned by a shifted Laplace multigrid method, J. Comp. Appl. Math. 236 (2011), 281–293.
  • I. Livshits, Multiple Galerkin adaptive algebraic multigrid algorithm for the Helmholtz equations, SIAM J. Sci. Comp. 37 (2015), S195–S215.
  • O.V. Nechaev, E.P. Shurina and M.A. Botchev, Multilevel iterative solvers for the edge finite element solution of the 3D Maxwell equation, Comp. Math. Appl. 55 (2008), 2346–2362.
  • V. Puzyrev, J. Puente and G. Houzeaux, Algebraic multigrid preconditioning within parallel finite-element solvers for 3D electromagnetic modelling problems in geophysics, Geophys. J. Int. 197 (2014), 1290–1298.
  • J. Shin, S. Kim, D. Lee and J. Kim, A parallel multigrid method of the Cahn-Hilliard equation, Comp. Mater. Sci. 71 (2013), 89–96.
  • R.N. Slaybaugh, T.M. Evans, G.G. Davidson and P.P.H. Wilson, Multigrid in energy preconditioner for Krylov solvers, J. Comp. Phys. 242 (2013), 405–419.
  • S.C. Stefano and T.P. Cristina, A note on algebraic multigrid methods for the discrete weighted Laplacian, Comp. Math. Appl. 60 (2010), 1290–1298.
  • C.C. Stolk, M. Ahmed and S.K. Bhowmik, A multigrid method for the Helmholtz equation with optimized coarse grid corrections, SIAM J. Sci. Comp. 36 (2014), A2819–A2841.
  • J.T. Tang, F.Y. Wang, Z.Y. Ren and R.W. Guo, 3D direct current resistivity forward modeling by adaptive multigrid finite element method, J. Central South Univ. 17 (2010), 587–592.
  • X. Wang, Z. Yang, G. Huang and B. Chen, A high-order compact difference scheme for 2D Laplace and Poisson equations in non-uniform grid systems, Comm. Nonlin. Sci. Numer. Simul. 14 (2009), 379–398.
  • Y. Wang and J. Zhang, Fast and robust sixth-order multigrid computation for the three-dimensional convection-diffusion equation, J. Comp. Appl. Math. 234 (2010), 3496–3506.
  • J. Wang, W.J. Zhong and J. Zhang, A general meshsize fourth-order compact difference discretization scheme for 3D Poisson equation, Appl. Math. Comp. 183 (2006), 804–812.