7.2
CiteScore
3.7
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
ABUNDANCE ESTIMATION IN AN ARID ENVIRONMENT
Case Study
Correspondence
Corrigendum
Editorial
Full Length Article
Invited review
Letter to the Editor
Original Article
Retraction notice
REVIEW
Review Article
SHORT COMMUNICATION
Short review
7.2
CiteScore
3.7
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
ABUNDANCE ESTIMATION IN AN ARID ENVIRONMENT
Case Study
Correspondence
Corrigendum
Editorial
Full Length Article
Invited review
Letter to the Editor
Original Article
Retraction notice
REVIEW
Review Article
SHORT COMMUNICATION
Short review
View/Download PDF

Translate this page into:

Original article
32 (
1
); 567-574
doi:
10.1016/j.jksus.2018.08.013

Sine-Gordon expansion method for exact solutions to conformable time fractional equations in RLW-class

Çankırı Karatekin University, Department of Mathematics, 18200 Çankırı, Turkey
Eskişehir Osmangazi University, Department of Computer and Mathematics, Eskişehir, Turkey
Department of Mathematics, Rasht Branch, Islamic Azad University, Rasht, Iran
Faculty of Engineering Technology, Amol University of Special Modern Technologies, Amol, Iran
Department of Mathematics, Faculty of Mathematical Sciences, University of Mazandaran, Babolsar, Iran

⁎Corresponding author. alperkorkmaz7@gmail.com (Alper Korkmaz)

Disclaimer:
This article was originally published by Elsevier and was migrated to Scientific Scholar after the change of Publisher.

Peer review under responsibility of King Saud University.

Abstract

The Sine-Gordon expansion method is implemented to construct exact solutions some conformable time fractional equations in Regularized Long Wave (RLW)-class. Compatible wave transform reduces the governing equation to classical ordinary differential equation. The homogeneous balance procedure gives the order of the predicted polynomial-type solution that is inspired from well-known Sine-Gordon equation. The substitution of this solution follows the previous step. Equating the coefficients of the powers of predicted solution leads a system of algebraic equations. The solution of resultant system for coefficients gives the necessary relations among the parameters and the coefficients to be able construct the solutions. Some solutions are simulated for some particular choices of parameters.

Keywords

35C07
35R11
35Q53
Sine-Gordon expansion method
Conformable time fractional RLW equation
Conformable time fractional modified RLW equation
Conformable time fractional symmetric-RLW equation
PubMed
1

1 Introduction

Even though some nonlinear partial differential equations are integrable, it may be not easy to integrate them. Instead, a predicted solution with parameters are assumed to be a solution of governing equations and the relations among the parameters are investigated. The logic is simply based on the similarity with exponential-type solutions to the ordinary differential equations with constant coefficients. These predicted solutions are of various forms covering exponential, hyperbolic, trigonometric or rational functions, and more. Moreover the interaction of various types of waves were also determined by using analytical techniques (Ma and Zhou, 2018; Zhang and Ma, 2017; Zhao and Ma, 2017; Yang et al., 2017; Ma et al., 2018). Parallel to the recent developments in computer algebra in the last four decades, a tendency has been observed to determine exact solutions to nonlinear PDEs by following the procedure that starts with a predicted solution. Recently, this tendency has focused on exact solutions to fractional nonlinear partial differential equations. Many of techniques implemented to nonlinear PDEs to find exact solutions have been adapted for fractional nonlinear PDES (Korkmaz and Hosseini, 2017; Hosseini and Ansari, 2017; Kumar et al., 2017; Korkmaz, 2017a,b; Hosseini et al., 2017a,b; Kumar et al., 2018a,b; Rezazadeh et al., 2018; Singh et al., 2017a,b; Guner et al., 2017a,b). We also derive exact solutions to some conformable fractional equations in RLW-class modeling various wave phenomena both in nature or technology implementations. Different from previous studies, we adapt Sine-Gordon expansion approach to determine exact solutions to governing equations in fractional RLW-class.

The first equation considered in this study is the RLW equation

(1)
D t α u + pu x + quu x + rD t α u xx = 0 , t 0 where p , q and r real coefficients, D t α conformable fractional differential operator and u = u ( x , t ) . The integer ordered form of the RLW equation describes formation and development of undular bore by a long wave in shallow water (Peregrine, 1966). The same study investigates the transition and interaction between still water and a uniform flow. The origins and similarities to the KdV Equation of the RLW equation are discussed in general terms (Benjamin et al., 1972). Moreover, Benjamin et al. developed an exact theory to the RLW equation and existence of classical solutions were proved. The Lagrangian density for the RLW equation was defined in Morrison et al. (1984). Significant concept covering peak positions, amplitude and widths of two solitary waves were studied by trial function approach in the same paper.

The second equation to discuss exact solutions here is the modified RLW (mRLW) equation of the form

(2)
D t α u + pu x + qu 2 u x + rD t α u xx = 0 , t 0 where D t α conformable fractional differential operator. This equation seemed in Gardners’ study (Gardner et al., 1997). B-spline finite elements based approximate solutions defining motion of solitary waves were investigated in that study. In the following decades, various numerical and exact solutions to different problems constructed on the mRLW equation are solved using diverse techniques (Soliman, 2017; Khalifa et al., 2008; Santarelli, 1978; Karakoc and Geyikli, 2013; Karakoc et al., 2013; Mei and Chen, 2012).

