Solution scheme development of the nonhomogeneous heat conduction equation in cylindrical coordinates with Neumann boundary conditions by finite difference method

Abstract

Partial differential heat conduction equations are typically used to determine temperature distribution within any solid domain. The difficulty and complexity of the solution of the equation depend on differential equation characteristics, boundary conditions, coordinate systems, and the number of dependent variables. In the current study, the numerical solution schemes were developed by the Explicit Finite Difference and the Implicit Method- the Crank-Nicolson techniques for the partial differential heat conduction equation including heat generation term described as one-dimensional, time-dependent with the Neumann boundary conditions. The solution schemes were, then, applied to the battery problem including highly varying heat generation. Besides, the solution of the problem was performed by using Matlab pdepe solver to verify the developed schemes. Results suggest that the Crank-Nicolson scheme is unconditionally stable, whereas the explicit scheme is only stable when the Courant-Friedrichs-Lewy condition requirement is less than 0.3404. Comparing the developed schemes to the results obtained from the pdepe solver, the schemes are as reliable as the pdepe solver with certain grid structures. Besides, the developed numerical schemes allow for shorter computational times than the pdepe solver at the same grid structures when considering CPU times.

Keywords

• Heat conduction
• heat generation
• Neumann boundary condition
• Crank-Nicolson
• pdepe solver

References

1. Çengel, Y. A. (2007). Heat and mass transfer: A practical approach (3rd ed.). McGraw-Hill Education.
2. Götschel, S., & Weiser, M. (2019). Compression challenges in large scale partial differential equation solvers. Algorithms, 12(9), 197.
3. Zada, L., Nawaz, R., Nisar, K. S., Tahir, M., Yavuz, M., Kaabar, M. K. A., & Martinez, F. (2021). New approximate-analytical solutions to partial differential equations via auxiliary function method. Partial Differential Equations in Applied Mathematics, 4, 100045.
4. Chiang, K.-T., Kung, K.-Y., & Srivastava, H. M. (2009). Analytic transient solutions of a cylindrical heat equation with a heat source. Applied Mathematics and Computation, 215, 2877–2885.
5. Cooper, J. (1998). Boundary value problems for heat equations. In Introduction to partial differential equations with Matlab (pp. 131–156). Birkhauser. https://doi.org/10.1007/978-1-4612-1754-1_4
6. Hanafi, L., Mardlijah, M., Utomo, D. B., & Amiruddin, A. (2021). Study numerical scheme of finite difference for solution partial differential equation of parabolic type to heat conduction problem. Journal of Physics: Conference Series, 1821, 012032. https://doi.org/10.1088/1742-6596/1821/1/012032
7. Taler, J., & Oclon, P. (2014). Finite element method in steady-state and transient heat conduction. In Encyclopedia of thermal stresses (pp. 1594–1600). Springer. https://doi.org/10.1007/978-94-007-2739-7_897
8. Li, W., Yu, B., Wang, X., Wang, P., & Sun, S. (2012). A finite volume method for cylindrical heat conduction problems based on local analytical solution. International Journal of Heat and Mass Transfer, 55, 5570–5582.

