About this Journal Submit a Manuscript Table of Contents
ISRN Mathematical Physics
Volumeย 2012ย (2012), Article IDย 234516, 11 pages
Research Article

Cubic Spline Iterative Method for Poissonโ€™s Equation in Cylindrical Polar Coordinates

1Department of Mathematics, Faculty of Mathematical Sciences, University of Delhi, Delhi 110 007, India
2Department of Mathematics, Deenbandhu Chhotu Ram University of Science & Technology, Murthal 131039, India

Received 4 October 2011; Accepted 16 November 2011

Academic Editors: J.-C.ย Wallet and H.ย Zhou

Copyright ยฉ 2012 R. K. Mohanty et al. 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.


Using nonpolynomial cubic spline approximation in ๐‘ฅ- and finite difference in ๐‘ฆ-direction, we discuss a numerical approximation of ๐‘‚(๐‘˜2+โ„Ž4) for the solutions of diffusion-convection equation, where ๐‘˜>0 and โ„Ž>0 are grid sizes in ๐‘ฆ- and ๐‘ฅ-coordinates, respectively. We also extend our technique to polar coordinate system and obtain high-order numerical scheme for Poissonโ€™s equation in cylindrical polar coordinates. Iterative method of the proposed method is discussed, and numerical examples are given in support of the theoretical results.

1. Introduction

We consider the two-dimensional elliptic equation of the form๐œ•2๐‘ข๐œ•๐‘ฅ2+๐œ•2๐‘ข๐œ•๐‘ฆ2=๐ท(๐‘ฅ)๐œ•๐‘ข๐œ•๐‘ฅ+๐‘”(๐‘ฅ,๐‘ฆ),(๐‘ฅ,๐‘ฆ)โˆˆฮฉ.(1.1) The Dirichlet boundary conditions are given by๐‘ข(๐‘ฅ,๐‘ฆ)=๐‘”(๐‘ฅ,๐‘ฆ),(๐‘ฅ,๐‘ฆ)โˆˆฮ“,(1.2) where ฮฉ={(๐‘ฅ,๐‘ฆ)โˆฃ0<๐‘ฅ,๐‘ฆ<1}is the solution domain and ฮ“ is its boundary. For ๐‘”(๐‘ฅ,๐‘ฆ)=0 and ๐ท(๐‘ฅ)=๐›ฝ, the above equation represents diffusion-convention equation. For ๐ท(๐‘ฅ)=โˆ’1/๐‘ฅ, the above equations represent the Poissonโ€™s equation with singular coefficients in rectangular coordinates. Similarly, for ๐ท(๐‘ฅ)=โˆ’1/๐‘ฅ and replacing the variables ๐‘ฅ,๐‘ฆ by ๐‘Ÿ,๐‘ง, we obtain a Poissonโ€™s equation in cylindrical polar coordinates. We will assume that the boundary conditions are given with sufficient smoothness to maintain the order of accuracy of the difference scheme and spline functions under consideration.

In this paper, we are interested in discussing a new approximation based on cubic spline polynomial for the solution of elliptic equation (1.1). In many practical problems, coefficients of the second derivatives term are small compared to the coefficients of the first derivatives term. These problems are called singular perturbation problems. Singularly perturbed elliptic boundary value problems are mathematically models of diffusion-convections process or related physical phenomenon. The diffusion term is the term involving the second-order derivative, and convective term is that involving the first-order derivative. During last three decades, several numerical schemes for the solution of elliptic partial differential equations have been developed by many researchers. First-Lynch and Rice [1] have discussed high-accuracy finite difference approximations to the solutions of elliptic partial differential equations. Boisvert [2] has discussed a class of high-order accurate discretization for the elliptic boundary value problems. Yavneh [3] has reported the analysis of fourth-order compact scheme for convection diffusion equation. A fourth-order difference method for elliptic equations with non-linear first derivative terms has been discussed by Jain et al. [4, 5]. In 1997, Mohanty [6] has derived order โ„Ž4 difference method for a class of 2D elliptic boundary value problems with singular coefficients. A new discretization method of order four for the numerical solution of 2D non-linear elliptic partial differential equations has been studied by Mohanty et al. [7โ€“9]. The use of cubic spline polynomial and its approximation plays an important role for the formation of stable numerical methods. In the past, many authors (see [10โ€“12]) have studied and analysed the use of cubic spline approximations in the solution of linear two-point boundary value problems. In 1983, Jain and Aziz [13] have developed a new method based on cubic spline approximations for the solution of two-point nonlinear boundary value problems. Later, Al-Said [14, 15] has discussed cubic spline methods for solving the system of second-order boundary value problems. Khan [16] has introduced parametric cubic spline approach for solving second order ordinary differential equations. Mohanty et al. [17, 18] have also reported a fourth-order accurate cubic spline alternating group explicit method for nonlinear singular two-point boundary value problems. Recently, Rashidinia et al. [19] have proposed a new cubic spline technique for two-point boundary value problems. Most recently, Mohanty and Dahiya [20] and Mohanty et al. [21] have used cubic spline polynomials and developed high order stable numerical methods for the solution of one space dimensional parabolic and hyperbolic partial differential equations. To the authors knowledge, no high-order method using cubic spline polynomial for the solution of 2D elliptic differential equations (1.1) has been discussed in the literature so far. In this paper, using nine-point compact cell (see Figure 1), we discuss a new compact cubic spline finite difference method of accuracy two in ๐‘ฆ- and four in ๐‘ฅ-coordinates for the solution of elliptic differential equation (1.1). In the next section, we discuss the derivation of the proposed cubic spline method. It has been experienced in the past that the cubic spline solutions for the Poissonโ€™s equation in polar coordinates usually deteriorate in the vicinity of the singularity. We overcome this difficulty by modifying the method in such a way that the solution retains its order and accuracy everywhere in the vicinity of the singularity. In Section 3, we discuss an iterative method. In Section 4, we compare the computed results with the results obtained by the method discussed in [8]. Concluding remarks are given in Section 5.

