Recent Advances in Function Spaces and its Applications in Fractional Differential Equations 2020
View this Special IssueResearch Article  Open Access
Yating Huang, Zhe Yin, "The Compact Finite Difference Method of TwoDimensional Cattaneo Model", Journal of Function Spaces, vol. 2020, Article ID 6301757, 12 pages, 2020. https://doi.org/10.1155/2020/6301757
The Compact Finite Difference Method of TwoDimensional Cattaneo Model
Abstract
In this paper, we propose and analyze the compact finite difference scheme of the twodimensional Cattaneo model. The stability and convergence of the scheme are proved by the energy method, the convergence orders are in time and in space. We also use the variables separation method to find the true solution of the problem. On this basis, the validity and accuracy of the scheme are verified by numerical experiments.
1. Introduction
In the diffusion phenomenon described by traditional Fick law and Fourier law, the moment a point passes by, it is disturbed at a great distance, the propagation speed of the disturbance seems to be infinite. However, this property is unphysical. To solve the problem, Cattaneo modified the constitutive equation by introducing a relaxation parameter which plays the role of relaxation time, it should be noted that the relaxation parameter is small and depends on the thermodynamic properties of the material. Thus, he proposed the Cattaneo model [1].
From a mathematical point of view, Cattaneo model transforms the traditional diffusion equation into a hyperbolic equation, the speed of propagation is finite, and improves the property of infinite propagation speed. On the other hand, because of the hyperbolic nature of the Cattaneo model, it may have oscillatory solutions and negative values.
From a physical point of view, the Cattaneo model describes a physical phenomenon which is called heat waves. Although this phenomenon can only be observed under special circumstances, materials, or states [2], it is still gradually accepted by the public. It can be used to describe not only ultrasonic wave propagation in certain diluted gases but also heat pulse propagation in some pure nonmetallic crystals [3]. Straughan et al. [4] studied thermal convection with the CattaneoChristov model in horizontal layers of an incompressible Newtonian fluid. Haddad [5] utilized the theory of CattaneoChristov to investigate Brinkman's porous media. Cattaneo model is widely used in extended irreversible thermodynamics, cosmological models, and crystalline solids [6–8].
Because the fractional models have the advantage of previous memory, scholars have studied various fractional models [9–13]. With the wide application of the Cattaneo model in physics and theoretical analysis, more and more people pay attention to the generalization and solution of the Cattaneo model. Compte and Metzler [14] generalized the Cattaneo model from three perspectives which called continuoustime random walks, nonlocal transport theory, and delayed fluxforce relation. Ferrillo et al. [15] compared the Cattaneo model and fractional Cattaneo model and studied the asymptotic behavior of solutions to the Cattaneo equations. Su et al. [16] studied solutions to the classical Cattaneo equation and compared with the solutions of a phaselagging heat transport equation. C ̆iegis [17] studied the numerical solutions of a class of hyperbolic heat conduction equations, proposed explicit and implicit Euler schemes for such equations, and verified that the convergence order of this scheme was . Huang and Yin [18] proposed a fourthorder compact difference scheme for the D Cattaneo model and verified that the convergence order of the scheme was through numerical experiments. Zhao and Sun [19] presented compact CrankNicolson schemes for a class fractional Cattaneo equation in an inhomogeneous medium and used the energy method to verify that the convergence order of the scheme is . Li and Cao [20] proposed an unconditionally stable scheme with convergence order of . Vong et al. [21] proposed a higherorder difference scheme for a class of generalized Cattaneo equations and proved the stability and convergence of the scheme by energy method.
There are many methods to solve partial differential equations, such as the finite element method [22, 23], finite volume element method [24, 25], and finite difference method [26]. Compared with other methods, the finite difference method has the advantages of a low requirement for grid nodes and simple calculation. In 1992, Lele [27] first proposed a general form of compact difference scheme. Compact difference scheme is a kind of finite difference scheme with high precision, which is widely used in solving parabolic equation and hyperbolic equation [28–33].
Previous scholars used different methods to study the integerorder Cattaneo model and the fractional Cattaneo model. The onedimensional Cattaneo model describes the heat conduction phenomenon of the uniform thin column with side insulation, while the twodimensional Cattaneo model needs to be studied for the heat conduction phenomenon of the thin plate with upper and lower bottom insulation. In this paper, we construct a compact difference scheme for a twodimensional Cattaneo model. The time derivative is discretized with a central difference, and the space derivative is discretized with compact operators, so that the convergence order of the scheme reaches . We use the energy method to verify the stability and convergence of the scheme. Finally, we verify the validity and accuracy of the scheme through numerical experiments. It should be noted that the above articles are all about constructing the exact solution of the problem, and then calculating the right term, so as to conduct numerical experiments. In this paper, our example is the definite problem corresponding to the homogeneous equation. The exact solution is obtained through the separation of variables method, and then the numerical experiment is conducted.
The arrangement of the article is as follows. In the second part, we give some lemmas and symbols to derive a compact difference scheme. In the third part, we give a strict proof of the stability and convergence of the scheme by energy method. In the fourth part, we obtain the true solution of the specific problem through the separation of variable method and use numerical experiments to prove the validity of the scheme. In the fifth part, we make a final summary of the article.
2. Construction of the Compact Finite Difference Scheme
In this chapter, we consider the following twodimensional Cattaneo model: where is a rectangular region in , is the boundary of the rectangular region, is a relaxation parameter, is the diffusion coefficient, and and are known functions. The above model is considered under homogeneous Dirichlet boundary conditions. The regularity assumption that problem (1) has a solution will be used in the following analysis.
In order to propose the compact finite difference scheme, we first introduce some basic symbols and preliminary lemmas. Take three positive integers , , and , then let , , . In this case, the spatial node can be expressed as , where ; ; , let , , , , , , . Let be a grid function space defined on . For any grid function , introduce the following notations:
The compact operators are defined as follows
Let be grid function space defined on . For any grid function , introduce the following notations:
The following lemmas also play an important role in constructing a scheme.
Lemma 1 [34]. Suppose , it holds that where is a constant.
Lemma 2 [34]. Suppose , it holds that
Lemma 3 [34]. Suppose , it holds that where is a constant.
Lemma 4 [34]. Suppose , it holds that
where , is a constant.
Define the grid functions as , . We discrete at and refer to the formula [34].
We can obtain
where
Utilizing Taylor formula results in
Using the above equation and Lemma 3 result in
where .
We apply compact operators and to both sides of Eq. (13) and apply Lemma 4, we get
where
We think about the differential equation (1) at point ,
By applying Lemma 1 and Lemma 2, we have
Substituting the above equation into (16) and applying to both sides, we have
where
Omitting small terms and . The numerical approximate solution is represented by the grid function . Replacing with it, the fourthorder compact difference scheme can be obtained
From the above analysis, the following theorems can be obtained.
Theorem 5. The truncation error of compact difference scheme (20) is
3. Analysis of the Compact Finite Difference Scheme
Let , we define the following inner products and the corresponding norms
Next, we introduce the Lemma that will be applied to the stability analysis of the scheme.
Lemma 6 [35]. For any grid function , it holds that Applying the Lemma, we have the following result.
Theorem 7. Suppose is the solution of the following difference scheme It holds that where , is a source item.
Proof. Multiplying Eq. (24) by , and summing up for from to , we have
Next, we deal with the two terms on the lefthand side of Eq. (29), we have
Then we deal with the first term on the righthand side of Eq. (29).
Analogously
Combining (32) and (33), we have
Substituting (30), (31), (34) into (29) and applying Lemma 6, CauchySchwarz inequality, we have
Sorting Eq. (35) and adding to both sides, we have
Let , .
Eq. (36) is written as
Summing over results in
Multiplying Eq. (25) by , summing up for from to , applying Lemma 4 and CauchySchwarz inequality, we have
where .
Sorting Eq. (39), we have
Sorting Eq. (40), we have
Substituting (41) into (38), we get
End of the proof.
Remark 8. According to the results, it can be seen that the derivative term of initial value affects stability and convergence, so this term in the Cattaneo model will be set as zero in general. For details, please refer to the literature [15].
Theorem 9. Suppose is the solution of Eq. (1), is the solution of the difference scheme (20). Let it holds that
Proof. Subtracting (14), (18), and (20), we can get the error system Using the Theorem 7 and noticing (47) and (48), we have End of the proof.
Remark 10. Through Theorem 7 and Theorem 9, we prove the stability and convergence of the scheme. Due to the basic relationship between seminorm and infinity norm, we can also obtain that the proposed scheme is stable and convergent in norm.
4. Numerical Experiments
In this chapter, we verify the numerical accuracy and validity of the proposed scheme.
In the following numerical experiment, for simplicity, the domain is set to the square domain. The domain is , and is equal to . , , , . In order to verify the numerical accuracy of the scheme, we need the following exact solution of the above problem. We can obtain the exact solution by using the separation of variables method,
Before we give the experimental results, we first give the definitions of norm of the error and norm of the error.
For all grid points, norm of the error is defined by
For all grid points, norm of the error is defined by where is the true solution, is the numerical solution, and is the last moment. Let , order of spatial convergence is denoted by and order of temporal convergence is denoted by , which are defined as follows
In Table 1, let , and , for and , take the corresponding values, respectively, the compact difference scheme has the following corresponding error, order of temporal convergence and order of spatial convergence. The third column is the order of spatial convergence, and the fourth column is the order of time convergence. We can find that when the spatial step is reduced to onehalf of the original and the temporal step is reduced to onequarter of the original, the error is reduced to about onesixteenth of the original. At the same time, the compact difference scheme keeps the fourthorder spatial rate of convergence and secondorder temporal rate of convergency. So we also verify the stability and convergence of the scheme.

In order to verify the applicability of the compact difference scheme more accurately, we give the image of convergence order defined the by norm. In Figure 1, we can see more clearly that the convergence order of the scheme (20) is , which also verifies that the scheme (20) proposed in this paper is stable and effective. And the relaxation parameter does not affect the convergence rate.
Figure 2 show the true solution and the numerical solution when and .
Figure 3 shows the images when is taken at different moments. It can be seen from the images that as time increases, the temperature also decreases gradually.
Next, the value of in the previous example is , and the temperature image of at different moments is also obtained. See Figure 4, it has the same property as Figure 3.
Finally, we use two examples to test the influence of relaxation parameters, as shown in Figure 5. From Figure 3–5, we can see that the larger is, the more obvious the fluctuation of temperature distribution is, and the slower the temperature change at the initial time.
In Figure 6, in order to observe the differences between the Cattaneo model and the traditional diffusion equation more conveniently, we give the temperature distribution images of the onedimensional Cattaneo model and the heat conduction equation at different moments at the midpoint of . It can be seen from the Figure 6 that the Cattaneo model has a slower temperature transformation at the initial moment than the heat conduction equation, but at the same time, the former has negative temperature and oscillation, which is also affected by the hyperbolic property of Cattaneo model.
5. Conclusion
In this paper, we propose a compact finite difference scheme for a twodimensional Cattaneo model and prove the stability and convergence of the scheme by energy method. A concrete example is given, and the exact solution is obtained by the method of separating variables, it can be seen from the figures that the fitting degree of the exact solution and numerical solution is higher, it is also verified by numerical experiments that the convergence order of the scheme is . In conclusion, the proposed scheme is simple, stable, and effective and has higher computational accuracy than some existing methods.
Data Availability
Some or all data, models, or code generated or used during the study are available from the corresponding author by request.
Conflicts of Interest
The authors declare that they have no competing interests.
Authors’ Contributions
Yating Huang carried out the main part of this article. All authors read and approved the final manuscript.
Acknowledgments
This work is supported by the National Science Foundation of China (Nos. 11501335), and the Natural Science Foundation of Shandong Province of China (Nos. ZR2017MA020, Nos. ZR2017MA003).
References
 C. Cattaneo, “Sulla conduzione del calore,” Atti del Seminario Matematico e Fisico dell'Università di Modena, vol. 3, pp. 83–101, 1948. View at: Google Scholar
 B. Straughan, Heat Waves, Springer Science Business, New York, 2011.
 K. V. Zhukovsky and H. M. Srivastava, “Analytical solutions for heat diffusion beyond Fourier low,” Applied Mathematics and Computation, vol. 293, pp. 423–437, 2017. View at: Publisher Site  Google Scholar
 B. Straughan, “Thermal convection with the Cattaneo–Christov model,” International Journal of Heat and Mass Transfer, vol. 53, no. 13, pp. 95–98, 2010. View at: Publisher Site  Google Scholar
 S. A. Haddad, “Thermal instability in Brinkman porous media with Cattaneo–Christov heat flux,” International Journal of Heat and Mass Transfer, vol. 68, no. 4, pp. 659–668, 2014. View at: Publisher Site  Google Scholar
 D. Jou, J. CasasVÁzquez, and G. Lebon, “Extended irreversible thermodynamics,” Reports on Progress in Physics, vol. 51, no. 8, pp. 1105–1179, 1988. View at: Publisher Site  Google Scholar
 M. Zakari and D. Jou, “Equations of state and transport equations in viscous cosmological models,” Physical Review D, vol. 48, no. 4, pp. 1597–1601, 1993. View at: Publisher Site  Google Scholar
 S. Godoy and L. S. GarcíaColín, “From the quantum random walk to classical mesoscopic diffusion in crystalline solids,” Physical Review D, vol. 53, pp. 5779–5785, 1996. View at: Publisher Site  Google Scholar
 C. Chen, H. Liu, X. Zheng, and H. Wang, “A twogrid MMOC finite element method for nonlinear variableorder timefractional mobile/immobile advection–diffusion equations,” Computers & Mathematcs with Applications, vol. 79, no. 9, pp. 2771–2783, 2020. View at: Publisher Site  Google Scholar
 Z. J. Zhou and W. Gong, “Finite element approximation of optimal control problems governed by time fractional diffusion equation,” Computers & Mathematcs with Applications, vol. 71, no. 1, pp. 301–318, 2016. View at: Publisher Site  Google Scholar
 Y. S. Li, H. Z. Chen, and H. Wang, “Amixedtype galerkin variational formulation and fast algorithms for variablecoefficient fractional diffusion equations,” Mathematical Methods in the Applied Sciences, vol. 40, no. 14, pp. 5018–5034, 2017. View at: Publisher Site  Google Scholar
 J. X. He, X. G. Zhang, L. S. Liu, Y. H. Wu, and Y. J. Cui, “A singular fractional Kelvin–Voigt model involving a nonlinear operator and their convergence properties,” Boundary Value Problems, vol. 2019, no. 1, Article ID 112, 2019. View at: Publisher Site  Google Scholar
 X. G. Zhang, L. S. Liu, and Y. H. Wu, “Multiple positive solutions of a singular fractional differential equation with negatively perturbed term,” Mathematical and Computer Modelling, vol. 55, no. 34, pp. 1263–1274, 2012. View at: Publisher Site  Google Scholar
 A. Compte and R. Metzler, “The generalized Cattaneo equation for the description of anomalous transport processes,” Journal of Physics A: Mathematical and General, vol. 30, no. 21, pp. 7277–7289, 1997. View at: Publisher Site  Google Scholar
 F. Ferrillo, R. Spigler, and M. Concezzi, “Comparing Cattaneo and fractional derivative models for heat transfer processes,” SIAM Journal on Applied Mathematics, vol. 78, no. 3, pp. 1450–1469, 2018. View at: Publisher Site  Google Scholar
 S. Su, W. Dai, P. M. Jordan, and R. E. Mickens, “Comparison of the solutions of a phaselagging heat transport equation and damped wave equation,” International Journal of Heat and Mass Transfer, vol. 48, no. 11, pp. 2233–2241, 2005. View at: Publisher Site  Google Scholar
 R. Čiegis, “Numerical solution of hyperbolic heat conduction equation,” Mathematical Modelling and Analysis, vol. 14, no. 1, pp. 11–24, 2009. View at: Publisher Site  Google Scholar
 Y. T. Huang and Z. Yin, “The compact finite difference method of Cattaneo model,” Journal of Shandong Normal University, vol. 34, no. 1, pp. 39–42, 2019. View at: Google Scholar
 X. Zhao and Z. Z. Sun, “Compact CrankNicolson schemes for a class of fractional Cattaneo equation in inhomogeneous medium,” Journal of Scientific Computing, vol. 62, no. 3, pp. 747–771, 2015. View at: Publisher Site  Google Scholar
 C. P. Li and J. X. Cao, “A finite difference method for timefractional telegraph equation,” in Proceedings of 2012 IEEE/ASME 8th IEEE/ASME International Conference on Mechatronic and Embedded Systems and Applications, pp. 314–318, Suzhou, China, 2012. View at: Publisher Site  Google Scholar
 S. W. Vong, H. K. Pang, and X. Q. Jin, “A highorder difference scheme for the generalized Cattaneo equation,” East Asian Journal on Applied Mathematics, vol. 2, no. 2, pp. 170–184, 2012. View at: Publisher Site  Google Scholar
 C. J. Chen, K. Li, Y. P. Chen, and Y. Q. Huang, “Twogrid finite element methods combined with CrankNicolson scheme for nonlinear Sobolev equations,” Advances in Computational Mathematics, vol. 45, no. 2, pp. 611–630, 2019. View at: Publisher Site  Google Scholar
 C. J. Chen, X. Y. Zhang, G. D. Zhang, and Y. Y. Zhang, “A twogrid finite element method for nonlinear parabolic integrodifferential equations,” International Journal of Computer Mathematics, vol. 96, no. 10, pp. 2010–2023, 2019. View at: Publisher Site  Google Scholar
 C. J. Chen and X. Zhao, “A posteriori error estimate for finite volume element method of the parabolic equations,” Numerical Methods for Partial Differential Equations, vol. 33, no. 1, pp. 259–275, 2017. View at: Publisher Site  Google Scholar
 C. J. Chen, W. Liu, and C. J. Bi, “A twogrid characteristic finite volume element method for semilinear advectiondominated diffusion equations,” Numerical Methods for Partial Differential Equations, vol. 29, no. 5, pp. 1543–1562, 2013. View at: Publisher Site  Google Scholar
 T. T. Wu and R. M. Xu, “An optimal compact sixthorder finite difference scheme for the Helmholtz equation,” Computers & Mathematics with Applications, vol. 75, no. 7, pp. 2520–2537, 2018. View at: Publisher Site  Google Scholar
 S. K. Lele, “Compact finite difference schemes with spectrallike resolution,” Journal of Computational Physics, vol. 103, no. 1, pp. 16–42, 1992. View at: Publisher Site  Google Scholar
 S. C. Dennid and J. D. Hundson, “Compact h^{4} finite difference approximations to operators of NavierStokes type,” Journal of Computational Physics, vol. 85, no. 2, pp. 390–416, 1989. View at: Publisher Site  Google Scholar
 Y. M. Wang and B. Y. Guo, “Fourthorder compact finite difference method for fourthorder nonlinear elliptic boundary value problems,” Journal of Computational and Applied Mathematics, vol. 221, no. 1, pp. 76–97, 2008. View at: Publisher Site  Google Scholar
 Z. D. Luo, S. J. Jin, and J. Chen, “A reducedorder extrapolation central difference scheme based on POD for twodimensional fourthorder hyperbolic equations,” Applied Mathematics and Computation, vol. 289, pp. 396–408, 2016. View at: Publisher Site  Google Scholar
 J. H. Jia and H. Wang, “A fast finite difference method for distributedorder spacefractional partial differential equations on convex domains,” Computers & Mathematcs with Applications, vol. 75, no. 6, pp. 2031–2043, 2018. View at: Publisher Site  Google Scholar
 X. Shen and A. L. Zhu, “A Crank–Nicolson linear difference scheme for a BBM equation with a time fractional nonlocal viscous term,” Advances in Difference Equations, vol. 2018, Article ID 351, 2018. View at: Publisher Site  Google Scholar
 L. Y. Li, Z. W. Jiang, and Z. Yin, “Fourthorder compact finite difference method for solving twodimensional convectiondiffusion equation,” Advances in Difference Equations, vol. 2018, no. 1, Article ID 234, 2018. View at: Publisher Site  Google Scholar
 Z. Z. Sun, Numerical Methods for Partial Differential Equations, Science Press, Beijing, 2012.
 Z. Z. Sun and G. H. Gao, Finite Difference Methods for Fractional Differential Equations, Science Press, Beijing, 2015.
Copyright
Copyright © 2020 Yating Huang and Zhe Yin. 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.