Lump, lumpoff and predictable rogue wave solutions to a dimensionally reduced Hirota bilinear equation
Wang Haifeng, Zhang Yufeng
School of Mathematics, China University of Mining and Technology, Xuzhou 221116, China

 

† Corresponding author. E-mail: zhangyfcumt@163.com

Project supported by the National Natural Science Foundation of China (Grant No. 11971475).

Abstract

We study a simplified (3+1)-dimensional model equation and construct a lump solution for the special case of z = y using the Hirota bilinear method. Then, a more general form of lump solution is constructed, which contains more arbitrary autocephalous parameters. In addition, a lumpoff solution is also derived based on the general lump solutions and a stripe soliton. Furthermore, we figure out instanton/rogue wave solutions via introducing two stripe solitons. Finally, one can better illustrate these propagation phenomena of these solutions by analyzing images.

1. Introduction

In recent years, the theme of lump waves in nonlinear science has attracted great interests because it was regarded as suitable prototypes to model rogue wave dynamics in both nonlinear optics[1] and oceanography.[2] In Ref. [3] the author proposed an effective method for constructing the lump solutions to the Kadomtsev–Petviashvili equation, and then all second-order lump solutions were systematically presented in Ref. [4] by Ma and Zhou. Henceforth, many nonlinear partial differential equations also admit lump solutions, such as the generalized Kadomtsev–Petviashvili–Boussinesq equation,[5] the Boussinesq–Kadomtsev–Petviashvili equation,[6] the (2+1)-dimensional Ito equation,[7] and the (2+1)-dimensional Burgers equation.[8] etc.[819] Many methods can help us study lump solutions, such as the multiple variable separation method[20] and the common expression method.[21] Among them the Hirota bilinear method is recognized as an effective approach to find lump solutions to nonlinear partial differential equations. The large amplitude wave related to the rogue wave is an extreme event that appears on the surface of the ocean. Such waves can be accompanied by deep grooves (holes) that appear before and/or after the maximum peak. Recently, special attention to theories and experiments has spread from oceanography into several other areas of research, such as nonlinear optical systems,[22,23] fluid dynamics and atmosphere.[2431] These developments help to understand the physical mechanisms of the flow phenomenon.

Jia and Lou[32] gave one special case of the lump solutions to the (2+1)-dimensional KP equation based on the general mathematical results established in Refs. [3,4] After that, the lumpoff solution and the instantan/rogue wave solution were derived. Based on the completed study, a series of related work was carried out.[3336] In this paper, we mainly consider a new (3+1)-dimensional model equation introduced by using a multivariate polynomial in Ref. [37] and then discuss the lump dynamics related to the new Hirota bilinear equation.[38] The form of the Hirota bilinear equation is

which is equivalent to

where f = f(x,y,z,t), and the D-operator is defined by[39,40]

Based on the theories of the binary Bell polynomials and the bilinear operator, one can introduce a variable transformation

then equation (2) is transformed into a (3+1)-dimensional model as follows:

The solution u of Eq. (3) will be analytical when the solution f of Eq. (2) is positive.

This paper is arranged as follows: In Section 2, we construct a general lump solution to the reduced (3+1)-dimensional model equation. On the basis of the general lump solution, a lumpoff solution is also derived with a stripe soliton in Section 3. In Section 4, we find that the lump will become a rough wave when a pair of solitons was determined by the lump and can only visible at a specific point in time.

2. Lump solution to reduction with z = y

We first take a special reduction z = y to construct the lump solutions to the (3+1)-dimensional model equation (3). Then equation (2) turns out to be

It follows that equation (3) is converted to

2.1. Lump solutions

To get the single lump solution of the equivalent equation (4), we take a general quadratic function as follows:

where x1 = x, x2 = y, x3 = t, x0 = 1, and aij, ij and f0 are constants to be determined. In Eq. (6), f has an expansion form as follows:

which contains eleven parameters to be determined. Substituting Eq. (7) into Eq. (4) and collecting the coefficients of x, y, t, we obtain ten equations as follows:

which can be used to determine the parameters of aij and f0. Solving Eqs. (8)–(17), one can get the following five constraints of ai3 (i = 0,1,2,3) and f0 as

where a01, a02, a11, a12, a22 are arbitrary constants, and a00 is canceled out by Eq. (22).

