# Line Soliton Interactions for Shallow Ocean Waves and Novel Solutions with Peakon, Ring, Conical, Columnar, and Lump Structures Based on Fractional KP Equation

Authors: Bo Xu; Yufeng Zhang; Sheng Zhang
Journal: Advances in Mathematical Physics
 (2021)

---

## Abstract

It is well known that the celebrated Kadomtsev-Petviashvili (KP) equation has many important applications. The aim of this article is to use fractional KP equation to not only simulate shallow ocean waves but also construct novel spatial structures. Firstly, the definitions of the conformable fractional partial derivatives and integrals together with a physical interpretation are introduced and then a fractional integrable KP equation consisting of fractional KPI and KPII equations is derived. Secondly, a formula for the fractionaln-soliton solutions of the derived fractional KP equation is obtained and fractional line one-solitons with bend, wavelet peaks, and peakon are constructed. Thirdly, fractional X-, Y- and 3-in-2-out-type interactions in the fractional line two- and three-soliton solutions of the fractional KPII equation are simulated for shallow ocean waves. Besides, a falling and spreading process of a columnar structure in the fractional line two-soliton solution is also simulated. Finally, a fractional rational solution of the fractional KP equation is obtained including the lump solution as a special case. With the development of time, the nonlinear dynamic evolution of the fractional lump solution of the fractional KPI equation can change from ring and conical structures to lump structure.


---

## Body

## 1. Introduction

As Ablowitz and Baldwin [1] pointed, the two-space and one-time dimensional equation governing unidirectional, maximally balanced, weakly nonlinear shallow water waves with weak transverse variation is the celebrated KP equation [2]:
(1)∂∂x1ghηt+ηx+32hηηx+h2γ2ηxxx+12ηyy=0,where g is gravity, h is constant mean height, η=ηx,y,t is wave high above h, γ=1−τ/3, τ=T/ρgh2 is dimensionless surface tension coefficient, and ρ is density. Obviously, if η is not related to y, then Equation (1) reduces to the Korteweg-de Vries (KdV) equation [3]. So, KP Equation (1) can be regarded as an extension of the KdV equation to two-dimensional space. In nondimensional form, KP Equation (1) can be rescaled as [1]:
(2)ut+6uux+uxxxx+3σuyy=0,where σ=1 and σ=−1, depending on the sign of γ, correspond to KPII and KPI equations, respectively. The KPI equation exists lump solutions rationally localized in all directions in the space, for example, those [4] obtained by Ma. At the same time, n-soliton solutions with a single phase, called line solitons, have been found for the KPII equation by the Hirota bilinear method [5].Recently, Abolwitz and Baldwin [1] obtained line two- and three-soliton solutions with X-, Y-, and more complex-type interactions of the KPII equation. Interestingly, such types of X, Y, and more complex (for example 3-in-2-out) interactions have been observed to appear frequently in shallow water on two relatively flat beaches. Three types (short stem, long stem with lower and higher incoming line solitons) of X- and typical Y-type ocean wave interactions observed by Abolwitz and Baldwin in Mexico and California occur daily, shortly before and after low tide. However, these complex interactions like the 3-in-2-out one with three incoming line solitons on one side and two outgoing line solitons on the other side are much less frequent than X- and Y-type interactions. If the distances in the open-ocean direction are large enough, the destructive tsunami waves can merge in similar ways to the initial formation of an X- or Y-type wave, such as the observations that happen in the tsunami induced by the 2011 Japanese Tohoku-Oki earthquake indicating there was a merging phenomenon from a cylindrical-wave-type interaction. For more details, refer to Ref. [1].The aim of this article is to derive not only the X-, Y-, and 3-in-2-out-type interactions but also other novel solutions forming ring, conical, columnar, and lump structures from a fractional version of KP Equation (2):
(3)ut,cα+6uux,aα+ux,a3αx,aα+3σuy,b2α=0,0<α=pq≤1,σ=±1,with the conformable fractional partial derivatives Dx,aα, Dy,bα, and Dt,cα (see the next section for Definitions 1 and 2). Some novel spatial structures obtained in fractional KP Equation (3) cannot be derived from corresponding integral-order KP Equation (2). This is due to the irreplaceable role of fractional calculus [6–23].The rest of this paper is organized as follows. In Section2, we present the conformable fractional partial derivatives and integrals and give a physical interpretation of the conformable fractional derivative. In Section 3, we derive fractional KP Equation (3) from the corresponding fractional Lax pair. In Section 4, we transform fractional KP Equation (3) into bilinear form and then not only construct its fractional line solitons with X-, Y-, and 3-in-2-out interactions but also simulate a falling and spreading process of a columnar structure in the fractional KPII equation. At the same time, we construct a fractional rational solution of fractional KP Equation (3) and employ it to form ring, conical, and lump structures of the fractional KPI equation. In Section 5, we conclude this paper.


## 2. Definitions and Physical Interpretation

