0
Research Papers

# An Accurate Jacobi Pseudospectral Algorithm for Parabolic Partial Differential Equations With Nonlocal Boundary Conditions

[+] Author and Article Information
E. H. Doha

Department of Mathematics,
Faculty of Science,
Cairo University,
Giza 12613, Egypt
e-mail: eiddoha@frcu.eun.eg

A. H. Bhrawy

Department of Mathematics,
Faculty of Science,
King Abdulaziz University,
Jeddah 21589, Saudi Arabia
Department of Mathematics,
Faculty of Science,
Beni-Suef University,
Beni-Suef 62511, Egypt
e-mail: alibhrawy@yahoo.co.uk

M. A. Abdelkawy

Department of Mathematics,
Faculty of Science,
Beni-Suef University,
Beni-Suef 62511, Egypt
e-mail: melkawy@yahoo.com

Contributed by the Design Engineering Division of ASME for publication in the Journal of Computational and Nonlinear Dynamics. Manuscript received December 30, 2013; final manuscript received February 20, 2014; published online January 12, 2015. Assoc. Editor: Ahmet S. Yigit.

J. Comput. Nonlinear Dynam 10(2), 021016 (Mar 01, 2015) (13 pages) Paper No: CND-13-1330; doi: 10.1115/1.4026930 History: Received December 30, 2013; Revised February 20, 2014; Online January 12, 2015

## Abstract

A new spectral Jacobi–Gauss–Lobatto collocation (J–GL–C) method is developed and analyzed to solve numerically parabolic partial differential equations (PPDEs) subject to initial and nonlocal boundary conditions. The method depends basically on the fact that an expansion in a series of Jacobi polynomials $Jn(θ,ϑ)(x)$ is assumed, for the function and its space derivatives occurring in the partial differential equation (PDE), the expansion coefficients are then determined by reducing the PDE with its boundary conditions into a system of ordinary differential equations (SODEs) for these coefficients. This system may be solved numerically in a step-by-step manner by using implicit the Runge–Kutta (IRK) method of order four. The proposed method, in contrast to common finite-difference and finite-element methods, has the exponential rate of convergence for the spatial discretizations. Numerical results indicating the high accuracy and effectiveness of this algorithm are presented.

<>

## References

