Next Article in Journal
The Generalized Gielis Geometric Equation and Its Application
Previous Article in Journal
An Ontology Based on the Timeline of Log2timeline and Psort Using Abstraction Approach in Digital Forensics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Exact Solutions of the Conformable Space-Time Sharma–Tasso–Olver Equation Using Two Reliable Methods

by
Sekson Sirisubtawee
1,2,*,
Sanoe Koonprasert
1,2 and
Surattana Sungnul
1,2
1
Department of Mathematics, Faculty of Applied Science, King Mongkut’s University of Technology North Bangkok, Bangkok 10800, Thailand
2
Centre of Excellence in Mathematics, CHE, Si Ayutthaya Road, Bangkok 10400, Thailand
*
Author to whom correspondence should be addressed.
Submission received: 18 March 2020 / Revised: 9 April 2020 / Accepted: 13 April 2020 / Published: 17 April 2020

Abstract

:
The major purpose of this article is to seek for exact traveling wave solutions of the nonlinear space-time Sharma–Tasso–Olver equation in the sense of conformable derivatives. The novel ( G G ) -expansion method and the generalized Kudryashov method, which are analytical, powerful, and reliable methods, are used to solve the equation via a fractional complex transformation. The exact solutions of the equation, obtained using the novel ( G G ) -expansion method, can be classified in terms of hyperbolic, trigonometric, and rational function solutions. Applying the generalized Kudryashov method to the equation, we obtain explicit exact solutions expressed as fractional solutions of the exponential functions. The exact solutions obtained using the two methods represent some physical behaviors such as a singularly periodic traveling wave solution and a singular multiple-soliton solution. Some selected solutions of the equation are graphically portrayed including 3-D, 2-D, and contour plots. As a result, some innovative exact solutions of the equation are produced via the methods, and they are not the same as the ones obtained using other techniques utilized previously.

1. Introduction

Many nonlinear physical phenomena, such as those found in solid-state physics, plasma physics, optical fibers, shallow water waves, fluid dynamics, and biology, have been modeled by nonlinear partial differential equations (NPDEs) of integer- or fractional-order. Therefore, finding explicit exact solutions and approximate analytical solutions of NPDEs is one of the most significant and active areas of investigation in pure and applied mathematics. Several powerful methods based on using symbolic software packages such as Mathematica or Maple 17 have been proposed and developed to analytically solve NPDEs. Some examples of approaches for obtaining approximate solutions in an analytical form to NPDEs are the Adomian decomposition method (ADM) [1,2], the variational iteration method (VIM) [3,4], the differential transform method (DTM) [5], and the homotopy perturbation method (HPM) [6,7]. Some algorithms for obtaining explicit exact solutions of NPDEs are, for instance, the ( G G ) -expansion method [8,9], the ( G G , 1 G ) -expansion method [10], the fractional Riccati expansion method [11,12], the improved extended tanh-coth method [13], the Kudryashov method [14,15], and the sub-equation method [16,17]. All of the above methods are based on the homogeneous balance principle.
In particular, the novel ( G G ) -expansion method and the generalized Kudryashov method have been widely applied to NPDEs for finding explicit exact solutions. The utilizations of the novel ( G G ) -expansion method for solving NPDEs are reviewed as follows. The novel ( G G ) -expansion method is a generalization of many kinds of the ( G G ) -expansion methods such as the basic ( G G ) -expansion method proposed by Wang et al. [18], the improved ( G G ) -expansion method introduced by Zhang et al. [19], and the generalized and improved ( G G ) -expansion method presented by Akbar et al. [20]. Later on, Hafez et al. [21] employed the novel ( G G ) -expansion method to obtain exact traveling wave solutions of the Klein–Gordon equation. Alam et al. obtained explicit exact solutions of the Boussinesq equation [22] and the ( 1 + 1 ) -dimensional modified Benjamin–Bona–Mahony equation [23] using the novel ( G G ) -expansion method. In 2015, Hafez et al. [24] employed the method to construct the exact solutions of the ( 2 + 1 ) -dimensional coupled integrable nonlinear Maccari system. Alam et al. [25] used the method to find explicit exact traveling wave solutions of the ( 1 + 1 ) -dimensional combined KdV-mKdV equation. They found that the solutions of the equation include solitary wave solutions, kink-type solutions, and periodic solutions. Moreover, Alam et al. [26] demonstrated applications of the novel ( G G ) -expansion method for obtaining novel exact traveling wave solutions of the complex coupled Higgs field equations. In 2016, Hafez [27] employed the novel ( G G ) -expansion method to construct some new traveling wave solutions of the ( 1 + 1 ) -dimensional cubic nonlinear Schrödinger’s equation. Exact solutions of the ( 1 + 1 ) -dimensional compound KdVB equation obtained by the novel ( G G ) -expansion method were studied by Alam and Muhammad Belgacem [28]. The novel ( G G ) -expansion method was implemented to find traveling wave solutions of the positive Gardner-KP equation by Akbar et al. [29]. The recent review of the applications of the generalized Kudryashov method is as follows. In 2015, the ( 3 + 1 ) -dimensional Kadomtsev–Petviashvili equation, the ( 3 + 1 ) -dimensional Jimbo–Miwa equation, and the ( 3 + 1 ) -dimensional Zakharov–Kuznetsov equation were studied by Islam et al. [30] to achieve their exact solutions through the generalized Kudryashov method. In 2016, Kaplan [31] utilized the method to seek for exact solutions of the nonlinear Jaulent–Miodek hierarchy and the ( 2 + 1 ) -dimensional Calogero–Bogoyavlenskii–Schiff equation. In 2017, Mahmud et al. [32] applied the generalized Kudryashov method to construct exact traveling wave solutions of the Fisher equation and the PHI-four equation. In 2018, Biswas et al. [33] obtained dark, bright, and singular soliton solutions of the perturbed nonlinear Schrödinger’s equation with fractional temporal evolution using the generalized Kudryashov method. Nestor et al. [34] investigated that the generalized Kudryashov method provides solitary wave solutions and hyperbolic function solutions for the dimensionless Schrödinger’s equation with quadratic-cubic nonlinearity. Houwe et al. [35] used the method to obtain dark and bright exact solutions for the perturbed nonlinear Schrödinger’s equation describing the dynamics of ultrashort optical solitons. Through applying the generalized Kudryashov and the novel G G -expansion methods to the nonlinear complex fractional generalized-Zakharov system, Lu et al. [36] obtained its new forms of analytical and solitary traveling wave solutions. In 2019, Demiray and Bulut [37] applied the generalized Kudryashov method to search for dark and bright soliton solutions of the coupled Higgs equation and the Nizhnik–Novikov–Veselov (NNV) system. In addition, Habib et al. [38] employed the generalized Kudryashov method to attain exact solutions of the Burgers-Huxley, the mKdV, and the first extended fifth-order nonlinear equations.
However, the objective of this paper is to apply the novel ( G G ) -expansion method and the generalized Kudryashov method to analytically solve the conformable space-time Sharma–Tasso–Olver (STO) equation in the sense of conformable derivatives to obtain a considerable number of exact solutions. The STO equation has been utilized in many scientific applications such as in nonlinear optics, dispersive wave phenomena, plasma physics, quantum field theory, and physical sciences (see [39,40,41,42,43] and references therein). To the best of the authors’ knowledge, there are no researchers who have utilized these methods to solve the STO equation with conformable derivatives for obtaining its exact solutions. Some of the exact solutions of the equation, which will be derived in the following sections, are new and reported in this paper for the first time. The rest of this paper is arranged as follows. In Section 2, the definition and some properties of the conformable derivative are given. The description of the two methods is included in this section as well. In Section 3, we illustrate the application of the used methods to the nonlinear conformable space-time Sharma–Tasso–Olver equation. In Section 4, we give some plots and their physical explanations of some chosen exact solutions of the equation. Some conclusions and discussions for the results obtained using the methods are given in Section 5.

2. Mathematical Preliminaries

In this section, we will provide fundamental concepts including the definition and vital properties of the conformable derivative and the algorithms of the novel ( G G ) -expansion method and the generalized Kudryashov method. They are required for constructing explicit exact solutions of the conformable space-time Sharma–Tasso–Olver equation using the methods.

2.1. Conformable Derivative and Its Properties

The definition of the conformable derivative and its important properties are provided as follows.
Definition 1.
Let f : [ 0 , ) R be a function. Then, the conformable derivative of f of order α is defined by [44,45]
D t α f ( t ) = lim ε 0 f ( t + ε t 1 α ) f ( t ) ε , f o r a l l t > 0 , 0 < α 1 .
The function f is α-conformable differentiable at a point t > 0 if the limit in Equation (1) exists.
Remark 1.
The derivative defined in Equation (1) was initially called the conformable fractional derivative and had been utilized in numerous applications of fractional differential equations (FDEs) [46,47,48]. Until 2018, Tarasov [49] showed that the conformable fractional derivative in Equation (1) does not provide innovative ideas in the spaces of differentiable functions and is not a fractional-order derivative. In particular, some of the following properties demonstrate that the conformable fractional derivative can be written in terms of an ordinary derivative. Throughout this work, we thus call the derivative in Equation (1) the conformable derivative.
Theorem 1.
[44,45] Let α ( 0 , 1 ] , and suppose functions f ( t ) , g ( t ) are α-conformable differentiable at a point t > 0 . Then,
D t α ( λ ) = 0 , where λ = constant , D t α ( t μ ) = μ t μ α , for all μ R , D t α ( a f ( t ) + b g ( t ) ) = a D t α f ( t ) + b D t α g ( t ) , for all a , b R , D t α ( f ( t ) g ( t ) ) = f ( t ) D t α g ( t ) + g ( t ) D t α f ( t ) , D t α f ( t ) g ( t ) = g ( t ) D t α f ( t ) f ( t ) D t α g ( t ) g ( t ) 2 .
Remark 2.
Conformable derivatives of some interesting functions are as follows [44].
(1) 
D t α ( e c t ) = c t 1 α e c t , c R .
(2) 
D t α ( sin b t ) = b t 1 α cos b t , b R .
(3) 
D t α ( cos b t ) = b t 1 α sin b t , b R .
(4) 
D t α ( 1 α t α ) = 1 .
(5) 
D t α ( f ( t ) ) = t 1 α d f ( t ) d t , provided that f ( t ) is differentiable.
Theorem 2.
[45,46] Suppose f : ( 0 , ) R is a function such that f is differentiable and α-conformable differentiable. Further, suppose that g is a differentiable function defined in the range of f. Then,
D t α ( f g ) ( t ) = t 1 α f ( g ( t ) ) g ( t ) ,
where the prime notation ( ) denotes the ordinary derivative.

