Abstract
In this paper, a high-order compact local one-dimensional (LOD) method is proposed to solve the three-dimensional reaction-diffusion equation with nonlinear reaction term. First, the governing equation is split into three one-dimensional equations. Then, the fourth-order compact difference formulas and the Crank-Nicolson method are employed to approximate the space and time derivatives, respectively. The nonlinear reaction term is linearized by utilizing the Taylor expansions. It only needs to solve three one-dimensional tridiagonal linear systems at each time layer, which makes the computation cost-effective. After conducting theoretical analysis, we find that the truncation error of the current LOD method is of order \(O({\tau ^2} + {h^4})\) for source spatially symmetry problems, but only of order \(O(\tau + \tau {h^2})\) for source spatially asymmetry problems. The convergence and stability of the proposed scheme are proved by the energy analysis method. Finally, some numerical experiments are carried out to demonstrate the accuracy and the efficiency of the presented method. Furthermore, our method is also exploited in the simulation of the quenching and blow-up phenomena.
Similar content being viewed by others
Data Availability
Data Availability The data that support this study are available from the corresponding author on reasonable request.
References
Ames WF (1965) Nonlinear partial differential equations in engineering. Academic Press, New York
Ahmed N, Korkmaz A, Rafiq M et al (2020) A novel time efficient structure-preserving splitting method for the solution of two-dimensional reaction–diffusion systems. Adv Differ Equ-ny 2020:1–26
Abia LM, Lpez-Marcos JC, Martnez J (1996) Blow-up for semidiscretizations of reaction–diffusion equations. Appl Numer Math 20:145–156
Beauregard MA, Sheng Q (2012) A semi-adaptive compact splitting method for the numerical solution of 2-dimensional quenching problems. Appl Math Comput 218:11240–11254
Bao A, Song X (2016) Bounds for the blowup time of the solution to a parabolic system with nonlocal factors in nonlinearities. Comput Math Appl 71:723–729
Cho CH (2013) On the computation of the numerical blow-up time. Jpn J Ind Appl Math 30:331–349
Deng K, Levine HA (1989) On the blow up of \(u_t\) at quenching. Proc Am Math Soc 106:1049–1056
Deng D, Zhang C (2015) Analysis and application of a compact multistep ADI solver for a class of nonlinear viscous wave equations. Appl Math Model 39:1033–1049
Frank-Kamenetskii DA (1969) Diffusion and heat transfer in chemical kinetics, 2nd edn. Plenum Press, New York
Ge Y, Cai Z, Sheng Q (2018) A compact adaptive approach for degenerate singular reaction–diffusion equations. Numer Methods Part D E 34:1166–1187
Gu Y, Liao W, Zhu J (2003) An effcient high-order algorithm for solving systems of 3-D reaction–diffusion equations. J Comput Appl Math 155:1–17
Hao A, Zhou J (2016) Blowup, extinction and non-extinction for a nonlocal p-biharmonic parabolic equation. Appl Math Lett 64:198–204
Kawarada H (1975) On solutions of initial-boundary problem for \(u_{t}=u_{xx} + 1/(1 - u)\). Publ Res I Math Sci 10:729–736
Karaa S (2005) An accurate LOD scheme for two-dimensional parabolic problems. Appl Math Comput 170:886–894
Levine HA (1990) Quenching, nonquenching and beyond quenching for solutions of some parabolic equations. Ann Mat Pura Appl 155:243–260
Li J, Chen Y, Liu G (2006) High-order compact ADI methods for parabolic equations. Comput Math Appl 52:1343–1356
Liao W, Zhu J, Khaliq AQM (2002) An efficient high-order algorithm for solving systems of reaction–diffusion equations. Numer Methods Part D E 18:340–354
Liao W (2012) A high-order ADI finite difference scheme for a 3D reaction–diffusion equation with neumann boundary condition. Numer Methods Part D E 29:778–798
Ma J, Jiang Y (2009) Moving mesh methods for blowup in reaction–diffusion equations with traveling heat source. J Comput Phys 228:6977–6990
Nakagawa T (1975) Blowing up of a finite difference solution to \({u_t} = u_{xx} + {u^2}\). Appl Math Optim 2:337–350
Qin J, Wang T (2011) A compact locally one-dimensional finite difference method for nonhomogeneous parabolic differential equations. Int J Numer Methods Biomed Eng 27:128–142
Samarskii A (1964) Local one-dimensional difference schemes for multi-dimensional hyperbolic equations in an arbitrary region. USSR Comput Math Math Phys 4:21–35
Sun ZZ (2005) Numerical methods of partial differential equations. Science Press, Beijing (in Chinese)
Wang YM (2011) Error and extrapolation of a compact LOD method for parabolic differential equations. J Comput Appl Math 235:1367–1382
Wu F, Cheng X, Li D, Duan J (2018) A two-level linearized compact ADI scheme for two-dimensional nonlinear reaction–diffusion equations. Comput Math Appl 75:2835–2850
Wang T, Wang YM (2014) A higher-order compact LOD method and its extrapolations for nonhomogeneous parabolic differential equations. Appl Math Comput 237:512–530
Wang K, Wang H (2018) Stability and error estimates of a new high-order compact ADI method for the unsteady 3D convection-diffusion equation. Appl Math Comput 331:140–159
Wang M, Zhou Y (1996) Exact solutions to the Chaffee-Infante reaction diffusion equation. J Lanzhou Univ (Natural Sciences) 32:26–30 (in Chinese)
Zhu X, Ge Y (2020) Adaptive ADI numerical analysis of 2D quenching-type reaction diffusion equation with convection term. Math Probl Eng 2020:1–19
Zhou H, Wu Y, Tian W (2012) Extrapolation algorithm of compact ADI approximation for two-dimensional parabolic equation. Appl Math Comput 219:2875–2884
Acknowledgements
We would like to thank the editors and anonymous reviewers for their constructive comments and suggestions that are helpful to improve the quality of the paper. This work is partially supported by the National Natural Science Foundation of China (12161067, 11772165, 11961054, 11902170), National Natural Science Foundation of Ningxia (2022AAC02023, 2020AAC03059), the Key Research and Development Program of Ningxia (2018BEE03007), National Youth Top-notch Talent Support Program of Ningxia, and the First Class Discipline Construction Project in Ningxia Universities: Mathematics.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of Interest
The authors declare no potential conflict of interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Wu, Y., Ge, Y. & Zhang, L. A high-order compact LOD method for solving the three-dimensional reaction–diffusion equation with nonlinear reaction term. Comp. Appl. Math. 42, 46 (2023). https://doi.org/10.1007/s40314-022-02178-6
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s40314-022-02178-6
Keywords
- Reaction–diffusion equation with nonlinear reaction
- High-order compact LOD scheme
- Convergence and stability
- Quenching
- Blow-up