Next Article in Journal
Wavefront Recovery for Multiple Sun Regions in Solar SCAO Scenarios with Deep Learning Techniques
Next Article in Special Issue
Universal Bifurcation Chaos Theory and Its New Applications
Previous Article in Journal
A New Approach to Artificial Intelligent Based Three-Way Decision Making and Analyzing S-Box Image Encryption Using TOPSIS Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Construction of Infinite Series Exact Solitary Wave Solution of the KPI Equation via an Auxiliary Equation Method

1
School of Economics and Management, Huainan Normal University, Huainan 232038, China
2
Institute of Plasma Physics, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei 230031, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(6), 1560; https://doi.org/10.3390/math11061560
Submission received: 3 March 2023 / Revised: 16 March 2023 / Accepted: 20 March 2023 / Published: 22 March 2023
(This article belongs to the Special Issue Advances in Nonlinear Dynamics and Chaos: Theory and Application)

Abstract

:
The KPI equation is one of most well-known nonlinear evolution equations, which was first used to described two-dimensional shallow water wavs. Recently, it has found important applications in fluid mechanics, plasma ion acoustic waves, nonlinear optics, and other fields. In the process of studying these topics, it is very important to obtain the exact solutions of the KPI equation. In this paper, a general Riccati equation is treated as an auxiliary equation, which is solved to obtain many new types of solutions through several different function transformations. We solve the KPI equation using this general Riccati equation, and construct ten sets of the infinite series exact solitary wave solution of the KPI equation. The results show that this method is simple and effective for the construction of infinite series solutions of nonlinear evolution models.

1. Introduction

In modern physics and engineering, many physics phenomena, including hydrodynamics, optics, plasma, condensed matter physics, elementary particle physics, material physics, ocean engineering, astrophysics, and biology [1,2,3,4,5,6], can be expressed as nonlinear evolution equations (NLEEs) varying with time and space. A lot of methods have been proposed to obtain the exact traveling wave solutions of NLEEs, such as the F-expansion method [7,8,9,10,11], the tanh method and its extension [12,13,14], the Jacobi elliptic function method [15,16,17,18], the auxiliary equation method [19,20,21,22,23], the homogenous balance method [24,25], the trial function method [26], the (G’/G)-expansion method and its extension [27,28,29,30], and so on.
The Kadomtsev–Petviashvili (KP) equation [31], shown in Equation (1), was proposed by Kadomtsev and Petviashvili in 1970 to study the physical model of a two-dimensional shallow water wave, which describes how the water wave mainly propagates along the X-direction while the movement is very small in the Y-direction.
u t + 6 · u · u x + u x x x x + σ 2 · u y y = 0
This NLEE has important applications in fluid mechanics [32,33,34], plasma ion acoustic wave [35], nonlinear optics [36,37], and other fields. The KP equation is renamed the KPI equation when σ 2 = 1 , and the KPII equation when σ 2 = 1 . Compared with the KPII equation, the KPI equation has many more types of solutions [38,39,40,41,42,43,44]. In this work, we only discuss the solutions of KPI equation, shown as Equation (2):
u t + 6 · u · u x + u x x x x u y y = 0
To study various complex and diverse evolutionary processes of the KPI equation, we introduce a method to construct infinite series solutions of the KPI equation by using of Riccati equation as an auxiliary equation. Firstly, the general Riccati equation is solved to obtain many new types of solutions through several different function transformations. Then, we use this general Riccati equation as an auxiliary equation to solve the KPI equation. Finally, the infinite series exact solitary wave solutions are obtained. In this work, a simple Riccati equation is used to construct the form of the KPI equation’s solution as f ξ in Section 2. In Section 3, three different function transformations of f ξ are used to obtain ten sets of infinite-series hyperbolic function solutions, which are used to construct the infinite series exact solitary wave solution of the KPI equation. In Section 4, the infinite series exact solitary wave solutions are discussed in detail. Finally, a conclusion is given in Section 5.

2. The Form of KPI Equation’s Solutions

We assume that Equation (2) has the following traveling wave transformation:
u x , y , t = u ξ ξ = k x · x + k y · y + c · t
where k x , k y and c are the wave number in X and Y directions and wave speed, respectively. Equation (2) can be converted to an ODE, shown as Equation (4):
k x 4 · u + c k x k y 2 · u + 3 k x 2 · u 2 = 0
where u represents d u / d ξ . u ξ can be assumed as a finite series of defined f ξ , shown as Equation (5):
ξ = j = 0 n a j · f n ξ
where a j are constants, f ( ξ ) is the solution to be constructed later, and n can be obtained by the homogeneous balance between u and u 2 in Equation (4).
The highest degree of d p u / d ξ p is taken as:
O d p u / d ξ p = n + p , p = 1 , 2 , 3
We can obtain n = 2 , and u ξ can be expressed as:
u ξ = a 0 + a 1 · f ξ + a 2 · f 2 ξ
a 0 , a 1 and a 2 in Equation (7) are constants to be determined later, and f ξ is defined by simple Riccati equation, shown as:
f ξ = f 2 ξ + μ
Substituting Equation (8) into Equation (4), and setting the coefficients of f n ξ to zero, we get the following equation:
6 k x 4 · a 2 + 3 k x 2 · a 2 2 = 0 2 k x 4 · a 1 + 6 k x 2 · a 1 a 2 = 0 8 k x 4 · μ a 2 + c k x · a 2 k y 2 · a 2 + 3 k x 2 · a 1 2 + 6 k x 2 · a 0 a 2 = 0 2 k x 4 · μ a 1 + c k x · a 1 k y 2 · a 1 + 6 k x 2 · a 0 a 1 = 0 2 k x 4 · μ 2 a 2 + c k x · a 0 k y 2 · a 0 + 3 k x 2 · a 0 2 = 0
By solving these equations, we obtain the two solutions, shown as Equation (10):
a 0 = 2 k x 2 · μ / 3 a 1 = 0 a 2 = 2 k x 2 c = k y 2 4 k x 4 · μ / k x k x 0 k y 0   or   a 0 = 2 k x 2 · μ a 1 = 0 a 2 = 2 k x 2 c = k y 2 + 4 k x 4 · μ / k x k x 0 k y 0
u ξ could be expressed as v ξ in Equation (11) and w ξ in Equation (12):
v ξ = 2 k x 2 · μ / 3 2 k x 2 · f 2 ξ ξ = k x · x + k y · y + k y 2 4 k x 4 · μ / k x · t
w ξ = 2 k x 2 · μ 2 k x 2 · f 2 ξ ξ = k x · x + k y · y + k y 2 + 4 k x 4 · μ / k x · t
where k x 0 and k y 0 . Obviously, in order to get the solution of u ξ , we need to supplement the definition of f ξ to obtain the form of f ξ and the value of μ .

3. The Definition of f ξ

According to Equations (11) and (12), different sets of solutions u ξ could be determined through different definitions of f ξ . In the following, we adopt three definitions of f ξ to obtain ten sets of infinite-series hyperbolic function solutions for the KPI equation.

3.1. The First Definition of f ξ