2.2. Description of the Methods

Consider the following nonlinear evolution partial differential equation for two independent variables x and t,
F ( u , D t α u , D x β u , D t 2 α u , D t α D x β u , D x 2 β u , ) = 0 ,
where D t α u , D x β u are the conformable derivatives of a dependent variable u with respect to variables t and x, respectively, with 0 < α , β < 1 . The function F is a polynomial of the unknown function u = u ( x , t ) , and its diverse partial derivatives in which the highest order derivatives and nonlinear terms are involved. There is a common step between the novel ( G G ) -expansion method and the generalized Kudryashov method which is a conversion from the partial differential equation in (3) to an ordinary differential equation (ODE) via using a traveling wave variable ξ [50,51,52]. We suppose that
U ( ξ ) = u ( x , t ) , ξ = x β β V t α α ,
where V is a non-zero arbitrary constant. Converting (3) via transformation (4) and then integrating the resulting equation with respect to ξ (if possible), Equation (3) is reduced to an ODE in the variable U = U ( ξ ) as
P ( U , U , U , U , ) = 0 ,
where P is a polynomial function of U ( ξ ) and its various derivatives. The prime notation (′) represents the ordinary derivative with respect to ξ .
Next, we provide the remaining steps of each of the methods.

2.2.1. Description of the Novel ( G G ) -Expansion Method

Step 1: Suppose that a solution of Equation (5) can be expressed in powers of ψ ( ξ ) as follows,
U ( ξ ) = j = N N a j ψ ( ξ ) j ,
where
ψ ( ξ ) = d + ϕ ( ξ ) with ϕ ( ξ ) = G ( ξ ) G ( ξ ) .
The unknown constants a N or a N may be zero, but both of them cannot be zero simultaneously. The constants a j ( j = 0 , ± 1 , ± 2 , , ± N ) and d are computed at a following step and the function G = G ( ξ ) satisfies the following nonlinear second-order ODE,
G G = λ G G + μ G 2 + v G 2 ,
where the prime notation (′) is the ordinary derivative with respect to ξ , and where λ , μ , and v are real parameters.
Using the Cole–Hopf transformation ϕ ( ξ ) = ln ( G ( ξ ) ) ξ = G ( ξ ) / G ( ξ ) , Equation (8) can be reduced into the generalized Riccati equation as follows,
ϕ ( ξ ) = μ + λ ϕ ( ξ ) + ( v 1 ) ϕ 2 ( ξ ) .
It has been discovered that Equation (9) has thirty-nine exact solutions (see [53,54,55] for the details).
Step 2: The value of the positive integer N in solution (6) can be computed by balancing the highest-order derivative term with the nonlinear terms of the highest order occurring in Equation (5). The degree formulas of some terms are given as
Deg d q U ( ξ ) d ξ q = N + q , Deg ( U ( ξ ) ) p d q U ( ξ ) d ξ q s = N p + s ( N + q ) ,
where N is the degree of U ( ξ ) , i.e., Deg [ U ( ξ ) ] = N .
Step 3: Substituting Equation (6) along with Equations (7) and (8) into Equation (5), we get polynomials in ( d + ϕ ( ξ ) ) k , k = 0 , ± 1 , ± 2 , , ± m , where m is some positive integer. Equating all coefficients of the same power of the resulting polynomials to zero, we obtain an over-determined set of algebraic equations for a j ( j = 0 , ± 1 , ± 2 , , ± N ) , d and V .
Step 4: Assuming that the unknown constants of the algebraic equations in Step 3 can be possibly obtained, we then substitute the values of these constants together with the solution ϕ ( ξ ) of Equation (9) into Equation (6) to get exact traveling wave solutions of (3) when ξ is defined in (4).

2.2.2. Description of the Generalized Kudryashov Method

Step 1: Assume that the solution of Equation (5) can be expressed in a rational form as
U ( ξ ) = i = 0 N a i Q i ( ξ ) j = 0 M b j Q j ( ξ ) ,
where a i ( i = 0 , 1 , 2 , N ) , b j ( j = 0 , 1 , 2 , M ) are constants to be determined at a later step such that a N 0 , b M 0 and the function Q = Q ( ξ ) is a solution of
Q ( ξ ) = Q 2 ( ξ ) Q ( ξ ) .
It is obviously found that the solution of Equation (12) is
Q ( ξ ) = 1 1 + C e ξ ,
where C is a constant of integration.
Step 2: We find the positive integers N and M in Equation (11) by employing the homogeneous balance method, i.e, equating between the highest order derivative and the highest power nonlinear term in Equation (5). The formulas in Equation (10) can also be used in this step.
Step 3: Substituting Equation (11) into Equation (5) along with Equation (12), we obtain a polynomial R ( Q ) of Q. Next, equating all of the coefficients of R ( Q ) to zero, a system of algebraic equations is obtained. Solving this system with the aid of symbolic software packages such as Maple, we can get the values of a i ( i = 0 , 1 , 2 , N ) , b j ( j = 0 , 1 , 2 , M ) . When we substitute these values and the function Q ( ξ ) in Equation (13) into Equation (11), then we attain the exact solutions of the reduced Equation (5). In consequence, exact solutions of (3) are finally obtained using ξ defined in (4).

3. Application of the Methods

The Sharma–Tasso–Olver equation of integer order is expressed as
u t + 3 ρ u x 2 + 3 ρ u 2 u x + 3 ρ u u x x + ρ u x x x = 0 ,
where ρ is an arbitrary real parameter. In this paper, we, however, consider the conformable space-time Sharma–Tasso–Olver equation described as
D t α u + 3 ρ D x β u 2 + 3 ρ u 2 D x β u + 3 ρ u D x 2 β u + ρ D x 3 β u = 0 , t > 0 , 0 < α , β 1 ,
where D s κ u represents the conformable derivative of the function u with respect to an independent variable s of order κ . This function u ( x , t ) is unknown and depends on the variables x and t, and ρ is an arbitrary real parameter. We will apply the novel G G -expansion method and the generalized Kudryashov method to produce exact traveling wave solutions of (15). However, the common step of both methods is to convert (15) into an ODE via the traveling wave transformation in the variables x and t as follows,
ξ = x β β V t α α .
Utilizing Theorem 2, we obtain D t α u = V U , D x β u = U , and D x 2 β u = D x β D x β u , D x 3 β u = D x β D x 2 β u . Then, Equation (15) is reduced into the ODE in the variable U = U ( ξ ) as
V U + ρ U + 3 ρ U 2 + 3 ρ U U + 3 ρ U 2 U = 0 ,
where the prime notation (′) denotes the derivative with respect to ξ . Integrating Equation (17) with respect to ξ along with algebraically manipulating some terms and then letting the constant of integration to be zero, we get the following ODE,
V U + ρ U + 3 ρ U U + ρ U 3 = 0 .

3.1. Obtaining Exact Solutions of Equation (15) Using the Novel ( G G ) -Expansion Method

