- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Abstract and Applied Analysis
Volume 2013 (2013), Article ID 741278, 12 pages
Fast Spectral Collocation Method for Solving Nonlinear Time-Delayed Burgers-Type Equations with Positive Power Terms
1Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
2Department of Mathematics, Faculty of Science, Beni-Suef University, Beni Suef 62511, Egypt
3Department of Mathematics, Faculty of Science, Um-Al-Qurah University, Makkah 21955, Saudi Arabia
Received 20 January 2013; Revised 6 June 2013; Accepted 6 June 2013
Academic Editor: Mustafa Bayram
Copyright © 2013 A. H. Bhrawy 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.
Since the collocation method approximates ordinary differential equations, partial differential equations, and integral equations in physical space, it is very easy to implement and adapt to various problems, including variable coefficient and nonlinear differential equations. In this paper, we derive a Jacobi-Gauss-Lobatto collocation method (J-GL-C) to solve numerically nonlinear time-delayed Burgers-type equations. The proposed technique is implemented in two successive steps. In the first one, we apply nodes of the Jacobi-Gauss-Lobatto quadrature which depend upon the two general parameters , and the resulting equations together with the two-point boundary conditions constitute a system of ordinary differential equations (ODEs) in time. In the second step, the implicit Runge-Kutta method of fourth order is applied to solve a system of ODEs of second order in time. We present numerical results which illustrate the accuracy and flexibility of these algorithms.
Spectral methods have emerged as powerful techniques used in applied mathematics and scientific computing to numerically solve linear and nonlinear differential equations [1–4] and integral equations [5–7]. Also, they have become increasingly popular for solving fractional differential equations [8–10]. The main idea of spectral methods is to put the solution of the problem as a sum of certain basic functions and then to choose the coefficients in the sum in order to minimize the difference between the exact solution and the approximate one as well as possible. The choice of test functions leads to the three well-known types of spectral methods, namely, the Galerkin, tau, and collocation methods [11–14]. Spectral collocation method has an exponential convergence rate, which is very useful in providing highly accurate solutions to nonlinear differential equations even using a small number of grids.
The Jacobi polynomials satisfy the orthogonality condition on the interval with respect to the weight function . There are many special cases of the Jacobi polynomials such as Gegenbauer, Legendre, Zernike, ultraspherical, and Chebyshev polynomials . In recent decades, the use of Jacobi polynomials for solving differential equations has gained increasing popularity due to obtaining the solution in terms of the Jacobi parameters and (see, e.g., [16, 17]).
Time-delay partial differential equations are a type of differential equations in which the derivative of the unknown function at a certain time is given in terms of the values of the function at previous times. They have a wide range of applications in science and engineering such as physical, chemical, and biological sciences [18, 19]. Also, time-delayed nonlinear partial differential equations described the propagation and transport phenomena or population dynamics.
The solution of delay differential equations has been presented by many authors, but we briefly review some of them. In , Ghasemi and Kajani developed a numerical approach of time-varying delay systems using Chebyshev wavelets. Hybrid functions together with Legendre polynomials are investigated in [21–23] to obtain efficient numerical solution of delay systems. Sedaghat et al.  presented a numerical scheme based on Chebyshev polynomials to treat the delay differential equations of pantograph type. The authors in  proposed a Bernoulli operational matrix method for solving generalized pantograph equation. Recently, Ali et al.  implemented spectral Legendre approach for solving pantograph-type differential and integral equations with studying the error analysis of the method. More recently, the work of Trif  discussed the application of the Tau method based on operational matrix of Chebyshev polynomials for solving DDEs of pantograph-type.
Recently, the authors of  and  presented some new travelling wave solutions of Burgers equation with finite transport memory and the Korteweg-de Vries-Burgers equation, respectively. Pandey et al.  investigated Du Fort-Frankel finite difference approach for solving Burgers equation in finite domain. Meanwhile, in , Sun and Wu presented and applied an efficient numerical solution for the Burgers equation based on a difference scheme in an unbounded domain. More recently, a differential quadrature scheme for Burgers equations was proposed in . The idea of commutative hypercomplex mathematics and the homotopy perturbation method were combined to investigate solutions of time-delayed Burgers equation by Rostamy and Karimi . The Darboux transformation was described in  to determine the exact solutions to the Burgers equation. Caglar and Ucar  proposed the nonpolynomial cubic spline scheme to develop a numerical solution of Burgers equation. Regarding the two-dimensional Burgers equation, Wang and Zhao  presented a novel combination of two-dimensional Haar wavelet functions based on tensorial products for solving two-dimensional Burgers equation.
Time-delayed Burgers equation has a wide range of application in many areas of applied sciences, for example, forest fire, population growth models, and Neolithic transitions [37, 38]. The time-delayed Burgers-Fisher equation is a very important model to forest fire, population growth, Neolithic transitions, the interaction between the reaction mechanism, convection effect and diffusion transport, and so forth . Fahmy et al.  used improved tanh function, variational iteration, and the Adomian decomposition methods to present some exact solutions and numerical solutions of the time-delayed Burgers equation. With the aid of a subsidiary high-order ODE, Zhang et al.  obtained the exact solutions of the generalized time-delayed Burgers-Fisher equation with positive fractional power terms. Jawad et al.  introduced the exact solution of time-delayed Burgers equation using complex tanh method. Moreover, ()-expansion method is applied by Kim and Sakthivel  to find the exact solutions of time-delayed Burgers equation.
The goal of this paper is to propose an orthogonal collocation scheme for solving three nonlinear time-delay partial differential equations based on Jacobi family in which the roots of the Jacobi orthogonal polynomials whose distributions can be tuned by two parameters and . Firstly, we apply the Jacobi-Gauss-Lobatto collocation (J-GL-C) method to the model equation for discretizing spatial derivatives, using nodes of the Jacobi-Gauss-Lobatto quadrature which depends upon the two general parameters , ; these equations together with the two-point boundary conditions constitute system of ordinary differential equations (ODEs) in time. Secondly, the Runge-Kutta method of fourth order is investigated for the time integration of the resulting system of nonlinear second-order ODEs.
Indeed, the main advantage of the proposed technique is that the Legendre and Chebyshev collocation methods and other methods can be obtained as special cases from our proposed technique. Comparison of the results obtained by various choices of Jacobi parameters and reveals that the present method is very effective and convenient for all choices of and . Finally, the accuracy of the proposed method is showed by test problems. From the results, these algorithms are extremely efficient and accurate for solving nonlinear time-delayed Burgers'-type equations.
The rest of this paper is structured as follows. In the next section, some basic properties of Jacobi polynomials, which are required in our paper, are given. Section 3 is devoted to the development of Gauss-Lobatto collocation technique for a general form of time-delay partial differential equation based on the Jacobi polynomials, and in Section 4 the proposed method is applied to obtain some numerical results for three problems of time-delay partial differential equations with known exact solutions. Finally, a brief conclusion and some remarks are provided in Section 5.
In this section, we briefly recall some properties of the Jacobi polynomials , , ), which are satisfying the following relations:
The th derivative of Jacobi polynomials of degree can be given by
These polynomials are the only polynomials arising as eigenfunctions of the following singular Sturm-Liouville equation:
Let ; then we define the weighted space as usual. The inner product and the norm of with respect to the weight function are defined as follows:
The set of Jacobi polynomials forms a complete -orthogonal system, and
3. Jacobi Spectral Collocation Method
The main objective of this section is to develop the J-GL-C method to numerically solve the nonlinear time-delayed Burgers-type equations: where
with the boundary-initial conditions where , , , and are real numbers and is the time delay. Now, suppose that the change of variables , , which will be used to transform problem (6)–(9) into another one in the classical interval for the space variable, to directly implement collocation method based on Jacobi family defined in , where
with the boundary-initial conditions
The aim of this work is to consider the advantage of the collocation point distribution in a specified domain using the roots of the Jacobi orthogonal polynomials whose distributions can be tuned by two parameters, and . Now, we outline the main step of the J-GL-C method for solving nonlinear time-delayed Burgers-type equation. Let us expand the dependent variable in a Jacobi series, and in virtue of (4)-(5), we deduce that
To evaluate the previous integral accurately, we present the Jacobi-Gauss-Lobatto quadrature. For any , where is the set of polynomials of degree less than or equal to , () and () are the nodes and the corresponding Christoffel numbers of the Jacobi-Gauss-Lobatto quadrature formula on the interval , respectively.
In accordance with (15), the coefficients in terms of the solution at the collocation points can be approximated by
Therefore, (13) can be rewritten as
Furthermore, if we differentiate (17) once (using (2)) and evaluate it at all Jacobi-Gauss-Lobatto collocation points, it is easy to compute the first spatial partial derivative in terms of the values at these collocation points as
or it can be shortened to where
Similar steps can be applied to the second spatial partial derivative to get where
In the proposed Jacobi-Gauss-Lobatto collocation method, the residual of (6) is set to zero at of Jacobi-Gauss-Lobatto points; moreover, the boundary conditions (8) will be enforced at the two collocation points and . Therefore, adopting (19)–(22) enables one to write (6)–(8) in the form where
This provides an system of second-order ordinary differential equations in the expansion coefficients , namely, where
subject to the initial values
Remark 1. We can replace the Jacobi polynomials by the Legendre, Chebyshev of the first, second, third, and fourth kinds, or Gegenbauer polynomials and replace the nodes of Jacobi-Gauss-Lobatto quadrature by the Legendre-, Chebyshev-, or Gegenbauer-Gauss-Lobatto quadrature (cf. [42–45]), just by taking the special cases , , , , , or respectively, in the resulting system of ordinary differential equations.
The SODEs (29) can be solved by using implicit Runge-Kutta method of fourth order, which is an important family of implicit and explicit iterative methods for the approximation of solution of system of ordinary differential equations. In the next section, we recall that the difference between the measured or inferred value of approximate solution and its actual value (absolute error) is given by where and are the exact solution and the approximate solution at the point , respectively. Moreover, the maximum absolute error is given by
4. Numerical Results and Discussions
This section considers three numerical examples to demonstrate the accuracy and applicability of the proposed method in the present paper. Comparison of the results obtained by various choices of Jacobi parameters and reveals that the present method is very effective and convenient for all choices of and .
We consider the following three examples.
Maximum absolute errors of (34) subject to (35) and (36) are introduced in Table 1 using J-GL-C method with various choices of , and , and ; this table indicates that the obtained results are very accurate for a small choice of , while the absolute errors of problem (34) are presented in Table 2 for , , , and with different values of in the interval .
In Figure 1, we see that the approximate solution and the exact solution for of problem (34) coincide with values of parameters listed in their caption. Moreover, the absolute error of problem (34) where and is displayed in Figure 2. In Figure 3, we plotted the approximate solution . This assertion shows that the obtained numerical results are accurate and compare favorably with the exact solution.
Example 3. Consider the nonlinear time-delayed one-dimensional generalized Burgers equation in the form
subject to the boundary conditions
and the initial conditions
The exact solution  of (39) is
Table 3 lists the maximum absolute errors of (39) subject to (40) and (41) in , using J-GL-C method for different values of ,, and . Moreover, in Tables 4 and 5, we evaluate the absolute errors of (39) for and , respectively, with different values of in the interval .
Figure 4 plots the approximate and the exact solutions at of (39) with the special values , , , and . Moreover, the absolute errors between the approximate and the exact solutions of (39) with , , , and are shown in Figure 5.
Maximum absolute errors of (44) subject to (45) and (46) are tabulated in Table 6 using J-GL-C method with various choices of , , and, while the absolute errors of (44) are presented in Table 7 for and with different values of in the interval .
In Figure 6, we see that the approximate solution and the exact solution of (44) coincide for different values of , and and , , , and . The approximate solution of (44) with , , , and is displayed in Figure 7, while the absolute error of (44) with , , , and is displayed in Figure 8.
5. Conclusions and Future Works
We have constructed in this paper an efficient spectral-collocation algorithm for solving nonlinear time-delayed Burgers-type equations with positive power terms subject to initial and boundary conditions. The Jacobi-Gauss-Lobatto collocation methods based upon two general parameters, and , are developed and applied to the time-delayed Burgers-type equations for reformulating the problem to a system of second-order ordinary differential equations.
In fact, Jacobi-Gauss-Lobatto collocation method is slightly more complicated to implement than other orthogonal collocation methods but is more efficient and can be applied to a wider class of problems. We presented some advantages of the presented algorithm, as well as numerical results which demonstrate its accuracy and flexibility.
While the new algorithm presented in this paper only applied to a specific class of nonlinear time-delayed equations, it is very efficient and accurate whenever applicable. Moreover, high accuracy in long computational intervals and the stability of the proposed method encourage us to apply a similar scheme for the numerical solution of coupled nonlinear partial differential equations and other applied mathematics problems (see, e.g, [46–49] and references therein) in the future.
- C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang, Spectral Methods: Fundamentals in Single Domains, Scientific Computation, Springer, Berlin, Germany, 2006.
- C. I. Gheorghiu, Spectral Methods for Differential Problems, T. Popoviciu Institute of Numerical Analysis, Cluj-Napoca, Romaina, 2007.
- E. H. Doha, A. H. Bhrawy, and R. M. Hafez, “On shifted Jacobi spectral method for high-order multi-point boundary value problems,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 10, pp. 3802–3810, 2012.
- E. H. Doha, W. M. Abd-Elhameed, and A. H. Bhrawy, “New spectral-Galerkin algorithms for direct solution of high even-order differential equations using symmetric generalized Jacobi polynomials,” Collectanea Mathematica, 2013.
- K. Zhang, J. Li, and H. Song, “Collocation methods for nonlinear convolution Volterra integral equations with multiple proportional delays,” Applied Mathematics and Computation, vol. 218, no. 22, pp. 10848–10860, 2012.
- S. K. Vanani and F. Soleymani, “Tau approximate solution of weakly singular Volterra integral equations,” Mathematical and Computer Modelling, vol. 57, no. 3-4, pp. 494–502, 2013.
- L. Zhu and Q. Fan, “Solving fractional nonlinear Fredholm integro-differential equations by the second kind Chebyshev wavelet,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 6, pp. 2333–2341, 2012.
- E. H. Doha, A. H. Bhrawy, and S. S. Ezz-Eldien, “A Chebyshev spectral method based on operational matrix for initial and boundary value problems of fractional order,” Computers & Mathematics with Applications, vol. 62, no. 5, pp. 2364–2373, 2011.
- E. H. Doha, A. H. Bhrawy, and S. S. Ezz-Eldien, “A new Jacobi operational matrix: an application for solving fractional differential equations,” Applied Mathematical Modelling, vol. 36, no. 10, pp. 4931–4943, 2012.
- D. Rostamy, K. Karimi, L. Gharacheh, and M. Khaksarfard, “Spectral method for fractional quadratic Riccati differential equation,” Journal of Applied Mathematics and Bioinformatics, vol. 2, pp. 85–97, 2012.
- A. H. Bhrawy and A. S. Alofi, “A Jacobi-Gauss collocation method for solving nonlinear Lane-Emden type equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 1, pp. 62–70, 2012.
- A. H. Bhrawy and M. A. Alghamdi, “A shifted Jacobi-Gauss-Lobatto collocation method for solving nonlinear fractional Langevin equation involving two fractional orders in different intervals,” Boundary Value Problems, vol. 2012, article 62, 2012.
- K. Maleknejad and M. Attary, “A Chebysheve collocation method for the solution of higher-order Fredholm-Volterra integro-differential equations system,” Scientific Bulletin, Universitatea Politehnica din Bucuresti Series A, vol. 74, no. 4, pp. 17–28, 2012.
- M. Maleki, I. Hashim, M. T. Kajani, and S. Abbasbandy, “An adaptive pseudospectral method for fractional order boundary value problems,” Abstract and Applied Analysis, vol. 2012, Article ID 381708, 19 pages, 2012.
- G. Szegő, Orthogonal Polynomials, Colloquium Publications. XXIII. American Mathematical Society, 1939.
- E. H. Doha, A. H. Bhrawy, and R. M. Hafez, “A Jacobi dual-Petrov-Galerkin method for solving some odd-order ordinary differential equations,” Abstract and Applied Analysis, vol. 2011, Article ID 947230, 21 pages, 2011.
- A. H. Bhrawy, E. H. Doha, and R. M. Hafez, “A Jacobi dual-Petrov-Galerkin method for solving some odd-order ordinary differential equations,” Abstract and Applied Analysis, vol. 2011, Article ID 947230, 16 pages, 2011.
- H. Kim and R. Sakthivel, “Travelling wave solutions for time-delayed nonlinear evolution equations,” Applied Mathematics Letters, vol. 23, no. 5, pp. 527–532, 2010.
- M. Dehghan and R. Salehi, “Solution of a nonlinear time-delay model in biology via semi-analytical approaches,” Computer Physics Communications, vol. 181, no. 7, pp. 1255–1265, 2010.
- M. Ghasemi and M. T. Kajani, “Numerical solution of time-varying delay systems by Chebyshev wavelets,” Applied Mathematical Modelling, vol. 35, no. 11, pp. 5235–5244, 2011.
- X. T. Wang, “Numerical solution of delay systems containing inverse time by hybrid functions,” Applied Mathematics and Computation, vol. 173, no. 1, pp. 535–546, 2006.
- X. T. Wang, “Numerical solutions of optimal control for time delay systems by hybrid of block-pulse functions and Legendre polynomials,” Applied Mathematics and Computation, vol. 184, no. 2, pp. 849–856, 2007.
- X. T. Wang, “Numerical solutions of optimal control for linear time-varying systems with delays via hybrid functions,” Journal of the Franklin Institute, vol. 344, no. 7, pp. 941–953, 2007.
- S. Sedaghat, Y. Ordokhani, and M. Dehghan, “Numerical solution of the delay differential equations of pantograph type via Chebyshev polynomials,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 12, pp. 4815–4830, 2012.
- E. Tohidi, A. H. Bhrawy, and Kh. Erfani, “A collocation method based on Bernoulli operational matrix for numerical solution of generalized pantograph equation,” Applied Mathematical Modelling, vol. 37, pp. 4283–4294, 2013.
- I. Ali, H. Brunner, and T. Tang, “Spectral methods for pantograph-type differential and integral equations with multiple delays,” Frontiers of Mathematics in China, vol. 4, no. 1, pp. 49–61, 2009.
- D. Trif, “Direct operatorial tau method for pantograph-type equations,” Applied Mathematics and Computation, vol. 219, no. 4, pp. 2194–2203, 2012.
- R. Sakthivel, C. Chun, and J. Lee, “New travelling wave solutions of burgers equation with finite transport memory,” Zeitschrift fur Naturforschung A, vol. 65, no. 8, pp. 633–640, 2010.
- R. Sakthivel, “Robust stabilization the Korteweg-de Vries-Burgers equation by boundary control,” Nonlinear Dynamics, vol. 58, no. 4, pp. 739–744, 2009.
- K. Pandey, L. Verma, and A. K. Verma, “Du Fort-Frankel finite difference scheme for Burgers equation,” Arabian Journal of Mathematics, vol. 2, no. 1, pp. 91–101, 2013.
- Z.-Z. Sun and X.-N. Wu, “A difference scheme for Burgers equation in an unbounded domain,” Applied Mathematics and Computation, vol. 209, no. 2, pp. 285–304, 2009.
- R. C. Mittal and R. Jiwari, “A differential quadrature method for numerical solutions of Burgers-type equations,” International Journal of Numerical Methods for Heat & Fluid Flow, vol. 22, no. 6-7, pp. 880–895, 2012.
- D. Rostamy and K. Karimi, “Hypercomplex mathematics and HPM for the time-delayed Burgers equation with convergence analysis,” Numerical Algorithms, vol. 58, no. 1, pp. 85–101, 2011.
- A. G. Kudryavtsev and O. A. Sapozhnikov, “Determination of the exact solutions to the inhomogeneous burgers equation with the use of the darboux transformation,” Acoustical Physics, vol. 57, no. 3, pp. 311–319, 2011.
- H. Caglar and M. F. Ucar, “Non-polynomial spline method for the solution of non-linear Burgers equation,” Chaos and Complex Systems, pp. 213–218, 2013.
- M. Wang and F. Zhao, “Haar Wavelet method for solving two-dimensional Burgers' equation,” in Proceedings of the 2nd International Congress on Computer Applications and Computational Science, Advances in Intelligent and Soft Computing, vol. 145, pp. 381–387, 2012.
- J. Zhang, P. Wei, and M. Wang, “The investigation into the exact solutions of the generalized time-delayed Burgers-Fisher equation with positive fractional power terms,” Applied Mathematical Modelling. Simulation and Computation for Engineering and Environmental Systems, vol. 36, no. 5, pp. 2192–2196, 2012.
- S. Rendine, A. Piazza, and L. L. Cavalli-Sforza, “Simulation and separation by principle components of multiple demic expansions in Europe,” The American Naturalist, vol. 128, pp. 681–706, 1986.
- E. S. Fahmy, H. A. Abdusalam, and K. R. Raslan, “On the solutions of the time-delayed Burgers equation,” Nonlinear Analysis. Theory, Methods & Applications, vol. 69, no. 12, pp. 4775–4786, 2008.
- A. J. M. Jawad, M. D. Petković, and A. Biswas, “Soliton solutions of Burgers equations and perturbed Burgers equation,” Applied Mathematics and Computation, vol. 216, no. 11, pp. 3370–3377, 2010.
- H. Kim and R. Sakthivel, “Travelling wave solutions for time-delayed nonlinear evolution equations,” Applied Mathematics Letters, vol. 23, no. 5, pp. 527–532, 2010.
- E. H. Doha and A. H. Bhrawy, “An efficient direct solver for multidimensional elliptic Robin boundary value problems using a Legendre spectral-Galerkin method,” Computers & Mathematics with Applications, vol. 64, no. 4, pp. 558–571, 2012.
- E. H. Doha, W. M. Abd-Elhameed, and M. A. Bassuony, “New algorithms for solving high evenorder differential equations using third and fourth Chebyshev-Galerkin methods,” Journal of Computational Physics, vol. 236, pp. 563–579, 2013.
- E. H. Doha and A. H. Bhrawy, “A Jacobi spectral Galerkin method for the integrated forms of fourth-order elliptic differential equations,” Numerical Methods for Partial Differential Equations, vol. 25, no. 3, pp. 712–739, 2009.
- E. H. Doha, A. H. Bhrawy, D. Baleanu, and S. S. Ezz-Eldien, “On shifted Jacobi spectral approximations for solving fractional differential equations,” Applied Mathematics and Computation, vol. 219, no. 15, pp. 8042–8056, 2013.
- R. Naz, I. Naeem, and F. M. Mahomed, “First integrals for two linearly coupled nonlinear Duffing oscillators,” Mathematical Problems in Engineering, vol. 2011, Article ID 831647, 14 pages, 2011.
- M. Asadzadeh, D. Rostamy, and F. Zabihi, “Discontinuous Galerkin and multiscale variational schemes for a coupled damped nonlinear system of Schrodinger equations,” Numerical Methods for Partial Differential Equations, 2013.
- R. A. Van Gorder and K. Vajravelu, “A general class of coupled nonlinear differential equations arising in self-similar solutions of convective heat transfer problems,” Applied Mathematics and Computation, vol. 217, no. 2, pp. 460–465, 2010.
- M. Huang and Z. Zhou, “Standing wave solutions for the discrete coupled nonlinear Schrodinger equations with unbounded potentials,” Abstract and Applied Analysis, vol. 2013, Article ID 842594, 6 pages, 2013.