Journal of Applied Mathematics

Journal of Applied Mathematics / 2014 / Article

Research Article | Open Access

Volume 2014 |Article ID 490540 | 6 pages | https://doi.org/10.1155/2014/490540

Richardson Cascadic Multigrid Method for 2D Poisson Equation Based on a Fourth Order Compact Scheme

Academic Editor: Song Cen
Received09 Sep 2013
Revised31 Dec 2013
Accepted02 Jan 2014
Published10 Mar 2014

Abstract

Based on a fourth order compact difference scheme, a Richardson cascadic multigrid (RCMG) method for 2D Poisson equation is proposed, in which the an initial value on the each grid level is given by the Richardson extrapolation technique (Wang and Zhang (2009)) and a cubic interpolation operator. The numerical experiments show that the new method is of higher accuracy and less computation time.

1. Introduction

Poisson equation is a partial differential equation (PDE) with broad applications in theoretical physics, mechanical engineering and other fields, such as groundwater flow [1, 2], fluid pressure prediction [3], electromagnetics [4], semiconductor modeling [5], and electrical power network modeling [6].

We consider the following two-dimensional (2D) Poisson equation: where is a rectangular domain or union of rectangular domains with Dirichlet boundary . The solution and the forcing function are assumed to be sufficiently smooth.

Multigrid (MG) method is one of the most effective algorithms to solve the large scale problem. In 1996, cascadic multigrid (CMG) method proposed by Bornemann and Deuflhard [7] and then analyzed by Shi et al. (see [811]) and Shaidurov (see [12]). In the recent years, there have been several theoretical analyses and the applications of these methods for the plate bending problems (see [13]), the parabolic problems (see [10]), the nonlinear problems (see [14, 15]), and the Stokes problems (see [16]). In order to improve the efficiency of the CMG, some new extrapolation formulas and extrapolation cascadic multigrid (EXCMG) methods are proposed by Chen et al. (see [1720]). These new methods can provide a better initial value for smoothing operator on the refined grid level to accelerate their convergence rate.

Based on the Richardson extrapolation technique, Wang and Zhang [21] presented a multiscale multigrid algorithm. Numerical experiments show that the new method is of higher accuracy solution and higher efficiency.

In this paper, in order to develop a more efficient CMG method, we use the Richardson extrapolation technique presented in [21] and a new extrapolation formula; a new Richardson extrapolation cascadic multigrid (RCMG) method for 2D Poisson equation is proposed.

The sections are arranged as follows: the fourth order compact difference scheme and Richardson extrapolation technique are given in Section 2. Chen’s new extrapolation formula and EXCMG method are introduced in Section 3. In Section 4, we present the RCMG method. In Section 5, the numerical experiments show the effectiveness of the new method.

2. Fourth Order Compact Difference Scheme and Richardson Extrapolation Technique

For convenience, we consider the rectangular domain . We discretize with uniform mesh sizes and in the and coordinate directions. The mesh points are with and , and , . Let's denote the mesh aspect ratio , and be the solution at the grid point , we can rewrite the fourth order compact difference scheme of (1) into the following form [22]: The coefficients in (2) are If the domain is subdivided into a sequence of grids (or ), with step length (namely, , by using the fourth order compact difference scheme (see (2)), a series of linear equations of the model problem (1) are given as follows