Figure 1: 9-point computational network.

2. The Approximation Based on Cubic Spline Polynomial

We consider our region of interest, a rectangular domain ฮฉ=[0,1]ร—[0,1]. We choose grid spacing โ„Ž>0 and ๐‘˜>0 in the directions ๐‘ฅ- and ๐‘ฆ- respectively, so that the mesh points (๐‘ฅ๐‘™,๐‘ฆ๐‘š) denoted by (๐‘™,๐‘š) are defined by ๐‘ฅ๐‘™=๐‘™โ„Ž and ๐‘ฆ๐‘š=๐‘š๐‘˜;๐‘™=0,1,โ€ฆ,๐‘+1;๐‘š=0,1,โ€ฆ,๐‘€+1, where ๐‘ and ๐‘€ are positive integers such that (๐‘+1)โ„Ž=1 and (๐‘€+1)๐‘˜=1. The mesh ratio parameter is denoted by ๐œ†=(๐‘˜/โ„Ž). The notations ๐‘ข๐‘™,๐‘šand๐‘ˆ๐‘™,๐‘š are used for the discrete approximation and the exact value of ๐‘ข(๐‘ฅ,๐‘ฆ) at the grid point (๐‘ฅ๐‘™,๐‘ฆ๐‘š), respectively. Similarly, we write ๐ท(๐‘ฅ๐‘™)=๐ท๐‘™ and ๐ท(๐‘ฅ๐‘™ยฑ1)=๐ท๐‘™ยฑ1.

At the grid point (๐‘ฅ๐‘™,๐‘ฆ๐‘š), we denote ๐‘Š๐‘Ž๐‘=๐œ•๐‘Ž+๐‘๐‘Š๐œ•๐‘ฅ๐‘™๐‘Ž๐œ•๐‘ฆ๐‘š๐‘,๐‘Š=๐‘ˆ,๐ทand๐‘”.(2.1)

Let ๐‘†๐‘š(๐‘ฅ) be the cubic spline interpolating polynomial of the function ๐‘ข(๐‘ฅ,๐‘ฆ๐‘š) between the grid point (๐‘ฅ๐‘™โˆ’1,๐‘ฆ๐‘š) and (๐‘ฅ๐‘™,๐‘ฆ๐‘š) and is given by๐‘†๐‘š๎€ท๐‘ฅ(๐‘ฅ)=๐‘™๎€ธโˆ’๐‘ฅ3๐‘€6โ„Ž๐‘™โˆ’1,๐‘š+๎€ท๐‘ฅโˆ’๐‘ฅ๐‘™โˆ’1๎€ธ3๐‘€6โ„Ž๐‘™,๐‘š+๎‚ต๐‘ˆ๐‘™โˆ’1,๐‘šโˆ’โ„Ž26๐‘€๐‘™โˆ’1,๐‘š๎‚ถ๎‚€๐‘ฅ๐‘™โˆ’๐‘ฅโ„Ž๎‚+๎‚ต๐‘ˆ๐‘™,๐‘šโˆ’โ„Ž26๐‘€๐‘™,๐‘š๎‚ถโ‹…๎‚€๐‘ฅโˆ’๐‘ฅ๐‘™โˆ’1โ„Ž๎‚,๐‘ฅ๐‘™โˆ’1โ‰ค๐‘ฅโ‰ค๐‘ฅ๐‘™;๐‘™=1,2,โ€ฆ,๐‘+1;๐‘š=0,1,2,โ€ฆ,๐‘€+1,(2.2) which satisfies at ๐‘šth-line parallel to ๐‘ฅ-axis the following properties:(i)๐‘†๐‘š(๐‘ฅ) coincides with a polynomial of degree three on each [๐‘ฅ๐‘™โˆ’1,๐‘ฅ๐‘™],๐‘™=1,2,โ€ฆ,๐‘+1and๐‘š=1,2,โ€ฆ๐‘€,(ii)๐‘†๐‘š(๐‘ฅ)โˆˆ๐ถ2[0,1], (iii)๐‘†๐‘š(๐‘ฅ๐‘™)=๐‘ˆ๐‘™,๐‘š,๐‘™=0,1,2,โ€ฆ,๐‘+1and๐‘š=1,2,โ€ฆ,๐‘€.

