# Research on Time-Space Fractional Model for Gravity Waves in Baroclinic Atmosphere

Authors: Yanwei Ren; Huanhe Dong; Xinzhu Meng; Hongwei Yang
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

The research of gravity solitary waves movement is of great significance to the study of ocean and atmosphere. Baroclinic atmosphere is a complex atmosphere, and it is closer to the real atmosphere. Thus, the study of gravity waves in complex atmosphere motion is becoming increasingly essential. Deriving fractional partial differential equation models to describe various waves in the atmosphere and ocean can open up a new window for us to understand the fluid movement more deeply. Generally, the time fractional equations are obtained to reflect the nonlinear waves and few space-time fractional equations are involved. In this paper, using multiscale analysis and perturbation method, from the basic dynamic multivariable equations under the baroclinic atmosphere, the integer order mKdV equation is derived to describe the gravity solitary waves which occur in the baroclinic atmosphere. Next, employing the semi-inverse and variational method, we get a new model under the Riemann-Liouville derivative definition, i.e., space-time fractional mKdV (STFmKdV) equation. Furthermore, the symmetry analysis and the nonlinear self-adjointness of STFmKdV equation are carried out and the conservation laws are analyzed. Finally, adopting theexp(-Φ(ξ)) method, we obtain five different solutions of STFmKdV equation by considering the different cases of the parameters (η,σ). Particularly, we study the formation and evolution of gravity solitary waves by considering the fractional derivatives of nonlinear terms.


---

## Body

## 1. Introduction

With the intercross and penetration of different knowledge, Rossby solitary waves have been applied to many fields successfully, such as physical oceanography, atmospheric science, hydraulic engineering, and communication engineering. Particularly, Rossby solitary waves have important theoretical significance and research value in marine atmospheric science. They have largely determined the impact of the oceans on the atmosphere and other climate change. On the time scale, the energy of Rossby waves determines the ocean energy spectrum, which makes the energy spread from east to west to form and maintain a strong ocean boundary flow, such as Kuroshio, Gulfstream, and East Australian flow. Great achievements have been made in this regard.As we know, Rossby solitary waves in the westerly shear flow were first found by Long [1]. He found that the amplitude of the Rossby waves satisfied KdV equation by the β-plane approximation(1)ut+μuux+δuxxx=0. With the development of Rossby waves theory, Wadati [2] derived the modified KdV equation(2)ut+μu2ux+δuxxx=0. In view of the barotropic fluid and stratified fluid model, the KdV model and the mKdV model are also generated to describe the generation and evolution of Rossby solitary waves by Redekopp [3]. Apart from the KdV model and the mKdV model, for other initial disturbance, employing a different time and space stretching transform, Boussinesq equation was derived by Meng and Lv [4]. Afterwards, the Rossby parameters β along with the changes of latitude were discussed by Luo [5] and generalized β-plane approximation was obtained. In recent years, in the theoretical study of Rossby waves, many new wave equations were obtained to describe the generation and evolution of various types of fluctuations in the ocean [6, 7], such as NLS equation, ILW-Burgers equation, and ZK-Burgers equation. In the past, predecessors studied wave equations in the barotropic atmospheric environment by using the β-plane approximation. But we know that the basic dynamic equations for describing the baroclinic atmospheric movement are more in line with the actual situation and are very complicated. And the baroclinic problem in real atmosphere is inevitable [8]. In this paper, starting with the basic equations adopting the Bousinesq approximation [9] and under the baroclinic atmospheric environment, using multiscale analysis and turbulence method, we get a new model (mKdV) to describe the Rossby solitary waves. The advantages of basic equations are as follows:( 1 ) The equations are multivariate, and the physical meaning of each variable is clear;( 2 ) The baroclinic atmosphere problem is considered to help us understand the generation and evolution of isolated waves in the ocean.In recent years, the study of integer partial differential equations has yielded many brilliant achievements [10–15]. Simultaneously, it has been found that fractional order partial differential equations also play an important role in many fields [16–22]. The fractional differentiation calculus [23, 24] was first developed by Liouville primary. Liouville expands the function into an exponential form and defines the q-order derivatives of this expanded form term by term. Afterwards, Riemann proposed a different definition that can be implemented to a power series with a negative power term. Finally, Ross and Oldham [25, 26] unified the two definitions, so that the application of fractional differential was further developed. Subsequently, a version of the Euler-Lagrange equations for problems of calculus of variation with fractional derivatives was formulated by Riewe in 1990s [27, 28]. Recently, Agrawal [29, 30] studied the fractional Euler-Lagrange equation deeper and a series of new methods have been put forward in his research, which provide a new idea for us to study fractional partial differential equations [31, 32]. The fact has shown that the new fractional equation is more suitable than the integer order equation due to the precise description of the nonlinear phenomena[33, 34]. At the same time, in the field of oceanography, the fractional partial differential equation can better describe the generation and evolution of solitary waves, which is more favorable for us to study the theory of fluctuation.Similar to the study of integer order differential equation, the conservation laws of the fractional differential equation are an important branch. As we know, if the fractional differential equation is an Euler-Lagrange equation, then conservation laws can be found using Noethem’s theorem by variational Lie point symmetries of this equation [35–37]. Lie symmetry analysis was proposed by Sophus Lie. The main idea of this method is that the infinitesimal transformation keeps the solution set of the partial differential equation unchanged. The Lie symmetry analysis offers an efficient and powerful tool for the study of conservation laws of fractional partial differential equation [38–42]. For this reason, the researchers are very interested in studying the symmetry analysis of fractional differential equations. As far as we know, in the past, the symmetric method was only used for time fraction partial differential equations (TFPDE), but has not been used to analyze space and time fraction partial differential equations (STFPDE) [43, 44]. In this paper, the Lie symmetry analysis was used to study the conservation law of the STFmKdV equation [45, 46].By studying the work of predecessors, we can find that several methods have been used to solve nonlinear partial differential equations, such as the trial equation method [47, 48], Hirota bilinear method [49], binary nonlinearization method [50], Darboux transformation method [51], Jacobi iteration method [52], (G′/G)-expansion method [53–55], exp-function method [48], sub-equation method [56], and others [57, 58]. Therefore, it is an important task to find an accurate and effective method to solve the fractional differential equation.This paper is organized as follows: In Section2, using multiscale analysis and turbulence method, from the basic dynamic multivariable equations under the baroclinic atmosphere [59, 60], the integer order mKdV equation is derived. In Section 3, we use the semi-inverse method to derive the Lagrangian form of the mKdV equation [61, 62] firstly. Then the Lagrangian space and time operator of the mKdV equation have been transformed into the fractional domain of the left Riemann-Liouville fractional differential operator. Finally, applying the variational method, we derive the STFmKdV from this Euler-Lagrangian equation. In Section 4, we first study the symmetry analysis of the fractional equation to obtain the corresponding infinitesimal generator of the equation. Then we discuss the nonlinear self-adjointness of the STFmKdV equation and finally get the conservation vectors of the equation. In Section 5, based on the STFmKdV equation, employing the exp⁡(-Φ(ξ)) method, and considering the different cases of the parameters(η,σ), we obtain five different solutions of the equation.