The last equation to derive the exact solutions is the symmetric RLW (sRLW) equation

(3)
D tt 2 α u + pu xx + quD t α u x + qu x D t α u + rD tt 2 α u xx = 0 , t 0 where u = u ( x , t ) , p , q and r real parameters, D tt 2 α is the fractional differential operator of order 2 α in conformable sense. This equation is an interesting model to describe ion-acoustic and space charge waves with weak non linearity (Seyler and Fenstermacher, 1984).

Before starting the solution procedure, we should give some significant properties of conformable fractional derivative. Thus, the next section focuses on conformable fractional derivative definition and some important properties. We explain the solution procedure in the third section. The following sections cover implementations of the proposed procedure to some conformable fractional PDEs in RLW-class.

2

2 Conformable fractional derivative

The conformable derivative of order α with respect to the independent variable t is defined as

(4)
D t α ( y ( t ) ) = lim τ 0 y ( t + τ t 1 - α ) - y ( t ) τ , t > 0 , α ( 0 , 1 ] . for a function y = y ( t ) : [ 0 , ) R (Khalil et al., 2014). This newly defined fractional derivative is capable of satisfying some well known required properties.
Theorem 1

Assume that the order of the derivative α ( 0 , 1 ] , and suppose that u = u ( t ) and y = y ( t ) are α -differentiable for all positive t. Then,

  • D t α ( c 1 u + c 2 y ) = c 1 D t α ( u ) + c 2 D t α ( y )

  • D t α ( t k ) = kt k - α , k R

  • D t α ( λ ) = 0 , for all constant function u ( t ) = λ

  • D t α ( uy ) = uD t α ( y ) + yD t α ( u )

  • D t α u y = yD t α ( u ) - uD t α ( y ) y 2

  • D t α ( u ) ( t ) = t 1 - α du dt

for c 1 , c 2 R (Atangana et al., 2015).

Conformable fractional differential operator satisfies some critical fundamental properties like the chain rule, Taylor series expansion and Laplace transform (Abdeljawad, 2015).

Theorem 2

Let u = u ( t ) be an α -conformable differentiable function and assume that y is differentiable and defined in the range of u. Then,

(5)
D t α ( u y ) ( t ) = t 1 - α y ( t ) u ( y ( t ) )

3

3 Sine-Gordon expansion method

The compatible fractional form u ( x , t ) = U ( ξ ) with ξ = a ( x - ν t α / α ) of the classical wave transform (Ma and Fuchssteiner, 1996; Ma and Lee, 2009; Ma et al., 2010; Korkmaz, 2018) reduces the fractional Sine–Gordon equation in one dimension of the form

(6)
2 u x 2 - D t 2 α u = m 2 sin u , m is constant to the ODE
(7)
d 2 U d ξ 2 = m 2 a 2 ( 1 - ν 2 ) sin U
where ν represents velocity of the traveling wave defined in the transform (Yan, 1996). Some simplifications lead
(8)
d ( U / 2 ) d ξ 2 = m 2 a 2 ( 1 - ν 2 ) sin 2 U / 2 + C
where C is constant of integration. C is assumed zero for simplicity. Let w ( ξ ) = U ( ξ ) / 2 and b 2 = m 2 / ( a 2 ( 1 - ν 2 ) ) . Then, (8) is converted to
(9)
d ( w ) d ξ = b sin w
Set b = 1 in (9). Then, (9) yields two significant relations
(10)
sin w ( ξ ) = 2 de ξ d 2 e 2 ξ + 1 d = 1 = sech ξ
or
(11)
cos w ( ξ ) = d 2 e 2 ξ - 1 d 2 e 2 ξ + 1 d = 1 = tanh ξ
where d is nonzero integral constant. The fractional PDE of the form
(12)
P ( u , D t α u , u x , D tt 2 α u , u xx , ) = 0
can be reduced to an ODE
(13)
P ̃ ( U , U , U , ) = 0
by using a compatible wave transform u ( x , t ) = U ( ξ ) where the transform variable ξ is defined as a ( x - ν t α / α ) . Then, the predicted solution to (13) of the form
(14)
U ( ξ ) = A 0 + i = 1 s tanh i - 1 ( ξ ) B i sech ξ + A i tanh ξ
can be written as
(15)
U ( w ) = A 0 + i = 1 s cos i - 1 ( w ) B i sin w + A i cos w
owing to (10)–(11). The procedure starts by determining index limit s by the assistance of homogenous balance of the terms in (13). Following the substitution of the predicted solution (15) into (13) the coefficients of powers of sin w cos w are assumed as zero. Next, the resultant algebraic system is tried to be solved for the coefficients A 0 , A 1 , B 1 , , a , ν . Then, the solutions are constructed, if exists, by using (10)–(11) and ξ

4

4 Solutions to the conformable time fractional RLW equation

