 Research note
 Open Access
 Published:
Accelerated nonstandard finite difference method for singularly perturbed BurgerHuxley equations
BMC Research Notes volume 14, Article number: 446 (2021)
Abstract
Objective
The main purpose of this paper is to present an accelerated nonstandard finite difference method for solving the singularly perturbed BurgerHuxley equation in order to produce more accurate solutions.
Results
The quasilinearization technique is used to linearize the nonlinear term. A nonstandard methodology of Mickens procedure is used in the spatial direction and also within the first order temporal direction that construct the firstorder finite difference approximation to solve the considered problem numerically. To accelerate the rate of convergence from first to secondorder, the Richardson extrapolation technique is applied. Numerical experiments were conducted to support the theoretical results.
Introduction
Singularly perturbed differential equations are typically characterized by a small perturbation parameter multiplied with the highest order derivative term. These types of equations appear in computational fluid dynamics, hydrodynamics, chemical reactor theory, financial modeling, mathematical biology [1,2,3]. The solutions of singularly perturbed problems exhibit layer. The conventional numerical methods on uniform meshes do not produce satisfactory numerical approximations for small values of the perturbation parameter. A uniformly convergent numerical method, which is a numerical method suitable for these problems and in which the error bound is independent of the size of the perturbation parameter (see [4, 5] and the references therein for more details).
In this work the singularly perturbed BurgerHuxley equations of the form:
where \(0 < \varepsilon < < 1\) is a perturbation parameter. The solution domain \(D = (0,1) \times (0,T]\), and \(\alpha \ge 1\), \(\beta \ge 0\), \(\gamma \in (0,1)\) are given constants.
BurgerHuxley equation describes the interaction between convection, diffusion, and reaction processes that have numerous fascinating phenomena such as busting oscillation, population genetics, bifurcation, and so on [6,7,8,9,10,11]. A numerical method for Eq. (1) has been developed only in [6] using a robust adaptive grid method, a uniformly convergent method, and a singular perturbation approach respectively. However, all these methods considered the standard backward Euler scheme. To solve Eq. (1) the idea of using classical numerical methods is not an efficient approach to produce an accurate solution. This should be develop an accelerated nonstandard finite difference method. As a general rule, higherorder convergent methods are preferred as they provide well numerical approximations with low computational cost. So, in this paper, our objective is to propose the accelerated nonstandard finite difference method for solving Eq. (1) by applying the nonstandard procedures of Mickens [12] in the space direction. We provide an error analysis for the method and prove that it is uniformly convergent with secondorder accuracy after applying the Richardson extrapolation technique. It is also shown that the method is computationally more efficient compared to some existing methods in the literature.
Main text
Numerical method
Consider equation under consideration in Eq. (1), rewritten as:
where \(F\left( {x,t,u,\frac{\partial u}{{\partial x}}} \right) = \beta (1  u)(u  \gamma )u  \alpha u\frac{\partial u}{{\partial x}}\).
Let us consider the homogenous part of Eq. (2):
We look for a solution to the dimensionless heat equation of the form of Eq. (3), when \(\varepsilon = 1\), using separation of variables, we get the solution to Eq. (3) as:
that satisfies both the initial and boundary conditions given in Eq. (2). The formula in Eq. (4) will be used to guess intial approximation in the linearization process. Thus, to linearize Eq. (2), by applying the quasilinearization technique on the nonlinear term, for the reasonable initial guess of the form of Eq. (4) is given by:
Thus, the nonlinear term \(F\left( {x,t,u,\frac{\partial u}{{\partial x}}} \right)\) can be linearized initially as:
Substituting Eq. (6) into Eq. (2) and inducing for iteration number \(i\), we obtain the linearized differential equation.
where \(\begin{aligned}a^{(i)} (x,t) &=  \left. {\,\,\frac{\partial F}{{\partial \left( {\frac{\partial u}{{\partial x}}} \right)}}} \right_{{\left( {x,\,t,\,u^{(i)} ,\,\,\frac{{\partial u^{(i)} }}{\partial x}} \right)}} , \ldots b^{(i)} (x,t) \\ & \quad =  \left. {\,\frac{\partial F}{{\partial u}}} \right_{{\,\left( {x,\,t,\,u^{(i)} ,\,\,\frac{{\partial u^{(i)} }}{\partial x}} \right)}} ,\end{aligned}\) \(\begin{aligned}f^{(i)} (x,t) &= F\left( {x,t,u^{(i)} ,\frac{{\partial u^{(i)} }}{\partial x}} \right) \\ & \quad  \,\left. {\,u^{(i)} \frac{\partial F}{{\partial u}}} \right_{{\left( {x,\,t,\,u^{(i)} ,\,\,\frac{{\partial u^{(i)} }}{\partial x}} \right)}}  \left. {\,\,\frac{{\partial u^{(i)} }}{\partial x}\frac{\partial F}{{\partial \left( {\frac{\partial u}{{\partial x}}} \right)}}} \right_{{\left( {x,\,t,\,u^{(i)} ,\,\,\frac{{\partial u^{(i)} }}{\partial x}} \right)}}\end{aligned}\).
Let N be a positive integer different from one, then discretize the interval \([0,\,T]\) on the temporal direction with uniform step length k. Hence, the interval \([0,\,T]\) is partitioned into N equal subintervals with each nodal point satisfies: \(0 = t_{0} < t_{1} < \, \cdots \, < \,t_{N} = T.\) Thus, the temporal nodal points are generated by
Using Taylor’s series expansion, we have:
\(u^{n + 1} (x) = u^{n} (x) + k\frac{{\partial u^{n} (x)}}{\partial t} + \frac{{k^{2} }}{2}\frac{{\partial^{2} u^{n} (x)}}{{\partial t^{2} }} + \frac{{k^{3} }}{6}\frac{{\partial^{3} u^{n} (x)}}{{\partial t^{3} }} + \cdots\), which gives:
where the truncation term is \(TE =  \frac{k}{2}\frac{{\partial^{2} u^{n} (x)}}{{\partial t^{2} }}\). This indicates, the error estimate:
for the arbitrary constant \(C = \frac{1}{2}\left\ {\frac{{\partial^{2} u^{n} (x)}}{{\partial t^{2} }}} \right\_{\infty }\). Substituting Eq. (9) into Eq. (7) gives the semidiscrete differential equation of the form:
Let us consider the singularly perturbed convection–diffusion differential equation of the form:
where \(\eta > 0\) and \(\theta \ge 0\) are constants. Then, Eq. (12) has two linearly independent solutions as:
for the constants \(\lambda_{1} = \frac{{  \eta + \sqrt {\eta^{2} + 4\varepsilon \theta } }}{2\varepsilon },\quad \,{\text{and}}\quad \lambda_{2} = \frac{{  \eta  \sqrt {\eta^{2} + 4\varepsilon \theta } }}{2\varepsilon }\). Let M be is a positive integer, and the interval \([0,\,1]\) on the spatial direction be divided into M equal subintervals through the nodes \(x_{m} = mh,\quad m = \frac{1}{M},\,\quad \,m = 0,\,1,\, \ldots ,M\), and the approximate solution to \(u(x)\) at the grid points \(x_{m} ^{\prime}s\) is denoted by \(u_{m} ^{\prime}s\). The theory of difference (see [12]), shows that the secondorder linear difference equation:
Substituting the values of \(\lambda_{1}\) and \(\lambda_{2}\) into Eq. (13), we get:
This Eq. (15) is the exact difference scheme of Eq. (12), in the wisdom that the difference Eq. (15) has the general solution \(u_{m} = C_{1} \exp (\lambda_{1} x_{m} ) + C_{2} \exp (\lambda_{2} x_{m} )\) as the differential Eq. (12). Since, \(\theta \ge 0\), we require the exact scheme for the reduced case. To this end, the exact scheme corresponding to \(\varepsilon \frac{{d^{2} u(x)}}{{dx^{2} }} + \eta \frac{du(x)}{{dx}} = 0.\,\,\) Thus, Eq. (15) is deduced to:
Multiplying both sides of Eq. (16) by \(\exp \left( {\frac{\eta h}{{2\varepsilon }}} \right)\), and incorporating the term \(u_{m + 1}  u_{m}\) into this equation, we get:
Consequently, Eq. (17) can be transformed into:
where \(\phi^{2} = \frac{h\varepsilon }{\eta }\left( {\exp (\frac{\eta h}{\varepsilon })  1} \right)\).
When we come to the BurgerHuxley equation, the differential equation under consideration the full discretization of \(\overline{D}_{M}^{N}\), denoting to approximation \(u(x_{m} ,t_{n} )\) by \(U_{m}^{n}\). Then, the nonstandard finite difference rules developed for ODE above can be extended for PDE is given by:
where \(i = 0,1,\,\,...\) iteration number, and \(\phi_{m}^{2} = \frac{h\varepsilon }{{a_{m}^{n + 1} }}\left( {\exp \left( {\frac{{ha_{m}^{n + 1} }}{\varepsilon }} \right)  1} \right)\).
Stability of the scheme
A partial differential equation is wellposed if the solution of the partial differential equation is exists, and depends continuously on the initial condition and boundary conditions. The Von Neumann stability technique is applied to investigate the stability of the developed scheme in Eq. (19), by assuming that the solution of Eq. (19) at the grid point \(\left( {x_{m} ,t_{n} } \right)\) is given by:
where \(i = \sqrt {  1} ,\,\,\,\theta\) is the real number and \(\xi\) is the amplitude factor.
Now, considering at the first iteration and putting Eq. (20) into the homogeneous scheme part of Eq. (19) gives:
Then, solving for the amplitude factor \(\xi\), yields:
The condition of stability is \(\xi \,\,\, \le 1\) and for sufficiently small k, we have \(\xi \,\,\, = 1\). Hence, the scheme given in Eq. (19) is stable. Thus, the scheme in Eq. (19) is unconditionally stable.
Consistency of the scheme
The local truncation error \(T(h,k)\) between the operator on the exact solution \(u_{m}^{n}\) to Eq. (7) and the approximate solution \(U_{m}^{n}\) to Eq. (19) at the fixed \(i = 0\) iteration is given by:
Using Taylor’s series expansion, we have:
Substituting Eqs. (22)–(24) into Eq. (21), which implies:
Now, the values of \(\phi_{m}^{2}\) provided in Eq. (19) can be expanded as:
Thus, the truncation error is
Therefore, the norm of the local truncation error can be written as:
where \(C_{1} = \frac{1}{2}\left\ {\frac{{\partial^{2} u_{m}^{n + 1} }}{{\partial t^{2} }}} \right\_{\infty } ,\,\,\,\,{\text{and}}\,\,C_{2} = \frac{1}{2}\left\ {a_{m}^{n + 1} \frac{{\partial^{2} u_{m}^{n + 1} }}{{\partial x^{2} }}} \right\_{\infty } .\)
To accelerate the rate of convergence, we apply the Richardson extrapolation techniques on Eq. (27). Assume that \(U(h,k)\) denote the approximate value of \(u(x_{m} ,t_{n} )\) with the mesh length of h and k. The approximate solution \(U\left( {\frac{h}{2},\frac{k}{2}} \right)\) also denotes the value of \(u(x_{m} ,t_{n} )\) obtained by using the same method with step length \(\frac{h}{2}\) and \(\frac{k}{2}\), then the order of convergence for the two approximate solutions can be written as:
where C is a constant independent of the perturbation and mesh parameters.
Eliminating C from Eq. (28), gives:
Let us denote the combination of the two approximate solutions in Eq. (29) by:
Then, it is rewritten as:
This indicates the method is accelerated to secondorder. Hence, we have
Therefore, the proposed method is secondorder uniformly convergent. Thus, the right hand side hand of Eq. (27) vanishes as \(k \to 0\,\,{\text{and}}\,\,h \to 0\) implies \(T(h,k) \to 0\). Hence, the scheme is consistent with the order of convergence \(O\left( {\mathop k\nolimits^{2} + \mathop h\nolimits^{2} } \right)\). Ended, the scheme developed in Eq. (19), is convergent. A consistent and stable finite difference method is convergent by Lax's equivalence theorem [2]. For further convergence analysis, one can see the works provided in [5, 11, 13, 14].
Numerical illustration and discussions
The exact solution for the considered examples is not available. Hence, the maximum absolute errors are calculated by the double mesh principle, [6], for before and after applying the Richardson extrapolation technique respectively by:
where \(U_{m}^{n}\) and \(U_{2m}^{2n}\) are approximate solutions evaluated on \(\overline{D}_{M}^{N}\) and \(\overline{D}_{2M}^{2N}\) respectively. Similarly, its extrapolated are induced. The corresponding rate of convergences is determined by:
Example 1.
Consider the singularly perturbed BurgersHuxley equation:
Example 2.
Consider the following singularly perturbed Burgers’ equation:
Tables 1 and 2 show the maximum absolute errors that demonstrate the validity of the present method and errors are monotonically decreasing behavior with increasing the number of intervals which confirm the convergence of the method. Table 3 validate that the corresponding rate of convergence. Thus, the proposed method is secondorder convergent. Furthermore, the method gives a more accurate solution than some existing methods in the literature.
Conclusion
The accelerated nonstandard finite difference method is presented for solving the singularly perturbed BurgerHuxley equation. The nonlinear terms are linearized by the quasilinearization technique. Firstorder finite difference approximation for the discretization of a time derivative and the nonstandard methodology of Micken’s procedure is applied for the spatial derivatives. To accelerate the convergence of the method, the Richardson extrapolation technique is applied. It is provided that from numerical results, the method gives a better accurate solution with a higher order of convergence than some existing methods in the literature. Therefore, the presented method is secondorder convergent and gives an accurate solution for solving the singularly perturbed BurgerHuxley equation.
Limitations
During the quasilinearization process we used only the first iteration. If more number of iterations were done, the proposed scheme can have more accurate solution than the existing results. Further, the scheme can more illustrate the physical behaviour of the problem under consideration.
Availability of data and materials
No additional data is used for this research work.
Abbreviations
 ODE:

