Rivers State University, Nigeria
* Corresponding author
Rivers State University, Nigeria

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

  1. 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
  2. Lele SK. Compact Finite Difference Schemes with Spectral-like resolution. Journal of Computational Physics, 1992; 103(1):16-42.
     Google Scholar
  3. 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
  4. Li J, Chen Y. Computational partial differential equations using MATLAB R . Crc Press, 2019.
     Google Scholar
  5. 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
  6. 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
  7. 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
  8. 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
  9. Ozbenli E, Vedula P. Construction of invariant compact finite-difference schemes. Physical Review E, 2020; 101(2):023303.
     Google Scholar
  10. 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
  11. 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
  12. 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
  13. 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
  14. 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
  15. 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
  16. Nwaigwe C. A new second-order solver for nonlinear convection- diffusion-reaction equations. 2020.
     Google Scholar
  17. Nwaigwe C. An unconditionally stable scheme for two- dimensional convection-diffusion-reaction equations. 2022.
     Google Scholar
  18. Guennebaud G, Jacob B, et al. Eigen v3. http://eigen.tuxfamily.org. 2010.
     Google Scholar
  19. 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
  20. 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
  21. 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
  22. 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
  23. 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