Skip to main content

ORIGINAL RESEARCH article

Front. Phys., 06 December 2019
Sec. Statistical and Computational Physics

New Solutions of Gardner's Equation Using Two Analytical Methods

  • 1Department of Engineering Science, Kermanshah University of Technology, Kermanshah, Iran
  • 2Department of Mathematics, Faculty of Engineering and Natural Sciences, Bahçeşehir University, Istanbul, Turkey
  • 3Department of Mathematics, Faculty of Arts and Sciences, Cankaya University, Ankara, Turkey
  • 4Department of Medical Research, China Medical University Hospital, China Medical University, Taichung, Taiwan
  • 5Institute of Space Sciences, Bucharest, Romania

This article introduces and applies new methods to determine the exact solutions of partial differential equations that will increase our understanding of the capabilities of applied models in real-world problems. With these new solutions, we can achieve remarkable advances in science and technology. This is the basic idea in this article. To accurately describe this, some exact solutions to the Gardner's equation are obtained with the help of two new analytical methods including the generalized exponential rational function method and a Jacobi elliptical solution finder method. A set of new exact solutions containing four parameters is reported. The results obtained in this paper are new solutions to this equation that have not been introduced in previous literature. Another advantage of these methods is the determination of the varied solutions involving various classes of functions, such as exponential, trigonometric, and elliptic Jacobian. The three-dimensional diagrams of some of these solutions are plotted with specific values for their existing parameters. By examining these graphs, the behavior of the solution to this equation will be revealed. Mathematica software was used to perform the computations and simulations. The suggested techniques can be used in other real-world models in science and engineering.

1. Introduction

It is difficult or impossible to determine the exact solution for many partial differential equations. In spite of these problems, in recent years a variety of efficient and practical methods have been proposed by mathematicians and physicists. Some of these methods are the exp-function method [1], the Darboux transformation [2], the Lie group analysis [3], the modified simple equation method [4], the homogeneous balance scheme [5], the sine-cosine method, and the tanh-coth method. Some new and effective attempts at determining solutions of partial differential equations can be found in [618].

The Gardner equation belongs to the category of integrable non-linear partial differential equations. The introduction of this equation is attributed to the famous mathematician Clifford Gardner in 1968 [19]. This equation can actually be generalized to the KdV equation. It is therefore sometimes referred to as the modified KdV equation. This equation is used in many areas of applications, such as hydrodynamics, plasma physics, and quantum field theory. This paper aims to employ two analytical methods to solve the following version of the integrable equation given by [20]

ut+k2uux+k3u2ux+k4uxxxx=0.    (1)

In this model, the dependent variable is u(x, t), and The independent variables x and t are the spatial and temporal variables, respectively. Abdul-Majid Wazwaz in [21] has obtained some multiple-soliton solutions for a variant of the equation called the Gardner-KP (GKP) equation. His approach is based on the Hirota's bilinear method. In [22] the authors have applied the mapping method to study the dynamics of solitary waves governed by Gardner's equation. This equation arises while studying the shallow water waves. The perturbed Gardner equation is also discussed in this article through the aid of He's semi-inverse variational approach. Very recently, a classification of Lie symmetries for the Gardner equation has been reported in [23]. They have also used the similarity transformation method to introduce the invariant solutions. Their solutions are of multisoliton, compacton, negaton, positon, and kink wave soliton types. Considering some suitable auxiliary dependent variables, the authors of [24, 25] have obtained some exact invariant solutions for the equation with non-local symmetries. By using the method of planar dynamical systems approach, in different parameter regions, the authors in [26] have constructed the bifurcation of phase portraits of a traveling wave system. The work of [27] presents the ill-posedness results for the initial value problem for the Gardner equation. In [28], a certain classification of single traveling wave solutions of the time-fraction Gardner equation is investigated. These forms of the Gardner equation can be utilized to model various physical phenomena, such as the non-linear propagation of ion acoustic waves in an unmagnetized plasma.

As can be seen, numerous numerical and analytical methods have been used to study this equation. That proves the importance of this equation. This is our main motivation for writing this article - to determine new solutions to this equation This paper is organized as follows. The analysis of the GERFM is outlined in section 2. The application of the method of solving (1) is presented in section 3. Also, to have a better insight into the resulting solutions, many numerical simulations are carried out in this section. Finally, some remarks are discussed in the last section.

2. The Analysis of the GERFM

The GERFM has recently been applied to solve many non-linear PDEs in some literature [2931]. The successful use of this method in solving different sets of equations has made it an efficient method for solving partial equations. In order to gain insight into the method, let us have a quick review of the method. The steps to apply this method include the following.

1. Consider the following general non-linear PDE as

N(ψ,ψx,ψt,ψxx,)=0.    (2)

For two unknown constants of μ, ν, we define the new variables of ψ = ψ(ϰ) and ϰ = μx − νt. then, Equation (28) can be reformulated as a non-linear ODE as

N(ψ,μψ,-νψ,μ2ψ,)=0.    (3)

2. Now, we take the solution Equation (29) into account for the following structure:

ψ(ϰ)=A0+k=1MAkΘ(ϰ)k+k=1MBkΘ(ϰ)-k.    (4)

where

Θ(ϰ)=r1es1ϰ+r2es2ϰr3es3ϰ+r4es4ϰ.    (5)

and ri, si(1 ≤ i ≤ 4), A0, Ak and Bk(1 ≤ kM) are unknown constants. Then, equating the two values of the amplitude, from (12) and (13), leads to the value of M.