We give, in this section, the definitions of the conformable fractional partial derivatives and integrals. Since these definitions are based on the similar definitions in Ref. [24], we add this reference when defining them.Definition 1 [24].
For the given functionux,t:R×R⟶C, arbitrary constant a∈R, and fractional-order α∈0,1, the conformable fractional partial derivative with respect to x has the definition:
(4)Dx,aαux,t=ux,aαx,t=limε→0ux+εx−a1−α,t−ux,tε,0<α=pq≤1,where p and q are coprime and positive integers guaranteeing x−aα can return a unique real or complex value. Here, if q is an even number, then x−aα is stipulated to locate in the first quadrant of the complex plane and has the smallest imaginary part for arbitrary x<a. Note that ux,aαx,tx=a=limx→a+ux,aαx,t is a supplement of Definition 1. As a special case of Definition 1, we have the following Definition 2.Definition 2 [24].
If functionux,t:R×R⟶R, arbitrary constant a∈R, and fractional-order α∈0,1 are given, then we define the conformable fractional partial derivative with respect to x:
(5)Dx,aαux,t=ux,aαx,t=limε→0ux+εx−ap1/p−1/q,t−ux,tε,0<α=pq≤1,where p and q are coprime and positive integers, and p is an even number. A supplement of Definition 2 is the stipulation ux,aαx,tx=a=limx→a+ux,aαx,t.Definition 3 [24].
The conformable fractional integral of functionux,t:R×R⟶C is defined:
(6)Ix,aαux,t=∂x,a−αux,t=∫axx−a1−αux,tdx,0<α=pq≤1,where a∈R is an arbitrary constant, and x−aα is stipulated to return a unique real or complex value.Definition 4 [24].
Suppose thatp is an even number and a∈R is an arbitrary constant, the conformable fractional integral of function ux,t:R×R⟶R is defined:
(7)Ix,aαux,t=∂x,a−αux,t=∫axx−ap1/p−1/qux,tdx,0<α=pq≤1.We may sayux,t is α-differentiable with respect to x in the case that Equation (4) or (5) holds. By comparisons, we can see that Definitions 1–4 include those in [24, 25] as special cases. The defined conformable partial derivatives and integrals in Definitions 1–4 have some basic properties, which are all similar to the ones in [24, 25]. Only one of them is listed here, and the rest are omitted, that is, Dx,aαux,t=x−a1−αuxx,t. In other words, Dx,aαux,t is the product of the rate of change of ux,t along the x-axis and x−a1−α, which is helpful to understand the physical meaning of conformable fractional partial derivatives.It is a challenging problem to explore the physical or geometric meaning of fractional derivatives and integrals. As Xue [26] pointed out, although fractional calculus was born more than 300 years ago, there are no widely accepted physical and geometric interpretations of fractional derivatives and fractional integrals. Inspired by Podlubny’s work [27], we would like to give a physical explanation of the conformable fractional derivative Tαft [25]. Suppose an object moves with variable acceleration, and its variable acceleration and velocity are at and vt at time t, respectively. Furthermore, it is assumed that the relationship between the cosmic time [27] (real time) T=gtτ and the individual time [27] (inaccurate timing) τ satisfies
(8)gtτ=t2−α+τ2−α−t2−α.With the above preparation, we can see that the real velocityvOt of the moving object with the variable acceleration aτ in the individual time measure at cosmic time T is
(9)vOt=∫0taτdgtτ=∫0tv′ττ1−αdτ,from which the variable accelerationaOt at time t can be derived as the following:
(10)aOt=v′Ot=ddt∫0tv′ττ1−αdτ=v′tt1−α=Tαvt.Therefore, theα-order conformable fractional derivative Tαvt of the velocity function vt is just the variable acceleration aOt from the observer’s perspective with accurate time. This is a physical interpretation of the conformable fractional derivative. For the special case α=1, Equations (8) and (10) reduce gtτ=τ and aOt=v′t=at, respectively. That is to say, when there is no difference between the cosmic time and the individual time, the corresponding variable accelerations are equal.


## 3. Lax Representation

For the Lax representation of fractional KP Equation (3), we have the following Theorem 5.Theorem 5.
Fractional KP Equation (3) is the Lax integrable, the Lax representation of which is
(11)L,A−∂t,cα=0,with the Lax pair
(12)L=σ∂y,bα+∂x,a2α+u,A=−4∂x,a3α−6u∂x,aα−3ux,aα+3σ∂x,a−αuy,bα.Proof.
It is easy to derive that(13)LA=−4∂x,a5α−6u∂x,a3α−3ux,aα∂x,a2α−u4∂x,a3α+6u∂x,aα+3ux,aα−12ux,aα∂x,a2α−12ux,a2α∂x,aα−3ux,a3α+σ−4ux,a3α−6u∂x,aα−3ux,aα+3σ∂x,a−αuy,bα∂y,bα+3σ∂x,a−αuy,bα∂x,a2α+u+3σ∂x,a−αuy,b2α,(14)AL=−4∂x,a5α−6u∂x,a3α−3ux,aα∂x,a2α−u4∂x,a3α+6u∂x,aα+3ux,aα−12ux,aα∂x,a2α−12ux,a2α∂x,aα−4ux,a3α−6uux,aα+σ−4∂x,a3α−6u∂x,aα−3ux,aα+3σ∂x,a−αuy,bα∂y,bα+3σ∂x,a−αuy,aα∂x,a2α+u.
On the other hand, we can verify the relation:(15)L,−Dt,cα=ut,cα.
Substituting Equations (13) and (14) into Equation (11) yields
(16)ut,cα+6uux,aα+ux,a3α+3σ∂x,a−αuy,b2α=0.
Further takingα-order conformable partial derivative with respect to x, we finally transform Equation (16) into Equation (3). Thus, we complete the proof.


## 4. Fractional Line Soliton Solutions and Their Interactions

There are four parts in this section for the fractional line soliton solutions and their interactions of fractional KP Equation (3): one part is the bilinearization of Equation (3), and the other three parts are the simulations for X-, Y-, and 3-in-2-out-type interactions and peakon, ring, conical, columnar, and lump structures in the obtained fractional line one-, two-, and three-soliton solutions.


### 4.1. Bilinearization and Formula of Soliton Solutions

We have the following Theorem6 for the bilinear form of fractional KP Equation (3).Theorem 6.
Let(17)u=2∂x,a2αlnfx,y,t,then fractional KP Equation (3) has the following bilinear form:
(18)Ht,cαHx,aα+Hx,a4α+3σHy,b2αf·f=0,where Hx,aα, Hy,bα, and Ht,cα are the fractional versions [28] of the Hirota bilinear operator [5]:
(19)Hx,amαHt,cnαf·g=∂x,aα−∂x′,aαm∂t,cα−∂t′,cαnfx,tgx′,t′x′=x,t′=t,m,n∈N.Proof.
Substituting Equation (17) into fractional KP Equation (3), we convert Equation (3) into
(20)∂t,cα∂x,aαf−∂t,cαf∂x,aαf+f∂x,a4αf−4∂x,aαf∂x,a3αf+3∂x,a2αf2+3σf∂y,b2αf−∂y,bαf2=0,which can be rewritten as Equation (18) by using the fractional Hirota’s bilinear operator (19). The proof ends.With the help of the fractional bilinear form (18), we then obtain fractional n-soliton solutions of fractional KP Equation (3) as follows:
(21)u=2∂x,a2αln∑μ=0,1e∑j=1nμjξj+∑1≤j<lnμjμlAjl,where Σμ=0,1 denotes all the possible combinations of each μj=0,1 for j=1,2,⋯,n,(22)ξj=kjx−aαα+pjy−bαα−kj2+3σpj2t−cαα+ξj0,while eAjl=kj−kl2−σpj−pl2/kj+kl2−σpj−pl2, kj, pj, and ξj0 are constants.In particular, whenn=1, Equation (21) gives the fractional one-soliton solution of fractional KP Equation (3):
(23)u=2∂x,a2αln1+eξ1=k122sech2ξ12,with
(24)ξ1=k1x−aαα+p1y−bαα−k12+3σp12t−cαα+ξ10.Compared with one-soliton of integer orderα=1, we can see from Figures 1–6 that the fractional ones of the corresponding fractional KPII equation have different profiles. Figure 1 shows there are different wave widths. But the essential difference is that every fractional soliton has a certain inclination in the opposite direction of its propagation, as if it is affected by the resistance of the carrier in the forward direction. Such a characteristic also appears in Figure 2. More interestingly, the fractional one-soliton in Figure 2 has a bend in the middle; such phenomena often occur on the edges of flat beaches. In addition to the same amplitude k12/2 and the similar bell profile known from Equation (23), to gain more insights in the dynamics of the one-solitons in Figure 1 mathematically, we insert a=0, b=2, c=3, k1=1, p1=1/2, σ=1, and ξ10=0 into Equation (24); then, the velocity along the x-axis of the one-solitons can be derived by the expression
(25)ẋ=74t−3α−142α2+74t−3α1/α−1,which shows that all the fractional one-solitons propagate along the x-axis with variable velocities depending on time t and fractional order α. However, the integral-order soliton propagates at a constant velocity of 7/4. At time t=0, the velocities corresponding to α=53/99,73/99,89/99 are about 0.62, 2.14, and 2.02, respectively. The simulation results show that when t=4, all these fractional and integral-order solitons propagate to almost the same position along the positive direction of the x-axis, but when t=8, their positions in the positive direction of x-axis are just opposite to those when t=0. If ẋ<0, such as t=−18, all these solitons, whether integer or fractional, propagate along the negative x-axis, and when t=−18, they all reach the negative x-axis.Figure 1
Profiles of fractional one-solitonu with different orders by selecting a=0, b=2, c=3, k1=1, p1=1/2, σ=1, ξ10=0, y=−40, and t=0.Figure 2
Fractional one-solitonu with α=73/99 by selecting a=1, b=2, c=3, k1=1, p1=1/2, σ=1, ξ10=0, and t=0.Figure 3
Profiles of fractional one-soliton∣u∣ with different orders by selecting a=1, b=2, c=3, k1=1, p1=1/2, σ=1, ξ10=0, y=−40, and t=0.Figure 4
Fractional one-soliton∣u∣ with α=53/98 by selecting a=1, b=1.5, c=0, k1=1, p1=1/2, σ=1, ξ10=0, and t=0.Figure 5
Peakon formed by the fractional one-solitonu with α=6/11 by selecting a=1, b=1.5, c=0, k1=1, p1=1/2, σ=1, ξ10=0, and t=0.Figure 6
Profiles of fractional one-solitonu with different orders by selecting a=1, b=1.5, c=0, k1=1, p1=1/2, σ=1, ξ10=0, y=−1, and t=0.Figure3 shows there are not only different wave widths but also amplitudes between the integer-order one-soliton and the fractional ones. Besides, a fractional one-soliton with α=53/98 in Figure 4 appears two obvious wavelet peaks in the negative direction of x-axis. These different wave widths, amplitudes, and wavelet peaks in Figures 3 and 4 benefit from the fact that both x−aα and y−bα in Equation (24) return complex values according to Definition 1 for arbitrary x<a and y<b when α=p/q and q is an even number. More specifically, it is known that sech2z≤4 always holds for any z∈R, but it does not hold when z=x+iyx,y≠0∈R is a complex number. It is because that sechz2=∣1/cosh2x−sin2y∣ can reach any given positive number as long as x and y are selected appropriately. For example, sech0.2+1.3i2>8.9.Figures5 and 6 indicate that the uplift in the middle of the fractional line one-solitons appears some features of peakons, the peaks of which have discontinuous first derivatives of integer order. For the fractional one-soliton u with α=6/11, the limit of uxx,y,t at the point 1,1.5,0 is infinity. The reason that the peakons appear in Figures 5 and 6 is due to Definition 2 which let Equation (24) perform the computation x−aα=x−ap1/q and y−bα=y−bp1/q when α=p/q and p is an even number.


### 4.2. Fractional X-, Y-, and 3-in-2-out-Type Interactions

Assignn=2 to Equation (21), we first write the fractional two-soliton solution of fractional KP Equation (3) as
(26)u=2∂x,a2αln1+eξ1+eξ2+eξ1+ξ2+A12,whereξ1 satisfies Equation (24), and
(27)ξ2=k2x−aαα+p2y−bαα−k22+3σp22t−cαα+ξ20,eA12=k1−k22−σp1−p22k1+k22−σp1−p22.Next, we use the fractional two-soliton solution (26) to simulate the X- and Y-type interactions [1] observed by Abolwitz and Baldwin. Since the phase shifts of the colliding 2-solitons in a two-dimensional KdV equation are determined by eA12 in Equation (27), more details can be found in [29], we choose appropriate values of eA12 to simulate the X- and Y-type interactions through the fractional KPII equation. Based on the selection of k1=1/2, k2=1/2, p1=2/3, p2=−2/3 [1], and σ=1 letting eA12≈2.3, Figure 7 shows an X-type interaction formed by the fractional two-soliton solution (26) with α=83/99. At the same time, it should be noted that a phase shift caused by the collision of two-solitons will form the so-called “stem,” so the change of eA12 will affect the length, height, and even disappearance of stem. By choosing k1=1/2, k2=1, p1=3/4, p2=1, and σ=1 so that eA12=0, Figure 8 simulates a Y-type interaction without a stem. In addition, two X-type interactions with long stem (eA12≈5.0×109) and short stem (eA12≈50.8) are shown in Figures 9 and 10, respectively. In Figure 11, we show another X-type interaction, the stem corresponding to eA12≈5.0×10−8 of which is lower than one incoming line soliton not as the ones having high stems in Figures 7, 9, and 10, while is almost equal to the other incoming line soliton in height.Figure 7
X-type interaction in the fractional line two-soliton solutionu with α=83/99, a=0, b=0, c=0, k1=1/2, k2=1/2, p1=2/3, p2=−2/3, σ=1, ξ10=0, and ξ20=0.Figure 8
Y-type interaction in the fractional line two-soliton solutionu with α=83/99, a=0, b=0, c=0, k1=1/2, k2=1, p1=3/4, p2=1, σ=1, ξ10=0, and ξ20=0.Figure 9
X-type interaction with a short stem in the fractional line two-soliton solutionu with α=83/99, a=0, b=0, c=0, k1=1/2, k2=1/2, p1=−1/4−10−2, p2=3/4, σ=1, ξ10=0, and ξ20=0.Figure 10
X-type interaction with a long stem in the fractional line two-soliton solutionu with α=83/99, a=0, b=0, c=0, k1=1/2, k2=1/2, p1=−1/2−10−10, p2=1/2, σ=1, ξ10=0, and ξ20=0.Figure 11
X-type interaction with a low stem in the fractional line two-soliton solutionu with α=83/99, a=0, b=0, c=0, k1=1, k2=1/2, p1=1/2−10−7, p2=0, σ=1, ξ10=0, and ξ20=0.Finally, we employ the fractional three-soliton solution of fractional KP Equation (3) to simulate not only X- and Y-type interactions but also 3-in-2-out-type interaction [1]. For this purpose, we write the fractional three-soliton solution by the substitution of n=3 into Equation (21):
(28)u=2∂x,a2αln1+eξ1+eξ2+eξ3+eξ1+ξ2+A12+eξ1+ξ3+A13+eξ2+ξ3+A23+eξ1+ξ2+ξ3+A12+A13+A23,where ξ1 satisfies Equation (24), eA12 and ξ2 are determined by Equation (27), and
(29)ξ3=k3x−aαα+p3y−bαα−k22+3σp32t−cαα+ξ30,(30)eA13=k1−k32−σp1−p32k1+k32−σp1−p32,eA23=k2−k32−σ2p2−p32k2+k32−σp2−p32.Substitutingk1=0.2, k2=1.2, k3=2.2, p1=−1/4, p2=−3/4, p3=−7/4, and σ=1 into Equations (24) and (27), we have eA12≈0.4, eA13≈0.5, and eA23≈4.2×10−17; then, an X-type interaction forming in the fractional line three-soliton solution (28) with α=83/99 is shown in Figure 12. Interestingly, the remarkable differences from Figure 11 are the stem in Figure 12 which is lower than all the two incoming line solitons and one curved outgoing line soliton. The appropriate parameters k1, k2, and k3, which are only slightly different from those in Figure 12, are selected; Figure 13 shows a Y-type structure (eA12≈0.56, eA13≈0.61, eA23=0) in the same fractional line three-soliton solution (28). In Figure 14, we show a 3-in-2 out-type interaction, one line soliton of which has high amplitude than all the incoming ones, where eA12≈0.08, eA13≈0.12, and eA23=0 are determined by selecting k1=1, k2=2, k3=2, p1=−1/3, p2=−2/3, p3=−5/3, and σ=1.Figure 12
X-type interaction in the fractional line three-soliton solutionu with α=83/99, a=0, b=0, c=0, k1=0.2, k2=1.2, k3=2.2, p1=−1/4, p2=−3/4, p3=−7/4, σ=1, ξ10=0, ξ20=0, and ξ30=0.Figure 13
Y-type interaction in the fractional line three-soliton solutionu with α=83/99, a=0, b=0, c=0, k1=0.13, k2=1.13, k3=2.13, p1=−1/4, p2=−3/4, p3=−7/4, σ=1, ξ10=0, ξ20=0, and ξ30=0.Figure 14
3-in-2-out-type interaction in the fractional line three-soliton solutionu with α=83/99, a=0, b=0, c=0, k1=1, k2=2, k3=2, p1=−1/3, p2=−2/3, p3=−5/3, σ=1, ξ10=0, ξ20=0, and ξ30=0.


### 4.3. Falling and Spreading Phenomena of a Columnar Structure

As mentioned earlier, all the values ofeA12, eA13, and eA23 in Figures 7–14 are nonnegative. From the forms of their compositions in Equations (27) and (30), we can see that eA12, eA13, and eA23 can also take negative values. In Figures 15(a)–15(l), we show a falling and spreading process of a columnar structure (eA12≈−0.3) in the fractional two-soliton solution (26) of the fractional KPII equation, where the parameters are selected as α=82/99, a=0, b=0, c=0, k1=−1/2, k2=−1/2, p1=1/4−10−4, p2=4/3, σ=1, ξ10=0, and ξ20=0, respectively. In this case, for any fixed time t, the limit of the amplitude u of this fractional two-soliton always tends to zero as x⟶±∞ and y⟶±∞. Figure 15(a) shows the columnar structure at the beginning t=0. With the passage of time, the top of the columnar began to fall vertically, forming a column in the opposite direction and widening continuously, until at t=18.83 the column suddenly shrinks into a very thin inverted cone, see Figures 15(b)–15(g). Then, the inverted cone becomes slender and splits into several inverted cones from about t=19.2. When t=20.5, the number of inverted cones began to decrease, and diamond pattern was formed on the background, see Figures 15(h) and 15(i). Finally, the number of inverted cones increases gradually, and the diamond background pattern gradually expands, see Figures 15(j)–15(l). In time inversion, the process of falling and spreading simulated in Figures 15 forms a long flat column, and the merging effect, as pointed out by Ablowitz and Baldwin [1], produced by multiple such columns may be useful to explain the generation of tsunami caused by earthquake.Figure 15
Falling and spreading phenomena of a columnar structure in the fractional two-soliton solutionu with α=82/99.
(a)
t=0(b)
t=5(c)
t=10(d)
t=15(e)
t=17(f)
t=18.82(g)
t=18.83(h)
t=20(i)
t=30(j)
t=50(k)
t=100(l)
t=500


### 4.4. Fractional Rational Solution and Lump Structure

Expanding Equation (18) yields
(31)fx,a;t,c2αf−fx,aαft,cα+fx,a4αf−4fx,a3αfx,aα+3fx,a2α2+3σfy,b2αf−fy,b2α2=0.Following Ma’s work [4], we suppose that
(32)f=g2+h2+a9,with
(33)g=a1x−aα+a2y−bα+a3t−cα+a4,(34)h=a5x−aα+a6y−bα+a7t−cα+a8,where aii=1,2,⋯,9 are undetermined constants. Substituting Equations (32)–(34) into Equation (31), collecting the coefficients of the same powers x−aαjy−bαlt−cαsj,l,s=0,1,2,⋯, we derive a set of algebraic equations about these undermined constants. Solving this set of algebraic equations, we can determine Equations (32)–(34) by
(35)a3=−3σ2a2a5a6+a1a22−a62a12+a52,(36)a7=−3σ−a22a5+2a1a2a6+a5a62a12+a52,a9=−a12+a523α2σa2a5−a1a62.In view of Equations (17) and (32)–(36), we obtain a rational solution of fractional KP Equation (3):
(37)u=4α2−x−aα+a2y−bα+3a22−a62t−cα+a422+a6y−bα+6a2a6t−cα+a822−α2/a62σx−aα+a2y−bα+3a22−a62t−cα2+a6y−bα+6a2a6t−cα2−α2/a62σ2by setting a1=1 and a5=0. Further inserting a4=0, a8=0, and σ=−1 into Equation (37), we reduce Equation (37) to a fractional lump solution of the fractional KPI equation:
(38)u=4α2−x−aα+a2y−bα+3a22−a62t−cα2+a6y−bα+6a2a6t−cα2+α2/a62x−aα+a2y−bα+3a22−a62t−cα2+a6y−bα+6a2a6t−cα2+α2/a622.In Figure16, the fractional lump solution (38) with α=91/99 of the fractional KPI equation is shown by setting a=0.5, b=0.1, c=−0.3, a2=−0.2, a6=1.2, and t=0. Here, we note that no matter whether x−a and y−b in Equation (38) are nonnegative or not, when α=p/q and q is a odd number, Definition 1 always returns a unique real number from x−aα and y−bα, respectively. Therefore, when we choose other different values of a and b, the solution (38) can still maintain such lump structure in the space coordinate plane; only the positions are different. In the next subsection, we will simulate different spatial evolution structures presented by solution (38) in different positions of the coordinate plane.Figure 16
Fractional lump solution∣u∣ with α=91/99 and a=0.5, b=0.1, c=−0.3, a2=−0.2, a6=1.2, and t=0.


### 4.5. Dynamic Evolutions from Ring and Conical Structures to Lump Structure

In order to simulate the different spatial evolution structures of solution (38) in the same evolution process, we pay attention to the positive and negative effects of x−a, y−b, and t−c, so that x−aα, y−bα, and t−cα are real numbers in a certain range and complex numbers in other ranges, see Figures 17–20 for examples. Figure 17 shows the dynamic evolution of the fractional lump solution (38) by the selection of α=51/98, a=0.2, b=0, c=0, a2=−0.6, and a6=0.8, where the spatial structure of ∣u∣ changes from ring to lump. Two profiles of the fractional lump solution (38) at time t=0 are shown in Figure 18. With a different fractional order α=61/98, Figure 19 shows another dynamic evolution of the fractional lump solution (38), where the spatial structure of ∣u∣ changes from cone to lump. Two different profiles of the fractional lump solution (38) at time t=0 are shown in Figure 20.Figure 17
Dynamic evolution of the fractional lump solution∣u∣ with α=51/98, a=0.2, b=0, c=0, a2=−0.6, and a6=0.8.
(a)
t=−20(b)
t=5(c)
t=20Figure 18
Profiles of the fractional lump solution∣u∣ with α=51/98, a=0.2, b=0, c=0, a2=−0.6, a6=0.8, and t=0.
(a)
x=−0.5(b)
y=−0.35Figure 19
Dynamic evolution of the fractional lump solution∣u∣ with α=61/98, a=0.2, b=0, c=0, a2=−0.6, and a6=0.8.
(a)
t=−20(b)
t=5(c)
t=20Figure 20
Profiles of the fractional lump solution∣u∣ with α=61/98, a=0.2, b=0, c=0, a2=−0.6, a6=0.8, and t=0.
(a)
x=−0.5(b)
y=−0.35


## 4.1. Bilinearization and Formula of Soliton Solutions

We have the following Theorem6 for the bilinear form of fractional KP Equation (3).Theorem 6.
Let(17)u=2∂x,a2αlnfx,y,t,then fractional KP Equation (3) has the following bilinear form:
(18)Ht,cαHx,aα+Hx,a4α+3σHy,b2αf·f=0,where Hx,aα, Hy,bα, and Ht,cα are the fractional versions [28] of the Hirota bilinear operator [5]:
(19)Hx,amαHt,cnαf·g=∂x,aα−∂x′,aαm∂t,cα−∂t′,cαnfx,tgx′,t′x′=x,t′=t,m,n∈N.Proof.
Substituting Equation (17) into fractional KP Equation (3), we convert Equation (3) into
(20)∂t,cα∂x,aαf−∂t,cαf∂x,aαf+f∂x,a4αf−4∂x,aαf∂x,a3αf+3∂x,a2αf2+3σf∂y,b2αf−∂y,bαf2=0,which can be rewritten as Equation (18) by using the fractional Hirota’s bilinear operator (19). The proof ends.With the help of the fractional bilinear form (18), we then obtain fractional n-soliton solutions of fractional KP Equation (3) as follows:
(21)u=2∂x,a2αln∑μ=0,1e∑j=1nμjξj+∑1≤j<lnμjμlAjl,where Σμ=0,1 denotes all the possible combinations of each μj=0,1 for j=1,2,⋯,n,(22)ξj=kjx−aαα+pjy−bαα−kj2+3σpj2t−cαα+ξj0,while eAjl=kj−kl2−σpj−pl2/kj+kl2−σpj−pl2, kj, pj, and ξj0 are constants.In particular, whenn=1, Equation (21) gives the fractional one-soliton solution of fractional KP Equation (3):
(23)u=2∂x,a2αln1+eξ1=k122sech2ξ12,with
(24)ξ1=k1x−aαα+p1y−bαα−k12+3σp12t−cαα+ξ10.Compared with one-soliton of integer orderα=1, we can see from Figures 1–6 that the fractional ones of the corresponding fractional KPII equation have different profiles. Figure 1 shows there are different wave widths. But the essential difference is that every fractional soliton has a certain inclination in the opposite direction of its propagation, as if it is affected by the resistance of the carrier in the forward direction. Such a characteristic also appears in Figure 2. More interestingly, the fractional one-soliton in Figure 2 has a bend in the middle; such phenomena often occur on the edges of flat beaches. In addition to the same amplitude k12/2 and the similar bell profile known from Equation (23), to gain more insights in the dynamics of the one-solitons in Figure 1 mathematically, we insert a=0, b=2, c=3, k1=1, p1=1/2, σ=1, and ξ10=0 into Equation (24); then, the velocity along the x-axis of the one-solitons can be derived by the expression
(25)ẋ=74t−3α−142α2+74t−3α1/α−1,which shows that all the fractional one-solitons propagate along the x-axis with variable velocities depending on time t and fractional order α. However, the integral-order soliton propagates at a constant velocity of 7/4. At time t=0, the velocities corresponding to α=53/99,73/99,89/99 are about 0.62, 2.14, and 2.02, respectively. The simulation results show that when t=4, all these fractional and integral-order solitons propagate to almost the same position along the positive direction of the x-axis, but when t=8, their positions in the positive direction of x-axis are just opposite to those when t=0. If ẋ<0, such as t=−18, all these solitons, whether integer or fractional, propagate along the negative x-axis, and when t=−18, they all reach the negative x-axis.Figure 1
Profiles of fractional one-solitonu with different orders by selecting a=0, b=2, c=3, k1=1, p1=1/2, σ=1, ξ10=0, y=−40, and t=0.Figure 2
Fractional one-solitonu with α=73/99 by selecting a=1, b=2, c=3, k1=1, p1=1/2, σ=1, ξ10=0, and t=0.Figure 3
Profiles of fractional one-soliton∣u∣ with different orders by selecting a=1, b=2, c=3, k1=1, p1=1/2, σ=1, ξ10=0, y=−40, and t=0.Figure 4
Fractional one-soliton∣u∣ with α=53/98 by selecting a=1, b=1.5, c=0, k1=1, p1=1/2, σ=1, ξ10=0, and t=0.Figure 5
Peakon formed by the fractional one-solitonu with α=6/11 by selecting a=1, b=1.5, c=0, k1=1, p1=1/2, σ=1, ξ10=0, and t=0.Figure 6
Profiles of fractional one-solitonu with different orders by selecting a=1, b=1.5, c=0, k1=1, p1=1/2, σ=1, ξ10=0, y=−1, and t=0.Figure3 shows there are not only different wave widths but also amplitudes between the integer-order one-soliton and the fractional ones. Besides, a fractional one-soliton with α=53/98 in Figure 4 appears two obvious wavelet peaks in the negative direction of x-axis. These different wave widths, amplitudes, and wavelet peaks in Figures 3 and 4 benefit from the fact that both x−aα and y−bα in Equation (24) return complex values according to Definition 1 for arbitrary x<a and y<b when α=p/q and q is an even number. More specifically, it is known that sech2z≤4 always holds for any z∈R, but it does not hold when z=x+iyx,y≠0∈R is a complex number. It is because that sechz2=∣1/cosh2x−sin2y∣ can reach any given positive number as long as x and y are selected appropriately. For example, sech0.2+1.3i2>8.9.Figures5 and 6 indicate that the uplift in the middle of the fractional line one-solitons appears some features of peakons, the peaks of which have discontinuous first derivatives of integer order. For the fractional one-soliton u with α=6/11, the limit of uxx,y,t at the point 1,1.5,0 is infinity. The reason that the peakons appear in Figures 5 and 6 is due to Definition 2 which let Equation (24) perform the computation x−aα=x−ap1/q and y−bα=y−bp1/q when α=p/q and p is an even number.


## 4.2. Fractional X-, Y-, and 3-in-2-out-Type Interactions

Assignn=2 to Equation (21), we first write the fractional two-soliton solution of fractional KP Equation (3) as
(26)u=2∂x,a2αln1+eξ1+eξ2+eξ1+ξ2+A12,whereξ1 satisfies Equation (24), and
(27)ξ2=k2x−aαα+p2y−bαα−k22+3σp22t−cαα+ξ20,eA12=k1−k22−σp1−p22k1+k22−σp1−p22.Next, we use the fractional two-soliton solution (26) to simulate the X- and Y-type interactions [1] observed by Abolwitz and Baldwin. Since the phase shifts of the colliding 2-solitons in a two-dimensional KdV equation are determined by eA12 in Equation (27), more details can be found in [29], we choose appropriate values of eA12 to simulate the X- and Y-type interactions through the fractional KPII equation. Based on the selection of k1=1/2, k2=1/2, p1=2/3, p2=−2/3 [1], and σ=1 letting eA12≈2.3, Figure 7 shows an X-type interaction formed by the fractional two-soliton solution (26) with α=83/99. At the same time, it should be noted that a phase shift caused by the collision of two-solitons will form the so-called “stem,” so the change of eA12 will affect the length, height, and even disappearance of stem. By choosing k1=1/2, k2=1, p1=3/4, p2=1, and σ=1 so that eA12=0, Figure 8 simulates a Y-type interaction without a stem. In addition, two X-type interactions with long stem (eA12≈5.0×109) and short stem (eA12≈50.8) are shown in Figures 9 and 10, respectively. In Figure 11, we show another X-type interaction, the stem corresponding to eA12≈5.0×10−8 of which is lower than one incoming line soliton not as the ones having high stems in Figures 7, 9, and 10, while is almost equal to the other incoming line soliton in height.Figure 7
X-type interaction in the fractional line two-soliton solutionu with α=83/99, a=0, b=0, c=0, k1=1/2, k2=1/2, p1=2/3, p2=−2/3, σ=1, ξ10=0, and ξ20=0.Figure 8
Y-type interaction in the fractional line two-soliton solutionu with α=83/99, a=0, b=0, c=0, k1=1/2, k2=1, p1=3/4, p2=1, σ=1, ξ10=0, and ξ20=0.Figure 9
X-type interaction with a short stem in the fractional line two-soliton solutionu with α=83/99, a=0, b=0, c=0, k1=1/2, k2=1/2, p1=−1/4−10−2, p2=3/4, σ=1, ξ10=0, and ξ20=0.Figure 10
X-type interaction with a long stem in the fractional line two-soliton solutionu with α=83/99, a=0, b=0, c=0, k1=1/2, k2=1/2, p1=−1/2−10−10, p2=1/2, σ=1, ξ10=0, and ξ20=0.Figure 11
X-type interaction with a low stem in the fractional line two-soliton solutionu with α=83/99, a=0, b=0, c=0, k1=1, k2=1/2, p1=1/2−10−7, p2=0, σ=1, ξ10=0, and ξ20=0.Finally, we employ the fractional three-soliton solution of fractional KP Equation (3) to simulate not only X- and Y-type interactions but also 3-in-2-out-type interaction [1]. For this purpose, we write the fractional three-soliton solution by the substitution of n=3 into Equation (21):
(28)u=2∂x,a2αln1+eξ1+eξ2+eξ3+eξ1+ξ2+A12+eξ1+ξ3+A13+eξ2+ξ3+A23+eξ1+ξ2+ξ3+A12+A13+A23,where ξ1 satisfies Equation (24), eA12 and ξ2 are determined by Equation (27), and
(29)ξ3=k3x−aαα+p3y−bαα−k22+3σp32t−cαα+ξ30,(30)eA13=k1−k32−σp1−p32k1+k32−σp1−p32,eA23=k2−k32−σ2p2−p32k2+k32−σp2−p32.Substitutingk1=0.2, k2=1.2, k3=2.2, p1=−1/4, p2=−3/4, p3=−7/4, and σ=1 into Equations (24) and (27), we have eA12≈0.4, eA13≈0.5, and eA23≈4.2×10−17; then, an X-type interaction forming in the fractional line three-soliton solution (28) with α=83/99 is shown in Figure 12. Interestingly, the remarkable differences from Figure 11 are the stem in Figure 12 which is lower than all the two incoming line solitons and one curved outgoing line soliton. The appropriate parameters k1, k2, and k3, which are only slightly different from those in Figure 12, are selected; Figure 13 shows a Y-type structure (eA12≈0.56, eA13≈0.61, eA23=0) in the same fractional line three-soliton solution (28). In Figure 14, we show a 3-in-2 out-type interaction, one line soliton of which has high amplitude than all the incoming ones, where eA12≈0.08, eA13≈0.12, and eA23=0 are determined by selecting k1=1, k2=2, k3=2, p1=−1/3, p2=−2/3, p3=−5/3, and σ=1.Figure 12
X-type interaction in the fractional line three-soliton solutionu with α=83/99, a=0, b=0, c=0, k1=0.2, k2=1.2, k3=2.2, p1=−1/4, p2=−3/4, p3=−7/4, σ=1, ξ10=0, ξ20=0, and ξ30=0.Figure 13
Y-type interaction in the fractional line three-soliton solutionu with α=83/99, a=0, b=0, c=0, k1=0.13, k2=1.13, k3=2.13, p1=−1/4, p2=−3/4, p3=−7/4, σ=1, ξ10=0, ξ20=0, and ξ30=0.Figure 14
3-in-2-out-type interaction in the fractional line three-soliton solutionu with α=83/99, a=0, b=0, c=0, k1=1, k2=2, k3=2, p1=−1/3, p2=−2/3, p3=−5/3, σ=1, ξ10=0, ξ20=0, and ξ30=0.


## 4.3. Falling and Spreading Phenomena of a Columnar Structure

As mentioned earlier, all the values ofeA12, eA13, and eA23 in Figures 7–14 are nonnegative. From the forms of their compositions in Equations (27) and (30), we can see that eA12, eA13, and eA23 can also take negative values. In Figures 15(a)–15(l), we show a falling and spreading process of a columnar structure (eA12≈−0.3) in the fractional two-soliton solution (26) of the fractional KPII equation, where the parameters are selected as α=82/99, a=0, b=0, c=0, k1=−1/2, k2=−1/2, p1=1/4−10−4, p2=4/3, σ=1, ξ10=0, and ξ20=0, respectively. In this case, for any fixed time t, the limit of the amplitude u of this fractional two-soliton always tends to zero as x⟶±∞ and y⟶±∞. Figure 15(a) shows the columnar structure at the beginning t=0. With the passage of time, the top of the columnar began to fall vertically, forming a column in the opposite direction and widening continuously, until at t=18.83 the column suddenly shrinks into a very thin inverted cone, see Figures 15(b)–15(g). Then, the inverted cone becomes slender and splits into several inverted cones from about t=19.2. When t=20.5, the number of inverted cones began to decrease, and diamond pattern was formed on the background, see Figures 15(h) and 15(i). Finally, the number of inverted cones increases gradually, and the diamond background pattern gradually expands, see Figures 15(j)–15(l). In time inversion, the process of falling and spreading simulated in Figures 15 forms a long flat column, and the merging effect, as pointed out by Ablowitz and Baldwin [1], produced by multiple such columns may be useful to explain the generation of tsunami caused by earthquake.Figure 15
Falling and spreading phenomena of a columnar structure in the fractional two-soliton solutionu with α=82/99.
(a)
t=0(b)
t=5(c)
t=10(d)
t=15(e)
t=17(f)
t=18.82(g)
t=18.83(h)
t=20(i)
t=30(j)
t=50(k)
t=100(l)
t=500


## 4.4. Fractional Rational Solution and Lump Structure

Expanding Equation (18) yields
(31)fx,a;t,c2αf−fx,aαft,cα+fx,a4αf−4fx,a3αfx,aα+3fx,a2α2+3σfy,b2αf−fy,b2α2=0.Following Ma’s work [4], we suppose that
(32)f=g2+h2+a9,with
(33)g=a1x−aα+a2y−bα+a3t−cα+a4,(34)h=a5x−aα+a6y−bα+a7t−cα+a8,where aii=1,2,⋯,9 are undetermined constants. Substituting Equations (32)–(34) into Equation (31), collecting the coefficients of the same powers x−aαjy−bαlt−cαsj,l,s=0,1,2,⋯, we derive a set of algebraic equations about these undermined constants. Solving this set of algebraic equations, we can determine Equations (32)–(34) by
(35)a3=−3σ2a2a5a6+a1a22−a62a12+a52,(36)a7=−3σ−a22a5+2a1a2a6+a5a62a12+a52,a9=−a12+a523α2σa2a5−a1a62.In view of Equations (17) and (32)–(36), we obtain a rational solution of fractional KP Equation (3):
(37)u=4α2−x−aα+a2y−bα+3a22−a62t−cα+a422+a6y−bα+6a2a6t−cα+a822−α2/a62σx−aα+a2y−bα+3a22−a62t−cα2+a6y−bα+6a2a6t−cα2−α2/a62σ2by setting a1=1 and a5=0. Further inserting a4=0, a8=0, and σ=−1 into Equation (37), we reduce Equation (37) to a fractional lump solution of the fractional KPI equation:
(38)u=4α2−x−aα+a2y−bα+3a22−a62t−cα2+a6y−bα+6a2a6t−cα2+α2/a62x−aα+a2y−bα+3a22−a62t−cα2+a6y−bα+6a2a6t−cα2+α2/a622.In Figure16, the fractional lump solution (38) with α=91/99 of the fractional KPI equation is shown by setting a=0.5, b=0.1, c=−0.3, a2=−0.2, a6=1.2, and t=0. Here, we note that no matter whether x−a and y−b in Equation (38) are nonnegative or not, when α=p/q and q is a odd number, Definition 1 always returns a unique real number from x−aα and y−bα, respectively. Therefore, when we choose other different values of a and b, the solution (38) can still maintain such lump structure in the space coordinate plane; only the positions are different. In the next subsection, we will simulate different spatial evolution structures presented by solution (38) in different positions of the coordinate plane.Figure 16
Fractional lump solution∣u∣ with α=91/99 and a=0.5, b=0.1, c=−0.3, a2=−0.2, a6=1.2, and t=0.


## 4.5. Dynamic Evolutions from Ring and Conical Structures to Lump Structure

In order to simulate the different spatial evolution structures of solution (38) in the same evolution process, we pay attention to the positive and negative effects of x−a, y−b, and t−c, so that x−aα, y−bα, and t−cα are real numbers in a certain range and complex numbers in other ranges, see Figures 17–20 for examples. Figure 17 shows the dynamic evolution of the fractional lump solution (38) by the selection of α=51/98, a=0.2, b=0, c=0, a2=−0.6, and a6=0.8, where the spatial structure of ∣u∣ changes from ring to lump. Two profiles of the fractional lump solution (38) at time t=0 are shown in Figure 18. With a different fractional order α=61/98, Figure 19 shows another dynamic evolution of the fractional lump solution (38), where the spatial structure of ∣u∣ changes from cone to lump. Two different profiles of the fractional lump solution (38) at time t=0 are shown in Figure 20.Figure 17
Dynamic evolution of the fractional lump solution∣u∣ with α=51/98, a=0.2, b=0, c=0, a2=−0.6, and a6=0.8.
(a)
t=−20(b)
t=5(c)
t=20Figure 18
Profiles of the fractional lump solution∣u∣ with α=51/98, a=0.2, b=0, c=0, a2=−0.6, a6=0.8, and t=0.
(a)
x=−0.5(b)
y=−0.35Figure 19
Dynamic evolution of the fractional lump solution∣u∣ with α=61/98, a=0.2, b=0, c=0, a2=−0.6, and a6=0.8.
(a)
t=−20(b)
t=5(c)
t=20Figure 20
Profiles of the fractional lump solution∣u∣ with α=61/98, a=0.2, b=0, c=0, a2=−0.6, a6=0.8, and t=0.
(a)
x=−0.5(b)
y=−0.35


## 5. Conclusions

Based on the definitions of the conformable fractional partial derivatives and integrals, fractional KP Equation (3) with Lax integrability is derived. Then, a formula (21) of the fractional n-soliton solutions and rational solution (37) of KP Equation (3) is obtained by employing the Hirota bilinear method with fractional bilinear operators. When the fractional order α=1, most of the simulated fractional X-, Y-, and 3-in-2-out-type interactions and fractional lump structures deform to the ones which can be constructed in [1, 4]. Even though, the selectivity of fractional order α results in the variation of velocity of soliton evolution propagation, the nonlinearity of motion trajectory, the inclination and steepness of appearance profile, the changeability of amplitude, the bend and wavelet peaks as well as peakon in the fractional line one-soliton, the falling and spreading process of a columnar structure in the fractional line two-soliton, and the ring and conical structures before lump structure generation are novel. This is due to the fractional partial derivatives contained in the fractional KP equation, which make the fractional soliton and lump solutions show different nonlinearity at different positions in the coordinate plane. In order to better understand the conformable fractional derivatives and integrals, this paper also attempts to give a physical explanation of the conformable fractional derivative from the perspective of variable acceleration. Extending fractional calculus to other important topics of nonlinear integrable systems is worth studying, such as soliton molecules [30], full reversal symmetric multiple soliton solutions [31], financial rogue wave [32], long-time asymptotic behavior [33], initial-boundary value problems [34], and high-dimensional hierarchies of evolution equations and their Hamiltonian structures [35].



---

*Source: 6664039-2021-02-27.xml*