Volume 9, Issue 2
A High-Accuracy Mechanical Quadrature Method for Solving the Axisymmetric Poisson’s Equation

Hu Li & Jin Huang

Adv. Appl. Math. Mech., 9 (2017), pp. 393-406.

Published online: 2018-05

Preview Full PDF 2 479
Export citation
  • Abstract

In this article, we consider the numerical solution for Poisson’s equation in axisymmetric geometry. When the boundary condition and source term are axisymmetric, the problem reduces to solving Poisson’s equation in cylindrical coordinates in the two-dimensional (r,z) region of the original three-dimensional domain S. Hence, the original boundary value problem is reduced to a two-dimensional one. To make use of the Mechanical quadrature method (MQM), it is necessary to calculate a particular solution, which can be subtracted off, so that MQM can be used to solve the resulting Laplace problem, which possesses high accuracy order O(h 3 max) and low computing complexities. Moreover, the multivariate asymptotic error expansion of MQM accompanied with O(h 3 i ) for all mesh widths hi is got. Hence, once discrete equations with coarse meshes are solved in parallel, the higher accuracy order of numerical approximations can be at least O(h 5 max) by the splitting extrapolation algorithm (SEA). Meanwhile, a posteriori asymptotic error estimate is derived, which can be used to construct self-adaptive algorithms. The numerical examples support our theoretical analysis.

  • Keywords

Mechanical quadrature method, splitting extrapolation algorithm, Poisson’s equation, a posteriori estimate.

  • AMS Subject Headings

65N38, 65R20

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • References
  • Hide All
    View All

@Article{AAMM-9-393, author = {}, title = {A High-Accuracy Mechanical Quadrature Method for Solving the Axisymmetric Poisson’s Equation}, journal = {Advances in Applied Mathematics and Mechanics}, year = {2018}, volume = {9}, number = {2}, pages = {393--406}, abstract = {

In this article, we consider the numerical solution for Poisson’s equation in axisymmetric geometry. When the boundary condition and source term are axisymmetric, the problem reduces to solving Poisson’s equation in cylindrical coordinates in the two-dimensional (r,z) region of the original three-dimensional domain S. Hence, the original boundary value problem is reduced to a two-dimensional one. To make use of the Mechanical quadrature method (MQM), it is necessary to calculate a particular solution, which can be subtracted off, so that MQM can be used to solve the resulting Laplace problem, which possesses high accuracy order O(h 3 max) and low computing complexities. Moreover, the multivariate asymptotic error expansion of MQM accompanied with O(h 3 i ) for all mesh widths hi is got. Hence, once discrete equations with coarse meshes are solved in parallel, the higher accuracy order of numerical approximations can be at least O(h 5 max) by the splitting extrapolation algorithm (SEA). Meanwhile, a posteriori asymptotic error estimate is derived, which can be used to construct self-adaptive algorithms. The numerical examples support our theoretical analysis.

}, issn = {2075-1354}, doi = {https://doi.org/10.4208/aamm.2015.m1287}, url = {http://global-sci.org/intro/article_detail/aamm/12155.html} }
Copy to clipboard
The citation has been copied to your clipboard