Skip to main content

Theory and Modern Applications

Analytic and numerical solutions of nonlinear diffusion equations via symmetry reductions

Abstract

In this article, the authors study analytic and numerical solutions of nonlinear diffusion equations of Fisher’s type with the help of classical Lie symmetry method. Lie symmetries are used to reduce the equations into ordinary differential equations (ODEs). Lie group classification with respect to time dependent coefficient and optimal system of one-dimensional sub-algebras is obtained. Then sub-algebras are used to construct symmetry reduction and analytic solutions. Finally, numerical solutions of nonlinear diffusion equations are obtained by using one of the differential quadrature methods.

1 Introduction

Consider the nonlinear diffusion equation

u t =D 2 u x 2 +g(u).
(1)

When g(u)=αu(1u), (1) stands for the Fisher equation, put forward by Fisher [1] as a model for the spatial and temporal propagation of a viral gene in an infinite medium. This equation also expresses a one-dimensional reaction-diffusion model for the evolution of the infected population. The equation is defined by

u t = 2 u x 2 +αu(1u)
(2)

in which α>0 is a parameter. Kawahara and Tanaka [2] found an exact solution describing the coalescence of two traveling wave fronts of (2). Applications of traveling wave fronts appear in biology, chemistry, and medicine [3]. Such wave fronts were studied by Fisher for the first time in 1930s by considering (2).

The Fisher equation (2) appears in chemical kinetics [4], in logistic population growth models [5], autocatalytic chemical reactions, branching Brownian motion processes, flame propagation, and neurophysiology. The reaction-diffusion equation (2) also expresses a model equation for the evolution of a neutron population in a nuclear reactor [6] and also arises in the study of chemical wave propagation [7]. This equation includes the effects of linear diffusion via u x x and nonlinear local multiplication or reaction via u(1u).

There is a large cycle of works on mathematical properties and discussion of the Fisher equation in the literature. Larson [8], Kawahara and Tanaka [2], and Brazhnik and Tyson [5] provided excellent summaries of the Fisher equation. Ablowitz and Zepetella [9] gave the presentation of explicit solutions of the Fisher equation for a special wave speed. For the generalized the Fisher equation, Wang [10] presented the exact and explicit solitary wave solutions. Then Wazwaz and Gorguis [11] considered an analytic study of the Fisher equation by using the Adomian decomposition method. But it was not before 1974 when numerical solutions of the Fisher equation were available in the literature. Gazdag and Canosa [12] were the first to study numerical solutions of the Fisher equation with a pseudo-spectral approach. Afterwards, a lot of researchers have studied numerical solutions on the Fisher equation. Hagstrom and Keller [13] presented asymptotic boundary conditions by using a centered finite-difference algorithm. Afterwards, Evans and Sahimi [14] used an alternating group explicit iterative method to solve (2) and obtained satisfactory results, of a qualitatively similar nature. The numerical scheme considered in [14] is quite complicated and it causes unexpected high-frequency oscillations, which must be filtered out at each time step. Next, Parekh and Puri [15] and Twizell et al. [16] developed implicit and explicit finite differences algorithms for numerically solving the Fisher equation. Then Tang and Weber [17] proposed a Galerkin finite element method. Mickens [18] put forward a best finite-difference scheme for the Fisher equation. Afterwards, Garey and Shen [19] used a least-squares finite element method and Qiu and Sloan [20] used a moving mesh method for numerical solution of the Fisher equation. Rizwan [21] compared the nodal integral method and non-standard finite-difference schemes, and Khaled [22] proposed the Sinc collocation method. Daniel et al. [23] proposed a pseudo-spectral method for the numerical solution of the Fisher equation. Mittal and Kumar [24] studied the Fisher equation by applying wavelet Galerkin method, while Jiwari et al. [2530] studied numerical solutions of some nonlinear evolution equations by using differential quadrature method. Finally, for numerical solutions of the nonlinear Fisher reaction-diffusion equation, Mittal and Jain [31] proposed a numerical method, based on collocation of modified cubic B-splines over finite elements.

The purpose of this paper is to present analytic and numerical solutions via symmetry reductions of nonlinear diffusion equations of Fisher’s type defined as:

F1: Fisher’s equation

u t u x x αu(1u)=0;
(3)

F2: Fisher’s type equation

u t . u x x u 2 (1u)=0;
(4)

F3: Fisher’s type equation