According to Eqs. (18)–(22), two nonzero conditions for a22 ≠ 0 and should be satisfied to find the solution to Eqs. (8)–(17).

Substituting the result Eqs. (18)–(22) into Eq. (7), one finds

Further, we get the lump solution u of Eq. (5) as follows:

where a01, a02, a11, a12, a22 are arbitrary constants.

Thus, the solution to Eqs. (23) and (24) containing ten parameters with five constraints and two nonzero conditions represents the general lump solutions to Eq. (3) with z = y. Next we will explore whether fewer constraints can get the similar lump solutions.

2.2. Lump solutions with more freedom

To contain more free parameters in the lump solution, we take the general quadratic function

where

with x being a column vector matrix, AR4 × 4 a symmetric matrix, and f0 a positive real parameter to be determined. Here f can be rewritten as

Different from Eq. (6), aij is redefined as

where

are n-dimensional vectors, and mλ,k (λ = 1,2,…,n and k = 1,2,3,4) are real scalar parameters to realize.

Substituting Eqs. (26)–(28) into Eq. (4) and eliminating the coefficients of x, y, t yield two constraints as follows:

According to the above results, it is not hard to find that Eqs. (26)–(28) include infinitely many parameters. Based on the results of Refs. [3,4], n = 3 is a more right choice for general lump solutions to Eq. (3). Therefore, one can rewrite constraint conditions by taking n = 3 as

where mλ,1, mλ,2 and mλ,3, λ = 1,2,3 are real constants. Thus, we obtain the general lump solution u of Eq. (5) as follows:

with

To precisely explore the moving path of lump waves expressed by Eq. (26), we give the critical point of the lump waves. Taking fx = fy = 0, the positions of the general single lump solution is

This means the lump waves move along the straight line

To further explore the propagation characteristics of lump wave (35) in detail, we present Fig. 1 shown with the arbitrary parameters chosen as

Fig. 1. Plots of the lump solution (35) for Eq. (5) with the parameters (39): (a) three-dimensional plot at time t = 0, (b) density plot, (b) the contour plot about the moving path described by the straight line (38).

Figure 1(a) describes the spatial structure of lump solution (35) with Eq. (39). Figure 1(b) presents the density plot of lump solution at t = 0. Figure 1(c) shows the relevant contour map of the lump solution when t = 0, t = 15 and t = 30, and the blue line expresses the propagation path (38), i.e. y = 2x.

3. Lumpoff solutions

The interaction between lump waves and stripe soliton waves formed the lumpoff solution which is also called a cutoff lump. It means that the lump solitons can be cut by stripe solitons before or after a specific time. Therefore, we consider the lumpoff solution to the reduced (3+1)-dimensional model (3) in this section. The lumpoff solution of the reduced (3+1)-dimensional model (3) has the form as follows:

where flump is given by Eq. (36), and the k, m01, m02, m03 and m04 are undetermined parameters.

Sulution (40) exhibits that the cutoff solutions composed of the lump solution and the exponential solutions. When m01 + m02x + m03y + m04 t > 0, the exponentiation part will become the dominant position. For this reason, the lump solution will only appear at m01 + m02x + m03y + m04t < 0.

Inserting Eq. (40) into Eq. (4), we have

where a11, a12 and a22 are defined in Eq. (27), with k and m01 being arbitrary constants.

The above results show that the parameters of the soliton parts m02, m03 and m04 are related to the parameters of the bulk solution (36), i.e., a11, a12 and a22. This means that the existence of the solitons is determined by the existence of a lump. Thus, the appearance of stripe soliton coincides with lump soliton. Once the soliton is induced, the lump will be invisible until the exponential part dominates. It implies that the lump is cut off by the solitons induced by itself. We further obtain the following lumpoff solution u of Eq. (3) with z = y by the transformation u = 2(ln f)x,

Then, a set of parameters are provided to better see the corresponding motion characteristics to lumpoff solutions to the dispersive relations in Eq. (41),

and the parameters of stripe soliton part can be given by

According to Eq. (40), we can find that the lump part is unchanged, thus the moving path of lump waves for Eq. (42) is invariable to expression in Eq. (38), one has y = 2x.

