Advances in Numerical Analysis

Volume 2013 (2013), Article ID 614508, 10 pages

http://dx.doi.org/10.1155/2013/614508

## Geometric Mesh Three-Point Discretization for Fourth-Order Nonlinear Singular Differential Equations in Polar System

^{1}Department of Mathematics, South Asian University, Akbar Bhawan, Chanakya Puri, New Delhi 110021, India^{2}Department of Mathematics, Faculty of Mathematical Sciences, University of Delhi, Delhi 110007, India

Received 9 June 2013; Revised 8 August 2013; Accepted 23 August 2013

Academic Editor: Rüdiger Weiner

Copyright © 2013 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

Numerical method based on three geometric stencils has been proposed for the numerical solution of nonlinear singular fourth-order ordinary differential equations. The method can be easily extended to the sixth-order differential equations. Convergence analysis proves the third-order convergence of the proposed scheme. The resulting difference equations lead to block tridiagonal matrices and can be easily solved using block Gauss-Seidel algorithm. The computational results are provided to justify the usefulness and reliability of the proposed method.

#### 1. Introduction

Consider the fourth-order boundary value problem: subject to the necessary boundary conditions:

where and , , , and are real constants and .

Or equivalently

subject to the natural boundary conditions:

Fourth-order differential equations occur in various areas of mathematics such as viscoelastic and inelastic flows, beam theory, Lifshitz point in phase transition physics (e.g., nematic liquid crystal, crystals, and ferroelectric crystals) [1], the rolls in a Rayleigh-Benard convection cell (two parallel plates of different temperature with a liquid in between) [2], spontaneous pattern formation in second-order materials (e.g., polymeric fibres) [3], the waves on a suspension bridge [4, 5], geological folding of rock layers [6], buckling of a strut on a nonlinear elastic foundation [7], traveling water waves in a shallow channel [8], pulse propagation in optical fibers [9], system of two reaction diffusion equation [10], and so forth.

The existence and uniqueness of the solution for the fourth and higher-order boundary value problems have been discussed in [11–14]. In the recent past, the numerical solution of fourth-order differential equations has been developed using multiderivative, finite element method, Ritz method, spline collocation, and finite difference method [15–18]. The determination of eigen values of self adjoint fourth-order differential equations was developed in [19] using finite difference scheme. The motivation of variable mesh technique for differential equations arises from the theory of electrochemical reaction-convection-diffusion problems in one-dimensional space geometry [20]. The geometric mesh method for self-adjoint singular perturbation problems using finite difference approximations was discussed in [21]. The use of geometric mesh in the context of boundary value problems was studied extensively in [22–24]. In this paper, we derive a geometric mesh finite difference method for the solution of fourth- and sixth-order differential boundary value problems with order of accuracy being three. The simplicity of the proposed method lies in its three-point discretization without any use of fictitious nodes. The scheme is compact and applicable to both singular and nonsingular problems. The resulting difference equations are solved using block Gauss-Seidel algorithm for linear case, and corresponding Newton’s method has been applied to nonlinear problems.

The paper is outlined in the following manner: in Section 2, the derivation of the method is discussed in detail. In Section 3, we define the procedure for numerical solution to singular problems in such a way that the method retains the order and accuracy even in the vicinity of singularity. In Section 4, algorithmic details are provided for the numerical solution of sixth-order differential equations. The convergence property has been discussed briefly in Section 4. The numerical illustrations based on geometric mesh as well as uniform mesh were provided in Section 5. The paper is concluded in the last section with future development and remarks.

#### 2. Derivation of the Numerical Scheme

We discretize the solution region such that . Let and be the nonuniform step size, and let be the geometric mesh ratio. Let and denote the exact solution values of and at the mesh and , be their approximate solutions, respectively.

Consider the following three-point geometric mesh discretizations for and : where

Define

It is easy to verify that where , , and so forth.

Now, let where and are free parameters to be determined.

With the help of (8) and (11), it follows that

Further, we define

With the help of (12)-(13), it follows that

With the help of approximations (10), (14), the difference equation (6) at each internal mesh , is approximated as follows: where

The difference scheme (15) to be of , the coefficients of in (16) must be zero, and hence we obtain

Thus the values associated with (17) are given by , and the local truncation errors given by (16) becomes , . However, for , the error reduces to .

#### 3. Application to Singular Problem

Consider the singular fourth-order linear differential equation in cylindrical polar coordinates: or equivalently

where , , and are the singular coefficients. The associated boundary conditions are given by (4). For , 1, or 2, the differential equation (18) shows planar, cylindrical, or spherical geometries (see, [25, 26]).

Applying the difference schemes (5) and (15) to (19) and (20), respectively, we obtain a system of coupled difference equations for : where

Note that the scheme (21) fails when the solution is to be determined at . We overcome this difficulty by modifying the scheme in such a way that the solutions retain order and accuracy even in the vicinity of singularity . We consider the following approximations:

Using the similar approximations of , and and neglecting terms, we can rewrite (21) with the help of compact operators as

where

The modified scheme (24) is free from the terms , hence easily solved for . The difference equation (24) along with the boundary conditions (4) gives a linear system of equations for the unknowns , , . The resulting block tridiagonal system can be easily solved using block Gauss-Seidel algorithm.

#### 4. Extension to Sixth-Order Differential Equations

The proposed method can be easily extended to the sixth-order differential equations: subject to the necessary boundary conditions:

or equivalently, subject to the natural boundary conditions:

We outline the similar algorithm for (28) as follows: where the values of are same as obtained in Section 2.

Then, the -approximations for (26) or (28) can be obtained by the following relations for :

The boundary conditions (27) are used to obtain values at for and , respectively. The numerical scheme may be implemented by neglecting terms from (31).

#### 5. Convergence Analysis

In this section, we derive the difference scheme of singular problem and investigate the convergence property of the proposed scheme. Consider the model problem or equivalently

along with the boundary conditions (2).

If and on , then the boundary value problem (33) has a unique solution. Under these assumptions (see, [27]), there are constants and such that , .

For the convergence, the coefficients , , and associated with (5) and (15) must be negative (see, [22]), from which we obtain the condition .

Now applying the methods (5) and (15) to (33) and using the similar technique discussed in Section 3 for singular coefficients and , we obtain the following system of difference equations: where Incorporating the boundary values , , , and , the system of difference equations (34) in the matrix-vector form can be written as where is the block tridiagonal matrix, , , and .

Let , , and , which satisfies

Let be the discretization error vector, and let , be the discretization errors at the node . Subtracting (37) from (36), we obtain the error equation: Also, we obtain

Thus for sufficiently small or equivalently as , we obtain the relations , and , . Hence, the graph of the matrix is strongly connected, and thus the matrix is irreducible (see, [28]).

Further, let be the sum of the th row sum of the matrix ; then we have the following.

For ,

For ,

For ,

This implies the following.

For ,

For ,

For ,

For sufficiently small value of , that is, in the limiting case as , we obtain

Hence we find that is monotone (see, [29, 30]). Consequently exists, and .

Let be the (, )th element of , and we define

From the theory of matrix, we know that

Thus the following bounds can be estimated with the help of series expansions.

For ,

For ,

For ,

With the help of (48), we obtain the following bounds:

From (38) and (52), we obtain the following error estimates:

This proves the third-order convergence of the proposed method. We generalize the above results in the following theorem.

Theorem 1. *The method given by (15) for the numerical solution of fourth-order singular differential equation (1) with sufficiently small and , gives a third-order convergent solution.*

#### 6. Computational Illustrations

To illustrate the geometric mesh finite difference method, we have solved both linear and nonlinear problems. The boundary conditions may be obtained from the analytical solution as a test procedure. The numerical accuracy of results are tested using maximum absolute errors and root mean square errors with the error tolerance being ≤10^{−15}. For the simplicity in computation, we choose , for and define the geometric mesh as follows ([24]):

The subsequent mesh spacing is determined by , . If the boundary value problems exhibit layer behaviour near the left boundary (see, [21]), 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 .

All the numerical computations are performed using long double length arithmetic in under Linux operating system with 2 GB operational memory.

*Example 1. *Consider the fourth-order linear problem (see, [31]) in the polar form:

The analytical solution is . The errors estimates for various values of are reported in Tables 1 and 2 for uniform mesh and geometric mesh , respectively.

*Example 2. *Consider the boundary value problems that arise from time-dependent Navier-Stokes equation (see, [32]) for axis symmetric flow of an incompressible fluid contained between infinite disks

The analytical solution is . The errors estimates are reported in Table 3 for various values of and .

*Example 3. *Consider the sixth-order linear singular problem:

The analytical solution is . The errors estimates for various values of are reported in Tables 4 and 5 for uniform mesh and geometric mesh , respectively.

*Example 4. *Consider the nonlinear problem

The analytical solution is . The errors estimates are reported in Table 6 for various values of and .

#### 7. Conclusion

The numerical results confirm that the proposed geometric mesh finite difference scheme converges and applicable to both singular and nonsingular differential equations. The numerical accuracy obtained using geometric mesh shows superiority over corresponding uniform mesh. The optimum mesh ratio parameter within the specified convergent region may be obtained by simulations. We have employed block Gauss-Seidel method to solve the block matrix systems. The method can be extended to general even-order nonlinear differential equations. Application to the proposed scheme to nonlinear singular elliptic problems is an open problem.

#### Acknowledgment

The authors are thankful to the Professor Rüdiger Weiner and unknown referee for their valuable suggestions which improve the quality of paper.

#### References

