- About this Journal ·
- Abstracting and Indexing ·
- 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 2014 (2014), Article ID 849682, 11 pages
Numerical Method Using Cubic Trigonometric B-Spline Technique for Nonclassical Diffusion Problems
1Department of Mathematics, University of Sargodha, Sargodha 40100, Pakistan
2School of Mathematical Sciences, Universiti Sains Malaysia (USM), 11800 Penang, Malaysia
3Department of Mathematics, Gomal University, Dera Ismail Khan 29050, Pakistan
Received 17 October 2013; Accepted 13 April 2014; Published 20 May 2014
Academic Editor: Sining Zheng
Copyright © 2014 Muhammad Abbas 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.
A new two-time level implicit technique based on cubic trigonometric B-spline is proposed for the approximate solution of a nonclassical diffusion problem with nonlocal boundary constraints. The standard finite difference approach is applied to discretize the time derivative while cubic trigonometric B-spline is utilized as an interpolating function in the space dimension. The technique is shown to be unconditionally stable using the von Neumann method. Several numerical examples are discussed to exhibit the feasibility and capability of the technique. The and error norms are also computed at different times for different space size steps to assess the performance of the proposed technique. The technique requires smaller computational time than several other methods and the numerical results are found to be in good agreement with known solutions and with existing schemes in the literature.
This study deals with the numerical solution of a nonclassical diffusion problem with two nonlocal boundary constraints using cubic trigonometric B-splines. This problem arises in several branches of science. In particular, electrochemistry , heat conduction process , thermoelasticity , plasma physics , semiconductor modeling , biotechnology , control theory, and inverse problems . The analysis, development, and implementation of numerical methods for the solution of such diffusion problems have received wide attention in the literature.
Consider an insulated rod of length located on the -axis of the interval . Let the rod have a source of heat. Let denote the temperature in the insulated rod with ends held at constant temperatures and , and the initial temperature distribution along the rod is . The problem is to study the flow of heat in the rod and in this paper the partial differential equation governing the flow of heat in the rod is given by the diffusion equation with specification of energy: with the initial constraint and the nonlocal boundary constraints where , are known constants, are known continuous functions, is the thermal diffusivity of the rod, and are prescribed functions. This problem has been studied by Dehghan , Martín-Vaquero and Vigo-Aguiar , Li and Wu , and Golbabai and Javidi . Several physical circumstances might be modeled by equation and constraints (1)–(3) and several examples of application in physics with comprehensive derivations of the above mentioned problem can be found in [3, 12–14].
There are several numerical methods in the literature that have been developed for solving the proposed problem (1) subject to initial and nonlocal boundary constraints (2)-(3). The methods were based, for instance, on the forward Euler method, the backward Euler approach or the Crank-Nicolson scheme [15, 16], Laplace transformation , and so forth. Dehghan  presented four finite difference approaches, namely, the BTCS (backward time centred space) scheme, the implicit Crandall’s formula, the 3-point FTCS (forward time centred space) two-level scheme, and the Dufort-Frankel three-level approach for the numerical solution of parabolic equation with nonlocal specification. Martín-Vaquero and Vigo-Aguiar  improved the order of convergence of the implicit Crandall’s formula proposed by Dehghan  and also improved the accuracy of the method. Martín-Vaquero and Vigo-Aguiar  developed an algorithm for the solution of the heat conduction equations with nonlocal constraints which reduced the CPU time and enhanced the accuracy of (3, 3) Crandall’s formula proposed in . Li and Wu  proposed an algorithm which was based on the transverse method of lines (TMOL) which can reduce a nonclassical diffusion equation to a series of ordinary differential equations (ODEs). Subsequently, the authors in  used an analytic reproducing kernel technique to solve ODEs with integral boundary constraints. Dehghan [19, 20], jointly with Tatari and Dehghan , has proposed several efficient techniques for the numerical solutions of partial differential equations subject to nonlocal boundary constraints. Golbabai and Javidi  introduced a Chebyshev spectral collocation method (CSCM) based on Chebyshev polynomials for solving a parabolic problem subject to nonlocal boundary constraints. For more details on other numerical methods for the solution of a one-dimensional heat equation subject to nonlocal boundary constraints in the literature, see [22–31].
The study of B-spline functions is a key element in computer-aided geometric design [32–35]. It has also attracted attention in the literature [36–51] to the numerical solution of various differential equations [38–40]. This is because they have important geometric properties and features that make them amenable to more detailed analysis. Numerical methods based on B-spline functions of various degrees have been utilized for solving initial and boundary value problems. As examples, a cubic B-spline collocation method was used to solve a nonlinear diffusion equation subject to certain initial and Dirichlet boundary constraints , a finite element method based on bivariate splines has been used for solving parabolic partial differential equation , and the combination of finite difference approach and cubic B-spline method was applied for the solution of a one-dimensional heat equation subject to local boundary constraints [43, 44]. Goh et al.  presented a comparison of cubic B-spline and extended cubic uniform B-spline based collocation methods for solving a one-dimensional heat equation with a nonlocal initial constraint and concluded that extended cubic uniform B-spline with an appropriate value of parameters gives better results than the cubic B-spline. A finite difference scheme based on cubic B-spline was also used for solving the one-dimensional wave equation , advection-diffusion equation , one-dimensional coupled viscous Burgers’ equation , system of strongly coupled reaction-diffusion equations , and one-dimensional hyperbolic problems .
In our present paper, a new two-time level implicit technique is developed to approximate the solution of the nonclassical diffusion problem (1) subject to initial constraints in (2) and nonlocal boundary constraints in (1)–(3). The technique is based on the cubic trigonometric B-spline functions. A finite difference approach and -weighted scheme are applied for the time and space discretization, respectively. Some researchers have considered the ordinary B-spline collocation method for solving the heat equation subject to local and nonlocal boundary constraints but, so far as we are aware, not with the cubic trigonometric B-spline collocation method. Cubic trigonometric B-spline is used as an interpolating function in the space dimension. The unconditional stability property of the method is proved by von Neumann method. The feasibility of the method is shown by test problems with , instead of smaller time step size and the approximated solutions are found to be in good agreement with the known exact solutions.
The outline of this study is as follows. A numerical solution of nonclassical diffusion problem is presented in Section 2. In Section 3, the cubic trigonometric B-spline is utilized as an interpolating function in the space dimension. The von Neumann approach is used to prove the stability of the method in Section 4. Numerical examples are considered in Section 5 to show the achievability of the proposed method. Finally, the concluding remarks of this study are given in Section 6.
2. Solution of Nonclassical Diffusion Problem
Consider a uniform mesh with grid points to discretize the grid region with , and , , . Here the quantities and are mesh space size and time step size, respectively. The time derivative can be approximated by using the standard finite difference formula: Using the approximation of (4), (1) becomes Using -weighted technique, the space derivatives of (5) can be written as where and the subscripts and are successive time levels. It is noted that the system becomes an explicit scheme when , a fully implicit scheme when , and a Crank-Nicolson scheme when [43, 49]. In this paper, we use the Crank-Nicolson approach. Hence, (6) becomes After simplification, (7) leads to The space derivatives are approximated by using cubic trigonometric B-spline and are discussed in the next section.
3. Cubic Trigonometric B-Spline Technique
In this section, we discuss the cubic trigonometric B-spline collocation method (CuTBSM) for the numerical solution of the nonclassical diffusion equation (1). Consider a mesh which is equally divided by knots into subintervals , where . Our approach for the nonclassical diffusion equation using collocation method with cubic trigonometric B-spline is to seek an approximate solution as  where are to be determined for the approximated solutions to the exact solutions , at the point . are twice continuously differentiable piecewise cubic trigonometric B-spline basis functions over the mesh defined by [49–51] where and where . The approximations at the point over subinterval can be defined as In order to obtain the approximations to the solutions, the values of and its derivatives at nodal points are required and these derivatives are tabulated in Table 1, where Using approximate functions (10) and (12) and following Mittal and Arora , the values at the knots of and their derivatives up to second order are determined in terms of time parameters as Substituting (12) into (8) gives the following equation: The system thus obtained on simplifying (15) consists of linear equations in unknowns at the time level . Equation (9) is applied to the boundary constraints (2) and (3) for two additional linear equations to obtain a unique solution of the resulting system: From (15), (16), and (17), the system can be written in the matrix vector form as follows: where and and are -dimensional matrix given by where Thus, the system (18) becomes a matrix system of dimension which is a tridiagonal system that can be solved by the Thomas Algorithm .
3.1. Initial State Vector
After the initial vectors have been computed from the initial constraints, the approximate solutions at a particular time level can be calculated repeatedly by solving the recurrence relation (15) .
The initial vectors can be obtained from the initial condition and boundary values of the derivatives of the initial condition as follows [40, 49]: Thus (22) yields a matrix system, of the form where The solution of this system can be found by the use of the Thomas Algorithm.
In this section, the von Neumann stability method is applied for investigating the stability of the proposed scheme. This approach has been used by many researchers [18, 40, 44, 45, 47–49]. Substituting the approximate solution and its derivatives at knots with , into (6) yields a difference equation with variables given by Substituting the values of into (25) we obtain where Simplifying it leads to where Now on inserting the trial solutions (one Fourier mode out of the full solution) at a given point into (28) and rearranging the equations, is the mode number, is the element size, and , we get Dividing (30) by and rearranging the equation, we get Let Therefore, (31) can be written as Equation (33) can be rewritten as where For stability, the maximum modulus of the eigenvalues of the matrix has to be less than or equal to one . Since , , and , we always have Thus, from (36), the proposed scheme for nonclassical diffusion equation (with term, ) is unconditionally stable and it is also unconditionally stable with a general term , since the modulus of the eigenvalues must be less than one . We recall Duhamel’s principle (, chapter 9); a scheme is stable for equation if it is stable for the equation . This means that there are no constraints on grid size and step size in time level , but we should prefer those values of and for which we obtain the best accuracy of the scheme.
5. Results and Discussions
In this section, the cubic trigonometric B-spline collocation method is employed to obtain the numerical solutions for one-dimensional nonclassical diffusion problem with nonlocal boundary constraints given in (1)–(3). Two numerical examples are discussed in this section to exhibit the capability and efficiency of the proposed trigonometric spline method. Numerical results are compared with existing methods in the literature and with the exact solution at the different nodal points for some time levels using some particular space step size and time step . In order to calculate the maximum errors and relative error norms of the proposed method numerically, we use the following formulas:
Example 1. Consider the nonclassical diffusion problem ((1)–(3)), with This test problem is from Dehghan  and Li and Wu  and the known solution is . We compare the maximum errors with TMOL  when they are considered with space size and errors are recorded with several values of time step , given in Table 2 and also shown in Figure 1. It is worth noting that the results obtained using CuTBSM are more accurate as compared to TMOL . We also compare the relative errors of numerical value of with different space step and with time step size , , instead of which was used in the BTCS , the implicit Crandall’s formula , the 3-point FTCS , and the Dufort-Frankel three-level approach  and TMOL , and these results are tabulated in Table 3. It is clearly shown from this table that the obtained results by using CuTBSM are more precise as compared to methods in [8, 10]. Figure 2 shows the approximate solution and exact solution for this example at different time levels with and .
Example 2. We consider another numerical test problem, with The exact solution of this equation is and this test problem has been taken from [8, 10, 11]. This problem is tested using different values of and to show the capability of the present method for solving nonclassical diffusion equation ((1)–(3)). The final time is taken . The maximum errors of the numerical method are calculated at different time levels with different time step size and it is observed that they are more accurate as compared to TMOL  and Chebyshev spectral collocation method (CSCM) based on Chebyshev polynomials . The numerical errors are tabulated in Tables 4 and 5 and are also depicted graphically in Figure 3. The relative errors of numerical value with different space step and with time step size , , instead of which was used in the BTCS , the implicit Crandall’s formula , the 3-point FTCS , and the Dufort-Frankel three-level approach  and TMOL  and they are recorded in Table 6. It is worth noting that numerical results are much better than the methods in [8, 10, 11]. A comparison of numerical solutions at different time levels with known solution is presented graphically in the Figure 4.
Example 3. Finally, we consider nonclassical diffusion problem ((1)-(3)), with This test problem has been taken from  and its exact solution is . The final time is taken as . The maximum errors of the proposed scheme are considered at with different time step sizes that are depicted graphically in Figure 5. The relative errors of numerical value are calculated with different space size step with , , and they are given in Table 7. It is worth noting that the numerical results are found to be in good agreement with exact solutions. A comparison of numerical solutions at different time levels with known solution is presented graphically in Figure 6.
6. Concluding Remarks
In this paper, a new two-time level implicit scheme based on cubic trigonometric B-spline has been used to solve the nonclassical diffusion problem with known initial and with nonlocal boundary constraints instead of the usual boundary constraints. A usual finite difference discretization is used for time derivatives and cubic trigonometric B-spline is applied for space derivatives. It is noted that the accuracy of solution may reduce as time increases due to the time truncation errors of time derivative term . The cubic trigonometric B-spline method used in this paper is simple and straightforward to apply. An advantage of using the cubic trigonometric B-spline method outlined in this paper is that it produces a spline function on each new time line which can be used to obtain the solutions at any intermediate point in the space direction whereas the finite difference approach yields the solution only at the selected points. The CuTBSM has approximated the solution with more accurate results for time step size , , as compared to some finite difference schemes with smaller time step size such as BTCS, Crandall’s formula, FTCS, the Dufort-Frankel scheme, and TMOL based on reproducing kernel. The proposed method is shown to be unconditionally stable. It is also evident from the examples that the approximate solution is very close to the exact solution.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors are grateful to the anonymous reviewers for their helpful, valuable comments, and suggestions to improve the paper. This study was fully supported by FRGS Grant no. 203/PMATHS/6711324 from the School of Mathematical Sciences, Universiti Sains Malaysia, Penang, Malaysia. The first author was supported by Post Doctorate Fellowship from School of Mathematical Sciences, Universiti Sains Malaysia, Penang, Malaysia, during a part of the time in which the research was carried out.
- Y. S. Choi and K.-Y. Chan, “A parabolic equation with nonlocal boundary conditions arising from electrochemistry,” Nonlinear Analysis: Theory, Methods & Applications, vol. 18, no. 4, pp. 317–331, 1992.
- A. Bouziani, “On a class of parabolic equations with a nonlocal boundary condition,” Académie Royale de Belgique. Bulletin de la Classe des Sciences. 6e Série, vol. 10, no. 1–6, pp. 61–77, 1999.
- W. A. Day, “Parabolic equations and thermodynamics,” Quarterly of Applied Mathematics, vol. 50, no. 3, pp. 523–533, 1992.
- A. A. Samarskiĭ, “Some problems of the theory of differential equations,” Differential Equations, vol. 16, no. 11, pp. 1925–1935, 1980.
- A. R. Bahadır, “Application of cubic B-spline finite element technique to the termistor problem,” Applied Mathematics and Computation, vol. 149, no. 2, pp. 379–387, 2004.
- J. H. Cushman, B. X. Hu, and F. Deng, “Nonlocal reactive transport with physical and chemical heterogeneity: localization errors,” Water Resources Research, vol. 31, no. 9, pp. 2219–2237, 1995.
- F. Kanca, “The inverse problem of the heat equation with periodic boundary and integral overdetermination conditions,” Journal of Inequalities and Applications, vol. 2013, article 108, 2013.
- M. Dehghan, “Efficient techniques for the second-order parabolic equation subject to nonlocal specifications,” Applied Numerical Mathematics, vol. 52, no. 1, pp. 39–62, 2005.
- J. Martín-Vaquero and J. Vigo-Aguiar, “A note on efficient techniques for the second-order parabolic equation subject to non-local conditions,” Applied Numerical Mathematics, vol. 59, no. 6, pp. 1258–1264, 2009.
- X. Li and B. Wu, “New algorithm for nonclassical parabolic problems based on the reproducing kernel method,” Mathematical Sciences, vol. 7, article 4, 2013.
- A. Golbabai and M. Javidi, “A numerical solution for non-classical parabolic problem based on Chebyshev spectral collocation method,” Applied Mathematics and Computation, vol. 190, no. 1, pp. 179–185, 2007.
- W. A. Day, “A decreasing property of solutions of parabolic equations with applications to thermoelasticity,” Quarterly of Applied Mathematics, vol. 40, no. 4, pp. 468–475, 1983.
- W. A. Day, “Extensions of a property of the heat equation to linear thermoelasticity and other theories,” Quarterly of Applied Mathematics, vol. 40, no. 3, pp. 319–330, 1982.
- W. A. Day, Heat Conduction within Linear Thermoelasticity, vol. 30 of Springer Tracts in Natural Philosophy, Springer, New York, NY, USA, 1985.
- J. C. Strikwerda, Finite Difference Schemes and Partial Differential Equations, Society for Industrial and Applied Mathematics, Philadelphia, Pa, USA, 2nd edition, 2004.
- Du. V. Rosenberg, Methods for Solution of Partial Differential Equations, vol. 113, American Elsevier, New York, NY, USA, 1969.
- W. T. Ang, “A method of solution for the one-dimensional heat equation subject to nonlocal conditions,” Southeast Asian Bulletin of Mathematics, vol. 26, no. 2, pp. 197–203, 2002.
- J. Martín-Vaquero and J. Vigo-Aguiar, “On the numerical solution of the heat conduction equations subject to nonlocal conditions,” Applied Numerical Mathematics, vol. 59, no. 10, pp. 2507–2514, 2009.
- M. Dehghan, “On the numerical solution of the diffusion equation with a nonlocal boundary condition,” Mathematical Problems in Engineering, vol. 2, pp. 81–92, 2003.
- M. Dehghan, “A computational study of the one-dimensional parabolic equation subject to nonclassical boundary specifications,” Numerical Methods for Partial Differential Equations, vol. 22, no. 1, pp. 220–257, 2006.
- M. Tatari and M. Dehghan, “On the solution of the non-local parabolic partial differential equations via radial basis functions,” Applied Mathematical Modelling, vol. 33, no. 3, pp. 1729–1738, 2009.
- J. H. Cushman, “Diffusion in fractal porous media,” Water Resources Research, vol. 27, no. 4, pp. 643–644, 1991.
- V. V. Shelukhin, “A non-local in time model for radionuclides propagation in Stokes fluid,” Siberian Branch of Russian Academy of Sciences, Institute of Hydrodynamics, no. 107, pp. 180–193, 1993.
- A. S. Vasudeva Murthy and J. G. Verwer, “Solving parabolic integro-differential equations by an explicit integration method,” Journal of Computational and Applied Mathematics, vol. 39, no. 1, pp. 121–132, 1992.
- C. V. Pao, “Numerical methods for nonlinear integro-parabolic equations of Fredholm type,” Computers & Mathematics with Applications, vol. 41, no. 7-8, pp. 857–877, 2001.
- C. V. Pao, “Reaction diffusion equations with nonlocal boundary and nonlocal initial conditions,” Journal of Mathematical Analysis and Applications, vol. 195, no. 3, pp. 702–718, 1995.
- J. R. Cannon and A. L. Matheson, “A numerical procedure for diffusion subject to the specification of mass,” International Journal of Engineering Science, vol. 31, no. 3, pp. 347–355, 1993.
- M. Dehghan, “Numerical solution of a parabolic equation with non-local boundary specifications,” Applied Mathematics and Computation, vol. 145, no. 1, pp. 185–194, 2003.
- G. Ekolin, “Finite difference methods for a nonlocal boundary value problem for the heat equation,” BIT. Numerical Mathematics, vol. 31, no. 2, pp. 245–261, 1991.
- G. Fairweather and J. C. López-Marcos, “Galerkin methods for a semilinear parabolic problem with nonlocal boundary conditions,” Advances in Computational Mathematics, vol. 6, no. 3-4, pp. 243–262, 1996.
- Y. Liu, “Numerical solution of the heat equation with nonlocal boundary conditions,” Journal of Computational and Applied Mathematics, vol. 110, no. 1, pp. 115–127, 1999.
- G. E. Farin, Curves and Surfaces for Computer-Aided Geometric Design: A Practical Code, Academic Press, 1996.
- J. Hoschek and D. Lasser, Fundamentals of Computer Aided Geometric Design, A K Peters, Boston, Mass, USA, 1993.
- Y. S. Lai, W. P. Du, and R. H. Wang, “The viro method for construction of piecewise algebraic hypersurfaces,” Abstract and Applied Analysis, vol. 2013, Article ID 690341, 7 pages, 2013.
- R. H. Wang, X. Q. Shi, Z. X. Luo, and Z. X. Su, Multivariate Spline and Its Applications, Science Press, Beijing, China; Kluwer Academic Publishers, New York, NY, USA, 2001.
- P. M. Prenter, Splines and Variational Methods, John Wiley & Sons, New York, NY, USA, 1989.
- C. de Boor, A Practical Guide to Splines, vol. 27 of Applied Mathematical Sciences, Springer, 1978.
- H. N. Caglar, S. H. Caglar, and E. H. Twizell, “The numerical solution of third-order boundary-value problems with fourth-degree -spline functions,” International Journal of Computer Mathematics, vol. 71, no. 3, pp. 373–381, 1999.
- H. N. Çaglar, S. H. Çaglar, and E. H. Twizell, “The numerical solution of fifth-order boundary value problems with sixth-degree -spline functions,” Applied Mathematics Letters, vol. 12, no. 5, pp. 25–30, 1999.
- I. Daǧ, D. Irk, and B. Saka, “A numerical solution of the Burgers' equation using cubic B-splines,” Applied Mathematics and Computation, vol. 163, no. 1, pp. 199–211, 2005.
- J. Rashidinia, M. Ghasemi, and R. Jalilian, “A collocation method for the solution of nonlinear one-dimensional parabolic equations,” Mathematical Sciences Quarterly Journal, vol. 4, no. 1, pp. 87–104, 2010.
- K. Qu, Z. Wang, and B. Jiang, “A finite element method by using bivariate splines for one dimensional heat equations,” Journal of Information & Computational Science, vol. 10, no. 12, pp. 3659–3666, 2013.
- J. Goh, A. A. Majid, and A. I. M. Ismail, “Numerical method using cubic B-spline for the heat and wave equation,” Computers and Mathematics with Applications, vol. 62, no. 12, pp. 4492–4498, 2011.
- J. Goh, A. A. Majid, and A. I. M. Ismail, “Cubic B-spline collocation method for one-dimensional heat and advection-diffusion equations,” Journal of Applied Mathematics, vol. 2012, Article ID 458701, 8 pages, 2012.
- J. Goh, A. A. Majid, and A. I. M. Ismail, “A comparison of some splines-based methods for the one-dimensional heat equation,” Proceedings of World Academy of Science, Engineering and Technology, vol. 70, pp. 858–861, 2010.
- R. C. Mittal and G. Arora, “Quintic B-spline collocation method for numerical solution of the Kuramoto-Sivashinsky equation,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 10, pp. 2798–2808, 2010.
- R. C. Mittal and G. Arora, “Numerical solution of the coupled viscous Burgers' equation,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 3, pp. 1304–1313, 2011.
- M. Abbas, A. A. Majid, A. I. M. Ismail, and A. Rashid, “Numerical method using cubic B-spline for a strongly coupled reaction-diffusion system,” PLoS ONE, vol. 9, no. 1, article e83265, 2014.
- M. Abbas, A. A. Majid, A. I. M. Ismail, and A. Rashid, “The application of cubic trigonometric B-spline to the numerical solution of the hyperbolic problems,” Applied Mathematics and Computation, vol. 239, pp. 74–88, 2014.
- A. Nikolis, “Numerical solutions of ordinary differential equations with quadratic trigonometric splines,” Applied Mathematics E-Notes, vol. 4, pp. 142–149, 1995.
- N. N. Abd Hamid, A. A. Majid, and A. I. M. Ismail, “Cubic trigonometric B-spline applied to linear two-point boundary value problems of order two,” World Academy of Science, Engineering and Technology, vol. 70, pp. 798–803, 2010.