The derivatives of cubic spline function ๐‘†๐‘š(๐‘ฅ) are given by๐‘†โ€ฒ๐‘šโˆ’๎€ท๐‘ฅ(๐‘ฅ)=๐‘™๎€ธโˆ’๐‘ฅ2๐‘€2โ„Ž๐‘™โˆ’1,๐‘š+๎€ท๐‘ฅโˆ’๐‘ฅ๐‘™โˆ’1๎€ธ2๐‘€2โ„Ž๐‘™,๐‘š+๐‘ˆ๐‘™,๐‘šโˆ’๐‘ˆ๐‘™โˆ’1,๐‘šโ„Žโˆ’โ„Ž6๎€บ๐‘€๐‘™,๐‘šโˆ’๐‘€๐‘™โˆ’1,๐‘š๎€ป,๐‘†๐‘š๎…ž๎…ž๎€ท๐‘ฅ(๐‘ฅ)=๐‘™๎€ธโˆ’๐‘ฅโ„Ž๐‘€๐‘™โˆ’1,๐‘š+๎€ท๐‘ฅโˆ’๐‘ฅ๐‘™โˆ’1๎€ธโ„Ž๐‘€๐‘™,๐‘š,(2.3) where๐‘€๐‘™,๐‘š=๐‘†๐‘š๎…ž๎…ž๎€ท๐‘ฅ๐‘™๎€ธ=๐‘ˆ๐‘ฅ๐‘ฅ๐‘™,๐‘š=โˆ’๐‘ˆ๐‘ฆ๐‘ฆ๐‘™,๐‘š+๐ท๐‘™๐‘ˆ๐‘ฅ๐‘™,๐‘š+๐‘”๐‘™,๐‘š,๐‘™=0,1,2,โ€ฆ,๐‘+1,๐‘—=1,2,โ€ฆ,๐ฝ,(2.4)๐‘š๐‘™,๐‘š=๐‘†๎…ž๐‘š๎€ท๐‘ฅ๐‘™๎€ธ=๐‘ˆ๐‘ฅ๐‘™,๐‘š=๐‘ˆ๐‘™,๐‘šโˆ’๐‘ˆ๐‘™โˆ’1,๐‘šโ„Ž+โ„Ž6๎€บ๐‘€๐‘™โˆ’1,๐‘š+2๐‘€๐‘™,๐‘š๎€ป,๐‘ฅ๐‘™โˆ’1โ‰ค๐‘ฅโ‰ค๐‘ฅ๐‘™,(2.5) and replacing โ„Ž by โ€œโˆ’โ„Žโ€, we get๐‘š๐‘™,๐‘š=๐‘†๎…ž๐‘š๎€ท๐‘ฅ๐‘™๎€ธ=๐‘ˆ๐‘ฅ๐‘™,๐‘š=๐‘ˆ๐‘™+1,๐‘šโˆ’๐‘ˆ๐‘™,๐‘šโ„Žโˆ’โ„Ž6๎€บ๐‘€๐‘™+1,๐‘š+2๐‘€๐‘™,๐‘š๎€ป,๐‘ฅ๐‘™โ‰ค๐‘ฅโ‰ค๐‘ฅ๐‘™+1,(2.6) Combining (2.5) and (2.6), we obtain๐‘š๐‘™,๐‘š=๐‘†๎…ž๐‘š๎€ท๐‘ฅ๐‘™๎€ธ=๐‘ˆ๐‘ฅ๐‘™,๐‘š=๐‘ˆ๐‘™+1,๐‘šโˆ’๐‘ˆ๐‘™โˆ’1,๐‘šโˆ’โ„Ž2โ„Ž๎€บ๐‘€12๐‘™+1,๐‘šโˆ’๐‘€๐‘™โˆ’1,๐‘š๎€ป,(2.7) Further, from (2.5), we have๐‘š๐‘™+1,๐‘š=๐‘†๎…ž๐‘š๎€ท๐‘ฅ๐‘™+1๎€ธ=๐‘ˆ๐‘ฅ๐‘™+1,๐‘š=๐‘ˆ๐‘™+1,๐‘šโˆ’๐‘ˆ๐‘™,๐‘šโ„Ž+โ„Ž6๎€บ๐‘€๐‘™,๐‘š+2๐‘€๐‘™+1,๐‘š๎€ป,(2.8) and from (2.6), ๐‘š๐‘™โˆ’1,๐‘š=๐‘†๎…ž๐‘š๎€ท๐‘ฅ๐‘™โˆ’1๎€ธ=๐‘ˆ๐‘ฅ๐‘™โˆ’1,๐‘š=๐‘ˆ๐‘™,๐‘šโˆ’๐‘ˆ๐‘™โˆ’1,๐‘šโ„Žโˆ’โ„Ž6๎€บ๐‘€๐‘™,๐‘š+2๐‘€๐‘™โˆ’1,๐‘š๎€ป,(2.9) We consider the following approximations:๐‘ˆ๐‘ฆ๐‘ฆ๐‘™,๐‘š=๎€ท๐‘ˆ๐‘™,๐‘š+1โˆ’2๐‘ˆ๐‘™,๐‘š+๐‘ˆ๐‘™,๐‘šโˆ’1๎€ธ๎€ท๐‘˜2๎€ธ=๐‘ˆ๐‘ฆ๐‘ฆ๐‘™,๐‘š๎€ท๐‘˜+๐‘‚2๎€ธ,๐‘ˆ๐‘ฆ๐‘ฆ๐‘™+1,๐‘š=๎€ท๐‘ˆ๐‘™+1,๐‘š+1โˆ’2๐‘ˆ๐‘™+1,๐‘š+๐‘ˆ๐‘™+1,๐‘šโˆ’1๎€ธ๎€ท๐‘˜2๎€ธ=๐‘ˆ๐‘ฆ๐‘ฆ๐‘™+1,๐‘š๎€ท๐‘˜+๐‘‚2๎€ธ,๐‘ˆ๐‘ฆ๐‘ฆ๐‘™โˆ’1,๐‘š=๎€ท๐‘ˆ๐‘™โˆ’1,๐‘š+1โˆ’2๐‘ˆ๐‘™โˆ’1,๐‘š+๐‘ˆ๐‘™โˆ’1,๐‘šโˆ’1๎€ธ๎€ท๐‘˜2๎€ธ=๐‘ˆ๐‘ฆ๐‘ฆ๐‘™โˆ’1,๐‘š๎€ท๐‘˜+๐‘‚2๎€ธ,๐‘ˆ๐‘ฅ๐‘™,๐‘š=๎€ท๐‘ˆ๐‘™+1,๐‘šโˆ’๐‘ˆ๐‘™โˆ’1,๐‘š๎€ธ(2โ„Ž)=๐‘ˆ๐‘ฅ๐‘™,๐‘š+โ„Ž26๐‘ˆ30๎€ทโ„Ž+๐‘‚4๎€ธ,๐‘ˆ๐‘ฅ๐‘™+1,๐‘š=๎€ท3๐‘ˆ๐‘™+1,๐‘šโˆ’4๐‘ˆ๐‘™,๐‘š+๐‘ˆ๐‘™โˆ’1,๐‘š๎€ธ(2โ„Ž)=๐‘ˆ๐‘ฅ๐‘™+1,๐‘šโˆ’โ„Ž23๐‘ˆ30๎€ทโ„Ž+๐‘‚3๎€ธ,๐‘ˆ๐‘ฅ๐‘™โˆ’1,๐‘š=๎€ทโˆ’3๐‘ˆ๐‘™โˆ’1,๐‘š+4๐‘ˆ๐‘™,๐‘šโˆ’๐‘ˆ๐‘™+1,๐‘š๎€ธ(2โ„Ž)=๐‘ˆ๐‘ฅ๐‘™โˆ’1,๐‘šโˆ’โ„Ž23๐‘ˆ30๎€ทโ„Ž+๐‘‚3๎€ธ.(2.10)