Ordinary differential equation
 TE:

Truncation errors
References
 1.
Bullo T, Duressa GF, Degla G. Accelerated fitted operator finite difference method for singularly perturbed parabolic reactiondiffusion problems. Computat Methods Diff Eq. 2021;9(3):886–98.
 2.
Bullo TA, Degla GA, Duressa GF. Uniformly convergent higherorder finite difference scheme for singularly perturbed parabolic problems with nonsmooth data. J Appl Math Comput Mech. 2021;20(1):5–16.
 3.
Mekonnen TB, Duressa GF. A computational method for singularly perturbed twoparameter parabolic convectiondiffusion problems. Cogent Math Stat. 2020;7(1):1829277.
 4.
Bullo TA, Duressa GF, Degla GA. Robust finite difference method for singularly perturbed twoparameter parabolic convectiondiffusion problems. Int J Comput Methods. 2021;18(2):2050034.
 5.
Woldaregay MM, Duressa GF. A uniformly convergent numerical method for singularly perturbed delay parabolic differential equations arising in computational neuroscience. Kragujevac J Math. 2022;46(1):65–84.
 6.
Liu LB, Liang Y, Zhang J, Bao X. A robust adaptive grid method for singularly perturbed BurgerHuxley equations. Elect Res Arch. 2020;28(4):1439.
 7.