u t . u x x u(ua)(1u)=0
(5)

in which a is a parameter.

Lie symmetries are used to reduce the equations to ordinary differential equations (ODEs). In the next section, lie group classification with respect to time dependent coefficient and optimal system of one-dimensional sub-algebras is obtained. Then, in Section 3, sub-algebras are used to construct symmetry reduction and analytic solutions. In Section 4, numerical solutions of nonlinear diffusion equations are obtained by using polynomial differential quadrature method. Finally, in Section 5, absolute, root mean square (RMS), and L errors are calculated.

2 Method of Lie symmetries

In this section, we recall the general procedure for determining symmetries for any system of partial differential equation. To begin, let us consider that the general case of a nonlinear system of partial differential equations of order n in p-independent and q-dependent variables is given as a system of equations,

Δ v ( x , u ( n ) ) =0,v=1,,l.
(6)

This entails x=( x 1 ,, x p ), u=( u 1 ,, u q ), and the derivatives of u with respect to x up to n, where u ( n ) expresses all the derivatives of u of all orders from 0 to n. We consider a one-parameter Lie group of infinitesimal transformation acting on the independent and dependent variables of the system (6)

( x ) i = x i + s ξ i ( x , u ) + O ( s 2 ) , i = 1 , , p , ( u ) j = u j + s η j ( x , u ) + O ( s 2 ) , j = 1 , , q ,
(7)

where s is the parameter of transformation and ξ i , η j are infinitesimals of transformation for the independent and dependent variables, respectively. The invariance of the system (6) under the infinitesimal transformation leads to the invariance conditions

p r ( n ) [ Δ v ( x , u ( n ) ) ] =0,v=1,,l whenever  Δ v ( x , u ( n ) ) =0,
(8)

where p r ( n ) is called the n th-order prolongation of infinitesimal generator given by

p r ( n ) =v+ α = 1 q ϕ α j ( x , u ( n ) ) u j α ,
(9)

where j=( j 1 ,, j k ), 1 j k p, 1kn, and the sum is over all orders of j. If j=k the coefficients ϕ j α of u j α will only depend on k th- and lower-order derivatives of u;

ϕ α j ( x , u ( n ) ) = D j ( ϕ α i = 1 p ξ i u i α ) + i = 1 p ξ i u j , i α ,
(10)

where u i α = u α x i and u j , i α = u j α x i .

One of the most important properties of these infinitesimal symmetries is that they form a Lie algebra under the usual Lie bracket.

3 Polynomial differential quadrature method

Differential quadrature method is a numerical technique to find the numerical solution of differential equations. The technique is used for the discretization of the spatial derivatives. The polynomial differential quadrature discretization of the first and the second derivatives at a point x i is given by the following equations:

u x ( x i ,t)= j = 1 N a i j u( x j ,t), u x x ( x i ,t)= j = 1 N b i j u( x j ,t),
(11)

where a i j and b i j express the weighting coefficients [32], i=1,2,,N. The following basis functions are used to obtain the weighting coefficients:

g k (x)= L ( x ) ( x x k ) L ( 1 ) ( x k ) ,k=1,2,,N,
(12)

where

L(x)=(x x 1 )(x x 2 )(x x N ),
(13)
L ( 1 ) ( x i )= k = 1 , k i N ( x i x k )
(14)

when x x i and then

g k (x)= x k ,k=1,2,,N
(15)

since both test functions given in (12) and (15) span the problem domain. The off-diagonal weighting coefficients for the first-order derivative are determined by using the set of basis functions given in (12) and the off-diagonal weighting coefficients of the first-order derivative are found as [32]

a i j = L ( 1 ) ( x i ) ( x i x j ) L ( 1 ) ( x j ) ,k=1,2,,N,ij.
(16)

The set of basis functions given in (15) are used to find the following diagonal weighting coefficients:

a i i = j = 1 , j i N a i j ,i=1,2,,N.
(17)

In the same way the weighting coefficients of the second-order derivative are derived by the same basis functions and the weighting coefficients are found as [32]

b i j =2 a i j ( a i i 1 x i x j ) ,i,j=1,2,,N,ij,
(18)
b i i = j = 1 , j i N b i j ,i=1,2,,N.
(19)

Similarly, Shu [32] proposed the weighting coefficients of higher-order derivatives in the explicit form