## 2. Derivation of the mKdV Equation

Using the sum of disturbance pressure gradient force and buoyancy force expressing the vertical pressure gradient force and gravity force, and adopting the Bousinesq approximation [60], the dimensionless basic dynamic equations of atmospheric motion are as follows [59]:(3)∂u′∂t′+UfLu′∂u′∂x′+v′∂u′∂y′+w′∂u′∂z′=-1ρs∂p′∂x′+v′,∂v′∂t′+UfLu′∂v′∂x′+v′∂v′∂y′+w′∂v′∂z′=-1ρs∂p′∂y′-u′,∂w′∂t′+UfLu′∂w′∂x′+v′∂w′∂y′+w′∂w′∂z′=gLδθDfUθ0-1ρs∂p′∂z′+θ′,∂θ′∂t′+UfLu′∂θ′∂x′+v′∂θ′∂y′+σUDfLδθw′=0,∂ρsu′∂x′+∂ρsv′∂y′+∂ρsw′∂z′=0.where u′, v′ are the level of the air speed, w′ is the vertical velocity, p′ is the atmospheric pressure, θ is the temperature field, and f is the Coriolis parameter. θ0 is the potential of environmental flow field and ρs is the density of environmental flow field; they are both the height functions.Because the second term of the fourth formula in the left side is lesser, we get the following approximation:(4)δθ~σUDfL,UfL~o1. Let the parameter ε=f2/N2(ε≪1), N2=gσ/θ0. By varying, (3) transforms to(5)∂u′∂t′+u′∂u′∂x′+v′∂u′∂y′+w′∂u′∂z′=-1ρs∂p′∂x′+v′,∂v′∂t′+u′∂v′∂x′+v′∂v′∂y′+w′∂v′∂z′=-1ρs∂p′∂y′-u′,∂w′∂t′+u′∂w′∂x′+v′∂w′∂y′+w′∂w′∂z′=ε-1-1ρs∂p′∂x′+θ′,∂θ′∂t′+u′∂θ′∂x′+v′∂θ′∂y′+w′=0,∂ρsu′∂x′+∂ρsv′∂y′+∂ρsw′∂z′=0.We introduce the multiscale variables (omitting the sign at the top right corner of the variables)(6)t=ε3/2t′,x=ε1/2x′,y=y′,z=z′, so long time and space scales are defined as(7)∂∂t=ε3/2∂∂t′,∂∂x=ε1/2∂∂x′,∂∂y=∂∂y′,∂∂z=∂∂z′.Further, according to the small parameter ε, u′,v′,w′,p′,θ′ in (5) can be expended:(8)u′=-∫0yUς,z-c+ελdς+ε1/2u0+εu1+ε3/2u2+ε2u3+⋯,v′=εv0+ε3/2v1+ε2v2+ε5/2v3+⋯,w′=εw0+ε3/2w1+ε2w2+ε5/2w3+⋯,θ′=Θy,z+ε1/2θ0+εθ1+ε3/2θ2+ε2θ3+⋯,p′=Py,z+ε1/2p0+εp1+ε3/2p2+ε2p3+⋯,where U, P, Θ are the function of y, z. U is the speed of the basic flow, P is the air pressure, and Θ is the temperature field. Obviously, the zonal flow is in the following forms:(9)U=Uy=0,as0≤y<y0,costant,asy≥y0, and the boundary conditions of (3) are(10)p′=0,asy=0,p′→0,asy→∞.Substituting (7) and (8) into (3), we get each order form for ε as follows:(11)
O
ϵ
0
:
1
ρ
s
∂
p
∂
y
-
∫
0
y
U
-
c
d
ς
=
0
,
1
ρ
s
∂
p
∂
z
-
Θ
=
0
.
∑
​
l
2 Assuming 1/ρs2≪1, we have(12)∂∂z∫0yUdς=Θy. Next, we write the first-, second-, and third-order approximation for ε as the following form:(13)∫0yU-cdς∂ui∂x+U-c+1vi+Θywi-1ρs∂pi∂x=Aui,1ρs∂pi∂y+ui=Avi,1ρs∂pi∂z-θi=Awi,∫0yU-cdς∂θi∂x-Θyvi-wi=Aθi,∂ρsui∂x+∂ρsvi∂y+∂ρswi∂z=0,i=0,1,2,3⋯ where(14)Au0=Av0=Aw0=Aθ0=0. By eliminating u0,v0,w0,θ0 in (13), we can obtain the equation for p0(15)Ly,z∂p0∂x=0,where (16)Ω=Uy-1+Uz2,Ωy=∂Ω∂y,Ωz=∂Ω∂z,Ly,z=∂2∂y2-Uy-1∂2∂z2+2Uz∂2∂y∂z+Uzz-ΩyΩ-UzΩzΩ∂∂y+Uy-1ΩzΩ-UzΩyΩ∂∂z-1UUzz-ΩyΩ-UzΩzΩ.Clearly, (15) is a variable separable equation. Assume its solution is(17)p0=p0~y,zAt,x, and under a certain definite solution condition, we can get p0~. Further, all the solutions of (15) can be obtained:(18)u0=u0~y,zAt,x,v0=v0~y,zAxt,x,w0=w0~y,zAxt,x,θ0=θ0~y,zAt,x.Further, we know that(19)Au1=u0∂u0∂x+v0∂u0∂y+w0∂u0∂z,Av1=0,Aw1=0,Aθ1=u0∂θ0∂x+v0∂θ0∂y. Similarly, we cannot get the equation for A(x,t). So, we eliminate u1,v1,w1,θ1 in (13) and (19), and we can obtain the equation of p1(20)Ly,z∂p1∂x=l1y,zAu1+l2y,zAθ1, where(21)l1y,z=--∂∂y-Uzz-Uz∂∂z-1U-c-1△△y+Uz△z,l2y,z=-1U-c∂∂y+Uzz+Uz∂∂z-1△△y+Uz△zSubstituting (19) into (20) and observing both ends of (20), we can get the solutions for the variables as follows:(22)u1=u1~y,zA2t,x,v1=v1~y,zAt,xAxt,x,w1=w1~y,zAx,tAxt,x,θ1=θ1~y,zA2t,x,p1=p1~y,zA2t,x.Next, we have a further discussion,(23)Au2=∂u0∂t+α∂u0∂x+u0∂u1∂x+u1∂u0∂x+v0∂u1∂y+v1∂u0∂y+w0∂u1∂z+w1∂u0∂z,Av2=∫0yU-cdς∂v0∂x,Aw2=0,Aθ2=∂θ0∂t+u0∂θ1∂x+u1∂θ2∂x+v0∂θ1∂y+v1∂θ0∂y, Simplify (13) and (23) to the following form:(24)Ly,z∂p2∂x≡l3y,z∂Av2∂x+l1y,zAu2+l2y,zAθ2, where(25)l1y,z=--∂∂y-Uzz-Uz∂∂z-1U-c-1△△y+Uz△zl2y,z=-1U-c∂∂y+Uzz+Uz∂∂z-1△△y+Uz△zl3y,z=-1U-c-Uz∂∂y+Uy-1∂∂z+1△Uz△y+Uy-1△z. We know that the homogeneous part in (24) is the same as (15). Substituting (18) and (22) into (19) and (23), and according to (13)×∂p0/∂x-(20)×∂p2/∂x, we can get(26)Au2=u0~At+λAx+3u0~u1~+v0~u1~y+v1~u0~y+w0~u1~z+w1~u0~zA2Ax,Av2=∫0yU-cdςv0~Axx,Aw2=0,Aθ2=θ0~At+2u0~θ1~+u1~θ0~+v0~θ1~y+v1~θ0~yA2Ax,(27)ddy∂p2∂xdp0~dy-p0~ddy∂p2∂x-Ly,zp0~∂p2∂x-Ly,z∂p2∂xp0~=-p0~l1y,zu0~+l2y,zθ0~At+l1y,zλAx+l1y,z3u0~u1~+v0~u1~y+v1~u0~y+w0~u1~z+w1~u0~z+l2y,z2u0~θ1~+u1~θ0~+v0~θ1~y+v1~θ0~yA2Ax+l3y,zU-cv0~Axxx. Integrating (27) over the domain [0,y0) leads to(28)∫p0~y0,zddy∂p2∂x-∂p2∂xdp0~y0,zdydz+∫∫0y0Lp0~y0,z∂p2∂x-L∂p2∂xp0~y0,zdydz=∫∫0y0-p0~l1y,zu0~+l2y,zθ0~At+l1y,zu0~λAx+l1y,z3u0~u1~+v0~u1~y+v1~u0~y+w0~u1~z+w1~u0~z+l2y,z2u0~θ1~+u1~θ0~+v0~θ1~y+v1~θ0~yA2Ax+l3y,zU-cv0~Axxxdydz. We know that the two ends of (28) are identically zero, so we can write the simple form as the following form:(29)At+a0λAx+a1A2Ax+a2Axxx=0, where(30)a0=l1y,zu0~l1y,zu0~+l2y,zθ0~,a1=l1y,z3u0~u1~+v0~u1~y+v1~u0~y+w0~u1~z+w1~u0~z+l2y,z2u0~θ1~+u1~θ0~+v0~θ1~y+v1~θ0~yl1y,zu0~+l2y,zθ0~,a2=l3y,zU-cv0~l1y,zu0~+l2y,zθ0~Remark. According to the above study, we obtain a new model in baroclinic atmosphere. Based on the nonlinear term A2Ax, the new model is called the generalized mKdV equation. Compared to the model which is obtained in the barotropic atmosphere, the new mKdV equation is more likely to describe the movement of solitary waves.


