International Journal of Computational MathematicsVolume 2014 (2014), Article ID 527924, 12 pageshttp://dx.doi.org/10.1155/2014/527924

Research Article

The Convergence of Geometric Mesh Cubic Spline Finite Difference Scheme for Nonlinear Higher Order Two-Point Boundary Value Problems

Navnit Jha,¹ R. K. Mohanty,1,2 and Vinod Chauhan²

¹Department of Mathematics, South Asian University, Akbar Bhawan, Chanakyapuri, New Delhi 110021, India²Department of Mathematics, University of Delhi, Delhi 110007, India

Received 31 March 2014; Accepted 10 July 2014; Published 23 July 2014

Academic Editor: Riccardo Dondi

Copyright © 2014 Navnit Jha 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.

Abstract.

An efficient algorithm for the numerical solution of higher (even) orders two-point nonlinear boundary value problems has been developed. The method is third order accurate and applicable to both singular and nonsingular cases. We have used cubic spline polynomial basis and geometric mesh finite difference technique for the generation of this new scheme. The irreducibility and monotone property of the iteration matrix have been established and the convergence analysis of the proposed method has been discussed. Some numerical experiments have been carried out to demonstrate the computational efficiency in terms of convergence order, maximum absolute errors, and root mean square errors. The numerical results justify the reliability and efficiency of the method in terms of both order and accuracy.

1. Introduction

Consider the following nonlinear two-point boundary value problems of order : where , , and , are finite real constants and .

The higher order two-point boundary value problems play an important role in various areas of mathematical physics and engineering. The mathematical modeling of geological folding of rock layers [1], theory of plates and shell [2], waves on a suspension bridge [3], reaction diffusion equation [4], astrophysical narrow convection layers bounded by stable layers [5], viscoelastic and inelastic flows, deformation of beam and plate deflection theory [6–8], and so forth are some of the modeling problems in mathematical physics.

The analytical solution of (1) for the arbitrary choice of is difficult and thus we attempt to develop an economical computational method. The existence and uniqueness of the solutions of higher order boundary value problems have been discussed by Agarwal and Krishnamoorthy [9], O’Regan [10], Aftabizadeh [11], and Wei [12]. In the past, the approximate solution for the second, fourth, and/or sixth order two-point boundary value problems has been discussed using homotopy analysis by Liang and Jeffrey [13], reproducing kernel space by Yao and Lin [14], spline solution by Siddiqi and Twizell [15], and the Sinc-Galerkin and Sinc-Collocation methods by Rashidinia and Nabati [16]. The monotone iterative technique and quasilinearization method for the higher order ordinary differential equations have been analysed by Koleva and Vulkov [17]. The geometric mesh technique gains its importance from the theory of electrochemical reaction-convection-diffusion problems in one-dimensional space geometry [18]. The formulation of the geometric mesh finite difference approximations for the two-point singular perturbation problems was developed by Jain et al. [19] and Iyengar et al. [20]. Later, the applications of geometric mesh in the context of second order two-point boundary value problems were studied extensively by Mohanty et al. [21] and Kadalbajoo and Kumar [22].

There are some higher order approximation methods available for the numerical solution of differential equations. The most important issues associated with higher order schemes are more computing time and more memory spaces. Therefore, the lower order numerical method with very low computing time is worthwhile. This paper aims to develop a third order accurate numerical method based on cubic polynomial spline basis and geometric mesh finite difference approximations for the numerical solution of second order two-point boundary value problems. The method can be easily extended to fourth, sixth, and even higher order problems. The simplicity of the proposed method lies in its three-point discretizations with evaluations at two neighbouring nodes and one central node. The scheme is inherently compact and minor modifications are required for the singular problems. The resulting systems of algebraic equations are solved using block gauss elimination method obtained from the discretizations of linear differential equations. Classical Newton’s method has been applied to the nonlinear coupled difference equations.

Section 2 discussed how the third order accurate method based on geometric mesh polynomial spline basis and finite difference discretizations for second order differential equations was derived. The convergence analysis has been discussed briefly in Section 3. Further in Section 4, the proposed method has been extended to higher (even) order two-point boundary value problems; however, we restrict our computations up to the sixth order problems. In Section 5, some examples are presented which show singular and nonsingular nature both in linear and nonlinear cases of second, fourth, and sixth order problems.

2. Geometric Mesh Cubic Spline Methods