The traveling wave transform u ( x , t ) U ( ξ ) , ξ = a ( x - ν t α / α ) reduces the time fractional RLW equation to

(16)
- a ν + ap U + 1 / 2 qaU 2 - r ν a 3 U = C where C is constant of integration, and ″ denotes differentiation wrt ξ . Here, we assume that C = 0 to reduce the complexity of the solutions. The balance between U 2 and U gives s = 2 . Thus, the predicted solution takes the form
(17)
U ( w ) = A 0 + B 1 sin w + A 1 cos w + B 2 cos w sin w + A 2 cos 2 w
Substituting this solution into (16) yields
(18)
- a 3 ν rB 2 + aqA 2 B 2 sin w ξ cos w ξ 3 + cos w ξ 3 aqA 1 A 2 + apA 2 - a ν A 2 + 1 / 2 aqA 1 2 + aqA 0 A 2 cos w ξ 2 + - a 3 ν rB 1 + aqA 1 B 2 + aqA 2 B 1 sin w ξ cos w ξ 2 + 1 / 2 aqB 2 2 + 4 a 3 ν rA 2 - 1 / 2 aqA 2 2 sin w ξ 2 cos w ξ 2 + 5 sin w ξ 3 cos w ξ a 3 ν rB 2 + 2 a 3 ν rA 1 + aqB 1 B 2 sin w ξ 2 cos w ξ + aqA 0 B 2 + aqA 1 B 1 - a ν B 2 + apB 2 sin w ξ cos w ξ + aqA 0 A 1 - a ν A 1 + apA 1 cos w ξ + 1 / 2 aqA 2 2 + 1 / 2 aqA 0 2 - a ν A 0 + apA 0 + aqA 0 B 1 - a ν B 1 + apB 1 sin w ξ + sin w ξ 3 a 3 ν rB 1 - 2 sin w ξ 4 a 3 ν rA 2 + 1 / 2 aqB 1 2 - 1 / 2 aqA 2 2 sin w ξ 2 = 0
Using some trigonometric identities and simplifications we find the following algebraic system of equations:
(19)
1 / 2 a - 4 a 2 ν rA 2 - qA 2 2 + qB 1 2 + 1 / 2 aqA 2 2 + 1 / 2 aqA 0 2 - a ν A 0 + apA 0 = 0 a 3 ν rB 1 + aqA 0 B 1 - a ν B 1 + apB 1 = 0 - a ν A 2 + apA 2 + 1 / 2 aqA 1 2 + aqA 0 A 2 + 6 a 3 ν rA 2 + 1 / 2 aqB 2 2 - 1 / 2 aqA 2 2 - 1 / 2 a - 4 a 2 ν rA 2 - qA 2 2 + qB 1 2 = 0 - 6 a 3 ν rA 2 - 1 / 2 aqB 2 2 + 1 / 2 aqA 2 2 = 0 - 2 a 3 ν rA 1 + aqA 1 A 2 - aqB 1 B 2 = 0 2 a 3 ν rA 1 + aqA 0 A 1 + aqB 1 B 2 - a ν A 1 + apA 1 = 0 - 6 a 3 ν rB 2 + aqA 2 B 2 = 0 - 2 a 3 ν rB 1 + aqA 1 B 2 + aqA 2 B 1 = 0 5 a 3 ν rB 2 + aqA 0 B 2 + aqA 1 B 1 - a ν B 2 + apB 2 = 0
Solution of this system for a , ν , A 0 , A 1 , A 2 , B 1 and B 2 gives various solution sets with a 0 and ν 0 .

Set 1:

(20)
ν = p ra 2 + 1 , A 0 = - 6 pa 2 r ( ra 2 + 1 ) q , A 1 = 0 , A 2 = 6 pa 2 r ( ra 2 + 1 ) q , B 1 = 0 , B 2 = 6 ra 2 pi ( ra 2 + 1 ) q , i = - 1

Set 2:

(21)
ν = p ra 2 + 1 , A 0 = - 6 pa 2 r ( ra 2 + 1 ) q , A 1 = 0 , A 2 = 6 pa 2 r ( ra 2 + 1 ) q , B 1 = 0 , B 2 = - 6 ra 2 pi ( ra 2 + 1 ) q , i = - 1

Set 3:

(22)
ν = - p ra 2 - 1 , A 0 = 4 pa 2 r ( ra 2 - 1 ) q , A 1 = 0 , A 2 = - 6 pa 2 r ( ra 2 - 1 ) q , B 1 = 0 , B 2 = 6 ra 2 pi ( ra 2 - 1 ) q , i = - 1

Set 4:

(23)
ν = - p ra 2 - 1 , A 0 = 4 pa 2 r ( ra 2 - 1 ) q , A 1 = 0 , A 2 = - 6 pa 2 r ( ra 2 - 1 ) q , B 1 = 0 , B 2 = - 6 ra 2 pi ( ra 2 - 1 ) q , i = - 1

Set 5:

(24)
ν = p 4 ra 2 + 1 , A 0 = - ( 8 ra 2 - 1 ) p 4 ra 2 + 1 + p q , A 1 = 0 , A 2 = 12 pa 2 r ( 4 ra 2 + 1 ) q , B 1 = 0 , B 2 = 0

