Elsevier

Chaos, Solitons & Fractals

Volume 40, Issue 4, 30 May 2009, Pages 1616-1622
Chaos, Solitons & Fractals

Solving the fractional BBM–Burgers equation using the homotopy analysis method

https://doi.org/10.1016/j.chaos.2007.09.042Get rights and content

Abstract

Based on the homotopy analysis method, a scheme is developed to obtain approximation solution of a fractional BBM–Burgers equation with initial condition, which is introduced by replacing some integer-order space derivatives by fractional derivatives. The fractional derivatives are described in the Caputo sense. So the traditional homotopy analysis method for differential equations of integer-order is directly extended to derive explicit and numerical solutions of the fractional differential equations. The solutions of our model equation are calculated in the form of convergent series with easily computable components.

Introduction

In the past decades, both mathematicians and physicists have devoted considerable effort to the study of explicit solutions to nonlinear integer-order differential equation. Many powerful methods have been presented [1], [2], [3], [4], [5], [6], [7], [8], [9], [10], [11], [12], [13], [14], [15], [16], [17], [18], [19], [20], [21], [22], [23], [24], [25], [26], [27], [28], [29], [30], [31], [32], [33], [34]. Among them, the homotopy analysis method (HAM) [22], [23], [24], [25], [26], [27] provides an effective procedure for explicit and numerical solutions of a wide and general class of differential systems representing real physical problems. Based on homotopy of topology, the validity of the HAM is independent of whether or not there exists small parameters in the considered equation. Therefore, the HAM can overcome the foregoing restrictions and limitations of perturbation techniques so that it provides us with a possibility to analyze strongly nonlinear problems. This method has been successfully applied to solve many types of nonlinear problems by others [28], [29], [30], [31], [32], [33], [34]. However, the application of HAM only circumscribes integer-order differential equation. Here, we introduce HAM to nonlinear fractional differential equation.

In recent years, considerable interest in fractional differential equation has been stimulated due to their numerous applications in the areas of physics and engineering [35]. Many important phenomena in electromagnetics, acoustics, viscoelasticity, electrochemistry and material science are well described by fractional differential equation [36], [37], [38]. The solution of fractional differential equation is much involved. In general, there exists no method that yields an exact solution for fractional differential equation. Only approximation solution can be derived using linearization or perturbation method.

The aim of this paper is to directly extend the HAM to consider the numerical solution of the fractional BBM–Burgers equation with space-fractional derivatives of the formut+Dxαu+uux-auxx-buxxt=0,t>0,0<α1,where a, b are constants and α is parameter describing the order of the fractional space-derivatives. The function u(x, t) is assumed to be a causal function of space, i.e. vanishing for t < 0 and x < 0. The fractional derivatives are considered in the Caputo sense [39]. The general response expression contains a parameter describing the order of the fractional derivative that can be varied to obtain various responses. In the case of α = 1, Eq. (1.1) reduces to the classical nonlinear BBM–Burgers equation. More important, above procedure is just an algebraic algorithm and can be easily applied in the symbolic computation system Maple.

Although there are a lot of studies for the classical BBM–Burgers equation and some profound results have been established, it seems that detailed studies of the nonlinear fractional differential equation are only beginning. According to our knowledge, the paper represents the first available numerical solution of the fractional BBM–Burgers equation with space-fractional derivatives.

The paper has been organized as follows. In Section 2, a brief review of the theory of fractional calculus will be given to fix notation and provide a convenient reference. In Section 3, we give analysis of the HAM. In Section 4, we extend the application of the HAM to construct numerical solution for the fractional BBM–Burgers equation. Conclusions are presented in Section 5.

Section snippets

Preliminaries and notations

In this section, let us recall essentials of fractional calculus first. The fractional calculus is a name for the theory of integrals and derivatives of arbitrary order, which unifies and generalizes the notions of integer-order differentiation and n-fold integration. There are many books [35], [36], [37], [38] that develop fractional calculus and various definitions of fractional integration and differentiation, such as Grünwald–Letnikov’s definition, Riemann–Liouville’s definition, Caputo’s

Analysis of the homotopy analysis method

