• Views 464
• Citations 0
• ePub 17
• PDF 216
`Journal of Computational Methods in PhysicsVolume 2014, Article ID 853198, 11 pageshttp://dx.doi.org/10.1155/2014/853198`
Research Article

## A Single Sweep AGE Algorithm on a Variable Mesh Based on Off-Step Discretization for the Solution of Nonlinear Burgers’ Equation

1Department of Applied Mathematics, Faculty of Mathematics & Computer Science, Akbar Bhawan, Chanakyapuri, New Delhi 110021, India
2Department of Mathematics, Faculty of Mathematical Sciences, University of Delhi, Delhi 110007, India

Received 4 October 2013; Accepted 1 December 2013; Published 16 January 2014

Copyright © 2014 R. K. Mohanty and Jyoti Talwar. 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

We discuss a new single sweep alternating group explicit iteration method, along with a third-order numerical method based on off-step discretization on a variable mesh to solve the nonlinear ordinary differential equation subject to given natural boundary conditions. Using the proposed method, we have solved Burgers’ equation both in singular and nonsingular cases, which is the main attraction of our work. The convergence of the proposed method is discussed in detail. We compared the results of the proposed iteration method with the results of the corresponding double sweep alternating group explicit iteration methods to demonstrate computationally the efficiency of the proposed method.

#### 1. Introduction

Consider the general nonlinear ordinary differential equation subject to essential boundary conditions where are finite constants.

We assume that for (i) is continuous,(ii) and exist and are continuous,(iii) and for some positive constant .

These conditions ensure that the boundary value problem (1) and (2) possesses a unique solution (see Keller [1]).

With the advent of parallel computers, scientists are focusing on developing finite difference methods with the property of parallelism. Working on this, in the early 1980s, Evans [2, 3] introduced the Group Explicit methods for large linear system of equations. Further he discussed the Alternating Group Explicit (AGE) method to solve periodic parabolic equations in a coupled manner. Mohanty and Evans applied AGE method along with various high order methods [4, 5] for the solution of two-point boundary value problems. Later, Sukon and Evans [6] introduced a Two-parameter Alternating Group Explicit (TAGE) method for the two-point boundary value problem with a lower order accuracy scheme. In 2003 Mohanty et al. [7] discussed the application of TAGE method for nonlinear singular two point boundary value problems using a fourth-order difference scheme. In 1990, Evans introduced the Coupled Alternating Group Explicit method [8] and applied it to periodic parabolic equations. Many scientists are applying these parallel algorithms to solve ordinary and partial differential equations [911].

Recently, Mohanty [12] has proposed a high order variable mesh method for nonlinear two-point boundary value problem. Mohanty and Khosla [13, 14] also devised a new third-order accurate arithmetic average variable mesh method for the solution of the boundary value problem (1) and (2), using three grid points, which is applicable to both singular and nonsingular problems. No special technique is required to handle singular problems. They also discussed the application of two-parameter double sweep alternating group explicit methods. Here, we discuss a new single sweep group explicit iteration method along with a third-order accurate variable mesh method based on two extra off-step grid points for the solution of the boundary value problem (1) and (2).

#### 2. Off-Step Discretization

Consider the solution interval with a nonuniform mesh such that . Let , , be the mesh size and let , , be the mesh ratio. Grid points are given by , . Let be the exact solution of at the grid point and approximated by . Let and .

The new third-order method is described as follows.

For , we first obtain the off-step discretization where , , and .

Now, let

Then, at each interior grid point , the proposed differential equation (1) is discretized by With the help of (3), from (5) it is easy to verify that , provided . However, for uniform mesh , the local truncation error associated with (5) becomes (see Chawla and Shivakumar [15]).

Note that the proposed method (5) is applicable to both singular and nonsingular problems. We do not require any special technique to handle singular problems (see Mohanty [13]). Further note that, the boundary values are given by and . If the differential equation (1) is linear we use the proposed single sweep AGE iterative method and in nonlinear case, we use the Newton-AGE iterative method to obtain the solution.

#### 3. Application to Singular Problems

##### 3.1. Linear Singular Problems

We now discuss the application of the proposed numerical method (5) to the linear differential equation with variable coefficients where represents a forcing function.

For and for and , the equation above represents singular equation in cylindrical and spherical symmetry, respectively.

Let us denote

Therefore applying the method (5) to the differential equation (6) and neglecting error term, we obtain a linear difference equation of the form where

The linear difference scheme (8) has a local truncation error of and is free from the terms and therefore can very easily be solved for in the region .

##### 3.2. Nonlinear Singular Problems

