Next Article in Journal
Asymmetrization of a Set of Degressively Proportional Allocations with Respect to Lexicographic Order. An Algorithmic Approach
Previous Article in Journal
Base Catalysis of Sodium Salts of [Ta6−xNbxO19]8− Mixed-Oxide Clusters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlocal Symmetry, Painlevé Integrable and Interaction Solutions for CKdV Equations

1
School of Computer Science, Shaanxi Normal University, Xi’an 710062, China
2
School of Information and Engineering, Xi’an University, Xi’an 710065, China
3
School of Mathematical Sciences, Liaocheng University, Liaocheng 252029, China
*
Authors to whom correspondence should be addressed.
Submission received: 12 June 2021 / Revised: 9 July 2021 / Accepted: 12 July 2021 / Published: 15 July 2021
(This article belongs to the Special Issue Differential and Difference Equations and Symmetries)

Abstract

:
In this paper, we provide a method to construct nonlocal symmetry of nonlinear partial differential equation (PDE), and apply it to the CKdV (CKdV) equations. In order to localize the nonlocal symmetry of the CKdV equations, we introduce two suitable auxiliary dependent variables. Then the nonlocal symmetries are localized to Lie point symmetries and the CKdV equations are extended to a closed enlarged system with auxiliary dependent variables. Via solving initial-value problems, a finite symmetry transformation for the closed system is derived. Furthermore, by applying similarity reduction method to the enlarged system, the Painlevé integral property of the CKdV equations are proved by the Painlevé analysis of the reduced ODE (Ordinary differential equation), and the new interaction solutions between kink, bright soliton and cnoidal waves are given. The corresponding dynamical evolution graphs are depicted to present the property of interaction solutions. Moreover, With the help of Maple, we obtain the numerical analysis of the CKdV equations. combining with the two and three-dimensional graphs, we further analyze the shapes and properties of solutions u and v.

1. Introduction and Motivation

Symmetry plays an important role in the construction of solutions for nonlinear (PDE), a continuous symmetry of a PDE system is a transformation that leaves invariant for the solution manifold of the system, i.e., it maps any solution of the system into a solution of the same system. As one of the main methods for studying differential equations, the Lie symmetry method has received extensive attention and has been developing vigorously since it was proposed by Sophus Lie [1]. The Lie symmetry method can not only obtain new solutions from old ones, but also reduce the dimensions of partial differential equations (PDEs) [2,3,4]. Once the reduction equations are derived, one can derive the Painlevé integral property through the reduced ODE and get several types of exact solutions. Nonlocal symmetry as a generalization of the symmetry was first studied by Vinogradov and Krasil’shchik in 1980 [5], which has close relationship with the integrable models. Compared with Lie point symmetry, the construction process of nonlocal symmetry is more complicated, but once the nonlocal symmetry of the equation is obtained, new types of symmetry reductions and exact solutions of some PDEs will also be obtained. Therefore, a variety of methods for constructing nonlocal symmetry have been proposed. For instance, Galas obtained the nonlocal symmetry from the pseudopotential [6]. Bluman, Euler et al. constructed nonlocal symmetry of PDEs by using potential system [7,8,9,10]. Lou and Hu started from the recursion operator and its inverse to construct the nonlocal symmetry for PDEs [11]. Lou, Hu and Chen derived the nonlocal symmetry from the Bäcklund transformation [12]. Lou and Reyes derived an infinite number of nonlocal symmetry of PDEs from a parameter dependent symmetry without using recursion operator [13,14]. Gao, Lou and Tang found that Painlevé analysis can be applied to obtain nonlocal symmetry corresponding to the residues with respect to the singular manifold of the truncated Painlevé expansion, which is also called residual symmetry [15]. Xin and Chen obtained the nonlocal symmetries from the auxiliary system (Lax pair) of PDEs [16,17].
This paper is motivated by Reference [16]; we start from the Lax pairs of the following CKdV equations
u t 1 2 u x x x + 3 2 v x x x 3 ( u v ) u x + 6 u v x = 0 , v t 1 2 v x x x + 3 2 u x x x 3 ( v u ) v x + 6 u x v = 0 ,
with u = u ( x , t ) , v = v ( x , t ) . By improving the classic Lie group method, along with altering the assumptions for symmetries, the local and nonlocal symmetries of Equation (1) are both derived firstly (Theorem 1). Since nonlocal symmetry cannot be applied to symmetry reduction directly, thus, a problem emerged of how to localize the nonlocal symmetry into local ones and how to use the obtained local symmetry to study the properties of the original PDEs.
Solving the above two problems is another innovative point of this article. Here, we give a method to localize the nonlocal symmetries by introducing the auxiliary potential variables, then the nonlocal symmetries are localized to Lie point symmetries of the extended system. When the prolonged system is regarded as a potential system, the nonlocal symmetries of the original equations can be obtained by calculating the Lie point symmetry of the entire system.
After localization the nonlocal symmetries, we employ the finite symmetry transformations theorem [6], and solving the initial value problems of the local symmetries, the corresponding finite symmetry transformations are derived (Theorem 2). That is to say, if we know the simple solution of Equation (1), then, by Theorem 2, other solutions for Equation (1) will be obtained. Moreover, we are interested in using standard Lie point symmetry approach to study the similarity reductions of the prolonged system. Then, one can not only study the Painlevé integral property of the CKdV Equation (1), but also construct some new interaction solutions for Equation (1) via localization procedures related with the nonlocal symmetry. This kind of solutions among different types of nonlinear excitation is hardly studied by other methods. This method not only can be applied in constant coefficient but also variable coefficient nonlinear systems [18,19,20,21,22,23,24,25,26].
The CKdV Equation (1) was first presented in Reference [27] to discuss the symmetry invariant and symmetry breaking soliton solutions of AB-KdV equations, which can be used to describe many physical phenomena, such as the internal-gravity-wave motion in a shallow stratified liquid [28], the atmospheric and oceanic blocking phenomena [29], etc. Equation (1) is integrable, which can be derived from the reduction of the Hirota-Satsuma systems [30,31,32,33], and the Lax pairs are as follows
ψ x x = P 1 ψ = u 0 0 v ψ ,
ψ t = P 2 ψ = 1 2 ( 3 v u ) x + ( u 3 v ) x 0 0 1 2 ( 3 u v ) x + ( v 3 u ) x ψ ,
where
ψ = ψ 1 ψ 2 ,
and the compatibility condition is
P 1 t P 2 x + [ P 1 , P 2 ] = 0 .
Enlightened by References [34,35], the rest contents are arranged as follows: In Section 2, we list some important definitions, theorems and provide a general method to construct nonlocal symmetry of nonlinear PDE. In Section 3, we apply the mehtod provided in Section 2 to study the nonlocal symmetries of the CKdV equations. In Section 4, by introducing potential variables to extend the original system, the nonlocal symmetries are transformed into the Lie point symmetries, and the corresponding finite transformation groups are given with Lie’s first theorem. In Section 5, two types of symmetry reductions of Equation (1) are discussed, the Painlevé property, interaction solutions and numerical analysis of Equation (1) are also derived. A short conclusion and further researches are included in the last Section.

2. Preliminaries

Firstly, we list some important definitions and theorems [7] which will be used later.
Definition 1.
A one-parameter Lie group of point transformations
( x * ) i = f i ( x , u ; ϵ ) , ( u * ) μ = g μ ( x , u ; ϵ ) ,
leaves the following PDE system
R σ ( x , u , u , , k u ) = 0 , ( x = ( x 1 , , x n ) , u ( x ) = ( u ( 1 ) ( x ) , , u ( m ) ( x ) ) ) ,
invariant if and only if its kth extension
X ( k ) = ξ i ( x , u ) x i + η μ ( x , u ) u μ + η i ( 1 ) μ ( x , u , u ) u i μ + + η i 1 i k ( k ) μ ( x , u , u , , k u ) u i 1 i k μ
leaves invariant the solution manifold of (5) in ( x , u , u , , k u ) , i.e., it maps any family of solution surfaces of PDE system (5). In this case, the one-parameter Lie group of point transformations (4) is called a point symmetry of the PDE system (5).
Definition 2.
A symmetry σ = σ ( u ) of a PDE system (5), is defined as a solution of its linearized equation,
R σ d d ϵ R ( u + ϵ σ ) | ϵ = 0 ,
which means (7) is invariant under the transformation
u u + ϵ σ ,
with the infinitesimal parameter ϵ.
Definition 3.
For PDE system (5), we define
V = ξ n ( x , u ) x n + η m ( x , u ) u m ,
be the infinitesimal generator of the Lie group of point transformations (4).
Definition 4.
For PDE system (5), if we derive the infinitesimal generators (9) of classical Lie point symmetries, then the corresponding finite symmetry transformations are given as the solution of the initial value problem
d x ^ d ϵ = ξ ( x ^ , t ^ , u ^ ) , x ^ ( ϵ = 0 ) = x , d t ^ d ϵ = τ ( x ^ , t ^ , u ^ ) , t ^ ( ϵ = 0 ) = t , d u i ^ d ϵ = η u i ( x ^ , t ^ , u ^ ) , u i ^ ( ϵ = 0 ) = u i .

Method for Seeking Nonlocal Symmetries

In this subsection, we give the concrete steps to construct nonlocal symmetry for system (5). For simplicity, we consider the case n = 2 , m = 1 , i.e., ( x 1 , x 2 ) = ( x , t ) .
Step 1. According to Definition 2, we write the symmetry equation of the system (5).
Step 2. Choosing the appropriate auxiliary systems. Usually, we can use the Lax pair, Bäcklund transformation, potential system, pseudo-potential, etc. and the general forms are as follows,
R α ( x , y , t , u , u x , u t , , ψ x , ψ t , ψ x x , ψ x t , ψ t t , , ψ λ x , ψ μ t ) = 0 , α Z +
where ψ = ( ψ 1 , ψ 2 , , ψ β ) denote β auxiliary variables and ψ λ x denote λ th-order partial derivatives of x, ψ μ t denote μ th-order partial derivatives of t.
Let U R be the space representing the single coordinate u, the space U 1 is isomorphic to R 2 with coordinates ( u x , u t ) . Similarly, U 2 R 3 has the coordinates representing the second order partial derivatives of u, and in general, U k R k + 1 , since there are k + 1 distinct k-th order partial derivatives of u. Finally, the space U ( k ) = U × U 1 × × U k with coordinates U ( k ) = ( u ; u x , u t ; u x x , u x t , u t t ; ) .
Step 3. Extending the original space X × U to the prolonged space X × U ( n ) with coordinates ( x , t , u , u , u x , u t , ) . V ˜ ( n ) denotes the n-th prolongation of V, which is a vector field on the n-jet space X × U ( n ) with the following form
V ˜ ( n ) = i = 1 2 ξ i x i + L η L u L .
then, we introduce the new variables into the coefficient function, and the transformed coefficient function ξ i , η L depend on the variables ( x , t , u , , ψ , ψ , ψ x , ψ t , , ψ λ x , ψ μ t ) , η 0 = η and η L have the form,
η L = D L u i = 1 2 u L D L ξ i .
Step 4. In order to construct the nonlocal symmetries, we need to solve the following equations
V ˜ ( n ) Δ v ( x , u ( n ) ) Δ v ( x , u ( n ) ) = 0 E q u a t i o n ( 3 ) = 0 .
Through the system (14), we obtain the coefficients of k u , ψ λ x , ψ μ t and let them equal to zero, which leads to a large number of determined equations for the coefficient functions ξ i and η L . Solving the determined equations with M a p l e , the local and nonlocal symmetry can be derived.
Next, take a CKdV equations as an example to illustrate the above processes. In this paper, we choose the Lax pair as the auxiliary system.

3. The Nonlocal Symmetry of the CKdV Equations

In order to construct nonlocal symmetry of the CKdV Equation (1), the first step is to write the symmetry equations by Definition 2.
The symmetry
σ σ 1 σ 2 ,
of CKdV (1) is defined as a solution of the following linearized system
σ 1 t 1 2 σ 1 x x x + 3 2 σ 2 x x x 3 σ 1 u x 3 u σ 1 x + 3 σ 2 u x + 3 v σ 1 x + 6 σ 1 v x + 6 u σ 2 x = 0 , σ 2 t 1 2 σ 2 x x x + 3 2 σ 1 x x x 3 σ 2 v x 3 v σ 2 x + 3 σ 1 v x + 3 u σ 2 x + 6 B σ 1 x + 6 σ 2 u x = 0 .
Namely, (15) is form invariant under the transformation.
u v u v   +   ϵ σ 1 σ 2
Similar to the standard Lie symmetry method, the symmetries σ 1 , σ 2 can be written as
σ 1 = X ¯ u x + T ¯ u t U ¯ , σ 2 = X ¯ v x + T ¯ v t V ¯ .
However, there exists the essential difference with Lie symmetry. According to the second step, we assume X ¯ , T ¯ , U ¯ , V ¯ are functions of x , t , u , v , ψ 1 , ψ 2 , ψ 1 x , ψ 2 x .
To proceed, using the third step, solving the determined equations, we derive the following nonlocal symmetry theorem of the CKdV Equation (1).
Theorem 1.
The CKdV Equation (1) possesses the following nonlocal symmetries
σ 1 = ( c 1 x + c 3 ) u x + ( 3 c 1 t + c 2 ) u t + 2 c 1 u 2 c 4 ψ 1 ψ 2 x 6 c 4 ψ 2 ψ 1 x , σ 2 = ( c 1 x + c 3 ) v x + ( 3 c 1 t + c 2 ) v t 6 c 4 ψ 1 ψ 2 x 2 c 4 ψ 2 ψ 1 x + 2 c 1 v ,
where c i ( i = 1 , , 4 ) is arbitrary constant.
Via the forth step, we give the specific proof of the Theorem 1 as follows.
Proof. 
Substituting Equation (18) into Equation (16), eliminating u t , v t , and ψ 1 x x , ψ 1 t , ψ 2 x x , ψ 2 t in terms of the Lax pairs (2) and (3), we derive a system of determining equations for the functions X ¯ , T ¯ , U ¯ , V ¯ , which can be solved by symbolic computing tool M a p l e to give
T ¯ = 3 c 1 t + c 2 , U ¯ = 2 c 4 ψ 1 ψ 2 x + 6 c 4 ψ 2 ψ 1 x 2 c 1 u , X ¯ = c 1 x + c 3 , V ¯ = 6 c 4 ψ 1 ψ 2 x + 2 c 4 ψ 2 ψ 1 x 2 c 1 v .
Substituting Equation (20) into Equation (18), we derive the exact form of σ 1 , σ 2 are expression (19). □
Remark 1.
The prolongation of vector fields shows that the symmetries σ 1 , σ 2 are neither classical Lie point symmetries nor Lie-Bäcklund symmetries because they depend not only on the variables x , t , u , v but also on the auxiliary variables and their high order partial derivatives. More results may be obtained if we assume the coefficients ξ i , η L have integral terms of the auxiliary variable i.e., they are the functions of ( x , t , u , , ψ d x , ) .
However, if we rewrite the symmetry Equation (19) as two parts
σ 11 = ( c 1 x + c 3 ) u x + ( 3 c 1 t + c 2 ) u t + 2 c 1 u , σ 21 = ( c 1 x + c 3 ) v x + ( 3 c 1 t + c 2 ) v t + 2 c 1 v ,
and
σ 12 = 2 c 4 ψ 1 ψ 2 x 6 c 4 ψ 2 ψ 1 x , σ 22 = 6 c 4 ψ 1 ψ 2 x 2 c 4 ψ 2 ψ 1 x .
we find that σ 11 , σ 21 in the expressions (22) denote the classical Lie symmetries and σ 12 , σ 22 denote the nonlocal symmetries of the CKdV equations respectivly. For simplicity, letting c 1 = 0 , c 2 = 0 , c 3 = 0 , c 4 = 1 in (19), we obtain the following nonlocal symmetries for Equation (1)
σ 1 = 2 ψ 1 ψ 2 x 6 ψ 2 ψ 1 x , σ 2 = 6 ψ 1 ψ 2 x 2 ψ 2 ψ 1 x .

4. Localization of the Nonlocal Symmetry

As we know, nonlocal symmetry cannot be applied to construct explicit solutions for PDEs directly, but once it is localized, especially as Lie point symmetry, we can further use the obtained results to explore the properties of the equations. According to this idea, we introduce two auxiliary variables ψ 3 = ψ 1 x , ψ 4 = ψ 2 x to extend the original system, the closed prolonged system of the CKdV equations become
u t 1 2 u x x x + 3 2 v x x x 3 ( u v ) u x + 6 u v x = 0 , v t 1 2 v x x x + 3 2 u x x x 3 ( v u ) v x + 6 u x v = 0 , ψ 1 x x = u ψ 1 , ψ 2 x x = v ψ 2 , ψ 1 t = 3 2 ψ 1 v x 1 2 ψ 1 u x + ( u 3 v ) ψ 1 x , ψ 2 t = 3 2 ψ 2 u x 1 2 ψ 2 u x + ( v 3 u ) ψ 2 x , ψ 3 = ψ 1 x , ψ 4 = ψ 2 x .
According to the Definition 2, the symmetries equations of system (24) are as follows
σ 3 x x + σ 1 ψ 1 + u σ 3 = 0 , σ 4 x x + σ 2 ψ 2 + v σ 4 = 0 , σ 3 t 3 2 σ 3 v x + 8 ψ 1 2 ψ 2 x x 8 ψ 1 ψ 1 x ψ 2 x + 1 2 σ 3 u x u σ 3 x + 3 v σ 3 x = 0 , σ 4 t 3 2 σ 4 u x + 8 ψ 2 2 ψ 1 x x 8 ψ 2 ψ 1 x ψ 2 x + 1 2 σ 4 v x B σ 4 x + 3 u σ 4 x = 0 , σ 5 σ 3 x = 0 , σ 6 σ 4 x = 0 ,
with σ 1 and σ 2 being given by (23), ψ 1 , ψ 2 , ψ 3 and ψ 4 satisfying the following transformations
ψ 1 ψ 1 + ϵ σ 3 , ψ 2 ψ 2 + ϵ σ 4 , ψ 3 ψ 3 + ϵ σ 5 , ψ 4 ψ 4 + ϵ σ 6 .
Solving Equations (16) and (25) leads to
σ 3 = ψ 1 f , σ 4 = ψ 2 f ,
with f satisfying
f x = 2 ψ 1 ψ 2 , f t = 2 v ψ 1 ψ 2 + 2 u ψ 1 ψ 2 + 8 ψ 1 x ψ 2 x ,
and the symmetry equations for f are
σ 7 x = 2 ( σ 3 ψ 2 + ψ 1 σ 4 ) , σ 7 t 2 ψ 1 ψ 2 ( σ 2 + σ 1 ) 2 v ( σ 3 ψ 2 + ψ 1 σ 4 ) 2 u ( σ 3 ψ 2 + ψ 1 σ 4 ) 8 ( σ 3 x ψ 2 x + ψ 1 x σ 4 x ) = 0 .
What is more interesting is that the auxiliary-dependent variable f satisfies the Schwartzian form of the CKdV equations
( C C x + C t ) ( S + C ) = 0 ,
with
C = f t f x , S = f x x x f x 3 2 ( f x x f x ) 2 .
The Schwarzian form (29) is invariant under the Möbious transformation
f a + b f c + d f ( a d b c ) .
That is to say, Equation (29) bears three symmetries σ f = d 1 , σ f = d 2 f and σ f = d 3 f 2 with arbitrary constants d 1 , d 2 and d 3 .
Here, if take
σ 7 = σ f = f 2 ,
one can see that the nonlocal symmetries (23) in the original spaces { x , t , A , B , ψ 1 , ψ 2 } can be localized to Lie point symmetries
σ 1 = 2 ψ 1 ψ 2 x 6 ψ 2 ψ 1 x , σ 2 = 6 ψ 1 ψ 2 x 2 ψ 2 ψ 1 x , σ 3 = ψ 1 f , σ 4 = ψ 2 f , σ 5 = ψ 3 f + 2 ψ 1 2 ψ 2 , σ 6 = ψ 4 f + 2 ψ 1 ψ 2 2 , σ 7 = f 2 .
In the extend spaces { x , t , u , v , ψ 1 , ψ 2 , ψ 3 , ψ 4 , f } , the corresponding vector is
V ^ = 2 ( ψ 1 ψ 2 x + 3 ψ 2 ψ 1 x ) u 2 ( 3 ψ 1 ψ 2 x + ψ 2 ψ 1 x ) v + ψ 1 f ψ 1 + ψ 2 f ψ 2 + ( ψ 3 f + 2 ψ 1 2 ψ 2 ) ψ 3 + ( ψ 4 f + 2 ψ 1 ψ 2 2 ) ψ 4 + f 2 f .
According to Definition 4, we obtain the following symmetry group transformation Theorem.
Theorem 2.
If u , v , ψ 1 , ψ 2 , ψ 3 , ψ 4 , f are solutions of the prolonged system, so are u ^ , v ^ , ψ 1 ^ , ψ 2 ^ , ψ 3 ^ , ψ 4 ^ , f ^
u ^ = u + 2 ϵ ψ 1 ψ 4 1 + ϵ f + 6 ϵ ψ 2 ψ 3 1 + ϵ f 8 ϵ 2 ψ 1 2 ψ 2 2 ( 1 + ϵ f ) 2 , v ^ = v + 6 ϵ ψ 1 ψ 4 1 + ϵ f + 2 ϵ ψ 2 ψ 3 1 + ϵ f 8 ϵ 2 ψ 1 2 ψ 2 2 ( 1 + ϵ f ) 2 , ψ 1 ^ = ψ 1 1 + ϵ f , ψ 2 ^ = ψ 2 1 + ϵ f , ψ 3 ^ = ψ 3 1 + ϵ f 2 ϵ ψ 1 2 ψ 2 ( 1 + ϵ f ) 2 , ψ 4 ^ = ψ 4 1 + ϵ f 2 ϵ ψ 1 ψ 2 2 ( 1 + ϵ f ) 2 , f ^ = f 1 + ϵ f ,
where ϵ is an arbitrary group parameter.
Proof. 
Letting ψ 3 = ψ 1 x , ψ 4 = ψ 2 x , via Definition 4, by solving the following initial value problems:
d u ^ ( ϵ ) d ϵ = 2 ( ψ ^ 1 ψ ^ 4 + 3 ψ ^ 2 ψ ^ 3 ) ,    u ^ ( 0 ) = u , d v ^ ( ϵ ) d ϵ = 2 ( 3 ψ ^ 1 ψ ^ 4 + ψ ^ 2 ψ ^ 3 ) ,    v ^ ( 0 ) = v , d ψ ^ 1 ( ϵ ) d ϵ = ψ ^ 1 f ^ ,    ψ 1 ( 0 ) = ψ 1 , d ψ ^ 2 ( ϵ ) d ϵ = ψ ^ 2 f ^ ,    ψ ^ 2 ( 0 ) = ψ 2 , d ψ ^ 3 ( ϵ ) d ϵ = ψ ^ 3 f ^ 2 ψ ^ 1 2 ψ 2 ,    ψ 3 ^ ( 0 ) = ψ 3 , d ψ ^ 4 ( ϵ ) d ϵ = ψ ^ 4 f ^ 2 ψ ^ 1 ψ ^ 2 2 ,    ψ 4 ^ ( 0 ) = ψ 4 , d f ^ ( ϵ ) d ϵ = f 2 ^ ( ϵ ) ,    f ^ ( 0 ) = f ,
one can easily derive the solutions of the system (33) given in Theorem 3. □
From Theorem 2, we can get a pair of new solutions u ^ , v ^ from given solutions u , v for Equation (1) through the finite symmetry transformations (32). It should be mentioned that the last expression in (32) is nothing but the Möbious transformation of the Schwartz form (29).
To proceed, we study the Lie point symmetries for the prolonged systems (24) and (27) of this section. By employing the classical Lie symmetry method, we first suppose the vector of the symmetries has the following form
V = X x + T t + U u + V v + Ψ 1 ψ 1 + Ψ 2 ψ 2 + Ψ 3 ψ 3 + Ψ 4 ψ 4 + F f ,
where X , T , U , V , Ψ 1 , Ψ 2 , Ψ 3 , Ψ 4 , F are functions of x , t , u , v , ψ 1 , ψ 2 , ψ 3 , ψ 4 and f. In other words, the closed system is invariant under the following infinitesimal transformations
( x , t , u , v , ψ 1 , ψ 2 , ψ 3 , ψ 4 , f ) ( x + ϵ X , t + ϵ T , u + ϵ U , v + ϵ V , ψ 1 + ϵ Ψ 1 , ψ 2 + ϵ Ψ 2 , ψ 3 + ϵ Ψ 3 , ψ 4 + ϵ Ψ 4 , f + ϵ f )
with
σ 1 = X u x + T u t U , σ 2 = X v x + T v t V , σ 3 = X ψ 1 x + T ψ 1 t Ψ 1 , σ 4 = X ψ 2 x + T ψ 2 t Ψ 2 , σ 5 = X ψ 3 x + T ψ 3 t Ψ 3 , σ 6 = X ψ 4 x + T ψ 4 t Ψ 4 , σ 7 = X f x + T f t F .
Inserting the above σ i (i = 1, …, 7) into the linearized equations of the prolonged systems, i.e., (16), (25) and (28), eliminating u t , v t , ψ 1 x x , ψ 1 t , ψ 2 x x , ψ 2 t , f x , f t , and collecting the coefficients of the independent partial derivatives of dependent variables u , v , we derive a system of overdetermined linear equations for the infinitesimals X , T , U , V , Ψ 1 , Ψ 2 , Ψ 3 , Ψ 4 , F . Solving the determining equations will lead to
X = c 1 x + c 3 , T = 3 c 1 t + c 2 , U = 2 ( c 4 ψ 1 ψ 4 + 3 c 4 ψ 2 ψ 3 + c 1 u ) , V = 2 ( 3 c 4 ψ 1 ψ 4 + c 4 ψ 2 ψ 3 + c 1 v ) , Ψ 1 = ( c 4 f + c 5 ) ψ 1 , Ψ 2 = ( c 4 f + c 6 ) ψ 2 , Ψ 3 = 2 c 4 ψ 1 2 ψ 2 + ( c 4 f c 1 + c 5 ) ψ 3 , Ψ 4 = 2 c 4 ψ 1 ψ 2 2 + ( c 4 f c 1 + c 6 ) ψ 4 , F = ( c 4 f + c 1 + c 5 + c 6 ) f + c 7 ,
where c i , ( i = 1 , , 7 ) is arbitrary constant.
For convenience, we restate and summarize the above content in theorem as follows.
Theorem 3.
The prolonged CKdV equations possess the following Lie point symmetries
σ 1 = ( c 1 x + c 3 ) u x + ( 3 c 1 t + c 2 ) u t + 2 ( c 4 ψ 1 ψ 4 + 3 c 4 ψ 2 ψ 3 + c 1 u ) , σ 2 = ( c 1 x + c 3 ) v x + ( 3 c 1 t + c 2 ) v t + 2 ( 3 c 4 ψ 1 ψ 4 + c 4 ψ 2 ψ 3 + c 1 v ) , σ 3 = ( c 1 x + c 3 ) ψ 1 x + ( 3 c 1 t + c 2 ) ψ 1 t ( c 4 f + c 5 ) ψ 1 , σ 4 = ( c 1 x + c 3 ) ψ 2 x + ( 3 c 1 t + c 2 ) ψ 2 t ( c 4 f + c 6 ) ψ 2 , σ 5 = ( c 1 x + c 3 ) ψ 3 x + ( 3 c 1 t + c 2 ) ψ 3 t 2 c 4 ψ 1 2 ψ 2 + ( c 4 f c 1 + c 5 ) ψ 3 , σ 6 = ( c 1 x + c 3 ) ψ 4 x + ( 3 c 1 t + c 2 ) ψ 4 t 2 c 4 ψ 1 ψ 2 2 + ( c 4 f c 1 + c 6 ) ψ 4 , σ 7 = ( c 1 x + c 3 ) f x + ( 3 c 1 t + c 2 ) f t ( c 4 f + c 1 + c 5 + c 6 ) f + c 7 .
When setting c i = 0 ( i = 1 , , 7 , i 4 ) , and c 4 = 1 , the above symmetries (36) degenerate into Equation (23).

5. Symmetry Reduction to the CKdV Equation

For nonlinear partial differential equation(NPDE), we are interested in studying its exact solutions by different methods, and then to analysis the physical properties which contains. Currently, there are lots of methods to construct exact solutions for NPDE [36,37,38,39], in this section, we mainly construct some kind of group invariant solutions for the CKdV equations by classical Lie symmetry method. Firstly, we need to solve the symmetry constraint conditions σ i = 0 ( i = 1 , , 7 ) defined by (35) with (36), which is equivalent to solving the following characteristic equations
d x X = d t T = d u U = d v V = d ψ 1 Ψ 1 = d ψ 2 Ψ 2 = d ψ 3 Ψ 3 = d ψ 4 Ψ 4 = d f F ,
we obtain two nontrivial similar reductions and several substantial invariant solutions listed below.
Next, we discuss the similarity reductions for (1) under the following two cases: c 1 0 and c 1 = 0 , then the Painlevé integrable property and two kinds of nontrivial group invariant solutions of Equation (1) listed in the follows.

5.1. Symmetry Reduction and Painlevé Integrable to the CKdV Equations

In this section, we will consider the symmetry reduction of Equation (1) in case of c 1 0 , and obtain the Painlevé integrability of Equation (1) through the Painlevé integrability of the reduced ordinary differential equation.
Reduction 1.
c 1 0 .
Without loss of generality, we assume c 2 = c 3 = c 5 = c 6 = 0 . By solving the characteristic Equation (38) and letting c 1 2 4 c 4 c 7 144 c 1 2 = k 2 , the following two sub-cases of k 0 and k = 0 need to be considered.
Case 1.
k 0 . The similarity solutions are as follows
u = U ( z ) t 2 3 + 2 c 4 9 c 1 2 k 2 t 2 3 exp ( 2 3 F ( z ) ) { c 4 Ψ 1 2 ( z ) Ψ 2 2 ( z ) 1 ( 2 cosh ( t 1 ) 2 1 ) 2 3 c 1 k tanh ( t 1 ) ( Ψ 1 ( z ) Ψ 4 ( z ) + 3 Ψ 2 ( z ) Ψ 3 ( z ) ) tanh ( t 1 ) 2 + 1 ) } , v = V ( z ) t 2 3 + 2 c 4 9 c 1 2 k 2 t 2 3 exp ( 2 3 F ( z ) ) { c 4 Ψ 1 2 ( z ) Ψ 2 2 ( z ) 1 ( 2 cosh ( t 1 ) 2 1 ) 2 3 c 1 k tanh ( t 1 ) ( 3 Ψ 1 ( z ) Ψ 4 ( z ) + Ψ 2 ( z ) Ψ 3 ( z ) ) tanh ( t 1 ) 2 + 1 ) } , ψ 1 = Ψ 1 ( z ) exp ( 1 6 F ( z ) ) t 1 6 ( 2 c o s h ( t 1 ) 2 1 ) , ψ 2 = Ψ 2 ( z ) exp ( 1 6 F ( z ) ) t 1 6 ( 2 c o s h ( t 1 ) 2 1 ) , f = c 1 ( 12 tanh ( 2 t 1 ) k + 1 ) 2 c 4 , ψ 3 = 1 3 exp ( 1 2 F ( z ) ) ( 6 Ψ 3 ( z ) t sech ( 2 t 1 ) + 2 c 4 Ψ 2 ( z ) Ψ 1 2 ( z ) sinh ( t 1 ) c o s h ( t 1 ) t k c 1 ( sinh ( t 1 ) 2 + c o s h ( t 1 ) 2 ) 2 ) , ψ 4 = 1 3 exp ( 1 2 F ( z ) ) ( 6 Ψ 4 ( z ) t sech ( 2 t 1 ) + 2 c 4 Ψ 1 ( z ) Ψ 2 2 ( z ) sinh ( t 1 ) c o s h ( t 1 ) t k c 1 ( sinh ( t 1 ) 2 + c o s h ( t 1 ) 2 ) 2 ) ,
with t 1 = k ( ln t + F ( z ) ) and z = x t 1 3 .
Substituting (39) into the extended equations, we obtain
Ψ 3 ( z ) = 1 6 exp ( 1 3 F ( z ) ) ( Ψ 1 ( z ) F z ( z ) 6 Ψ 1 z ( z ) , Ψ 4 ( z ) = 1 c 4 Ψ 1 2 ( z ) exp ( 2 3 F ( z ) ) c 1 k 2 ( Ψ 1 ( z ) F z 2 ( z ) 6 Ψ 1 z ( z ) F z ( z ) + 6 Ψ 1 ( z ) F z z ( z ) ) , U ( z ) = 1 36 ( F z 2 ( z ) 6 F z z ( z ) ) 1 Ψ 1 ( z ) ( 4 k 2 Ψ 1 ( z ) F z 2 ( z ) + Ψ 1 z z ( z ) ) + 1 3 Ψ 1 ( z ) Ψ 1 z ( z ) F z ( z ) ) , V ( z ) = 1 36 ( ( 144 k 2 + 1 ) F z 2 ( z ) + 18 F z z ( z ) ) + 1 3 Ψ 1 ( z ) ( Ψ 1 z ( z ) F z ( z ) + 3 Ψ 1 z z ) 2 Ψ 1 z F z ( z ) Ψ 1 2 ( Ψ 1 z ( z ) F z ( z ) Ψ 1 ( z ) F z z ( z ) ) , Ψ 1 ( z ) = 6 F z ( z ) c 1 c 4 k exp ( Q 1 ( z ) d z ) , Ψ 2 ( z ) = 6 F z ( z ) c 1 c 4 k exp ( Q ( z ) d z ) ,
where F 1 ( z ) = F z ( z ) , Q 1 ( z ) = Q z ( z ) satisfy the following equations
6 F 1 z z F 1 ( z ) 48 k 2 F 1 4 ( z ) + 36 Q 1 2 ( z ) F 1 2 ( z ) 2 z F 1 2 ( z ) 9 F 1 z 2 + 6 F 1 ( z ) = 0 , 6 Q 1 z z ( z ) F 1 ( z ) + Q 1 ( z ) ( 4 z F 1 ( z ) 48 Q 1 2 ( z ) F 1 ( z ) 9 ) = 0 .
Solving the second equation of system (41) leads to
F 1 ( z ) = 9 2 Q 1 ( z ) 24 Q 1 3 ( z ) 2 z Q 1 ( z ) 3 Q 1 z z ( z ) .
Substituting the expression (42) into the first equation of system (41), we derive a fourth-order ordinary differential equation of Q ( z ) . Then, we use the Ablowitz–Ramani–Segur (ARS) algorithm to test the Painlevé property of Equation (41). The dominant behavior for the first equation of (41) around ( z z 0 ) is Q 1 ( z ) Q 0 ( z ) ( z z 0 ) r , where Q 0 ( z ) and z 0 are arbitrary constants. Via the Leading terms analysis leads to the following two branches:
Q 1 ( z ) 1 2 ( z z 0 ) r , Q 1 ( z ) 1 2 ( z z 0 ) r ,
and the resonant points { 1 , 1 , 4 , 5 } . After boring calculation, it indicates the system (41) is Painlevé integrable. By the relationship of partial differential equation and the corresponding reduction equation [40], we know that Equation (1) is Painlevé integrable.
Case 2.
k = 0 . The similarity solutions are as follows
u = U ( z ) t 2 3 8 c 4 2 Ψ 1 2 ( z ) Ψ 2 2 ( z ) 9 t 2 3 c 1 2 t 1 2 + 2 c 4 ( 3 Ψ 3 ( z ) Ψ 2 ( z ) + Ψ 4 ( z ) Ψ 1 ( z ) ) 3 t 2 3 c 1 t 1 , v = V ( z ) t 2 3 8 c 4 2 Ψ 1 2 ( z ) Ψ 2 2 ( z ) 9 t 2 3 c 1 2 t 1 2 + 2 c 4 ( Ψ 3 ( z ) Ψ 2 ( z ) + 3 Ψ 4 ( z ) Ψ 1 ( z ) ) 3 t 2 3 c 1 t 1 , ψ 1 = Ψ 1 ( z ) t 1 6 t 1 , ψ 2 = Ψ 2 ( z ) t 1 6 t 1 , f = c 1 ( t 1 + 6 ) 2 c 4 t 1 , ψ 3 = 2 c 4 Ψ 1 2 ( z ) Ψ 2 ( z ) 3 c 1 t 1 2 t + Ψ 3 ( z ) t 1 t , ψ 4 = 2 c 4 Ψ 1 ( z ) Ψ 2 2 ( z ) 3 c 1 t 1 2 t + Ψ 4 ( z ) t 1 t .
Substituting (43) into the extended equations, we obtain
U ( z ) = Ψ 1 z z ( z ) ψ 1 ( z ) , V ( z ) = Ψ 2 z z ( z ) ψ 2 ( z ) , Ψ 1 ( z ) = 6 c 1 F 1 ( z ) e x p ( Q 1 ( z ) d z ) 2 c 4 , Ψ 2 ( z ) = 6 c 1 F 1 ( z ) e x p ( Q 1 ( z ) d z ) 2 c 4 , Ψ 3 ( z ) = Ψ 1 z ( z ) , Ψ 4 ( z ) = Ψ 2 z ( z ) ,
where F 1 ( z ) = F z ( z ) and Q 1 ( z ) = Q z ( z ) satisfy the following equations
6 F 1 ( z ) F 1 z z ( z ) + 36 F 1 2 ( z ) Q 1 2 ( z ) 2 z F 1 2 ( z ) 9 F 1 z 2 ( z ) + 6 F 1 ( z ) = 0 , 6 F 1 ( z ) Q 1 z z ( z ) 48 Q 1 ( z ) F 1 ( z ) Q 1 2 ( z ) + 4 z F 1 ( z ) 9 Q 1 ( z ) = 0 .
For the system (44), applying the same method as subcase 1, we also derive the following two branches
Q 1 ( z ) 1 2 ( z z 0 ) r , Q 1 ( z ) 1 2 ( z z 0 ) r ,
and the resonant points { 1 , 1 , 4 , 5 } . Furthermore, the Painlevé integrable of system (44) can be obtained and the same property of Equation (1) is also guaranteed.

5.2. Symmetry Reduction and Group Invariant Solutions to the CKdV Equation

In this section, we consider the symmetry reduction of Equation (1) in case of c 1 = 0 , and get two kinds of group invariant solutions.
Reduction 2.
c 1 = 0 .
For simplicity, we take c 2 = 1 and let ( c 5 + c 6 ) 2 4 c 4 c 7 = l 2 . Then, the following two subcases need to be discussed.
Case 3.
l 0 . It leads to the following similarity solutions
u = U ( z ) 2 c 4 l ( Ψ 1 ( z ) Ψ 4 ( z ) + 3 Ψ 2 ( z ) Ψ 3 ( z ) ) tanh t 1 8 c 4 2 l 2 Ψ 1 2 Ψ 2 2 sech 2 t 1 , v = V ( z ) 2 c 4 l ( 3 Ψ 1 ( z ) Ψ 4 ( z ) + Ψ 2 ( z ) Ψ 3 ( z ) ) tanh t 1 8 c 4 2 l 2 Ψ 1 2 Ψ 2 2 sech 2 t 1 , ψ 1 = Ψ 1 ( z ) exp ( t ( c 5 c 6 ) 2 ) sech t 1 , ψ 2 = Ψ 2 ( z ) exp ( ( c 5 c 6 ) t 2 ) sech t 1 , ψ 3 = 2 exp ( ( c 5 c 6 ) t 2 ) l 2 c 4 Ψ 1 2 ( z ) Ψ 2 ( z ) tanh ( t 1 ) + 2 exp ( t ( c 5 c 6 ) 2 ) Ψ 3 sech t 1 , f = 1 2 c 4 ( 2 l tanh t 1 c 5 c 6 ) ,
where t 1 = l ( F ( z ) + t ) , and z = x c 3 t are similarity variables.
Substituting (45) into the extended system, we derive
U ( z ) = 4 c 4 2 Ψ 1 2 ( z ) Ψ 2 2 ( z ) l 2 Ψ 1 z z ( z ) Ψ 1 ( z ) , V ( z ) = 4 c 4 2 Ψ 1 2 ( z ) Ψ 2 2 ( z ) l 2 Ψ 2 z z ( z ) Ψ 2 ( z ) , Ψ 1 ( z ) = l 2 2 F 1 ( z ) c 4 exp ( Q 1 ( z ) d z ) , Ψ 2 ( z ) = l 2 2 F 1 ( z ) c 4 exp ( Q 1 ( z ) d z ) , Ψ 3 ( z ) = Ψ 1 z ( z ) , Ψ 4 ( z ) = Ψ 2 z ( z ) , F 1 ( z ) = 6 Q 1 ( z ) 32 Q 1 3 ( z ) 8 c 3 Q 1 ( z ) 4 Q 1 z z ( z ) + c 5 c 6 ,
where Q 1 ( z ) = Q z ( z ) satisfies the following elliptic equation
Q 1 z 2 ( z ) = a 0 + a 1 Q 1 ( z ) + a 2 Q 1 2 ( z ) + a 3 Q 1 3 ( z ) + a 4 Q 1 4 ( z ) ,
with { a 0 = 4 l 2 a 3 2 , a 1 = 1 2 ( c 5 c 6 ) , a 2 = 2 c 3 , a 4 = 4 } or { a 0 = 0 , a 1 = c 6 c 5 4 , a 2 = 2 c 3 , a 3 = c 5 c 6 2 c 3 + 2 l 2 c 3 ( c 5 c 6 ) } .
Since Q 1 ( z ) in Equation (47) can be expressed as different types of Jacobian elliptic functions, we can derive the interactions solutions between soliton and many kinds of cnoidal periodic waves through the (45) and (46). In order to show this soliton and cnodial waves more clearly, we take a special solution of Equation (47) as
Q 1 ( z ) = μ 1 + μ 2 JacobiSN ( z , m ) .
Substituting Equation (48) into (47) leads to
a 3 = 16 μ 1 , c 3 = 2 μ 2 2 12 μ 1 2 + 1 2 , c 6 = 32 μ 1 3 16 μ 1 μ 2 2 + c 5 4 μ 1 , l = 8 μ 1 4 μ 1 4 + 4 μ 1 2 μ 2 2 + μ 1 2 μ 2 2 , m = 2 μ 2 .
Substituting expressions (49) into (47) and (46), we obtain the exact forms of F 1 ( z ) , Ψ 1 ( z ) , Ψ 2 ( z ) , Ψ 3 ( z ) , Ψ 4 ( z ) , U ( z ) , V ( z ) , and then we put them into the first two expressions of Equation (45), we will get the exact solutions of u and v. The result is omitted here because of its prolixity, but the corresponding images of u and v are as follows with the parameters { μ 1 = 2 , μ 2 = 1 , c 3 = 2 , c 5 = 1 , c 6 = 1 3 , l = 10 }.
At the same time, we give the numerical solution of u at x = 0 as follows.
With the help of M a p l e , we give the numerical analysis of u at x = 0 . By taking the value of time t at some key points, the result of numerical analysis is almost completely consistent with the two-dimensional graph Figure 1a of the exact solution. From Table 1, we can see that the value of u changes periodically from t = 7 to t = 0.2516 , the value at the highest point is approximately 8.34 and the value at the lowest point is approximately 2.21. When t = 0.5049 , the value of u suddenly becomes smaller, and then it varies periodically. Combining the two-dimensional and three-dimensional graph, we know that the soliton and the periodic wave interact at this moment to form a kink soliton solution, and the collision between soliton and periodic wave is an elastic collision. Namely, the shape of the wave does not change before and after collision.
At the same time, we give the numerical solution of v at x = 0 as follows.
Similar to Table 1, the numerical analysis of v at x = 0 also be derived. By taking the value of time t at some key points, the numerical analysis is almost completely consistent with the two-dimensional graph Figure 2a of the exact solution. From Table 2, we can see that the value of v reaches the maximum at t = 2 , which is close to 29.07. That is the result of the interaction between the soliton and the elliptical periodic wave. The collision is an elastic collision. Since, when t = 2 to 1.2234 and t = 2 to 5.7031, the values of v present a periodic characteristic, the maximum value is 6.70, and the minimum value is 1.36.
Similarly, we can get the change of u with time t at different positions of x, and the change rule is consistent with that at x = 0 .
Case 4.
l = 0 . It leads to the following similarity solutions
u = U ( z ) 8 c 4 2 Ψ 1 2 ( z ) Ψ 2 2 ( z ) t 1 2 + 2 c 4 t 1 ( Ψ 1 ( z ) Ψ 4 ( z ) + 3 Ψ 2 ( z ) Ψ 3 ( z ) ) , v = V ( z ) 8 c 4 2 Ψ 1 2 ( z ) Ψ 2 2 ( z ) t 1 2 + 2 c 4 t 1 ( 3 Ψ 1 ( z ) Ψ 4 ( z ) + Ψ 2 ( z ) Ψ 3 ( z ) ) , ψ 1 = Ψ 1 ( z ) t 1 exp ( t ( c 5 c 6 ) 2 ) , ψ 2 = Ψ 2 ( z ) t 1 exp ( t ( c 6 c 5 ) 2 ) , ψ 3 = 2 c 4 Ψ 1 2 ( z ) Ψ 2 ( z ) exp ( t ( c 5 c 6 ) 2 ) t 1 2 + exp ( t ( c 5 c 6 ) 2 ) Ψ 3 ( z ) t 1 , ψ 4 = 2 c 4 Ψ 1 ( z ) Ψ 2 2 ( z ) exp ( t ( c 6 c 5 ) 2 ) t 1 2 + exp ( t ( c 6 c 5 ) 2 ) Ψ 4 ( z ) t 1 , f = ( c 5 + c 6 ) ( t + F ( z ) ) + 2 2 c 4 t 1 ,
where t 1 = t + F ( z ) and z = x c 3 t .
Substituting the expressions (50) into the prolonged system, we derive
U ( z ) = ψ 1 z z Ψ 1 , V ( z ) = ψ 2 z z Ψ 2 , Ψ 1 ( z ) = 1 2 2 F 1 ( z ) c 4 exp ( Q 1 ( z ) d z ) , Ψ 2 ( z ) = 1 2 2 F 1 ( z ) c 4 exp ( Q 1 ( z ) d z ) , Ψ 3 ( z ) = Ψ 1 z ( z ) , Ψ 4 ( z ) = Ψ 2 z ( z ) , F 1 ( z ) = 6 Q 1 ( z ) 32 Q 1 3 ( z ) 8 c 3 Q 1 ( z ) 4 Q 1 z z ( z ) + c 5 c 6 ,
and Q 1 ( z ) satisfies elliptic equation
Q 1 z 2 ( z ) = b 0 + b 1 Q 1 ( z ) + b 2 Q 1 2 ( z ) + b 3 Q 1 3 ( z ) + b 4 Q 1 4 ( z ) ,
with { a 0 = 0 , a 1 = c 6 c 5 4 , a 2 = 2 c 3 , a 3 = c 5 c 6 2 c 3 , a 4 = 4 } , or { a 0 = 0 , a 1 = c 5 c 6 2 , a 2 = 2 c 3 , a 3 = a 3 , a 4 = 4 } .
For this subcase, according to the method of subcase 3, one can also get the interaction solutions of u and v, which reflect the interactions between Jacobian elliptic function and rational function.

6. Conclusions and Further Researches

In this paper, we provide the basic method to construct the nonlocal symmetry of nonlinear PDE, employ this method to derive the local and nonlocal symmetries of Equation (1), give the corresponding finite symmetry transformations theorem, and discuss the similarity reductions of the prolonged system. Furthermore, the Painlevé integral property and two types of new exact solutions are presented, which include the special interaction solutions between the soliton and the cnoidal periodic wave, and the soliton with rational function solution. This kind of solution can be easily applicable to the analysis of physically interesting processes. Besides, we also give the numerical analysis of the CKdV equations, which can verify the properties of exact solution.
For the provided method in this paper, we only consider the Lax pairs as the auxiliary system. However, if we choose the Bäcklund transformation, potential system, pseudo-potential as the auxiliary system, maybe, we can obtain different nonlocal symmetries.
Meanwhile, there is an other interesting thing that the CKdV Equation (1) have close relationship with the AB-KdV. If we let u = A , v = B = A P ^ T ^ = A ( x , t ) (where P ^ and T ^ are the usual parity and time-reversal operators respectively) and put the above assumption into Equation (1), then (1) will be reduced to the following integrable nonlocal AB-KdV equations
A t 1 2 A x x x + 3 2 B x x x 3 ( A B ) A x + 6 A B x = 0 , B = A P ^ T ^ = A ( x , t ) .
Therefore, for the symmetry reduction 1 and 2, Q ( z ) = 0 leads to Ψ 1 = Ψ 2 , then we can further deduce u = v and derive the PT (Parity and Time) symmetric exact solution for the AB-KdV equation by the transformations of x = x , t = t .
Based on the above facts, we consider whether all of the coupled equations can be reduced to a corresponding nonlocal AB equation. If so, whether we can get some new types of solutions of nonlocal equations by solving the corresponding local equations will be discussed in our future work.

Author Contributions

Investigation, Y.X., R.Y., Y.L., Software & Computing, Y.X., Methodology, X.X., Writing-original draft, Y.X., Analysis, Y.L., Validation, Y.X., Review & editing, R.Y. All authors have read and agreed to the published version of the manuscript.

Funding

The project is supported by the National Natural Science Foundation of China (Grant Nos. 12001424, 11471004), the Natural Science Basic research program of Shaanxi Province (No. 2021JZ-21), The Chinese Post doctoral Science Foundation (No. 2020M673332), Liaocheng University level science and technology research fund (No. 318012018), Research Award Foundation for Outstanding Young Scientists of Shandong Province (No. BS2015SF009), the Fundamental Research Funds for the Central Universities (No. 2020CBLY013).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lie, S. Sophus Lie’s 1880 Transformation Group Paper; Ackerman, M., Hermann, R., Eds.; Mathematical Sciene Press: Brookline, MA, USA, 1975. [Google Scholar]
  2. Freire, I.L.; Torrisi, M. Symmetry methods in mathematical modeling of Aedes aegypti dispersal dynamics. Nonlinear Anal-Real 2013, 14, 1300. [Google Scholar] [CrossRef]
  3. Hussain, A.; Bano, S.; Khan, I.; Baleanu, D.; Nisar, K.S. Lie Symmetry Analysis, Explicit Solutions and Conservation Laws of a Spatially Two-Dimensional Burgers-Huxley Equation. Symmetry 2020, 12, 170. [Google Scholar] [CrossRef] [Green Version]
  4. Aliyu, A.I.; Inc, M.; Yusuf, A.; Baleanu, D. Symmetry Analysis, Explicit Solutions, and Conservation Laws of a Sixth-Order Nonlinear Ramani Equation. Symmetry 2018, 10, 341. [Google Scholar] [CrossRef] [Green Version]
  5. Krasil’shchik, I.S.; Vinogradov, A.M. A method of calculating higher symmetries of nonlinear evolutionary equations, and nonlocal symmetries. Dokl. Akad. Nauk SSSR 1980, 253, 1289. [Google Scholar]
  6. Galas, F. New nonlocal symmetries with pseudopotentials. J. Phys. A Math. Gen. 1992, 25, L981. [Google Scholar] [CrossRef]
  7. Bluman, G.W.; Cheviakov, A.F.; Anco, S.C. Applications of Symmetry Methods to Partial Differential Equations; Springer: New York, NY, USA, 2010. [Google Scholar]
  8. Bluman, G.W.; Cheviakov, A.F. Framework for potential systems and nonlocal symmetries: Algorithmic approach. J. Math. Phys. 2005, 46, 123506. [Google Scholar] [CrossRef] [Green Version]
  9. Euler, N.; Euler, M. On Nonlocal Symmetries, Nonlocal Conseration Laws and Nonlocal Transformations of Evolution Equation: Two linearisable Hierarchies. J. Nonlinear Math. Phys. 2009, 16, 489. [Google Scholar] [CrossRef] [Green Version]
  10. Euler, M.; Euler, N.; Reyes, E.G. Multipotentializations and nonlocal symmetries: Kupershmidt, Kaup-Kupershmidt and Sawada-Kotera equations. J. Nonlinear Math. Phys. 2017, 24, 303. [Google Scholar] [CrossRef]
  11. Lou, S.Y.; Hu, X.B. Nonlocal Lie-Bäclund symmetries and Olver symmetries of the KdV equation. Chin. Phys. Lett. 1993, 10, 577. [Google Scholar] [CrossRef]
  12. Lou, S.Y.; Hu, X.R.; Chen, Y. Nonlocal symmetries related to Bäcklund transformation and their applications. J. Phys A Math. Gen. 2012, 45, 155205. [Google Scholar] [CrossRef] [Green Version]
  13. Lou, S.Y. Negative Kadomtsev-Petviashvili hierarchy. Phys. Scr. 1998, 57, 481. [Google Scholar] [CrossRef]
  14. Reyes, E.G. The modified Camassa-Holm equation. Int. Math. Res. Notices 2015, 12, 2617. [Google Scholar]
  15. Gao, X.N.; Lou, S.Y.; Tang, X.Y. Bosonization, singularity analysis, nonlocal symmetry reductions and exact solutions of supersymmetric KdV equation. J. High Energy Phys. 2013, 5, 029. [Google Scholar] [CrossRef] [Green Version]
  16. Xin, X.P.; Miao, Q.; Chen, Y. Nonlocal symmetry, optimal systems, and explicit solutions of the mKdV equation. Chin. Phys. B 2014, 23, 010203. [Google Scholar] [CrossRef]
  17. Miao, Q.; Xin, X.P.; Chen, Y. Nonlocal symmetries and explicit solutions of the AKNS system. Appl. Math. Lett. 2014, 28, 7. [Google Scholar] [CrossRef]
  18. Ren, B.; Cheng, X.P.; Lin, J. The (2+ 1)-dimensional Konopelchenko-Dubrovsky equation: Nonlocal symmetries and interaction solutions. Nonlinear Dyn. 2016, 86, 1855. [Google Scholar] [CrossRef]
  19. Huang, L.L.; Chen, Y. Nonlocal symmetry and similarity reductions for a (2 + 1)-dimensional Korteweg-de Vries equation. Nonlinear Dyn. 2018, 92, 221. [Google Scholar] [CrossRef]
  20. Xia, Y.R.; Xin, X.P.; Zhang, S.L. Residual symmetry, interaction solutions, and conservation laws of the (2+1)-dimensional dispersive long-wave system. Chin. Phys. B 2017, 26, 030202. [Google Scholar] [CrossRef]
  21. Xia, Y.R.; Xin, X.P.; Zhang, S.L. Nonlinear Self-Adjointness, Conservation Laws and Soliton-Cnoidal Wave Interaction Solutions of (2+1)-dimensional Modified Dispersive Water-Wave System. Commun. Theor. Phys. 2017, 67, 15. [Google Scholar] [CrossRef]
  22. Xin, X.P.; Liu, H.Z.; Zhang, L.L.; Wang, Z.G. High order nonlocal symmetries and exact interaction solutions of the variable coefficient KdV equation. Appl. Math. Lett. 2019, 88, 132. [Google Scholar] [CrossRef]
  23. Feng, Y.; Bilige, S.; Wang, X. Diverse exact analytical solutions and novel interaction solutions for the (2+1)-dimensional Ito equation. Phys. Scr. 2020, 95, 095201. [Google Scholar] [CrossRef]
  24. Kumar, S.; Kumar, A.; Kharbanda, H. Lie symmetry analysis and generalized invariant solutions of (2+1)-dimensional dispersive long wave (DLW) equations. Phys. Scr. 2020, 95, 065207. [Google Scholar] [CrossRef]
  25. Xin, X.P.; Zhang, L.L.; Xia, Y.R.; Liu, H.Z. Nonlocal symmetries and exact solutions of the (2+1)-dimensional generalized variable coefficient shallow water wave equation. Appl. Math. Lett. 2019, 94, 112. [Google Scholar] [CrossRef] [Green Version]
  26. Xia, Y.R.; Yao, R.X.; Xin, X.P. Nonlocal symmetries and group invariant solutions for the coupled variable-coefficient Newell-Whitehead system. J. Nonlinear Math. Phys. 2020, 27, 581. [Google Scholar] [CrossRef]
  27. Lou, S.Y. Alice-Bob systems, P ^ T ^ C ^ symmetry invariant and symmetry breaking soliton solutions. J. Math. Phys. 2018, 59, 083507. [Google Scholar] [CrossRef]
  28. Gear, J.A.; Grimshaw, R. Weak and srong interactions between interal solitary waves. Stud. Appl. Math. 1984, 70, 235. [Google Scholar] [CrossRef]
  29. Lou, S.Y.; Tong, B.; Hu, H.C.; Tang, X.Y. CKdV equations derived from two-layer fluids. J. Phys. A Math. Gen. 2006, 39, 513. [Google Scholar] [CrossRef] [Green Version]
  30. Hirota, R.; Satsuma, J. Soliton solutions of a couple Korteweg-de-Vries equation. Phys. Lett. A 1981, 85, 407. [Google Scholar] [CrossRef]
  31. Ramani, A.; Dorizzi, B.; Grammaticos, B. Integrability of Hirota-Satsuma equations: Two tests. Phys. Lett. A 1983, 99, 411. [Google Scholar] [CrossRef]
  32. Xu, M.H.; Jia, M. Exact Solutions, Symmetry Reductions, Painlevé Test and Bäcklund Transformations of A CKdV Equation. Commun. Theor. Phys. 2017, 68, 417. [Google Scholar]
  33. Parra, H.P.; Cisneros-Ake, L.A. The direct method for multisolitons and two-hump solitons in the Hirota-Satsuma system. Phys. Lett. A 2020, 384, 126471. [Google Scholar] [CrossRef]
  34. Milovanovic, G.V.; Parmar, R.K.; Rathie, A.K. Certain Laplace transforms of convolution type integrals involving product of two special pFp functions. Demonstr. Math. 2018, 51, 264. [Google Scholar] [CrossRef]
  35. Fitri, S.; Marjono, M.; Thomas, D.K.; Wibowo, R.B.E. Coefficient inequalities for a subclass of Bazilevic functions. Demonstr. Math. 2020, 53, 27. [Google Scholar] [CrossRef]
  36. Kumar, H.; Kumar, A.; Ch, F.; Singh, R.M.; Gautam, M.S. Construction of new traveling and solitary wave solutions of a nonlinear PDE characterizing the nonlinear low-pass electrical transmission lines. Phys. Scr. 2021, 96, 085215. [Google Scholar] [CrossRef]
  37. Riaz, M.B.; Atangana, A.; Jhangeer, A.; Rehman, M.J.U. Some exact explicit solutions and conservation laws of Chaffee-Infante equation by Lie symmetry analysis. Phys. Scr. 2021, 96, 084008. [Google Scholar] [CrossRef]
  38. Shen, Y.; Tian, B.; Liu, S.H.; Yang, D.Y. Bilinear Bäcklund transformation, soliton and breather solutions for a (3+1)-dimensional generalized Kadomtsev-Petviashvili equation in fluid dynamics and plasma physics. Phys. Scr. 2021, 96, 075212. [Google Scholar] [CrossRef]
  39. Fokou, M.; Kofane, T.C.; Mohamadou, A.; Yomba, E. Lump periodic wave, soliton periodic wave, and breather periodic wave solutions for third-order (2+1)-dimensional equation. Phys. Scr. 2021, 96, 055223. [Google Scholar] [CrossRef]
  40. Ablowitz, M.J.; Ramani, A.; Segur, H. A connection between nonlinear evolution equations and ordinary differential equations of P-type.I. J. Math. Phys. 1980, 21, 715. [Google Scholar] [CrossRef]
Figure 1. The kink soliton+cnoidal periodic wave solution for u. (a) The profile of the special structure with x = 0 . (b) The profile of the special structure with t = 0 . (c) The 3D soliton-cnodial periodic wave interaction solution to u.
Figure 1. The kink soliton+cnoidal periodic wave solution for u. (a) The profile of the special structure with x = 0 . (b) The profile of the special structure with t = 0 . (c) The 3D soliton-cnodial periodic wave interaction solution to u.
Symmetry 13 01268 g001
Figure 2. The bright soliton+cnoidal periodic wave solution for v. (a) The profile of the special structure with x = 0 , (b) The profile of the special structure with t = 0 , (c) The 3D soliton-cnodial periodic wave interaction solution to v.
Figure 2. The bright soliton+cnoidal periodic wave solution for v. (a) The profile of the special structure with x = 0 , (b) The profile of the special structure with t = 0 , (c) The 3D soliton-cnodial periodic wave interaction solution to v.
Symmetry 13 01268 g002
Table 1. Numerical analysis of u at x = 0 .
Table 1. Numerical analysis of u at x = 0 .
x000000000
t−7−6.2354−6−5.2777−4−3.5919−2.8640−2−1.9063
u2.213348.339768.317072.213445.062702.213528.339542.383912.21331
x000000000
t−1.1783−1−0.25160.504911.24011.474222.1255
u8.339438.452732.212355.359370.017791.3420 0.0171175.073705.35921
x000000000
t2.68812.925833.15993.811344.37154.52504.8492
u0.018161.342081.108750.017615.359245.026720.017801.341880.01737
x000000000
t55.497166.05736.30306.5315
u1.820855.359290.989800.017041.330290.01664
Table 2. Numerical analysis of v at x = 0 .
Table 2. Numerical analysis of v at x = 0 .
x000000000
t−2−2.0035−1.5320−1−0.9694−0.34120.30510.77051
u5.590655.549691.356796.555276.697611.356056.699811.463562.96166
x000000000
t1.22341.853422.08832.14182.49002.977033.6751
u6.696932.095754.6965229.071610.0393926.532381.356761.28513 6.64508
x000000000
t44.14794.615555.26195.7031
u2.456041.382546.696794.681811.509116.69698
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xia, Y.; Yao, R.; Xin, X.; Li, Y. Nonlocal Symmetry, Painlevé Integrable and Interaction Solutions for CKdV Equations. Symmetry 2021, 13, 1268. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13071268

AMA Style

Xia Y, Yao R, Xin X, Li Y. Nonlocal Symmetry, Painlevé Integrable and Interaction Solutions for CKdV Equations. Symmetry. 2021; 13(7):1268. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13071268

Chicago/Turabian Style

Xia, Yarong, Ruoxia Yao, Xiangpeng Xin, and Yan Li. 2021. "Nonlocal Symmetry, Painlevé Integrable and Interaction Solutions for CKdV Equations" Symmetry 13, no. 7: 1268. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13071268

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