Skip to main content


Front. Appl. Math. Stat., 13 February 2023
Sec. Numerical Analysis and Scientific Computation
Volume 9 - 2023 |

A NSFD method for the singularly perturbed Burgers-Huxley equation

  • 1Department of Mathematics, Adama Science and Technology University, Adama, Ethiopia
  • 2Department of Mathematics and Applied Mathematics, University of the Western Cape, Bellville, South Africa

This article focuses on a numerical solution of the singularly perturbed Burgers-Huxley equation. The simultaneous presence of a singular perturbation parameter and the nonlinearity raise the challenge of finding a reliable and efficient numerical solution for this equation via the classical numerical methods. To overcome this challenge, a nonstandard finite difference (NSFD) scheme is developed in the following manner. The time variable is discretized using the backward Euler method. This gives rise to a system of nonlinear ordinary differential equations which are then dealt with using the concept of nonlocal approximation. Through a rigorous error analysis, the proposed scheme has been shown to be parameter-uniform convergent. Simulations conducted on two numerical examples confirm the theoretical result. A comparison with other methods in terms of accuracy and computational cost reveals the superiority of the proposed scheme.

1. Introduction

The time-dependent singularly perturbed Burgers-Huxley equation models a large class of physical phenomena such as the interaction between convection effect, reaction mechanism and diffusion transport. In one dimensional space this equation has the form [1, 2]