Canuto, C., Hussaini, M. Y., Quarteroni, A., and Zang, T. A., 2006, Spectral Methods: Fundamentals in Single Domains, Springer-Verlag, New York.
Dehghan, M., and Sabouri, M., 2013, “A Legendre Spectral Element Method on a Large Spatial Domain to Solve the Predator-Prey System Modeling Interacting Populations,” Appl. Math. Model., 37, pp. 1028–1038.
Doha, E. H., Bhrawy, A. H., Baleanu, D., and Hafez, R. M., 2014, “A New Jacobi Rational-Gauss Collocation Method for Numerical Solution of Generalized Pantograph Equations,” Appl. Numer. Math., 77, pp. 43–54.
Bhrawy, A. H., 2013, “A Jacobi-Gauss-Lobatto Collocation Method for Solving Generalized Fitzhugh-Nagumo Equation With Time-Dependent Coefficients,” Appl. Math. Comput., 222, pp. 255–264.
Wan, Z., Chen, Y., and Huang, Y., 2009, “Legendre Spectral Galerkin Method for Second-Kind Volterra Integral Equations,” Front. Math. China, 4, pp. 181–193.
Rahmoune, A., 2013, “Spectral Collocation Method for Solving Fredholm Integral Equations on the Half-Line,” Appl. Math. Comput., 219, pp. 9254–9260.
Ahmadian, A., Suleiman, M., Salahshour, S., and Baleanu, D., 2013, “A Jacobi Operational Matrix for Solving a Fuzzy Linear Fractional Differential Equation,” Adv. Difference Equations, 2013, p. 104.
Bhrawy, A. H., and Alghamdi, M., 2012, “A Shifted Jacobi-Gauss-Lobatto Collocation Method for Solving Nonlinear Fractional Langevin Equation Involving Two Fractional Orders in Different Intervals,” Boundary Value Problems, 2012, p. 62.
Doha, E. H., Bhrawy, A. H., Hafez, R. M., and Abdelkawy, M. A., 2014, “A Chebyshev-Gauss-Radau Scheme for Nonlinear Hyperbolic System of First Order,” Appl. Math. Info. Sci., 8(2), pp. 535–544.
Bhrawy, A. H., and Alofi, A. S., 2012, “A Jacobi-Gauss Collocation Method for Solving Nonlinear Lane-Emden Type Equations,” Commun. Nonlin. Sci. Numer. Simul., 17, pp. 62–70.
Bialecki, B., and Karageorghis, A., 2008, “Spectral Chebyshev-Fourier Collocation for the Helmholtz and Variable Coefficient Equations in a Disk,” J. Comput. Phys., 227(19), pp. 8588–8603.
Yuzbasi, S., Sezer, M., and Kemanci, B., 2013, “Numerical Solutions of Integro-differential Equations and Application of a Population Model With an Improved Legendre Method,” Appl. Math. Model., 37, pp. 2086–2101.
Yuzbasi, S., and Sahin, N., 2013, “Numerical Solutions of Singularly Perturbed One-Dimensional Parabolic Convection. Diffusion Problems by the Bessel Collocation Method,” Appl. Math. Comput., 220, pp. 305–315.
Esmaeili, S., Shamsi, M., and Luchko, Y., 2011, “Numerical Solution of Fractional Differential Equations With a Collocation Method Based on Mntz Polynomials,” Comput. Math. Appl., 62(3), pp. 918–929.
Doha, E. H., Bhrawy, A. H., Abdelkawy, M. A., and Hafez, R. M., 2014, “A Jacobi Collocation Approximation for Nonlinear Coupled Viscous Burgers' Equation,” Cent. Eur. J. Phys., 12, pp. 111–122.
Eslahchi, M. R., Dehghan, M., and Parvizi, M., 2014, “Application of the Collocation Method for Solving Nonlinear Fractional Integro-differential Equations,” J. Comput. Appl. Math., 257, pp. 105–128.
Capsso, V., and Kunisch, K., 1988, “A Reaction-Diffusion System Arising in Modeling Man-Environment Diseases,” Q. Appl. Math., 48, pp. 431–449.
Yurchuk, N. I., 1986, “Mixed Problem With an Integral Condition for Certain Parabolic Equations,” J. Differential Equations, 22, pp. 2117–2126.
Cannon, J. R., Lin, Y., and Matheson, A. L., 1993, “The Solution of the Diffusion Equation in Two-Space Variables Subject to the Specification of Mass,” Applicable Anal., 50, pp. 1–19.
Bouziani, A., 1999, “On a Class of Parabolic Equations With a Nonlocal Boundary Condition,” Acad. R. Belg., Bull. Classe Sci., 10, pp. 61–77.
Diaz, J. I., 1997, “On a Nonlocal Elliptic Problem Arising in the Magnetic Confinement of a Plasma in a Stellarator,” Nonlin. Anal., 30, pp. 3963–3974.
Farengo, R., Lee, Y. C., and Guzdar, P. N., 1983, “An Electromagnetic Integral Equation: Application to Microtearing Modes,” Phys. Fluids, 26, pp. 3515–3523.
Ashyralyev, A., and Aggez, N., 2004, “A Note on the Difference Schemes of the Nonlocal Boundary Value Problems for Hyperbolic Equations,” Numer. Funct. Anal. Opt., 25, pp. 439–462.
Sajavičius, S.-U., 2012, “Stability of the Weighted Splitting Finite-Difference Scheme for a Two-Dimensional Parabolic Equation With Two Nonlocal Integral Conditions,” Comput. Math. Appl., 64, pp. 3485–3499.
Martín-Vaquero, J., and Vigo-Aguiar, J., 2009, “On the Numerical Solution of the Heat Conduction Equations Subject to Nonlocal Conditions,” Appl. Numer. Math., 59, pp. 2507–2514.
Dehghan, M., 2004, “Numerical Schemes for One-Dimensional Parabolic Equations With Nonstandard Initial Condition,” Appl. Math. Comput., 147, pp. 321–331.
Golbabai, A., and Javidi, M., 2007, “Numerical Solution for Non-classical Parabolic Problem Based on Chebyshev Spectral Collocation Method,” Appl. Math. Comput., 190, pp. 179–185.
Liu, Y., 1999, “Numerical Solution of the Heat Equation With Nonlocal Boundary Conditions,” J. Comput. Appl. Math., 110, pp. 115–127.
Olmstead, W. E., and Roberts, C. A., 1997, “The One-Dimensional Heat Equation With a Nonlocal Initial Condition,” Appl. Math. Lett., 10, pp. 89–94.
Mu, L., and Du, H., 2008, “The Solution of a Parabolic Differential Equation With Non-local Boundary Conditions in the Reproducing Kernel Space,” Appl. Math. Comput., 202, pp. 708–714.
Ang, W.-T., 2006, “Numerical Solution of a Non-classical Parabolic Problem: An Integro-differential Approach,” Appl. Math. Comput., 175, pp. 969–979.
Lardner, R. W., 1990, “Stability of the Numerical Solution of a Parabolic System With Integral Subsidiary Conditions,” Comput. Math. Appl., 19, pp. 41–46.
Slodička, M., and Dehilis, S., 2009, “A Numerical Approach for a Semilinear Parabolic Equation With a Nonlocal Boundary Condition,” J. Comput. Appl. Math., 231, pp. 715–724.
Doha, E. H., Bhrawy, A. H., and Hafez, R. M., 2011, “A Jacobi-Jacobi Dual-Petrov-Galerkin Method for Third- and Fifth-Order Differential Equations,” Math. Comput. Model., 53, pp. 1820–1932.
Tchiotsop, D., Wolf, D., Louis-Dorr, V., and Husson, R., 2007, “ECG Data Compression Using Jacobi Polynomials,” Proceedings of the 29th Annual International Conference of the IEEE EMBS, pp. 1863–1867.
Gottlieb, D., and Shu, C.-W., 1997, “On the Gibbs Phenomenon and Its Resolution,” SIAM Rev., 29, pp. 644–668.
Doha, E. H., and Bhrawy, A. H., 2009, “A Jacobi Spectral Galerkin Method for the Integrated Forms of Fourth-Order Elliptic Differential Equations,” Numer. Methods Partial Differential Equations, 25(3), pp. 712–739.
El-Kady, M., 2012, “Jacobi Discrete Approximation for Solving Optimal Control Problems,” J. Korean Math. Soc., 49, pp. 99–112.
Szegö, G., 1939, Orthogonal Polynomials, Colloquium Publications, XXIII, American Mathematical Society, MR 0372517G, Providence, RI.
Doha, E. H., Bhrawy, A. H., Abdelkawy, M. A., and Van Gorder, R. A., 2014, “Jacobi-Gauss-Lobatto Collocation Method for the Numerical Solution of 1+1 Nonlinear Schrodinger Equations,” J. Comput. Phys., 261, pp. 244–255.
Butcher, J. C., 1987, The Numerical Analysis of Ordinary Differential Equations, Runge-Kutta and General Linear Methods, Wiley, New York.
Butcher, J. C., 1964, “Implicit Runge-Kutta Processes,” Math. Comput., 18, pp. 50–64.
Doha, E. H., 1983, “On the Chebyshev Method to Linear Parabolic P. D. E. s With Inhomogeneous Mixed Boundary Conditions,” Arab J. Math., 4, pp. 31–57.
Mittal, R. C., and Jain, R. K., 2012, “Cubic B-Splines Collocation Method for Solving Nonlinear Parabolic Partial Differential Equations With Neumann Boundary Conditions,” Commun. Nonlin. Sci. Numer. Simul., 17, pp. 4616–4625.
Dehghan, M., 2004, “Numerical Solution of a Parabolic Equation Subject to Specification of Energy,” Appl. Math. Comput., 149, pp. 31–45.
Zhoua, Y., Cui, M., and Lin, Y., 2009, “Numerical Algorithm for Parabolic Problems With Non-classical Conditions,” J. Comput. Appl. Math., 230, pp. 770–780.

