Research Article | Open Access
Chebyshev Collocation Method for Parabolic Partial Integrodifferential Equations
An efficient technique for solving parabolic partial integrodifferential equation is presented. This technique is based on Chebyshev polynomials and finite difference method. A priori error estimate for the proposed technique is deduced. Some examples are presented to illustrate the validity and efficiency of the presented method.
Partial integrodifferential equations (PIDEs) are the equations that combine partial differentiation and integration of the unknown function. They are used in modeling several phenomena where the effect of memory must be considered.
In this work, we consider the parabolic PIDE of the following formsubject to the following initial and boundary conditions:where , , , and are continuous functions.
This class of equations appears in various fields of physics and engineering such as heat conduction , compression of poroviscoelastic media , reaction diffusion problems , and nuclear reactor dynamics .
Chebyshev polynomials were proposed by the Russian mathematician Chebyshev. Many authors depend on Chebyshev collocation method to solve different types of equations such as linear differential equations , systems of high-order linear differential equations with variable coefficients , systems of high-order linear Fredholm-Volterra integrodifferential equations , fourth-order Sturm-Liouville problems , Troeschs problem , nonlinear differential equations , linear partial differential equations [19, 20], and nonlinear Fredholm-Volterra integrodifferential equations .
This paper is organized into six sections. Section 2 lists some notations and definitions of Chebyshev polynomials. Section 3 presents a description for the technique of solution to problems (1) and (2). In Section 4, we deduce a priori error estimate for the approximate solution. Section 5 is devoted to present some examples that illustrate the proposed technique of solution. Finally, Section 6 presents the conclusions of this study.
2. Fundamental Relations
On a general interval , the shifted Chebyshev polynomials take the following form:and the Chebyshev collocation points are given byThe function on is approximated using truncated shifted Chebyshev series in the form
The symbol sum with single prime indicates that the summation involves rather than . Similarly, the derivatives , , are written asThe function and its derivatives can be written in matrix form as where Matrix can be obtained from by the following relation : where where Then, can be written in the form
3. Description of Chebyshev Collocation Method
The partial time-derivative in (1) is discretized using finite difference method. Denote the time step by and the value of the function at time by . Then we have
The integral is handled numerically using the composite trapezoidal rule, and (1) is discretized into the following form:which is simplified towhere The approximate solution of (15) is computed using truncated Chebyshev series.
Theorem 1. If the assumed approximate solution of problem (15) is described by (5), then the discrete Chebyshev system is given by where are the Chebyshev collocation points. The fundamental matrix form of the discrete Chebyshev system is given bywhere The boundary conditions are integrated into system of (18) in the following form: The final form of the system is given byWe construct the two matrices and by replacing the first row and last row of the matrix by the corresponding row of boundary conditions.
4. Error Analysis
In (1), the time is discretized using finite difference method. From Taylor series, we have Approximation of integral term using composite trapezoidal integration rule leads to an error of order . So (15) can be written in the formEquation (23) shows that the truncation error due to time discretization of (1) is of order .
If is the exact solution of (1), then is the exact solution of (23). Let and be the Chebyshev series solution of (15) and (23), respectively. If , then the discrete Chebyshev system is given byBy subtracting (18) from (24), the following relation is obtained:Let be the Lagrange interpolating polynomial of order of on the grid of Chebyshev collocation points. Then is given bywhere
Theorem 2 (see ). If , then for any grid of nodes and for any the interpolation error at is given bywhere and are the nodal polynomial given byLet . Then is the solution of the following equation:where can be written as . So the discrete Chebyshev system of (30) is given byby subtracting (18) from (32), we obtain
Proof. Adding and subtracting the polynomial yield The upper bound for the second term can be found using some properties of normThe upper bound for the second term can be found by adding and subtracting the polynomial following the work in . Summing up these two upper bounds yields (34).
5. Numerical Examples
Example 1 (see ). Consider the parabolic PIDE (1) with the initial condition and with boundary conditions and with , , , and The exact solution of this example is , where and .
The maximum absolute error is tabulated in Table 1 for Chebyshev polynomial together with the results of  at and . The graphs of exact and approximate solutions for are illustrated in Figure 1.
The maximum absolute error is tabulated in Table 2 for Chebyshev polynomial at and .
Example 2. Consider the parabolic PIDE (1) with the initial condition and with boundary conditions and with , , , and . The exact solution of this example is given by
The maximum absolute error is tabulated in Table 3 for Chebyshev polynomial at different times. The graphs of exact and approximate solutions for and are illustrated in Figure 2.
Example 3 (see ). Consider the parabolic PIDE (1) with the initial condition and with boundary conditions and with , , , and which can be chosen such that the exact solution in this example is .
The maximum absolute error is tabulated in Table 4 for Chebyshev polynomial together with the results of . The graphs of exact and approximate solutions for are illustrated in Figure 3.
Example 4. Consider the parabolic PIDE (1) with the initial condition and with boundary conditions and with , , , and which can be chosen such that the exact solution in this example is .
The maximum absolute error is tabulated in Table 5 for Chebyshev polynomial at different values of and . The graphs of exact and approximate solutions for , at , and are illustrated in Figures 4 and 5, respectively.
Chebyshev collocation method is successfully used for solving parabolic partial integrodifferential equation. This method reduced the considered problem into linear system of algebraic equations that can be solved successively to obtain a numerical solution at varied time levels. Numerical examples show that the results of above scheme are in good agreement with the exact ones. Comparisons with the results obtained by using radial basis functions and backward-Euler scheme show that the Chebyshev method yields good results with fewer number of iterations. The method also renders good results for problems where solution exhibits fast changes but with larger number of iterations as anticipated by the error estimate. Moreover, the above scheme can be developed to solve nonlinear parabolic partial integrodifferential equation.
The authors declare that there are no competing interests regarding the publication of this paper.
- M. Dehghan and F. Shakeri, “Solution of parabolic integro-differential equations arising in heat conduction in materials with memory via He's variational iteration technique,” International Journal for Numerical Methods in Biomedical Engineering, vol. 26, no. 6, pp. 705–715, 2010.
- G. J. Habetler and R. L. Schiffman, “A finite difference method for analyzing the compression of poro-viscoelastic media,” Computing, vol. 6, pp. 342–348, 1970.
- H. Engler, “On some parabolic integro-differential equations: existence and asymptotics of solutions,” in Equadiff 82, H. W. Knobloch and K. Schmitt, Eds., vol. 1017 of Lecture Notes in Mathematics, pp. 161–167, 1983.
- C. V. Pao, “Solution of a nonlinear integrodifferential system arising in nuclear reactor dynamics,” Journal of Mathematical Analysis and Applications, vol. 48, pp. 470–492, 1974.
- F. Fakhar-Izadi and M. Dehghan, “The spectral methods for parabolic Volterra integro-differential equations,” Journal of Computational and Applied Mathematics, vol. 235, no. 14, pp. 4032–4046, 2011.
- C. E. Greenwell-Yanik and G. Fairweather, “Analyses of spline collocation methods for parabolic and hyperbolic problems in two space variables,” SIAM Journal on Numerical Analysis, vol. 23, no. 2, pp. 282–296, 1986.
- J.-P. Kauthen, “The method of lines for parabolic partial integro-differential equations,” Journal of Integral Equations and Applications, vol. 4, no. 1, pp. 69–81, 1992.
- S. Larsson, V. Thomée, and L. B. Wahlbin, “Numerical solution of parabolic integro-differential equations by the discontinuous Galerkin method,” Mathematics of Computation, vol. 67, no. 221, pp. 45–71, 1998.
- J. Ma, “Finite element methods for partial Volterra integro-differential equations on two-dimensional unbounded spatial domains,” Applied Mathematics and Computation, vol. 186, no. 1, pp. 598–609, 2007.
- Z. Avazzadeh, Z. Beygi Rizi, F. M. Maalek Ghaini, and G. B. Loghmani, “A numerical solution of nonlinear parabolic-type Volterra partial integro-differential equations using radial basis functions,” Engineering Analysis with Boundary Elements, vol. 36, no. 5, pp. 881–893, 2012.
- A. F. Soliman, A. M. El-asyed, and M. S. El-Azab, “On the numerical solution of partial integro-differential equations,” Mathematical Sciences Letters, vol. 1, no. 1, pp. 71–80, 2012.
- A. Filiz, “Numerical Solution of Parabolic Volterra integro-differential equations via backward Euler scheme,” Journal of Computational and Applied Mathematics, vol. 3, no. 6, pp. 277–282, 2013.
- M. Sezer and M. Kaynak, “Chebyshev polynomial solutions of linear differential equations,” International Journal of Mathematical Education in Science and Technology, vol. 27, no. 4, pp. 607–611, 1996.
- A. Akyüz and M. Sezer, “Chebyshev polynomial solutions of systems of high-order linear differential equations with variable coefficients,” Applied Mathematics and Computation, vol. 144, no. 2-3, pp. 237–247, 2003.
- A. Akyüz and M. Sezer, “A Chebyshev collocation method for the solution of linear integro-differential equations,” International Journal of Computer Mathematics, vol. 72, no. 4, pp. 491–507, 1999.
- M. El-Gamel and M. Sameeh, “An efficient technique for finding the eigenvalues of fourth-order Sturm-Liouville problems,” Applied Mathematics, vol. 3, no. 8, pp. 920–925, 2012.
- M. El-Gamel and M. Sameeh, “A Chebyshev collocation method for solving Troesch's problem,” International Journal of Mathematics and Computer Applications Research, vol. 3, pp. 23–32, 2013.
- A. Akyüz-Daşcıoğlu and H. Çerdik-Yaslan, “The solution of high-order nonlinear ordinary differential equations by Chebyshev series,” Applied Mathematics and Computation, vol. 217, no. 12, pp. 5658–5666, 2011.
- C. Keşan, “Chebyshev polynomial solutions of second-order linear partial differential equations,” Applied Mathematics and Computation, vol. 134, no. 1, pp. 109–124, 2003.
- G. Yuksel, O. R. Isik, and M. Sezer, “Error analysis of the Chebyshev collocation method for linear second-order partial differential equations,” International Journal of Computer Mathematics, vol. 92, no. 10, pp. 2121–2138, 2015.
- H. Çerdk-Yaslan and A. Akyüz-Yaslan, “Chebyshev polynomial solution of nonlinear Fredholm Volterra integro-differential equations,” Journal of Science and Arts, vol. 6, pp. 89–101, 2006.
- A. Quarteroni, R. Sacco, and F. Saleri, Numerical Mathematics, vol. 37 of Texts in Applied Mathematics, Springer, Berlin, Germany, 2nd edition, 2007.
Copyright © 2016 M. Sameeh and A. Elsaid. 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.