Assume the fourth order accurate solutions and on the grid and the grid are given, respectively (Figure 1). In 2009, Wang and Zhang [21] applied the Richardson extrapolation (where to get a sixth order accurate solution on .

The above extrapolation operator is rewritten as the following iterative operator .

Algorithm 1. Consider .
Step  1. Set , .
Step  2. Update every (even, even) grid point on by Richardson extrapolation formula (see (5)); then use direct interpolation to get . Consider Step  3. Update every (odd, odd) grid point on . From (2), for each (odd, odd) point , the updated solution is Here, represents the right-hand side part of (2).
Step  4. Update every (odd, even) grid point on . From (2), for each (odd, even) grid point, the updated value is
Step  5. Update every (even, odd) grid point on . From (2), the idea is similar to the (odd, even) grid point. Let .
Step  6. If or , stop. Else, let and return to Step  3.

3. New Extrapolation Formula and EXCMG Method

Based on an asymptotic expansion of finite element method, a new extrapolation formula and an extrapolation cascadic multigrid (EXCMG) method are proposed by Chen et al. (see [1720]). The numerical experiments show that the EXCMG method is of high accuracy and efficiency. Now we rewrite the new extrapolation formula as follows.

Let us denote the above new extrapolation formula by operator

Now let , on , denote the exact solutions, the EXCMG method is as following:

Algorithm 2 (EXCMG). For , consider the following
Step  1. Extrapolate by using the new extrapolation formula (see (10))
Step  2. Compute the initial value on by using quadratic interpolation operator .
Step  3. Smooth times to get the iterative solution on by using some classical iterative operator .
Step  4. Return to Step  1 if , until you get the final iterative solution on the finest grid .

4. Richardson Cascadic Multigrid Method

One of the main tasks in cascadic multigrid method is constructing a suitable interpolation. Based on a new extrapolation-interpolation formula, Chen [1720] proposed the following extrapolation cascadic multigrid (EXCMG) method, in which the new extrapolation and quadratic interpolation are used to provide a better initial value on refined grid.

In this section, we use RET operator and a cubic interpolation to interpolate the initial guess on the refined grid . Then a classical iterative operator (such as conjugate gradient method) is used as a smoothing operator to compute the high accuracy solution on the fine grid . Similar to the standard CMG method, we propose the following Richardson cascadic multigrid (RCMG) method.

Algorithm  3 [RCMG]

Step  1. Exactly solve the equation on coarsest grid .

Step  2. Run Algorithm 1; we have

Step  3. Use a cubic interpolation operator to have the initial value on the gird level .

Step  4. Smoothing times by using the classical iterative operator , on the level . Set ;

Step  5. Return to Step  2, if .

The difference between RCMG method and EXCMG method is that

5. Numerical Experiment and Comparison

Numerical experiments are conducted to solve a 2D Poisson equation (1) on the unit square domain .

Example 4. The exact solution ; the forcing function

Example 5. The exact solution ; the forcing function
We use the conjugate gradient (CG) method as a smoothing iterative operator in EXCMG method and RCMG method. In EXCMG method, the number of iterations on each grid level has to increase from finer to coarser grids; in this paper let . And in RCMG, we set the number of iteration (Step  2) and (Step  4) be . We set of RET in the RCMG method (on Step  2).

5.1. Comparison of the Initial Errors

Assume that the exact solutions of the difference equation on grids and are given. We compare EXCMG method with RCMG method for the initial error on grid .

From Figure 2, the accuracy of the initial error on the next grid of RCMG method is higher than EXCMG method. Namely, a better initial value on the fine grid can be got by using RCMG method. Based on the results of the literature [1720], the RCMG method can obtain good convergence rate.

5.2. Comparison between EXCMG Method and RCMG Method

Let denote the maximum absolute error between the computed solution and the exact solution on the finest grid points. The “” denotes the computing time (unit: second) of EXCMG method and RCMG method.

From Figures 3 and 4 and Tables 1 and 2, we see that, under the same conditions, the RCMG method can obtain higher computational precision and spend less computing time than EXCMG method.


EXCMG RCMG
cpu cpu

3
128 0.37 0.14
256 1.12 0.31
512 4.10 1.12

4
128 0.34 0.08
256 1.19 0.42
512 4.38 1.29

5
128 0.51 0.08
256 1.06 0.31
512 4.15 1.15


EXCMG RCMG
cpu cpu

3
128 0.39 0.16
256 1.22 0.50
512 4.17 1.23

4
128 0.25 0.11
256 1.00 0.30
512 4.23 1.17

5
128 0.28 0.09
256 1.01 0.31
512 4.07 1.06

6. Conclusion

In this paper, based on a fourth order compact scheme, we present a Richardson cascadic multigrid method for 2D Poisson problem by using Richardson technique presented by [21]. The numerical results show that RCMG method has higher computational accuracy and higher efficiency.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grant no. 11161014), the National Natural Science Foundation of Yunnan Province (Grant no. 2012FD054), and Scientific Research Starting Foundation for Master or Ph.D. of Honghe University (Grant no. XJ1S0925).