## 3. Formulation of the STFmKdV Equation

First we introduce the definition of Riemann-Liouville fractional derivatives and Caputo fractional derivatives[41, 63].Definition 1 (Riemann-Liouville fractional derivative [41, 63]).
f ( t ) is a function defined in the [a,b], for any nonnegative real α, satisfying n-1≤α<n,(31)aRDtαft≔DnaRDtα-nft=1Γn-αdndtn∫atfτt-τα+1-ndτ,∀t∈a,b, and this formula is called the α-order left Riemann-Liouville fractional derivative. And(32)tRDbαft≔-1nDntRDbα-nft=-1nΓn-αdndtn∫tbfττ-tα+1-ndτ,∀t∈a,b, and this formula is the α-order right Riemann-Liouville fractional derivative.Remark 2.
We note that if the functionf(t) is n-order continuous derivable on the interval of [a,b], when α tends to n, the left fractional derivative is the traditional n-order derivative. In addition, if the function f(t) is n-order continuous derivable on the interval of [a,b], when α tends to n, the right fractional derivative is the traditional n(n-1)-order derivative multiplied by (-1)n.Definition 3 (Caputo fractional derivative [41, 63]).
f ( t ) is a function defined in the [a,b], for any nonnegative real α, satisfying n-1≤α<n,(33)aCDtαft≔DaRtα-nDnft=1Γn-α∫atfnτt-τα+1-ndτ,∀t∈a,b, and this formula is the α-order left Caputo fractional derivative.(34)tCDbαft≔-1ntRDbα-nDnft=-1nΓn-α∫tbfnττ-tα+1-ndτ,∀t∈a,b, and this formula is the α-order right Caputo fractional derivative.Lemma 4 (see [61, 62]).
Riemann-Liouville fractional derivative and Caputo fractional derivative have the following relationship:(35)aRDtαft=DaCtαft+∑j=0n-1fiat-ai-αΓ1+i-α.In this section, based on the generalized mKdV equation obtained in Section2, we use semi-inverse method and variational method[63] to construct the generalized STFmKdV equation under the Riemann-Liouville derivative definition. Let a0,3=1, so that (29) transforms to the following form:(36)At+λAx+a1A2Ax+a2Axxx=0, where a1,a2 are arbitrary constants, A(x,t) denotes the amplitude of the Rossby waves, x∈R is the space variable in the propagation of the field, and t∈T(=[0,T0]) is the time variable. The main steps are arranged as follows[62].First of all, we introduceu(x,t) as a potential function. Let A(x,t)=ux(x,t), so that (36) can be written as(37)uxt+λuxx+a1ux2uxx+a2uxxx=0. The Lagrangian form of (36) can be defined using the semi-inverse method. μA is considered as a fixed function. The functional of the potential equation (37) can be represented by(38)Ju=∫Rdx∫Tdtux,tc1uxt+c2λuxx+c3a1ux2uxx+c4a2uxxxx, where c1,c2,c3, and c4 are Lagrangian multipliers to be determined later. Integrating (38) by parts and taking utR=uxR=uxT=uxxR=uxxxR=0 lead to(39)Ju=∫Rdx∫Tdt-c1uxut-c2λux2-13c3a1ux4+c4a2uxx2.Secondly, by applying the variation of this functional with respect tou(x,t), and integrating each term by parts, optimizing the variation δJ(u)=0, we have(40)2c1uxt+2c2λuxx+4c3a1ux2uxx+2c4a2uxxxx=0. Equation (40) is equivalent to (37), so we can get(41)c1=12,c2=12,c3=14,c4=12. A functional relation (39) is given to produce the direct Lagrange form of the mKdV equation.(42)Lut,ux,uxx,⋯=-12uxut-12λux2-112a1ux4+12a2uxx2.Thirdly, according to the Lagrangian form of the integer order, we can obtain the Lagrangian form of the space-time fractional order similarly.(43)FDtα∗Dxβu,Dxβu,Dx2βu,⋯=-12Dtαu∗Dxβu-12λDxβu2-112a1Dxβu4+12a2Dx2βu2, where 0≤α,β<1. The fractional derivative Dtαu(x,t) or Dxβu(x,t) in terms of the left Riemann-Liouville fractional derivative is defined by(44)Dςγfς=1Γk-γdkdςk∫dsς-sk-γ-1fs,k-1≤γ≤k,ς=ςt,x.Then, the STFmKdV-Burgers equation takes the following form:(45)JFu=∫Rdxβ∫TdtαF∗Dtα∗Dxβu,Dxβu,Dx2βu. The variation of (45) with respect to u(x,t) leads to(46)δJFu=∫Rdxβ∫Tdtα∂F∗∂uδu+∂F∗∂DtαuδDtαu+∂F∗∂DxβuδDxβu+∂F∗∂Dx2βuδDx2βu. Adopting the fractional integration rule and the right Riemann-Liouville fractional derivative, δJF∗(u) is written as(47)δJFu=∫Rdxβ∫Tdtα∂F∂u-Dtα∂F∂Dtαu-Dxβ∂F∂Dxβu+Dx2β∂F∂Dx2βuδu.Optimizing the variation of the functional, i.e.,δJF∗(u)=0, the Euler-Lagrange form for the STFmKdV equation leads to(48)∂F∂u-Dtα∂F∂Dtαu-Dxβ∂F∂Dxβu+Dx2β∂F∂Dx2βu=0. Substituting (43) into (48), we have(49)DtαDxβu+λDxβDxβu+a1Dxβu2Dx2βu+a2Dx2βDx2βu=0. Substituting Dxβu(x,t)=A(x,t) into (49), we have the STFmKdV equation(50)DtαA+λDxβA+a1A2DxβA+a2Dx3βA=0. In this paper, in order to make the content more complete, then we will study the conservation laws and the solutions of the STFmKdV equation.


