Solution of Moving Boundary Space-Time Fractional Burger’s Equation

1 Department of Mathematics, Faculty of Science, Assiut University, New Valley Branch, El-Kharja 72511, Egypt 2 Department of Mathematics, Faculty of Science, Northern Border University, Arar 91431, Saudi Arabia 3 Department of Applied Mathematics, Faculty of Mathematical Sciences, University of Khartoum, 11111 Khartoum, Sudan 4Department of Mathematics, Faculty of Science, Sudan University of Science and Technology, 11115 Khartoum, Sudan 5Department of Mathematics, Faculty of Applied and Pure Science, International University of Africa, 14415 Khartoum, Sudan


Introduction
In recent years, fractional differential equations (FDEs) have become one of the most exciting and extremely active areas of research because of their potential applications in physics and engineering.These include electromagnetic, fluid flow, dynamical process in self-similar and porous structures, electrical networks, probability and statistics control theory of dynamical systems, chemical physics, optics, acoustic, viscoelasticity, electrochemistry, and material science [1][2][3][4][5][6][7][8][9].
There are different kinds of fractional integration and differentiation operators.The most famous one is the Riemann-Liouville definition [10][11][12], which has been used in various fields of science and engineering successfully, but this definition leads to the result that constant function differentiation is not zero.Caputo put definitions which give zero value for fractional differentiation of constant function, but these definitions require that the function should be smooth and differentiable [10][11][12].Recently, Jumarie derived definitions for the fractional integral and derivative called modified Riemann-Liouville [13][14][15], which are suitable for continuous and nondifferentiable functions and give differentiation of a constant function equal to zero.The modified Riemann-Liouville fractional definitions are used effectively in many different problems [16][17][18][19][20].
In the literature, there are many effective methods to treat FDEs such as the Adomian decomposition method, the variational iteration method, the homotopy perturbation method, the differential transform method, the finite difference method, the finite element method, the exponential function method [21], the fractional subequation method [22], the (  /)-expansion method [23,24], and the first integral method [25].Based on Jumarie's modified Riemann-Liouville derivative and the fractional Riccati equation    () =  + () 2 , Zhang and Zhang in [26] introduced the subequation method for solving nonlinear time fractional biological population model and (4 + 1)dimensional space-time fractional Fokas equation.Guo et al. in [27] improved the subequation method; they obtained the analytical solutions of the space-time fractional Whitham-Broer-Kaup and generalized Hirota-Satsuma coupled KdV equations by introducing a new general ansätz.Recently, by extending the fractional Riccati equation in [26] to the more general form    () =  + () 2 , Abdel-Salam and Yousif [28] presented the fractional Riccati expansion method to obtain exact solutions of the space-time fractional Kortewegde Vries equation, the spacetime fractional RLW equation, the space-time fractional Boussinesq equation, and the spacetime fractional Klein-Gordon equation.In addition, Li et al. in [29] extended fractional Riccati expansion method for solving the time fractional Burger's equation and the spacetime fractional Cahn-Hilliard equation.
The moving boundary problem is a nonlinear initialboundary value problem that requires extra boundary conditions to determine the motion of the boundary [30,31].Moving boundary problems arise in applications that involve time varying geometries.They have numerous applications in areas of physics, engineering, and biology.For example, coating flows, viscous sintering in industrial processes, melting and solidification problems, heat flow and diffusion with phase change, the formation of earth's crust according to plate tectonics theory, and the phenomenon of oxygen diffusion from blood into oxygen consuming tissue give rise to a moving boundary [32][33][34][35][36]. Generally, the solution differs according to the boundary conditions, coordinate dimensions, and the method of solution.Numerous numerical methods have been considered for solving this problem; examples include front-tracking, the front-fixing, the domain-fixing, finite-difference, finite-element, moving grid, phase field, and a fully implicit ghost-cell immersed boundary [37][38][39][40][41][42][43][44][45].Atkinson [46] studied time fractional diffusion for the motion of planar boundaries as well as cylindrical and spherical ones, the solution obtained for slow growth by the method of matched asymptotic expansion.Yin and Xu [47] used two-parameter regular perturbation technique and Fourier-Laplace transform method in finding analytical solution given in terms of the Wright function for a problem of a drug released from a polymeric matrix that can be dissolved into a solvent.Kushwaha and Kumar [48] approximated the solution of a model governed by spacetime fractional derivative for a moving boundary problem which occurs in fluviodeltaic sedimentation process on earth surface.Few researchers investigated some exact analytical solutions using the similarity transformation method and matched asymptotic expansion method.Also, the analytical and numerical solutions of the fractional single phase moving boundary problem were obtained by Li et al. [49].
Burger's equation is a classical nonlinear differential equation which was firstly introduced by Burger in 1948.It is used as a model for many nonlinear physical phenomena such as acoustics, continuous stochastic processes, dispersive water waves, gas dynamics, heat conduction, longitudinal elastic waves in an isotropic solid, number theory, shock waves, and turbulence [50][51][52][53][54][55].The space-time fractional Burger's equation, which is a transformed generalization of the Burger's equation, is defined as follows: where  = (, ),,  are arbitrary constants and  is the fractional order derivative.The initial-boundary conditions of (1) are where () ≤  < ∞,  ≥ 0, ℎ(), (), and () are given functions of their arguments; also, () is continuous and bounded.The function () is assumed to be smooth and it describes the motion of the boundary.It is important to note that when  = 1 and the motion of the boundary is unknown, one obtains a Stefan problem [56][57][58].The aim of the paper is to find the analytic solution of the moving boundary spacetime Burger's equation ( 1) via variable coefficient fractional Riccati expansion method.To solve ( 1)-( 2), a linear transformation is used to convert a constant coefficient FDE with moving boundary into a variable coefficient FDE with a fixed boundary.
The rest of this paper is organized as follows: the description of the fractional Riccati expansion method with variable coefficients is presented in Section 2. In Section 3, the solution of the moving boundary space-time Burger's equation is studied.In Section 4, discussion and conclusion are presented.

