Laplace equation, numerical methods
Methods that replace the original boundary value problem by a discrete problem containing a finite number of unknows, such that if one finds a solution of the latter with suitable accuracy, this enables one to determine the solution of the original problem with given accuracy
;
depends on
and tends to infinity as
.
In the case of spatial variables
, the Laplace equation has the form
![]() |
and is the homogeneous Poisson equation. Boundary value problems for the Laplace equation are special cases of boundary value problems for the Poisson equation and more general equations of elliptic type (see [1]), and numerical methods for solving boundary value problems for equations of elliptic type (see [1], [2]) comprise many numerical methods for the Laplace equation. The specific character of the Laplace equation makes it possible to construct and use methods that have essentially better characteristics than methods for more general equations, although in practice one often prefers the simplicity of using general methods on a computer, rather than these possibilities.
The main numerical methods for equations of elliptic type are: projection-grid methods (finite-element methods) and difference methods. Both classes of methods are connected with the approximation of the original domain by a grid domain
containing
nodes of the grid and the construction of a system of algebraic equations
![]() | (*) |
for the values of the function defined at these nodes. In projection-grid methods, which are special cases of variational and projection methods, one uses the idea of approximation of the space of functions under consideration that contains the solution of the original problem by certain special finite-dimensional subspaces with given basis functions, and in the system (*) the vector consists of the coefficients of the expansion of the resulting approximation of the required solution with respect to the chosen basis. On the assumption that the solution of the original problem in a bounded domain
of the plane has the form
![]() |
where ,
, where
denotes a Sobolev space, and the functions
are specified and reflect the asymptotic behaviour of
close to singular points (corner points of the boundary and points where the type of boundary condition changes), for many types of domains
and mixed boundary value problems these methods make it possible, for example, to find the solution
to within
in
in
arithmetic operations (see [3]), and in a number of more special cases the estimates of the computational work are lowered to
and
.
In difference methods one commonly uses in one form or another an approximation of the derivatives by differences, and in the system (*) the vector consists of components that approximate the values of the solution at the nodes of the grid . The characteristics of the methods mentioned have been extensively studied for boundary value problems in bounded domains
of the plane. For example, for the Dirichlet condition
, where
is the boundary of
and the
are sufficiently smooth, on the basis of an improvement of the differentiability properties of the solution of the Laplace equation as one moves away from
one can construct the system (*) in such a way that the number
has equal order as the number
of points on
used to specify
to within
, and
can be found to within
in
arithmetic operations, and this makes it possible to find a solution of the original problem to within
at any fixed point of a strictly interior subdomain in a finite number of operations (see [4]). Methods of this type are asymptotically optimal; in the case when one uses, for example, simpler methods with a rectangular grid that have accuracy
, the number of operations required to find
to within
is
;
(see [4]). Estimates of the error in the grid method for the Laplace equation have been studied in more detail (see [4], [5]); if there are singular points on
it is advisable to use the special structure of the grid close to these points (see [5]). One often uses difference methods based on the approximation of certain integral characteristics for the Laplace equation (see [6]).
Fairly infrequently one uses collocation methods, in which the system (*) is obtained as a consequence of the original equation being satisfied at nodes of the grid and the assumption that an approximation to a solution of the original problem is sought in some finite-dimensional subspace. A special class of numerical methods for solving boundary value problems for the Laplace equation is based on the reduction of these problems to singular integral equations and the subsequent solution of the resulting integral equations by numerical methods (see [7], [8]).
References
[1] | J.P. Aubin, "Approximation of elliptic boundary-value problems" , Wiley (1972) |
[2] | P.G. Ciarlet, "The finite element method for elliptic problems" , North-Holland (1978) |
[3] | E.G. Dyakonov, "Minimization of computational work. Asymptotically optimal algorithms for elliptic problems" , Moscow (1989) (In Russian) |
[4] | N.S. Bakhvalov, "About optimization of numerical methods" , Proc. Internat. Congress Mathematicians Nice, 1970 , 3 , Gauthier-Villars (1971) pp. 289–295 |
[5] | E.A. Volkov, "On the smoothness of solutions of the Dirichlet problem and the composite mesh method on polyhedra" Proc. Steklov Inst. Math. , 150 (1979) pp. 71–104 Trudy Mat. Inst. Steklov. , 150 (1979) pp. 67–98 |
[6] | A.A. Samarskii, I.V. Fryazinov, "Difference approximation methods for problems of mathematical physics" Russian Math. Surveys , 31 : 6 (1976) pp. 179–213 Uspekhi Mat. Nauk , 31 : 6 (1976) pp. 167–197 |
[7] | P.K. Banerjee, R. Butterfield, "Boundary element methods in engineering science" , McGraw-Hill (1981) |
[8] | S.L. Grouch, A.M. Starfield, "Boundary element methods in solid mechanics" , Allen & Unwin (1983) |
Comments
References
[a1] | J.R. Rices, R.F. Boisvert, "Solving elliptic problems using ELLPACK" , Springer (1984) |
[a2] | G. Birkhoff, R.E. Lynch, "Numerical solution of elliptic problems" , SIAM (1984) |
[a3] | R.S. Varga, "Matrix iterative analysis" , Prentice-Hall (1962) |
[a4] | E.L. Wachspress, "Iterative solution of elliptic systems, and applications to the neutron diffusion equations of nuclear physics" , Prentice-Hall (1966) |
[a5] | D.M. Young, "Iterative solutions of large linear system" , Acad. Press (1971) |
[a6] | G.J. Fix, "An analyse of the finite element method" , Prentice-Hall (1973) |
[a7] | A.R. Mitchell, D.F. Griffiths, "The finite difference method in partial differential equations" , Wiley (1980) |
[a8] | C. Johnson, "Numerical solutions of partial differential equations by the finite element method" , Cambridge Univ. Press (1987) |
Laplace equation, numerical methods. Encyclopedia of Mathematics. URL: http://encyclopediaofmath.org/index.php?title=Laplace_equation,_numerical_methods&oldid=17252