## 4. Lie Symmetry Analysis and Conservation Laws of the STFmKdV Equation

In this section, we employ Lie symmetry analysis to discuss the conservation laws[64, 65] of STFmKdV equaiton which does not contain dissipation item. The details are as follows.


### 4.1. Lie Symmetry Analysis

The STFPDE of a functionf(x,t) with two independent variables considered in the Riemann-Liouville sense is defined as the following form:(51)Dtαu=∂αu∂tα=1Γn-α∂n∂tn∫0tt-sn-α-1fx,sds,n-1<α<n∈N,∂nfx,t∂tn,α=n∈N.∑​l First, we define a Lie group of point transformations under one parameter(52)x∗=x+εξx,t,A+oε2,t∗=t+ετx,t,A+oε2,u∗=u+εηx,t,A+oε2, where ξ,τ,η are infinitesimal functions and ε is a small continuous parameter. The associated Lie algebra is spanned by(53)X=ξx,t,A∂β∂xβ+τx,t,A∂α∂tα+ηx,t,A∂∂A. The prolonged generator can be defined(54)Prα,β,2X=ξx,t,A∂β∂xβ+τx,t,A∂α∂tα+ηx,t,A∂∂A+ηα,t∂∂DtαA+ηβ,x∂∂DxβA+ηβ,xx∂∂Dx2βA. On the basis of the infinitesimal invariance criterion, one can get(55)PrnXFF=0=0, where(56)F=DtαA+λDxβA+a1A2DxβA+a2Dx3βA. The operators ηα,t,ηβ,x,ηβ,xx are fractional extended symmetry operators defined as follows[64]:(57)ηα,t=Dtαη+ξDtαAx-DtαξAx+DtαADtτ-Dtα+1τA+τDtα+1A,ηβ,x=Dxβη+DxβADxξ-Dxβ+1ξA+ξDxβ+1A+τDxβAt-DxβτAt,ηβ,xx=Dxβηβ,x-AxxDxβξ-AxtDxβτ, where the symbols Dt,Dx represent the total derivative operators defined by(58)Dt=∂t+At∂A+Att∂At+Axt∂Ax+⋯,Dx=∂x+Ax∂A+Axx∂Ax+Axt∂At+⋯.Second, the conserved vectors of the STFmKdV equation are investigated as follows. Applying the second prolongationPr(α,β,2)X, we get(59)ηα,t+ληβ,x+2a1AηDxβA+a1A2ηβ,x+a2ηβ,xxx=0 Substituting (57) and (58) into (59), and equating the coefficients of alike partial derivatives of u, we can obtain the determining equations(60)ξt=ξA=τx=τA=ηAA=0,ατt-3βξx=0,αn∂tαηA-αn+1Dtn+1τ=0,βn∂xβηA-βn+1Dtn+1ξ=0,∂tαη-A∂tαηA+λ∂tαη-A∂tαηA+a0A2ηx+a1ηxxx=0.Solving these equations, the infinitesimals can be derived in the following form:(61)ξ=c1βx+c2,τ=3c1αt+c3,η=3c1β-1βA. where c1,c2,c3 are arbitrary constants. Thus, the corresponding infinitesimal generator can be written(62)X=xβ∂∂x+3tα∂∂t+3Aβ-1β∂∂A