Description of the Fractional Riccati Expansion Method with Variable Coefficients
In this section, we present the fractional Riccati expansion method with variable coefficients to find exact analytical solutions of nonlinear FDEs.The fractional derivatives are described in sense of the modified Riemann-Liouville derivative defined by Jumarie [59][60][61] as which has merits over the original one; for example, the order derivative of a constant is zero.Some properties of the Jumarie's modified Riemann-Liouville derivative are The above properties play an important role in the fractional Riccati expansion method with variable coefficients.
Consider a nonlinear FDE in two variables  and  to be where  = (, ) is an unknown function,     and     are Jumarie's modified Riemann-Liouville derivatives of , and  is a polynomial in  and its various partial derivatives, in which the highest order derivatives and nonlinear terms are involved.
Suppose that () can be expressed by a finite power series of () as where  =  + (),  is arbitrary constant and (),   () ( = 0, 1, 2, . . ., ) are arbitrary functions of  to be determined later,  is a positive integer determined by balancing the linear term of the highest order with the nonlinear term in (5), and  = () satisfies the following fractional Riccati equation: where  and  are constants.Using the Mittag-Leffler function in one parameter   () = ∑ ∞ =0 (  /Γ(1 + )) ( > 0), we obtain the following solutions of (7).
where the generalized hyperbolic and trigonometric functions are defined as Determining the integer  and substituting ( 6) with ( 7) into ( 5) and collecting all terms with the same order of (), then setting each coefficient of () to zero yields a system of overdetermined nonlinear fractional equations for (),   (), and .Solving this system for (),   (), and , we get explicit expressions for   () and .Using the solutions of ( 7) into (6), we obtain exact solutions of FDE (5).
Remark 1.It can be easily found that if   () are constants and  is linear function of  and , then (6) becomes the same as (6) constructed in [28].Therefore, we may get more general exact solutions of the FDE (5).