Set 6:

(25)
ν = - p 4 ra 2 - 1 , A 0 = - - ( 8 ra 2 - 1 ) p 4 ra 2 - 1 + p q , A 1 = 0 , A 2 = - 12 pa 2 r ( 4 ra 2 - 1 ) q , B 1 = 0 , B 2 = 0

Using these solution sets of algebraic equations, we construct the solutions to (1) as u 1 ( x , t ) = - 6 a 2 r ν q + 6 a 2 r ν q i tanh a x - ν t α α sech a x - ν t α α + 6 a 2 r ν q tanh 2 a x - ν t α α u 2 ( x , t ) = - 6 a 2 r ν q + 6 a 2 r ν q i tanh a x - ν t α α sech a x - ν t α α - 6 a 2 r ν q tanh 2 a ( x - ν t α α ) for ν = p / ( ra 2 + 1 ) , u 3 ( x , t ) = - 4 a 2 r ν q + 6 a 2 r ν q i tanh a x - ν t α α sech a x - ν t α α - 6 a 2 r ν q tanh 2 a x - ν t α α u 4 ( x , t ) = - 4 a 2 r ν q + 6 a 2 r ν q i tanh a x - ν t α α sech a x - ν t α α + 6 a 2 r ν q tanh 2 a x - ν t α α for ν = - p / ( ra 2 - 1 ) , u 5 ( x , t ) = ( 8 ra 2 - 1 ) ν + p q + 12 a 2 r ν q tanh 2 a x - ν t α α for ν = p / ( 4 ra 2 + 1 ) . u 6 ( x , t ) = - ( 8 ra 2 - 1 ) ν + p q + 12 a 2 r ν q tanh 2 a ( x - ν t α α ) for ν = - p / ( 4 ra 2 + 1 ) .

A particular form of u 5 ( x , t ) is depicted for various values of α { 0.25 , 0.5 , 0.75 , 1 } in a finite domain in Fig. 1(a)–(d). It is observed that the propagation of the initial pulse propagates along the x -axis as time proceeds by preserving its shape and amplitude in all cases. α affects only propagation velocity. The propagation is faster in small times but later it gets slower when α is less than 1. α = 1 choice gives a constant propagation velocity to the pulse, Fig. 1(d).

The solution u 5 ( x , t ) for { a = 1 / 4 , q = - 1 , r = 1 , p = 10 } .
Fig. 1
The solution u 5 ( x , t ) for { a = 1 / 4 , q = - 1 , r = 1 , p = 10 } .

5

5 Solutions to the conformable time fractional mRLW equation

The traveling wave transform u ( x , t ) U ( ξ ) , ξ = a ( x - ν t α / α ) reduces the mRLW Eq. (2) to