### 4.2. Nonlinear Self-Adjointness

The concept of nonlinear self adjoint is proposed in the application of new conservation theorem to the conservation laws of equations. This concept is extended to the space-time fractional partial differential equations. The lagrangian form of nondissipative STFmKdV equation is given by(63)L=vx,tDtαA+λDxβA+a1A2DxβA+a2Dx3βA, where v(x,t) is a new dependent variable. The adjoint equation of the STFmKdV equation is defined by(64)F∗≡δLδA=0, where δ/δA is the Euler-Lagrange operator. According to (63) and (64), the adjoint equation of STFmKdV equation can be obtained as follows:(65)F∗=Dtα∗v+a1A2-λDxβ∗v-a2Dx3β∗v=0. Here, (Dtα)∗,(Dxβ)∗ are the adjoint operators of Dtα,Dxβ. For the Riemann-Liouville fractional differential operators, the corresponding adjoint operators[41] have the following form:(66)​0Dtα∗=-1tnITn-αDtn=DtCTα,​0Dxβ∗=-1xmIRm-βDxm=xCDRβ, where ITn-α,IRm-β are the right fractional integral operators of order n-α,m-β, defined by(67)ITn-αfx,t=1Γn-α∫tTfx,ττ-tα+1-ndτ,n=α+1,IRm-βfx,t=1Γm-β∫xRfς,tς-xβ+1-mdς,m=β+1. For nonlinear self-adjointness, let us assume v=ψ(x,t,A), where ψ(x,t,A)≠0.Substitutingv=ψ(x,t,A) into (65)(68)Dtα∗ψ+a1A2-λDxβ∗ψ-a2Dx3β∗ψ+3Dx2β∗ψAAx+3Dxβ∗ψAAAx2+3Dxβ∗ψAxAxx+3ψAAAxAxx+ψAAxxx+ψAAAAx3≡F∗. By using the method of undetermined coefficients, we obtain the following cases[64, 65]:(69)case  1:0<β<1,vx,t=C0forDtα=0Dtα,case2:0<β<1,vx,t=C1x+C2forDtα=0cDtα. Here, C0,C1,C2 are arbitrary constants. These solutions of v(x,t) are introduced in the Lagrangian form for the construction of conserved vectors in the next subsection.