We introduce a finite set of geometric grids for the solution region . Let , be the nonuniform step size and let be the geometric mesh ratio. Let be the exact solution values and let be the approximate solution of at the mesh , respectively. Let be the cubic interpolating polynomial which interpolates the solution of the second order differential equation at defined as follows: where the coefficients , , , and are obtained using the following relations: We obtain via algebraic calculations the following expressions: Using the continuity of the first order derivatives , , we obtain Now, we define the following finite difference geometric mesh approximations so as to obtain the third order of accuracy for the second order differential equations with significant first order derivative: With the help of series expansion, we obtain Define With the help of (7)–(10), we obtain where .

Let , where is a free parameter to be determined. With the help of (7) and (12), we find Therefore, if we choose , we achieve .

Now consider the spline relation (2), for ; we can easily obtain the following approximations: It is easy to verify that We define Thus, we obtain Analogous to spline relation (5), with the help of (17), consider the following relationship: If the coefficients of , in the relation (18) vanish, then the right hand side of (18) has local truncation error of and hence we obtain Finally, we conclude the derivation for the approximate solution of second order differential equation with significant first order derivative and the corresponding geometric mesh spline approximation is given by The difference equation (21) has local truncation error of for . However, for the uniform mesh (), the local truncation errors become . If contains any singularity in the given domain , then we use finite Taylor’s expansion for the terms containing and then neglecting terms from (21) and incorporating the boundary values , , we obtain a system of recurrence relations with tridiagonal nature. If is linear, then the resulting matrix system can be solved by gauss elimination method; otherwise we use Newton’s solver for computing the numerical solutions. For the convergence of the difference equations (21), must be positive [23].

3. Convergence Analysis

In this section, we discuss the convergence theory for the geometric mesh cubic spline method (21) discussed for the numerical solution of second order two-point boundary value problems (20). At , (20) is written as Then, the geometric mesh cubic spline scheme (21) is given by where Scheme (23) in the matrix vector notation is written as where , .

Our aim is to find an approximation for by solving system From (26) and (27), we obtain Let and   be the discretization error vector.

Let us define where With the help of mean value theorem, (30) can be expressed as where , .

Now, from (24), we find that Equivalently, in the matrix notation, we have where is the tridiagonal matrix with the following entries: From (28) and (35), we obtain It is easy to see that for sufficiently small , that is, and .

Hence the graph of matrix is strongly connected and thus is irreducible [24].

Now let Let and .

Then, we obtain Thus, for sufficiently small values of , we find that Hence is monotone [25]; consequently we find that exists and [26]. If denotes the th element of and we define the matrix-vector norm as Also, from the theory of matrix, we know that With the help of Taylor’s expansion, we obtain Therefore, from (37) and (41), we obtain This proves the third order convergence of the geometric mesh spline method. The convergence theory may be now formally stated.

Theorem 1. The geometric mesh cubic spline finite difference method given by (21) for the numerical solution of second order two-point boundary value problems (20) with sufficiently small and gives a third order convergent solution.

4. Geometric Mesh Spline Algorithms for Higher Order Boundary Value Problem

The geometric mesh spline method, developed in Section 2, can be easily extended to higher order (even) boundary value problem (1) and it can be written as a system of second order boundary value problem: Then, the third order accurate algorithm based on the geometric mesh cubic spline method for the numerical solution of (45) may be written as follows: The boundary values are used to obtain values at for and , respectively. The numerical scheme may be implemented by neglecting terms from the system of (47). The resulting difference equations in case of linear boundary value problems give a block tridiagonal system of equations for the unknowns , , and can be easily solved using block gauss elimination method.

5. Numerical Results

To illustrate the efficiency of the cubic spline geometric mesh finite difference method, we have computed both nonsingular and singular equations for linear and nonlinear two-point boundary value problems of second, fourth, and sixth order. The boundary conditions are obtained using the analytical solution as a test procedure. The numerical accuracy of results is presented using maximum absolute errors (), root mean square errors (), and computational order of convergence () for interpolating th order derivative of with the error tolerance being : For the simplicity in computation, we have chosen , for , and computed the geometric mesh as follows: The subsequent mesh spacing is determined by , . If the boundary value problems exhibit layer behaviour near the left boundary, the solution value can be captured by choosing . If the layer occurs at the right boundary, we choose . If the layer occurs in the interior region, then mesh in the first half of the interval may be arranged by choosing and second half of the interval by choosing [27]. All the numerical computations are performed using long double length arithmetic in C under Linux operating system with 2 GB operational memory.

Example 1 (see [28]). Consider the linear nonsingular convection diffusion problem The analytical solution is . The numerical results are computed for = 10, 100 and 500 in Tables 1, 2, and 3 for various values of .

Table 1: Error estimates and order of convergence for Example 1 at .

20

—

—

40

4.0

3.9

80

3.0

3.9

Table 2: Error estimates and order of convergence for Example 1 at .