Since the derivative values of ๐‘†๐‘š(๐‘ฅ) defined by (2.4), (2.7), (2.8), and (2.9) are not known at each grid point (๐‘ฅ๐‘™,๐‘ฆ๐‘š), we use the following approximations for the derivatives of ๐‘†๐‘š(๐‘ฅ).


Then at each grid point (๐‘ฅ๐‘™,๐‘ฆ๐‘š), a cubic spline finite difference method of Numerov type with accuracy of ๐‘‚(๐‘˜2+โ„Ž4) for the solution of differential equation (1.1) may be written as๐œ†2๐›ฟ2๐‘ฅ๐‘ˆ๐‘™,๐‘š+๐‘˜2๎‚ƒ12๐‘ˆ๐‘ฆ๐‘ฆ๐‘™+1,๐‘š+๐‘ˆ๐‘ฆ๐‘ฆ๐‘™โˆ’1,๐‘š+10๐‘ˆ๐‘ฆ๐‘ฆ๐‘™,๐‘š๎‚„=๐‘˜2๎‚ธ๐ท12๐‘™+1๐‘ˆ๐‘ฅ๐‘™+1,๐‘š+๐ท๐‘™โˆ’1๐‘ˆ๐‘ฅ๐‘™โˆ’1,๐‘š+10๐ท๐‘™๎๐‘ˆ๐‘ฅ๐‘™,๐‘š๎‚น+๐‘˜2๎€บ๐‘”12๐‘™+1,๐‘š+๐‘”๐‘™โˆ’1,๐‘š+10๐‘”๐‘™,๐‘š๎€ป+๐‘‡๐‘™,๐‘š,(2.12) where the local truncation error ๐‘‡๐‘™,๐‘š=๐‘‚(๐‘˜4+๐‘˜2โ„Ž4). Note that the method (2.12) is of ๐‘‚(๐‘˜2+โ„Ž4) for the numerical solution of (1.1). However, the method (2.12) fails to compute at ๐‘™=1, when ๐ท(๐‘ฅ) and/or ๐‘”(๐‘ฅ,๐‘ฆ) contains the singular terms like 1/๐‘ฅ,1/๐‘ฅ2, and so forth. For example, if ๐ท(๐‘ฅ)=1/๐‘ฅ, this implies ๐ท๐‘™โˆ’1=1/๐‘ฅ๐‘™โˆ’1, which cannot be evaluated at ๐‘™=1 (since ๐‘ฅ0=0). We modify the method (2.12) in such a manner, so that the method retains its order and accuracy everywhere in the vicinity of the singularity.