## Figures

Fig. 1

The numerical solution u˜(x,t) of the problem in Eq. (4.1) where θ = ϑ = 0 and N = 12

Fig. 2

The absolute error of the problem in Eq. (4.1), where θ = ϑ = 0 and N = 12

Fig. 3

The curves of numerical and exact solutions of the problem in Eq. (4.1) for the different values of t = 0.5, 1.0, and 1.5 where θ = ϑ = 0 and N = 12

Fig. 4

The curves of numerical and exact solutions of the problem in Eq. (4.1) for the different values of x = 0.5, 1.0, and 1.5 where θ = ϑ = 0 and N = 12

Fig. 5

The curve of absolute error of the problem in Eq. (4.1) at x = 0, where θ = ϑ = 0 and N = 12

Fig. 6

The numerical solution u˜(x,t) of the problem in Eq. (4.7), where -θ = ϑ = 1/2 and N = 10

Fig. 7

The curves of numerical and exact solutions of the problem in Eq. (4.7) for the different values of t = 0.1, 0.5, and 0.9 where -θ = ϑ = 1/2 and N = 10

Fig. 8

The curves of numerical and exact solutions of the problem in Eq. (4.7) for the different values of x = 0.1, 0.2, and 0.3 where θ = ϑ = 0,ɛ = 2-4 and N = 16

