This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.
Introduction
Many nonlinear phenomena are described by differential equations, especially by partial differential equations (PDEs). Symmetry group analysis of a differential equation appears as a powerful method to analyse PDEs [11] and fractional PDEs [13, 17]. Among its many applications, we highlight the fact that they allow us to obtain exact solutions of a PDE, directly or by using similarity solutions [1, 14]; classify invariant equations; reduce the number of independent variables and construct conservation laws [4, 7, 12, 21].
The symmetry group of a PDE is defined as the largest group of transformations acting on independent and dependent variables which transforms solutions of the equation into other solutions. Probably the most famous method used to obtain symmetries is the Lie classical method. The basic idea of Lie method is that, when a differential equation is invariant under a Lie group of transformations, there is a transformation which reduces the number of independent variables, in the case of a PDE; either the order of the equation, in the case of an ordinary differential equation (ODE).
Nevertheless, symmetry reductions for many PDEs are unobtainable by using classical symmetries. Accordingly, several generalizations of the Lie classical method have been established, for instance, the nonclassical method of Bluman and Cole [2]. This method generalizes the Lie classical method, however it is much more difficult to implement due to it leads to a system of nonlinear determining equations. This, couple with the fact that for some equations, such as the Korteweg-de Vries, the nonclassical method does not lead to new symmetries, it is appropriate to apply previously the Lie classical method. Many authors have used the nonclassical method to solve PDEs. In [10] Clarkson and Mansfield proposed an algorithm for calculating the determining equations associated to the nonclassical method. A different procedure for finding nonclassical symmetries was proposed by Bilă and Niesen in [3]. In [4] Bruzón and Gandarias extended the algorithm described for Bilă and Niesen to determine the nonclassical symmetries of a PDE for the case ξp = 0.
In the last decades travelling wave solutions of nonlinear PDEs have been studied [5]. A very successful method to obtain exact travelling wave solutions of numerous nonlinear PDEs is the method of simplest equation, especially its version called modified method of simplest equation. Simplest equation method is based on a procedure analogous to the first step of the test for the Painlevè property. In the modified simplest equation method, this procedure is substituted by the concept of balance equations [15, 16].
In this paper, we consider a generalized Gardner equation given by
where n is a positive constant, a and b are not simultaneously equal to zero, c≠ 0, d, e and f are arbitrary constants.
The Gardner equation, also known as combined KdV-mKdV equation, is widely used in various areas of physics, such that plasma physics, fluid dynamics, quantum field theory, and it is a useful model for the description of a great variety of wave phenomena in plasma and solid state. The Gardner equation has been recently considered by different authors [6, 12, 18, 20].
The aim of this paper is to study equation (1) from the point of view of symmetry reductions in PDEs. We obtain the classical and nonclassical symmetries of equation (1). Using the optimal system, we get the similarity variables and the similarity solutions which allow us to transform our equation into an ordinary differential equation. From these reductions we derive exact travelling waves solutions by using the simplest equation method given by Kudryashov. Some concluding remarks will end the paper.
Classical symmetries
In order to obtain the classical symmetries of equation (1) we apply the Lie classical method. This method is based on the determination of the symmetry group of a differential equation. For equation (1), a general element of the symmetry group is given by
The invariance of equation (1) under the infinitesimal generator (2) leads us to a set of determining equations for the unknown infinitesimals τ(x,t,u), ξ(x,t,u) and η(x,t,u) [19]. Simplifying this system we obtain that τ = τ(t), ξ = ξ(x,t) and η = η(x,t,u) must satisfy the following conditions:
$$\begin{array}{}
\displaystyle
{\bf v}_8= \displaystyle \left( 2 b x+ 4 b d t-a^2 t \right) \partial_x+ 6 b t \partial_t
- \left( 2 b u +a \right) \partial_u.
\end{array}$$
In the previous cases, k ≠ 0 represents an arbitrary constant.
Suppose that 𝒜 is an r-dimensional Lie algebra, and v1,v2,...vr, form a basis for 𝒜 . If we consider that two subalgebras are related by a transformation of the group of symmetries, invariant solutions calculated from them will be related by the same transformation. We construct the optimal system of subalgebras to obtain those invariant solutions that cannot be derived from others. By using the optimal system, we calculate the similarity variables and the similarity solutions. This allows us to transform equation (1) into an ODE, solving the characteristic system
In Table 1, we show the elements of the optimal system for each case in the symmetry classification along with their corresponding similarity variables and similarity solutions. Furthermore, in Table 2 we present the corresponding reduced ODEs.
Similarity solutions and similarity variables of equation (1)
λ e exp(ez) (h′ + f) - μ (e2x + kh′ - ehk + fk - de) = 0
2.1. a)
cμ3h‴ + kμhh′ + (dμ - λ)h′ + f = 0
2.1. b)
h′ + μd + (λ + μkz) f = 0
2.1. c)
3 c h‴ + 3k h h′-h′ z - 2 h = 0
2.2. a)
cμ3h‴ + bμh2nh′ + (dμ - λ)h′ = 0
2.2. b)
3 c n h‴ + 3 k n h2nh′- n h′ z - h = 0
2.3. a)
cμ3h‴ + μahnh′ + (dμ - λ)h′ = 0
2.3. b)
3 c n h‴ + 3 a n hn h′ - nh′ z - 2 h = 0
2.4. a)
cμ 3h‴ + bμh2h′ + aμhh′ + (dμ - λ )h′ = 0
2.4. b)
3 c h‴ + 3 b h2 h′ - h′ z - h = 0
Nonclassical symmetries
We apply the nonclassical method of Bluman and Cole [2] to get nonclassical symmetries of equation (1). Let us consider a s-th order PDE with p independent variables, x = (x1,...,xp), and one dependent variable, u = u(x),
where pr(s)v is the s-th prolongation of the vector field (14). This yields an overdetermined nonlinear system of equations for the infinitesimals.
To calculate the determining equations we implement the algorithm described in [9]. The application of this algorithm involves tedious, mechanical computations. Therefore, we make use of the Macsyma program symmgrp.max [8]. We can distinguish two different cases: τ ≠ 0 and τ = 0.
Case 1: τ ≠ 0
In the case that τ ≠ 0, without loss of generality, we may set τ(x,t,u) = 1, and we obtain a set of nine determining equations for the infinitesimals ξ (x,t,u) and η(x,t,u)
If $\begin{array}{}
\displaystyle
\displaystyle n\neq \frac{1}{2}
\end{array}$, 1, a = 0
$$\begin{array}{}
\displaystyle
\xi ={{c x+2 c d t+3 { k_1}+c}\over{3 c \left(t+{ k_1} \right)}}, \qquad \eta =-{{u}\over{3 n( t+ { k_1})}}.
\end{array}$$
If $\begin{array}{}
\displaystyle
\displaystyle n=\frac{1}{2}
\end{array}$, a = 0,
$$\begin{array}{}
\displaystyle
\xi ={{c x+2 c d t+3 { k_1}+c}\over{3 c \left(t+{ k_1}
\right)}}, \qquad \eta =-{{2 u}\over{3 \left(t+{ k_1}\right)}}.
\end{array}$$
If n = 1, b ≠ 0,
$$\begin{array}{}
\displaystyle
\xi ={{2 b c x+4 b c d t-a^2 c t+2 b c+6 k_1 b}\over{
6 b c \left(t+{k_1}\right)}}, \qquad
\eta ={{-2 b u-a}\over{6 b \left(t+{ k_1}\right)}}.
\end{array}$$
In the above cases, k1 represents an arbitrary constant.
Case 2: τ = 0
In the case τ = 0, without loss of generality, we may set ξ = 1 and the determining equation for the infinitesimal η are:
The complexity of this equation is the reason why we cannot solve (17) in general. Thus, we proceed by making ansatz on the form of η(x,t,u).
For $\begin{array}{}
\displaystyle
n=\displaystyle \frac{1}{2}
\end{array}$, a = 0, b = k or n = 1, a = k, b = 0, k ≠ 0 an arbitrary constant, choosing η = η(t), we find that the infinitesimal generator takes the form:
It is easy to check that generator (18) does not satisfy the Lie classical determining equations. From (18) the similarity variables for equation (1) have the form
with α and β unknown constants which must be determined. Taking the homogeneous balance between the highest order derivative h‴ and the nonlinear term of highest order h2h′ we obtain N = 2. Therefore, the solution of (25) takes the following form
Inserting (28) and the derivatives h′, h‴,..., into (25) we get a polynomial in Y(z) and its derivatives. Requiring the vanishing of the coefficients of the different powers of the function Y(z), we obtain an overdetermined system of equations
Equation (27) with β = 0 is the Bernoulli equation, so we can obtain the corresponding solution h of the ODE (25) in terms of this equation. As a result, the solution of the Bernoulli equation is
By using transformation (23) we can obtain a solution of equation (1) with n = 1.
Conclusions
In this paper, we have studied the classical and nonclassical symmetries admitted by the generalized Gardner equation (1). We have established a symmetry classification of equation (1) in terms of the arbitrary constants a, b, c, d, e, f and n. We have proved that the nonclassical method applied to this equation leads to new symmetries. Furthermore, similarity variables and similarity solutions of equation (1) have been obtained from the optimal system of subalgebras. Taking into account the elements of the optimal system of subalgebras we reduce equation (1) to an ordinary differential equation. Finally, we have constructed some travelling wave solutions by using the modified simplest equation method.