w i j ( r ) =r [ a i j w i i ( r 1 ) w i i ( r 1 ) x i x j ] ,for ij,i,j=1,2,,N;r=3,4,,N1,
(20)
w i i ( r ) = j = 1 , j i N w i j ( r ) ,for i=j,
(21)

where a i j and w i j ( r ) are the weighting coefficients of the first-order derivative and r th-order derivative, respectively.

4 Lie classical analysis for nonlinear diffusion Fisher’s type equations

In this section, we study the infinitesimal transformations and reductions by one-dimensional sub-algebras of (3)-(5) by applying the classical Lie symmetry method [33] one by one. Consider the one-parameter group of infinitesimal transformations in (x,t,u) given by

x = x + ε ξ ( x , t , u ) + O ( ε 2 ) , t = t + ε τ ( x , t , u ) + O ( ε 2 ) , u = u + ε η ( x , t , u ) + O ( ε 2 ) ,
(22)

where ε is the group parameter. The functions ξ, τ, η are the infinitesimals of the transformations for the variables x, t, and u, respectively. We shall denote the infinitesimals for u t , u x x by η x , η x x . The infinitesimals are as follows:

η x ( η ξ u x τ u t ) x + ξ u x x + τ u t x , η x x ( η ξ u x τ u t ) x x + ξ u x x x + τ u t x x .
(23)

Using these various extensions, the infinitesimal criterion for the invariance of (3)-(5) under the group (22) is given by

VH | H = 0 0,
(24)

where H=( x , t , v , v x , v t , v x x ) and v=θ( x , t ) are also solutions of (3)-(5). In (24), the prolongation of the tangent vector field V is given by

V=ξ x +τ t + η x u x + η t u t + η x x u x x +.
(25)

Now, substitute (23) and (25) into (24). Then we collect together the coefficients of u, u x , u t , u x x , u t t and set all of them to zero. Finally, we get a system of linear partial differential equations from which we can find ξ, τ, and η in practice.

Next we will use the above method to find the Lie symmetry group of the Fisher type equations (3)-(5). The prolongation of the tangent vector field of (3) is given by

V=ξ x +τ t + η x u x + η t u t + η x x u x x + η t t u t t .
(26)

After substituting (23) and (26) into (24), we have

η t η x x η+2uη=0.
(27)

Now, substituting (23) in (27), we get a system of linear partial differential equations. We set all the coefficients of u, u x , u t , u x x to zero. Then we obtain

ξ=a,τ=b,η=0,

where a and b are arbitrary constants. The associated vector fields for the one-parameter Lie group of infinitesimal transformations are V 1 and V 2 as given by

V 1 = x , V 2 = t .

Firstly, we construct an optimal system to classify the group-invariant solutions of (3) and the problem of finding an optimal subgroup is equivalent to that of finding an optimal system of sub-algebras. Here, by using the method presented in [33], we will construct an optimal system of one-dimensional sub-algebras of (22). The general one-parameter group of symmetries can be obtained by considering the linear combination V= a 1 V 1 + a 2 V 2 of given vector fields. But the explicit formulas for the above transformations are very complicated. Factually, they can be expressed uniquely in the form

g=exp( ε 1 V 1 )+exp( ε 2 V 2 ).
(28)

Our task is to simplify as many of the coefficients ε i , i=1,2 as possible though judicious applications of adjoint maps to V. By taking a 2 0, a 2 =1, we have V 2 + a 1 V 1 .

4.1 Reductions by one-dimensional sub-algebras

Case 1: Thus making use of this group transformation, the similarity variable and similarity solution are given by

ζ=xct,u(x,t)=F(ζ).
(29)

After substituting these similarity variables in (3), the equation reduced into ordinary differential equation (ODE) given by

c F (ζ) F (ζ)αF(ζ) ( 1 F ( ζ ) ) =0.
(30)

The exact solution corresponding to this ODE is given by

c= 5 α 6 ,F(ζ)= 1 4 1 2 tanh ( c 1 + 1 12 6 α ζ ) + 1 4 tanh ( c 1 + 1 12 6 α ζ ) 2 ,
(31)

where c 1 is a constant.

Thus, the exact solution of Fisher equation (3) is given by

u(x,t)= 1 4 1 2 tanh ( c 1 + 1 12 6 α ( x c t ) ) + 1 4 tanh ( c 1 + 1 12 6 α ( x c t ) ) 2 .
(32)

In this similar way, we can apply the Lie classical method to (4) and (5) and can study the exact solutions of these equations.