From Fig. 2, one can find the process of propagation in the space for lumpoff waves at t = −2, t = 0, t = 5, t = 20 with the selection of Eq. (44) {, , }. Additionally, Fig. 2 illustrates that the lumpoff wave is cut by the produced soliton. The lump wave appears when m01 + m02x + m03y + m04t < 0, soon afterwards the lump is cutoff by the induced soliton and eventually disappears. Figures 3(a)3(c) display the contour plots at t = −2, t = 0, t = 5 and t = 20. Figures 3(d) and 3(e) show the relevant density plots. Moreover, the similar phenomena can be analyzed when we choose {, , .

Fig. 2. Profiles of the lumpoff waves (42) for Eq. (3) with the parameters (43) and the selection of Eq. (44) {, , } at different times: (a) t = −2, (b) t = 0, (c) t = 5, (d) t = 20.
Fig. 3. Relevant contour plots of the lumpoff waves (42) for Eq. (3) with the parameters (43) at different times: (a) t = −2, (b) t = 0, (c) t = 5, (d) t = 20; and density plots of lumpoff waves (42) for Eq. (3) with the parameters (43) at times: (e) t = −2, (f) t = 0, (g) t = 5, (h) t = 20.
4. Instanton/rogue wave solution

An instanton/rogue wave is a localized wave decayed in all space and time directions. The motivation of studying predictable rogue wave is that the emerging time and place of this special rogue wave can be predicted. Thus, we focus on the instanton wave solutions to the reduced (3+1)-dimensional model (3) in this section. The instanton wave solutions of the reduced (3+1)-dimensional model (3) has the form as follows:

where flump is given by Eq. (36), the parameters m02, m03 and m04 are provided by Eq. (41), h and m01 are arbitrary constants, and l is another restricted condition that will be identified by later. Based on the previous analysis, one can find that the lumpoff soliton only appears at m01 + m02x + m03y + m04t < 0, and becomes invisible for m01 + m02x + m03y + m04 t > 0. However, the lump waves are both invisible for m01 + m02x + m03y + m04t < 0 and m01 + m02x + m03y + m04t > 0 due to the existence of a pair of exponential solutions. The particular rogue waves will appear only when the lump waves move to the line m01 + m02x + m03y + m04 t ∼ 0. Substituting Eq. (45) into Eq. (4), collecting the coefficients and considering the restricted conditions (41), one has

It follows that we obtain the particular rogue solution u of Eq. (3) with z = y as

Our results show that a pair of solitons are caused by block waves. The lump wave always disappears because the cosh part being predominant, which results in the lump visible only when it moves to the line m01 + m02x + m03y + m04t ∼ 0. Then, the lump continue to move until achieve the peak at the center line m01 + m02x + m03y + m04t = 0 of a pair of resonance stripe soliton, the lump soliton becomes rogue soliton. We therefore find that the time roughly appears to read

and the position

A set of parameters is selected to explain the phenomenon of predictable rogue wave (47) as

Substituting it into Eq. (44) yields the moving path

and the travelling direction

Thus, the time of the rogue waves arise at approximately t = −0.17 in (x = −4.31, y = −8.62) with arriving its peak based on Eqs. (48) and (49).

Figures 4(a)4(e) indicate the process of propagation in the special rogue wave at different times t = −5, t = −3, t = 0, t = 2, t = 10, respectively, and panel (f) is the wave propagation along the y-axis at t = −5, t = −3, t = 0, t = 2, t = 10. Figure 5 shows the corresponding density maps of the predictable rogue wave at t = −5, t = −3, t = 0, t = 2, t = 10.

Fig. 4. Profiles of the special instanton (47) for Eq. (3) with the parameters (50) at times: (a) t = −5, (b) t = −3, (c) t = 0, (d) t = 2, (e) t = 10; (f) the relevant wave height of predictable rogue waves in y = 0 plane at t = −3 (blue), t = 0 (red) and t = 2 (black).
Fig. 5. Density plots of the special instanton (47) for Eq. (3) with the parameters (50) at times: (a) t = −5, (b) t = −3, (c) t = 0, (d) t = 2, (e) t = 10.
5. Conclusion

In summary, we have studied the (3+1)-dimensional model (3) with z = y, and constructed the lump solution (6) by solving the bilinear equation (4) with z = y. It follows that the more general form of lump solution is also found to contain seven arbitrary independent parameters and four constraint conditions. On the basis of the general lump solution, we find the lumpoff and instanton/rogue wave solutions. The soliton is induced by the lump as a result of lump determines a soliton, which leads to the fact that the lump will be cut and then become lumpoff. Additionally, the lump becomes a rough wave when a pair of solitons are determined by the lump, and can only be visible at a specific point of time. It is possible for us to predict the exitance of the rogue wave thanks to the existence of the visible soliton.

Reference
[1] Solli D R Ropers C Koonath P Jalali B 2007 Nature 450 1054
[2] Müller P Garrett C Osborne A 2005 Oceanography 18 66
[3] Ma W X 2015 Phys. Lett. 379 1975
[4] Ma W X Zhou Y 2018 J. Differ. Equ. 264 2633
[5] X Chen S T Ma W X 2016 Nonlinear Dyn. 86 523
[6] Yang J Y Ma W X 2016 Int. J. Mod. Phys. 30 1640028
[7] Yang J Y Ma W X Qin Z Y 2018 Anal. Math. Phys. 8 427
[8] Wang H 2018 Appl. Math. Lett. 85 27
[9] X Tian B Qi F H 2012 Nonlinear Anal-Real 13 1130
[10] Manakov S V Zakharov V E Bordag L A Its A R Matveev V B 1977 Phys. Lett. 63 205
[11] Ma W X 2018 J. Geom. Phys. 133 10
[12] Ma W X Yong X Zhang H Q 2018 Comput. Math. Appl. 75 289
[13] Ma W X 2013 Rep. Math. Phys. 72 41
[14] Ma W X 2019 Mod. Phys. Lett. 33 1950457
[15] Ma W X 2019 J. Appl. Anal. Comput. 9 1319
[16] Ma W X 2019 Front. Math. Chin. 14 619
[17] Wang H F Zhang Y F 2019 Symmetry 11 1365
[18] Zhang Y Liu Y P Tang X Y 2018 Comput. Math. Appl. 76 592
[19] Zhang Y F Ma W X 2015 Z. Naturforsch. 70 263
[20] Tang X Y Lou S Y Zhang Y 2002 Phys. Rev. 66 046601
[21] Yang J Y Ma W X 2017 Nonlinear Dyn. 89 1539
[22] Shats M Punzmann H Xia H 2010 Phys. Rev. Lett. 104 104503
[23] Ginzburg N S Rozental R M Sergeev A S Fedotov A E Zotova I V Tarakanov V P 2017 Phys. Rev. Lett. 119 034801
[24] Moslem W M Shukla P K Eliasson B 2011 Europhys. Lett. 96 25002
[25] Höhmann R Kuhl U Stöckmann H J Kaplan L Heller E J 2010 Phys. Rev. Lett. 104 093901
[26] Satsuma J Ablowitz M J 1979 J. Math. Phys. 20 1496
[27] Osman M S Machado J A T 2018 Nonlinear Dyn. 93 733
[28] Sun Y Tian B Xie X Y Chai J Yin H M 2018 Wave. Random Complex 28 544
[29] Villarroel J Prada J Estévez P G 2009 Stud. Appl. Math. 122 395
[30] Liu J G Yang X J Cheng M H Fenf Y Y Wang Y D 2019 Comput. Math. Appl. 78 1947
[31] Akhmediev N Ankiewicz A Taki M 2009 Phys. Lett. 373 675
[32] Jia M Lou S Y 2018 arXiv:1803.01730v1[nlin.SI]
[33] Wu P X Zhang Y F 2019 Phys. Lett. 383 1755
[34] Mao J J Tian S F Zou L Zhang T T Yan X J 2019 Nonlinear Dyn. 95 3005
[35] Peng W Q Tian S F Zhang T T 2018 Phys. Lett. 382 2701
[36] Zheng P F Jia M 2018 Chin. Phys. 27 120201
[37] Gao L N Zhao X Y Zi Y Y Yu J X 2016 Comput. Math. Appl. 72 1225
[38] X Ma W X 2016 Nonlinear Dyn. 85 1217
[39] Hirota R 1971 Phys. Rev. Lett. 27 1192
[40] Hirota R 2004 Direct Method Soliton Theory Cambridge Cambridge University Press 155