In the paper, we apply the homotopy analysis method [22], [23], [24], [25], [26], [27] to the discussed problem. We extend Liao’ basic ideas to the fractional differential equation.

Let us consider the fractional differential equationFD(u(x,t))=0,where FD is a fractional differential operator, x and t denote independent variables, u(x, t) is an unknown function. For simplicity, we ignore all boundary or initial conditions, which can be treated in the same way.

Based on the constructed zero-order

Application

To demonstrate the effectiveness of the method, we consider Eq. (1.1) with the initial conditionu(x,0)=ρx2,ρC

We choose the linear integer-order operatorL[U(x,t;q)]=U(x,t;q)t.

Furthermore, Eq. (1.1) suggests to define the nonlinear fractional differential operatorNFD[U(x,t;q)]=Ut(x,t;q)+DxαU(x,t;q)+U(x,t;q)Ux(x,t;q)-aUxx(x,t;q)-bUxxt(x,t;q).

Using above definition, we construct the zeroth-order deformation equation(1-q)L[U(x,t;q)-u0(x,t)]=qhNFD[U(x,t;q)].

Obviously, when q = 0 and q = 1,U(x,t;0)=u0(x

Conclusions

Here, we solve fractional BBM–Burgers equation with space-fractional derivatives. We can also apply the HAM to solve fractional BBM–Burgers equation with time-fractional derivatives,Dtαu+ux+uux-auxx-buxxt=0,t>0,0<α1,

At the same time, we choose linear integer-order operator L[U(x,t;q)]=Ux(x,t;q), Uxx(x, t; q), or Uxxt(x, t; q). If we choose linear fractional operator, there are difficulties in solving process. How to conquer the difficulties by modifying the method need us study further.

The

Acknowledgement

The work is partially supported by the National Key Basic Research Project of China under the Grant No. 2004CB318000.

References (39)

  • S.K. Liu et al.

    Phys Lett A

    (2001)
  • J.H. He et al.

    Chaos Solitons & Fractals

    (2006)
  • E.G. Fan et al.

    Phys Lett A

    (2001)
  • Z.Y. Yan

    Phys Lett A

    (2001)
  • Z.S.

    Chaos Solitons & Fractals

    (2003)
  • B. Li et al.

    Chaos Solitons & Fractals

    (2003)
  • Y. Chen et al.

    Appl Math Comput

    (2004)
  • F.D. Xie et al.

    Chaos Solitons & Fractals

    (2005)
  • Q. Wang et al.

    Chaos Solitons & Fractals

    (2006)
  • L.N. Song et al.

    Appl Math Comput

    (2006)
  • S.J. Liao

    Int J Non-linear Mech

    (1999)
  • S.J. Liao

    Appl Math Comput

    (2004)
  • S.J. Liao

    Int J Heat Mass Transfer

    (2005)
  • S.J. Liao

    Appl Math Comput

    (2005)
  • S.J. Liao et al.

    Int J Heat Mass Transfer

    (2006)
  • S. Abbasbandy

    Phys Lett A

    (2006)
  • S. Abbasbandy

    Phys Lett A

    (2007)
  • M.J. Ablowitz et al.

    Solitons, nonlinear evolution equations and inverse scattering

    (1991)
  • M. Wadat

    J Phys Soc Jpn

    (1972)
  • Cited by (56)

    • A new modified definition of Caputo–Fabrizio fractional-order derivative and their applications to the Multi Step Homotopy Analysis Method (MHAM)

      2019, Journal of Computational and Applied Mathematics
      Citation Excerpt :

      This technique has been applied to solve linear and non-linear fractional partial differential equations [25,26]. For instance, the fractional KdV–Burgers–Kuramoto equation [27], the non-linear fractional Riccati equation [28], fractional initial value problems (fIVPs) [29], the Benjamin–Bona–Mahony–Burgers equation [30], the fractional Klein–Gordon equation [31,32], this method has also been used to solve linear homogeneous one and two-dimensional fractional heat-like partial differential equations subject to the Neumann boundary conditions [33], linear and non-linear homogeneous fractional diffusion-wave equations [34] or linear and non-linear systems of FDEs [35]. Other interesting applications of HAM can be found in [36–40].

    View all citing articles on Scopus
    View full text