Now we consider the application to nonlinear differential equation (1). Neglecting the error term, we may rewrite the nonlinear difference equation (5) as Let We denote The Jacobian may be represented as

#### 4. Single Sweep AGE Method

##### 4.1. Description of the Method

In this section we discuss the single sweep AGE iteration method. Using the boundary conditions , , the linear difference equation (8) may be written in the matrix form as where

To implement the single sweep AGE iterative method, we split the coefficient matrix into two submatrices , where and satisfy the following conditions.(i) and are nonsingular for suitable choice of and .(ii)For any vectors and and , , it is “convenient” to solve the system explicitly; that is, and for vectors and , respectively.We will be concerned here with the situation where and are small block systems.

Now we discuss the case when is even (with , ).

Let So that the system (14) can be rewritten as Then a two-parameter AGE method for solving the afore mentioned system may be written as where is an intermediate vector.

Eliminating and combining (18) and (19), we obtain the iterative method or where

The new iterative method (20) or (21) is called the two-parameter single sweep AGE iterative method and the matrix is called the iteration matrix.

Now we discuss the single sweep AGE algorithm, when is even.

For simplicity, we denote and we define for .

By carrying out the necessary algebra in (20), we obtain the following algorithm.

For , For .

Let and then Finally, for , Similarly, we can write the single sweep AGE algorithm when is odd.

##### 4.2. Convergence Analysis

The single sweep AGE iteration method is given by or where The matrix is called the iteration matrix.

To prove the convergence of the method, we need to prove that , where denotes the spectral radius of .

Lemma 1. Let be sufficiently small. Then, the eigenvalues of and are all real.

Proof. Consider Therefore, we have , for .
Let , be the eigenvalues of . Then ’s are the roots of the quadratic equation Simplifying, we get The discriminants of the quadratic equations are Hence, the eigenvalues of are real.
In a similar manner, we can show that the eigenvalues of are real.

Now we give the sufficient condition for the convergence of the method.

Theorem 2. Let and , be the eigenvalues of and , respectively. If then the iterative method is convergent for the system (14).

Proof. Let
Since the off-diagonal entries of are negative, therefore , . Therefore, the diagonal entries of are positive.
The iteration matrix is given by Define where and , and It is easy to verify that and are symmetric. Therefore, the matrices and are also symmetric. Hence, Also, is symmetric, and therefore Therefore, we have
From (35) and (36), we have and for .
Hence, Also, from (37), we have Hence, we conclude that Thus, Similarly, we can prove that Using (44), (48), and (49), we get .

Hence, the convergence of the method (20) follows.

##### 4.3. Application of Single Sweep Newton-AGE Method

Now we discuss the single sweep Newton-AGE iterative method for the nonlinear difference equation (11). We follow the approaches given by Evans [8, 16].

Let us define Then the Jacobian of can be written as the th-order tridiagonal matrix Now with any initial vector , we define where is the solution of the nonlinear system

For the single sweep Newton-AGE method, we consider the case when is even. We split the matrix as , where and then we write single sweep Newton-AGE method as where are relaxation parameters and and are nonsingular.

Since and consist of submatrices, they can be easily inverted aswith , , , and .

The matrices and can be evaluated in a manner suitable for parallel computing. In order for this method to converge, it is sufficient that the initial vector be close to the solution.

Similarly, we can write the single sweep Newton-AGE algorithm when is odd.

#### 5. Results and Observations

We have applied the methods to the following three examples, whose exact solutions are known to us, and have compared the results with the corresponding double sweep AGE and Newton-AGE method [14]. For single sweep Newton-AGE method, we use the technique given by Evans [16]. The right hand side function and boundary conditions may be obtained using the exact solutions. Here, we have taken constant, . The value of the first mesh spacing on the left is given by

Therefore, given the value of and , we can calculate from the above relation and the remaining mesh points are determined by , . The initial vector is used in all iterative methods and the iterations were stopped when the absolute error tolerance was achieved.

Example 1 (linear singular problem). Consider The exact solution is . The root mean square errors (RMSE), the maximum absolute errors (MAE), and number of iterations both for single and double sweep AGE methods are tabulated in Table 1 for , , , , , and , . The graph of Example 1 of exact solution versus numerical solution is plotted in Figure 1.

Table 1: Example 1.
Figure 1: Graph of exact solution and the numerical solution for alpha = 2.0, sigma = 1.1, and for Example 1.

Example 2 (Burgers’ equation). Consider The exact solution is given by , where . The RMSE, MAE, and the number of iterations for both single and double sweep AGE methods are tabulated in Table 2 for , , , , and , . The graph of Example 2 of exact solution versus numerical solution for , , is plotted in Figure 2.