Fig. 9

The absolute error of the problem in Eq. (4.7), at x = 0, where -θ = ϑ = 1/2 and N = 10

Fig. 10

The absolute error of the problem in Eq. (4.7) at x = 0, where -θ = ϑ = 1/2 and N = 10

Fig. 11

The numerical solution u˜(x,t) of the problem in Eq. (4.12), where θ = ϑ = -1/2 and N = 8

Fig. 12

The absolute error of the problem in Eq. (4.12), where θ = ϑ = -1/2 and N = 8

Fig. 13

The curves of numerical and exact solutions of the problem in Eq. (4.12) for the different values of t = 1.0, 1.5, and 2.0 where θ = ϑ = -1/2 and N = 8

Fig. 14

The curves of numerical and exact solutions of the problem in Eq. (4.12) for the different values of x = 0.1, 0.2, and 0.3 where θ = ϑ = -1/2 and N = 8

Fig. 15

The absolute error of the problem in Eq. (4.12) at x = 0, where θ = ϑ = -1/2 and N = 8

Fig. 16

The absolute error of the problem in Eq. (4.12) at t = 0, where θ = ϑ = -1/2 and N = 8

Fig. 17

The numerical solution u˜(x,t) of the problem in Eq. (4.15), where θ = ϑ = 1/2 and N = 8

Fig. 18

The absolute error of the problem in Eq. (4.15), where θ = ϑ = 1/2 and N = 8

Fig. 19

The curves of numerical and exact solutions of the problem in Eq. (4.15) for the different values of t = 0.1, 0.2, and 0.3 where θ = ϑ = 1/2 and N = 8

Fig. 20

The curves of numerical and exact solutions of the problem in Eq. (4.15) for the different values of t = 0.1, 0.2, and 0.3 where θ = ϑ = 1/2 and N = 8

Fig. 21

The numerical solution u˜(x,t) of the problem in Eq. (4.18), where θ = ϑ = -1/2 and N = 10

Fig. 22

The curves of numerical and exact solutions of the problem in Eq. (4.18) for the different values of t = 0.1, 0.5, and 1.0 where θ = ϑ = -1/2 and N = 10

Fig. 23

The curves of numerical and exact solutions of the problem in Eq. (4.18) for the different values of x = 0.1, 0.2, and 0.3 where θ = ϑ = -1/2 and N = 10

Fig. 24

The absolute error of the problem in Eq. (4.18) at x = 0, where θ = ϑ = -1/2 and N = 10

## Discussions

Some tools below are only available to our subscribers or users with an online account.

### Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related Proceedings Articles
Related eBook Content
Topic Collections