Case 2: Now, by applying the Lie classical method to (4), the following vector fields are obtained:

V 1 = x , V 2 = t .

The optimal system of these vector fields is

(i) V 2 + a 2 V 1 ,

where a 2 is an arbitrary constant.

(ii) V 1 .

The similarity variable and similarity solution corresponding to the basic vector field V 2 + a 2 V 1 are given by

ζ=xlt,u(x,t)=F(ζ),

where l= 1 a 2 , a 2 0. On using this similarity variable and similarity solution into (4), this equation reduces to

l F (ζ) F (ζ) F 2 (ζ) ( 1 F ( ζ ) ) =0.
(33)

The solution of the above ODE is given by

l= 1 2 ,F(ζ)= 1 2 1 2 tanh ( c 2 + 1 4 2 ζ ) .
(34)

Thus, the exact solution of the Fisher equation (4) is given by

u(x,t)= 1 2 1 2 tanh ( c 2 + 1 4 2 ( x 1 2 t ) ) ,
(35)

where c 2 is a constant.

Corresponding to the basic vector field V 1 in the optimal system, we can obtain only a constant solution.

Case 3: After utilizing the Lie classical method for (5), the following vector fields are obtained:

V 1 = x , V 2 = t .

The optimal systems of these vector fields contain the following sub-algebras:

(i) V 2 + a 3 V 1 ,

where a 3 is an arbitrary constant.

(ii) V 1 .

For the sub-algebra, V 2 + a 3 V 1 , the similarity variable and similarity solution are

ζ=xmt,u(x,t)=F(ζ),where m= 1 a 3 , a 3 0.

Using these similarity variables in the nonlinear diffusion equation (5), the reduced ODE is

m F (ζ) F (ζ)F(ζ) ( 1 F ( ζ ) ) ( F ( ζ ) a ) =0.
(36)

The exact solution of the above ODE is given as

m= 1 2 (1+a),F(ζ)= 1 2 + 1 2 a+ ( 1 2 1 2 a ) tanh ( c 3 + 1 4 2 ( 1 + a ) ζ ) ,
(37)

where c 3 is an arbitrary constant.

Now, the exact solution of the nonlinear diffusion equation (5) is given by

u(x,t)= 1 2 + 1 2 a+ ( 1 2 1 2 a ) tanh ( c 3 + 1 4 2 ( 1 + a ) ( x m t ) ) .
(38)

For the sub-algebra V 1 , only a constant can be obtained.

5 Comparative study of exact and numerical solutions of Fisher’s type equations

In this section, a comparative study of exact and numerical solutions is made by finding a numerical solution of the equations with the help of the polynomial differential quadrature method. For this purpose, the exact solutions obtained in Section 4 are used for initial and boundary conditions to find the numerical solution by polynomial differential quadrature method (PDQM). Discretizing the spatial derivatives of equations (3)-(5) by using PDQM at the point x i , we have

d u i d t = j = 1 N b i j u j +α u i (1 u i ),i=1,2,,N,
(39)
d u i . d t = j = 1 N b i j u j + u i 2 (1 u i ),i=1,2,,N,
(40)
d u i . d t = j = 1 N b i j u j + u i (1 u i )( u i a),i=1,2,,N,
(41)

where α, a are parameters, b i j are weighting coefficients of the second-order partial derivative and u i =u( x i ,t). Equations (39)-(40) are systems of first-order nonlinear differential equations. The initial and boundary conditions are taken from the analytic solutions obtained by Lie symmetry method in Section 4. Finally, the systems of initial and boundary value problems are solved by Pike and Roe’s fourth-stage RK4 [34].

5.1 Numerical experiments and discussion

In this subsection, three particular numerical examples are considered with fixed values of the arbitrary constants occurring in the solutions of (3)-(5) and numerical solutions are obtained by using the polynomial differential quadrature method. The whole computation work is done by the MATLAB and DEV C++ software. Absolute errors L , root mean square error (RMS), and L 2 are computed according to the following formulas:

L Error = max 1 i n | e i | , RMS Error = ( i = 1 n e i 2 n ) 1 / 2 , and L 2 Error = ( i = 1 n e i 2 ) 1 / 2 ,

where e i =( u i U i ), u i are approximated solutions and U i are exact solutions.

Example 1 Consider (3) over the domain [0,1] with the following initial conditions:

u(x,0)= 1 4 1 2 tanh ( c 1 + 1 12 6 α ( x ) ) + 1 4 tanh ( c 1 + 1 12 6 α ( x ) ) 2

and boundary conditions:

u ( 0 , t ) = 1 4 1 2 tanh ( c 1 + 1 12 6 α ( c t ) ) + 1 4 tanh ( c 1 + 1 12 6 α ( c t ) ) 2 , u ( 1 , t ) = 1 4 1 2 tanh ( c 1 + 1 12 6 α ( 1 c t ) ) + 1 4 tanh ( c 1 + 1 12 6 α ( 1 c t ) ) 2 ,

where c= 5 α 6 and c 1 is an arbitrary constant.

The exact solution of the equation is taken from (32). Numerical solutions in the form of errors of the example are given in Table 1. The table shows that the errors are small and negligible. Figure 1 compares the numerical and exact solutions in 3D form and it is concluded that the solutions are very similar.

Figure 1
figure 1

Numerical (left) and exact (right) solutions of Example 1 for c 1 =α=1 up to time t=1 .

Table 1 L , RMS, and L 2 errors of Example 1 at different times t and for different constants

Example 2 In this example, we have considered the Fisher equation (4) over the domain [0,1] with the following initial conditions:

u(x,t)= 1 2 1 2 tanh ( c 2 + 1 4 2 ( x ) )

and boundary conditions:

u ( 0 , t ) = 1 2 1 2 tanh ( c 2 + 1 4 2 ( 1 2 t ) ) , u ( 1 , t ) = 1 2 1 2 tanh ( c 2 + 1 4 2 ( 1 1 2 t ) ) .

The exact solution of the equation is taken from (35). Numerical solutions of the example are given in Table 2 and Figure 2. Table 2 presents the errors at different times and parameter values which are small and negligible. Figure 2 depicts a comparison of numerical and exact solutions in 3D form and it is concluded that the solutions are very similar.

Figure 2
figure 2

Numerical (left) and exact (right) solutions of Example 2 for c 2 =1 up to time t=1 .

Table 2 L , RMS, and L 2 errors of Example 2 at different times t and for different constants

Example 3 Consider the Fisher equation (5) over the domain [0,1] with the following initial and boundary conditions:

u ( x , 0 ) = 1 2 + 1 2 a + ( 1 2 1 2 a ) tanh ( c 3 + 1 4 2 ( 1 + a ) ( x ) ) , u ( 0 , t ) = 1 2 + 1 2 a + ( 1 2 1 2 a ) tanh ( c 3 + 1 4 2 ( 1 + a ) ( m t ) ) , u ( 1 , t ) = 1 2 + 1 2 a + ( 1 2 1 2 a ) tanh ( c 3 + 1 4 2 ( 1 + a ) ( 1 m t ) ) .

In this example, the exact solution is taken from (38). Table 3 and Figure 3 present the numerical results of the example. Table 3 presents the errors at different times and parameter values which are small and negligible. Figure 3 depicts a comparison of the numerical and exact solutions in 3D form and it is concluded that the solutions are very similar.

Figure 3
figure 3

Numerical (left) and exact (right) solutions of Example 3 for c 3 =1 , a=0.5 up to time t=1 .

Table 3 L , RMS, and L 2 errors of Example 3 at different times t and for different constants

6 Conclusion

In this article, the authors studied analytic and numerical solutions of the Fisher type equations with the help of the classical Lie symmetry method and the polynomial differential quadrature method. The Lie symmetry method is utilized to investigate the symmetries and invariant solutions of the equations. By determining the transformation group under which a given system is invariant, information about the invariants and symmetries of that equation is obtained. This information, in turn, is used to determine similarity variables that reduce the number of independent variables. The vector fields of the optimal system lead to a reduction of the nonlinear system of partial differential equations to ordinary differential equations. The infinitesimal generators in the optimal system are used for reductions and exact solutions. Finally, the polynomial differential quadrature method is used to find the numerical solutions of the Fisher type equations with the help of initial and boundary conditions taken from the analytic solutions obtained by the classical Lie symmetry method. It is concluded that the numerical solutions are in good agreement with the analytical solutions. L , RMS, and L 2 errors are calculated for each equation with particular values of arbitrary constants, which are small and negligible.