3. Putting Equation (30) into Equation (29) and collecting all terms, the left-hand side of Equation (29) give us an algebraic equation P(Z1, Z2, Z3, Z4) = 0 in terms of Zi=esiϰ for i = 1, …, 4. Zeroing each coefficient of P, we get a system of non-linear equations in terms of ri, si(1 ≤ i ≤ 4), and μ, l, A0, Ak and Bk(1 ≤ kM).

4. Any symbolic computation software can be utilized to solve this system to determine the values of ri, si(1 ≤ i ≤ 4), A0, Ak, and Bk(1 ≤ kM). Using these results will direct us to soliton solutions of the main non-linear PDE.

3. Application of the Method

Below, we present a detailed presentation of the solution of Equation (1). To this end, let us consider the following new definitions

u(x,t)=U(ϰ),ϰ=μx-νt,    (6)

where μ and ν are arbitrary unknown parameters. Utilizing the wave transformation (36) converts Equation (1) into the following single NODE:

(k1ν-μ)U+k2νUU+k3νU2U+k4ν3U=0,    (7)

Performing the integral with respect to ϰ and with c = 0, the last equation becomes

(k1ν-μ)U+13k2νU2+13k3νU3+k4ν3U=0.    (8)

Then, equating the two values of 3M and M + 2, corresponding to U3 and U in Equation (8), leads to the value of M = 1. Using Equation (5) together with M = 1, we have

U(ϰ)=A0+A1Θ(ϰ)+B1Θ(ϰ).    (9)

Proceeding as outlined in the second section and depending on the values of the parameters we obtain in the solitary wave solutions.

Set 1:

One obtains r = [−3, −1, 1, 1] along with s = [2, 0, 2, 0], so (5) turns to

Θ(ϰ)=-3e2ϰ-1e2ϰ+1.    (10)

In this case we obtain two exact solutions, as:

I.

   μ=-k22k2-6k3k472k32k4,ν=k2-6k3k412k3k4,A0=k22k3,A1=0,B1=3k22k3.

Putting these values in Equations (10) and (37) yields a solitary wave solution for Equation (1) as:

u1(x,t)=-k2k3(3e2ϰ+1),    (11)

where

ϰ=-k32-6k3k472k32k4x-k2-6k3k46k3k4t.

II.

   μ=-k22k2-6k3k472k32k4,ν=k2-6k3k412k3k4,A0=-3k22k3,A1=0,B1=-3k22k3.

Putting these values in Equations (10) and (37) yields a solitary wave solution for Equation (1) as:

u2(x,t)=-3k2e2ϰk3(3e2ϰ+1),    (12)

where

ϰ=-k32-6k3k436k32k4x-k2-6k3k46k3k4t.

Set 2:

One obtains r = [−1, 3, 1, −1] along with s = [1, −1, 1, −1], so (5) turns to

Θ(ϰ)=cosh(ϰ)-2sinh(ϰ)sinh(ϰ).    (13)

In this case we obtain two exact solutions, as:

I.

μ=-k22k2-6k3k4576k32k4,ν=k2-6k3k424k3k4,A0=-k2k3,                                                  A1=-k24k3,B1=-3k24k3.

Now, from Equations (13) and (37) we will reach to a solitary wave solution for Equation (1) as:

u3(x,t)=-k22k3(sinh(2ϰ)-4sinh2(ϰ)),    (14)

where

ϰ=-k32-6k3k4576k32k4x-k2-6k3k424k3k4t.

II.

    μ=k2-3k3k4(7+4)(406k22-56k227)1944k32k4(57-7)(7+1)2,    ν=k2-3k3k4(7+4)36k3k4,A0=-k2(57+11)6k3(7+1),A1=-k2(7+1)12k3,B1=-59k27-119k2k3(7+1)3(57-7).

Equations (13) and (37) for these values will introduce a solitary wave solution for Equation (1) as:

u4(x,t)=2k23k3(1167112)cosh2(ϰ)(29728)sinh(2ϰ)577+231(577)(7+1)3(sinh(2ϰ)4sinh2(ϰ)),    (15)

where

ϰ=   k23k3k4(7+4)(406k22+14047k1k356k227+756k1k3)1944k32k4(577)(7+1)2x               k23k3k4(7+4)36k3k4t.

Set 3:

One obtains r = [3, 2, 1, 1] along with s = [1, 0, 1, 0], so (5) turns to

Θ(ϰ)=3eϰ+2eϰ+1.    (16)

In this case we obtain an exact solution, as:

I.

   μ=-k23-6k3k44500k32k4,ν=k2-6k3k430k3k4,A0=-k2k3,A1=k25k3,B1=6k25k3.

Putting these values in Equations (16) and (37) yields a solitary wave solution for Equation (1) as:

u5(x,t)=-k2eϰ5k3(1+eϰ)(3eϰ+2),    (17)

where

ϰ=(150k1k3-k22)k2-6k3k44500k32k4x-k2-6k3k430k3k4t.

Set 5:

One obtains r = [1, 1, 1, −1] along with s = [2, 0, 2, 0], so (5) turns to

Θ(ϰ)=e2ϰ+1e2ϰ-1.    (18)

In this case we obtain an exact solution, as:

I.

   μ=-k32-6k3k4144k32k4,ν=k2-6k3k424k3k4,A0=-k22k3,A1=-k24k3,B1=-k24k3.

For these solutions in Equations (18) and (37) yields a solitary wave solution for Equation (1) as:

u6(x,t)=-k2e4ϰk3(e4ϰ-1),    (19)

where

ϰ=-k32-6k3k436k32k4x-k2-6k3k46k3k4t.

Set 6:

One obtains r = [−2−i, 2−i, −1, 1] along with s = [i, −i, i, −i], so (5) turns to

Θ(ϰ)=cos(ϰ)+2sin(ϰ)sin(ϰ).    (20)

In this case we obtain an exact solution, as:

I.

   μ=k32-6k3k4576k32k4,ν=k2-6k3k424k3k4,A0=-k2k3,A1=k24k3,B1=5k24k3.

Inserting these values in Equations (20) and (37) yields a solitary wave solution for Equation (1) as:

u7(x,t)=k22k3(sin(2ϰ)+4sin2(ϰ)),    (21)

where

ϰ=k23-6k3k4576k32k4x-k2-6k3k424k3k4t.

Set 7:

One obtains r = [−3, −1, 1, 1] along with s = [1, −1, 1, −1], so (5) turns to

Θ(ϰ)=-2cosh(ϰ)-sinh(ϰ)cosh(ϰ).    (22)

In this case we obtain an exact solution, as:

I.

   μ=-k32-6k3k472k32k4,ν=k2-6k3k412k3k4,A0=-3k23k3,A1=0,B1=-3k22k3.

Putting these values in Equations (22) and (37) yields a solitary wave solution for Equation (1) as:

u8(x,t)=-3k2(cosh(ϰ)+sinh(ϰ))2k3(2cosh(ϰ)+sinh(ϰ)),    (23)

where

ϰ=-k32-6k3k472k32k4x-k2-6k3k412k3k4t.

Set 8:

One obtains r = [1 + i, 1 − i, 1, 1] along with s = [i, −i, i, −i], so (5) turns to

Θ(ϰ)=-sin(ϰ)+cos(ϰ)cos(ϰ).    (24)

In this case we obtain an exact solution, as:

I.

    μ=k2(-374k22-6k2217)-3k3k4(17+9)9216k32(17+1)2k4,    ν=k2-3k3k4(17+9)48k3k4,A0=-k2(517+13)8k3(17+1),A1=k2(17+1)16k3,B1=k2(17+1)8k3.

Using the above solutions in Equations (24) and (37) yields a solitary wave solution for Equation (1) as:

u9(x,t)=-k28k3(617-10)cos3(ϰ)-4(17+1)cos(ϰ)-(17+9)sin(ϰ)(17+1)(2cos(ϰ)3-cos(ϰ)),    (25)

where

ϰ=k2(-374k22-6k2217)-3k3k4(17+9)9216k32(17+1)2k4x    -k2-3k3k4(17+9)48k3k4t.

Set 9:

One obtains r = [−1, −2, 1, 1] along with s = [1, 0, 1, 0], so (5) turns to

Θ(ϰ)=-eϰ-2eϰ+1.    (26)

In this case we obtain an exact solution, as:

I.

   μ=k2(771k2273-2263k22)-3k3k4(373+41)147456k32(73+3)2k4,   ν=k2-3k3k4(373+41)96k3k4,A0=-k2(2573+171)32k3(73+3),A1=-k2(73+3)32k3,B1=-k2(73+3)16k3.

Inserting these values in Equations (26) and (37) yields a solitary wave solution for Equation (1) as:

u10(x,t)=-k2((773-75)e2ϰ+(2773-143)eϰ+1473-150)32k3(73+3)(1+eϰ)(eϰ+2),    (27)

where