(26)
- a ν + ap U + qa 3 U 3 - r ν a 3 U = C where C is constant of integration. The balance between U 3 and U leads s = 1 . Thus, the predicted solution takes the form
(27)
U ( w ) = A 0 + A 1 cos w + B 1 sin w
Substituting the predicted solution (27) into (26) gives
(28)
1 / 3 cos w ξ 3 aqA 1 3 - a 3 ν rB 1 + aqA 1 2 B 1 sin w ξ cos w ξ 2 + cos w ξ 2 aqA 0 A 1 2 + 2 a 3 ν rA 1 + aqA 1 B 1 2 sin w ξ 2 cos w ξ + 2 sin w ξ cos w ξ aqA 0 A 1 B 1 + aqA 0 2 A 1 - a ν A 1 + apA 1 cos w ξ - C - a ν A 0 + apA 0 + 1 / 3 aqA 0 3 + a 3 ν rB 1 + 1 / 3 aqB 1 3 sin w ξ 3 + sin w ξ 2 aqA 0 B 1 2 + aqA 0 2 B 1 - a ν B 1 + apB 1 sin w ξ = 0
Following implementation of some trigonometric identities and some simplifications, we equate the coefficients of powers of sin , cos and multiplications of them to zero to give the algebraic system of equations
(29)
aqA 0 B 1 2 - C - a ν A 0 + apA 0 + 1 / 3 aqA 0 3 = 0 aqA 0 2 B 1 - a ν B 1 + apB 1 + a 3 ν rB 1 + 1 / 3 aqB 1 3 = 0 aqA 0 A 1 2 - aqA 0 B 1 2 = 0 - 2 a 3 ν rB 1 + aqA 1 2 B 1 - 1 / 3 aqB 1 3 = 0 - 2 a 3 ν rA 1 - aqA 1 B 1 2 + 1 / 3 aqA 1 3 = 0 2 aqA 0 A 1 B 1 = 0 2 a 3 ν rA 1 + aqA 0 2 A 1 + aqA 1 B 1 2 - a ν A 1 + apA 1 = 0
The solution of this system for ν , A 0 , A 1 , B 1 and C forces C to be zero. Thus, the solutions of this system can be summarized as
(30)
ν = - - 2 pq qra 2 - 2 q , A 0 = 0 , A 1 = 3 r ν 2 q a , B 1 = - 3 r ν 2 q a ν = - - 2 pq qra 2 - 2 q , A 0 = 0 , A 1 = 3 r ν 2 q a , B 1 = - - 3 r ν 2 q a ν = - - 2 pq qra 2 - 2 q , A 0 = 0 , A 1 = - 3 r ν 2 q a , B 1 = - 3 r ν 2 q a ν = - - 2 pq qra 2 - 2 q , A 0 = 0 , A 1 = - 3 r ν 2 q a , B 1 = - - 3 r ν 2 q a ν = p ra 2 + 1 , A 0 = 0 , A 1 = 0 , B 1 = 6 r ν q a ν = p ra 2 + 1 , A 0 = 0 , A 1 = 0 , B 1 = - 6 r ν q a ν = - p 2 ra 2 - 1 , A 0 = 0 , A 1 = - 6 r ν q a , B 1 = 0 ν = - p 2 ra 2 - 1 , A 0 = 0 , A 1 = - - 6 r ν q a , B 1 = 0
for arbitrary a. When a is considered as non arbitrary, we have some more solutions such that
(31)
a = 1 r , ν = 2 p , A 0 = 0 , A 1 = 3 p q , B 1 = - 3 p q a = - 1 r , ν = 2 p , A 0 = 0 , A 1 = 3 p q , B 1 = - 3 p q a = 1 r , ν = 2 p , A 0 = 0 , A 1 = - 3 p q , B 1 = - 3 p q a = - 1 r , ν = 2 p , A 0 = 0 , A 1 = - 3 p q , B 1 = - 3 p q a = 1 r , ν = 2 p , A 0 = 0 , A 1 = 3 p q , B 1 = - - 3 p q a = - 1 r , ν = 2 p , A 0 = 0 , A 1 = 3 p q , B 1 = - - 3 p q a = 1 r , ν = 2 p , A 0 = 0 , A 1 = - 3 p q , B 1 = - - 3 p q a = - 1 r , ν = 2 p , A 0 = 0 , A 1 = - 3 p q , B 1 = - - 3 p q
Thus, the solutions to mRLW Eq. (2) are represented by
(32)
u 7 ( x , t ) = 3 r ν 2 q a tanh a x - ν t α α + - 3 r ν 2 q a sech a x - ν t α α , ν = - - 2 pq qra 2 - 2 q u 8 ( x , t ) = 3 r ν 2 q a tanh a x - ν t α α - - 3 r ν 2 q a sech a x - ν t α α , ν = - - 2 pq qra 2 - 2 q u 9 ( x , t ) = - 3 r ν 2 q a tanh a x - ν t α α + - 3 r ν 2 q a sech a x - ν t α α , ν = - - 2 pq qra 2 - 2 q u 10 ( x , t ) = - 3 r ν 2 q a tanh a x - ν t α α - - 3 r ν 2 q a sech a x - ν t α α , ν = - - 2 pq qra 2 - 2 q u 11 ( x , t ) = 6 r ν q a sech a x - ν t α α , ν = p ra 2 + 1 u 12 ( x , t ) = - 6 r ν q a sech a x - ν t α α , ν = p ra 2 + 1 u 13 ( x , t ) = - 6 r ν q a tanh a x - ν t α α , ν = - - p 2 ra 2 - 1 u 14 ( x , t ) = - - 6 r ν q a tanh a x - ν t α α , ν = - - p 2 ra 2 - 1
when ν is dependent on a. Whenever ν and a are independent upon each other, the solution are represented by
(33)
u 15 ( x , t ) = 3 p q tanh a x - ν t α α + - 3 p q sech a x - ν t α α , ν = 2 p , a = 1 r u 16 ( x , t ) = 3 p q tanh a x - ν t α α + - 3 p q sech a x - ν t α α , ν = 2 p , a = - 1 r u 17 ( x , t ) = - 3 p q tanh a x - ν t α α + - 3 p q sech a x - ν t α α , ν = 2 p , a = 1 r u 18 ( x , t ) = - 3 p q tanh a x - ν t α α + - 3 p q sech a x - ν t α α , ν = 2 p , a = - 1 r u 19 ( x , t ) = 3 p q tanh a x - ν t α α - - 3 p q sech a x - ν t α α , ν = 2 p , a = 1 r u 20 ( x , t ) = 3 p q tanh a x - ν t α α - - 3 p q sech a x - ν t α α , ν = 2 p , a = - 1 r u 21 ( x , t ) = - 3 p q tanh a x - ν t α α - - 3 p q sech a x - ν t α α , ν = 2 p , a = 1 r u 22 ( x , t ) = - 3 p q tanh a x - ν t α α - - 3 p q sech a x - ν t α α , ν = 2 p , a = - 1 r
The solution u 11 ( x , t ) is simulated for various values of α in Fig. 2(a)–(d) by using the parameter set { a = 1 , q = - 1 , r = 1 / 3 , p = 10 } . The solution propagates along x -axis as time proceeds with preserved shape and amplitude. The propagation is not linear when α is less than 1, Fig. 2(a)–(c) as we observe a linear propagation with constant speed with α = 1 , Fig. 2(d).
The solution u 11 ( x , t ) for { a = 1 , q = - 1 , r = 1 / 3 , p = 10 } .
Fig. 2
The solution u 11 ( x , t ) for { a = 1 , q = - 1 , r = 1 / 3 , p = 10 } .