### 4.3. Conservation Laws

It is well known thatC=(Ct,Cx) is called the conserved vector that satisfies the conservation equation(70)DtCt+DxCx=0,Noether’s theorem is widely adopted to construct the conserved vectors. Using this method, the conserved vectors are obtained by using the Noether operators for the Lagrangian form. According to[40], the fractional Noether operator for the variable x,t has been given.Before constructing the conserved vectors, we should consider the Lie characteristic function for the Lie symmetry infinitesimal generatorX=ξ(∂β/∂xβ)+τ(∂α/∂tα)+η(∂/∂A). The Lie characteristic function is defined by(71)W=η-ξAx-τAt=3Aβ-1β-xβAx-3tαAt,Then, the fractional Noether operator fort of the STFmKdV-Burgers equation is(72)Ct=∑k=0n-1-1kDtα-1-kWDtk∂F∗∂DtαA, Similarly, the fractional Noether operator for x of the STFmKdV-Burgers equation is(73)Cx=∑k=0m-1-1kDxβ-1-kWDxk∂F∗∂DxβA, Substituting the Lagrangian form (65) and (71) into (72) with case 1 in (69), let C0=1; the t-component of the conserved vectors can be obtained:(74)Ct=3β-1βIt1-αA-xβIt1-αAx-3αIt1-αtAt,as0<α<1,(75)Ct=3β-1βDtα-1A-xβDtα-1Ax-3αDtα-1tAt,as1<α<2, The x-component of the conserved vectors is(76)Cx=3β-1βDxβ-1A-1βDxβ-1xAx-3tαDxβ-1At,as1<α<2,In case 2 in (69), let C1=1,C2=0; the t-component of the conserved vectors can be obtained as follows:(77)Ct=x×3β-1βIt1-αA-xβIt1-αAx-3αIt1-αtAt,as0<α<1,(78)Ct=x×3β-1βDtα-1A-xβDtα-1Ax-3αDtα-1tAt,as1<α<2,The x-component of the conserved vectors is(79)Cx=x×3β-1βDxβ-1A-1βDxβ-1xAx-3tαDxβ-1At-3β-1βIx2-βA-1βIx2-βxAx-3tαIx2-βAt.


## 4.1. Lie Symmetry Analysis

The STFPDE of a functionf(x,t) with two independent variables considered in the Riemann-Liouville sense is defined as the following form:(51)Dtαu=∂αu∂tα=1Γn-α∂n∂tn∫0tt-sn-α-1fx,sds,n-1<α<n∈N,∂nfx,t∂tn,α=n∈N.∑​l First, we define a Lie group of point transformations under one parameter(52)x∗=x+εξx,t,A+oε2,t∗=t+ετx,t,A+oε2,u∗=u+εηx,t,A+oε2, where ξ,τ,η are infinitesimal functions and ε is a small continuous parameter. The associated Lie algebra is spanned by(53)X=ξx,t,A∂β∂xβ+τx,t,A∂α∂tα+ηx,t,A∂∂A. The prolonged generator can be defined(54)Prα,β,2X=ξx,t,A∂β∂xβ+τx,t,A∂α∂tα+ηx,t,A∂∂A+ηα,t∂∂DtαA+ηβ,x∂∂DxβA+ηβ,xx∂∂Dx2βA. On the basis of the infinitesimal invariance criterion, one can get(55)PrnXFF=0=0, where(56)F=DtαA+λDxβA+a1A2DxβA+a2Dx3βA. The operators ηα,t,ηβ,x,ηβ,xx are fractional extended symmetry operators defined as follows[64]:(57)ηα,t=Dtαη+ξDtαAx-DtαξAx+DtαADtτ-Dtα+1τA+τDtα+1A,ηβ,x=Dxβη+DxβADxξ-Dxβ+1ξA+ξDxβ+1A+τDxβAt-DxβτAt,ηβ,xx=Dxβηβ,x-AxxDxβξ-AxtDxβτ, where the symbols Dt,Dx represent the total derivative operators defined by(58)Dt=∂t+At∂A+Att∂At+Axt∂Ax+⋯,Dx=∂x+Ax∂A+Axx∂Ax+Axt∂At+⋯.Second, the conserved vectors of the STFmKdV equation are investigated as follows. Applying the second prolongationPr(α,β,2)X, we get(59)ηα,t+ληβ,x+2a1AηDxβA+a1A2ηβ,x+a2ηβ,xxx=0 Substituting (57) and (58) into (59), and equating the coefficients of alike partial derivatives of u, we can obtain the determining equations(60)ξt=ξA=τx=τA=ηAA=0,ατt-3βξx=0,αn∂tαηA-αn+1Dtn+1τ=0,βn∂xβηA-βn+1Dtn+1ξ=0,∂tαη-A∂tαηA+λ∂tαη-A∂tαηA+a0A2ηx+a1ηxxx=0.Solving these equations, the infinitesimals can be derived in the following form:(61)ξ=c1βx+c2,τ=3c1αt+c3,η=3c1β-1βA. where c1,c2,c3 are arbitrary constants. Thus, the corresponding infinitesimal generator can be written(62)X=xβ∂∂x+3tα∂∂t+3Aβ-1β∂∂A