f ξ is defined as Equation (13):
f ξ = k 0 + k 1 · g n ξ / g n ξ + r
where k 0 , k 1 , r are constants to be determined later.
While g n ξ is defined as Equation (14), and p n and q n are constants in these equations.
g n ξ 2 = p n · g n 2 ξ + q n
When n = 0 , g 0 ξ 2 = p 0 · g 0 2 ξ + q 0 . Obviously, this equation has three solutions:
g 01 ξ = s i n h ξ ( p 0 = 1 , q 0 = 1 ) g 02 ξ = c o s h ( ξ ) ( p 0 = 1 , q 0 = 1 ) g 03 ξ = s i n h ( ξ ) + ε · c o s h ( ξ ) ( p 0 = 1 , q 0 = 0 , ε 2 = 1 )
By substituting Equation (13) into Equation (8), using Equation (14), and then setting each coefficient of g n i ξ ( i = 0 , 1 , 2 , ) to zero, the following equations can be obtained.
k 0 = 0 , k 1 2 p n + μ = 0 , k 1 p n r = 2 μ r , k 1 q n = k 1 2 q n + μ r 2 ;   or   q n = 0 , k 0 2 + k 1 2 p n ± 2 k 0 k 1 p n + μ = 0 , k 1 p n r = 2 k 0 2 r ± 2 k 0 k 1 r p n + 2 μ r , k 0 2 r 2 + μ r 2 = 0 .
Solving them, we can obtain Equation (17).
k 0 = 0 , k 1 = 1 2 , μ = p n 4 , r = ± q n p n ;   or   q n = 0 , k 0 = ± 1 2 p n k 1 = 1 , μ = p n 4 , r 0 .
The constants p n and q n in Equation (14) are the key to obtain the values of constants k 0 , k 1 and r in Equation (13), so as to determine the f ξ . Therefore, two recursive relations of g n ξ are introduced below:
Case 1:
g n ξ = g n 1 2 ξ + a n 1
For n = 1 :
g 1 ξ 2 = p 1 · g 1 2 ξ + q 1 g 1 ξ = g 0 2 ξ + a 0
We obtain the following numerical relationship:
4 p 0 = p 1 , 4 q 0 = 2 p 1 a 0 , p 1 a 0 2 + q 1 = 0
Solving these equations, we get:
p 1 = 4 p 0 q 1 = q 0 2 / p 0 a 0 = q 0 / 2 · p 0
g 11 ξ = sinh 2 ( ξ ) + 1 2 ( p 0 = 4 , q 0 = 1 ) g 12 ξ = cosh 2 ξ 1 2 ( p 0 = 4 , q 0 = 1 ) g 13 ξ = [ sinh ( ξ ) + ε · c o s h ( ξ ) ] 2 ( p 0 = 4 , q 0 = 0 , ε 2 = 1 )
For n 2 :
Due to s i n h 2 ξ + 1 2 = c o s h 2 ξ 1 2 , the infinite series hyperbolic function solutions for g n 1 ξ are the same as that for g n 2 ξ . It is easy to obtain mathematical recurrence formula satisfying Equation (14), shown in Equation (23):
g n ξ = g n 1 2 ( ξ ) + q n 1 2 p n 1 ,   ( p n = 4 p n 1 , q n = q n 1 2 p n 1 ) .
We obtain two sets of infinite-series hyperbolic function. The first set is shown as Equation (24):
f 1 n ξ = g n 1 ξ 2 g n 1 ξ ± q n 1 / p n 1 μ = p n 1 / 4 , n = 1 , 2 , 3 , g n ξ = g n 1 2 ξ + q n 1 / 2 · p n 1 p n = 4 · p n 1 , q n = q n 1 2 / p n 1
The corresponding equations of n = 0 are g 0 ξ = s i n h ξ when p 0 = 1 and q 0 = 1 and g 0 ξ = c o s h ξ when p 0 = 1 and q 0 = 1 .
The second set is shown as Equation (25):
f 2 n ξ = ± 1 2 ν n ε ν n [ sinh ( ξ ) + ε · c o s h ξ ] ν n [ sinh ( ξ ) + ε · c o s h ξ ] ν n + r n = 1 , 2 , 3 ,
where μ = ν 2 · n 2 / 4 , r 0 , ν 0 , ε 2 = 1 .
Case 2:
g n ξ = g n 1 ξ + a n 1
For n = 1 :
g 1 ξ 2 = p 1 · g 1 2 ξ + q 1 g 1 ξ = g 0 ξ + a 0
We obtain the following numerical relationship:
p 0 / 4 = p 1 2 p 1 a 0 + q 1 = 0 q 0 / 4 = p 1 a 0 2 + q 1 a 0
Solving these equations, we get:
p 1 = p 0 / 4 q 1 = ± 2 p 0 q 0 a 0 = q 0 / p 0
g 14 ξ = s i n h ( ξ ) + ε ( p 1 = 1 / 4 , q 1 = ε / 2 , ε 2 = 1 ) g 15 ξ = c o s h ( ξ ) + ε p 1 = 1 / 4 , q 1 = ε / 2 , ε 2 = 1 g 16 ξ = s i n h ( ξ ) ± c o s h ( ξ ) ( p 1 = 1 4 , q 1 = 0 )
For n 2 :
It is easy to obtain mathematical recurrence formula satisfying Equation (14), shown as Equation (31):
g n ξ = g n 1 ξ ± q n 1 / p n 1 , p n = p n 1 / 4 , q n = p n 1 · q n 1 / 2
Besides one set of infinite series hyperbolic function solutions being the same to Equation (25), a new set of infinite series hyperbolic function solutions is provided in this case.
f 3 n ξ = g n 1 ξ 2 g n 1 ξ ± q n 1 / p n 1 μ = p n 1 / 4 , n = 1 , 2 , 3 , g n ξ = g n 1 ξ ± q n 1 / p n 1 p n = p n 1 / 4 , q n = 2 p n 1 · q n 1
The corresponding equations of n = 0 are g 0 ξ = s i n h ξ when p 0 = 1 and q 0 = 1 and g 0 ξ = c o s h ξ when p 0 = 1 and q 0 = 1 .

3.2. The Second Definition of f ξ

f ξ is defined as Equation (20):
f ξ = k 0 + k 1 · g n ξ · g n ξ / g n 2 ξ + r
where k 0 , k 1 , r are constants to be determined later. Similar to Section 3.1, substituting Equation (33) into Equation (8) and use Equation (14), we get Equation (34), which provides two sets solutions shown as Equation (35):
k 0 = 0 k 1 2 p n + μ = 0 k 1 q n + 2 k 1 p n r = k 1 2 q n + 2 μ r k 1 q n r = μ r 2   or   q n = 0 k 0 2 + k 1 2 p n ± 2 k 0 k 1 p n + μ = 0 2 k 1 p n r = 2 k 0 2 r ± 2 k 0 k 1 r p n + 2 μ r k 0 2 r 2 + μ r 2 = 0
k 0 = 0 k 1 = 1 μ = p n r = q n / p n   or   q n = 0 k 0 = ± p n k 1 = 2 μ = p n r 0
Corresponding to the first definition of g n ξ in Section 3.1, shown as Equation (18), two new sets of infinite-series hyperbolic function solutions are obtained, shown as Equations (36) and (37):
f 4 n ξ = g n 1 ξ g n 1 ξ g n 1 2 ξ + q n 1 / p n 1 μ = p n 1 , n = 1 , 2 , 3 , g n ξ = g n 1 2 ξ + q n 1 / 2 · p n 1 p n = 4 · p n 1 , q n = q n 1 2 / p n 1
The corresponding equations of n = 0 are g 0 ξ = s i n h ξ when p 0 = 1 and q 0 = 1 and g 0 ξ = c o s h ξ when p 0 = 1 and q 0 = 1 .
f 5 n ξ = ± ν n 2 ε ν n [ sinh ( ξ ) + ε · c o s h ξ ] 2 ν n [ sinh ( ξ ) + ε · c o s h ξ ] 2 ν n + r n = 1 , 2 , 3 ,
where μ = ν 2 n 2 , r 0 , ν 0 , ε 2 = 1 .
Corresponding to the second definition of g n ξ in Section 3.1, shown as Equation (26), the infinite series hyperbolic function solutions are same to the f 3 n ξ and f 5 n ξ , we do not get new form of solution.

3.3. The Third Definition of f ξ

Obviously, the function a 1 / f ξ also satisfies Equation (8) when a 1 = μ . We can easily obtain five new sets of infinite-series hyperbolic function solutions, which are shown below:
f 6 n ξ = p n 1 g n 1 ξ + r 2 g n 1 ξ μ = p n 1 / 4 , r = ± q n 1 / p n 1 g n ξ = g n 1 2 ξ + q n 1 / 2 · p n 1 p n = 4 · p n 1 , q n = q n 1 2 / p n 1
where n = 1 , 2 , 3 , . The corresponding equations of n = 0 are g 0 ξ = s i n h ξ when p 0 = 1 and q 0 = 1 and g 0 ξ = c o s h ξ when p 0 = 1 and q 0 = 1 .
f 7 n ξ = ν n [ sinh ( ξ ) + ε · c o s h ξ ] ν n + ν n r ( ± 2 4 ε ) [ sinh ( ξ ) + ε · c o s h ξ ] ν n ± 2 r n = 1 , 2 , 3 ,
where μ = ν 2 · n 2 / 4 , r 0 , ν 0 , ε 2 = 1 .
f 8 n ξ = p n 1 g n 1 ξ + r 2 g n 1 ξ μ = p n 1 / 4 , r = ± q n 1 / p n 1 g n ξ = g n 1 ξ ± q n 1 / p n 1 p n = p n 1 / 4 , q n = 2 p n 1 · q n 1
where n = 1 , 2 , 3 , . The corresponding equations of n = 0 are g 0 ξ = s i n h ξ when p 0 = 1 and q 0 = 1 and g 0 ξ = c o s h ξ when p 0 = 1 and q 0 = 1 .
f 9 n ξ = p n 1 [ g n 1 2 ξ + r ] g n 1 ξ g n 1 ξ μ = p n 1 , r = q n 1 / p n 1 g n ξ = g n 1 2 ξ + q n 1 / 2 · p n 1 p n = 4 · p n 1 , q n = q n 1 2 / p n 1
where n = 1 , 2 , 3 , . The corresponding equations of n = 0 are g 0 ξ = s i n h ξ when p 0 = 1 and q 0 = 1 and g 0 ξ = c o s h ξ when p 0 = 1 and q 0 = 1 .
f X n ξ = ν n [ sinh ( ξ ) + ε · c o s h ξ ] 2 ν n + ν n r ( ± 1 2 ε ) [ sinh ( ξ ) + ε · c o s h ξ ] 2 ν n ± r n = 1 , 2 , 3 ,
where μ = ν 2 n 2 , r 0 , ν 0 , ε 2 = 1 .
Through the definition of three different functional forms of f ξ , we obtain ten sets of infinite-series hyperbolic function solutions satisfying the Riccati equation (Equation (8)), which provide the base for the solution for the KPI equation shown as Equation (7).

4. The Discussion of u ξ

Due to the space limitation of the article, we only discuss n 4 of the ten sets of infinite-series hyperbolic function solutions v ξ and w ξ shown as Equations (11) and (12). Because the wave described by The KPI equation mainly propagates along the X-direction, we only consider v ξ and w ξ vs. x , t by assuming y = 0 .

4.1. The First Set of Solutions

According to the Equation (15), g 01 ξ = s i n h ξ when p 0 = 1 and q 0 = 1 and g 02 ξ = c o s h ξ when p 0 = 1 and q 0 = 1 . When n 1 , g n 1 ξ are the same as g n 2 ξ . So, we only discuss the case of g 02 ξ = c o s h ξ here, and discard the imaginary number solutions. When n = 1 , we get Equations (43) and (44).
f 11 ξ = sinh ξ 2 c o s h ξ + r μ = 1 / 4 , r = ± 1 g 1 ξ = c o s h 2 ξ 1 / 2 p 1 = 4 , q 1 = 1
v 11 ξ = k x 2 6 k x 2 · s i n h 2 ξ 2 · c o s h ξ + r 2 ξ = k x · x + k y · y + k y 2 + k x 4 / k x · t   and   w 11 ξ = k x 2 2 k x 2 · s i n h 2 ξ 2 · c o s h ξ + r 2 ξ = k x · x + k y · y + k y 2 k x 4 / k x · t
For the convivence of description, we assume k x = 1 , k y = 1 and r = 1 . Figure 1 shows bright solitary wave evolution of v 11 ξ (left figure) and w 11 ξ (right figure) vs. x , t . The structure of these two solitary waves is roughly the same, except for the amplitude and the position of the wave crest. The difference between v 11 ξ and w 11 ξ is caused by the constant terms and the coefficient of t , which will not affect the general structure of solitary waves.
When n = 2 , we get Equations (45) and (46).
f 12 ξ = sinh ξ · c o s h ξ c o s h 2 ξ 1 / 2 + r μ = 1 , r = ± 1 / 2 g 2 ξ = c o s h 4 ξ c o s h 2 ξ + 1 / 8 p 2 = 16 , q 2 = 1 / 4
v 12 ξ = 2 k x 2 3 2 k x 2 · s i n h 2 ξ · c o s h 2 ξ c o s h 2 ξ 1 / 2 + r 2 ξ = k x · x + k y · y + k y 2 + 4 k x 4 / k x · t   and   w 12 ξ = 2 k x 2 2 k x 2 · s i n h 2 ξ · c o s h 2 ξ c o s h 2 ξ 1 / 2 + r 2 ξ = k x · x + k y · y + k y 2 4 k x 4 / k x · t
Bright solitary wave evolutions of v 12 ξ (left figure) and w 12 ξ (right figure) are shown as Figure 2.
When n = 3 , we get Equations (47)–(49). Bright solitary wave evolutions of v 13 ξ (left figure) and w 13 ξ (right figure) are shown as Figure 3:
f 13 ξ = 2 · s i n h ξ · c o s h 3 ξ sinh ξ · c o s h ξ c o s h 4 ξ c o s h 2 ξ + 1 / 8 + r g 3 ξ = c o s h 8 ξ 2 · c o s h 6 ξ + 5 4 · c o s h 4 ξ 1 4 · c o s h 2 ξ + 1 128 μ = 4 , r = ± 1 8 , p 3 = 64 , q 2 = 1 256
v 13 ξ = 8 k x 2 3 2 k x 2 · 4 c o s h 8 ξ 8 c o s h 6 ξ + 5 c o s h 4 ξ c o s h 2 ξ c o s h 4 ξ c o s h 2 ξ + 1 8 + r 2 ξ = k x · x + k y · y + k y 2 + 16 k x 4 / k x · t
w 13 ξ = 8 k x 2 2 k x 2 · 4 c o s h 8 ξ 8 c o s h 6 ξ + 5 c o s h 4 ξ c o s h 2 ξ c o s h 4 ξ c o s h 2 ξ + 1 8 + r 2 ξ = k x · x + k y · y + k y 2 16 k x 4 / k x · t
When n = 4 , we get Equations (50)–(52). Solitary wave evolutions of v 14 ξ (left figure) and w 14 ξ (right figure) are shown in Figure 4.
f 14 ξ = 8 · s i n h ξ · c o s h 7 ξ 12 · s i n h ξ · c o s h 5 ξ + 5 · s i n h ξ · c o s h 3 ξ s i n h ξ · c o s h ξ / 2 2 · c o s h 8 ξ 2 · c o s h 6 ξ + 5 4 · c o s h 4 ξ 1 4 · c o s h 2 ξ + 1 128 + r g 4 ξ = c o s h 8 ξ 2 · c o s h 6 ξ + 5 4 · c o s h 4 ξ 1 4 · c o s h 2 ξ + 1 128 2 1 32768 μ = 16 , r = ± 1 128 , p 4 = 256 , q 4 = 1 4194304
v 14 ξ = 32 · k x 2 3 k x 2 · 32 c o s h 8 ξ 2 c o s h 6 ξ + 5 4 c o s h 4 ξ 1 4 c o s h 2 ξ 7 64 2 1 512 c o s h 8 ξ 2 c o s h 6 ξ + 5 4 c o s h 4 ξ 1 4 c o s h 2 ξ 7 64 + r 2 ξ = k x · x + k y · y + k y 2 + 64 k x 4 / k x · t
w 14 ξ = 32 · k x 2 k x 2 · 32 c o s h 8 ξ 2 c o s h 6 ξ + 5 4 c o s h 4 ξ 1 4 c o s h 2 ξ 7 64 2 1 512 c o s h 8 ξ 2 c o s h 6 ξ + 5 4 c o s h 4 ξ 1 4 c o s h 2 ξ 7 64 + r 2 ξ = k x · x + k y · y + k y 2 64 k x 4 / k x · t
With the increase of n in v 1 n ξ and w 1 n ξ , the amplitude of v 1 n ξ and w 1 n ξ increases, and the bandwidth in space–time decreases.

4.2. The Second Set of Solution

According to f 2 n ξ given by Equation (25), the second set of solution u 2 n ξ is given by Equation (53).
u 2 n ξ = 1 6 ± 1 6 ν 2 n 2 k x 2 + 2 λ k x 2 ν 2 n 2 [ s i n h ( ξ ) + ε · c o s h ξ ] ν n [ s i n h ( ξ ) + ε · c o s h ξ ] ν n + r 2 k x ν n [ s i n h ( ξ ) + ε · c o s h ξ ] ν n [ s i n h ( ξ ) + ε · c o s h ξ ] ν n + r 2 ξ = k x x + k y y + c t , c = k y 2 ν 2 n 2 k x 4 / k x , μ = ν 2 · n 2 4 , λ 2 = 1 , ε 2 = 1
Here, we only consider the case of λ = 1 , ν = 1 , ε = 1 and r = 1 . The expressions of v ξ and w ξ are shown as Equations (54) and (55):
ν 2 n ξ = 1 3 n 2 k x 2 + 2 k x 2 n 2 [ s i n h ( ξ ) + c o s h ξ ] n [ s i n h ( ξ ) + c o s h ξ ] n + 1 2 k x n [ s i n h ( ξ ) + c o s h ξ ] n [ s i n h ( ξ ) + c o s h ξ ] n + 1 2 ξ = k x · x + k y · y + k y 2 + n 2 k x 4 / k x · t
w 2 n ξ = 2 k x 2 n 2 [ s i n h ( ξ ) + c o s h ξ ] n [ s i n h ( ξ ) + c o s h ξ ] n + 1 2 k x n [ s i n h ( ξ ) + c o s h ξ ] n [ s i n h ( ξ ) + c o s h ξ ] n + 1 2 ξ = k x · x + k y · y + k y 2 n 2 k x 4 / k x · t
Figure 5, Figure 6, Figure 7 and Figure 8 show the structure of bright solitary wave of Equations (54) and (55) when n = 1 to 4. Under these conditions, with the increase of n , the amplitude of the bright solitary wave increases and the bandwidth in space–time decreases. The solitary wave structure is relatively stable under low- n conditions. Under high- n conditions, ν 2 n ξ and w 2 n ξ increase sharply in a small space–time range and becomes unstable.

4.3. The Third Set of Solutions

Here, we only consider the examples corresponding to g 0 ξ = c o s h ξ . According to the recurrence relation of Equation (32), we obtain the expressions of f 3 n ξ and g n ξ corresponding to different n . These expressions are brought into Equations (11) and (12) to find the two solutions v ξ and w ξ . The cases of n 4 are shown below.
When n = 1 , Equation (56) is the same as Equation (44), in which r = ± 1 . The solutions are shown as Figure 1, which will not be repeated here.
v 31 ξ = k x 2 6 k x 2 · s i n h 2 ξ 2 · c o s h ξ + r 2 ξ = k x · x + k y · y + k y 2 + k x 4 / k x · t   and   w 31 ξ = k x 2 2 k x 2 · s i n h 2 ξ 2 · c o s h ξ + r 2 ξ = k x · x + k y · y + k y 2 k x 4 / k x · t
When n = 2 , the solutions are shown as Equation (57), in which r = ± 2 2 . The corresponding solitary wave evolution is shown as Figure 9 with k x = 1 , k y = 1 and y = 0 .
v 32 ξ = k x 2 24 k x 2 · c o s h ξ 7 · k x 2 8 c o s h ξ + 1 + r 2 ξ = k x · x + k y · y + k y 2 + k x 4 4 / k x · t   and   w 32 ξ = k x 2 8 k x 2 · c o s h ξ 7 · k x 2 8 c o s h ξ + 1 + r 2 ξ = k x · x + k y · y + k y 2 k x 4 4 / k x · t
When n = 3 , the solutions are shown as Equation (58), in which r = ± 4 2 4 . The corresponding solitary wave evolution is shown as Figure 10 with k x = 1 , k y = 1 and y = 0 .
v 33 ξ = k x 2 96 k x 2 · c o s h ξ + 1 15 2 · k x 2 32 · c o s h ξ + 1 + 2 2 + r 2 ξ = k x · x + k y · y + k y 2 + k x 4 16 / k x · t   and   w 33 ξ = k x 2 32 k x 2 · c o s h ξ + 1 15 2 · k x 2 32 · c o s h ξ + 1 + 2 2 + r 2 ξ = k x · x + k y · y + k y 2 k x 4 16 / k x · t
When n = 4 , the solutions are shown as Equation (59), in which r = ± 4 2 2 4 . The corresponding solitary wave evolution is shown as Figure 11 with k x = 1 , k y = 1 and y = 0 .
v 34 ξ = k x 2 384 k x 2 · c o s h ξ + 1 + 2 2 7 2 4 · k x 2 128 · c o s h ξ + 1 + 2 2 + 4 2 4 + r 2 ξ = k x · x + k y · y + k y 2 + k x 4 64 / k x · t   and   w 34 ξ = k x 2 128 k x 2 · c o s h ξ + 1 + 2 2 7 2 4 · k x 2 128 · c o s h ξ + 1 + 2 2 + 4 2 4 + r 2 ξ = k x · x + k y · y + k y 2 k x 4 64 / k x · t
In this set of solutions, with the increase in n , the amplitude of the bright solitary wave decreases and the bandwidth in space–time increases. It would be easy to conclude that the solitary waves have relatively stable structures for different n .

4.4. The Fourth Set of Solutions

When n = 1 , the solutions are shown as Equation (60). The corresponding solitary wave evolution is shown as Figure 12 with k x = 1 , k y = 1 and y = 0 .
v 41 ξ = 2 · k x 2 3 2 · k x 2 · c o s h 2 ξ s i n h 2 ξ ξ = k x · x + k y · y + k y 2 + 4 k x 4 / k x · t   and   w 41 ξ = 2 · k x 2 2 · k x 2 · c o s h 2 ξ s i n h 2 ξ ξ = k x · x + k y · y + k y 2 4 k x 4 / k x · t
When n = 2 , the solutions are shown as Equation (61). The corresponding solitary wave evolution is shown as Figure 13 with k x = 1 , k y = 1 and y = 0 .
v 42 ξ = 8 · k x 2 3 32 · k x 2 · [ c o s h 2 ξ 1 2 ] 2 4 [ c o s h 2 ξ 1 2 ] 2 1 ξ = k x · x + k y · y + k y 2 + 16 k x 4 / k x · t   and   w 42 ξ = 8 · k x 2 32 · k x 2 · [ c o s h 2 ξ 1 2 ] 2 4 [ c o s h 2 ξ 1 2 ] 2 1 ξ = k x · x + k y · y + k y 2 16 k x 4 / k x · t
When n = 3 , the solutions are shown as Equation (62). The corresponding solitary wave evolution is shown as Figure 14 with k x = 1 , k y = 1 and y = 0 .
v 43 ξ = 32 k x 2 3 512 k x 2 · c o s h 4 ξ c o s h 2 ξ + 1 8 2 16 c o s h 4 ξ c o s h 2 ξ + 1 8 2 1 4 ξ = k x · x + k y · y + k y 2 + 64 k x 4 / k x · t and w 43 ξ = 32 k x 2 512 k x 2 · c o s h 4 ξ c o s h 2 ξ + 1 8 2 16 c o s h 4 ξ c o s h 2 ξ + 1 8 2 1 4 ξ = k x · x + k y · y + k y 2 64 k x 4 / k x · t
When n = 4 , the solutions are shown as Equation (63). The corresponding solitary wave evolution is shown as Figure 15 with k x = 1 , k y = 1 and y = 0 .
v 44 ξ = 128 · k x 2 3 8192 · k x 2 · c o s h 8 ξ 2 c o s h 6 ξ + 5 4 c o s h 4 ξ 1 4 c o s h 2 ξ 7 64 2 64 c o s h 8 ξ 2 c o s h 6 ξ + 5 4 c o s h 4 ξ 1 4 c o s h 2 ξ 7 64 2 1 256 ξ = k x · x + k y · y + k y 2 + 256 k x 4 / k x · t and w 44 ξ = 128 · k x 2 8192 · k x 2 · c o s h 8 ξ 2 c o s h 6 ξ + 5 4 c o s h 4 ξ 1 4 c o s h 2 ξ 7 64 2 64 c o s h 8 ξ 2 c o s h 6 ξ + 5 4 c o s h 4 ξ 1 4 c o s h 2 ξ 7 64 2 1 256 ξ = k x · x + k y · y + k y 2 256 k x 4 / k x · t
Figure 12, Figure 13, Figure 14 and Figure 15 show the dark solitary wave evolution for n = 1 to 4, respectively. The space–time width of the solitary wave decreases with the increase in n . However, under certain conditions, this set of infinite solutions has singularities in some specific space–time positions, which makes the solitary wave structure unstable.

4.5. The Fifth Set of Solutions

According to f 5 n ξ given by Equation (37), the second set of solution u 5 n ξ is given by Equation (64):
u 5 n ξ = 2 ± 2 3 ν 2 n 2 k x 2 + 8 λ k x 2 ν 2 n 2 [ s i n h ( ξ ) + ε · c o s h ξ ] 2 ν n [ s i n h ( ξ ) + ε · c o s h ξ ] 2 ν n + r 8 k x 2 ν n [ s i n h ( ξ ) + ε · c o s h ξ ] 2 ν n [ s i n h ( ξ ) + ε · c o s h ξ ] 2 ν n + r 2 ξ = k x x + k y y + c t , c = k y 2 4 ν 2 n 2 k x 4 / k x , λ 2 = 1 , ε 2 = 1
Here, we only consider the case of λ = 1 , ν = 1 , ε = 1 and r = 1 . The expressions of v ξ and w ξ are shown as Equation (65):
ν 5 n ξ = 4 3 n 2 k x 2 + 8 k x 2 n 2 [ s i n h ( ξ ) + c o s h ξ ] 2 n [ s i n h ( ξ ) + c o s h ξ ] 2 n + 1 8 k x 2 n [ s i n h ( ξ ) + c o s h ξ ] 2 n [ s i n h ( ξ ) + c o s h ξ ] 2 n + 1 2 ξ = k x · x + k y · y + k y 2 + 4 n 2 k x 4 / k x · t   and   w 5 n ξ = 8 k x 2 n 2 [ s i n h ( ξ ) + c o s h ξ ] 2 n [ s i n h ( ξ ) + c o s h ξ ] 2 n + 1 8 k x 2 n [ s i n h ( ξ ) + c o s h ξ ] 2 n [ s i n h ( ξ ) + c o s h ξ ] 2 n + 1 2 ξ = k x · x + k y · y + k y 2 4 n 2 k x 4 / k x · t
Respectively, Figure 16, Figure 17, Figure 18 and Figure 19 show the dark solitary wave evolution for n = 1 to 4, when k x = 1 , k y = 1 and y = 0 . While n increases, the amplitude of the solitary wave increases, and the bandwidth decreases.

4.6. The Sixth Set of Solutions

Here, we only give the cases of n 4 for u 6 n ξ for Equation (3). When n = 1 , the solutions are shown as Equation (66), in which r = ± 1 . The corresponding solitary wave evolution is shown as Figure 20 with k x = 1 , k y = 1 and y = 0 .
v 61 ξ = k x 2 6 k x 2 · c o s h ξ + r 2 2 · s i n h 2 ξ ξ = k x · x + k y · y + k y 2 + k x 4 / k x · t   and   w 61 ξ = k x 2 2 k x 2 · c o s h ξ + r 2 2 · s i n h 2 ξ ξ = k x · x + k y · y + k y 2 k x 4 / k x · t
When n = 2 , the solutions are shown as Equation (67), in which r = ± 1 2 . The corresponding solitary wave evolution is shown as Figure 21 with k x = 1 , k y = 1 and y = 0 .
v 62 ξ = 2 k x 2 3 8 k x 2 · [ c o s h 2 ξ 1 2 + r ] 2 4 [ c o s h 2 ξ 1 2 ] 2 1 ξ = k x · x + k y · y + k y 2 + 4 k x 4 / k x · t   and   w 62 ξ = 2 k x 2 8 k x 2 · [ c o s h 2 ξ 1 2 + r ] 2 4 [ c o s h 2 ξ 1 2 ] 2 1 ξ = k x · x + k y · y + k y 2 4 k x 4 / k x · t
When n = 3 , the solutions are shown as Equation (68), in which r = ± 1 8 . The corresponding solitary wave evolution is shown as Figure 22 with k x = 1 , k y = 1 and y = 0 .
v 63 ξ = 8 k x 2 3 128 · k x 2 · c o s h 4 ξ c o s h 2 ξ + 1 8 + r 2 4 c o s h 8 ξ 8 c o s h 6 ξ + 5 c o s h 4 ξ c o s h 2 ξ ξ = k x · x + k y · y + k y 2 + 16 k x 4 / k x · t   and   w 63 ξ = 8 k x 2 128 · k x 2 · c o s h 4 ξ c o s h 2 ξ + 1 8 + r 2 4 c o s h 8 ξ 8 c o s h 6 ξ + 5 c o s h 4 ξ c o s h 2 ξ ξ = k x · x + k y · y + k y 2 16 k x 4 / k x · t
When n = 4 , the solutions are shown as Equation (69), in which r = ± 1 128 . The corresponding solitary wave evolution is shown as Figure 23 with k x = 1 , k y = 1 and y = 0 .
v 64 ξ = 32 · k x 2 3 2048 · k x 2 · c o s h 8 ξ 2 c o s h 6 ξ + 5 4 c o s h 4 ξ 1 4 c o s h 2 ξ 7 64 + r 2 64 · c o s h 8 ξ 2 c o s h 6 ξ + 5 4 c o s h 4 ξ 1 4 c o s h 2 ξ 7 64 2 1 256 ξ = k x · x + k y · y + k y 2 + 64 k x 4 / k x · t   and   w 64 ξ = 32 · k x 2 2048 · k x 2 · c o s h 8 ξ 2 c o s h 6 ξ + 5 4 c o s h 4 ξ 1 4 c o s h 2 ξ 7 64 + r 2 64 · c o s h 8 ξ 2 c o s h 6 ξ + 5 4 c o s h 4 ξ 1 4 c o s h 2 ξ 7 64 2 1 256 ξ = k x · x + k y · y + k y 2 64 k x 4 / k x · t
Figure 20, Figure 21 and Figure 22 represent the dark solitary wave solutions of the KPI equation. Due to the existence of singularities, the solitary wave structure of each solution is completely unstable.

4.7. The Seventh Set of Solutions

According to f 7 n ξ given by Equation (39), the seventh set of solution u 7 n ξ is given by Equation (70).
u 7 n ξ = 1 3 ± 1 6 ν 2 n 2 k 2 2 k 2 ν n [ s i n h ( ξ ) + ε · c o s h ξ ] ν n + ν n r 2 λ + 4 ε [ s i n h ( ξ ) + ε · c o s h ξ ] ν n + 2 λ r 2 ξ = k x x + k y y + c t , c = k y 2 ν 2 n 2 k x 4 / k x , λ 2 = 1 , ε 2 = 1 .
Here, we only consider the case of λ = 1 , ν = 1 , ε = 1 and r = 1 . The expressions of v ξ and w ξ are shown as Equation (71).
ν 7 n ξ = n 2 k x 2 6 2 k x 2 n [ s i n h ( ξ ) + c o s h ξ ] n + n 6 · [ s i n h ( ξ ) + c o s h ξ ] n + 2 2 ξ = k x · x + k y · y + k y 2 + n 2 k x 4 / k x · t   and   w 7 n ξ = n 2 k x 2 2 2 k x 2 n [ s i n h ( ξ ) + c o s h ξ ] n + n 6 · [ s i n h ( ξ ) + c o s h ξ ] n + 2 2 ξ = k x · x + k y · y + k y 2 n 2 k x 4 / k x · t
The solutions of n = 1   to   4 are shown as Figure 24, Figure 25, Figure 26 and Figure 27, respectively. We find that the structures of these solutions have an abrupt change in their amplitude, which maybe implies a bistable mode, including those with fractional derivatives [45], which is of particular interest in optical technology, power networks, and so on [46,47]. Different values of n can make the amplitude, space–time width, kink direction, and singularity position of the solitary wave different. In this case, all the solitary waves have relatively stable structures.

4.8. The Eighth Set of Solutions

The cases of n 4 for u 8 n ξ for Equation (3) is given below. When n = 1 , the solutions are shown as Equation (72), in which r = ± 1 . These solutions are to the same as v 61 and w 61 given in Section 4.6, which will not be repeated here.
v 81 ξ = k x 2 6 k x 2 · c o s h ξ + r 2 2 · s i n h 2 ξ ξ = k x · x + k y · y + k y 2 + k x 4 / k x · t   and   w 81 ξ = k x 2 2 k x 2 · c o s h ξ + r 2 2 · s i n h 2 ξ ξ = k x · x + k y · y + k y 2 k x 4 / k x · t
When n = 2 , the solutions are shown as Equation (73), in which r = ± 2 2 . The corresponding solitary wave evolution is shown as Figure 28 with k x = 1 , k y = 1 and y = 0 .
v 82 ξ = k x 2 24 k x 2 · c o s h ξ + 1 + r 2 8 · c o s h ξ 56 ξ = k x · x + k y · y + k y 2 + k x 4 4 / k x · t   and   w 82 ξ = k x 2 8 k x 2 · c o s h ξ + 1 + r 2 8 · c o s h ξ 56 ξ = k x · x + k y · y + k y 2 k x 4 4 / k x · t
When n = 3 , the solutions are shown as Equation (74), in which r = ± 4 2 4 . The corresponding solitary wave evolution is shown as Figure 29 with k x = 1 , k y = 1 and y = 0 .
v 83 ξ = k x 2 96 k x 2 · c o s h ξ + 1 + 2 2 + r 2 32 · c o s h ξ + 1 480 2 ξ = k x · x + k y · y + k y 2 + k x 4 16 / k x · t   and   w 83 ξ = k x 2 32 k x 2 · c o s h ξ + 1 + 2 2 + r 2 32 · c o s h ξ + 1 480 2 ξ = k x · x + k y · y + k y 2 k x 4 16 / k x · t
When n = 4 , the solutions are shown as Equation (75), in which r = ± 4 2 2 4 . The corresponding solitary wave evolution is shown as Figure 30 with k x = 1 , k y = 1 and y = 0 .
v 84 ξ = k x 2 384 k x 2 · c o s h ξ + 1 + 2 2 + 4 2 4 + r 2 128 · c o s h ξ + 1 + 2 2 3584 2 4 ξ = k x · x + k y · y + k y 2 + k x 4 64 / k x · t   and   w 84 ξ = k x 2 128 k x 2 · c o s h ξ + 1 + 2 2 + 4 2 4 + r 2 128 · c o s h ξ + 1 + 2 2 3584 2 4 ξ = k x · x + k y · y + k y 2 k x 4 64 / k x · t
It can be seen that the shape of the solitary wave and the space–time position of singular points will change with the value of n . v 81 ξ and w 81 ξ represent dark solitary waves. With the increase in n , it gradually evolves into a traveling wave with two rows of singular points. The solitary wave structure becomes unstable under high- n conditions.

4.9. The Ninth Set of Solutions

When n = 1 , the solutions are shown as Equation (76). The corresponding solitary wave evolution is shown as Figure 31 with k x = 1 , k y = 1 and y = 0 .
v 91 ξ = 2 k x 2 3 2 k x 2 · s i n h 2 ξ c o s h 2 ξ ξ = k x · x + k y · y + k y 2 + 4 k x 4 / k x · t   and   w 91 ξ = 2 k x 2 2 k x 2 · s i n h 2 ξ c o s h 2 ξ ξ = k x · x + k y · y + k y 2 4 k x 4 / k x · t
When n = 2 , the solutions are shown as Equation (77). The corresponding solitary wave evolution is shown as Figure 32 with k x = 1 , k y = 1 and y = 0 .
v 92 ξ = 8 k x 2 3 2 k x 2 · 4 · [ c o s h 2 ξ 1 2 ] 2 1 [ c o s h 2 ξ 1 2 ] 2 ξ = k x · x + k y · y + k y 2 + 16 k x 4 / k x · t   and   w 92 ξ = 8 k x 2 2 k x 2 · 4 · [ c o s h 2 ξ 1 2 ] 2 1 [ c o s h 2 ξ 1 2 ] 2 ξ = k x · x + k y · y + k y 2 16 k x 4 / k x · t
When n = 3 , the solutions are shown as Equation (78). The corresponding solitary wave evolution is shown as Figure 33 with k x = 1 , k y = 1 and y = 0 .
v 93 ξ = 32 k x 2 3 2 k x 2 16 c o s h 2 ξ 1 2 2 1 8 2 1 4 c o s h 2 ξ 1 2 2 1 8 2 ξ = k x · x + k y · y + k y 2 + 64 k x 4 / k x · t   and w 93 ξ = 32 k x 2 2 k x 2 16 c o s h 2 ξ 1 2 2 1 8 2 1 4 c o s h 2 ξ 1 2 2 1 8 2 ξ = k x · x + k y · y + k y 2 64 k x 4 / k x · t
When n = 4 , the solutions are shown as Equation (79). The corresponding solitary wave evolution is shown as Figure 34 with k x = 1 , k y = 1 and y = 0 .
v 94 ξ = 128 k x 2 3 2 k x 2 · 64 c o s h 8 ξ 2 c o s h 6 ξ + 5 4 c o s h 4 ξ 1 4 c o s h 2 ξ 7 64 2 1 256 c o s h 8 ξ 2 c o s h 6 ξ + 5 4 c o s h 4 ξ 1 4 c o s h 2 ξ 7 64 2 ξ = k x · x + k y · y + k y 2 + 256 k x 4 / k x · t   and   w 94 ξ = 128 k x 2 2 k x 2 · 64 c o s h 8 ξ 2 c o s h 6 ξ + 5 4 c o s h 4 ξ 1 4 c o s h 2 ξ 7 64 2 1 256 c o s h 8 ξ 2 c o s h 6 ξ + 5 4 c o s h 4 ξ 1 4 c o s h 2 ξ 7 64 2 ξ = k x · x + k y · y + k y 2 256 k x 4 / k x · t
In this case, the increase in n reduces the space–time width of the solitary wave and changes the space–time position of the singularity. v 91 ξ and w 91 ξ represent bright solitary waves. With the increase in n , it gradually evolves into a traveling wave with many singular points. In this case, stable solitary wave structures only exist under low- n conditions, while they become unstable under high- n conditions.

4.10. The Tenth Set of Solutions

According to f X n ξ given by Equation (42), the second set of solution u X n ξ is given by Equation (80).
u X n ξ = 4 3 ± 2 3 ν 2 n 2 k 2 2 k 2 ν n [ s i n h ( ξ ) + ε · c o s h ξ ] 2 ν n + ν n r λ 2 ε [ s i n h ( ξ ) + ε · c o s h ξ ] 2 ν n ± r 2 ξ = k x x + k y y + c t , c = k y 2 4 ν 2 n 2 k x 4 / k x , μ = ν 2 · n 2 4 , λ 2 = 1 , ε 2 = 1 .
Here, we only consider the cases of λ = 1 , ν = 1 , ε = 1 and r = 1 . The expressions of v ξ and w ξ are shown as Equation (81).
ν X n ξ = 2 3 n 2 k x 2 2 k 2 n [ s i n h ( ξ ) + c o s h ξ ] 2 n + n [ s i n h ( ξ ) + c o s h ξ ] 2 n + 1 2 ξ = k x · x + k y · y + k y 2 + 4 n 2 k x 4 / k x · t   and   w X n ξ = 2 n 2 k x 2 2 k 2 n [ s i n h ( ξ ) + c o s h ξ ] 2 n + n [ s i n h ( ξ ) + c o s h ξ ] 2 n + 1 2 ξ = k x · x + k y · y + k y 2 4 n 2 k x 4 / k x · t
Figure 35, Figure 36, Figure 37 and Figure 38 represent the 10th set of solitary wave solutions of the KPI equation. The increase in n can change the amplitude, space–time width, singularity position, and propagation direction of the solitary wave. In this case, under low- n conditions, there are very stable solitary wave structures, while under high- n conditions, the wave structures become unstable.

5. Summary

In nature, most nonlinear phenomena can be expressed as corresponding mathematical models under specific approximate conditions, and be simplified as NLEEs varying with time and space. In the study of such processes, the accuracy and stability of the solution is important. Many researchers have developed various methods to obtain the analytical solutions of NLEEs. The KPI equation is a well-known NLEE. In recent years, its applications in nonlinear optics, plasma ion acoustic mechanics, and other fields have made significant progress. Acquiring precise answers to the KPI equation is a crucial research tool in these fields.
In this paper, a general Riccati equation is treated as an auxiliary equation. Through different function transformations of f ξ shown as Equation (13), (33) and a 1 / f ξ , ten sets of infinite-series solitary wave solutions of the KPI equation have been constructed via the auxiliary equation method, which would be helpful to study the nonlinear problems described by the KPI equation. Unlike some other methods which can only provide a limited number of solutions, we can construct infinite-series exact solitary wave solutions with a relatively simple method.
The images of the first four solutions of each set have been plotted to provide intuitive structure of waves. Among these infinite-series exact wave solutions of the KPI equation, we find different structures of waves, including bright solitary waves, dark solitary waves, kink solitary waves and, traveling waves with singular points. Some of these solutions are stable, and the others are relatively unstable. A simple and intuitive introduction is given in Section 4. These results not only include known solutions, but also put forward some new forms of solutions, such as kink mode and traveling wave with many singular points, which may imply a new evolution mode. The detailed analysis should be combined with specific research objects, which will be discussed in future research papers. The results show that this method is simple and effective for the construction of infinite series solutions of nonlinear evolution models.

Author Contributions

Conceptualization, G.W. and Y.G.; methodology, G.W.; software, F.P.; validation, G.W. and Y.G.; formal analysis, F.P.; investigation, F.P.; resources, G.W. and Y.G.; data curation, F.P.; writing—original draft preparation, F.P.; writing—review and editing, G.W.; visualization, F.P.; supervision, G.W. and Y.G.; project administration, G.W. and Y.G.; funding acquisition, Y.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by National Natural Science Foundation of China with Contract Nos. 12275307, 11575238, 11475222, 11505228, 10975159, 11275232, and 11735016, the National MCF Energy R&D Program (2019YFE03040000, 2019YFE03030000), the Research Project of Huainan Normal University (2020XJZD011).

Data Availability Statement

No data, models, or code were generated or used during the study.

Acknowledgments

Most numerical calculations in this paper have been done on the supercomputing system in the Supercomputing Center of University of Science and Technology of China, and the ShenMa High Performance Computing Cluster in Institute of Plasma Physics, Chinese Academy of Sciences.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Guo, H.D.; Xia, T.C.; Hu, B.B. High-order lumps, high-order breathers and hybrid solutions for an extended (3+1)-dimensional Jimbo-Miwa equation in fluid dynamics. Nonlinear Dyn. 2020, 100, 601–614. [Google Scholar] [CrossRef]
  2. Lan, Z.Z.; Guo, B.L. Nonlinear waves behaviors for a coupled generalized nonlinear Schrodinger-Boussinesq system in a homogeneous magnetized plasma. Nonlinear Dyn. 2020, 100, 3771–3784. [Google Scholar] [CrossRef]
  3. Biswas, A.; Ekici, M.; Sonmezoglu, A.; Belic, M.R. Solitons in optical fiber Bragg gratings with dispersive reflectivity by extended trial function method. Optik 2019, 182, 88–94. [Google Scholar] [CrossRef]
  4. Seadawy, A.R.; Lu, D.C.; Nasreen, N.; Nasreen, S. Structure of optical solitons of resonant Schrodinger equation with quadratic cubic nonlinearity and modulation instability analysis. Phys. A Stat. Mech. Its Appl. 2019, 534, 122155. [Google Scholar] [CrossRef]
  5. Abdoud, M.A.; Owyed, S.; Abdel-Aty, A.; Raffan, B.M.; Abdel-Khalek, S. Optical soliton solutions for a space-time fractional perturbed nonlinear Schrödinger equation arising in quantum physics. Results Phys. 2020, 16, 102895. [Google Scholar] [CrossRef]
  6. Peng, W.Q.; Tian, S.F.; Zhang, T.T. Dynamics of the soliton waves, breather waves, and rogue waves to the cylindrical Kadomtsev-Petviashvili equation in pair-ion-electron plasma. Phys. Fluids 2019, 31, 102107. [Google Scholar] [CrossRef]
  7. Liu, J.B.; Yang, K.Q. The extended F-expansion method and exact solutions of nonlinear PDEs. Chaos Solitons Fractals 2004, 22, 111. [Google Scholar] [CrossRef]
  8. Zhang, S. Application of Exp-function method to a KdV equation with variable coefficients. Phys. Lett. A 2007, 365, 448. [Google Scholar] [CrossRef]
  9. Zhang, S. The periodic wave solutions for the (2+1)-dimensional Konopelchenko–Dubrovsky equations. Chaos Solitons Fractals 2006, 30, 1213. [Google Scholar] [CrossRef]
  10. Zhang, S. The periodic wave solutions for the (2+1)-dimensional dispersive long water equations. Chaos Solitons Fractals 2007, 32, 847. [Google Scholar] [CrossRef]
  11. Abdou, M.A. The extended F-expansion method and its application for a class of nonlinear evolution equations. Chaos Solitons Fractals 2007, 31, 95. [Google Scholar] [CrossRef]
  12. Wazwaz, A.M. The extended tanh method for new solitons solutions for many forms of the fifth-order KdV equations. Appl. Math. Comput. 2007, 184, 1002–1014. [Google Scholar] [CrossRef]
  13. Wazwaz, A.M. The tanh–coth method for solitons and kink solutions for nonlinear parabolic equations. Appl. Math. Comput. 2007, 188, 1467–1475. [Google Scholar] [CrossRef]
  14. Fan, E. Extended tanh-function method and its applications to nonlinear equations. Phys. Lett. A 2000, 277, 212. [Google Scholar] [CrossRef]
  15. Liu, S.K.; Fu, Z.T.; Liu, S.D.; Zhao, Q. Jacobi elliptic function expansion method and periodic wave solutions of nonlinear wave equations. Phys. Lett. A 2001, 289, 69. [Google Scholar] [CrossRef]
  16. Liu, S.K.; Fu, Z.T.; Liu, S.D.; Zhao, Q. New Jacobi elliptic function expansion and new periodic solutions of nonlinear wave equations. Phys. Lett. A 2001, 290, 72. [Google Scholar] [CrossRef]
  17. Elgarayhi, A. New periodic wave solutions for the shallow water equations and the generalized Klein–Gordon equation. Commun. Nonlinear Sci. Numer. Simul. 2008, 13, 877–888. [Google Scholar] [CrossRef]
  18. Wu, G.; Guo, Y. New Complex Wave Solutions and Diverse Wave Structures of the (2+1)-Dimensional Asymmetric Nizhnik–Novikov–Veselov Equation. Fractal Fract. 2023, 7, 170. [Google Scholar] [CrossRef]
  19. Wu, G.; Han, J.; Zhang, W.; Zhang, M. New periodic wave solutions to nonlinear evolution equations by the extended mapping method. Physical D-Nonlinear Phenom. 2007, 229, 116. [Google Scholar] [CrossRef]
  20. Sirendaoreji; Jiong, S. Auxiliary equation method for solving nonlinear partial differential equations. Phys. Lett. A 2003, 309, 387. [Google Scholar] [CrossRef]
  21. Sirendaoreji. New exact travelling wave solutions for the Kawahara and modified Kawahara equations. Chaos Solitons Fractals 2004, 19, 147. [Google Scholar] [CrossRef]
  22. Zhu, X.; Cheng, J.; Chen, Z.; Wu, G. New Solitary-Wave Solutions of the Van der Waals Normal Form for Granular Materials via New Auxiliary Equation Method. Mathematics 2022, 10, 2560. [Google Scholar] [CrossRef]
  23. Wu, G.; Guo, Y. Construction of New Infinite-Series Exact Solitary Wave Solutions and Its Application to the Korteweg–De Vries Equation. Fractal Fract. 2023, 7, 75. [Google Scholar] [CrossRef]
  24. Wang, M.L. Solitary wave solutions for variant Boussinesq equations. Phys. Lett. A 1995, 199, 169–172. [Google Scholar] [CrossRef]
  25. Wang, M.L.; Zhou, Y.B.; Li, Z.B. Application of a homogenous balance method to exact solitons of nonlinear equations in mathematical physics. Phys. Lett. A 1996, 216, 67–75. [Google Scholar] [CrossRef]
  26. Otwinowski, M.; Paul, R. Laidlaw WG. Exact travelling wave solutions of a class of nonlinear diffusion equations by reduction to a quadrature. Phys. Lett. A 1988, 128, 483–487. [Google Scholar] [CrossRef]
  27. Wang, M.; Li, X.; Zhang, J. The (G′/G)-expansion method and travelling wave solutions of nonlinear evolution equations in mathematical physics. Phys. Lett. A 2008, 372, 417. [Google Scholar] [CrossRef]
  28. Zayed, E.M.E.; Gepreel, K.A. The G′/G-expansion method for finding the traveling wave solutions of nonlinear partial differential equations in mathematical physics. J. Math. Phys. 2009, 50, 013502. [Google Scholar] [CrossRef]
  29. Guo, S.; Zhou, Y. The extended G′/G -expansion method and its applications to the Whitham–Broer–Kaup–Like equations and coupled Hirota–Satsuma KdV equations. Appl. Math. Comput. 2010, 215, 3214. [Google Scholar] [CrossRef]
  30. Islam, M.S.; Khan, K.; Akbar, M.A. An analytical method for finding exact solutions of modified Korteweg–de Vries equation. Results Phys. 2015, 5, 131. [Google Scholar] [CrossRef] [Green Version]
  31. Kadomtsev, B.B.; Petviashvili, V.I. On the stability of solitary waves in weakly dispersing media. Sov. Phys. Dokl. 1970, 15, 539–554. [Google Scholar]
  32. Segur, H.; Finkel, A. An analytical model of periodic waves in shallow water. Stud. Appl. Math. 1985, 73, 183–220. [Google Scholar] [CrossRef]
  33. Hammack, J.; Scheffner, N.; Segur, H. Two-dimensional periodic waves in shallow water. J. Fluid Mech. 1989, 209, 567–589. [Google Scholar] [CrossRef]
  34. Hammack, J.; Scheffner, N.; Segur, H. Two-dimensional periodic waves in shallow water. Part 2. Asymmetric waves. J. Fluid Mech. 1995, 285, 95–122. [Google Scholar] [CrossRef]
  35. Infeld, E.; Rowlands, G. Nonlinear Waves, Solitons and Chaos; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  36. Pelinovsky, D.E.; Stepanyants, Y.A.; Kivshar, Y.S. Self-focusing of plane dark solitons in nonlinear defocusing media. Phys. Rev. E 1995, 51, 5016–5026. [Google Scholar] [CrossRef] [PubMed]
  37. Baronio, F.; Wabnitz, S.; Kodama, Y. Optical Kerr spatio-temporal dark-lump dynamics of hydrodynamic origin. Phys. Rev. Lett. 2016, 116, 173901. [Google Scholar] [CrossRef] [Green Version]
  38. Tajiri, M.; Murakami, Y. Two-Dimensional Multisoliton Solutions: Periodic Soliton Solutions to the Kadomtsev-Petviashvili Equation with Positive Dispersion. J. Phys. Soc. Jpn. 1989, 58, 3029–3032. [Google Scholar] [CrossRef]
  39. Manakov, S.; Zakharov, V.; Bordag, L.; Its, A.; Matveev, V. Two-dimensional solitons of the Kadomtsev-Petviashvili equation and their interaction. Phys. Lett. A 1977, 63, 205–206. [Google Scholar] [CrossRef]
  40. Johnson, R.; Thompson, S. A solution of the inverse scattering problem for the Kadomtsev-Petviashvili equation by the method of separation of variables. Phys. Lett. A 1978, 66, 279–281. [Google Scholar] [CrossRef]
  41. Tajiri, M.; Fujimura, Y.; Murakami, Y. Resonant Interactions between Y-Periodic Soliton and Algebraic Soliton: Solutions to the Kadomtsev-Petviashvili Equation with Positive Dispersion. J. Phys. Soc. Jpn. 1992, 61, 783–790. [Google Scholar] [CrossRef]
  42. Murakami, Y.; Tajiri, M. Resonant Interaction between Line Soliton and Y-Periodic Soliton: Solutions to the Kadomtsev-Petviashvili Equation with Positive Dispersion. J. Phys. Soc. Jpn. 1992, 61, 791–805. [Google Scholar] [CrossRef]
  43. Klein, C.; Saut, J.-C. Numerical Study of Blow up and Stability of Solutions of Generalized Kadomtsev–Petviashvili Equations. J. Nonlinear Sci. 2012, 22, 763–811. [Google Scholar] [CrossRef] [Green Version]
  44. Chakravarty, S.; Zowada, M. Multi-lump wave patterns of KPI via integer partitions. Phys. D Nonlinear Phenom. 2023, 446, 133644. [Google Scholar] [CrossRef]
  45. Kim, V.A.; Parovik, R.I.; Rakhmonov, Z.R. Implicit Finite-Difference Scheme for a Duffing Oscillator with a Derivative of Variable Fractional Order of the Riemann-Liouville Type. Mathematics 2023, 11, 558. [Google Scholar] [CrossRef]
  46. Loi, W.S.; Ong, C.T. Optical solitons: Mathematical model and simulations. In Proceedings of the 2013 IEEE Business Engineering and Industrial Applications Colloquium (BEIAC), Langkawi, Malaysia, 7–9 April 2013; 2013; pp. 901–905. [Google Scholar] [CrossRef]
  47. Fedoruk, M.P. Mathematical modeling of optical communication lines with dispersion management. In Computational Science and High Performance Computing. Notes on Numerical Fluid Mechanics and Multidisciplinary Design (NNFM); Krause, E., Shokin, Y.I., Resch, M., Shokina, N., Eds.; Springer: Berlin/Heidelberg, Germany, 2005; Volume 88. [Google Scholar] [CrossRef]
Figure 1. Bright solitary wave evolution of v 11 ξ (left figure) and w 11 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 and y = 0 .
Figure 1. Bright solitary wave evolution of v 11 ξ (left figure) and w 11 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 and y = 0 .
Mathematics 11 01560 g001
Figure 2. Bright solitary wave evolution of v 12 ξ (left figure) and w 12 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 / 2 and y = 0 .
Figure 2. Bright solitary wave evolution of v 12 ξ (left figure) and w 12 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 / 2 and y = 0 .
Mathematics 11 01560 g002
Figure 3. Bright solitary wave evolution of v 13 ξ (left figure) and w 13 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 / 8 and y = 0 .
Figure 3. Bright solitary wave evolution of v 13 ξ (left figure) and w 13 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 / 8 and y = 0 .
Mathematics 11 01560 g003
Figure 4. Solitary wave evolution of v 14 ξ (left figure) and w 14 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 / 128 and y = 0 .
Figure 4. Solitary wave evolution of v 14 ξ (left figure) and w 14 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 / 128 and y = 0 .
Mathematics 11 01560 g004
Figure 5. Bright solitary wave evolution of v 21 ξ (left figure) and w 21 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 and y = 0 .
Figure 5. Bright solitary wave evolution of v 21 ξ (left figure) and w 21 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 and y = 0 .
Mathematics 11 01560 g005
Figure 6. Bright solitary wave evolution of v 22 ξ (left figure) and w 22 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 and y = 0 .
Figure 6. Bright solitary wave evolution of v 22 ξ (left figure) and w 22 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 and y = 0 .
Mathematics 11 01560 g006
Figure 7. Bright solitary wave evolution of v 23 ξ (left figure) and w 23 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 and y = 0 .
Figure 7. Bright solitary wave evolution of v 23 ξ (left figure) and w 23 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 and y = 0 .
Mathematics 11 01560 g007
Figure 8. Bright solitary wave evolution of v 24 ξ (left figure) and w 24 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 and y = 0 .
Figure 8. Bright solitary wave evolution of v 24 ξ (left figure) and w 24 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 and y = 0 .
Mathematics 11 01560 g008
Figure 9. Bright solitary wave evolution of v 32 ξ (left figure) and w 32 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = ± 2 2 and y = 0 .
Figure 9. Bright solitary wave evolution of v 32 ξ (left figure) and w 32 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = ± 2 2 and y = 0 .
Mathematics 11 01560 g009
Figure 10. Bright solitary wave evolution of v 33 ξ (left figure) and w 33 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = ± 4 2 4 and y = 0 .
Figure 10. Bright solitary wave evolution of v 33 ξ (left figure) and w 33 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = ± 4 2 4 and y = 0 .
Mathematics 11 01560 g010
Figure 11. Bright solitary wave evolution of v 34 ξ (left figure) and w 34 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = ± 4 2 2 4 and y = 0 .
Figure 11. Bright solitary wave evolution of v 34 ξ (left figure) and w 34 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = ± 4 2 2 4 and y = 0 .
Mathematics 11 01560 g011
Figure 12. Dark solitary wave evolution of v 41 ξ (left figure) and w 41 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 12. Dark solitary wave evolution of v 41 ξ (left figure) and w 41 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g012
Figure 13. Dark solitary wave evolution of v 42 ξ (left figure) and w 42 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 13. Dark solitary wave evolution of v 42 ξ (left figure) and w 42 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g013
Figure 14. Dark solitary wave evolution of v 43 ξ (left figure) and w 43 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 14. Dark solitary wave evolution of v 43 ξ (left figure) and w 43 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g014
Figure 15. Solitary wave evolution of v 44 ξ (left figure) and w 44 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 15. Solitary wave evolution of v 44 ξ (left figure) and w 44 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g015
Figure 16. Wave evolution of v 51 ξ (left figure) and w 51 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 16. Wave evolution of v 51 ξ (left figure) and w 51 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g016
Figure 17. Wave evolution of v 52 ξ (left figure) and w 52 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 17. Wave evolution of v 52 ξ (left figure) and w 52 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g017
Figure 18. Wave evolution of v 53 ξ (left figure) and w 53 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 18. Wave evolution of v 53 ξ (left figure) and w 53 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g018
Figure 19. Wave evolution of v 54 ξ (left figure) and w 54 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 19. Wave evolution of v 54 ξ (left figure) and w 54 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g019
Figure 20. Dark solitary wave evolution of v 61 ξ (left figure) and w 61 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 and y = 0 .
Figure 20. Dark solitary wave evolution of v 61 ξ (left figure) and w 61 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 and y = 0 .
Mathematics 11 01560 g020
Figure 21. Dark solitary wave evolution of v 62 ξ (left figure) and w 62 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 / 2 and y = 0 .
Figure 21. Dark solitary wave evolution of v 62 ξ (left figure) and w 62 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 / 2 and y = 0 .
Mathematics 11 01560 g021
Figure 22. Dark solitary wave evolution of v 63 ξ (left figure) and w 63 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 / 8 and y = 0 .
Figure 22. Dark solitary wave evolution of v 63 ξ (left figure) and w 63 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 / 8 and y = 0 .
Mathematics 11 01560 g022
Figure 23. Solitary wave evolution of v 64 ξ (left figure) and w 64 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 / 128 and y = 0 .
Figure 23. Solitary wave evolution of v 64 ξ (left figure) and w 64 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 1 / 128 and y = 0 .
Mathematics 11 01560 g023
Figure 24. Wave evolution of v 71 ξ (left figure) and w 71 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 24. Wave evolution of v 71 ξ (left figure) and w 71 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g024
Figure 25. Wave evolution of v 72 ξ (left figure) and w 72 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 25. Wave evolution of v 72 ξ (left figure) and w 72 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g025
Figure 26. Wave evolution of v 73 ξ (left figure) and w 73 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 26. Wave evolution of v 73 ξ (left figure) and w 73 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g026
Figure 27. Wave evolution of v 74 ξ (left figure) and w 74 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 27. Wave evolution of v 74 ξ (left figure) and w 74 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g027
Figure 28. Solitary wave evolution of v 82 ξ (left figure) and w 82 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 2 2 and y = 0 .
Figure 28. Solitary wave evolution of v 82 ξ (left figure) and w 82 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 2 2 and y = 0 .
Mathematics 11 01560 g028
Figure 29. Solitary wave evolution of v 83 ξ (left figure) and w 83 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 4 2 4 and y = 0 .
Figure 29. Solitary wave evolution of v 83 ξ (left figure) and w 83 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 4 2 4 and y = 0 .
Mathematics 11 01560 g029
Figure 30. Solitary wave evolution of v 84 ξ (left figure) and w 84 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 4 2 2 4 and y = 0 .
Figure 30. Solitary wave evolution of v 84 ξ (left figure) and w 84 ξ (right figure) vs. x , t when k x = 1 , k y = 1 , r = 4 2 2 4 and y = 0 .
Mathematics 11 01560 g030
Figure 31. Wave evolution of v 91 ξ (left figure) and w 91 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 31. Wave evolution of v 91 ξ (left figure) and w 91 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g031
Figure 32. Wave evolution of v 92 ξ (left figure) and w 92 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 32. Wave evolution of v 92 ξ (left figure) and w 92 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g032
Figure 33. Wave evolution of v 93 ξ (left figure) and w 93 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 33. Wave evolution of v 93 ξ (left figure) and w 93 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g033
Figure 34. Wave evolution of v 94 ξ (left figure) and w 94 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 34. Wave evolution of v 94 ξ (left figure) and w 94 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g034
Figure 35. Wave evolution of v X 1 ξ (left figure) and w X 1 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 35. Wave evolution of v X 1 ξ (left figure) and w X 1 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g035
Figure 36. Wave evolution of v X 2 ξ (left figure) and w X 2 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 36. Wave evolution of v X 2 ξ (left figure) and w X 2 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g036
Figure 37. Wave evolution of v X 3 ξ (left figure) and w X 3 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 37. Wave evolution of v X 3 ξ (left figure) and w X 3 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g037
Figure 38. Wave evolution of v X 4 ξ (left figure) and w X 4 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Figure 38. Wave evolution of v X 4 ξ (left figure) and w X 4 ξ (right figure) vs. x , t when k x = 1 , k y = 1 and y = 0 .
Mathematics 11 01560 g038
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pei, F.; Wu, G.; Guo, Y. Construction of Infinite Series Exact Solitary Wave Solution of the KPI Equation via an Auxiliary Equation Method. Mathematics 2023, 11, 1560. https://doi.org/10.3390/math11061560

AMA Style

Pei F, Wu G, Guo Y. Construction of Infinite Series Exact Solitary Wave Solution of the KPI Equation via an Auxiliary Equation Method. Mathematics. 2023; 11(6):1560. https://doi.org/10.3390/math11061560

Chicago/Turabian Style

Pei, Feiyun, Guojiang Wu, and Yong Guo. 2023. "Construction of Infinite Series Exact Solitary Wave Solution of the KPI Equation via an Auxiliary Equation Method" Mathematics 11, no. 6: 1560. https://doi.org/10.3390/math11061560

APA Style

Pei, F., Wu, G., & Guo, Y. (2023). Construction of Infinite Series Exact Solitary Wave Solution of the KPI Equation via an Auxiliary Equation Method. Mathematics, 11(6), 1560. https://doi.org/10.3390/math11061560

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