We use the following approximation:๐ท๐‘™+1=๐ท๐‘™+โ„Ž๐ท10+โ„Ž22๐ท20๎€ทโ„Ž+๐‘‚3๎€ธ,๐ท๐‘™โˆ’1=๐ท๐‘™โˆ’โ„Ž๐ท10+โ„Ž22๐ท20๎€ทโ„Žโˆ’๐‘‚3๎€ธ,๐‘”๐‘™+1,๐‘š=๐‘”๐‘™,๐‘š+โ„Ž๐‘”10+โ„Ž22๐‘”20๎€ทโ„Ž+๐‘‚3๎€ธ,๐‘”๐‘™โˆ’1,๐‘š=๐‘”๐‘™,๐‘šโˆ’โ„Ž๐‘”10+โ„Ž22๐‘”20๎€ทโ„Žโˆ’๐‘‚3๎€ธ.(2.13)

Substituting the approximations (2.13) into (2.12) and neglecting higher-order terms and local truncation error, we get๎€ท12+๐›ฟ2๐‘ฅ๎€ธ๐›ฟ2๐‘ฆ๐‘ข๐‘™,๐‘š+12๐œ†2๐›ฟ2๐‘ฅ๐‘ข๐‘™,๐‘š=๐œ†2โ„Ž2๎€ท12๐ท00+โ„Ž2๐ท20๎€ธ๎€ท2๐œ‡๐‘ฅ๐›ฟ๐‘ฅ๎€ธ๐‘ข๐‘™,๐‘š+๐œ†2โ„Ž2๎€ท๐ท10โˆ’๐ท200๎€ธ๐›ฟ2๐‘ฅ๐‘ข๐‘™,๐‘šโˆ’โ„Ž2๐ท10๐›ฟ2๐‘ฆ๐‘ข๐‘™,๐‘š+โ„Ž๐ท002๎€ท๐›ฟ2๐‘ฆ2๐œ‡๐‘ฅ๐›ฟ๐‘ฅ๎€ธ๐‘ข๐‘™,๐‘š+๐‘˜2๎€บ12๐‘”00+โ„Ž2๎€ท๐‘”20+๐ท10๐‘”00โˆ’๐ท00๐‘”10,๎€ธ๎€ป(2.14) where ๐œ‡๐‘ฅ๐‘ข๐‘™=(1/2)(๐‘ข๐‘™+1/2+๐‘ข๐‘™โˆ’1/2) and ๐›ฟ๐‘ฅ๐‘ข๐‘™=(๐‘ข๐‘™+1/2โˆ’๐‘ข๐‘™โˆ’1/2). Note that the cubic spline method (2.14) is of ๐‘‚(๐‘˜2+โ„Ž4) for the numerical solution of elliptic equation (1.1), which is also free from the terms 1/๐‘ฅ๐‘™ยฑ1 and hence can be computed for ๐‘™=1(1)๐‘,๐‘š=1(1)๐‘€.

3. Iterative Method

Now consider the convection-diffusion equation ๐œ•2๐‘ข๐œ•๐‘ฅ2+๐œ•2๐‘ข๐œ•๐‘ฆ2=๐›ฝ๐œ•๐‘ข๐œ•๐‘ฅ,0<๐‘ฅ,๐‘ฆ<1,(3.1) where ๐›ฝ>0 is a constant, and magnitude of ๐›ฝ determines the ratio of convection to diffusion term. Substituting ๐ท(๐‘ฅ)=๐›ฝ and ๐‘”(๐‘ฅ,๐‘ฆ)=0 into the difference scheme (2.14) and simplifying, we obtain a nine-point cubic spline difference scheme of ๐‘‚(๐‘˜2+โ„Ž4) accuracy for the solution of the convection-diffusion equation (3.1) given by๐›ผ0๐‘ข๐‘™,๐‘š+๐›ผ1๐‘ข๐‘™+1,๐‘š+๐›ผ2๐‘ข๐‘™โˆ’1,๐‘š+๐›ผ3๐‘ข๐‘™,๐‘š+1+๐›ผ4๐‘ข๐‘™,๐‘šโˆ’1+๐›ผ5๐‘ข๐‘™+1,๐‘š+1+๐›ผ6๐‘ข๐‘™+1,๐‘šโˆ’1+๐›ผ7๐‘ข๐‘™โˆ’1,๐‘š+1+๐›ผ8๐‘ข๐‘™โˆ’1,๐‘šโˆ’1=0,(3.2) where ๐‘…=๐›ฝโ„Ž/2 is called the cell Reynold number, and the coefficients ๐›ผ๐‘—,๐‘—=0,1,2,โ€ฆ,8 are defined by๐›ผ0=20+24๐œ†2๎‚ต๐‘…1+23๎‚ถ,๐›ผ1=โˆ’12๐œ†2๎‚ต๐‘…1โˆ’๐‘…+23๎‚ถ๐›ผ+2(1โˆ’๐‘…),2=โˆ’12๐œ†2๎‚ต๐‘…1+๐‘…+23๎‚ถ+2(1+๐‘…),๐›ผ3=๐›ผ4๐›ผ=โˆ’10,5=๐›ผ6=โˆ’(1โˆ’๐‘…),๐›ผ7=๐›ผ8=โˆ’(1+๐‘…).(3.3)