{ut-εuxx+a(u)ux+b(u)u=0, (x,t)DΩ×(0,T](0,1)×(0,T],u(x,0)=u0(x), xΩ,u(0,t)=h1(t), u(1,t)=h2(t), t[0,T],    (1.1)


a(u)=αu,  b(u)=-β(1-u)(u-γ),

α ≥ 0, β ≥ 0, γ ∈ (0, 1) are real parameters and u0(x), h1(t),and h2(t) are sufficiently smooth prescribed functions. Here, 0 < ε ≪ 1 denotes the perturbation parameter.

The field of numerical methods for singularly perturbed problems has flourished significantly over the last couple of decades. The studies focused on the design and analysis of numerical methods that were parameter-uniformly convergent. In other words, research was mostly motivated by the challenge posed by the presence of the singular perturbation parameter. A vast majority of researchers discussed linear problems. Nonlinear singularly perturbed equations, and in particular Equation (1.1), received little attention from scholars in this research area. Nonlinearity constitutes an another layer of difficulties of singularly perturbed problems.

In general, the solution of the singularly perturbed Burgers-Huxley problem (1.1) with arbitrary initial and boundary conditions (IBCs) cannot be expressed in terms of a finite number of elementary functions [2, 3]. Thus, several scholars seek approximate analytic solutions using different transformation techniques [46]. However, these techniques still hold for only some specific parameters and initial boundary conditions. Hence, for general cases, there is a need to develop numerical methods to obtain an approximate solution of the problem.

For the last two decades, various numerical methods have been studied to solve the Burgers-Huxley equations of type (1.1) for ε = 1 and specific IBCs in the framework of NSFD methods. For example, A.R. Appadu et al. [7] developed two novel nonstandard finite difference schemes, and explicit exponential finite difference method and a fully implicit exponential finite difference method.

For the singularly perturbed case, that is when 0 < ε ≪ 1, various parameter-uniformly convergent methods have been developed to solve (1.1) [1, 2, 8, 9]. It is to be noted that all these works used the Newton's quasilinearization to deal with the nonlinearity.

In this article, a nonstandard finite difference method for the singularly perturbed case is proposed. To deal with the nonlinearity, rather than using quasilinearization, for the first time, nonlinear terms are approximated in a nonlocal manner following one of Micken's rules [10]. The resulting method preserves the properties of the continuous solution and provides reliable numerical results. The method is proved to be first order parameter-uniform convergent in time and space.

The rest of this article is structured in the following manner. In the next section, we study a priori estimates of the analytic solution to the problem. Section 3 is about the proposed NSFD method and its parameter-uniform convergence. Section 4 deals with the implementation of the method to confirm the theoretical results and compare with other methods. Section 5 concludes the present work and provides a direction for future work.

2. Some analytical results: A priori estimates

Throughout this article, we assume that u0, h1, and h2 are sufficiently smooth functions and that a and b satisfy

a(u)p>0, b(u)q0.    (2.1)

Under these assumptions the problem (1.1) has a unique solution u(x,t)C2,1(D̄) which exhibits a boundary layer near x = 1 [1].

In the rest of the paper, C denotes a generic constant independent of the parameters and mesh sizes, and ||.||D denotes the maximum norm over D. For the solution u(x, t) of (1.1), we have the following bounds.

Lemma 2.1. [Uniform stability estimate for continuous problem]

Let u(x, t) be the exact solution of (1.1) on D̄, then


where Di = {(x, t):t = 0, x ∈ [0, 1]}, ∂D = DiDLDR, DL = {(x, t):x = 0, t ∈ [0, T]}, and DR = {(x, t):x = 1, t ∈ [0, T]}.

Proof: The proof can be seen in [9].

Remark 2.2. From our assumption in (2.1) and Lemma 2.1, one can observe that the solution u(x, t) satisfies

pαu(x,t)C||u0||Di+||u||D, (x,t)D̄.    (2.2)

3. Proposed scheme and its convergence analysis

This section is dedicated to the construction of a new scheme and to the analysis of its parameter-uniform convergence. The first step is to discretize the solution domain D̄ and provide the definition of the NSFD methods given in [11] (a revised version of [10]).

Let M and N be positive integers and D̄M,N=Ω̄M×[0,T]N be uniform grid discretization of the solution domain D̄=Ω̄×[0,T] such that




where h = 1/M and Δt = T/N are spatial and temporal step sizes, respectively.

Definition 3.1. A discrete scheme to determine approximate solutions umn to the solution u(x, t) of the problem (1.1) is called a NSFD method if at least one of the following conditions is satisfied [10]:

1. The classical denominator h or h2 of the discrete first or second order derivative is replaced by a nonnegative function ψ such that

ψ(h)=h+O(h)  or  ψ2(h)=h2+O(h2).

For example, denominator functions that satisfy the above conditions are

ψ(h)=h, sin(h), eβh-1β, hεa(exp(ahε)-1)    (3.1)

and so on.

2. Nonlinear terms that occur in the differential equation are approximated in a nonlocal way. For example,

umn+1umn,   (umn)2umnvm+1n,  (umn)2umnum-1n+umn+vm+1n3,(umn+1)32(umn)3-(umn)2umn+1, (umn)3um-1numnum+1n,    (3.2)

and so on.

Definition 3.2. Assume that the solution of (1.1) satisfies some property P. The difference equation of (1.1) in umn is called (qualitatively) stable with respect to P if, for any values of the mesh sizes Δt and h, solution of the difference equation replicates the property P.

Remark 3.3. In [10], Mickens set five rules for the constructions of the finite difference models that can replicate the properties of the exact solution. Definition 3.1 satisfies only two of these rules.

The remaining rules are expressed in terms of definition 3.2. For example, the schemes under consideration in this paper is qualitatively stable with respect to the order of the differential equation, they do satisfy positivity and uniform boundedness.

To construct the scheme, first we semidiscretize the problem (1.1) in time direction and then in space direction.

3.1. Semidiscrete scheme

Denote un(x) as the approximation of u(x, tn) at time level tn, 0 ≤ nN. Now, we apply backward Euler finite difference method to discretize the continuous problem (1.1) in the temporal direction and obtain the following semidiscrete scheme:

{u0(x)=u(x,0)=u0(x),(I+ΔtLεN)un+1(x)=un(x),xΩ,un+1(0)=h1(tn+1), un+1(1)=h2(tn+1),  n=0,1,,N-1,    (3.3)



Here the nonlinearity a(u)ux and b(u)u approximated in a nonlocal way as

a(un+1)dun+1dxa(un)dun+1dx  and b(un+1)un+1b(un)un+1.

and thus this approximation satisfies the condition (2) in definition 3.1.

3.2. Convergence analysis for the semidiscrete scheme

The local truncation error en+1 of the semidiscrete scheme (3.3) is given by en+1=u(x,tn+1)-ûn+1(x), where ûn+1(x) is the solution of

{(I+ΔtLεN)u^n+1=u(x,tn), xΩ,u^n+1(0)=h1(tn+1),  u^n+1(1)=h2(tn+1),  n=0,1,,N-1.    (3.4)

That is, ûn+1(x) is the solution obtained after one step of semidiscrete scheme (3.3) by taking the exact value u(x, tn) instead of un(x) as the starting data.

In order to analyze the uniform convergence of the solution un(x) of (3.3) to the exact solution u(x, tn), we will perform the stability analysis and establish the consistency result. First, let us consider the semidiscrete maximum principle for the operator I+ΔtLεN.

Lemma 3.4. Let Φn+1(x) be a function such that Φn+1(0) ≥ 0, Φn+1(1) ≥ 0 and (I+ΔtLεN)Φn+1(x)0 for all x ∈ Ω. Then Φn+1(x) ≥ 0 for all xΩ̄.

Proof: Suppose there is x*Ω̄ such that Φn+1(x*)=minxΩ̄Φn+1(x)<0. From the given hypothesis and second derivative test, we have x*0,x*1, Φxn+1(x*)=0 and Φxxn+1(x*)>0. Then, from (3.3) we have

(I+ΔtεN)Φn+1(x*)=Φn+1(x*)+Δt( εΦxxn+1+a(Φn)Φxn+1                                           +b(Φn)Φn+1 )(x*)                                         =Φn+1(x*)+Δt( εΦxxn+1(x*)                                          +a(Φn(x*))Φxn+1(x*)+b(Φn(x*))Φn+1(x*) )                                         =Φn+1(x*)+Δt( εΦxxn+1(x*)                                          +b(Φn(x*))Φn+1(x*) ).

Assumption (2.1) leads to (I+ΔtLεN)Φn+1(x*)<0, which contradicts the given assumption in Lemma 3.4 and thus Φn+1(x) ≥ 0 for all xΩ̄.

This maximum principle leads to the following stability result

||(I+ΔtLεN)-1||C.    (3.5)

Lemma 3.5. (Local error estimate)

Estimate of the local error en+1 is given by

||en+1||CΔt2.    (3.6)

Proof: See [2].

The global error En+1 associated to the semidiscrete scheme (3.3) at (n+1)-th time level is given by En+1=u(x,tn+1)-un+1(x)=i=1n+1ei. Using the local error estimate (Lemma 3.5) and triangular inequality the following global error estimate follows.

Theorem 3.6. (Global error estimate)

The global error En+1 of the time discretisation at the n + 1 time step satisfies

||En+1||CΔt,  nΔtT.    (3.7)

Therefore, the semidiscrete scheme (3.3) is a first order uniformly convergent in the time direction.

The following lemma provides the asymptotic estimates of the exact solution un+1 of (3.3) and its derivatives. These estimates will be used in the convergence analysis of the fully discrete scheme

Lemma 3.7. If un+1(x) is the solution of the problem (3.3), then there exists a constant C such that

||iun+1xi||Ω̄C(1+ε-iexp(-p(1-x)ε)), 0i4, xΩ̄.    (3.8)

Proof: See the proof of Lemma 4.1 in [12]

3.3. The spatial discretisation

In this subsection, we totally discretise the semidiscrete scheme (3.3) on a uniform mesh Ω̄M. Let the approximation of un(x) at xm be denoted by umn, 0 ≤ mM. Similarly, let the approximations of a(umn) and b(umn) be denoted respectively by amn and bmn.

Using the upwind finite difference scheme and nonstandard finite difference methodology of Mickens [13], the semidiscrete problem (3.3) can be fully-discretised as

{(I+ΔtLεM,N)umn+1=un(xm),  1mM-1,u0n+1=h1(tn+1),  uMn+1=h2(tn+1),    (3.9)


LεM,Numn+1-ε[um+1n+1-2umn+1+um-1n+1(ψmn)2]+amn[umn+1-um-1n+1h]                     +bmnumn+1,

and the denominator function ψmn is


The denominator function is derived systematically to replicate the dissipativity properties of the solution of the differential equations. Interested readers may refer to [14] for details. Observe that


and thus this function satisfies the condition (1) in definition 3.1.

The method in Eq. (3.9) is a linear nonstandard finite difference (LNSFD) and can be written as:

Emn+1um-1n+1+Fmn+1umn+1+Gmn+1um+1n+1=Hmn+1,  1mM-1,    (3.10)


Emn+1=-Δt(ε(ψmn)2+amnh),    (3.11)
Fmn+1=1+Δt(2ε(ψmn)2+amnh+bmn),    (3.12)
Gmn+1=-εΔt(ψmn)2,    (3.13)

Equation (3.10) leads to the tridiagonal system which can be solved with Thomas Algorithm [15, 16].

Remark 3.8. The developed scheme is a linear difference equation even if the original equation (1.1) is nonlinear. This is one of the feature of NSFD.

3.4. Convergence analysis of the space discretisation

The difference operator I+ΔtLεM,N of (3.9) satisfies the following maximum principle. Hence (3.9) has a unique discrete solution umn.

Lemma 3.9. Let Φmn+1,m=0,1,,M be fully discrete mesh functions. If Φ0n+10, ΦMn+10 and (I+ΔtLεM,N)Φmn+10 for 1 ≤ mM−1 then Φmn+10 for 0 ≤ mM.

Proof: Let us proof by contradiction. Assume there is m* such that Φm*n+1=min0mMΦmn+1<0. Thus, we have, m*∉{0, M}, Φm*+1n+1-Φm*n+10, and Φm*-1n+1-Φm*n+10. Now, from (3.9) we have

(I+ΔtεM,N)Φm*n+1=Φm*n+1+Δt [ ε(Φm*+1n+12Φm*n+1+Φm*1n+1(ψm*n)2)                                       +am*n(Φm*n+1Φm*1n+1h)+bm*nΦm*n+1 ]                                   =εΔt(Φm*+1n+1Φm*n+1+Φm*1n1Φm*n+1(ψm*n)2)                                     +am*nΔt(Φm*n+1Φm*1n+1h)                                     +(1+Δtbm*n)Φm*n+1<0,                                

which contradicts the given assumption (I+ΔtLεM,N)Φmn+10 for 1 ≤ mM−1 and hence Φmn+10 for 0 ≤ mM.

This leads immediately to the following stability result, analogous to the continuous result.

Lemma 3.10. [Uniform stability estimate for discrete problem]

Let Φmn+1, 0mM, be any mesh functions such that Φ0n+1=ΦMn+1=0, then

|Φmn+1|11+qΔtmax1iM-1|(I+ΔtLεM,N)Φin+1|,  0mM.

Proof: Define


which implies

[ Ψ± ]0n+1=11+qΔtmax1iM1|(I+ΔtεM,N)Φin+1|± Φ0n+1               =11+qΔtmax1iM1|(I+ΔtεM,N)Φin+1|  0,[ Ψ± ]Mn+1=11+qΔtmax1iM1|(I+ΔtεM,N)Φin+1|± ΦMn+1               =11+qΔtmax1iM1|(I+ΔtεM,N)Φin+1|  0,

and, for m = 1, 2, ⋯ , M−1,


Since bn(xm)q, one has (I+ΔtLεM,N)[Ψ±]mn+10. Thus, by discrete maximum principle given in Lemma 3.9, one obtains [Ψ±]mn+10. This gives the required result

|Φmn+1|1qmax1iM-1|(I+ΔtLεM,N)Φin+1|,  0mM.

Remark 3.11. Lemmas 3.9 and 3.10 show the developed nonstandard finite difference method replicates the positivity and boundedness of the solution, respectively. Hence, the proposed method is qualitatively stable.

Theorem 3.12. (Error estimate in the spatial direction)

Let un+1(x) be the solution of the semidiscrete problem (3.4) and umn+1 be the solution of the full discretisation (3.9). Then, the error estimate is given by

|un+1(xm)-umn+1|Ch, 0mM.

Proof: The truncation error of the complete discretization (3.9) is given by

(I+ΔtLεM,N)(un+1(xm)-umn+1)    (3.14)
=Δt[-ε(un+1)(xm)+an(xm)(un+1)(xm)+bn(xm)un+1(xm)]-Δt[-εh2(ψmn)2D+D-un+1(xm)+an(xm)D-un+1(xm)+bn(xm)un+1(xm)]=εΔt(h2(ψmn)2D+D--d2dx2)un+1(xm)+an(xm)Δt(ddx-D-)un+1(xm),    (3.15)



Applying absolute values and using triangle inequalities in (3.15) leads to

| (I+ΔtεM,N)(un+1(xm)umn+1) |            εΔt| (h2(ψmn)21)D+Dun+1(xm) |          +εΔt| (D+Dd2dx2)un+1(xm) |           +an(xm)Δt| (ddxD)un+1(xm) |.                      (3.16)

From the use the Taylor series expansions of some terms in the above equation, on obtains

h2(ψmn)2-1=-an(η1)h2ε+ha(η1),    (3.17)
D+D-un+1(xm)=d2un+1dx2(η2),    (3.18)
(D+D--d2dx2)un+1(xm)=h216d4un+1dx4(η3),    (3.19)
(ddx-D-)un+1(xm)=h2d2un+1dx2(η4),    (3.20)

for some ηi such that xm−1 ≤ ηixm+1, i = 1, 2, 3, 4.

Substituting equations (3.17)–(3.20) in (3.16) and from the boundedness of an(x) in (3.8), one has

| (I+ΔtεM,N)(un+1(xm)umn+1) |εΔtan(η1)h2ε+han(η1)|d2un+1dx2(η2)|+εΔth216|d4un+1dx4(η3)|+an(xm)Δth2|d2un+1dx2(η4)|,    (3.21)

and this gives

| (I+ΔtM,N)(un+1(xm)umn+1) |  Ch(1+ε2exp(p(1xm)ε))  +Ch2(ε+ε3exp(p(1xm)ε))  Ch(1+ε3exp(p(1xm)ε)), for m=1,2,,M1.

Applying Lemma 7 of [17] gives

limε0|(I+ΔtLM,N)(un+1(xm)-umn+1)|Ch    (3.22)

which, upon use of Lemma 3.10, leads to

limε0|un+1(xm)-umn+1|Ch.    (3.23)


Combining the schemes (3.3) and (3.9), one obtains the following fully-discrete scheme on the mesh Ω̄M×[0,T]N.

{(I+ΔtLεM,N)umn+1=umn, 1mM-1, 0nN-1,un+1(0)=h1(tn+1),un+1(1)=h2(tn+1), 0nN-1.    (3.24)

The temporal and spatial error estimates in Theorem 3.6 and 3.12, respectively, give the the following main result of this paper.

Theorem 3.13. If u(x, t) is the exact solution of the continuous problem (1.1) and umn is the solution of the fully-discrete scheme (3.24), then

|u(xm,tn)-umn|C(Δt+h), 0mM,0nN.

Therefore, the presented discrete scheme is ε-uniform convergent of order one both in time and space.

4. Numerical implementation

In this section, two test examples are provided to demonstrate the efficiency and applicability of the proposed numerical method. The first example is taken from a recent article [7] for ε = 1, and we modified it by multiplying the highest derivative term by ε to make the problem singularly perturbed. It is the first time to consider this example for the singularly perturbed case. The second example is taken from [2] for different initial and boundary conditions to satisfy our assumption (2.1). All the computations are carried out by Intel Coreå i5-4210M CPU @2.60GHz×4 with MATLAB 2017.

Example 4.1. Consider the following singularly perturbed Burgers-Huxley problem:

{ut-ε2ux2+αuux-β(1-u)(u-γ)u=0,  (x,t)(0,1)×(0,T],    u(x,0)=γ2+γ2tanh(A1x), 0x1,    u(0,t)=γ2+γ2tanh(-A1A2t), u(1,t)=γ2+γ2tanh(A1(1-A2t)),                                                                         0tT,


A1=-α+α2+8β8, A2=γα2-(2-γ)(-α+α2+8β)8.

For ε = 10−4, α = β = 1 and γ = 0.5, we plotted the numerical solution of Example 4.1 using the proposed method in Figure 1. As we see, the method resolves the boundary layer at x = 1. Thus, our proposed method replicates the property of the continuous solution.


Figure 1. Numerical solutions of Example 4.1 using the proposed method for T = 1, M = 64, N = 40 and ε = 10−4. A layer is observed near x = 1.

Since the exact solution of this problem is not known, for each ε, the maximum pointwise error is calculated using double mesh principle [18] given by


where um;Mn;N and u2m;2M2n;2N are approximate solutions to problem (3.9) on DM, N and D2M, 2N, respectively. The corresponding order of convergence is given by


In Tables 1, 2, we compute the maximum absolute errors and the corresponding order of convergences using the proposed method for α = β = 1, γ = 0.5 and various values of ε with fixed values of M and N, respectively. From these results, we can observe that the method is ε-uniform convergent of order one both in time and space. This confirms the theoretical error results.


Table 1. Maximum absolute errors for Example 4.1 for M = 32 at the number of intervals N.


Table 2. Maximum absolute errors for Example 4.1 for N = 10 at the number of intervals M.

Table 3 compares absolute errors at some values of x and t using the proposed method and methods in [7, 19] and the nonlinear nonstandard finite difference (NNSFD) method given by

umn+1-umnϕ(Δt,ε)-εum+1n+1-2umn+1+um-1n+1(ψmn+1)2(h,ε)+amn+1umn+1-um-1n+1h              +bmn+1umn+1=0,    (4.1)

where the denominator functions

ϕ(Δt,ε)=eεΔt-1ε,  (ψmn+1)2(h,ε)=hεamn+1(exp(amn+1hε)-1).

The nonlinear system (4.1) can be easily solved using Newton's method (see in [7]). The solution at the previous time-step is taken as the initial estimate. The methods considered in [7] and [19] are NSFD methods and variational iteration method (VIM), respectively. From the table, we observe that our proposed method gives more accurate results.


Table 3. Comparison of absolute errors for Example 4.1 using the proposed scheme (LNSFD) with results of [7, 19] and NNSFD for ε = 1, α = β = 1, γ = 0.001 at some values of x and t.

Also, in Table 4, we compare the CPU times required to compute the solutions of Example 4.1 for M = 10, N = 2000, α = β = 1, γ = 0.001, T = 10 using the LNSFD and NNSFD methods. As one can observe in the table, the CPU time for NNSFD is larger than the LNSFD.


Table 4. CPU time for Example 4.1.

Example 4.2. Consider the following SPBHE [2]:

{ut-ε2ux2+uux-(1-u)(u-0.5)u=0,  (x,t)(0,1)×(0,1],    u(x,0)=1+sin(πx/2), 0<x<1,    u(0,t)=1,u(1,t)=2, 0t1.

5. Conclusion

In this paper, a nonstandard finite difference method for a singularly perturbed Burgers-Huxley equation has been developed. First, the backward-Euler scheme was applied to discretize the problem (1.1) with respect to time derivative and the upwind nonstandard finite difference scheme on uniform mesh to approximate the spatial derivative. Then, the presented method was proved to be first-order convergent in both the spatial and temporal variables. Numerical results are given in Figures 1, 2 and in Tables 16 for two test examples to confirm the theoretical results and to compare with recent results. It has been observed from these figures and tables that the numerical results are in agreement with the theoretical findings. For ε = 1, in Table 3, comparisons with the NNSFD and the VIM of [19] and the NSFD method of [7] reveal that the proposed NSFD method gives more accurate results. In addition, the present method is also applicable when 0 < ε < 1. Thus, in all cases, the present method produces more accurate results than the existing schemes.


Figure 2. Numerical solutions of Example 4.2 using the proposed method for T = 1, M = 64, N = 40 and ε = 10−4. A layer is observed near x = 1.


Table 5. Maximum absolute errors for Example 4.2 for M = 32 at the number of intervals N.


Table 6. Maximum absolute errors for Example 4.2 for N = 10 at the number of intervals M.

For future work, one may proceed with a higher order scheme for the problem under consideration for example by using the Crank-Nicolson method or those presented in [20, 21]. We are currently working in this direction. Also to note is that higher order parameter-uniform numerical methods for Burgers-Huxley equations are quasi-absent. The design of direct higher order parameter-uniform convergent methods is thus an interesting direction to explore.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

ED constructed the scheme and implemented the MATLAB code under JM's supervision. All authors participated in the analysis of the theoretical results and in the writing of the manuscript.


ED acknowledges the support received from the Department of Applied Mathematics of the Adama Science and Technology University which ensured that the author is registered for his Ph.D. study at the University for 4 years.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.


1. Kaushik A, Sharma MD. A uniformly convergent numerical method on non-uniform mesh for singularly perturbed unsteady Burger-Huxley equation. Appl Math Comput. (2008) 195:688–706. doi: 10.1016/j.amc.2007.05.067

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Liu LB, Liang Y, Zhang J, Bao X. A robust adaptive grid method for singularly perturbed Burger-Huxley equations. Elect Res Arch. (2020) 28:1439. doi: 10.3934/era.2020076

CrossRef Full Text | Google Scholar

3. Kadalbajoo MK, Sharma KK, Awasthi A. A parameter-uniform implicit difference scheme for solving time-dependent Burgers' equations. Appl Math Comput. (2005) 170:1365–93. doi: 10.1016/j.amc.2005.01.032

CrossRef Full Text | Google Scholar

4. Wang X, Zhu Z, Lu Y. Solitary wave solutions of the generalised Burgers-Huxley equation. J Phys A Math Gen. (1990) 23:271. doi: 10.1088/0305-4470/23/3/011

CrossRef Full Text | Google Scholar

5. Wazwaz AM. Travelling wave solutions of generalized forms of Burgers, Burgers-KdV and Burgers-Huxley equations. Appl Math Comput. (2005) 169:639–56. doi: 10.1016/j.amc.2004.09.081

CrossRef Full Text | Google Scholar

6. Yefimova OY, Kudryashov N. Exact solutions of the Burgers-Huxley equation. Appl Math Mech. (2004) 3:413–20. doi: 10.1016/S0021-8928(04)00055-3

CrossRef Full Text | Google Scholar

7. Appadu AR, Inan B, Tijani YO. Comparative study of some numerical methods for the Burgers-Huxley equation. Symmetry. (2019) 11:1333. doi: 10.3390/sym11111333

CrossRef Full Text | Google Scholar

8. Derzie EB, Munyakazi JB, Dinka TG. Parameter-uniform fitted operator method for singularly perturbed Burgers-Huxley equation. J Math Model. (2022) 10:515–34. doi: 10.22124/jmm.2022.21484.1883

CrossRef Full Text | Google Scholar

9. Gupta V, Kadalbajoo MK. A singular perturbation approach to solve Burgers-Huxley equation via monotone finite difference scheme on layer-adaptive mesh. Commun Nonlinear Sci Numer. (2011) 16:1825–44. doi: 10.1016/j.cnsns.2010.07.020

CrossRef Full Text | Google Scholar

10. Mickens RE. Nonstandard finite difference models of differential equations. World Scientific. (1994). doi: 10.1142/2081

CrossRef Full Text | Google Scholar

11. Anguelov R, Lubuma JMS. Contributions to the mathematics of the nonstandard finite difference method and applications. Numer Methods Partial Differ Equ. (2001) 17:518–43. doi: 10.1002/num.1025

CrossRef Full Text | Google Scholar

12. Kadalbajoo MK, Gupta V. Hybrid finite difference methods for solving modified Burgers and Burgers-Huxley equations. Neural Parallel Sci Comput. (2010) 18:409–22.

Google Scholar

13. Mickens RE. Advances in the applications of nonstandard finite difference schemes. World Scientific. (2005). doi: 10.1142/5884

CrossRef Full Text | Google Scholar

14. Mickens RE. Calculation of denominator functions for nonstandard finite difference schemes for differential equations satisfying a positivity condition. Numer Methods Partial Differ Equ. (2007) 23:672–91. doi: 10.1002/num.20198

CrossRef Full Text | Google Scholar

15. Kumari A, Kukreja V. Robust septic Hermite collocation technique for singularly perturbed generalized Hodgkin-Huxley equation. Int J Comput Math. (2022) 95:24. doi: 10.1080/00207160.2021.1939317

CrossRef Full Text | Google Scholar

16. Kumari A, Kukreja V. Error bounds for septic Hermite interpolation and its implementation to study modified Burgers' equation. Numer Algorithms. (2022) 89:22. doi: 10.1007/s11075-021-01173-y

CrossRef Full Text | Google Scholar

17. Munyakazi JB, Patidar KC. Limitations of Richardson's extrapolation for a high order fitted mesh method for self-adjoint singularly perturbed problems. J Appl Math Comput. (2010) 32:219–36. doi: 10.1007/s12190-009-0245-6

CrossRef Full Text | Google Scholar

18. Doolan EP, Miller JJ, Schilders WH. Uniform Numerical Methods for Problems with Initial and Boundary Layers. Dublin: Boole Press. (1980).

Google Scholar

19. Batiha B, Noorani M, Hashim I. Application of variational iteration method to the generalized Burgers-Huxley equation. Chaos, Solitons & Fractals. (2008) 36:660–3. doi: 10.1016/j.chaos.2006.06.080

CrossRef Full Text | Google Scholar

20. Kovacs E, A. class of new stable, explicit methods to solve the non-stationary heat equation. Numer Methods Partial Differential Eq. (2021) 37:2469–89. doi: 10.1002/num.22730

CrossRef Full Text | Google Scholar

21. Nagy A, Saleh M, Omle I, Humam K, Kovacs E. New stable, explicit, shifted-hopscotch algorithms for the heat equation. Math Comput Appl. (2021) 26:61. doi: 10.3390/mca26030061

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Burgers-Huxley equation, singularly perturbed problem, nonlinear equation, nonstandard finite difference methods, parameter-uniform convergence, error analysis, error bounds

Citation: Derzie EB, Munyakazi JB and Dinka TG (2023) A NSFD method for the singularly perturbed Burgers-Huxley equation. Front. Appl. Math. Stat. 9:1068890. doi: 10.3389/fams.2023.1068890

Received: 13 October 2022; Accepted: 27 January 2023;
Published: 13 February 2023.

Edited by:

Lei Zhang, Shanghai Jiao Tong University, China

Reviewed by:

Youssri Hassan Youssri, Cairo University, Egypt
V. K. Kukreja, Sant Longowal Institute of Engineering and Technology, India

Copyright © 2023 Derzie, Munyakazi and Dinka. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Eshetu B. Derzie, yes