Hermite Interpolation Approach to High-Order Approximation of Heat Equations
Article Main Content
It is usually desirable to approximate the solution of mathemati- cal problems with high-order of accuracy and preferably using com- pact stencils. This work presents an approach for deriving high-order compact discretization of heat equation with source term. The key contribution of this work is the use of Hermite polynomials to reduce second order spatial derivatives to lower order derivatives. This does not involve the use of the given equation, so it is universal. Then, Tay- lor expansion is used to obtain a compact scheme for first derivatives. This leads to a fourth-order approximation in space. Crank-Nicholson scheme is then applied to derive a fully discrete scheme. The result- ing scheme coincides with the fourth-order compact scheme, but our derivation follows a different philosophy which can be adapted for other equations and higher order accuracy. Two numerical experiments are provided to verify the fourth-order accuracy of the approach.
References
-
Nwaigwe C, Orji CU. Second-order non-oscillatory scheme for simulating a pressure-driven flow. Journal of the Nigerian Association of Mathematical Physics, 2019; 52(1):53-58.
Google Scholar
1
-
Lele SK. Compact Finite Difference Schemes with Spectral-like resolution. Journal of Computational Physics, 1992; 103(1):16-42.
Google Scholar
2
-
Gao G, Sun Z. A Compact Finite Difference Scheme for the fractional sub-diffusion equations. Journal of Computational Physic, 2011; 230(3):586-595.
Google Scholar
3
-
Li J, Chen Y. Computational partial differential equations using MATLAB R . Crc Press, 2019.
Google Scholar
4
-
Dlamini PG, Motsa SS, Khumalo M. Higher order compact finite difference schemes for unsteady boundary layer flow problems. Advances in Mathematical Physics, 2013.
Google Scholar
5
-
Bhatt HP, Khaliq AM. Fourth-order compact schemes for the numerical simulation of coupled burgers equation. Computer Physics Communications, 2016; 200:117–138.
Google Scholar
6
-
Chen J, Ge Y. High order locally one-dimensional methods for solving two-dimensional parabolic equations. Advances in Difference Equations, 2018; 2018(1):1–17.
Google Scholar
7
-
Ran M, Luo T, Zhang L. Unconditionally stable compact theta schemes for solving the linear and semi-linear fourth-order diffusion equations. Applied Mathematics and Computation, 2019; 342:118– 129.
Google Scholar
8
-
Ozbenli E, Vedula P. Construction of invariant compact finite-difference schemes. Physical Review E, 2020; 101(2):023303.
Google Scholar
9
-
Shukla RK, Zhong X. Derivation of high-order compact finite difference schemes for non-uniform grid using polynomial inter- polation. Journal of Computational Physics, 2005; 204(2):404–429.
Google Scholar
10
-
Hanh TTH, Thanh DNH, et al. Monotone finite-difference schemes with second order approximation based on regularization approach for the dirichlet boundary problem of the gamma equation. IEEE Access, 2020; 8:45119–45132.
Google Scholar
11
-
Matus P, Gaspar F, Tuyen VTK, et al. Monotone difference schemes for weakly coupled elliptic and parabolic systems. Computational Methods in Applied Mathematics, 2017; 17(2):287–298.
Google Scholar
12
-
Matus P, Vulkov LG, et al. Analysis of secondorder difference schemes on nonuniform grids for quasilinear parabolic equations. Journal of Computational and Applied Mathematics, 2017; 310:186–199.
Google Scholar
13
-
Matus PP, Poliakov DB. Consistent two-sided estimates for the solutions of quasilinear parabolic equations and their approximations. Differential Equations, 2017; 53(7):964–973.
Google Scholar
14
-
Matus P. On convergence of difference schemes for ibvp for quasi- linear parabolic equations with generalized solutions. Computational Methods in Applied Mathematics, 2014; 14(3):361–371.
Google Scholar
15
-
Nwaigwe C. A new second-order solver for nonlinear convection- diffusion-reaction equations. 2020.
Google Scholar
16
-
Nwaigwe C. An unconditionally stable scheme for two- dimensional convection-diffusion-reaction equations. 2022.
Google Scholar
17
-
Guennebaud G, Jacob B, et al. Eigen v3. http://eigen.tuxfamily.org. 2010.
Google Scholar
18
-
Nwaigwe C. Analysis and application of a convergent difference scheme to nonlinear transport in a brinkman flow. International Journal of Numerical Methods for Heat & Fluid Flow, 2020; 30(10):4453-4473.
Google Scholar
19
-
Nwaigwe C, Makinde OD. Finite difference investigation of a polluted non-isothermal non-newtonian porous media flow. Accepted by Defect and Diffusion Forum special issue on Engineering Fluid Flows and Heat Transfer Analysis. 2019; 26(4):145-156.
Google Scholar
20
-
Nwaigwe C, Weli A, Makinde OD.Computational analysis of porous channel flow with cross-diffusion. American Journal of Computational and Applied Mathematics, 2019; 9(5):119–132.
Google Scholar
21
-
Nwaigwe C, Benedict DN. Generalized banach fixed-point theorem and numerical discretization for nonlinear volterra-fredholm equations. Journal of Computational and Applied Mathematics, 2023; 425: 115019.
Google Scholar
22
-
Nwaigwe C, Weli A. Analysis of two finite difference schemes for a channel flow problem. Asian Research Journal of Mathematics, 2019; 15:1-14.
Google Scholar
23