Research Article  Open Access
A CrankNicolson Difference Scheme for Solving a Type of Variable Coefficient Delay Partial Differential Equations
Abstract
A linearized CrankNicolson difference scheme is constructed to solve a type of variable coefficient delay partial differential equations. The difference scheme is proved to be unconditionally stable and convergent, where the convergence order is two in both space and time. A numerical test is provided to illustrate the theoretical results.
1. Introduction
In the past few years, many scholars pay their attention to the theory of delay differential equations (DDEs) [1, 2]. There are many research results on delay ordinary differential equations [3, 4]; however, only few scholars focus on studies of delay partial differential equations. As we know, since, in most cases, DDEs’ exact solutions cannot be computed analytically, efficient numerical methods are needed to solve such equations.
In this paper, the numerical solutions of the following variable coefficient delay partial differential equations are considered: where is the constant diffusion coefficient, is the delay term, , , and . In the case of , numerical solutions of (1)–(3) have been considered in [5–7]. A CrankNicolson scheme and a linearized compact difference scheme have been proposed by Zhang and Sun in [5] and [6], respectively. Q. F. Zhang and C. J. Zhang considered a new linearized compact multisplitting scheme in [7]. We will construct a CrankNicolson scheme for solving (1)–(3). The unconditional stability and convergence will be shown in this paper, where the convergence order is two in both space and time. To testify the theoretical results, a numerical test is provided.
The paper is organized as follows. In Section 2, a linearized CrankNicolson scheme is constructed to solve (1)–(3). Section 3 considers the solvability, stability, and convergence of the CrankNicolson scheme. In Section 4, a numerical test is provided to illustrate the theoretical results. Section 5 gives a brief discussion of this paper.
2. Construction of the Linearized CrankNicolson Scheme
In this subsection, a linearized CrankNicolson scheme for solving (1)–(3) is constructed. In this paper, we make the following assumptions:(H1)assume that (1)–(3) had a unique solution , , and its partial derivatives are bounded by a constant ;(H2) has second derivatives, and we denote where , and are constants.
Two positive integers and are taken; then let , , , and . Define , where and , . Denote , , . Let be the grid function space defined on . Introduce the following notations:
Considering (1) at the point , we have
From Taylor expansion, where . Substituting (8) into (7) and denoting , we obtain where Discretizing the initial and boundary conditions of (2) and (3), we obtain Replacing by and omitting , we obtain the following CrankNicolson scheme:
3. The Solvability, Convergence, and Stability of the CrankNicolson Scheme
Define the following grid function space on :
If , introducing the following notations: The following two inequalities are satisfied [8]: For the analysis of the difference scheme, the following Lemma is needed.
Lemma 1 (see [8]). Let be nonnegative sequence and satisfy then where and are nonnegative constants.
Theorem 2. The difference scheme (13) has a unique solution, under the condition that and are small enough.
Proof. From the positive definiteness of the coefficient matrix of the scheme (13), we can easily obtain the results of Theorem 2 by the mathematical induction method.
Denoting , , , subtracting (13) from (9), (11), and (12), respectively, we obtain the following error equations:
Theorem 3. Letting and be small enough, one has where is independent of and .
Proof. Multiplying (20) by and summing up for from 1 to , we obtain
where
The mathematical induction method will be used to prove Theorem 3. From (21), we have , for . Suppose that (23) is true for , we will prove that (23) is also valid for .
From the inductive assumption, we have
In the following, each term of (24) will be estimated. Consider
From (H1) and (H2), we have
Using the above inequality, we have
Inserting (27)–(29) into (24), we obtain
Taking , we have
The above inequality has the following form:
Summing up (32) for , noticing (21), and exploiting (17), we have
By Lemma 1, we have
where is a constant which depends on , , , , and . From (16), we obtain
By the inductive principle, this completes the proof.
Remark 4. Theorem 3 shows that the convergence order of the variable coefficient delay partial differential equations (1) is . However, for the constant coefficient delay partial differential equations ( in (1)), a CrankNicolson scheme with ) convergence is constructed in [5], and a new difference scheme with ) convergence is constructed in [9].
To discuss the stability of the difference scheme (13), we consider the following problem: The following difference scheme solving for (36) can be obtained: where is a perturbation of .
Similar to the proof of Theorem 3, the following stability result can be obtained.
Theorem 5. Denote Then, there exist constants and such that under the condition that and are small enough and .
Remark 6. Under the condition of assumptions (H1) and (H2) and , for small and , we can get the stability results of Theorem 5 (which can be referred to in [8, 10, 11]), where the difficulty is that ; the proof can be referred to in the proof of Theorem 3.
4. Numerical Test
In this section, a numerical example is considered to validate the algorithm provided in this paper, and the numerical solutions of the example are obtained by exploiting scheme (13). Define
Consider the following problem: where . The exact solution of (41) is .
Table 1 provides some numerical results of difference scheme (13) solving for (41) with step size . Table 2 gives the maximum absolute errors between numerical solutions and exact solutions with different step sizes. From Table 2, we can see that when both the space step size and the time step size are reduced by a factor of 1/2, then the maximum absolute errors are reduced by a factor of approximately 1/4.