6

6 Solutions to the conformable time fractional sRLW equation

The traveling wave transform u ( x , t ) U ( ξ ) , ξ = a ( x - ν t α / α ) reduces the sRLW Eq. (2) to

(34)
a 2 ν 2 + a 2 p U + r ν 2 a 4 U - q ν a 2 UU = C The balance between U and UU gives s = 2 . Thus, we seek a solution of the form
(35)
U ( w ) = A 0 + A 1 cos w + B 1 sin w + B 2 cos w sin w + A 2 cos 2 w
Substitution of this solution into (34) gives
(36)
a 4 ν 2 rB 2 sin w ξ cos w ξ 4 + - 8 a 4 ν 2 rA 2 + 2 a 2 ν qA 2 2 - a 2 ν qB 2 2 sin w ξ 2 cos w ξ 3 + a 4 ν 2 rB 1 - a 2 ν qA 1 B 2 - a 2 ν qA 2 B 1 sin w ξ cos w ξ 3 + - 18 a 4 ν 2 rB 2 + 4 a 2 ν qA 2 B 2 sin w ξ 3 cos w ξ 2 + - 4 a 4 ν 2 rA 1 + 3 a 2 ν qA 1 A 2 - 2 a 2 ν qB 1 B 2 sin w ξ 2 cos w ξ 2 + - a 2 ν qA 0 B 2 - a 2 ν qA 1 B 1 + a 2 ν 2 B 2 + a 2 pB 2 sin w ξ cos w ξ 2 + 16 a 4 ν 2 rA 2 + a 2 ν qB 2 2 sin w ξ 4 cos w ξ + - 5 a 4 ν 2 rB 1 + 2 a 2 ν qA 1 B 2 + 2 a 2 ν qA 2 B 1 sin w ξ 3 cos w ξ + 2 a 2 ν qA 0 A 2 + a 2 ν qA 1 2 - a 2 ν qB 1 2 - 2 a 2 ν 2 A 2 - 2 a 2 pA 2 sin w ξ 2 cos w ξ + - a 2 ν qA 0 B 1 + a 2 ν 2 B 1 + a 2 pB 1 sin w ξ cos w ξ + 5 sin w ξ 5 a 4 ν 2 rB 2 + 2 a 4 ν 2 rA 1 + a 2 ν qB 1 B 2 sin w ξ 4 + a 2 ν qA 0 B 2 + a 2 ν qA 1 B 1 + a 2 ν qA 2 B 2 - a 2 ν 2 B 2 - a 2 pB 2 sin w ξ 3 + a 2 ν qA 0 A 1 - a 2 ν 2 A 1 - a 2 pA 1 sin w ξ 2 - a 2 ν qA 2 B 2 sin w ξ - C = 0
Substitution of some trigonometric identities and some simplifications leads the following system of algebraic equations from the equating the coefficients of powers of cos , sin functions and their multiplications:
(37)
- 3 a 2 ν 2 a 2 ν rA 1 - qA 1 A 2 + qB 1 B 2 + 8 a 4 ν 2 rA 1 + a 2 ν qA 0 A 1 - 3 a 2 ν qA 1 A 2 + 4 a 2 ν qB 1 B 2 - a 2 ν 2 A 1 - a 2 pA 1 - C = 0 5 a 4 ν 2 rB 2 + a 2 ν qA 0 B 2 + a 2 ν qA 1 B 1 - a 2 ν 2 B 2 - a 2 pB 2 = 0 24 a 4 ν 2 rB 2 - 4 a 2 ν qA 2 B 2 = 0 6 a 4 ν 2 rB 1 - 3 a 2 ν qA 1 B 2 - 3 a 2 ν qA 2 B 1 = 0 - 28 a 4 ν 2 rB 2 - 2 a 2 ν qA 0 B 2 - 2 a 2 ν qA 1 B 1 + 3 a 2 ν qA 2 B 2 + 2 a 2 ν 2 B 2 + 2 a 2 pB 2 = 0 - 5 a 4 ν 2 rB 1 - a 2 ν qA 0 B 1 + 2 a 2 ν qA 1 B 2 + 2 a 2 ν qA 2 B 1 + a 2 ν 2 B 1 + a 2 pB 1 = 0 a 2 24 a 2 ν 2 rA 2 - 2 ν qA 2 2 + 2 ν qB 2 2 = 0 6 a 4 ν 2 rA 1 - 3 a 2 ν qA 1 A 2 + 3 a 2 ν qB 1 B 2 = 0 a 2 - 16 a 2 ν 2 rA 2 - 2 ν qA 0 A 2 - ν qA 1 2 + ν qB 1 2 - ν qB 2 2 + 2 ν 2 A 2 + 2 pA 2 - a 2 24 a 2 ν 2 rA 2 - 2 ν qA 2 2 + 2 ν qB 2 2 = 0 - 14 a 4 ν 2 rA 1 - a 2 ν qA 0 A 1 + 6 a 2 ν qA 1 A 2 - 7 a 2 ν qB 1 B 2 + a 2 ν 2 A 1 + a 2 pA 1 + 3 a 2 ν 2 a 2 ν rA 1 - qA 1 A 2 + qB 1 B 2 = 0 - a 2 - 16 a 2 ν 2 rA 2 - 2 ν qA 0 A 2 - ν qA 1 2 + ν qB 1 2 - ν qB 2 2 + 2 ν 2 A 2 + 2 pA 2 = 0
Solution of this system gives
(38)
A 0 = - 5 a 2 ν 2 r - ν 2 - p ν q , A 1 = 0 , A 2 = 12 a 2 ν r q , B 1 = 0 , B 2 = 0 A 0 = - 5 a 2 ν 2 r - ν 2 - p ν q , A 1 = 0 , A 2 = 6 a 2 ν r q , B 1 = 0 , B 2 = 6 a 2 r ν i q A 0 = - 5 a 2 ν 2 r - ν 2 - p ν q , A 1 = 0 , A 2 = 6 a 2 ν r q , B 1 = 0 , B 2 = - 6 a 2 r ν i q
for arbitrarily chosen ν , a and C = 0 . One should note that the solution of this system is determinable whenever C = 0 . Thus, the solutions to (3) are constructed as
(39)
u 23 ( x , t ) = - 5 a 2 ν 2 r - ν 2 - p ν q + 6 a 2 r ν q i tanh a ( x - ν t α α ) sech a ( x - ν t α α ) + 6 a 2 r ν q tanh 2 a ( x - ν t α α ) u 24 ( x , t ) = - 5 a 2 ν 2 r - ν 2 - p ν q + 6 a 2 r ν q i tanh a ( x - ν t α α ) sech a ( x - ν t α α ) - 6 a 2 r ν q tanh 2 a ( x - ν t α α ) u 25 ( x , t ) = - 8 a 2 ν 2 r - ν 2 - p ν q + 12 a 2 ν r q tanh 2 a ( x - ν t α α )
The solution u 25 ( x , t ) is depicted in Fig. 3(a)–(d) for various α values in some finite domain of independent variables. The initial pulse propagates along the space axis without changing its shape and amplitude in all cases. α affects only propagation velocity due to being multiplier of time variable. Even though the pulse propagates non linearly for smaller α values, the propagation is linear when α is 1.
The solution u 25 ( x , t ) for { a = 1 , q = 1 , r = - 1 / 3 , p = 1 , ν = 5 } .
Fig. 3
The solution u 25 ( x , t ) for { a = 1 , q = 1 , r = - 1 / 3 , p = 1 , ν = 5 } .