References

  1. Fisher RA: The wave of advance of advantageous genes. Annu. Eugen. 1937, 7: 355-369. 10.1111/j.1469-1809.1937.tb02153.x

    Article  Google Scholar 

  2. Kawahara T, Tanaka M: Interactions of traveling fronts: an exact solution of a nonlinear diffusion equation. Phys. Lett. A 1983, 97(8):311-314. 10.1016/0375-9601(83)90648-5

    Article  MathSciNet  Google Scholar 

  3. Civan F, Sliepcevich CM: Differential quadrature for multi-dimensional problems. J. Math. Anal. Appl. 1984, 101: 423-443. 10.1016/0022-247X(84)90111-2

    Article  MathSciNet  MATH  Google Scholar 

  4. Tang S, Weber RO: Numerical study of Fisher’s equation by a Petrov-Galerkin finite element method. J. Aust. Math. Soc. Ser. B, Appl. Math 1991, 33: 27-38. 10.1017/S0334270000008602

    Article  MathSciNet  MATH  Google Scholar 

  5. Brazhnik P, Tyson J: On traveling wave solutions of Fisher’s equation in two spatial dimensions. SIAM J. Appl. Math. 1999, 2: 371-391.

    MathSciNet  Google Scholar 

  6. Canosa J: On a nonlinear diffusion equation describing population growth. IBM J. Res. Dev. 1973, 17: 307-313.

    Article  MathSciNet  MATH  Google Scholar 

  7. Arnold R, Showalter K, Tyson JJ: Propagation of chemical reactions in space. J. Chem. Educ. 1987, 64: 740-742. 10.1021/ed064p740

    Article  Google Scholar 

  8. Larson DA: Transient bounds and time-asymptotic behavior of solutions to nonlinear equations of Fisher’s type. SIAM J. Appl. Math. 1978, 34: 93-103. 10.1137/0134008

    Article  MathSciNet  MATH  Google Scholar 

  9. Ablowitz M, Zepetella A: Explicit solution of Fisher’s equation for a special wave speed. Bull. Math. Biol. 1979, 41: 835-840.

    Article  MathSciNet  MATH  Google Scholar 

  10. Wang XY: Exact and explicit solitary wave solutions for the generalized Fisher’s equation. Phys. Lett. A 1988, 131(4/5):277-279.

    Article  MathSciNet  Google Scholar 

  11. Wazwaz AM, Gorguis A: An analytic study of Fisher’s equation by using Adomian decomposition method. Appl. Math. Comput. 2004, 154: 609-620. 10.1016/S0096-3003(03)00738-0

    Article  MathSciNet  MATH  Google Scholar 

  12. Gazdag J, Canosa J: Numerical solution of Fisher’s equation. J. Appl. Probab. 1974, 11: 445-457. 10.2307/3212689

    Article  MathSciNet  MATH  Google Scholar 

  13. Hagstrom T, Keller HB: The numerical calculation of travelling wave solutions of nonlinear parabolic equations. SIAM J. Sci. Comput. 1986, 7: 978-988.

    Article  MathSciNet  MATH  Google Scholar 

  14. Evans DJ, Sahimi MS: The alternating group explicit iterative method to solve parabolic and hyperbolic partial differential equations. 2. Annual Review of Numerical Fluid Mechanics and Heat Transfer 1989, 283-389.

    Google Scholar 

  15. Parekh N, Puri S: A new numerical scheme for the Fisher’s equation. J. Phys. A 1990, 23: 1085-1091. 10.1088/0305-4470/23/7/015

    Article  MathSciNet  Google Scholar 

  16. Twizell EH, Wang Y, Price WG: Chaos free numerical solutions of reaction-diffusion equations. Proc. R. Soc. Lond., Math. Phys. Sci. 1990, 430: 541-576.

    Article  MathSciNet  MATH  Google Scholar 

  17. Tang S, Weber RO: Numerical study of Fisher’s equation by a Petrov-Galerkin finite element method. J. Aust. Math. Soc. Ser. B, Appl. Math 1991, 33: 27-38. 10.1017/S0334270000008602

    Article  MathSciNet  MATH  Google Scholar 

  18. Mickens RE: A best finite-difference scheme for Fisher’s equation. Numer. Methods Partial Differ. Equ. 1994, 10: 581-585. 10.1002/num.1690100505

    Article  MathSciNet  MATH  Google Scholar 

  19. Garey GF, Shen Y: Least-squares finite element approximation of Fisher’s reaction-diffusion equation. Numer. Methods Partial Differ. Equ. 1995, 11: 175-186. 10.1002/num.1690110206

    Article  Google Scholar 

  20. Qiu Y, Sloan DM: Numerical solution of Fisher’s equation using a moving mesh method. J. Comput. Phys. 1998, 146: 726-746. 10.1006/jcph.1998.6081

    Article  MathSciNet  MATH  Google Scholar 

  21. Rizwan U: Comparison of the nodal integral method and non-standard finite-difference scheme for the Fisher’s equation. SIAM J. Sci. Comput. 2001, 22: 1926-1942. 10.1137/S1064827597325463

    Article  MATH  Google Scholar 

  22. Khaled KA: Numerical study of Fisher’s reaction-diffusion equation by the Sinc collocation method. J. Comput. Appl. Math. 2001, 137: 245-255. 10.1016/S0377-0427(01)00356-9

    Article  MathSciNet  Google Scholar 

  23. Daniel O, Bernie D, Shizgal A: A pseudospectral method of solution of Fisher’s equation. J. Comput. Appl. Math. 2006, 193: 219-242. 10.1016/j.cam.2005.06.028

    Article  MathSciNet  MATH  Google Scholar 

  24. Mittal RC, Kumar S: Numerical study of Fisher’s equation by wavelet Galerkin method. Int. J. Comput. Math. 2006, 83: 287-298. 10.1080/00207160600717758

    Article  MathSciNet  MATH  Google Scholar 

  25. Mittal RC, Jiwari R: Numerical study of Fisher’s equation by using differential quadrature method. Int. J. Inf. Syst. Sci. 2008, 5(1):143-160.

    MathSciNet  Google Scholar 

  26. Jiwari R: Haar wavelet quasilinearization approach for numerical simulation of Burgers’ equation. Comput. Phys. Commun. 2012, 183: 2413-2423. 10.1016/j.cpc.2012.06.009

    Article  MathSciNet  MATH  Google Scholar 

  27. Jiwari R, Pandit S, Mittal RC: Numerical simulation of two-dimensional sine-Gordon solutions by differential quadrature method. Comput. Phys. Commun. 2012, 183: 600-616. 10.1016/j.cpc.2011.12.004

    Article  MathSciNet  MATH  Google Scholar 

  28. Mittal RC, Jiwari R: Differential quadrature method for numerical solution of coupled viscous Burgers’ equations. Int. J. Comput. Methods Eng. Sci. Mech. 2012, 13: 1-5.

    Article  MathSciNet  Google Scholar 

  29. Jiwari R, Mittal RC, Sharma KK: A numerical scheme based on weighted average differential quadrature method for the numerical solution of Burgers’ equation. Appl. Math. Comput. 2013, 219: 6680-6691. 10.1016/j.amc.2012.12.035

    Article  MathSciNet  MATH  Google Scholar 

  30. Mittal RC, Jiwari R, Sharma KK: A numerical scheme based on differential quadrature method to solve time dependent Burgers’ equation. Eng. Comput. 2013, 30(1):117-131.

    Article  Google Scholar 

  31. Mittal RC, Jain RK: Numerical solutions of nonlinear Fisher’s reaction-diffusion equation with modified cubic B -spline collocation method. Math. Sci. 2013, 7(12):1-10.

    MathSciNet  Google Scholar 

  32. Shu C: Differential Quadrature and Its Application in Engineering. Springer, London; 2000.

    Book  MATH  Google Scholar 

  33. Olver PJ: Applications of Lie Groups to Differential Equations. Springer, New York; 1986.

    Book  MATH  Google Scholar 

  34. Pike J, Roe PL: Accelerated convergence of Jameson’s finite volume Euler scheme using Van Der Houwen integrators. Comput. Fluids 1985, 13: 223-236. 10.1016/0045-7930(85)90027-1

    Article  MathSciNet  MATH  Google Scholar 

Download references

Acknowledgements

The authors are very thankful to the reviewers for their valuable suggestions to improve the quality of the paper.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ram Jiwari.

Additional information

Competing interests

The authors declare that they have no competing interest.

Authors’ contributions

The authors contributed equally in the writing of this paper. All authors read and approved final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

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 https://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Verma, A., Jiwari, R. & Koksal, M.E. Analytic and numerical solutions of nonlinear diffusion equations via symmetry reductions. Adv Differ Equ 2014, 229 (2014). https://doi.org/10.1186/1687-1847-2014-229

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-1847-2014-229

Keywords