9. Peiro, J., & Sherwin, S. (2005). Finite difference, finite element and finite volume methods for partial differential equations. In Handbook of materials modeling (pp. 2417–2431). Springer. https://doi.org/10.1007/978-1-4020-3286-8_127
10. Chapra, S. C., & Canale, R. P. (2015). Numerical methods for engineers (7th ed.). McGraw-Hill Education.
11. Huang, P., Feng, X., & Liu, D. (2013). A stabilized finite element method for the time-dependent Stokes equations based on Crank-Nicolson scheme. Applied Mathematical Modelling, 37, 1910–1919.
12. Mojumder, M., Haque, M., & Alam, M. (2023). Efficient finite difference methods for the numerical analysis of one-dimensional heat equation. Journal of Applied Mathematics and Physics, 11, 3099–3123. https://doi.org/10.4236/jamp.2023.1110204
13. Suarez-Carreno, F., & Rosales-Romero, L. (2021). Convergency and stability of explicit and implicit schemes in the simulation of the heat equation. Applied Sciences, 11, 4468. https://doi.org/10.3390/app11104468
14. Han, F., & Dai, W. (2013). New higher-order compact finite difference schemes for 1D heat conduction equations. Applied Mathematical Modelling, 37, 7940–7952.
15. Whole, A., Lobo, M., & Ginting, K. B. (2021). The application of finite difference method on 2-D heat conductivity problem. Journal of Physics: Conference Series, 2017, 012009.
16. Rieth, A., Kovacs, R., & Fülöp, T. (2018). Implicit numerical schemes for generalized heat conduction equation. International Journal of Heat and Mass Transfer, 126, 1177–1182.
17. Yosaf, A., Rehman, S. U., Ahmad, F., Ullah, M. Z., & Alshomrani, A. S. (2016). Eighth-order compact finite difference scheme for 1D heat conduction equation. Advances in Numerical Analysis, 2016, Article ID 8376061, 12 pages. https://doi.org/10.1155/2016/8376061
18. Dai, W. (2010). A new accurate finite difference scheme for Neumann (insulated) boundary condition of heat conduction. International Journal of Thermal Science, 49, 571–579.
19. Kumar, D., Kumar, V., & Singh, V. P. (2010). Mathematical modeling of brown stock washing problems and their numerical solution using Matlab. Computers & Chemical Engineering, 34, 9–16.
20. Li, Y.-X., Mishra, S. R., Pattnaik, P. K., Baag, S., Li, Y.-M., Khan, M. I., Khan, N. B., Alaoui, M. K., & Khan, S. U. (2022). Numerical treatment of time dependent magnetohydrodynamic nanofluid flow of mass and heat transport subject to chemical reaction and heat source. Alexandria Engineering Journal, 61, 2484–2491.
21. Courtier, N. E., Richardson, G., & Foster, J. M. (2018). A fast and robust numerical scheme for solving models of charge carrier transport and ion vacancy motion in perovskite solar cells. Applied Mathematical Modelling, 68, 329–348.
22. Yudianto, D., & Yuebo, X. (2010). A comparison of some numerical method in solving 1-D steady-state advection dispersion reaction equation. Civil Engineering and Environmental Systems, 27(2), 155–172.
23. MathWorks. (n.d.). Partial differential equations. MathWorks. Retrieved February 8, 2024, from https://www.mathworks.com/help/matlab/math/partial-differential-equations.html
24. Hwang, F. S., Confrey, T., Scully, S., Callaghan, D., Nolan, C., Kent, N., & Flannery, B. (2020). Modeling of heat generation in an 18650 lithium-ion battery cell under varying discharge rates. In 5th Thermal and Fluids Engineering Conference (pp. 333–341). New Orleans, LA, USA.
25. Wang, Z., Ma, J., & Zhang, L. (2017). Finite element thermal model and simulation for a cylindrical Li-ion battery. IEEE Access, 5, 15372–15379.
26. Gümüşsu, E., Ekici, Ö., & Köksal, M. (2017). 3D CFD modeling and experimental testing of thermal behavior of Li-ion battery. Applied Thermal Engineering, 120, 484–495.
27. Rostamy, D., & Abdollahi, N. (2019). Stability analysis for some numerical schemes of partial differential equation with extra measurements. Hacettepe Journal of Mathematics & Statistics, 48(5), 1324–1335.
28. Trivellato, F., & Castelli, M. R. (2014). On the Courant-Friedrichs-Lewy criterion of rotating grids in 2D vertical-axis wind turbine analysis. Renewable Energy, 62, 53–62.
29. Liu, W., & Wu, B. (2022). Unconditional stability and optimal error estimates of a Crank-Nicolson Legendre-Galerkin method for the two-dimensional second-order wave equation. Numerical Algorithms, 90, 137–158. https://doi.org/10.1007/s11075-021-01182-x
30. Majchrzak, E., & Mochnacki, B. (2017). Implicit scheme of the finite difference method for 1D dual-phase lag equation. Journal of Applied Mathematics and Computational Mechanics, 16(3), 37–46.
31. Khalifa, I. (2020). Comparing numerical methods for solving the Fisher equation (Master’s thesis, Lappeenranta-Lahti University of Technology LUT, Finland).