Utilizing the formulas (10) to balance the highest order derivative U with the nonlinear term of the highest order U 3 in Equation (18), we obtain N = 1 . By Equation (6) in Section 2.2.1, the solution of Equation (18) has the following form,
U ( ξ ) = a 1 ψ ( ξ ) 1 + a 0 + a 1 ψ ( ξ ) ,
where ψ ( ξ ) = d + ϕ ( ξ ) and ϕ ( ξ ) = G G is a solution of the generalized Riccati equation in Equation (9). Substituting Equation (19) into Equation (18), the left hand side of Equation (18) is converted into polynomials of ψ ( ξ ) j = ( d + ϕ ( ξ ) ) j , where j = 0 , ± 1 , ± 2 , ± 3 .
Equating the coefficients of the same power of the resulting polynomials to zero, we have the following set of nonlinear algebraic equations,
( ψ ( ξ ) ) 3 : 2 ρ v 2 a 1 + 3 ρ v a 1 2 + ρ a 1 3 4 ρ v a 1 3 ρ a 1 2 + 2 ρ a 1 = 0 , ( ψ ( ξ ) ) 2 : 6 d ρ v 2 a 1 6 d ρ v a 1 2 + 12 d ρ v a 1 + 6 d ρ a 1 2 + 3 λ ρ v a 1 + 3 λ ρ a 1 2 + 3 ρ v a 0 a 1 + 3 ρ a 0 a 1 2 6 d ρ a 1 3 λ ρ a 1 3 ρ a 0 a 1 = 0 , ψ ( ξ ) : 6 d 2 ρ v 2 a 1 + 3 d 2 ρ v a 1 2 12 d 2 ρ v a 1 3 d 2 ρ a 1 2 6 d λ ρ v a 1 3 d λ ρ a 1 2 6 d ρ v a 0 a 1 + 6 d 2 ρ a 1 + 6 d λ ρ a 1 + 6 d ρ a 0 a 1 + λ 2 ρ a 1 + 3 λ ρ a 0 a 1 + 2 μ ρ v a 1 + 3 μ ρ a 1 2 + 3 ρ a 1 a 1 2 + 3 ρ a 0 2 a 1 2 μ ρ a 1 V a 1 = 0 ,
( ψ ( ξ ) ) 0 : 2 d 3 ρ v 2 a 1 + 4 d 3 ρ v a 1 + 3 d 2 λ ρ v a 1 + 3 d 2 ρ v a 0 a 1 2 d 3 ρ a 1 3 d 2 λ ρ a 1 3 d 2 ρ a 0 a 1 d λ 2 ρ a 1 3 d λ ρ a 0 a 1 2 d μ ρ v a 1 2 d ρ v 2 a 1 + 2 d μ ρ a 1 + 4 d ρ v a 1 + λ μ ρ a 1 + λ ρ v a 1 + 3 μ ρ a 0 a 1 3 ρ v a 1 a 0 + 6 ρ a 1 a 0 a 1 + ρ a 0 3 2 d ρ a 1 λ ρ a 1 + 3 ρ a 1 a 0 V a 0 = 0 , ( ψ ( ξ ) ) 1 : 6 d 2 ρ v 2 a 1 12 d 2 ρ v a 1 6 d λ ρ v a 1 + 6 d ρ v a 1 a 0 + 6 d 2 ρ a 1 + 6 d λ ρ a 1 6 d ρ a 1 a 0 + λ 2 ρ a 1 3 λ ρ a 1 a 0 + 2 μ ρ v a 1 3 ρ v a 1 2 + 3 ρ a 1 2 a 1 + 3 ρ a 1 a 0 2 2 μ ρ a 1 + 3 ρ a 1 2 V a 1 = 0 , ( ψ ( ξ ) ) 2 : 6 d 3 ρ v 2 a 1 + 12 d 3 ρ v a 1 + 9 d 2 λ ρ v a 1 3 d 2 ρ v a 1 a 0 6 d 3 ρ a 1 9 d 2 λ ρ a 1 + 3 d 2 ρ a 1 a 0 3 d λ 2 ρ a 1 + 3 d λ ρ a 1 a 0 6 d μ ρ v a 1 + 6 d ρ v a 1 2 + 6 d μ ρ a 1 6 d ρ a 1 2 + 3 λ μ ρ a 1 3 λ ρ a 1 2 3 μ ρ a 1 a 0 + 3 ρ a 1 2 a 0 = 0 , ( ψ ( ξ ) ) 3 : 2 d 4 ρ v 2 a 1 4 d 4 ρ v a 1 4 d 3 λ ρ v a 1 + 2 d 4 ρ a 1 + 4 d 3 λ ρ a 1 + 2 d 2 λ 2 ρ a 1 + 4 d 2 μ ρ v a 1 3 d 2 ρ v a 1 2 4 d 2 μ ρ a 1 + 3 d 2 ρ a 1 2 4 d λ μ ρ a 1 + 3 d λ ρ a 1 2 + 2 μ 2 ρ a 1 3 μ ρ a 1 2 + ρ a 1 3 = 0 .
Using the symbolic computation software Maple 17 to solve the above system (20), we obtain three independent cases of the unknown constants a 0 , a 1 , a 1 , d , and V. For the sake of convenience, we set
Δ = λ 2 4 μ ( v 1 ) .
Consequently, the exact solutions of Equation (15), depending upon the following three cases of the unknowns in Equation (20) and the families of the solution ϕ ( ξ ) as shown in Appendix of [53], are expressed as follows.
Case 1: The first set of the unknown constants is
a 1 = 0 , a 0 = 2 d v 2 d λ , a 1 = 2 v + 2 , d = d , V = ρ Δ ,
where μ , λ , v , d , and ρ are arbitrary constants. Substituting Equation (22) into Equation (19) and then using Equations (16) and (21), we obtain the following exact solutions in which ξ = x β β ρ Δ t α α .
Family 1: When Δ > 0 and λ ( v 1 ) 0 (or μ ( v 1 ) 0 ), the solutions of Equation (15) written in terms of the hyperbolic functions are as follows,
u 1 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d λ + Δ tanh 1 2 Δ ξ 2 ( v 1 ) ,
u 2 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d λ + Δ coth 1 2 Δ ξ 2 ( v 1 ) ,
u 3 , 4 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d λ + Δ tanh Δ ξ ± i sech Δ ξ 2 ( v 1 ) ,
u 5 , 6 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d λ + Δ coth Δ ξ ± csch Δ ξ 2 ( v 1 ) ,
u 7 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d 2 λ + Δ tanh 1 4 Δ ξ + coth 1 4 Δ ξ 4 ( v 1 ) ,
u 8 , 9 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d + 1 2 ( v 1 ) λ + ± A 2 + B 2 Δ A Δ cosh Δ ξ A sinh Δ ξ + B ,
u 10 , 11 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d + 1 2 ( v 1 ) λ ± A 2 + B 2 Δ + A Δ sinh Δ ξ A cosh Δ ξ + B ,
where A and B are two non-zero real constants that satisfy B 2 A 2 > 0 ,
u 12 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d + 2 μ cosh 1 2 Δ ξ Δ sinh 1 2 Δ ξ λ cosh 1 2 Δ ξ ,
u 13 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d 2 μ sinh 1 2 Δ ξ λ sinh 1 2 Δ ξ Δ cosh 1 2 Δ ξ ,
u 14 , 15 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d + 2 μ cosh Δ ξ Δ sinh Δ ξ λ cosh Δ ξ ± i Δ ,
u 16 , 17 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d + 2 μ sinh Δ ξ λ sinh Δ ξ + Δ cosh Δ ξ ± Δ ,
u 18 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d + 4 μ sinh 1 4 Δ ξ cosh 1 4 Δ ξ 2 λ sinh 1 4 Δ ξ cosh 1 4 Δ ξ + 2 Δ cosh 1 4 Δ ξ 2 Δ .
Family 2: When Δ < 0 and λ ( v 1 ) 0 (or μ ( v 1 ) 0 ), the exact solutions of Equation (15) expressed in terms of the trigonometric functions are as follows,
u 19 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d + λ + Δ tan 1 2 Δ ξ 2 ( v 1 ) ,
u 20 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d λ + Δ cot 1 2 Δ ξ 2 ( v 1 ) ,
u 21 , 22 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d + λ + Δ tan Δ ξ ± sec Δ ξ 2 ( v 1 ) ,
u 23 , 24 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d λ + Δ cot Δ ξ ± csc Δ ξ 2 ( v 1 ) ,
u 25 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d + 2 λ + Δ tan 1 4 Δ ξ cot 1 4 Δ ξ 4 ( v 1 ) ,
u 26 , 27 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d + 1 2 ( v 1 ) λ + ± A 2 B 2 Δ A Δ cos Δ ξ A sin Δ ξ + B ,
u 28 , 29 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d + 1 2 ( v 1 ) λ ± A 2 B 2 Δ + A Δ cos Δ ξ A sin Δ ξ + B ,
where A and B are two non-zero real constants satisfying A 2 B 2 > 0 ,
u 30 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d 2 μ cos 1 2 Δ ξ Δ sin 1 2 Δ ξ + λ cos 1 2 Δ ξ ,
u 31 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d + 2 μ sin 1 2 Δ ξ λ sin 1 2 Δ ξ + Δ cos 1 2 Δ ξ ,
u 32 , 33 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d 2 μ cos Δ ξ Δ sin Δ ξ + λ cos Δ ξ ± Δ ,
u 34 , 35 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d + 2 μ sin Δ ξ λ sin Δ ξ + Δ cos Δ ξ ± Δ ,
u 36 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2
× d + 4 μ sin 1 4 Δ ξ cos 1 / 4 Δ ξ 2 λ sin 1 4 Δ ξ cos 1 4 Δ ξ + 2 Δ cos 1 4 Δ ξ 2 Δ .
Family 3: When μ = 0 and λ ( v 1 ) 0 , the exact solutions of Equation (15) written as the hyperbolic functions are as follows,
u 37 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d λ c 1 v 1 c 1 + cosh λ ξ sinh λ ξ ,
u 38 1 ( x , t ) = 2 d v 2 d λ + 2 v + 2 d λ cosh λ ξ + sinh λ ξ v 1 c 1 + cosh λ ξ + sinh λ ξ ,
where c 1 is an arbitrary constant.
Family 4: When μ = λ = 0 and v 1 0 , the rational function solution of Equation (15) is as follows,
u 39 1 ( x , t ) = 2 d v 2 d + 2 v + 2 d 1 v 1 ξ + c 2 ,
where c 2 is an arbitrary constant.
Case 2: The second set of the unknown constants is
a 1 = λ 2 4 μ v + 4 μ 2 ( v 1 ) , a 0 = 0 , a 1 = 2 v + 2 , d = λ 2 ( v 1 ) , V = 4 ρ Δ ,
where μ , λ , v , and ρ are arbitrary constants. Substituting Equation (51) into Equation (19), and then using Equations (16) and (21), we obtain the following explicit exact solutions in which ξ = x β β 4 ρ Δ t α α .
Family 1: When Δ > 0 and λ ( v 1 ) 0 (or μ ( v 1 ) 0 ), the exact solutions of Equation (15) expressed as the hyperbolic functions are
u 1 2 ( x , t ) = Δ tanh 1 2 Δ ξ + Δ tanh 1 2 Δ ξ ,
u 2 2 ( x , t ) = Δ coth 1 2 Δ ξ + Δ coth 1 2 Δ ξ ,
u 3 , 4 2 ( x , t ) = Δ tanh Δ ξ ± i sech Δ ξ + Δ tanh Δ ξ ± i sech Δ ξ ,
u 5 , 6 2 ( x , t ) = Δ coth Δ ξ ± csch Δ ξ + Δ coth Δ ξ ± csch Δ ξ ,
u 7 2 ( x , t ) = 2 Δ tanh 1 4 Δ ξ + coth 1 4 Δ ξ + Δ tanh 1 4 Δ ξ + coth 1 4 Δ ξ 2 ,
u 8 , 9 2 ( x , t ) = Δ ± A 2 + B 2 Δ A Δ cosh Δ ξ A sinh Δ ξ + B ± A 2 + B 2 Δ A Δ cosh Δ ξ A sinh Δ ξ + B ,
u 10 , 11 2 ( x , t ) = Δ ± A 2 + B 2 Δ + A Δ sinh Δ ξ A cosh Δ ξ + B + ± A 2 + B 2 Δ + A Δ sinh Δ ξ A cosh Δ ξ + B ,
where A and B are non-zero real constants satisfying B 2 A 2 > 0 ,
u 12 2 ( x , t ) = Δ 2 v 1 λ 2 v 2 + 2 μ cosh 1 2 Δ ξ Δ sinh 1 2 Δ ξ λ cosh 1 2 Δ ξ 2 v 1 × λ 2 v 2 + 2 μ cosh 1 2 Δ ξ Δ sinh 1 2 Δ ξ λ cosh 1 2 Δ ξ ,
u 13 2 ( x , t ) = Δ 2 v 1 λ 2 v 2 2 μ sinh 1 2 Δ ξ λ sinh 1 2 Δ ξ Δ cosh 1 2 Δ ξ 2 v 1 × λ 2 v 2 2 μ sinh 1 2 Δ ξ λ sinh 1 2 Δ ξ Δ cosh 1 2 Δ ξ ,
u 14 , 15 2 ( x , t ) = Δ 2 v 1 λ 2 v 2 + 2 μ cosh Δ ξ Δ sinh Δ ξ λ cosh Δ ξ ± i Δ 2 v 1 × λ 2 v 2 + 2 μ cosh Δ ξ Δ sinh Δ ξ λ cosh Δ ξ ± i Δ ,
u 16 , 17 2 ( x , t ) = Δ 2 v 1 λ 2 v 2 + 2 μ sinh Δ ξ λ sinh Δ ξ + Δ cosh Δ ξ ± Δ 2 v 1 × λ 2 v 2 + 2 μ sinh Δ ξ λ sinh Δ ξ + Δ cosh Δ ξ ± Δ ,
u 18 2 ( x , t ) = Δ 2 v 1 λ 2 v 2 + 4 μ sinh 1 4 Δ ξ cosh 1 4 Δ ξ 2 λ sinh 1 4 Δ ξ cosh 1 4 Δ ξ + 2 Δ cosh 1 4 Δ ξ 2 Δ 2 v 1 × λ 2 v 2 + 4 μ sinh 1 4 Δ ξ cosh 1 4 Δ ξ 2 λ sinh 1 4 Δ ξ cosh 1 4 Δ ξ + 2 Δ cosh 1 4 Δ ξ 2 Δ .
Family 2: When Δ < 0 and λ ( v 1 ) 0 (or μ ( v 1 ) 0 ), the exact solutions of Equation (15) written in terms of the trigonometric functions are shown as follows,
u 19 2 ( x , t ) = Δ Δ tan 1 2 Δ ξ Δ tan 1 2 Δ ξ ,
u 20 2 ( x , t ) = Δ Δ cot 1 2 Δ ξ + Δ cot 1 2 Δ ξ ,
u 21 , 22 2 ( x , t ) = Δ Δ tan Δ ξ ± sec Δ ξ Δ tan Δ ξ ± sec Δ ξ ,
u 23 , 24 2 ( x , t ) = Δ Δ cot Δ ξ ± csc Δ ξ + Δ cot Δ ξ ± csc Δ ξ ,
u 25 2 ( x , t ) = 2 Δ Δ tan 1 4 Δ ξ cot 1 4 Δ ξ Δ tan 1 4 Δ ξ cot 1 4 Δ ξ 2 ,
u 26 , 27 2 ( x , t ) = Δ ± A 2 B 2 Δ A Δ cos Δ ξ A sin Δ ξ + B ± A 2 B 2 Δ A Δ cos Δ ξ A sin Δ ξ + B ,
u 28 , 29 2 ( x , t ) = Δ ± A 2 B 2 Δ + A Δ cos Δ ξ A sin Δ ξ + B + ± A 2 B 2 Δ + A Δ cos Δ ξ A sin Δ ξ + B ,
where A and B are two non-zero real constants such that A 2 B 2 > 0 ,
u 30 2 ( x , t ) = Δ 2 v 1 λ 2 v 2 2 μ cos 1 2 Δ ξ Δ sin 1 2 Δ ξ + λ cos 1 2 Δ ξ 2 v 1 × λ 2 v 2 2 μ cos 1 2 Δ ξ Δ sin 1 2 Δ ξ + λ cos 1 2 Δ ξ ,
u 31 2 ( x , t ) = Δ 2 v 1 λ 2 v 2 + 2 μ sin 1 2 Δ ξ λ sin 1 2 Δ ξ + Δ cos 1 2 Δ ξ 2 v 1 × λ 2 v 2 + 2 μ sin 1 2 Δ ξ λ sin 1 2 Δ ξ + Δ cos 1 2 Δ ξ ,
u 32 , 33 2 ( x , t ) = Δ 2 v 1 λ 2 v 2 2 μ cos Δ ξ Δ sin Δ ξ + λ cos Δ ξ ± Δ 2 v 1 × λ 2 v 2 2 μ cos Δ ξ Δ sin Δ ξ + λ cos Δ ξ ± Δ ,
u 34 , 35 2 ( x , t ) = Δ 2 v 1 λ 2 v 2 + 2 μ sin Δ ξ λ sin Δ ξ + Δ cos Δ ξ ± Δ 2 v 1 × λ 2 v 2 + 2 μ sin Δ ξ λ sin Δ ξ + Δ cos Δ ξ ± Δ ,
u 36 2 ( x , t ) = Δ 2 v 1 λ 2 v 2 + 4 μ sin 1 4 Δ ξ cos 1 4 Δ ξ 2 λ sin 1 4 Δ ξ cos 1 4 Δ ξ + 2 Δ cos 1 4 Δ ξ 2 Δ 2 v 1 × λ 2 v 2 + 4 μ sin 1 4 Δ ξ cos 1 4 Δ ξ 2 λ sin 1 4 Δ ξ cos 1 4 Δ ξ + 2 Δ cos 1 4 Δ ξ 2 Δ .
Family 3: When μ = 0 and λ ( v 1 ) 0 , Equation (15) has the hyperbolic function solutions as follows,
u 37 2 ( x , t ) = λ 2 2 v 1 λ 2 v 2 λ c 1 v 1 c 1 + cosh λ ξ sinh λ ξ 2 v 1 × λ 2 v 2 λ c 1 v 1 c 1 + cosh λ ξ sinh λ ξ ,
u 38 2 ( x , t ) = λ 2 2 v 1 λ 2 v 2 λ cosh λ ξ + sinh λ ξ c 1 + cosh λ ξ + sinh λ ξ v 1 2 v 1 × λ 2 v 2 λ cosh λ ξ + sinh λ ξ c 1 + cosh λ ξ + sinh λ ξ v 1 ,
where c 1 is an arbitrary constant.
Family 4: When μ = λ = 0 and v 1 0 , Equation (15) has the rational function solution as follows,
u 39 2 ( x , t ) = 2 ξ + c 2 ,
where c 2 is an arbitrary constant.
Case 3: The last set of the unknown constants is expressed as
a 1 = 2 d 2 v 2 d 2 2 d λ + 2 μ , a 0 = 2 d v + 2 d + λ , a 1 = 0 , d = d , V = ρ Δ ,
where μ , λ , v , d , and ρ are arbitrary constants. Substituting Equation (79) into Equation (19), and then using Equations (16) and (21), we obtain the following exact solutions in which ξ = x β β ρ Δ t α α .
Family 1: When Δ > 0 and λ ( v 1 ) 0 (or μ ( v 1 ) 0 ), the exact solutions of Equation (15) written as the hyperbolic function solutions are as follows,
u 1 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d 1 2 ( v 1 ) λ + Δ tanh 1 2 Δ ξ 2 d v + 2 d + λ ,
u 2 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d 1 2 ( v 1 ) λ + Δ coth 1 2 Δ ξ 2 d v + 2 d + λ ,
u 3 , 4 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d 1 2 ( v 1 ) λ + Δ tanh Δ ξ ± i sec h Δ ξ 2 d v + 2 d + λ ,
u 5 , 6 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d 1 2 ( v 1 ) λ + Δ coth Δ ξ ± csc h Δ ξ 2 d v + 2 d + λ ,
u 7 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d 1 4 ( v 1 ) 2 λ + Δ tanh 1 4 Δ ξ + coth 1 4 Δ ξ 2 d v + 2 d + λ ,
u 8 , 9 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d + 1 2 ( v 1 ) λ + ± A 2 + B 2 Δ A Δ cosh Δ ξ A sinh Δ ξ + B 2 d v + 2 d + λ ,
u 10 , 11 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d + 1 2 ( v 1 ) λ ± A 2 + B 2 Δ + A Δ sinh Δ ξ A cosh Δ ξ + B 2 d v + 2 d + λ ,
where A and B are two non-zero real constants satisfying the condition B 2 A 2 > 0 ,
u 12 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d + 2 μ cosh 1 2 Δ ξ Δ sinh 1 2 Δ ξ λ cosh 1 2 Δ ξ 2 d v + 2 d + λ ,
u 13 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d 2 μ sinh 1 2 Δ ξ λ sinh 1 2 Δ ξ Δ cosh 1 2 Δ ξ 2 d v + 2 d + λ ,
u 14 , 15 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d + 2 μ cosh Δ ξ Δ sinh Δ ξ λ cosh Δ ξ ± i Δ 2 d v + 2 d + λ ,
u 16 , 17 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d + 2 μ sinh Δ ξ λ sinh Δ ξ + Δ cosh Δ ξ ± Δ 2 d v + 2 d + λ ,
u 18 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d + 4 μ sinh 1 4 Δ ξ cosh 1 4 Δ ξ 2 λ sinh 1 4 Δ ξ cosh 1 4 Δ ξ + 2 Δ cosh 1 4 Δ ξ 2 Δ 2 d v + 2 d + λ .
u 19 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d + 1 2 ( v 1 ) λ + Δ tan 1 2 Δ ξ 2 d v + 2 d + λ ,
u 20 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d 1 2 ( v 1 ) λ + Δ cot 1 2 Δ ξ 2 d v + 2 d + λ ,
u 21 , 22 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d + 1 2 ( v 1 ) λ + Δ tan Δ ξ ± sec Δ ξ 2 d v + 2 d + λ ,
u 23 , 24 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d 1 2 ( v 1 ) λ + Δ cot Δ ξ ± csc Δ ξ 2 d v + 2 d + λ ,
u 25 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d + 1 4 ( v 1 ) 2 λ + Δ tan 1 4 Δ ξ cot 1 4 Δ ξ 2 d v + 2 d + λ ,
u 26 , 27 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d + 1 2 ( v 1 ) λ + ± ( A 2 B 2 ) Δ A Δ cos Δ ξ A sin Δ ξ + B 2 d v + 2 d + λ ,
u 28 , 29 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d + 1 2 ( v 1 ) λ ± ( A 2 B 2 ) Δ + A Δ cos Δ ξ A sin Δ ξ + B 2 d v + 2 d + λ ,
where A and B are two non-zero real constants satisfying the condition A 2 B 2 > 0 ,
u 30 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d 2 μ cos 1 2 Δ ξ Δ sin 1 2 Δ ξ + λ cos 1 2 Δ ξ 2 d v + 2 d + λ ,
u 31 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d + 2 μ sin 1 2 Δ ξ λ sin 1 2 Δ ξ + Δ cos 1 2 Δ ξ 2 d v + 2 d + λ ,
u 32 , 33 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d 2 μ cos ( Δ ξ ) Δ sin ( Δ ξ ) + λ cos ( Δ ξ ) ± Δ 2 d v + 2 d + λ ,
u 34 , 35 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d + 2 μ sin ( Δ ξ ) λ sin ( Δ ξ ) + Δ cos ( Δ ξ ) ± Δ 2 d v + 2 d + λ ,
u 36 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ + 2 μ d + 4 μ sin 1 4 Δ ξ cos 1 4 Δ ξ 2 λ sin 1 4 Δ ξ cos 1 4 Δ ξ + 2 Δ cos 1 4 Δ ξ 2 Δ 2 d v + 2 d + λ .
Family 3: When μ = 0 and λ ( v 1 ) 0 , Equation (15) has the hyperbolic function solutions as follows,
u 37 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ d λ c 1 ( v 1 ) c 1 + cosh ( λ ξ ) sinh ( λ ξ ) 2 d v + 2 d + λ ,
u 38 3 ( x , t ) = 2 d 2 v 2 d 2 2 d λ d λ ( cosh ( λ ξ ) + sinh ( λ ξ ) ) ( v 1 ) c 1 + cosh ( λ ξ ) + sinh ( λ ξ ) 2 d v + 2 d + λ ,
where c 1 is an arbitrary constant.
Family 4: When μ = λ = 0 and v 1 0 , the exact solution of Equation (15) expressed as the rational function is
u 39 3 ( x , t ) = 2 d 2 v 2 d 2 d 1 ( v 1 ) ξ + c 2 2 d v + 2 d ,
where c 2 is an arbitrary constant.