- R. M. Hornreich, M. Luban, and S. Shtrikman, “Critical behavior at the onset of k→-space instability on the
*λ*line,”*Physical Review Letters*, vol. 35, no. 25, pp. 1678–1681, 1975. View at Publisher · View at Google Scholar · View at Scopus - M. C. Cross and P. C. Hohenberg, “Pattern formation outside of equilibrium,”
*Reviews of Modern Physics*, vol. 65, no. 3, pp. 851–1112, 1993. View at Publisher · View at Google Scholar · View at Scopus - A. Leizarowitz and V. J. Mizel, “One dimensional infinite-horizon variational problems arising in continuum mechanics,”
*Archive for Rational Mechanics and Analysis*, vol. 106, no. 2, pp. 161–194, 1989. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. C. Lazer and P. J. Mckenna, “Large-amplitude periodic oscillations in suspension bridges. Some new connections with nonlinear analysis,”
*SIAM Review*, vol. 32, no. 4, pp. 537–578, 1990. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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–1004, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. J. Amick and J. F. Toland, “Homoclinic orbits in the dynamic phase-space analogy of an elastic strut,”
*European Journal of Applied Mathematics*, vol. 3, no. 2, pp. 97–114, 1992. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. Buffoni, A. R. Champneys, and J. F. Toland, “Bifurcation and coalescence of a plethora of homoclinic orbits for a Hamiltonian system,”
*Journal of Dynamics and Differential Equations*, vol. 8, no. 2, pp. 221–279, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - N. N. Akhmediev, A. V. Buryak, and M. Karlsson, “Radiationless optical solitons with oscillating tails,”
*Optics Communications*, vol. 110, no. 5-6, pp. 540–544, 1994. View at Google Scholar · View at Scopus - A. Doelman and V. Rottschä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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. P. Agarwal and P. R. Krishnamoorthy, “Boundary value problems for $n$th order ordinary differential equations,”
*Bulletin of the Institute of Mathematics*, vol. 7, no. 2, pp. 211–230, 1979. View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Schröder, “Numerical error bounds for fourth order boundary value problems, simultaneous estimation of $u(x)$ and ${u}^{\u2033}(x)$,”
*Numerische Mathematik*, vol. 44, no. 2, pp. 233–245, 1984. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - V. Shanthi and N. Ramanujam, “A numerical method for boundary value problems for singularly perturbed fourth-order ordinary differential equations,”
*Applied Mathematics and Computation*, vol. 129, no. 2-3, pp. 269–294, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. P. Agarwal and Y. M. Chow, “Iterative methods for a fourth order boundary value problem,”
*Journal of Computational and Applied Mathematics*, vol. 10, no. 2, pp. 203–217, 1984. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - W. K. Zahra, “A smooth approximation based on exponential spline solutions for nonlinear fourth order two point boundary value problems,”
*Applied Mathematics and Computation*, vol. 217, no. 21, pp. 8447–8457, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Rashidinia and M. Ghasemi, “B-spline collocation for solution of two-point boundary value problems,”
*Journal of Computational and Applied Mathematics*, vol. 235, no. 8, pp. 2325–2342, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. A. Usmani and P. J. Taylor, “Finite difference methods for solving ${\left(p\left(x\right){y}^{\prime \prime}\right)}^{\prime \prime}+q\left(x\right)y=r\left(x\right)$,”
*International Journal of Computer Mathematics*, vol. 14, no. 3-4, pp. 277–293, 1983. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Britz,
*Digital Simulation in Electrochemistry*, vol. 66 of*Lecture Notes in Physics*, Springer, Berlin, Germany, 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Google Scholar · View at Scopus - R. K. Mohanty, “A class of non-uniform mesh three point arithmetic average discretization for ${y}^{\u2033}=f(x,y,{y}^{\prime})$ and the estimates of y′,”
*Applied Mathematics and Computation*, vol. 183, no. 1, pp. 477–485, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at MathSciNet - S. R. K. Iyengar and P. Jain, “Spline finite difference methods for singular two point boundary value problems,”
*Numerische Mathematik*, vol. 50, no. 3, pp. 363–376, 1986. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. D. Russell and L. F. Shampine, “Numerical methods for singular boundary value problems,”
*SIAM Journal on Numerical Analysis*, vol. 12, pp. 13–36, 1975. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - W. Gautschi,
*Numerical Analysis*, Birkhause, 2011. - R. S. Varga,
*Matrix Iterative Analysis*, vol. 27 of*Springer Series in Computational Mathematics*, Springer, Berlin, 2000. View at Publisher · View at Google Scholar · View at MathSciNet - P. Henrici,
*Discrete Variable Methods in Ordinary Differential Equations*, John Wiley & Sons, New York, NY, USA, 1962. View at MathSciNet - D. M. Young,
*Iterative Solution of Large Linear Systems*, Academic Press, New York, NY, USA, 1971. View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. R. Elcrat, “On the radial flow of a viscous fluid between porous disks,”
*Archive for Rational Mechanics and Analysis*, vol. 61, no. 1, pp. 91–96, 1976. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet