Journal of Function Spaces
Volume 2020 (2020), Article ID 6301757, 12 pages
https://doi.org/10.1155/2020/6301757

Research Article

The Compact Finite Difference Method of Two-Dimensional Cattaneo Model

Yating Huang and Zhe Yin

School of Mathematics and Statistics, Shandong Normal University, Jinan, Shandong, China

Correspondence should be addressed to Zhe Yin; zyin_sdnu@163.com

Received 1 March 2020; Accepted 15 April 2020; Published 28 May 2020

Guest Editor: Chuanjun Chen

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.

Abstract.

In this paper, we propose and analyze the compact finite difference scheme of the two-dimensional 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 Cattaneo-Christov model in horizontal layers of an incompressible Newtonian fluid. Haddad [5] utilized the theory of Cattaneo-Christov 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 continuous-time random walks, nonlocal transport theory, and delayed flux-force 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 phase-lagging 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 fourth-order 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 Crank-Nicolson 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 higher-order 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 integer-order Cattaneo model and the fractional Cattaneo model. The one-dimensional Cattaneo model describes the heat conduction phenomenon of the uniform thin column with side insulation, while the two-dimensional 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 two-dimensional 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 two-dimensional 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 fourth-order 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 left-hand side of Eq. (29), we have

Then we deal with the first term on the right-hand side of Eq. (29).

Analogously

Combining (32) and (33), we have

Substituting (30), (31), (34) into (29) and applying Lemma 6, Cauchy-Schwarz 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 Cauchy-Schwarz 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 semi-norm 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 one-half of the original and the temporal step is reduced to one-quarter of the original, the error is reduced to about one-sixteenth of the original. At the same time, the compact difference scheme keeps the fourth-order spatial rate of convergence and second-order temporal rate of convergency. So we also verify the stability and convergence of the scheme.

Table 1: Error and convergence order of compact difference scheme.

—
—
7.1854e-03
3.5927e-03

3.517
1.759
5.8076e-04
2.9038e-04

4.005
2.003
3.6204e-05
1.8102e-05

3.999
2.000
2.2637e-06
1.3185e-06

4.000
2.000
1.4151e-07
7.0753e-08

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 1: Order of spatial and temporal convergence.

Figure 2 show the true solution and the numerical solution when and .

Figure 2: The exact solution and the numerical at fixed , , , 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.

Figure 3: Temperature images at different times when , , , and .

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.

Figure 4: Temperature images at different times when , , , and .

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.

Figure 5: Temperature images with different relaxation parameters when , and .

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 one-dimensional 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.

Figure 6: Temperature distribution images at .

5. Conclusion

In this paper, we propose a compact finite difference scheme for a two-dimensional 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.

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.

B. Straughan, “Thermal convection with the Cattaneo–Christov model,” International Journal of Heat and Mass Transfer, vol. 53, no. 1-3, pp. 95–98, 2010.

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.

D. Jou, J. Casas-VÁzquez, and G. Lebon, “Extended irreversible thermodynamics,” Reports on Progress in Physics, vol. 51, no. 8, pp. 1105–1179, 1988.

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.

S. Godoy and L. S. García-Colín, “From the quantum random walk to classical mesoscopic diffusion in crystalline solids,” Physical Review D, vol. 53, pp. 5779–5785, 1996.

C. Chen, H. Liu, X. Zheng, and H. Wang, “A two-grid MMOC finite element method for nonlinear variable-order time-fractional mobile/immobile advection–diffusion equations,” Computers & Mathematcs with Applications, vol. 79, no. 9, pp. 2771–2783, 2020.

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.

Y. S. Li, H. Z. Chen, and H. Wang, “Amixed-type galerkin variational formulation and fast algorithms for variable-coefficient fractional diffusion equations,” Mathematical Methods in the Applied Sciences, vol. 40, no. 14, pp. 5018–5034, 2017.

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.

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. 3-4, pp. 1263–1274, 2012.

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.

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.

S. Su, W. Dai, P. M. Jordan, and R. E. Mickens, “Comparison of the solutions of a phase-lagging heat transport equation and damped wave equation,” International Journal of Heat and Mass Transfer, vol. 48, no. 11, pp. 2233–2241, 2005.

R. Čiegis, “Numerical solution of hyperbolic heat conduction equation,” Mathematical Modelling and Analysis, vol. 14, no. 1, pp. 11–24, 2009.

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.

X. Zhao and Z. Z. Sun, “Compact Crank-Nicolson schemes for a class of fractional Cattaneo equation in inhomogeneous medium,” Journal of Scientific Computing, vol. 62, no. 3, pp. 747–771, 2015.

C. P. Li and J. X. Cao, “A finite difference method for time-fractional 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.

S. W. Vong, H. K. Pang, and X. Q. Jin, “A high-order difference scheme for the generalized Cattaneo equation,” East Asian Journal on Applied Mathematics, vol. 2, no. 2, pp. 170–184, 2012.

C. J. Chen, K. Li, Y. P. Chen, and Y. Q. Huang, “Two-grid finite element methods combined with Crank-Nicolson scheme for nonlinear Sobolev equations,” Advances in Computational Mathematics, vol. 45, no. 2, pp. 611–630, 2019.

C. J. Chen, X. Y. Zhang, G. D. Zhang, and Y. Y. Zhang, “A two-grid finite element method for nonlinear parabolic integro-differential equations,” International Journal of Computer Mathematics, vol. 96, no. 10, pp. 2010–2023, 2019.

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.

C. J. Chen, W. Liu, and C. J. Bi, “A two-grid characteristic finite volume element method for semilinear advection-dominated diffusion equations,” Numerical Methods for Partial Differential Equations, vol. 29, no. 5, pp. 1543–1562, 2013.

T. T. Wu and R. M. Xu, “An optimal compact sixth-order finite difference scheme for the Helmholtz equation,” Computers & Mathematics with Applications, vol. 75, no. 7, pp. 2520–2537, 2018.

S. K. Lele, “Compact finite difference schemes with spectral-like resolution,” Journal of Computational Physics, vol. 103, no. 1, pp. 16–42, 1992.

S. C. Dennid and J. D. Hundson, “Compact h⁴ finite difference approximations to operators of Navier-Stokes type,” Journal of Computational Physics, vol. 85, no. 2, pp. 390–416, 1989.

Y. M. Wang and B. Y. Guo, “Fourth-order compact finite difference method for fourth-order nonlinear elliptic boundary value problems,” Journal of Computational and Applied Mathematics, vol. 221, no. 1, pp. 76–97, 2008.

Z. D. Luo, S. J. Jin, and J. Chen, “A reduced-order extrapolation central difference scheme based on POD for two-dimensional fourth-order hyperbolic equations,” Applied Mathematics and Computation, vol. 289, pp. 396–408, 2016.

J. H. Jia and H. Wang, “A fast finite difference method for distributed-order space-fractional partial differential equations on convex domains,” Computers & Mathematcs with Applications, vol. 75, no. 6, pp. 2031–2043, 2018.

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.

L. Y. Li, Z. W. Jiang, and Z. Yin, “Fourth-order compact finite difference method for solving two-dimensional convection-diffusion equation,” Advances in Difference Equations, vol. 2018, no. 1, Article ID 234, 2018.

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.