3.2. Obtaining Exact Solutions of Equation (15) Using the Generalized Kudryashov Method

Substituting solution form (11) into (18) and then applying the homogeneous balance principle to the resulting equation, we have
N M + 2 = 3 N 3 M .
If we choose M = 1 , then N = 2 . Using Equation (11) in Section 2.2.2, the exact solution of Equation (18) takes the form
U ( ξ ) = a 0 + a 1 Q + a 2 Q 2 b 0 + b 1 Q ,
where Q = Q ( ξ ) satisfies Equation (12). The parameters a 0 , a 1 , a 2 , b 0 , and b 1 are determined at the next step. Substituting Equation (108) into Equation (18), and utilizing Equation (12) and then setting all coefficients of the functions Q k to zero, we get
Q 0 : V a 0 b 0 2 + ρ a 0 3 = 0 , Q 1 : 2 V a 0 b 0 b 1 V a 1 b 0 2 + 3 ρ a 0 2 a 1 + 3 ρ a 0 2 b 1 3 ρ a 0 a 1 b 0 ρ a 0 b 0 b 1 + ρ a 1 b 0 2 = 0 , Q 2 : V a 0 b 1 2 2 V a 1 b 0 b 1 V a 2 b 0 2 + 3 ρ a 0 2 a 2 3 ρ a 0 2 b 1 + 3 ρ a 0 a 1 2 + 3 ρ a 0 a 1 b 0 + 3 ρ a 0 a 1 b 1 6 ρ a 0 a 2 b 0 + 3 ρ a 0 b 0 b 1 + ρ a 0 b 1 2 3 ρ a 1 2 b 0 3 ρ a 1 b 0 2 ρ a 1 b 0 b 1 + 4 ρ a 2 b 0 2 = 0 , Q 3 : V a 1 b 1 2 2 V a 2 b 0 b 1 + 6 ρ a 0 a 1 a 2 3 ρ a 0 a 1 b 1 + 6 ρ a 0 a 2 b 0 2 ρ a 0 b 0 b 1 ρ a 0 b 1 2 + ρ a 1 3 + 3 ρ a 1 2 b 0 9 ρ a 1 a 2 b 0 + 2 ρ a 1 b 0 2 + ρ a 1 b 0 b 1 10 ρ a 2 b 0 2 + 3 ρ a 2 b 0 b 1 = 0 , Q 4 : V a 2 b 1 2 + 3 ρ a 0 a 2 2 + 3 ρ a 1 2 a 2 + 9 ρ a 1 a 2 b 0 3 ρ a 1 a 2 b 1 6 ρ a 2 2 b 0 + 6 ρ a 2 b 0 2 9 ρ a 2 b 0 b 1 + ρ a 2 b 1 2 = 0 , Q 5 : 3 ρ a 1 a 2 2 + 3 ρ a 1 a 2 b 1 + 6 ρ a 2 2 b 0 3 ρ a 2 2 b 1 + 6 ρ a 2 b 0 b 1 3 ρ a 2 b 1 2 = 0 , Q 6 : ρ a 2 3 + 3 ρ a 2 2 b 1 + 2 ρ a 2 b 1 2 = 0 .
Solving the above algebraic system with the aid of the Maple package program, we get the following results.
Case 1:
a 0 = 0 , a 1 = 0 , a 2 = b 1 , b 0 = b 1 2 , b 1 = b 1 , V = 4 ρ ,
where b 1 is an arbitrary constant. From Equations (13), (108), and (110), we obtain the simplified exact solution of Equation (15) as follows,
u ( x , t ) = 2 C 2 e 2 ξ 1 ,
where ξ = x β β 4 ρ t α α .
Case 2:
a 0 = b 0 , a 1 = 2 b 0 , a 2 = b 1 , b 0 = b 0 , b 1 = b 1 , V = ρ ,
where b 0 , b 1 are arbitrary constants. From Equations (13), (108), and (112), the simplified exact solution of Equation (15) can be obtained as
u ( x , t ) = b 0 C 2 e 2 ξ b 0 b 1 1 + C e ξ b 0 C e ξ + b 0 + b 1 ,
where ξ = x β β ρ t α α .
Case 3:
a 0 = b 0 2 , a 1 = b 0 + b 1 2 , a 2 = b 1 , b 0 = b 0 , b 1 = b 1 , V = ρ 4 ,
where b 0 , b 1 are arbitrary constants. From Equations (13), (108), and (114), the simplified exact solution of Equation (15) can be expressed as
u 6 ( x , t ) = Ce ξ 1 2 C e ξ + 2 .
where ξ = x β β ρ t α 4 α .
Case 4:
a 0 = b 0 , a 1 = 2 b 0 + b 1 , a 2 = 2 b 1 , b 0 = b 0 , b 1 = b 1 , V = ρ ,
where b 0 , b 1 are arbitrary constants. From Equations (13), (108), and (116), we obtain the simplified exact solution of Equation (15) as follows,
u 3 ( x , t ) = Ce ξ 1 1 + Ce ξ ,
where ξ = x β β ρ t α α .
Case 5:
a 0 = b 1 , a 1 = 2 b 1 , a 2 = 2 b 1 , b 0 = b 1 2 , b 1 = b 1 , V = 4 ρ ,
where b 1 is an arbitrary constant. From Equations (13), (108), and (118), the simplified exact solution of Equation (15) can be obtained as follows,
u 4 ( x , t ) = 2 C 2 e 2 ξ + 2 C 2 e 2 ξ 1
where ξ = x β β 4 ρ t α α .
Case 6:
a 0 = 2 b 0 , a 1 = 4 b 0 , a 2 = 2 b 0 , b 0 = b 0 , b 1 = 2 b 0 , V = 4 ρ ,
where b 0 is an arbitrary constant. From Equations (13), (108), and (120), we get the simplified exact solution of Equation (15) as follows,
u 8 ( x , t ) = 2 C 2 e 2 ξ C 2 e 2 ξ 1 ,
where ξ = x β β 4 ρ t α α .

4. Graphical Representations of Some Exact Solutions and Their Physical Explanations

In this section, we will give some graphical representations of the above-determined exact solutions of the conformable space-time Sharma–Tasso–Olver Equation (15), which were obtained using the novel ( G G ) -expansion method and the generalized Kudryashov method in Section 3. Here, we take ρ = 2 , and use the following sets of the fractional orders, { β = 1 , α = 1 } , { β = 0.8 , α = 0.2 } , and { β = 0.5 , α = 0.5 } , for the equation. Some selected explicit exact solutions will be plotted as two- and three-dimensional graphs in which the used domain is 0 x 100 and 0 t 10 . Their corresponding contours are drawn as well. Furthermore, their physical explanations of the solutions are included.
Figure 1, Figure 2 and Figure 3 demonstrate the graphical representations of some chosen exact solutions of the problem obtained using the novel ( G G ) -expansion method. They are described below. In Figure 1, we show different plots of the exact solution u 9 1 ( x , t ) in Equation (28) using the following parameter values, μ = 0.5 , λ = 1 , v = 0.5 , d = 1 , A = 0.5 , and B = 1 . In particular, Figure 1a–c shows the 3-D plot, the 2-D plot with t fixed at t = 1 , and the contour plot of solution (28), respectively, when the set of the fractional orders { β = 1 , α = 1 } is used. Using the same parameter values as shown above except using { β = 0.8 , α = 0.2 } , the 3-D graph, the 2-D graph while t is held fixed at t = 1 and the contour graph of solution (28) are plotted in Figure 1d–f, respectively. Proceeding in a manner analogous to the above plots except using the fractional order set { β = 0.5 , α = 0.5 } , we obtain the 3-D plot, the 2-D plot with t = 1 , and the contour plot of solution (28) in Figure 1g–i, respectively. By characterizing the shapes of the plots in Figure 1, solution (28) behaves as a singular kink-type solution.
Distinct plots of the exact solution u 27 2 ( x , t ) in Equation (69) are portrayed in Figure 2 using the following parameter values, μ = 1 , λ = 0.5 , v = 1.5 , A = 1 , and B = 0.5 . Specifically, Figure 2a–c describes the 3-D plot, the 2-D plot as t is fixed at t = 1 , and the contour plot of solution (69), respectively, when the set of the fractional orders { β = 1 , α = 1 } is used. Utilizing the same parameter values as shown above, except using { β = 0.8 , α = 0.2 } , the 3-D graph, the 2-D graph when t is held fixed at t = 1 , and the contour graph of solution (69) are plotted in Figure 2d–f, respectively. Proceeding in a like manner to the above plots, except using the fractional order set { β = 0.5 , α = 0.5 } , we have the 3-D plot, the 2-D plot with t = 1 , and the contour graph of solution (69) in Figure 2g–i, respectively. By characterizing the shapes of the plots in Figure 2, solution (69) displays the behavior of a singularly periodic traveling wave solution.
The selected exact solution u 38 3 ( x , t ) in Equation (105) is plotted in Figure 3 using the three sets of the fractional orders and the following parameter values, λ = 0.5 , v = 1.6 , d = c 1 = 1 . In particular, Figure 3a–c shows the 3-D plot, the 2-D plot when t = 1 , and the contour plot of solution (105), respectively, when we use { β = 1 , α = 1 } . Employing the same parameter values as shown above except using { β = 0.8 , α = 0.2 } , the 3-D graph, the 2-D graph as t is held fixed at t = 1 , and the contour graph of solution (105) are plotted in Figure 3d–f, respectively. Proceeding in a manner analogous to the above plots, except using the fractional order set { β = 0.5 , α = 0.5 } , we obtain the 3-D plot, the 2-D plot with t = 1 , and the contour graph of solution (105) in Figure 3g–i, respectively. From these plots, solution (105) is characterized as a kink-type solution.
Figure 4 and Figure 5 display graphical representations of some chosen exact solutions of Equation (15) obtained using the generalized Kudryashov method. The following details graphically explain the obtained outcomes. Distinct plots of the exact solution (113) are depicted in Figure 4 using the following parameter values, b 0 = b 1 = C = 1 . Particularly, Figure 4a–c describes the 3-D graph, the 2-D graph with t = 2 , and the contour graph of solution (113), respectively, when the set of the fractional orders { β = 1 , α = 1 } is employed. Using the same parameter values as mentioned above, except utilizing { β = 0.8 , α = 0.2 } , the 3-D plot, the 2-D plot as t is held fixed at t = 2 , and the contour plot of solution (113) are plotted in Figure 4d–f, respectively. Proceeding in like manner to the above plots, except using the fractional-order set { β = 0.5 , α = 0.5 } , we have the 3-D plot, the 2-D plot with t = 2 , and the contour graph of solution (113) in Figure 4g–i, respectively. By characterizing the shapes of the plots in Figure 4, solution (113) shows the behavior of a kink-type solution.
Using the second method, the exact solution u ( x , t ) in Equation (119) is graphically shown in Figure 5 with the three different sets of the fractional orders and the constant C = 1 . In particular, Figure 5a–c displays the 3-D plot, the 2-D plot with t = 2 , and the contour plot of solution (119), respectively, when { β = 1 , α = 1 } is used. Employing the same parameter values above, except using { β = 0.8 , α = 0.2 } , the 3-D graph, the 2-D graph with t = 2 , and the contour graph of solution (119) are plotted in Figure 5d–f, respectively. Proceeding in a manner analogous to the previous plots, except using the fractional-order set { β = 0.5 , α = 0.5 } , we obtain the 3-D plot, the 2-D plot with t = 2 , and the contour graph of solution (119) in Figure 5g–i, respectively. By observing the shapes of the plots in Figure 5, solution (119) is classified into a singular multiple-soliton solution.