References

  1. R. E. Alcouffe, A. Brandt, J. E. Dendy et al., “The multigrid method or the diffusion equation with strongly discontinuous coefficients,” SIAM Journal on Scientific and Statistical Computing, no. 2, pp. 430–454, 1981. View at: Google Scholar | MathSciNet
  2. C. E. Kees, C. T. Miller, E. W. Jenkins, and C. T. Kelley, “Versatile two-level Schwarz preconditioners for multiphase flow,” Computational Geosciences, vol. 7, no. 2, pp. 91–114, 2003. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  3. C. Vuik, A. Segal, and J. A. Meijerink, “An efficient preconditioned CG method for the solution of a class of layered problems with extreme contrasts in the coefficients,” Journal of Computational Physics, vol. 152, no. 1, pp. 385–403, 1999. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  4. B. Heise and M. Kuhn, “Parallel solvers for linear and nonlinear exterior magnetic field problems based upon coupled FE/BE formulations,” Computing, vol. 56, no. 3, pp. 237–258, 1996. View at: Google Scholar | Zentralblatt MATH | MathSciNet
  5. R. K. Coomer and I. G. Graham, “Massively parallel methods for semiconductor device modelling,” Computing, vol. 56, no. 1, pp. 1–27, 1996. View at: Google Scholar | Zentralblatt MATH | MathSciNet
  6. V. E. Howle and S. A. Vavasis, “An iterative method for solving complex-symmetric systems arising in electrical power modeling,” SIAM Journal on Matrix Analysis and Applications, vol. 26, no. 4, pp. 1150–1178, 2005. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  7. F. A. Bornemann and P. Deuflhard, “The cascadic multigrid method for elliptic problems,” Numerische Mathematik, vol. 75, no. 2, pp. 135–152, 1996. View at: Google Scholar | Zentralblatt MATH | MathSciNet
  8. Z.-C. Shi, X.-J. Xu, and Y.-Q. Huang, “Economical cascadic multigrid method (ECMG),” Science in China A, vol. 50, no. 12, pp. 1765–1780, 2007. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  9. Z.-C. Shi and X. Xu, “Cascadic multigrid method for elliptic problems,” East-West Journal of Numerical Mathematics, vol. 7, no. 3, pp. 199–209, 1999. View at: Google Scholar | Zentralblatt MATH | MathSciNet
  10. Z.-C. Shi and X.-J. Xu, “Cascadic multigrid for parabolic problems,” Journal of Computational Mathematics, vol. 18, no. 5, pp. 551–560, 2000. View at: Google Scholar
  11. Z. Shi and X. Xu, “A new cascadic multigrid,” Science in China A, vol. 44, no. 1, pp. 21–30, 2001. View at: Google Scholar
  12. V. V. Shaidurov, “Some estimates of the rate of convergence for the cascadic conjugate-gradient method,” Computers and Mathematics with Applications, vol. 31, no. 4-5, pp. 161–171, 1996. View at: Google Scholar
  13. R. Stevenson, “Nonconforming finite elements and the cascadic multi-grid method,” Numerische Mathematik, vol. 91, no. 2, pp. 351–387, 2002. View at: Publisher Site | Google Scholar
  14. Y. Huang, Z. Si, T. Tang, and W. Xue, “A multilevel successive iteration method for nonlinear elliptic problems,” Mathematics of Computation, vol. 73, no. 246, pp. 525–539, 2004. View at: Publisher Site | Google Scholar
  15. G. Timmermann, “A cascadic multigrid algorithm for semilinear elliptic problems,” Numerische Mathematik, vol. 86, no. 4, pp. 717–731, 2000. View at: Google Scholar
  16. D. Braess and W. Dahmen, “A cascadic multigrid algorithm for the Stokes equations,” Numerische Mathematik, vol. 82, no. 2, pp. 179–191, 1999. View at: Google Scholar
  17. C. Chen, Z.-C. Shi, and H. Hu, “On extrapolation cascadic multigrid method,” Journal of Computational Mathematics, vol. 29, no. 6, pp. 684–697, 2011. View at: Publisher Site | Google Scholar
  18. C. M. Chen, Z. Q. Xie, C. L. Li et al., “Study of a new extrapolation multigrid method,” Journal of Natural Science of Hunan Normal University, vol. 30, no. 2, pp. 1–5, 2007. View at: Google Scholar
  19. C. M. Chen, H. L. Hu, Z. Q. Xie, and C. Li, “Analysis of extrapolation cascadic multigrid method (EXCMG),” Science in China A, vol. 51, no. 8, pp. 1349–1360, 2008. View at: Publisher Site | Google Scholar
  20. C. M. Chen, H. L. Hu, Z. Q. Xie et al., “L2-error of extrapolation cascadic multigrid (EXCMG),” Acta Mathematica Scientia, vol. 29, no. 3, pp. 539–551, 2009. View at: Publisher Site | Google Scholar
  21. Y. Wang and J. Zhang, “Sixth order compact scheme combined with multigrid method and extrapolation technique for 2D poisson equation,” Journal of Computational Physics, vol. 228, no. 1, pp. 137–146, 2009. View at: Publisher Site | Google Scholar
  22. J. Zhang, “Multigrid method and fourth-order compact scheme for 2D Poisson equation with unequal mesh-size discretization,” Journal of Computational Physics, vol. 179, no. 1, pp. 170–179, 2002. View at: Publisher Site | Google Scholar

Copyright © 2014 Li Ming and Li Chen-Liang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

838 Views | 482 Downloads | 1 Citation
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.