The scheme (3.2) may be written in matrix form๐€๐ฎ=๐›,(3.4) where ๐€ is a square matrix of order ๐‘๐‘€ร—๐‘๐‘€, ๐ฎ is the solution vector, and ๐› is the right hand side vector consisting of boundary values.

The coefficient matrix ๐€ has a block tridiagonal structure๐€=tri[โˆ’๐‹,๐ƒ,โˆ’๐”], with the submatrices โˆ’๐‹,๐ƒ, and โˆ’๐” given by๎€บ๐›ผโˆ’๐‹=tri8,๐›ผ4,๐›ผ6๎€ป๎€บ๐›ผ=โˆ’๐”,๐ƒ=tri2,๐›ผ0,๐›ผ1๎€ป.(3.5)

We focus on line stationary iterative methods for solving the linear system (3.5). The coefficient matrix ๐€ can be written as ๐€=๐ƒโˆ’๐‹โˆ’๐”, where ๐ƒ is block tri-diagonal matrix of ๐€, โˆ’๐‹ is strictly block lower triangular part, and โˆ’๐” is strictly block upper triangular part of matrix ๐€. The iteration matrices of the block Jacobi and block Gauss-Seidel methods are described by ๐†๐ฝ=๐ƒโˆ’1(๐‹+๐”),๐†GS=(๐ƒโˆ’๐‹)โˆ’1๐”.(3.6)

The matrix ๐€ has block tri-diagonal form and hence is block consistently ordered (see Varga [22]).

It can be verified that ๐›ผ0>0 and ๐›ผ๐‘—<0 for ๐‘—=1,2,โ€ฆ,8 provided |๐‘…|<1. One can also verify that ๐›ผ0=8๎“๐‘–=1||๐›ผ๐‘–||,(3.7) which implies that ๐€ is weakly diagonally dominant. Since ๐€ is reducible, we conclude that it is also an ๐Œ-matrix (see Varga [22]).

Applying the Jacobi iterative method to the scheme (3.2), we get the iterative scheme(1โˆ’๐‘…)๐‘ข(๐‘˜)๐‘™+1,๐‘š+1+๎‚ธ12๐œ†2๎‚ต๐‘…1โˆ’๐‘…+23๎‚ถ๎‚น๐‘ขโˆ’2(1โˆ’๐‘…)(๐‘˜)๐‘™+1,๐‘š+(1โˆ’๐‘…)๐‘ข(๐‘˜)๐‘™+1,๐‘šโˆ’1+10๐‘ข(๐‘˜)๐‘™,๐‘š+1โˆ’๎‚ธ20+24๐œ†2๎‚ต๐‘…1+23๐‘ข๎‚ถ๎‚น(๐‘˜+1)๐‘™,๐‘š+10๐‘ข(๐‘˜)๐‘™,๐‘šโˆ’1+(1+๐‘…)๐‘ข(๐‘˜)๐‘™โˆ’1,๐‘š+1+๎‚ธ12๐œ†2๎‚ต๐‘…1+๐‘…+23๎‚ถ๎‚น๐‘ขโˆ’2(1+๐‘…)(๐‘˜)๐‘™โˆ’1,๐‘š+(1+๐‘…)๐‘ข(๐‘˜)๐‘™โˆ’1,๐‘šโˆ’1=0.(3.8)

The propagating factors for the Jacobi iterative methods are given by๐œ‰๐‘—=15+6๐œ†2๎€ท๎€ท๐‘…1+2โŽกโŽขโŽขโŽฃโŽ›โŽœโŽœโŽ๎ƒŽ/3๎€ธ๎€ธ5cos(๐œ‹๐‘˜)+6๐œ†2๎‚ต๐‘…1+๐‘…+23๎‚ถโˆ’โŽžโŽŸโŽŸโŽ ร—โŽ›โŽœโŽœโŽ๎ƒŽ(1+๐‘…)(1โˆ’cos๐‘˜๐œ‹)6๐œ†2๎‚ต๐‘…1โˆ’๐‘…+23๎‚ถโŽžโŽŸโŽŸโŽ โŽคโŽฅโŽฅโŽฆ,โˆ’(1โˆ’๐‘…)(1โˆ’cos๐‘˜๐œ‹)cos๐œ‹โ„Ž(3.9) where (๐‘€+1)โ„Ž=1,(๐‘+1)๐‘˜=1. Consequently, the spectral radii ๐œŒ of the Jacobi and Gauss-Seidel matrices are related by ๐œŒ๎€ท๐†GS๎€ธ=๐œŒ(๐†๐ฝ)2.(3.10)