20

—

—

40

6.0

2.6

80

6.0

3.5

Table 3: Error estimates and order of convergence for Example 1 at .

20

—

—

40

5.5

0.7

80

0.1

1.3

Example 2 (see [29]). Consider the linear nonsingular fourth order stiff two-point boundary value problem: The theoretical solution is given by . We know that and are the eigenvalues of this equation and hence the problem is stiff for large values of . We have solved the problem for small as well as large values of and the behavior of solution is sufficiently smooth for in case of both uniform mesh and geometric mesh. We have computed the approximate solution with . Tables 4-5 present the accuracy of solution and computational order of convergence in case of uniform mesh and geometric mesh.

Table 4: Error estimates and order of convergence for Example 2 at .

20

—
—

40

3.9
3.9

80

3.9
3.9

Table 5: Error estimates for Example 2 at .

20

—
—

40

4.4
4.4

80

5.5
4.6

Example 3 (see [30]). Consider the linear nonsingular sixth order stiff two-point boundary value problem: The theoretical solution is given by . We know that and are the eigenvalues of this equation and hence the problem is stiff for large values of . We have solved the problem for small as well as large values of and the behavior of solution is sufficiently smooth in case of both uniform mesh and geometric mesh. We have computed the approximate solution with . Tables 6-7 present the accuracy of solution and computational order of convergence in case of uniform mesh and geometric mesh.

Table 6: Error estimates for Example 3 at .

20

—

40

3.9

80

3.9

Table 7: Error estimates for Example 3 at .

20

—

40

4.0

80

4.5

Example 4 (see [21]). Consider the linear singular second order problem: The above problem for represents Cartesian, cylindrical, and spherical symmetry, respectively. The unknown function is chosen in such a manner that is the theoretical solution. The accuracy of numerical and theoretical solutions is obtained for with different values of . The results are reported with geometric mesh as well as with uniform mesh. Tables 8-9 present the maximum absolute errors, root mean square errors, and computational order of convergence computed for various values of .

Table 8: Error estimates and order of convergence for Example 4 at .

20

—

—

40

4.8

3.8

80

6.7

3.9

Table 9: Error estimates and order of convergence for Example 4 at .

20

—

—

40

4.5

3.8

80

2.2

3.9

Example 5 (see [31]). Consider the linear singular fourth order problem: For , the above problem represents Cartesian, cylindrical, and spherical symmetry, respectively. The theoretical solution is given by . The accuracy of numerical and theoretical solutions is obtained for = 2 with different values of . The results are reported with geometric mesh as well as with uniform mesh. Tables 10-11 present the maximum absolute errors, root mean square errors, and computational order of convergence. The improvement in the accuracy of solutions is evident with geometric mesh as compared to spline method with uniform mesh.

Table 10: Error estimates for Example 5 at .

20

—
—

40

3.8
3.0

80

3.8
3.0

Table 11: Error estimates for Example 5 at .

20

—
—

40

5.6
4.6

80

1.3
0.5

Example 6 (see [32]). Consider the sixth order linear singular problem: The theoretical solution is given by . The problem has been solved numerically in the spherical polar coordinate (, ). The errors estimates for the various values of are reported in Tables 12-13 for uniform mesh as well as for geometric mesh .

Table 12: Error estimates for Example 6 at .

20

—

40

3.8

80

3.9

Table 13: Error estimates for Example 6 at .

20

—

40

3.5

80

3.4

Example 7 (see [33]). Consider the nonlinear nonsingular system that describes heat and mass transfer and an exothermal chemical reaction in a flat plate: The test solution is taken as and it has been observed that the behavior of the solution with uniform mesh deteriorates; however, the geometric mesh shows smooth behaviour of the solution. The accuracy of the solutions is computed for , , and in Table 14.

Table 14: Error estimates and order of convergence for Example 7.

20

—
Overflow
Overflow
—

40

3.0
Overflow
Overflow
—

80

2.0
Overflow
Overflow
—

Example 8 (see [17]). Consider the nonlinear nonsingular boundary value problems related to Fisher-Kolmogorov and Swift-Hohenberg equation: The analytical solution is . The errors estimates are reported in Tables 15-16 for , and for various values of using uniform mesh and geometric mesh, respectively.

Table 15: Error estimates for Example 8 at .

20

—
—

40

3.9
3.9

80

3.9
3.9

Table 16: Error estimates for Example 8 at .

20

—
—

40

3.3
2.8

80

2.9
2.0

Example 9 (see [34]). Consider the sixth order nonlinear nonsingular problem: The theoretical solution is given by . The accuracy of solutions has been reported in Tables 17-18 for various values of with uniform mesh and geometric mesh, respectively.