Figure 1 provides us with the error curves of numerical solutions for (41) at by using scheme (13). Figures 2 and 3 give the error surface of the numerical solutions with step sizes and , respectively.
Generally speaking, from the results of the tables and the figures provided, we can see that the numerical results are coincident with the theoretical results.
5. Conclusion
In this paper, a type of variable coefficient delay partial differential equations is considered. A linearized CrankNicolson scheme is constructed and is proved to be unconditionally stable and convergent. Finally, a numerical test is provided to illustrate the theoretical results.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work is supported by the Project Sponsored by the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry (no. 2013693), the National Natural Science Foundation of PR China (nos. 71301166, 11301544, 11201487, and 11101184), and the Science Foundation for Young Scientists of Jilin Province (20130522101JH).
References
 D. F. Li and C. J. Zhang, “Nonlinear stability of discontinuous Galerkin methods for delay differential equations,” Applied Mathematics Letters, vol. 23, no. 4, pp. 457–461, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 N. V. Azbelev, V. P. Maksimov, and L. F. Rakhmatullina, Introduction to the Theory of Functional Differential Equations, Nauka, Moscow, Russia, 1991. View at: MathSciNet
 J. K. Hale, Theory of Functional Differential Equations, Springer, New York, NY, USA, 1993.
 O. Diekmann, S. A. van Gils, S. V. Lunel, and H.O. Walther, Delay Equations: Functional, Complex, and Nonlinear Analysis, Springer, New York, NY, USA, 1995. View at: Publisher Site  MathSciNet
 Z.B. Zhang and Z.Z. Sun, “A CrankNicolson scheme for a class of delay nonlinear parabolic differential equations,” Journal on Numerical Methods and Computer Applications, vol. 31, no. 2, pp. 131–140, 2010. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 Z.Z. Sun and Z.B. Zhang, “A linearized compact difference scheme for a class of nonlinear delay partial differential equations,” Applied Mathematical Modelling, vol. 37, no. 3, pp. 742–752, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 Q. F. Zhang and C. J. Zhang, “A new linearized compact multisplitting scheme for the nonlinear convectionreactiondiffusion equations with delay,” Communications in Nonlinear Science and Numerical Simulation, vol. 18, no. 12, pp. 3278–3288, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 Z.Z. Sun, The Numerical Methods for Partial Equations, Science Press, Beijing, China, 2005. (Chinese).
 Y. F. Jin, J. X. Jiang, C. M. Hou, and D. H. Guan, “New difference scheme for general delay parabolic equations,” Journal of Information and Computational Science, vol. 9, article 18, pp. 5579–5586, 2012. View at: Google Scholar
 C. M. Huang and S. Vandewalle, “An analysis of delaydependent stability for ordinary and partial differential equations with fixed and distributed delays,” SIAM Journal on Scientific Computing, vol. 25, no. 5, pp. 1608–1632, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. M. Huang and S. Vandewalle, “Unconditionally stable difference methods for delay partial differential equations,” Numerische Mathematik, vol. 122, no. 3, pp. 579–601, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2014 Wei Gu and Peng Wang. 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.