Hence, the associated iteration๐ฎ(๐‘˜+1)=๐†๐ฎ(๐‘˜)+๐œ(3.11) converges for any initial guess, where ๐† is either Jacobi or Gauss-Seidel iteration matrix. The Jacobi and Gauss-Seidel splitting are regular for both the line and point versions, and hence, they converge for any initial guess.

4. Numerical Results

If we replace the partial derivatives in (1.1) by the central difference approximations at the grid point (๐‘™,๐‘š), we obtain a central difference scheme (CDS) which is of ๐‘‚(๐‘˜2+โ„Ž2). We now solve the following two benchmark problems whose exact solutions are known. The right hand side homogeneous function and boundary conditions may be obtained by using the exact solution as a test procedure. We use block Gauss-Seidel iterative method (see [22โ€“24]) to solve the proposed scheme (2.14). In all cases, we have considered ๐ฎ(0)=0 as the initial guess, and the iterations were stopped when the absolute error tolerance |๐ฎ(๐‘˜+1)โˆ’๐ฎ(๐‘˜)|โ‰ค10โˆ’12 was achieved. In all cases, we have calculated maximum absolute errors (๐‘™โˆž-norm) for different grid sizes. All computations were performed using double-precision arithmetic.

Example 4.1 (Convection-Diffusion Equation). The problem is to solve (3.1) in the solution region 0<๐‘ฅ,๐‘ฆ<1 whose exact solution is given by ๐‘ข(๐‘ฅ,๐‘ฆ)=๐‘’๐›ฝ๐‘ฅ/2sin๐œ‹๐‘ฆ๎€บsinh๐œŽ2๐‘’โˆ’๐›ฝ/2๎€ปsinh๐œŽ๐‘ฅ+sinh๐œŽ(1โˆ’๐‘ฅ),(4.1) where ๐œŽ2=๐œ‹2+๐›ฝ2/4, ๐›ฝ>0.
The maximum absolute errors for ๐‘ข are tabulated in Tables 1 and 2.

Table 1: Example 4.1: The maximum absolute errors.
Table 2: Example 4.1: The maximum absolute errors (๐‘˜/โ„Ž2=64).

Example 4.2 (Poissonโ€™s Equation in Polar Cylindrical Coordinates). Consider the following: ๐œ•2๐‘ข๐œ•๐‘Ÿ2+๐œ•2๐‘ข๐œ•๐‘ง2+1๐‘Ÿ๐œ•๐‘ข๐œ•๐‘Ÿ=๐‘“(๐‘Ÿ,๐‘ง),0<๐‘Ÿ,๐‘ง<1.(4.2)
The exact solution is given by ๐‘ข(๐‘Ÿ,๐‘ง)=๐‘Ÿ2sinh๐‘Ÿcosh๐‘ง. The maximum absolute errors for ๐‘ข are tabulated in Tables 3 and 4.

Table 3: Example 4.2: The maximum absolute errors.
Table 4: Example 4.2: The maximum absolute errors (๐‘˜/โ„Ž2=64).

5. Concluding Remarks