Table 17: Error estimates for Example 9 at .

20

—

40

3.9

80

3.9

Table 18: Error estimates for Example 9 at .

20

—

40

4.3

80

5.5

Example 10 (see [28]). Consider nonlinear singular Burger’s equations: The analytical solution is given by . The errors estimates for various values of are reported in Tables 19-20 for uniform mesh as well as for geometric mesh with high Reynold’s number .

Table 19: Error estimates and order of convergence for Example 10 at .

20

—

—

40

4.7

3.1

80

3.2

3.5

Table 20: Error estimates and order of convergence for Example 10 at .

20

—

—

40

4.0

3.1

80

5.4

3.5

Example 11 (see [35]). Consider the fourth order nonlinear singular problem which arises in the study of spreading of oil on water: The theoretical solution is given by . The errors estimates are reported in Tables 21-22 for various values of with constant mesh and geometric mesh, respectively. The computational order of convergence has very good agreement with the theoretical order of convergence.

Table 21: Error estimates for Example 11 at .

20

—
—

40

3.0
3.0

80

3.0
3.0

Table 22: Error estimates for Example 11 at .

20

—
—

40

4.6
4.5

80

7.4
4.8

6. Conclusion

We have developed an -accurate algorithm based on cubic spline finite difference approximations using geometric mesh for the numerical treatment of two-point boundary value problems of even higher orders. In this particular case, if the geometric mesh ratio parameter is unity, the method achieves fourth order of accuracy. The proposed method shows its unique importance in the equation governing heat and mass transfer and exothermal chemical reaction, where the classical uniform mesh fails to converge, whereas with small rearrangement of grids it shows superiority in terms of order and accuracy. The numerical results confirm that the proposed geometric mesh cubic spline finite difference scheme converges uniformly and is applicable to both singular and nonsingular differential equations. The optimum mesh ratio parameter near unity may be obtained by simulations. Applications of the proposed scheme to the parabolic problems are an open problem.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

C. J. Budd, G. W. Hunt, and M. A. Peletier, “Self-similar fold evolution under prescribed end-shortening,” Mathematical Geology, vol. 31, no. 8, pp. 989–1005, 1999.

S. Timoshenko and S. W. Krieger, Theory of Plates and Shell, McGraw-Hill, New York, NY, USA, 1959.

Y. Chen and P. J. McKenna, “Traveling waves in a nonlinearly suspended beam: theoretical results and numerical observations,” Journal of Differential Equations, vol. 136, no. 2, pp. 325–355, 1997.