## 4.2. Nonlinear Self-Adjointness

The concept of nonlinear self adjoint is proposed in the application of new conservation theorem to the conservation laws of equations. This concept is extended to the space-time fractional partial differential equations. The lagrangian form of nondissipative STFmKdV equation is given by(63)L=vx,tDtαA+λDxβA+a1A2DxβA+a2Dx3βA, where v(x,t) is a new dependent variable. The adjoint equation of the STFmKdV equation is defined by(64)F∗≡δLδA=0, where δ/δA is the Euler-Lagrange operator. According to (63) and (64), the adjoint equation of STFmKdV equation can be obtained as follows:(65)F∗=Dtα∗v+a1A2-λDxβ∗v-a2Dx3β∗v=0. Here, (Dtα)∗,(Dxβ)∗ are the adjoint operators of Dtα,Dxβ. For the Riemann-Liouville fractional differential operators, the corresponding adjoint operators[41] have the following form:(66)​0Dtα∗=-1tnITn-αDtn=DtCTα,​0Dxβ∗=-1xmIRm-βDxm=xCDRβ, where ITn-α,IRm-β are the right fractional integral operators of order n-α,m-β, defined by(67)ITn-αfx,t=1Γn-α∫tTfx,ττ-tα+1-ndτ,n=α+1,IRm-βfx,t=1Γm-β∫xRfς,tς-xβ+1-mdς,m=β+1. For nonlinear self-adjointness, let us assume v=ψ(x,t,A), where ψ(x,t,A)≠0.Substitutingv=ψ(x,t,A) into (65)(68)Dtα∗ψ+a1A2-λDxβ∗ψ-a2Dx3β∗ψ+3Dx2β∗ψAAx+3Dxβ∗ψAAAx2+3Dxβ∗ψAxAxx+3ψAAAxAxx+ψAAxxx+ψAAAAx3≡F∗. By using the method of undetermined coefficients, we obtain the following cases[64, 65]:(69)case  1:0<β<1,vx,t=C0forDtα=0Dtα,case2:0<β<1,vx,t=C1x+C2forDtα=0cDtα. Here, C0,C1,C2 are arbitrary constants. These solutions of v(x,t) are introduced in the Lagrangian form for the construction of conserved vectors in the next subsection.


## 4.3. Conservation Laws

It is well known thatC=(Ct,Cx) is called the conserved vector that satisfies the conservation equation(70)DtCt+DxCx=0,Noether’s theorem is widely adopted to construct the conserved vectors. Using this method, the conserved vectors are obtained by using the Noether operators for the Lagrangian form. According to[40], the fractional Noether operator for the variable x,t has been given.Before constructing the conserved vectors, we should consider the Lie characteristic function for the Lie symmetry infinitesimal generatorX=ξ(∂β/∂xβ)+τ(∂α/∂tα)+η(∂/∂A). The Lie characteristic function is defined by(71)W=η-ξAx-τAt=3Aβ-1β-xβAx-3tαAt,Then, the fractional Noether operator fort of the STFmKdV-Burgers equation is(72)Ct=∑k=0n-1-1kDtα-1-kWDtk∂F∗∂DtαA, Similarly, the fractional Noether operator for x of the STFmKdV-Burgers equation is(73)Cx=∑k=0m-1-1kDxβ-1-kWDxk∂F∗∂DxβA, Substituting the Lagrangian form (65) and (71) into (72) with case 1 in (69), let C0=1; the t-component of the conserved vectors can be obtained:(74)Ct=3β-1βIt1-αA-xβIt1-αAx-3αIt1-αtAt,as0<α<1,(75)Ct=3β-1βDtα-1A-xβDtα-1Ax-3αDtα-1tAt,as1<α<2, The x-component of the conserved vectors is(76)Cx=3β-1βDxβ-1A-1βDxβ-1xAx-3tαDxβ-1At,as1<α<2,In case 2 in (69), let C1=1,C2=0; the t-component of the conserved vectors can be obtained as follows:(77)Ct=x×3β-1βIt1-αA-xβIt1-αAx-3αIt1-αtAt,as0<α<1,(78)Ct=x×3β-1βDtα-1A-xβDtα-1Ax-3αDtα-1tAt,as1<α<2,The x-component of the conserved vectors is(79)Cx=x×3β-1βDxβ-1A-1βDxβ-1xAx-3tαDxβ-1At-3β-1βIx2-βA-1βIx2-βxAx-3tαIx2-βAt.