5. Conclusions

In this article, we have constructed explicit exact solutions of the (1+1)-dimensional conformable space-time Sharma–Tasso–Olver equation expressed in Equation (15) using the novel G G -expansion method and the generalized Kudryashov method with the help of the fractional complex transform and the symbolic computation package Maple 17. The obtained results have revealed that the methods are straightforward, reliable, and powerful. In particular, the novel G G -expansion method gives hyperbolic, trigonometric, and rational function solutions for the equation; however, the generalized Kudryashov method provides fractional solutions of the exponential functions, which are possibly converted into hyperbolic function solutions. Some of these obtained exact solutions have been graphically characterized into a variety of various physical structures such as the single-kink wave solution, the singular periodic wave solution, and the singular multiple-soliton solution. The applications of these exact solutions have been discovered in several physical phenomena such as plasma waves and optical fibers. All calculations in this investigation have been made and verified using the Maple package program. Previously, many authors had tackled the fractional Sharma–Tasso–Olver equation in different approaches. For example, finding exact solutions of the fractional Sharma–Tasso–Olver equation in the sense of the modified Riemann-Liouville derivative using the ( G / G , 1 / G ) -expansion method [56], the tanh ansatz method [57], the exp ( Φ ( ξ ) ) method [58], the sub-equation method [59], and the improved extended tanh-coth method [13]. In addition, constructing exact solutions of the nonlinear conformable time Sharma–Tasso–Olver equation via conformable derivatives was done using the simplest equation method [60] and the direct algebraic method [61]. The Sharma–Tasso–Olver equation as shown in the above articles involves with both of only the conformable time derivative and the conformable space-time derivatives. However, in this paper, we construct exact solutions of Equation (15) in the sense of conformable derivatives with respect to x and t. To the best of our knowledge, in [62], the authors analytically solved (15) using the Exp-function method and only two exact solutions, which are expressed in terms of the exponential function solutions and characterized as kink-type solutions, were constructed. Comparing our results obtained using the novel G G -expansion method with the results in [62], not only the number of our exact solutions is considerably more than the number of their results, but our solutions are also classified into more different types. This is because the novel G G -expansion method is generalized from many similar methods such as the G G -expansion method, improved G G -expansion method, and the generalized and improved G G -expansion method. In consequence, some of our exact solutions of (15) via the two methods are novel and reported here for the first time. According to the mentioned advantages of the methods for obtaining exact traveling wave solutions, they could be applied efficiently for a wide range of nonlinear conformable partial differential equations or other fractional-order PDEs, which appear in several branches of the applied sciences and engineering.

Author Contributions

All authors contributed equally to the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by King Mongkut’s University of Technology North Bangkok under contract no. KMUTNB-61-KNOW-037.

Acknowledgments