The available numerical methods based on spline approximations for the numerical solution of 2D Poissonโ€™s equation are of ๐‘‚(๐‘˜2+โ„Ž2) accurate, which require nine grid points. In this paper, using the same number of grid points, we have discussed a new stable compact nine point cubic spline finite difference method of ๐‘‚(๐‘˜2+โ„Ž4) accuracy for the solution of Poissonโ€™s equation in polar cylindrical coordinates. For a fixed parameter ๐›พ=๐‘˜/โ„Ž2, the proposed method behaves like a fourth-order method, which is exhibited from the computed results.


  1. R. E. Lynch and J. R. Rice, โ€œHigh accuracy finite difference approximation to solutions of elliptic partial differential equations,โ€ Proceedings of the National Academy of Sciences of the United States of America, vol. 75, no. 6, pp. 2541โ€“2544, 1978. View at Scopus
  2. R. F. Boisvert, โ€œFamilies of High order Accurate discretization of some elliptic problems,โ€ SIAM Journal on Scientific Computing, vol. 2, pp. 268โ€“285, 1981.
  3. I. Yavneh, โ€œAnalysis of a fourth-order compact scheme for convection-diffusion,โ€ Journal of Computational Physics, vol. 133, no. 2, pp. 361โ€“364, 1997. View at Scopus
  4. M.K. Jain, R.K. Jain, and R.K. Mohanty, โ€œA fourth order difference method for elliptic equations with Non Linear first derivative terms,โ€ Numerical Methods for Partial Differential Equations, vol. 5, no. 2, pp. 87โ€“95, 1989.
  5. M. K. Jain, R. K. Jain, and R. K. Mohanty, โ€œFourth order difference methods for the system of 2-D nonlinear elliptic partial differential equations,โ€ Numerical Methods for Partial Differential Equations, vol. 7, pp. 227โ€“244, 1991.
  6. R. K. Mohanty, โ€œOrder h4 difference methods for a class of singular two space elliptic boundary value problems,โ€ Journal of Computational and Applied Mathematics, vol. 81, no. 2, pp. 229โ€“247, 1997. View at Scopus
  7. R. K. Mohanty and S. Dey, โ€œA new finite difference discretization of order four for (u/n) for two-dimensional quasi-linear elliptic boundary value problem,โ€ International Journal of Computer Mathematics, vol. 76, no. 4, pp. 505โ€“516, 2001. View at Scopus
  8. R. K. Mohanty, S. Karaa, and U. Arora, โ€œFourth order nine point unequal mesh discretization for the solution of 2D nonlinear elliptic partial differential equations,โ€ Neural, Parallel and Scientific Computations, vol. 14, no. 4, pp. 453โ€“470, 2006. View at Scopus
  9. R. K. Mohanty and S. Singh, โ€œA new fourth order discretization for singularly perturbed two dimensional non-linear elliptic boundary value problems,โ€ Applied Mathematics and Computation, vol. 175, no. 2, pp. 1400โ€“1414, 2006. View at Publisher ยท View at Google Scholar ยท View at Scopus
  10. W. G. Bickley, โ€œPiecewise cubic interpolation and two-point boundary problems,โ€ Computer Journal, vol. 11, no. 2, pp. 206โ€“208, 1968. View at Publisher ยท View at Google Scholar ยท View at Scopus
  11. D. J. Fyfe, โ€œThe use of cubic splines in the solution of two point boundary value problems,โ€ The Computer Journal, vol. 12, pp. 188โ€“192, 1969.
  12. R. P. Tewarson, โ€œOn the use of splines for the numerical solution of nonlinear two-point boundary value problems,โ€ BIT, vol. 20, no. 2, pp. 223โ€“232, 1980. View at Publisher ยท View at Google Scholar ยท View at Scopus
  13. M. K. Jain and T. Aziz, โ€œCubic spline solution of two-point boundary value problems with significant first derivatives,โ€ Computer Methods in Applied Mechanics and Engineering, vol. 39, no. 1, pp. 83โ€“91, 1983. View at Scopus
  14. E. A. Al-Said, โ€œSpline methods for solving system of second-order boundary-value problems,โ€ International Journal of Computer Mathematics, vol. 70, no. 4, pp. 717โ€“727, 1999. View at Scopus
  15. E. A. Al-Said, โ€œThe use of cubic splines in the numerical solution of a system of second-order boundary value problems,โ€ Computers and Mathematics with Applications, vol. 42, no. 6-7, pp. 861โ€“869, 2001. View at Publisher ยท View at Google Scholar ยท View at Scopus
  16. A. Khan, โ€œParametric cubic spline solution of two point boundary value problems,โ€ Applied Mathematics and Computation, vol. 154, no. 1, pp. 175โ€“182, 2004. View at Publisher ยท View at Google Scholar ยท View at Scopus
  17. R. K. Mohanty and D. J. Evans, โ€œA fourth order accurate cubic spline alternating group explicit method for non-linear singular two point boundary value problems,โ€ International Journal of Computer Mathematics, vol. 80, no. 4, pp. 479โ€“492, 2003. View at Publisher ยท View at Google Scholar ยท View at Scopus
  18. R. K. Mohanty, P. L. Sachdev, and N. Jha, โ€œAn O(h4) accurate cubic spline TAGE method for nonlinear singular two point boundary value problems,โ€ Applied Mathematics and Computation, vol. 158, no. 3, pp. 853โ€“868, 2004. View at Publisher ยท View at Google Scholar ยท View at Scopus
  19. J. Rashidinia, R. Mohammadi, and R. Jalilian, โ€œCubic spline method for two-point boundary value problems,โ€ International Journal of Engineering Science, vol. 19, pp. 39โ€“43, 2008.
  20. R. K. Mohanty and V. Dahiya, โ€œAn O(k2+kh2+h4) accurate two-level implicit cubic spline method for one space dimensional quasi-linear parabolic equations,โ€ American Journal of Computational Mathematics, vol. 1, no. 1, pp. 11โ€“17, 2011.
  21. R.K. Mohanty, R. Kumar, and V. Dahiya, โ€œCubic spline method for 1D wave equation in polar coordinates,โ€ ISRN Computational Mathematics, vol. 2012, Article ID 302923, 6 pages, 2012. View at Publisher ยท View at Google Scholar
  22. R.S. Varga, Matrix Iterative Analysis, Springer, New York, NY, USA, 2000.
  23. Y. Saad, Iterative methods for Sparse Linear Systems, SIAM, Philadelphia, Pa, USA, 2nd edition, 2003.
  24. L. A. Hageman and D. M. Young, Applied Iterative Methods, Dover, New York, NY, USA, 2004.