7

7 Conclusion

In the paper, exact solutions of some conformable fractional equations in the RLW-class are investigated by using Sine-Gordon expansion approach. Using compatible wave transform, the equations are reduced to some ODEs. Then, the predicted solutions are substituted into the resultant ODE. Equating the coefficients of cosine and sine functions and their multiplications to zero leads to some algebraic system of equation. Solving this system gives the relations among the parameters. The method differs from many classical method since the terms in the series are multiplication of powers of different hyperbolic functions. In conclusion, some real and complex solutions that are combinations of powers of hyperbolic tangent and hyperbolic secant functions are determined explicitly. Graphical representations of some real valued solutions are depicted in some finite domains to comprehend the effects of α .

References

  1. , . On conformable fractional calculus. J. Comput. Appl. Math.. 2015;279:57-66.
    [Google Scholar]
  2. , , , . New properties of conformable derivative. Open Math.. 2015;13(1):1-10.
    [Google Scholar]
  3. , , , . Model equations for long waves in nonlinear dispersive systems. Philos. Trans. R. Soc. London A: Math., Phys. Eng. Sci.. 1972;272(1220):47-78.
    [Google Scholar]
  4. , , , , . Approximations of solitary waves of the MRLW equation by B-spline finite element. Arab. J. Sci. Eng.. 1997;22:183-193.
    [Google Scholar]
  5. , , , . Tanh-type and sech-type solitons for some space-time fractional PDE models. Eur. Phys. J. Plus. 2017;132(92):1-12.
    [Google Scholar]
  6. , , , . Dark solutions of space-time fractional Sharma-Tasso-Olver and potential Kadomtsev-Petviashvili Equations. Commun. Theor. Phys.. 2017;67(2):182-188.
    [Google Scholar]
  7. , , . New exact solutions of nonlinear conformable time-fractional Boussinesq equations using the modified Kudryashov method. Waves Random Complex Media 2017:1-9.
    [Google Scholar]
  8. , , , . Bright and singular soliton solutions of the conformable time-fractional Klein-Gordon equations with different nonlinearities. Waves Random Complex Media 2017:1-9.
    [Google Scholar]
  9. , , , . Modified Kudryashov method for solving the conformable time-fractional Klein-Gordon equations with quadratic and cubic nonlinearities. Optik-Int. J. Light Electron Opt.. 2017;130:737-742.
    [Google Scholar]
  10. , , . Petrov-Galerkin finite element method for solving the MRLW equation. Math. Sci.. 2013;7(1):25.
    [Google Scholar]
  11. , , , . Numerical approximation to a solution of the modified regularized long wave equation using quintic B-splines. Boundary Value Prob.. 2013;2013(1):27.
    [Google Scholar]
  12. , , , . A collocation method with cubic B-splines for solving the MRLW equation. J. Comput. Appl. Math.. 2008;212(2):406-418.
    [Google Scholar]
  13. , , , , . A new definition of fractional derivative. J. Comput. Appl. Math.. 2014;264:65-70.
    [Google Scholar]
  14. , . Exact solutions of space-time fractional EW and modified EW equations. Chaos, Solitons Fractals. 2017;96:132-138.
    [Google Scholar]
  15. , . Exact solutions to (3 + 1) conformable time fractional Jimbo-Miwa, Zakharov-Kuznetsov and modified Zakharov-Kuznetsov equations. Commun. Theor. Phys.. 2017;67(5):479-482.
    [Google Scholar]
  16. , . Complex wave solutions to mathematical biology models I: Newell-Whitehead-Segel and Zeldovich equations. . 2018;13(8):081004-1. 1
    [Google Scholar]
  17. , , . Exact solutions of a nonlinear conformable time-fractional parabolic equation with exponential nonlinearity using reliable methods. Opt. Quant. Electron.. 2017;49(8):278.
    [Google Scholar]
  18. , , , . A new analysis for fractional model of regularized long wave equation arising in ion acoustic plasma waves. Math. Methods Appl. Sci.. 2017;40(15):5642-5653.
    [Google Scholar]
  19. , , , . Analysis of regularized long-wave equation associated with a new fractional operator with Mittag-Leffler type kernel. Physica A. 2018;492:155-167.
    [Google Scholar]
  20. , , , . A new numerical algorithm for fractional Fitzhugh-Nagumo equation arising in transmission of nerve impulses. Nonlinear Dyn.. 2018;91(1):307-317.
    [Google Scholar]
  21. , , . Explicit and exact solutions to a Kolmogorov-Petrovskii-Piskunov equation. Int. J. Non-Linear Mech.. 1996;31(3):329-338.
    [Google Scholar]
  22. , , . A transformed rational function method and exact solutions to the 3 + 1 dimensional Jimbo-Miwa equation. Chaos, Solitons Fractals. 2009;42(3):1356-1363.
    [Google Scholar]
  23. , , . Lump solutions to nonlinear partial differential equations via Hirota bilinear forms. J. Diff. Equ.. 2018;264(4):2633-2659.
    [Google Scholar]
  24. , , , . A multiple exp-function method for nonlinear differential equations and its application. Phys. Scr.. 2010;82(6) 065003
    [Google Scholar]
  25. , , , . Diversity of interaction solutions to the (2 + 1)-dimensional Ito equation. Comput. Math. Appl.. 2018;75:289-295.
    [Google Scholar]
  26. , , . Numerical solutions of RLW equation using Galerkin method with extrapolation techniques. Comput. Phys. Commun.. 2012;183(8):1609-1616.
    [Google Scholar]
  27. , , , . Scattering of regularized-long-wave solitary waves. Physica D. 1984;11(3):324-336.
    [Google Scholar]
  28. , . Calculations of the development of an undular bore. J. Fluid Mech.. 1966;25(02):321-330.
    [Google Scholar]
  29. , , , , , . Traveling wave solution of conformable fractional generalized reaction Duffing model by generalized projective Riccati equation method. Opt. Quant. Electron.. 2018;50(3):150.
    [Google Scholar]
  30. , . Numerical analysis of the regularized long-wave equation: anelastic collision of solitary waves. Il Nuovo Cimento B (1971–1996). 1978;46(1):179-188.
    [Google Scholar]
  31. , , . A symmetric regularized-long-wave equation. Phys. Fluids (1958-1988). 1984;27(1):4-7.
    [Google Scholar]
  32. , , , , . A novel numerical approach for a nonlinear fractional dynamical model of interpersonal and romantic relationships. Entropy. 2017;19(7):375.
    [Google Scholar]
  33. , , , . A hybrid computational approach for Jeffery-Hamel flow in non-parallel walls. Neural Comput. Appl. 2017:1-7.
    [Google Scholar]
  34. , . Collocation method using quartic B-splines for solving the modified RLW equation. Indian J. Sci. Technol.. 2017;10(31)
    [Google Scholar]
  35. , . J. Comput. Nonlinear Dyn.. 1996;224:77-84.
  36. , , , . Lump and lump-soliton solutions to the (2 + 1)-dimensional Ito equation. Anal. Math. Phys. 2017:1-10.
    [Google Scholar]
  37. , , . Mixed lump-kink solutions to the BKP equation. Comput. Math. Appl.. 2017;74(3):591-596.
    [Google Scholar]
  38. , , . Mixed lump-kink solutions to the KP equation. Comput. Math. Appl.. 2017;74(6):1399-1405.
    [Google Scholar]
Show Sections