The authors are grateful to anonymous referees for the valuable comments which have significantly improved this paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Hashim, I. Adomian decomposition method for solving BVPs for fourth-order integro-differential equations. J. Comput. Appl. Math. 2006, 193, 658–664. [Google Scholar] [CrossRef] [Green Version]
  2. Basak, K.C.; Ray, P.C.; Bera, R.K. Solution of non-linear Klein-Gordon equation with a quadratic non-linear term by Adomian decomposition method. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 718–723. [Google Scholar] [CrossRef]
  3. Gorder, R.A.V. The variational iteration method is a special case of the homotopy analysis method. Appl. Math. Lett. 2015, 45, 81–85. [Google Scholar] [CrossRef]
  4. Wazwaz, A. The variational iteration method: A reliable analytic tool for solving linear and nonlinear wave equations. Comput. Math. Appl. 2007, 54, 926–932. [Google Scholar] [CrossRef] [Green Version]
  5. Nazari, D.; Shahmorad, S. Application of the fractional differential transform method to fractional-order integro-differential equations with nonlocal boundary conditions. J. Comput. Appl. Math. 2010, 234, 883–891. [Google Scholar] [CrossRef] [Green Version]
  6. Bera, P.; Sil, T. Homotopy perturbation method in quantum mechanical problems. Appl. Math. Comput. 2012, 219, 3272–3278. [Google Scholar] [CrossRef]
  7. Nik, H.S.; Effati, S.; Shirazian, M. An approximate-analytical solution for the Hamilton-Jacobi-Bellman equation via homotopy perturbation method. Appl. Math. Model. 2012, 36, 5614–5623. [Google Scholar]
  8. Qin, Z.; Mu, G.; Ma, H. The ( G G )-expansion method for the fifth-order forms of KdV-Sawada-Kotera equation. Appl. Math. Comput. 2013, 222, 29–33. [Google Scholar] [CrossRef]
  9. Aslan, İ. Exact solutions for fractional DDEs via auxiliary equation method coupled with the fractional complex transform. Math. Methods Appl. Sci. 2016, 39, 5619–5625. [Google Scholar] [CrossRef] [Green Version]
  10. Guner, O.; Bekir, A.; Ünsal, Ö. Two reliable methods for solving the time fractional Clannish Random Walkers Parabolic equation. Opt. Int. J. Light Electron Opt. 2016, 127, 9571–9577. [Google Scholar] [CrossRef]
  11. Abdel-Salam, E.A.B.; Gumma, E.A. Analytical solution of nonlinear space-time fractional differential equations using the improved fractional Riccati expansion method. Ain Shams Eng. J. 2015, 6, 613–620. [Google Scholar] [CrossRef] [Green Version]
  12. Abdel-Salam, E.A.B.; Yousif, E.A.; El-Aasser, M.A. Analytical solution of the space-time fractional nonlinear Schrödinger equation. Rep. Math. Phys. 2016, 77, 19–34. [Google Scholar] [CrossRef]
  13. Gómez, S.; Cesar, A. A nonlinear fractional Sharma–Tasso–Olver equation: New exact solutions. Appl. Math. Comput. 2015, 266, 385–389. [Google Scholar] [CrossRef]
  14. Hubert, M.B.; Betchewe, G.; Doka, S.Y.; Crepin, K.T. Soliton wave solutions for the nonlinear transmission line using the Kudryashov method and the ( G G )-expansion method. Appl. Math. Comput. 2014, 239, 299–309. [Google Scholar] [CrossRef]
  15. Ryabov, P.N.; Sinelshchikov, D.I.; Kochanov, M.B. Application of the Kudryashov method for finding exact solutions of the high order nonlinear evolution equations. Appl. Math. Comput. 2011, 218, 3965–3972. [Google Scholar] [CrossRef] [Green Version]
  16. Mohyud-Din, S.T.; Nawaz, T.; Azhar, E.; Akbar, M.A. Fractional sub-equation method to space-time fractional Calogero-Degasperis and potential Kadomtsev-Petviashvili equations. J. Taibah Univ. Sci. 2017, 11, 258–263. [Google Scholar] [CrossRef] [Green Version]
  17. Abdel-Salam, E.A.; Hassan, G.F. Solutions to class of linear and nonlinear fractional differential equations. Commun. Theor. Phys. 2016, 65, 127. [Google Scholar] [CrossRef]
  18. Wang, M.; Li, X.; Zhang, J. The ( G G )-expansion method and travelling wave solutions of nonlinear evolution equations in mathematical physics. Phys. Lett. A 2008, 372, 417–423. [Google Scholar] [CrossRef]
  19. Zhang, J.; Jiang, F.; Zhao, X. An improved ( G G )-expansion method for solving nonlinear evolution equations. Int. J. Comput. Math. 2010, 87, 1716–1725. [Google Scholar] [CrossRef]
  20. Akbar, M.A.; Ali, N.H.M.; Zayed, E.M.E. A generalized and improved ( G G )-expansion method for nonlinear evolution equations. Math. Probl. Eng. 2012, 2012, 459879. [Google Scholar] [CrossRef] [Green Version]
  21. Hafez, M.; Alam, M.N.; Akbar, M.A. Exact traveling wave solutions to the Klein-Gordon equation using the novel ( G G )-expansion method. Results Phys. 2014, 4, 177–184. [Google Scholar] [CrossRef] [Green Version]
  22. Alam, M.N.; Akbar, M.A.; Mohyud-Din, S.T. A novel ( G G )-expansion method and its application to the Boussinesq equation. Chin. Phys. B 2014, 23, 020203. [Google Scholar] [CrossRef]
  23. Alam, M.N.; Akbar, M.A. Traveling wave solutions of the nonlinear (1+1)-dimensional modified Benjamin-Bona-Mahony equation by using novel ( G G )-expansion method. Phys. Rev. Res. Int. 2014, 4, 147–165. [Google Scholar]
  24. Hafez, M.G.; Zheng, B.; Akbar, M. Exact travelling wave solutions of the coupled nonlinear evolution equation via the Maccari system using novel ( G G )-expansion method. Egypt. J. Basic Appl. Sci. 2015, 2, 206–220. [Google Scholar] [CrossRef] [Green Version]
  25. Alam, M.N.; Belgacem, F.B.M.; Akbar, M.A. Analytical Treatment of the Evolutionary (1+ 1)-Dimensional Combined KdV-mKdV Equation via the Novel ( G G )-Expansion Method. J. Appl. Math. Phys. 2015, 3, 1571. [Google Scholar] [CrossRef] [Green Version]
  26. Alam, M.N.; Hafez, M.; Belgacem, F.B.M.; Akbar, M.A. Applications of the novel ( G G )-expansion method to find new exact traveling wave solutions of the nonlinear coupled Hsiggs field equation. Nonlinear Stud. 2015, 22, 613–633. [Google Scholar]
  27. Hafez, M. New travelling wave solutions of the (1+1)-dimensional cubic nonlinear Schrödinger equation using novel ( G G )-expansion method. Beni-Suef Univ. J. Basic Appl. Sci. 2016, 5, 109–118. [Google Scholar] [CrossRef] [Green Version]
  28. Alam, M.N.; Belgacem, F.B.M. Exact Traveling Wave Solutions for the (1+ 1)-Dimensional Compound KdVB Equation via the Novel ( G G )-Expansion Method. Int. J. Mod. Nonlinear Theory Appl. 2016, 5, 28. [Google Scholar] [CrossRef] [Green Version]
  29. Akbar, M.A.; Alam, M.N.; Hafez, M.G. Application of the novel ( G G )-expansion method to construct traveling wave solutions to the positive Gardner–KP equation. Indian J. Pure Appl. Math. 2016, 47, 85–96. [Google Scholar] [CrossRef]
  30. Islam, M.S.; Khan, K.; Arnous, A.H. Generalized Kudryashov method for solving some (3+ 1)-dimensional nonlinear evolution equations. New Trends Math. Sci. 2015, 3, 46. [Google Scholar]
  31. Kaplan, M.; Bekir, A.; Akbulut, A. A generalized Kudryashov method to some nonlinear evolution equations in mathematical physics. Nonlinear Dyn. 2016, 85, 2843–2850. [Google Scholar] [CrossRef]
  32. Mahmud, F.; Samsuzzoha, M.; Akbar, M.A. The generalized Kudryashov method to obtain exact traveling wave solutions of the PHI-four equation and the Fisher equation. Results Phys. 2017, 7, 4296–4302. [Google Scholar] [CrossRef]
  33. Biswas, A.; Sonmezoglu, A.; Ekici, M.; Mirzazadeh, M.; Zhou, Q.; Moshokoa, S.P.; Belic, M. Optical soliton perturbation with fractional temporal evolution by generalized Kudryashov’s method. Optik 2018, 164, 303–310. [Google Scholar] [CrossRef]
  34. Nestor, S.; Justin, M.; Douvagai; Betchewe, G.; Doka, S.Y.; Kofane, T. New Jacobi elliptic solutions and other solutions with quadratic–Cubic nonlinearity using two mathematical methods. Asian-Eur. J. Math. 2018, 13, 2050043. [Google Scholar] [CrossRef]
  35. Houwe, A.; Justin, M.; Doka, S.Y.; Crepin, K.T. New traveling wave solutions of the perturbed nonlinear Schrödingers equation in the left–handed metamaterials. Asian-Eur. J. Math. 2018, 13, 2050022. [Google Scholar] [CrossRef]
  36. Lu, D.; Seadawy, A.R.; Khater, M.M. Structure of solitary wave solutions of the nonlinear complex fractional generalized Zakharov dynamical system. Adv. Differ. Equ. 2018, 2018, 266. [Google Scholar] [CrossRef] [Green Version]
  37. Demiray, S.T.; Bulut, H. Soliton solutions of some nonlinear evolution problems by GKM. Neural Comput. Appl. 2019, 31, 287–294. [Google Scholar] [CrossRef]
  38. Habib, M.; Ali, H.S.; Miah, M.M.; Akbar, M.A. The generalized Kudryashov method for new closed form traveling wave solutions to some NLEEs. AIMS Math. 2019, 4, 896. [Google Scholar] [CrossRef]
  39. Wazwaz, A.M. New solitons and kinks solutions to the Sharma-Tasso-Olver equation. Appl. Math. Comput. 2007, 188, 1205–1213. [Google Scholar] [CrossRef]
  40. Pan, J.T.; Chen, W.Z. A new auxiliary equation method and its application to the Sharma-Tasso-Olver model. Phys. Lett. A 2009, 373, 3118–3121. [Google Scholar]
  41. Bekir, A.; Boz, A. Exact solutions for nonlinear evolution equations using Exp-function method. Phys. Lett. A 2008, 372, 1619–1625. [Google Scholar] [CrossRef]
  42. He, Y.; Li, S.; Long, Y. Exact solutions to the Sharma-Tasso-Olver equation by using improved ( G G )-expansion method. J. Appl. Math. 2013, 2013, 247234. [Google Scholar] [CrossRef] [Green Version]
  43. Uğurlu, Y.; Kaya, D. Analytic method for solitary solutions of some partial differential equations. Phys. Lett. A 2007, 370, 251–259. [Google Scholar] [CrossRef]
  44. Khalil, R.; Horani, M.A.; Yousef, A.; Sababheh, M. A new definition of fractional derivative. J. Comput. Appl. Math. 2014, 264, 65–70. [Google Scholar] [CrossRef]
  45. Kumar, D.; Seadawy, A.R.; Joardar, A.K. Modified Kudryashov method via new exact solutions for some conformable fractional differential equations arising in mathematical biology. Chin. J. Phys. 2018, 56, 75–85. [Google Scholar] [CrossRef]
  46. Rezazadeh, H.; Korkmaz, A.; Eslami, M.; Vahidi, J.; Asghari, R. Traveling wave solution of conformable fractional generalized reaction Duffing model by generalized projective Riccati equation method. Opt. Quantum Electron. 2018, 50, 150. [Google Scholar] [CrossRef]
  47. Korkmaz, A. Explicit exact solutions to some one-dimensional conformable time fractional equations. Waves Random Complex Media 2019, 29, 124–137. [Google Scholar] [CrossRef]
  48. Ekici, M.; Mirzazadeh, M.; Eslami, M.; Zhou, Q.; Moshokoa, S.P.; Biswas, A.; Belic, M. Optical soliton perturbation with fractional-temporal evolution by first integral method with conformable fractional derivatives. Optik 2016, 127, 10659–10669. [Google Scholar] [CrossRef]
  49. Tarasov, V.E. No nonlocality. No fractional derivative. Commun. Nonlinear Sci. Numer. Simul. 2018, 62, 157–163. [Google Scholar] [CrossRef] [Green Version]
  50. He, J.H.; Elagan, S.; Li, Z. Geometrical explanation of the fractional complex transform and derivative chain rule for fractional calculus. Phys. Lett. A 2012, 376, 257–259. [Google Scholar] [CrossRef]
  51. Su, W.H.; Yang, X.J.; Jafari, H.; Baleanu, D. Fractional complex transform method for wave equations on Cantor sets within local fractional differential operator. Adv. Differ. Equ. 2013, 2013, 97. [Google Scholar] [CrossRef] [Green Version]
  52. Güner, O.; Bekir, A.; Cevikel, A.C. A variety of exact solutions for the time fractional Cahn–Allen equation. Eur. Phys. J. Plus 2015, 130, 146. [Google Scholar] [CrossRef]
  53. Sirisubtawee, S.; Koonprasert, S.; Khaopant, C.; Porka, W. Two Reliable Methods for Solving the (3 + 1)-Dimensional Space-Time Fractional Jimbo-Miwa Equation. Math. Probl. Eng. 2017, 2017, 1–30. [Google Scholar] [CrossRef] [Green Version]
  54. Zhu, S.D. The generalizing Riccati equation mapping method in non-linear evolution equation: Application to (2+1)-dimensional Boiti-Leon-Pempinelle equation. Chaos Solitons Fractals 2008, 37, 1335–1342. [Google Scholar]
  55. Zheng, C.L. Comments on “The generalizing Riccati equation mapping method in nonlinear evolution equation: Application to (2 + 1)-dimensional Boiti-Leon-Pempinelle equation”. Chaos Solitons Fractals 2009, 39, 1493–1495. [Google Scholar] [CrossRef]
  56. Uddin, M.H.; Khan, M.A.; Akbar, M.A.; Haque, M.A. Multi-solitary wave solutions to the general time fractional Sharma–Tasso–Olver equation and the time fractional Cahn–Allen equation. Arab. J. Basic Appl. Sci. 2019, 26, 193–201. [Google Scholar] [CrossRef]
  57. Guner, O.; Korkmaz, A.; Bekir, A. Dark soliton solutions of space-time fractional Sharma–Tasso–Olver and potential Kadomtsev–Petviashvili equations. Commun. Theor. Phys. 2017, 67, 182. [Google Scholar] [CrossRef]
  58. Kaplan, M.; Bekir, A. A novel analytical method for time-fractional differential equations. Optik 2016, 127, 8209–8214. [Google Scholar] [CrossRef]
  59. Jafari, H.; Tajadodi, H.; Baleanu, D.; Al-Zahrani, A.; Alhamed, Y.; Zahid, A. Fractional sub-equation method for the fractional generalized reaction Duffing model and nonlinear fractional Sharma-Tasso-Olver equation. Open Phys. 2013, 11, 1482–1486. [Google Scholar] [CrossRef]
  60. Chen, C.; Jiang, Y.L. Simplest equation method for some time-fractional partial differential equations with conformable derivative. Comput. Math. Appl. 2018, 75, 2978–2988. [Google Scholar] [CrossRef]
  61. Rezazadeh, H.; Khodadad, F.S.; Manafian, J. New structure for exact solutions of nonlinear time fractional Sharma–Tasso–Olver equation via conformable fractional derivative. Appl. Appl. Math. Int. J. 2017, 12, 13–21. [Google Scholar]
  62. Yaslan, H.Ç.; Girgin, A. Exp-function method for the conformable space-time fractional STO, ZKBBM and coupled Boussinesq equations. Arab. J. Basic Appl. Sci. 2019, 26, 163–170. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Associated plots of u 9 1 ( x , t ) in Equation (28): (ac) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 1 , α = 1 } ; (df) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.8 , α = 0.2 } ; (gi) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.5 , α = 0.5 } .