Solution of Space-Time Fractional Burger's Equation
To solve the space-time fractional Burger's equation ( 1) with moving boundary (2), firstly, we apply the following transformation to (1): This transformation maps equation (1) into the following forced space-time Burger's equation: with initial-boundary conditions where   () = / and (0) = 0, which implies the compatibility condition (0) = (0, 0).The moving boundary problem ( 1)-( 2) of the fractional space-time Burger's equation is transformed into a fixed boundary problem with forced nonlinear fractional differential equation ( 16) in the quarter plane  ≥ 0,  ≥ 0 with boundary condition (17).In order to solve ( 16) by the fractional Riccati expansion method, we use the following transformation: In particular, when() =  then  represents a travelling wave with velocity .By applying ( 18) into ( 16), we obtain the following nonlinear fractional ordinary differential equation: where   = /.Balancing  2   with     gives  = 1.Therefore, the solution of ( 16) is expressed as Substituting ( 20) into ( 19) using ( 7) and setting the coefficients of () to zero, we get where  0 and  are arbitrary constants.Therefore, general formula of the solution of the forced space-time fractional Burger's equation ( 16) is From the solutions of (7), by choosing the special values of , , and the corresponding function (), we get the following solutions of ( 16): The remaining solutions can be obtained in a similar manner.When  = 1, we obtain the following forced Burger's equation (see [56][57][58]): Thus, the solutions ( 23) are reduced to the solutions of (24), which are given by Meanwhile, from solutions ( 25)-( 27), we see that the velocity of the wave is related to the moving boundary function ().To the best of our knowledge some of the obtained solutions appear for the first time concerning this problem.In order to understand the evolution of the solutions ( 25)-( 27), the main soliton features of them are investigated by using direct computer simulations with the accuracy 10 −9 .We discuss the effect of the moving boundary condition in some special cases, such as linear and periodic moving boundary to the kink soliton solution (26).
Firstly, we study the features of the solution (26) for linear moving boundary.In order to do that we choose Taking  =  = −1,  0 =  = 1, then, from ( 17) and ( 26), the initial-boundary condition is given by  () = 1 + tanh ((V − 2) ) , ℎ() = 1 + tanh () .(29) Note that in the static case when V → 0, solution (26) becomes kink solution with velocity 2 0 .Figure 1 shows the solution (26) of the forced Burger's equation (24) in the moving frame of reference for a linear moving boundary (28), with V = 100, −100, 0.001, −0.001.The kink behavior appears clearly when the absolute value of V is sufficiently small as shown in Figures 1(c) and 1(d).The velocity of the boundary in the moving frame moves in a concave shape as shown in Figure 1(a), whilst moves in a convex shape as shown in Figure 1(b) in the -direction.Also, the velocity of the boundary in the moving frame is shown as a convex shape in the -direction in Figures 1(c) and 1(d).
In the special case of a linear moving boundary, the solution of the Burger's equation can be found directly from (24) via a change of dependent variable.In fact, if () = V, then (24) becomes By taking the dependent variable to be then (30) becomes the Burger's equation with the initial-boundary conditions Secondly, we discuss the features of the kink soliton solution (26) for the periodic moving boundary.Now let us choose  () = 10 sin (V) . ( The initial-boundary conditions have the form Figure 2 presents the surface plot of solution (26) for the forced Burger's equation (24) with periodic moving boundary (34) for V = 1 and V = 10.Because of the periodic property of the moving boundary (), the kink-type soliton propagation of the solution (26) shows the periodic property.

Conclusions
In this paper, the fractional Riccati expansion method with variable coefficients has been successfully applied to the forced space-time fractional Burger's equation.A number of new analytical solutions have been obtained.Figures 1   and 2 show the solutions of the forced Burger's equation with linear and periodic moving boundaries at  = 1.
The variable coefficient fractional Riccati expansion method can be applied to other nonlinear FDEs with variable coefficients.