## 5. The Exact Solutions of the STFmKdV Equation

In this part, we deal with the exact solutions [66–68] of (50). Using exp(-Φ(ξ)) method, the main steps of this method to solve the space time fractional partial equation can be given as follows.Firstly, we introduce the fractional complex transform:(80)Ax.t=Uξ,ξ=xβΓ1+β-νtαΓ1+α, where ν is the wave speed. Substitute (80) into the STFmKdV equation as follows:(81)DtαA+λDxβA+a1A2DxβA+a2Dx3βA=0 Equation (81) transforms to an ordinary differential equation(82)-νU′+λU′+a1U2U′+a2U′′′=0. Integrating (82) with respect to ξ and setting the integration constant to zero, we have the following equation:(83)-ν+λU+a13U3+a2U′′=0.Secondly, balancing the highest order derivative term and the highest order nonlinear term in (83), we get the balancing number n=1. Thus the solution of (83) takes the form(84)Uξ=k0+k1e-Φξ, where k0,k1 are constants to be determined later. And Φ(ξ) satisfies the following auxiliary ordinary differential equation:(85)Φ′ξ=exp⁡-Φξ+ηexp⁡Φξ+σ.Substitute (84) and (81) into (83) and collect all terms with the same degree of exp(-Φ(ξ))n together. Equating each coefficient of the same degree of exp(-Φ(ξ))n to zero, a set of algebraic equations for k0,k1,η,σ,ν can be obtained.(86)e0Φξ:λ-νa0+ησa12+a043=0,e-Φξ:λ-νa1+a03a1+σ2+2ηa12=0,e-2Φξ:a062a12+3a12σ=0,e-3Φξ:a0a133+2a12=0. Solving the algebraic equation system, we can get the solutions. Let ρ=±σ; the solutions of (86) are expressed by(87)a0=3ρ,a1=23ρ,ν=λ+34η+ρ4ρ, Different solutions of the auxiliary equation (85) have been given in [68], so we get the solutions of (81). Different cases are discussed as follows.Case 1 (hyperbolic function solutions).
Whenσ2-4η>0,η≠0,(88)A1x,t=3ρ+2ησ2-4ηtanh⁡σ2-4η/2xβ/Γ1+β-λ+34η+ρ4/ρtα/Γ1+α+c+σCase 2 (trigonometric function solutions).
Whenσ2-4η<0,η≠0,(89)A2x,t=3ρ+2η4η-σ2tan4η-σ2/2xβ/Γ1+β-λ+34η+ρ4/ρtα/Γ1+α+c-σ,Case 3 (hyperbolic function solutions).
Whenσ2-4η>0,η=0,σ≠0,(90)A3x,t=3ρ+σcoshσxβ/Γ1+β-λ+34η+ρ4/ρtα/Γ1+α+c+sinhσxβ/Γ1+β-λ+34η+ρ4/ρtα/Γ1+α+c-1,Case 4 (rational function solutions).
Whenσ2-4η=0,η≠0,σ≠0,(91)A4x,t=3ρ-σ2xβ/Γ1+β-λ+34η+ρ4/ρtα/Γ1+α+c2σxβ/Γ1+β-λ+34η+ρ4/ρtα/Γ1+α+c+4,Case 5.
Whenσ2-4η=0,η=0,σ=0,(92)A5x,t=1xβ/Γ1+β-λ+34η+ρ4/ρtα/Γ1+α+c.We take the solution of Case1 as an example to discuss the influence of fractional derivatives α and β on the wave, respectively.Figure1 shows that when the space and time derivative are determined, the amplitude of the wave decreases as time increases.Figure 1
The effect of the variation of the dissipation coefficient on the amplitude if the space and time derivative are determined.
(a)
α = 0.25, β=0.25 (b)
μ = 0.5, β=0.25 (c)
μ = 0.85, β=0.25In Figures2 and 3, we investigate how the parameters α, β affect the nonlinear solitary waves; the corresponding physical interpretations can be given as follows:Figure 2
The amplitude of solitary waveA(x,10) at t=10,β=0.5, for different value of α.Figure 3
The amplitude of solitary waveA(0,t) at x=0,α=0.5, for different value of β.( 1 ) From Figure 2, we can see that the amplitude of solitary wave increases with increase of β value, while the width of the wave decreases.( 2 ) Figure 3 shows that when the values of α increase, the amplitude of solitary waves has an increasing tread; however, the amplitude of the wave declines more and more quickly.


## 6. Conclusions

In this paper, using multiscale analysis and turbulence method, from the basic dynamic multivariable equations under the baroclinic atmosphere, the integer order mKdV equation is derived. In Section3, we use the semi-inverse method and variational method to derive the STFmKdV equation under the Riemann-Liouville definition. In Section 4, we extend the symmetry analysis of the fractional equation to obtain the corresponding infinitesimal generator of the equation. Then we discuss the nonlinear self-adjointness of the STFmKdV equation and finally get the conservation vectors of the equation. In Section 5, based on the STFmKdV equation, employing the exp(-Φ(ξ)) method, and considering the different cases of the parameters (η,σ), we obtain five different solutions of the equation.Note 5.
In this paper, we only study the space-time fractional order equation under the Riemann-Liouville derivative definition. In future studies, we can also consider the fractional equation under the Caputo derivative definition.Note 6.
In future studies, we can also consider the conservation laws of fractional equations under the Caputo derivative definition.



---

*Source: 1345346-2018-10-24.xml*