ϰ=k2((921673+(771732263)k22)3k3k4(373+41)147456k32(73+3)2k4x   k23k3k4(373+41)96k3k4t.

It is worth mentioning that the necessary condition to establish the existence of the acquired solutions u1(x, t) − u10(x, t) is k3k4 < 0.

4. A Jacobi Elliptical Solutions Finder Method

In this part, we are going to obtain new exact soliton solutions to the equation under investigation, using a newly proposed method [32]. To this end, we will briefly review the steps of using the method.

1. The main purpose of this method is to solve an equation as follows:

N(ϕ,ϕx,ϕt,ϕxx,)=0.    (28)

2. Defining ϕ = ϕ(ϰ) and ϰ = μxlt, Equation (28) is converted to

N(ϕ,ϕ,ϕ,)=0,    (29)

where μ and l are two constants.

3. At this point, the symbolic form of the Equation (29) can be formulated as follows:

ϕ(ϰ)=α0+k=12NαkΘ(ϰ)kβ0+k=12NβkΘ(ϰ)k,    (30)

where the values of constants A0, B0 and Ak, Bk(1 ≤ k ≤ 2N) are so that (30) is a solution to the Equation (29).

4. The value of N in Equation (30) is obtained using the balance principles and Θ(ϰ) satisfies the following non-linear ODE:

Θ(ϰ)2=h0+h2Θ(ϰ)2+h4Θ(ϰ)4+h6Θ(ϰ)6,Θ(ϰ) =h2Θ(ϰ)+2h4Θ(ϰ)3+3h6Θ(ϰ)5,    (31)

where hi(i = 0, 2, 4, 6) are real constants.

5. The solution of the Equation (31) should be as follows

Θ(ϰ)=Φ(ϰ)fΦ(ϰ)2+g,    (32)

where Φ(ϰ)2 + g > 0, and Φ(ϰ) is the solution of the Jacobian elliptic equation

Φ(ϰ)2=l0+l2Φ(ϰ)2+l4Φ(ϰ)4,    (33)

and lj(j = 0, 2, 4) are constants need to be calculated, The relationships for f and g will also be as follows:

f=h4(l2-h2)(l2-h2)2+3l0l4-2l2(l2-h2),g=3h4l0(l2-h2)2+3l0l4-2l2(l2-h2),    (34)

under the constraint condition

h24(l2-h2)[9l0l4-(l2-h2)(2l2+h2)]+3h6[3l0l4-(l22-h22)]2=0.    (35)

6. It is known that solutions of Equation (33) are in terms of Jacobi elliptic solutions. Inserting both (33) and (32) into Equation (30), one gets the optical solutions of Equation (28). It should be noted that by using the limits in Table 2, the Jacobian elliptic functions used in the solutions reduce to the known triangular functions.

5. The Application of the Method

In this section, to begin solving the equation, we first introduce the following new variables

ϕ=U(ϰ),ϰ=μx-νt.    (36)

Then we will consider the balancing principles in Equation (8). So, one gets N = 1. So, the Equation (30) can be rewritten as follows

U(ϰ)=α0+α1Θ(ϰ)+α2Θ2(ϰ)β0+β1Θ(ϰ)+β2Θ2(ϰ).    (37)

The following results will be obtained using the method presented in section 4 of this article.

Set 11: We attain

  μ=-2νk2227k3,ν=ν,α0=-β0k23k3,α1=α1,α2=0,    β0=β0,β1=0,β2=0,h0=h0,h2=k2227ν2k3k4,h4=-α12k36ν2β02k4,h6=0.    (38)

Using No. 1 in Table 1 we have

U(ϰ)=-k23k3+α16-(sn(ϰ,m))2ν2k4α12k3(m4-m2-k24729ν4k32k42+1)-3+(m2+k2227ν2k3k4+1)(sn(ϰ,m))2,

provided that

(27m2ν2k3k4+27ν2k3k4+k22)(27m2ν2k3k4-54ν2k3k4+k22)(54m2ν2k3k4-27ν2k3k4-k22)=0.

The exact soliton solution to the equation will thus be determined as follows

u11(x,t)=-k23k3+α16-(sn(ϰ,m))2ν2k4α12k3(m4-m2-k24729ν4k32k42+1)-3+(m2+k2227ν2k3k4+1)(sn(ϰ,m))2,    (39)

where

ϰ=-2νk2227k3x-νt.
TABLE 1
www.frontiersin.org

Table 1. Jacobi elliptic solutions of Equation (33).

TABLE 2
www.frontiersin.org

Table 2. Jacobi elliptic functions and their limits.

Using No. 2 in Table 1 we have

U(ϰ)=-k23k3+α16-k4ν2((sn(ϰ,m))2-1)(m4-m2-k24729ν4k32k42+1)α12k3(2(cn(ϰ,m))2m2-(cn(ϰ,m))2k2227ν2k3k4-(cn(ϰ,m))2-3m2+3),

provided that

(27m2ν2k3k4+27ν2k3k4+k22)(27m2ν2k3k4-54ν2k3k4+k22)(54m2ν2k3k4-27ν2k3k4-k22)=0.

The exact soliton solution to the equation will thus be determined as follows

u12(x,t)=-k23k3+α16-k4ν2((sn(ϰ,m))2-1)(m4-m2-k24729ν4k32k42+1)α12k3(2(cn(ϰ,m))2m2-(cn(ϰ,m))2k2227ν2k3k4-(cn(ϰ,m))2-3m2+3),    (40)

where

ϰ=-2νk2227k3x-νt.

Set 12: We attain

μ=-2νk2227k3,ν=ν,α0=α0,α1=α1,α2=0,β0=0,β1=-3k3α1k2,β2=0,                       h0=-α02k2254α12k3ν2k4,h2=k2227ν2k3k4,h4=h4,h6=0.    (41)

Using No. 1 in Table 1 we have

U(ϰ)=-k23k3α1(α1(sn(ϰ,m))2h4(m4-m2-k24729ν4k32k42+1)(-3+(m2+k2227ν2k3k4+1)(sn(ϰ,m))2)+α0)(sn(ϰ,m))2h4(m4-m2-k24729ν4k32k42+1)(-3+(m2+k2227ν2k3k4+1)(sn(ϰ,m))2),

provided that

(27m2ν2k3k4-54ν2k3k4+k22)(54m2ν2k3k4-27ν2k3k4-k22)(27m2ν2k3k4+27ν2k3k4+k22)=0.

The exact soliton solution to the equation will thus be determined as follows

u13(x,t)=-k23k3α1(α1(sn(ϰ,m))2h4(m4-m2-k24729ν4k32k42+1)(-3+(m2+k2227ν2k3k4+1)(sn(ϰ,m))2)+α0)(sn(ϰ,m))2h4(m4-m2-k24729ν4k32k42+1)(-3+(m2+k2227ν2k3k4+1)(sn(ϰ,m))2),    (42)

where

ϰ=-2νk2227k3x-νt.

Using No. 4 in Table 1 we have

U(ϰ)=k2α1((sn(ϰ,m))2m21)(m4m2k24729ν4k32k42+1)(h4(m2k2227ν2k3k41)(dn(ϰ,m))2+3m2h4cn2(ϰ,m))+α03k3α1((sn(ϰ,m))2m21)(m4m2k24729ν4k32k42+1)(h4(m2k2227ν2k3k41)(dn(ϰ,m))2+3m2h4(cn(ϰ,m))2),    (43)

provided that

(27m2ν2k3k4-54ν2k3k4+k22)(54m2ν2k3k4-27ν2k3k4-k22)(27m2ν2k3k4+27ν2k3k4+k22)=0.

The exact soliton solution to the equation will thus be determined as follows

u14(x,t)=-k2α1((sn(ϰ,m))2m2-1)(m4-m2-k24729ν4k32k42+1)(h4(-m2-k2227ν2k3k4-1)(dn(ϰ,m))2+3m2h4cn2(ϰ,m))+α03k3α1((sn(ϰ,m))2m2-1)(m4-m2-k24729ν4k32k42+1)(h4(-m2-k2227ν2k3k4-1)(dn(ϰ,m))2+3m2h4(cn(ϰ,m))2),    (44)

where

ϰ=-2νk2227k3x-νt.

Using No. 7 in Table 1 we have

U(ϰ)=-k23k3α1(α1(sn(ϰ,m))2(m4-m2-k24729ν4k32k42+1)((m2+k2227ν2k3k4-2)(sn(ϰ,m))2-3(cn(ϰ,m))2)h4+α0)(sn(ϰ,m))2(m4-m2-k24729ν4k32k42+1)((m2+k2227ν2k3k4-2)(sn(ϰ,m))2-3(cn(ϰ,m))2)h4,    (45)

provided that

(27m2ν2k3k4-54ν2k3k4+k22)(27m2ν2k3k4+27ν2k3k4+k22)(54m2ν2k3k4-27ν2k3k4-k22)=0.

The exact soliton solution to the equation will thus be determined as follows

u15(x,t)=-k23k3α1(α1(sn(ϰ,m))2(m4-m2-k24729ν4k32k42+1)((m2+k2227ν2k3k4-2)(sn(ϰ,m))2-3(cn(ϰ,m))2)h4+α0)(sn(ϰ,m))2(m4-m2-k24729ν4k32k42+1)((m2+k2227ν2k3k4-2)(sn(ϰ,m))2-3(cn(ϰ,m))2)h4,    (46)

where

ϰ=-2νk2227k3x-νt.

Set 13: We attain

μ=-2νk2227k3,ν=ν,α0=α0,α1=-β1k23k3,α2=α2,β0=0,β1=β1,β2=0,h0=-α02k36k4β12ν2,                                          h2=-27α0α2k32-β12k2227β12ν2k3k4,h4=-α22k36k4β12ν2,h6=0.    (47)

Using No. 1 in Table 1 we have

U(ϰ)=(α0-18(sn(ξ,m))2ν2k4(m4-m2-α02k2281β02ν4k42+1)k2(-3+(m2-α0k29β0ν2k4+1)(sn(ξ,m))2))(1+3β16-(sn(ξ,m))2ν2k4k3(m4-m2-α02k2281β02ν4k42+1)β12k22(-3+(m2-α0k29β0ν2k4+1)(sn(ξ,m))2)),

provided that

(9m2β0ν2k4-18β0ν2k4-α0k2)(9m2β0ν2k4+9β0ν2k4-α0k2)(18m2β0ν2k4-9β0ν2k4+α0k2)=0.

The exact soliton solution to the equation will thus be determined as follows

u16(x,t)=(α0-18(sn(ξ,m))2ν2k4(m4-m2-α02k2281β02ν4k42+1)k2(-3+(m2-α0k29β0ν2k4+1)(sn(ξ,m))2))(1+3β16-(sn(ξ,m))2ν2k4k3(m4-m2-α02k2281β02ν4k42+1)β12k22(-3+(m2-α0k29β0ν2k4+1)(sn(ξ,m))2)),    (48)

where

ϰ=-2νk2227k3x-νt.

Using No. 8 in Table 1 we have

U(ϰ)=-k23k3α1(α1(sn(ϰ,m))2h4(m4-m2-k24729ν4k32k42+1)(-3+(m2+k2227ν2k3k4+1)(sn(ϰ,m))2)+α0)(sn(ϰ,m))2h4(m4-m2-k24729ν4k32k42+1)(-3+(m2+k2227ν2k3k4+1)(sn(ϰ,m))2)

provided that

(9m2β0ν2k4-18β0ν2k4-α0k2)(9m2β0ν2k4+9β0ν2k4-α0k2)(18m2β0ν2k4-9β0ν2k4+α0k2)=0.

The exact soliton solution to the equation will thus be determined as follows

u17(x,t)=-k23k3α1(α1(sn(ϰ,m))2h4(m4-m2-k24729ν4k32k42+1)(-3+(m2+k2227ν2k3k4+1)(sn(ϰ,m))2)+α0)(sn(ϰ,m))2h4(m4-m2-k24729ν4k32k42+1)(-3+(m2+k2227ν2k3k4+1)(sn(ϰ,m))2),    (49)

where

ϰ=ν(27k1k3-2k22)x27k3-νt.

Set 14: We attain

              μ=ν(432ν4β04h62k3k42+4ν2β02β22h6k22k4-4ν2β0β23h4k22k4+β24k1k22)k22β24,ν=ν,                          α0=36ν2β03h6k4β22k2,α1=0,α2=-36ν2β02h6k4β2k2,β0=β0,β1=0,β2=β2,h0=β02(2β0h6+β2h4)β23,h2=-β0(216ν2β03h62k3k4-β0β22h6k22-2β23h4k22)k22β24,h4=h4,h6=h6.    (50)

Using No. 1 in Table 1 we have

U(ϰ)=-36ν2h6k4(((m4-m2-Δ2+1)β2-h4(m2-Δ+1))(sn(ξ,m))2+3h4)β22k2(((m4-m2-Δ2+1)β2+h4(m2-Δ+1))(sn(ξ,m))2-3h4),

where Δ=β0(216ν2β03h62k3k4-β0β22h6k22-2β23h4k22)k22β24, provided that one of following conditions holds

(216ν2β04h62k3k4+m2k22β24+β02β22h6k22               +2β0β23h4k222k22β24)=0,(216ν2β04h62k3k4+m2k22β24            +β02β22h6k22+2β0β23h4k22+k22β24)=0,(216ν2β04h62k3k4+2m2k22β24β02β22h6k22               2β0β23h4k22k22β24)h42=0.

The exact soliton solution to the equation will thus be determined as follows

u18(x,t)=36ν2h6k4(((m4m2Δ2+1)β2    h4(m2Δ+1))(sn(ξ,m))2+3h4)β22k2(((m4m2Δ2+1)β2  +h4(m2Δ+1))(sn(ξ,m))23h4),    (51)

where

ϰ=ν(432ν4β04h62k3k42+4ν2β02β22h6k22k44ν2β0β23h4k22k4+β24k1k22)xk22β24ν t.

Using No. 5 in Table 1 we have

U(ϰ)=36ν2h6k4(3(cn(ξ,m))2m2h4+(m4m2Δ2+1)β2+h4(2m2+Δ1))β22k2(3(cn(ξ,m))2m2h4+(m4m2Δ2+1)β2+(2m2Δ+1)h4),

where Δ=β0(216ν2β03h62k3k4-β0β22h6k22-2β23h4k22)k22β24, provided that one of following conditions holds

(216 ν2β04h62k3k4+m2k22β24+β02β22h6k22              +2β0β23h4k222k22β24)=0,(216 ν2β04h62k3k4+m2k22β24+β02β22h6k22              +2β0β23h4k22+k22β24)=0,(216 ν2β04h62k3k4+2m2k22β24β02β22h6k22              2β0β23h4k22k22β24)h42=0.

The exact soliton solution to the equation will thus be determined as follows

u18(x,t)=         36ν2h6k4(3(cn(ξ,m))2m2h4+(m4m2Δ2+1)β2+h4(2m2+Δ1))β22k2(3(cn(ξ,m))2m2h4+(m4m2Δ2+1)β2  +(2m2Δ+1)h4)    (52)

where

ϰ=ν(432 ν4β04h62k3k42+4ν2β02β22h6k22k44ν2β0β23h4k22k4+β24k1k22)xk22β24ν t.

Likewise, other new families of solutions are obtained by following steps similar to the above using the following sets of parameters.

Set 15: We attain

   μ=-2α22k22ν162ν2β12h4k4+27α22k3,ν=ν,α0=0,α1=0,α2=α2,β0=0,β1=β1,β2=-18ν2β12h4k4+3α22k32α2k2,h0=0,h2=-2α22k2227k4(6ν2β12h4k4+α22k3)ν2,h4=h4,h6=0.    (53)

Set 16: We attain

 μ=ν(-3456ν4β04h42k22k3k42-48ν2β02β12h4k24k4)(216ν2β02h4k4k3+β12k22)2,ν=ν,α0=-576β12ν2β03h4k4k23(216ν2β02h4k4k3+β12k22)2,α1=-144ν2β02β1h4k4k2216ν2β02h4k4k3+β12k22,α2=0,β0=β0,β1=β1,β2=216ν2β02h4k4k3+β12k224k22β0,    (54)
h0=144β04h4(5184ν4β04h42k32k42144ν2β02β12h4k22k3k4+β14k24)k24(216 ν2β02h4k4k3+β12k22)4,                 h2=24h4(72ν2β02h4k4k3-β12k22)k22β02(216ν2β02h4k4k3+β12k22)2.

Set 17: We attain

   μ=-2νk2227k3,ν=ν,α0=-β0k23k3,α1=α1,α2=3(96ν2β02h4k4-α12k3)8β0k2,β0=β0,β1=0,β2=-9k3(96ν2β02h4k4-α12k3)8k22β0,h0=64k24h4β0481k32(96ν2β02h4k4-α12k3)2,h2=-k22(48ν2β02h4k4+α12k3)27k3k4(96ν2β02h4k4-α12k3)ν2,h4=h4,h6=0.    (55)

As can be seen, many varied sets of soliton solutions to Gardner's equation will be obtained by applying this method. In the structure of these solutions, rational, hyperbolic, trigonometric, exponential and Jacobi elliptical functions are used. The correctness of all the obtained answers has been carefully examined. All of these soliton solutions are new findings presented for the first time in this article.

6. Graphical Representation

We aimed to find new solutions for a given problem in Equation (1), and these new solutions should be described graphically. Thus, we present a graphical representation of some obtained solutions with the help of Mathematica in Figures 19. From these plots, some interesting and important physics phenomena can be observed.

FIGURE 1
www.frontiersin.org

Figure 1. Dynamic behavior of u1(x, t), u2(x, t) for k2 = 1.2, k3 = 2, k4 = −2. (A) u1(x, t). (B) u2(x, t).

FIGURE 2
www.frontiersin.org

Figure 2. Dynamic behavior of u5(x, t), u8(x, t) for k2 = 1.2, k3 = 2, k4 = −2. (A) u5(x, t). (B) u8(x, t).

FIGURE 3
www.frontiersin.org

Figure 3. Dynamic behavior of u9(x, t), and u10(x, t) for k2 = 1.1, k3 = −1.2, k4 = 0.5. (A) u9(x, t). (B) u10(x, t).

FIGURE 4
www.frontiersin.org

Figure 4. Dynamic behavior of |u11(x, t)| for k2 = 0.8, k3 = −0.5, k4 = 0.5., and α1 = 0.1. (A) m = 0.5 and ν=-2k23-6k3k4. (B) m = 0.8 and ν=5102k2306k3k4.

FIGURE 5
www.frontiersin.org

Figure 5. Dynamic behavior of |u12(x, t)| for k2 = 0.3, k3 = −0.1, k4 = 0.2. (A) m = 0.2 and ν=-53k263k3k4. (B) m = 0.5 and ν=-221k263k3k4.

FIGURE 6
www.frontiersin.org

Figure 6. Dynamic behavior of |u15(x, t)| for k2 = 0.2, k3 = −0.9, k4 = 1. (A) m = 0.3 and ν=-10k23-246k3k4. (B) m = 0.7 and ν=-10k23-6k3k4.

FIGURE 7
www.frontiersin.org

Figure 7. Dynamic behavior of |u16(x, t)| for k2 = 0.1, k3 = 0.7, k4 = −1.3, and α2 = 0.1, β1 = 1. (A) m = 0.1 and ν=56k3k4(27α0α2k32-β12k22)63k3k4β1. (B) m = 0.9 and ν=10-357k3k4(27α0α2k32-β12k22)1071k3k4β1.

FIGURE 8
www.frontiersin.org

Figure 8. Dynamic behavior of |u18(x, t)| for k2 = 0.1, k3 = 0.7, k4 = −1.3, and β2 = 0.1, h6 = 1. (A) m = 0.3 and ν=3k3k4(50β02h6+100β0β2h4+41β22)k2β2180β02h6k3k4. (B) m = 0.8 and ν=6k3k4(25β02h6+50β0β2h4-7β22)k2β2180β02h6k3k4.

FIGURE 9
www.frontiersin.org

Figure 9. Dynamic behavior of |u19(x, t)| for k2 = 0.7, k3 = 0.4, k4 = −0.5, and β2 = h6 = 1. (A) m = 0.1 and ν=3k3k4(50β02h6+100β0β2h4+49β22)k2β2180k3k4h6β02. (B) m = 0.7 and ν=3k3k4(50β02h6+100β0β2h4+β22)k2β2180k3k4h6β02.

7. Conclusion

In this manuscript, we have studied the Gardner equation with the help of two exact solution finder methods. A set of new exact solutions, including bright, kink, multi-soliton solutions, and singular solitons were found corresponding to four parameters, namely k1, k2, k3, and k4. The dynamic behavior of the acquired solutions was also demonstrated to deeply understand the features of the non-linear model. In order to better their properties, we have drawn some 3-D graphs. To the best of the authors knowledge, all the acquired results are novel findings, and cannot be found in the previous works. This result verifies the power of two suggested methods. The main advantages of the method are that they are very simple and quite efficient for the estimation of the optical solutions of PDES. Moreover, the proposed approaches represent efficient methodologies to investigate the exact solutions of the non-linear PDEs.

Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

Funding

This work was supported by JMDC (100%.2019.FPHY.827).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

1. Biazar J, Ayati Z. Extension of the Exp-function method for systems of two-dimensional Burgers equations. Comput Math Appl. (2009) 58:2103–6. doi: 10.1016/j.camwa.2009.03.003

CrossRef Full Text | Google Scholar

2. Zhaqilao, Qiao Z. Darboux transformation and explicit solutions for two integrable equations. J Math Anal Appl. (2011) 380:794–806. doi: 10.1016/j.jmaa.2011.01.078

CrossRef Full Text | Google Scholar

3. Yang S, Hua C. Lie symmetry reductions and exact solutions of a coupled KdV–Burgers equation. Appl Math Comput. (2014) 234:579–83. doi: 10.1016/j.amc.2014.01.044

CrossRef Full Text | Google Scholar

4. Younis M. A new approach for the exact solutions of nonlinear equations of fractional order via modified simple equation method. Appl Math. (2014) 5:1927–32. doi: 10.4236/am.2014.513186

CrossRef Full Text | Google Scholar

5. Wang M. Exact solutions for a compound KdV-Burgers equation. Phys Lett A. (1996) 213:279–87. doi: 10.1016/0375-9601(96)00103-x

CrossRef Full Text | Google Scholar

6. Djilali S. Impact of prey herd shape on the predator-prey interaction. Chaos Solit Fract. (2019) 120:139–48. doi: 10.1016/j.chaos.2019.01.022

CrossRef Full Text | Google Scholar

7. Djilali S, Bentout S. Spatiotemporal patterns in a diffusive predator-prey model with prey social behavior. Acta Appl Math. (2019). doi: 10.1007/s10440-019-00291-z. [Epub ahead of print].

CrossRef Full Text | Google Scholar

8. Djilali S. Herd behavior in a predator–prey model with spatial diffusion: bifurcation analysis and Turing instability. J Appl Math Comput. (2017) 58:125–49. doi: 10.1007/s12190-017-1137-9

CrossRef Full Text | Google Scholar

9. Djilali S, Touaoula TM, Miri SEH. A heroin epidemic model: very general non linear incidence, treat-age, and global stability. Acta Appl Math. (2017) 152:171–94. doi: 10.1007/s10440-017-0117-2

CrossRef Full Text | Google Scholar

10. Goufo EFD, Kumar S, Mugisha SB. Similarities in a fifth-order evolution equation with and with no singular kernel. Chaos Solit Fract. (2020) 130:109467. doi: 10.1016/j.chaos.2019.109467

CrossRef Full Text | Google Scholar

11. Odibat Z, Kumar S. A robust computational algorithm of homotopy asymptotic method for solving systems of fractional differential equations. J Comput Nonlin Dyn. (2019) 14:081004. doi: 10.1115/1.4043617

CrossRef Full Text | Google Scholar

12. El-Ajou A, Oqielat MN, Al-Zhour Z, Kumar S, Momani S. Solitary solutions for time-fractional nonlinear dispersive PDEs in the sense of conformable fractional derivative. Chaos. (2019) 29:093102. doi: 10.1063/1.5100234

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Cattani C, Rushchitskii YY. Cubically nonlinear elastic waves: wave equations and methods of analysis. Int Appl Mech. (2003) 39:1115–45. doi: 10.1023/b:inam.0000010366.48158.48

CrossRef Full Text | Google Scholar

14. Yang AM, Zhang YZ, Cattani C, Xie GN, Rashidi MM, Zhou YJ, et al. Application of local fractional series expansion method to solve Klein-Gordon equations on cantor sets. Abstr Appl Anal. (2014) 2014:1–6. doi: 10.1155/2014/372741

CrossRef Full Text | Google Scholar

15. Cattani C. Haar wavelet-based technique for sharp jumps classification. Math Comput Model. (2004) 39:255–78. doi: 10.1016/s0895-7177(04)90010-6

CrossRef Full Text | Google Scholar

16. Cattani C. Harmonic wavelet solutions of the Schrodinger equation. Int J Fluid Mech Res. (2003) 30:463–72. doi: 10.1615/interjfluidmechres.v30.i5.10

CrossRef Full Text | Google Scholar

17. Cattani C, Sulaiman TA, Baskonus HM, Bulut H. On the soliton solutions to the Nizhnik-Novikov-Veselov and the Drinfel'd-Sokolov systems. Opt Quant Electr. (2018) 50:138. doi: 10.1007/s11082-018-1406-3

CrossRef Full Text | Google Scholar

18. Avazzadeh Z, Heydari MH, Cattani C. Legendre wavelets for fractional partial integro-differential viscoelastic equations with weakly singular kernels. Eur Phys J Plus. (2019) 134:368. doi: 10.1140/epjp/i2019-12743-6

CrossRef Full Text | Google Scholar

19. Griffiths GW, Shiesser WE. Traveling Wave Analysis of Partial Differential Equations. Cambridge, MA: Elsevier (2012). doi: 10.1016/c2009-0-64536-0

CrossRef Full Text | Google Scholar

20. Wazwaz AM. Soliton solutions for two three plus one dimensional non-integrable KdV-type equations. Math Comput Model. (2012) 55:1845–8. doi: 10.1016/j.mcm.2011.11.082

CrossRef Full Text | Google Scholar

21. Wazwaz AM. Solitons and singular solitons for the Gardner–KP equation. Appl Math Comput. (2008) 204:162–9. doi: 10.1016/j.amc.2008.06.011

CrossRef Full Text | Google Scholar

22. Krishnan EV, Triki H, Labidi M, Biswas A. A study of shallow water waves with Gardner's equation. Nonlin Dyn. (2011) 66:497–507. doi: 10.1007/s11071-010-9928-7

CrossRef Full Text | Google Scholar

23. Kumar M, Tanwar DV. On Lie symmetries and invariant solutions of (2+1)–dimensional Gardner equation. Commun Nonlin Sci Numer Simul. (2019) 69:45–57. doi: 10.1016/j.cnsns.2018.09.009

CrossRef Full Text | Google Scholar

24. Fei J, Cao W, Ma Z. Nonlocal symmetries and explicit solutions for the Gardner equation. Appl Math Comput. (2017) 314:293–8. doi: 10.1016/j.amc.2017.07.002

CrossRef Full Text | Google Scholar

25. kai Liu Y, Li B. Nonlocal symmetry and exact solutions of the (2+1)-dimensional Gardner equation. Chin J Phys. (2016) 54:718–23. doi: 10.1016/j.cjph.2016.05.014

CrossRef Full Text | Google Scholar

26. Betchewe G, Victor KK, Thomas BB, Crepin KT. New solutions of the Gardner equation: analytical and numerical analysis of its dynamical understanding. Appl Math Comput. (2013) 223:377–88. doi: 10.1016/j.amc.2013.08.028

CrossRef Full Text | Google Scholar

27. Alejo MA. On the ill-posedness of the Gardner equation. J Math Anal Appl. (2012) 396:256–60. doi: 10.1016/j.jmaa.2012.06.018

CrossRef Full Text | Google Scholar

28. Cao D. The classification of the single traveling wave solutions to the time-fraction Gardner equation. Chin J Phys. (2019) 59:379–92. doi: 10.1016/j.cjph.2019.03.003

CrossRef Full Text | Google Scholar

29. Ghanbari B, Inc M. A new generalized exponential rational function method to find exact special solutions for the resonance nonlinear Schrödinger equation. Eur Phys J Plus. (2018) 133:142. doi: 10.1140/epjp/i2018-11984-1

CrossRef Full Text | Google Scholar

30. Osman MS, Ghanbari B, Machado JAT. New complex waves in nonlinear optics based on the complex Ginzburg-Landau equation with Kerr law nonlinearity. Eur Phys J Plus. (2019) 134:20. doi: 10.1140/epjp/i2019-12442-4

CrossRef Full Text | Google Scholar

31. Ghanbari B, Baleanu D, Qurashi MA. New exact solutions of the generalized Benjamin–Bona–Mahony equation. Symmetry. (2018) 11:20. doi: 10.3390/sym11010020

CrossRef Full Text | Google Scholar

32. Ghanbari B, Baleanu D. A novel technique to construct exact solutions for nonlinear partial differential equations. Eur Phys J Plus. (2019) 134:506. doi: 10.1140/epjp/i2019-13037-9

CrossRef Full Text | Google Scholar

Keywords: soliton solutions, generalized exponential rational function method, analytical solutions, PDE, computational, solitons, Gardner's equation

Citation: Ghanbari B and Baleanu D (2019) New Solutions of Gardner's Equation Using Two Analytical Methods. Front. Phys. 7:202. doi: 10.3389/fphy.2019.00202

Received: 10 October 2019; Accepted: 13 November 2019;
Published: 06 December 2019.

Edited by:

Carlo Cattani, Università degli Studi della Tuscia, Italy

Reviewed by:

Salih Djilali, University of Chlef, Algeria
Sunil Kumar, National Institute of Technology, Jamshedpur, India
Aliyu Isa Aliyu, Sun Yat-sen University, China

Copyright © 2019 Ghanbari and Baleanu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Behzad Ghanbari, b.ghanbary@yahoo.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.