Premium
On higher order multigrid methods with application to a geothermal reservoir model
Author(s) -
Gary John
Publication year - 1982
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/fld.1650020104
Subject(s) - multigrid method , mathematics , partial differential equation , extrapolation , dirichlet boundary condition , elliptic partial differential equation , iterative method , elliptic curve , separable space , convergence (economics) , mathematical optimization , boundary value problem , mathematical analysis , economics , economic growth
This paper considers the multigrid iterative method applied to the solution of finite difference approximations to a linear second‐order self‐adjoint elliptic equation. It represents an extension of work by Dinar and Brandt. We compare two methods to obtain fourth‐order convergence. The first is local error extrapolation developed by Brandt, the second is iterative improvement developed by Lindberg. This work considers non‐separable problems, but only on a rectangular domain with Dirichlet boundary conditions. We consider test cases with non‐smooth (i.e. discontinuous second derivatives) as well as smooth solutions. We also apply the multigrid method to an elliptic equation with non‐separable coefficients which occurs in a geothermal model. In this case an analysis of the error fails to show any advantage in a fourth‐order difference scheme over a second‐order scheme. However, we do demonstrate that the multigrid iteration performs well on this problem. Also, this example shows that the multigrid iteration can be combined with iterative improvement to create an efficient fourth‐order method for a non‐separable elliptic equation which is coupled with a marching equation. Other work has found an advantage in this fourth‐order scheme for a similar geothermal model.