Table 2: Example 2.
Figure 2: Graph of the exact solution and the numerical solution for , sigma = 1.0, and for Example 2.

Example 3 (Burgers’ equation in polar coordinates). Consider The exact solution is given by . The RMSE, MAE, and the number of iterations for both single and double sweep AGE methods are tabulated in Table 3 for various values of , and . The graph of Example 3 of exact solution versus numerical solution for , , , is plotted in Figure 3.

Table 3: Example 3.
Figure 3: Graph of the exact solution and the numerical solution for alpha = 2.0, , sigma = 0.9, and for Example 3.

#### 6. Discussion and Conclusion

We have discussed a new single sweep AGE iterative method and three-point off-step method of accuracy on a variable mesh for the solution of nonlinear two point boundary value problems. To demonstrate the efficiency of the method, three examples including two nonlinear and singular cases are presented. The results obtained are compared with the corresponding double sweep AGE method and show superiority over the latter. The double sweep AGE method requires two sweeps to solve a problem, whereas the single sweep AGE method requires only one sweep to solve the problem. Experimentally, as compared to the double sweep method the corresponding single sweep method requires much less number of iterations as it uses less intermediate variables. The method can be extended to solve multidimensional problems and is suitable for use on parallel computers.

#### Conflict of Interests

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

#### Acknowledgment

This research was supported by the Council of Scientific and Industrial Research under research Grant no. 09/045(0836)2009-EMR-I.

#### References

1. H. B. Keller, Numerical Methods for Two Point Boundary Value Problems, Blaisdell, 1968.
2. D. J. Evans, “Parallel strategies for linear systems of equations,” International Journal of Computer Mathematics, vol. 81, no. 4, pp. 417–446, 2004.
3. D. J. Evans, “Group explicit methods for solving large linear systems,” International Journal of Computer Mathematics, vol. 17, pp. 81–108, 1985.
4. D. J. Evans and R. K. Mohanty, “Alternating group explicit method for the numerical solution of non-linear singular two-point boundary value problems using a fourth order finite difference method,” International Journal of Computer Mathematics, vol. 79, no. 10, pp. 1121–1133, 2002.
5. 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.
6. K. S. Sukon and D. J. Evans, “Two parameter AGE (TAGE) method for the solution of a tri-diagonal linear system of equations,” International Journal of Computer Mathematics, vol. 60, no. 3-4, pp. 265–278, 1996.
7. R. K. Mohanty, P. L. Sachdev, and N. Jha, “TAGE Method for nonlinear singular two point boundary value problems using a fourth order difference scheme,” Neural, Parallel and Scientific Computations, vol. 11, pp. 281–287, 2003.
8. D. J. Evans, “The solution of periodic parabolic equations by the coupled alternating group explicit (CAGE) iterative method,” International Journal of Computer Mathematics, vol. 34, pp. 227–235, 1990.
9. D. J. Evans and Pragya Jain, “The coupled reduced alternating group explicit (CRAGE) method,” Parallel Algorithms and Applications, vol. 2, pp. 193–208, 1993.
10. Q. Feng, “Explicit finite difference method for convection-diffusion equations,” in Proceedings of the World Congress on Engineering, vol. 2, pp. 1094–1097, London, UK, 2009.
11. B. Zheng and Q. Feng, “Parallel finite difference method for diffusion equations,” in Proceedings of the 15th American Conference on Applied Mathematics, pp. 60–62, 2009.
12. R. K. Mohanty, “A family of variable mesh methods for the estimates of (du/dr) and solution of non-linear two point boundary value problems with singularity,” Journal of Computational and Applied Mathematics, vol. 182, no. 1, pp. 173–187, 2005.
13. R. K. Mohanty, “A class of non-uniform mesh three point arithmetic average discretization for ${y}^{″}=f\left(x,y,{y}^{\prime }\right)$ and the estimates of ${y}^{\prime }$,” Applied Mathematics and Computation, vol. 183, no. 1, pp. 477–485, 2006.
14. R. K. Mohanty and N. Khosla, “Application of TAGE iterative algorithms to an efficient third order arithmetic average variable mesh discretization for two-point non-linear boundary value problems,” Applied Mathematics and Computation, vol. 172, no. 1, pp. 148–162, 2006.
15. 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.
16. D. J. Evans, “Iterative methods for solving non-linear two point boundary value problems,” International Journal of Computer Mathematics, vol. 72, no. 3, pp. 395–401, 1999.