Figure 1. Associated plots of u 9 1 ( x , t ) in Equation (28): (ac) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 1 , α = 1 } ; (df) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.8 , α = 0.2 } ; (gi) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.5 , α = 0.5 } .
Symmetry 12 00644 g001
Figure 2. Associated plots of u 27 2 ( x , t ) in Equation (69): (ac) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 1 , α = 1 } ; (df) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.8 , α = 0.2 } ; (gi) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.5 , α = 0.5 } .
Figure 2. Associated plots of u 27 2 ( x , t ) in Equation (69): (ac) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 1 , α = 1 } ; (df) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.8 , α = 0.2 } ; (gi) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.5 , α = 0.5 } .
Symmetry 12 00644 g002
Figure 3. Associated plots of u 38 3 ( x , t ) in Equation (105): (ac) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 1 , α = 1 } ; (df) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.8 , α = 0.2 } ; (gi) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.5 , α = 0.5 } .
Figure 3. Associated plots of u 38 3 ( x , t ) in Equation (105): (ac) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 1 , α = 1 } ; (df) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.8 , α = 0.2 } ; (gi) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.5 , α = 0.5 } .
Symmetry 12 00644 g003
Figure 4. Associated plots of u(x, t) in Equation (113): (ac) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 1 , α = 1 } ; (df) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.8 , α = 0.2 } ; (gi) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.5 , α = 0.5 } .
Figure 4. Associated plots of u(x, t) in Equation (113): (ac) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 1 , α = 1 } ; (df) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.8 , α = 0.2 } ; (gi) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.5 , α = 0.5 } .
Symmetry 12 00644 g004
Figure 5. Associated plots of u(x, t) in Equation (119): (ac) 3-D plot, 2-D plot, and contour plot, respectively, when { { β = 1 , α = 1 } }; (df) 3-D plot, 2-D plot, and contour plot, respectively, when { { β = 0.8 , α = 0.2 } }; (gi) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.5 , α = 0.5 } .
Figure 5. Associated plots of u(x, t) in Equation (119): (ac) 3-D plot, 2-D plot, and contour plot, respectively, when { { β = 1 , α = 1 } }; (df) 3-D plot, 2-D plot, and contour plot, respectively, when { { β = 0.8 , α = 0.2 } }; (gi) 3-D plot, 2-D plot, and contour plot, respectively, when { β = 0.5 , α = 0.5 } .
Symmetry 12 00644 g005

Share and Cite

MDPI and ACS Style

Sirisubtawee, S.; Koonprasert, S.; Sungnul, S. New Exact Solutions of the Conformable Space-Time Sharma–Tasso–Olver Equation Using Two Reliable Methods. Symmetry 2020, 12, 644. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12040644

AMA Style

Sirisubtawee S, Koonprasert S, Sungnul S. New Exact Solutions of the Conformable Space-Time Sharma–Tasso–Olver Equation Using Two Reliable Methods. Symmetry. 2020; 12(4):644. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12040644

Chicago/Turabian Style

Sirisubtawee, Sekson, Sanoe Koonprasert, and Surattana Sungnul. 2020. "New Exact Solutions of the Conformable Space-Time Sharma–Tasso–Olver Equation Using Two Reliable Methods" Symmetry 12, no. 4: 644. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12040644

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop