In this article, we have considered for numerical solution of a Poisson and Laplace equation in a domain. we have presented a novel finite difference method for solving the system of the boundary value problems subject to Dirichlet boundary conditions. We have derived the solution of the Poisson and Laplace equations in a two-dimensional finite region. We present numerical experiments to demonstrate the efficiency of the method.
In this article, we propose a finite difference method for solving system of elliptic boundary value problems in a square region Ω = [a,b]×[a,b],
where ∂Ω is boundary of Ω.
Many physical problems in engineering and physics can be modeled mathematically. These modeled problems in general governed by partial differential equations. We study the Poisson/Laplace equation a specific partial differential equation and these equations describe the behavior of electric, gravitational, and fluid potentials in the study of electrodynamics and electrostatics, gravitation and fluid dynamics respectively. It is difficult to get analytical solution of the most of the partial differential equations that arise in mathematical models of physical phenomena. So we have to use numerical methods to approximate the solution of such problems. Finite-difference, finite element and finite volume method are three important methods to numerically solve partial differential equations. A powerful and oldest method for solving Poisson**** or Laplace*** equation subject to conditions on boundary is the finite difference method, which makes use of finite-difference approximations. The Finite difference method is very simple and effective. The emphasis in this article is on the development of numerical method not to prove theoretical concepts of convergence and existence. Thus existence and uniqueness of the solution to problems (1) is assumed. We further assume that problems (1) is well posed. We will not consider the specific assumption to ensure existence and uniqueness of the solution to problems (1). Though it is important to prove theorems on uniqueness, existence and convergence and that can be developed in the future.
There are different methods depends on the geometry of the domain to solve Poisson or Laplace equation. The finite element and finite difference methods are mostly preferred method respectively for irregular and regular domain. There also exist several studies in the literature in which the standard finite difference method applied to solve elliptic PDEs in irregular domain too , , , . We are considering Poisson and Laplace equation in a regular domain. For example, we know that the electric potential is related to the charge density by Poisson’s equation and in a charge-free region of space i.e as charge density becomes negligible, this becomes Laplace equation. Consider a problem in a domain where these two equations arise simultaneously and we wish to determine electric potential of this problem. This type of the problem may be considered as an obstacle boundary value problem and numerically solved by following the ideas in [5, 6].
In this article, we consider a finite difference method for numerical solution of system of elliptic equations (1) in a square domain. This is the usual second-order accurate scheme. The derivation of the method depends on difference approximation of the derivative on discrete mesh points in region of interest. So much research have been reported on numerical solution of elliptic partial differential equations, many of them are excellent work. But to best of our knowledge, neither similar problems (1) nor numerical method for the solution of problems (1) has been discussed in literature so for.
The present work is organized as follows. In section 2: we present novel finite difference approximation for system of elliptic equations. A novel finite difference method is presented so that the resulting difference equation need satisfies the boundary conditions exactly. A derivation of the present method will discuss in section 3. A local truncation error in the method and convergence will discuss in section 4 and 5. Finally, the application of the developed method will be present together with illustrative numerical results has been produced to show the efficiency of the method in section 6. A discussion and conclusion on the performance of the method will present in section 7.
2 The Finite Difference Method
Consider the square domain Ω = [a,b]×[a,b] for the solution of problems (1). Let h =(b-a)÷ N be the uniform mesh size in the x and y directions of the Cartesian coordinate system parallel to coordinate axes. Generate mesh poinIn this article, we consider a finite differencets (xi, yj), xi = a + ih,i = 0,1,2,⋅⋅⋅⋅ N and yj = a + jh, j = 0,1,2,⋅⋅⋅⋅ N. Let denote the mesh point (xi, yj) by (i,j). Let us denote the exact, approximate values of the solution of the problems (1) and source function f(x, y) at mesh point (i,j) by Ui,j, ui,j and fi,j respectively. Also we define mesh point
where Fi,j =
3 Derivation of the Method
If we write system of equations given by (3) at each mesh point, we will obtain N × N linear system of equations if source functions f(x, y) and g(x, y) are linear otherwise nonlinear system of equations in
4 Local Truncation Error
Thus we have obtained an expression of local truncation error term in case i = j = 1 i.e. first equation of (3). Similarly we can compute local error terms for different values and cases of i and j i.e. in other equations in (3) as follows,
5 Convergence of the method
We next discuss the convergence of the proposed method and under suitable conditions . We prove that the order of the convergence is O(h2). For each i,j = 1(1)N, let us define
Let define matrix S(N–1)2 × 1,
where M(N–1)2 ×(N–1)2 is block tridiagonal matrix and defined as:
in which blocks are
But U is exact solution of (3), which in matrix form can be written as
The matrix M is row diagonally dominant. Let M be the adjacency matrix of some graph G. We may easily prove that graph G is connected. From this fact it follows that adjacency matrix M is irreducible . By the row sum criterion it follows that M is monotone . Thus positive M–1 exist. Moreover M is regular, thus system of equations (3) can be solved by Gaussian elimination method by preserving its band structure. Thus from (12), we have
Thus the proposed difference method (3) converges and the order of the convergence is quadratic.
6 Numerical Experiments
The validity and effectiveness of the proposed finite difference method on uniform mesh have been tested on linear model problem. In our experiment we took a square as a region of integration and covered it with a uniform mesh of size h. We have used iterative method Gauss Seidel to solve the system of linear equations. In table, we have shown maximum absolute error MAU, computed for different values of N, using following formula,
All the computations in the experiment were performed on MS Window 2007 professional operating system in the GNU FORTRAN environment version -99 compiler(2.95 of gcc) running on Intel Duo core 2.20 Ghz PC . The stopping condition for iteration was either error of order 10–10 or number of iterations 103⋅
Consider a linear elliptic boundary value problems in a square region Ω = [0,1]×[0,1] which, when solving consists of
where a0,a1 and a2 are constants and determined  by the considering the continuity of u, ux and uy points on ∂Ω = ∂Ω1 ∩ ∂Ω2. In table 1, we have presented the time Etime in seconds and the computed MAU for a0 = 0.5, a1 = 1.0–
Maximum absolute errors in u(x, y) in Ω1 and Ω2 for problem 1.
In this article we have discussed finite difference method for numerical solution of system of elliptic equations. The method is theoretically second order accurate and Dirichlet boundary conditions incorporated into the method. In considered model problem we find that the computational performance of the method and accuracy depends on the constructed exact solution. The results obtained in model problem is not in good agreement to the estimated order of the method. In future we wish to improve proposed method. Though method developed on rectangular domain and used equally spaced grid mesh systems, it has good potential for efficient application to many problems on different geometries and dimensions, work in this specific direction is in progress.
We are interested in computing approximate numerical solution of the elliptic obstacle problem at the mesh points (xi, yj) in two dimension space. The development of the method in this article may be applied for problems in higher dimensions and different geometry. Fortunately we able to find some idea about the convergence of finite difference method under an appropriate condition but unable to measure maximum absolute error in terms of the mesh size h in numerical experiment. In numerical experiment we have considered exact solution and its first order partial derivatives are continuous at the boundary points
Communicated by Juan L.G. Guirao
Pandey, P. K., The Numerical Solution of Third Order Differential Equation Containing the First Derivative. Neural Parallel & Scientific Comp. 13, 297-304 (2005).
Pandey, P. K., Solving system of second order boundary value problems by Numerov type method. Journal of Science and Arts 15, No. 2(31), 143-150 (2015).
Jain, M.K., Iyenger, S. R. K. and Jain, R. K., Numerical Methods for Scientific and Engineering Computation (2/e). Willey Eastern Limited, New Delhi, 1987.
Henrici, P., Discrete Variable Methods in Ordinary Differential Equations, John Wiley and Sons, New York, (1962).
Avetik, A., Rafayel B. and Michael, P., Numerical Solution of the Two-Phase Obstacle Problem by Finite Difference Method. Armenian Journal of Mathematics, Vol. 7, Number 2, 164-182 (2015).