Mohanty RK, Dai W, Liu D. Operator compact method of accuracy two in time and four in space for the solution of timedependent BurgersHuxley equation. Numerical Algorithm. 2015;70(3):591–605.
 8.
Mohammadi R. Bspline collocation algorithm for the numerical solution of the generalized Burger’sHuxley equation. Numer Methods Partial Diff Equ. 2013;29(4):1173–91.
 9.
Khattak AJ. A computational meshless method for the generalized Burger’s–Huxley equation. Appl Math Model. 2009;33(9):3718–29.
 10.
Mittal RC, Tripathi A. Numerical solutions of generalized BurgersFisher and generalized BurgersHuxley equations using collocation of cubic Bsplines. Int J Comput Math. 2015;92(5):1053–77.
 11.
Kaushik A, Sharma MD. A uniformly convergent numerical method on nonuniform mesh for singularly perturbed unsteady BurgerHuxley equation. Appl Math Comput. 2008;195(2):688–706.
 12.
Mickens RE. Advances in the applications of nonstandard finite difference schemes. World Scientific; 2005.
 13.
Jima KM, File DG. Implicit Finite Difference Scheme for Singularly Perturbed BurgerHuxley Equations. 2022;35(1):87–100.
 14.
Kabeto MJ, Duressa GF. Robust numerical method for singularly perturbed semilinear parabolic differential difference equations. Math Comput Simul. 2021;188:537–47.
Acknowledgements
The authors thank Jimma University for their material supports.
Funding
No funding organization for this research work.
Author information
Affiliations
Contributions
MJK carried out scheme development, MATLAB coding and numerical experimentation. GFD formulated the problem, design, and draft the manuscript. Both authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Kabeto, M.J., Duressa, G.F. Accelerated nonstandard finite difference method for singularly perturbed BurgerHuxley equations. BMC Res Notes 14, 446 (2021). https://doi.org/10.1186/s13104021058584
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13104021058584
Keywords
 Singularly perturbed BurgerHuxley equation
 Accelerated nonstandard method
 Accurate solution
Mathematics Subject Classification
 Primary 65M06
 65M12
 Secondary 65M15