Abstract
In this work, we investigate numerically a diffusion equation with nonlinear reaction, defined spatially over a closed and bounded interval of the real line. The partial differential equation is expressed in an equivalent logarithmic form, and initial and Dirichlet boundary data are imposed upon the problem. An implicit finite-difference discretization of this logarithmic model is proposed then. We show that the numerical scheme is capable of preserving the constant solutions of the continuous model. Moreover, we establish the existence of positive and bounded numerical solutions using analytical arguments. The theoretical analysis of the numerical model is carried out also. In particular, we establish that the method is a consistent technique, we provide a priori bounds for the numerical solutions, and we prove the stability and the convergence of the scheme by applying a suitable Gronwall-type inequality. As one of the consequences of stability, we show not only that the solutions of the numerical model exist, but also that they are unique.
Similar content being viewed by others
References
K. Al-Khaled, Numerical study of Fisher’s reaction–diffusion equation by the Sinc collocation method. J. Comput. Appl. Math. 137(2), 245–255 (2001)
K. Baba, M. Tabata, On a conservation upwind finite element scheme for convective diffusion equations. RAIRO-Analyse Numérique 15(1), 3–25 (1981)
C. Bo-Kui, M. Song-Qiang, W. Bing-Hong, Exact solutions of generalized Burgers–Fisher equation with variable coefficients. Commun. Theor. Phys. 53(3), 443 (2010)
G. Carey, Y. Shen, Least-squares finite element approximation of Fisher’s reaction–diffusion equation. Numer. Methods Partial Differ. Equ. 11(2), 175–186 (1995)
P. Das, S. Rana, H. Ramos, Homotopy perturbation method for solving caputo type fractional order Volterra–Fredholm integro-differential equations. Comput. Math. Methods 1, e1047 (2019)
X. Deng, Travelling wave solutions for the generalized Burgers–Huxley equation. Appl. Math. Comput. 204(2), 733–737 (2008)
S.A. El Morsy, M.S. El-Azab, Logarithmic finite difference method applied to KdVB equation. Am. Acad. Sch. Res. J. (AASRJ) 4(2), 41–48 (2012)
V. Ervin, J. Macías-Díaz, J. Ruiz-Ramírez, A positive and bounded finite element approximation of the generalized Burgers–Huxley equation. J. Math. Anal. Appl. 424(2), 1143–1160 (2015)
R.A. Fisher, The wave of advance of advantageous genes. Ann. Eug. 7(4), 355–369 (1937)
D. Furihata, Finite-difference schemes for nonlinear wave equation that inherit energy conservation property. J. Comput. Appl. Math. 134(1), 37–57 (2001)
D. Furihata, T. Matsuo, Discrete Variational Derivative Method: A Structure-Preserving Numerical Method for Partial Differential Equations (CRC Press, New York, 2010)
G. Hariharan, K. Kannan, K. Sharma, Haar wavelet method for solving Fisher’s equation. Appl. Math. Comput. 211(2), 284–292 (2009)
A.L. Hodgkin, A.F. Huxley, A quantitative description of membrane current and its application to conduction and excitation in nerve. J. Physiol. 117(4), 500–544 (1952)
B. İnan, A logarithmic finite difference technique for numerical solution of the generalized Huxley equation, in Book of Abstracts, 7th International Eurasian Conference on Mathematical Sciences and Applications, Kiev, Ukraine (2018), pp. 100–101
E.W. Jenkins, C. Paribello, N.E. Wilson, Discrete mass conservation for porous media saturated flow. Numer. Methods Partial Differ. Equ. 30(2), 625–640 (2014)
A.N. Kolmogorov, Étude de l’équation de la diffusion avec croissance de la quantité de matière et son application à un problème biologique. Bull. Univ. Moskow, Ser. Internat. Sec. A 1, 1–25 (1937)
J.A. López-Campos, A. Segade, E. Casarejos, J.R. Fernández, Behavior characterization of viscoelastic materials for the finite element method calculation applying Prony series. Comput. Math. Methods 1(1), e1014 (2019)
J. Macías-Díaz, A. Puri, On the propagation of binary signals in damped mechanical systems of oscillators. Physica D 228(2), 112–121 (2007)
J.E. Macías-Díaz, Numerical study of the transmission of energy in discrete arrays of sine-Gordon equations in two space dimensions. Phys. Rev. E 77(1), 016,602 (2008)
J.E. Macías-Díaz, Sufficient conditions for the preservation of the boundedness in a numerical method for a physical model with transport memory and nonlinear damping. Comput. Phys. Commun. 182(12), 2471–2478 (2011)
J.E. Macías-Díaz, A Mickens-type monotone discretization for bounded travelling-wave solutions of a Burgers–Fisher partial differential equation. J. Differ. Equ. Appl. 19(11), 1907–1920 (2013)
J.E. Macías-Díaz, Conciliating efficiency and dynamical consistency in the simulation of the effects of proliferation and motility of transforming growth factor \(\beta\) on cancer cells. Commun. Nonlinear Sci. Numer. Simul. 40, 173–188 (2016)
J.E. Macías-Díaz, A structure-preserving method for a class of nonlinear dissipative wave equations with Riesz space-fractional derivatives. J. Comput. Phys. 351, 40–58 (2017)
J.E. Macías-Díaz, Persistence of nonlinear hysteresis in fractional models of Josephson transmission lines. Commun. Nonlinear Sci. Numer. Simul. 53, 31–43 (2017)
J.E. Macías-Díaz, Numerical simulation of the nonlinear dynamics of harmonically driven Riesz-fractional extensions of the Fermi–Pasta–Ulam chains. Commun. Nonlinear Sci. Numer. Simul. 55, 248–264 (2018)
J.E. Macías-Díaz, On the numerical and structural properties of a logarithmic scheme for diffusion–reaction equations. Appl. Numer. Math. 140, 104–114 (2019)
J.E. Macías-Díaz, F.J. Avelar-González, R.S. Landry, On an efficient implementation and mass boundedness conditions for a discrete Dirichlet problem associated with a nonlinear system of singular partial differential equations. J. Differ. Equ. Appl. 21(11), 1021–1043 (2015)
J.E. Macías-Díaz, A. González, A convergent and dynamically consistent finite-difference method to approximate the positive and bounded solutions of the classical Burgers–Fisher equation. J. Comput. Appl. Math. 318, 604–615 (2017)
J.E. Macías-Díaz, A.S. Hendy, R.H. De Staelen, A compact fourth-order in space energy-preserving method for Riesz space-fractional nonlinear wave equations. Appl. Math. Comput. 325, 1–14 (2018)
J.E. Macías-Díaz, A.S. Hendy, R.H. De Staelen, A pseudo energy-invariant method for relativistic wave equations with Riesz space-fractional derivatives. Comput. Phys. Commun. 224, 98–107 (2018)
J.E. Macías-Díaz, B. İnan et al., Structural and numerical analysis of an implicit logarithmic scheme for diffusion equations with nonlinear reaction. Int. J. Mod. Phys. C (IJMPC) 30(09), 1–14 (2019)
J.E. Macías-Díaz, I.E. Medina-Ramírez, A. Puri, Numerical treatment of the spherically symmetric solutions of a generalized Fisher–Kolmogorov–Petrovsky–Piscounov equation. J. Comput. Appl. Math. 231(2), 851–868 (2009)
J.E. Macías-Díaz, A. Puri, A boundedness-preserving finite-difference scheme for a damped nonlinear wave equation. Appl. Numer. Math. 60(9), 934–948 (2010)
J.E. Macías-Díaz, A. Puri, An explicit positivity-preserving finite-difference scheme for the classical Fisher–Kolmogorov–Petrovsky–Piscounov equation. Appl. Math. Comput. 218(9), 5829–5837 (2012)
F. Martínez, I. Martínez-Vidal, S. Paredes, Conformable Euler’s Theorem on homogeneous functions. Comput. Math. Methods 1, e1048 (2019)
T. Matsuo, D. Furihata, Dissipative or conservative finite-difference schemes for complex-valued nonlinear partial differential equations. J. Comput. Phys. 171(2), 425–447 (2001)
R.E. Mickens, Dynamic consistency: a fundamental principle for constructing nonstandard finite difference schemes for differential equations. J. Differ. Equ. Appl. 11, 645–653 (2005)
R.E. Mickens, Dynamic consistency: a fundamental principle for constructing nonstandard finite difference schemes for differential equations. J. Differ. Equ. Appl. 11(7), 645–653 (2005)
R. Mittal, R. Jiwari, Numerical study of Fisher’s equation by using differential quadrature method. Int. J. Inf. Syst. Sci. 5(1), 143–160 (2009)
R.C. Mittal, R.K. Jain, Numerical solutions of nonlinear Fisher’s reaction–diffusion equation with modified cubic B-spline collocation method. Math. Sci. 7(1), 12 (2013)
M.D. Morales-Hernández, I.E. Medina-Ramírez, F.J. Avelar-González, J.E. Macías-Díaz, An efficient recursive algorithm in the computational simulation of the bounded growth of biological films. Int. J. Comput. Methods 9(04), 1250,050 (2012)
D. Prakasha, P. Veeresha, H.M. Baskonus, Two novel computational techniques for fractional Gardner and Cahn–Hilliard equations. Comput. Math. Methods 1(2), e1021 (2019)
H. Ramos, Development of a new Runge–Kutta method and its economical implementation. Comput. Math. Methods 1(2), e1016 (2019)
M. Ran, C. Zhang, A linearly implicit conservative scheme for the fractional nonlinear Schrödinger equation with wave operator. Int. J. Comput. Math. 93(7), 1103–1118 (2016)
J. Ruiz-Ramírez, J.E. Macías-Díaz, A finite-difference scheme to approximate non-negative and bounded solutions of a FitzHugh–Nagumo equation. Int. J. Comput. Math. 88(15), 3186–3201 (2011)
R. Santos, J.P. Martins, M. Felgueiras, S. Ferreira, Revisiting random polygonal lines iteratively generated on the plane. Comput. Math. Methods 1(4), e1036 (2019)
W. Tang, J. Zhang, Symplecticity-preserving continuous-stage Runge–Kutta–Nyström methods. Appl. Math. Comput. 323, 204–219 (2018)
T. Wang, B. Guo, Analysis of some finite difference schemes for two-dimensional Ginzburg–Landau equation. Numer. Methods Partial Differ. Equ. 27(5), 1340–1363 (2011)
X. Wang, Z. Zhu, Y. Lu, Solitary wave solutions of the generalised Burgers–Huxley equation. J. Phys. A Math. Gen. 23(3), 271 (1990)
B. Wongsaijai, T. Mouktonglang, N. Sukantamala, K. Poochinapan, Compact structure-preserving approach to solitary wave in shallow water modeled by the Rosenau-RLW equation. Appl. Math. Comput. 340, 84–100 (2019)
Acknowledgements
The corresponding author wishes to thank the anonymous reviewers for their comments and criticisms. All of their comments were taken into account in the revised version of the paper, resulting in a substantial improvement with respect to the original submission. The first author wishes to acknowledge the financial support from the National Council for Science and Technology of Mexico (CONACYT) through the Grant A1-S-45928.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Macías-Díaz, J.E., Hendy, A.S. On the stability and convergence of an implicit logarithmic scheme for diffusion equations with nonlinear reaction. J Math Chem 58, 735–749 (2020). https://doi.org/10.1007/s10910-020-01103-8
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10910-020-01103-8
Keywords
- Reaction–diffusion equations
- Implicit logarithmic scheme
- Structure-preserving method
- Numerical efficiency analysis