A. Doelman and V. Rottsch{\"a}fer, “Singularly perturbed and nonlocal modulation equations for systems with interacting instability mechanisms,” Journal of Nonlinear Science, vol. 7, no. 4, pp. 371–409, 1997.

J. Toomre, J. R. Zahn, J. Latour, and E. A. Spiegel, “Stellar convection theory II: Single mode study of the second convection zone in A-type stars,” Astrophysics Journal, vol. 207, pp. 545–563, 1976.

D. B. Tien and R. A. Usmani, “Solving boundary value problems in plate deflection theory,” Simulation, vol. 37, no. 6, pp. 195–206, 1981.

A. R. Davies, A. Karageorghis, and T. N. Phillips, “Spectral Galerkin methods for the primary two point boundary value problems in modeling viscoelastic flows,” International Journal for Numerical Methods in Engineering, vol. 26, no. 3, pp. 647–662, 1988.

R. P. Agarwal and G. Akrivis, “Boundary value problems occurring in plate deflection theory,” Journal of Computational and Applied Mathematics, vol. 8, no. 3, pp. 145–154, 1982.

R. P. Agarwal and P. R. Krishnamoorthy, “Boundary value problems for th order ordinary differential equations,” Bulletin of the Institute of Mathematics. Academia Sinica, vol. 7, no. 2, pp. 211–230, 1979.

D. O'Regan, “Solvability of some fourth (and higher) order singular boundary value problems,” Journal of Mathematical Analysis and Applications, vol. 161, no. 1, pp. 78–116, 1991.

A. R. Aftabizadeh, “Existence and uniqueness theorems for fourth-order boundary value problems,” Journal of Mathematical Analysis and Applications, vol. 116, no. 2, pp. 415–426, 1986.

Z. Wei, “A necessary and sufficient condition for 2nth order singular super-linear m-point boundary value problems,” Journal of Mathematical Analysis and Applications, vol. 327, no. 2, pp. 930–947, 2007.

S. Liang and D. J. Jeffrey, “An efficient analytical approach for solving fourth order boundary value problems,” Computer Physics Communications, vol. 180, no. 11, pp. 2034–2040, 2009.

H. Yao and Y. Lin, “Solving singular boundary-value problems of higher even-order,” Journal of Computational and Applied Mathematics, vol. 223, no. 2, pp. 703–713, 2009.

S. S. Siddiqi and E. H. Twizell, “Spline solution of linear sixth order boundary value problems,” International Journal of Computer Mathematics, vol. 60, pp. 295–304, 1996.

J. Rashidinia and M. Nabati, “Sinc-Galerkin and sinc-collocation methods in the solution of nonlinear two-point boundary value problems,” Computational & Applied Mathematics, vol. 32, no. 2, pp. 315–330, 2013.

M. N. Koleva and L. G. Vulkov, “Two-grid quasilinearization approach to ODEs with applications to model problems in physics and mechanics,” Computer Physics Communications, vol. 181, no. 3, pp. 663–670, 2010.

D. Britz, Digital Simulation in Electrochemistry, vol. 666 of Lecture Notes Physics, Springer, Berlin, Germany, 2005.

M. K. Jain, S. R. K. Iyengar, and G. S. Subramanyam, “Variable mesh methods for the numerical solution of two-point singular perturbation problems,” Computer Methods in Applied Mechanics and Engineering, vol. 42, no. 3, pp. 273–286, 1984.

S. R. K. Iyengar, G. Jayaraman, and V. Balasubramanian, “Variable mesh difference schemes for solving a nonlinear Schrödinger equation with a linear damping term,” Computers & Mathematics with Applications, vol. 40, no. 12, pp. 1375–1385, 2000.

R. K. Mohanty, D. J. Evans, and N. Khosla, “An non-uniform mesh cubic spline {TAGE} method for non-linear singular two-point boundary value problems,” International Journal of Computer Mathematics, vol. 82, no. 9, pp. 1125–1139, 2005.

M. K. Kadalbajoo and D. Kumar, “Geometric mesh FDM for self-adjoint singular perturbation boundary value problems,” Applied Mathematics and Computation, vol. 190, no. 2, pp. 1646–1656, 2007.

M. M. Chawla and P. N. Shivakumar, “An efficient finite difference method for two-point boundary value problems,” Neural, Parallel & Scientific Computations, vol. 4, no. 3, pp. 387–395, 1996.

R. S. Varga, Matrix Iterative Analysis, vol. 27 of Springer Series in Computational Mathematics, Springer, Berlin , Germany, 2000.

P. Henrici, Discrete Variable Methods in Ordinary Differential Equations, John Wiley & Sons, New York, NY, USA, 1962.

D. M. Young, Iterative Solution of Large Linear Systems, Academic Press, New York, NY, USA, 1971.

N. Jha, “A fifth order accurate geometric mesh finite difference method for general nonlinear two point boundary value problems,” Applied Mathematics and Computation, vol. 219, no. 16, pp. 8425–8434, 2013.

R. K. Mohanty, “A class of non-uniform mesh three point arithmetic average discretization for and the estimates of ,” Applied Mathematics and Computation, vol. 183, no. 1, pp. 477–485, 2006.

S. D. Conte, “The numerical solution of linear boundary value problems,” SIAM Review, vol. 8, pp. 309–321, 1966.

M. A. Noor and S. T. Mohyud-Din, “Homotopy perturbation method for solving sixth-order boundary value problems,” Computers & Mathematics with Applications, vol. 55, no. 12, pp. 2953–2972, 2008.

J. Talwar and R. K. Mohanty, “A class of numerical methods for the solution of fourth-order ordinary differential equations in polar coordinates,” Advances in Numerical Analysis, vol. 2012, Article ID 626419, 20 pages, 2012.

N. Jha, R. K. Mohanty, and V. Chauhan, “Geometric mesh three-point discretization for fourth-order nonlinear singular differential equations in polar system,” Advances in Numerical Analysis, vol. 2013, Art. ID 614508, 10 pages, 2013.

M. Kubicek and V. Hlavacek, Numerical Solutions of Nonlinear Boundary Value Problems with Applications, Prentice Hall, Englewood Cliffs, NJ, USA, 1983.

M. A. Ramadan, I. F. Lashien, and W. K. Zahra, “A class of methods based on a septic non-polynomial spline function for the solution of sixth-order two-point boundary value problems,” International Journal of Computer Mathematics, vol. 85, no. 5, pp. 759–770, 2008.

Y. W. Wang, “Effect of spreading of material on the surface of a fluid—an exact solution,” International Journal of Non-Linear Mechanics, vol. 6, no. 2, pp. 255–262, 1971.