Bifurcation and travelling wave solutions for a (2+1)-dimensional KdV equation

This work aims to study a new equation that is recently introduced in (Phys. Lett. A. 383: 728–731, 2019). By using the method of dynamical systems, we examine the bifurcation and construct exact travelling wave solutions for a (2+1) KdV equation. Exact parametric representations of all wave solutions are introduced and they are clariﬁed graphically.


Introduction
The construction of the new travelling wave solutions for nonlinear evolution equations (NELLEs) is substantial and significant in various sides for the majority of phenomena physical and mathematics. The nonlinear wave phenomena arise in numerous branches of engineering and science, such as ocean engineering fluid mechanics, chemical kinematics, biology, solid-state physics, chemical physics, geochemistry, and meteorology [1][2][3][4]. The nonlinear wave phenomena of dissipation, reaction, diffusion, dispersion and convection are very meaningful in nonlinear wave equations. Thus, the search for the exact solutions for those equations has long been one of the fundamental topics of the perennial interest in mathematics and physics. In the procedures of seeking for exact solutions for those equations, numerous methods have been introduced and some of them have been developed, such as homogeneous balance method [5], Darboux transformations [6,7], extended tanh method [8,9], inverse scattering transform [10] and so on. It is well known that the traditional (1 + 1)KdV equation is a famous one that is utilized to characterize the waves on shallow water surfaces. Several authors have studied numerous versions of the KdV equation with distinct procedures and techniques from diverse points of view (see, e.g. [11][12][13][14][15][16][17][18]) and the references therein.
Recently in 2019, a new (2+1)-dimensional KdV equation has been proved based on the extended Lax pair and has been announced in [19]. It takes the following form: u t + 6u 2 (u y − u x ) − u xxx + u yy + 3u xxy − 3u xyy = 0. (2) To our knowledge, Equation (2) was not previously studied in other works from the point of view of constructing the travelling wave solutions, which may be helpful in understanding the physical interpretation for it. This motivates us to utilize the combination of the qualitative theory for differential equation and bifurcation for planar dynamical system associated with Equation (2) to construct new exact travelling wave solutions based on different values of the parameters.
Notice that, if u does not depend on x independently, i.e. u = u(y, t), Equation (2) reduces immediately u t + 6u 2 u y + u yyy = 0. ( Re-expressing the variable y as x we get the traditional mKdV equation. This paper is organized as follows: Section two contains a popular transformation which is utilized to convert the given time-evolution Equation (2) a planar conservative Hamiltonian system. Section three, we study the phase portrait for the deduced planar system in section two. In section four, based on the formula of the total energy corresponding to the deduced planar system, we construct several types of travelling wave solution and illustrate them graphically.

Travelling wave transformation
In this subsection, we apply the transformation: to Equation (2) to be able to construct a travelling wave solution, we obtain cu + 3λu 2 u +λu = 0 where denotes derivatives with respect to ξ while where k 1 and k 2 are arbitrary constants and c is the speed of the wave. Integrating Equation (5) with respect to ξ and equating the integration constant to zero, we obtain Equation (7) can be expressed as a system of first-order differential equations in the following form Preforming the transformation to the dynamical system (8), we obtain It is worth noting that the dynamical system (10) is a conservative Hamiltonian system that describes physically the motion of a particle in the Euclidean plane under the influence of the potential forces. The Hamiltonian function (Total energy) takes the following form where h is an arbitrary constant denoting the value of the energy. The Hamiltonian system (10) is onedimensional integrable system. Thus, the constant of motion, the total energy, is given by (11) and is used to find the explicit solution of the Hamiltonian system (10) and so, we can construct a travelling wave solution for Equation (2). For more details about the integrable systems with more than one dimension (see, e.g. [20][21][22][23]).

Bifurcation and phase portrait of Equation (2)
In this section, we are interested in investigating all the possible bifurcations and phase portrait for the dynamical system (10), taking into account the transformation (9). It is clear that the equilibrium points corresponding the Hamiltonian system (10) laying on the axis z = 0 are the zeros of f (u) = cu + λu 2 . Thus, the Hamiltonian system (10) has E 1 = (0, 0) as a unique equilibrium point if cλ ≥ 0, while if cλ < 0, it has three equilibrium points that are E 1 = (0, 0) and E 2,3 = ± √ (−c/λ), 0 . The determinant of the Jacobi matrix corresponding to the Hamiltonian system (10) admits the following form where (u 0 , 0) is any equilibrium point for the system (10). According to the theory of planar dynamical system, the equilibrium point (u 0 , 0) of the planar integrable dynamical system (10) The determinants of the Jacobi matrix (12) calculated at the equilibrium points Ei, i = 1, 2, 3 are We have three bifurcation curves: Those curves divide the space of the parameters into 10 regions that are defined as Taking into account the above restriction on the parameters and utilizing the software Maple for symbolic computations, we can present the phase portrait for the system (10) in the plane (u, z).
Case I: When (λ, c, λ) ∈ R 1 , the dynamical system (10) has a unique equilibrium point E 1 = (0, 0) and it is a saddle point because the determinant of the Jacobi matrix (14) is negative. The phase portrait for this case is clarified in Figure 1(a).
Case II: When (λ, c, λ) ∈ R 2 , E 1 = (0, 0) is a unique equilibrium point, and the determinant of the Jacobi matrix (14) is negative , so it is a saddle point. The phase portrait for this case is outlined in Figure 1(a).
Case III: When (λ, c, λ) ∈ R 3 , the point, E 1 = (0, 0) is a unique equilibrium point for a system (10). It is a centre since the determinant of the Jacobi matrix (14) is positive. The phase portrait for this case is illustrated in Figure 1 Case IV: When (λ, c, λ) ∈ R 4 , the equilibrium point E 1 is unique and it is a centre since the determinant of the Jacobi matrix (14) is positive. The phase portrait for this case is outlined in Figure 1 Case V: Assuming (λ, c, λ) ∈ R 5 , the system (10) has a unique equilibrium point E 1 = (0, 0) which is a centre point. The phase portrait for this case is clarified in Figure 1 Case VI: For(λ, c, λ) ∈ R 6 , the dynamical system (10) has a unique equilibrium point E 1 = (0, 0) that is a saddle point. The phase portrait for the present case is outlined in Figure 1(a).
Case VII: For (λ, c, λ) ∈ R 7 , the dynamical system (10) has three equilibrium points E 1 and E 2,3 . As a result of that, the determinant of the Jacobi matrix (14) is positive, E 1 it is a centre point. While the determinant of the Jacobi matrix (15) is negative, the equilibrium points E 2,3 are saddle points. The phase portrait is clarified in Figure 1(c).
Case VIII: If (λ, c, λ) ∈ R 8 , the dynamical system (10) has three equilibrium points E 1 and E 2,3 . The equilibrium point E 1 is a saddle point, while the other two are centres as it is outlined in Figure 1(d).
Case IX: When (λ, c, λ) ∈ R 9 the dynamical system has three equilibrium points E 1 is a saddle point due to the determinant of the Jacobi matrix (14) is negative, while the other two are centre points as a result of that, the determinant of the Jacobi matrix (15) is positive Figure 1(d).
Case X: If (λ, c, λ) ∈ R 10 , the dynamical system (10) has three equilibrium points, E 1 and E 2,3 and is a center point, while the other two equilibrium points are saddle points. The phase portrait for this case is outlined in Figure 1(c).

Remark 1:
We notice that there is more than a case can be described by the same figure that represents the phase portrait. These figures are not completely identical. Let us clarify that: For the case VII, the green colour curves define orbits in the phase space for the dynamical system (10) on a certain level of the energy h > (−c2/4λ), while the similar curve in case X represents the orbits for the dynamical system (10) on the level h < (−c2/4λ) of the energy. The red colour curves represent the orbits in the phase space for the system (10) on a zero level of the energy for both cases. For case VII, the blue curves represent the orbits in the phase space for the system (10) on a zero level of the energy, but similar curves represent the orbits of the system (10) on a level h > (−c2/4λ) of the energy for case X. Similar conclusion can be done for the other figures. Moreover, this will be illustrated in more detail in the following section.

Exact travelling wave solution
Assume that u(ξ ) is a continuous solution for the is named a solitary wave solution, and otherwise, it is called a kink (or anti-kink) wave solution. It is well known that homoclinic, heteroclinic and periodic orbits of the dynamical system (10) correspond to solitary, kink (or anti-kink) and periodic wave solutions, respectively. Therefore to examine all the bifurcations of solitary waves, kink waves and periodic waves for Equation (2), we must find all periodic annuli of the dynamical system (10) depending on the parameter space. The exact parametric representation for all bounded functions u(ξ ) is equivalent to construct exact travelling wave solutions for Equation (2). We utilize the energy integral (11) to find such solutions. It gives where P 4 (u) is a polynomial of degree four in u and it takes the following form Depending on the values of c,λ, λ and on a certain level of energy h, we can evaluate the integral in Equation (18) and construct a travelling wave solution for Equation (2).

Suppose (λ, c, λ) ∈ R 1
As it is outlined in Figure 1(a), there are three families of orbits for different values of the parameters h. The explicit representation for the travelling wave solution of Equation (2) is introduced by using Equation (18) for different values of the parameter h: • On a zero-level of the total energy (11), a family of orbits pass through the origin. This family appears with a red colour in Figure 1(a). Thus, the polynomial (19) is P 4 (u) = −λ 4 u 2 2c λ + u 2 . The expression (18) gives where δ is an arbitrary integration constant. The solution (20) is clarified in Figure 2(a).
• On a positive level of the energy h > 0, there is a family of orbits with blue colour as it is outlined in Figure 1(a). These orbits do not intersect the u-axis (z = 0) and so the polynomial (19) does not have real zeroes, i.e. P 4 (u) = − λ 4 (u 2 1 + u 2 )(u 2 2 + u 2 ). The expression (18) gives where sn(u, k), cn(u, k) and dn(u, k) are the Jacobi elliptic functions [24]. The travelling wave solution (21) is illustrated in Figure 2(b).
• On a negative level of the energy, the Hamiltonian system (10) has a family of orbits with green colour as it is clarified in Figure 1(a). These orbits intersect the u-axis in two points, so the polynomial (19) has two real zeroes, i.e. P 4 (u) = − λ 4 (u 2 − u 2 1 )(u 2 + u 2 2 ). The expression (18) implies the following.
The travelling wave solution (22) is illustrated in Figure 2(c).

. Suppose (λ, c, λ) ∈ R 2
There are three families of orbits for the Hamiltonian system (10). These families are defined by H(u, z) = h for different values of h. In Figure 1(a), the family of orbits with red colour is defined on a zero level of the constant h see, Figure 1(a) and the parametric representation of the travelling wave solution is given by equation (20) and it is clarified in Figure 2(a). On a negative level of h, the orbits with blue colour are defined by H(u, z) = h ∈]−∞, 0[, see, Figure 1(b). The parametric representation of the travelling wave solution is obtained by Equation (21) and it is outlined in Figure 2(b). The orbits with green colour correspond to H(u, z) = h, where h > 0, see, Figure 1(a). The parametric representation of the travelling wave solution reads as Equation (22) and it is illustrated by Figure 2(c).

Suppose (λ, c, λ) ∈ R 3
There are a family of periodic orbits around the origin that is defined by H(u, z) = h, where as it is outlined in Figure 1(b). The existence of this type of orbits refers to the existence of periodic travelling wave solutions. Any orbit of this family intersects the u-axis (z = 0), i.e. the polynomial (19) has two real zeroes and so where u 2 2 > u 2 1 . Consequently, the expression (18) leads to It is clear that the travelling wave solution (23) is periodic with period where K(k) is a complete elliptic integral of the first type. Figure 2(d) illustrates the solution (23).

Suppose (λ, c, λ) ∈ R 4
This case is similar to the previous case. The explicit representation of the travelling wave solution is given by (23) and it is clarified by Figure 2(d).

Suppose (λ, c, λ) ∈ R 5
There is a family of periodic orbits around the origin as it is outlined in Figure 1(b). This family corresponds to the curves defined by H(u, z) = h, where h > 0. Following Figure 1(b), any orbit of this family intersects the u-axis (z = 0) in two points, so the polynomial P 4 (u) has two real zeroes, i.e.
where u 2 1 = 4h λ . Consequently, the expression (18) leads to Notice that the solution (25) does not contain t explicitly due to c = 0 in the region R 5 . This type of solution is named static solution and it is is clarified by Figure 2(e).

Suppose (λ, c, λ) ∈ R 6
There is a family of orbits as outlined in Figure 1(a) with green colour. This family is defined by H(u, z) = h, where h < 0. For this case, any orbit of this family intersects the u-axis (z = 0) in two points and having a green colour as outlined in Figure 1(a). Thus, the polynomial (19) has two real zeroes and it takes the following form The expression (18) implies to The solution (26) represents a static solution for Equation (2). Figure 2(f) illustrates solution (26). (λ, c, λ)

Suppose
For this case, there are four families of orbits that are defined by the curves H(u, z) , as outlined by Figure 1(c). Let us study each case individually: • When h = − c 2 4λ , there is a family of orbits connecting the two saddle points and they are defined by H(u, z) = − c 2 4λ . These orbits are characterized by red colour as it is outlined by Figure 1(c). Any orbits of this family intersect the u-axis (z = 0) in two points that are saddle points, so the polynomial (19) has two repeated real zeroes and it takes the following form • When h = 0, there is a family of orbits that are outlined in Figure 1(c) with a blue colour. Any orbit of this family does not intersect the zaxis (u = 0). The polynomial (19) reads as P 4 (u) = − λ 4 u 2 u 2 + 2c λ and under the expression (18), the travelling wave solution for equation (2) admits the following form • On the level, h > − c 2 4λ , there is a family of orbits having a green colour in Figure 1(c). These orbits do not intersect the u-axis (z = 0) and the polynomial (19) has complex roots. The polynomial (19) has the form P 4 (u) = − λ 4 (u 2 + u 2 1 )(u 2 + u 2 2 ), where u 2 1,2 = a + ib. Using the expression (18), we obtain a travelling wave solution in the form The solution (29) is illustrated in Figure 2(i).
• On the level, 0 < h < − c 2 4λ , there is a family of periodic orbits around the origin. This family is characterized by H(u, z) = h and it is outlined by cyan colour in Figure 1(c). Any orbit of this family cuts the u-axis (z = 0) in four points and so the polynomial (18) has four real zeroes and so it becomes . The expression (18) gives an explicit travelling wave solution for Equation (2) in the following form

Suppose (λ, c, λ) ∈ R 8
For this case there are different families of orbits for the Hamiltonian system (10) and these orbits are defined by H(u, z) = h for different values of the constant h, where h > 0, h = 0, h ∈ −c 2 4λ , 0 . Let us study and examine each case individually: • On the level h > 0, there is a family of orbits which intersect the u-axis (z = 0) in two points, as outlined by red colour in Figure 1(d). Thus, the polynomial (19) has two real zeroes, so it becomes P 4 (u) = λ 4 (u 2 1 − u 2 )(u 2 2 + u 2 ). By virtue of expression (18), we can construct a travelling wave solution for Equation (2) as given by Equation (23) and it is outlined by Figure 2(d).
• On a zero level of the energy, there is an orbit passing through the origin which is a saddle point and returns to it again, as illustrated by blue curves in Figure 1(d). This type of orbits is named homoclinic orbit, which refers usually to the existence of solitary wave solution. It is clear from Figure  1(d), such orbit intersects the u-axis (z = 0) in three points and so, the polynomial (19) becomes P 4 (u) = λ 4 u 2 − 2c λ − u 2 . Thanks the expression (18), we obtain Solution (31) is clarified by Figure 2(k).
• For h ∈ −c 2 4λ , 0 , equation H(u, z) = h defines periodic orbits around the two equilibrium points E 2,3 . Any orbit of this family intersects the u− axis (z = 0) in four points, so the polynomial (19) has four real zeroes, so it can be written as P 4 (u) = λ 4 (u 2 1 − u 2 )(u 2 2 − u 2 ). The periodic travelling wave solution for this case can be obtained from Equation (30) by replacing λ by −λ and it is illustrated by Figure 2(j). (λ, c, λ)

Suppose
There are different families of orbits for the Hamiltonian system (10) that are defined by H(u, z) They are outlined in Figure 1(d) in different colours. Let us study them individually.
• There is a family of periodic orbits around the equilibrium points. These orbits are defined by H(uz) = h, where h < 0, and they are clarified by red colour in Figure 1(d). Any orbit of this family intersects the u-axis (z = 0) in two points. This clarifies the polynomial (18) has two real zeroes, so it reads as P 4 (u) = − λ 4 (u 2 − u 2 1 )(u 2 + u 2 2 ). Thanks to expression (18), the travelling wave solution for Equation (2) can be obtained from Equation (22) by replacing u 2 by u 1 . This solution is illustrated in Figure 2(c).
• On a zero level of the energy h, there is a heteroclinic orbit that is defined by H(u, z) = 0 and it is outlined by blue colour in Figure 1(d). The travelling wave solution for Equation (2) is given by Equation (28) and it is illustrated by Figure 2(h).
defines two families of periodic orbits around the two centre points E 2,3 and they are illustrated by Figure 1(d) by green colour. Figure 1(d) illustrates that the polynomial (19) has four real zeroes, so it can be rewritten as . Using expression (18), the travelling wave solution for Equation (2) takes the following form Solution (32) is clarified by Figure 2(l).

Suppose (λ, c, λ) ∈ R 10
There are different families of orbits for the Hamiltonian system (10) and they are defined by H(u, z) = h for different values of h as it is outlined in Figure 1(c) with different colours.
• On the level h = − c 2 4λ , the corresponding family of heteroclinic orbits is illustrated by Figure 1(c) with red colour. The travelling wave solution for equation (2) can be obtained by equation (27) and it is clarified by Figure 2(g).
• When h > 0, the travelling wave solution for Equation (2) can be obtained by Equation (22) and it is illustrated in Figure 2(c). • When h ∈ − c 2 4λ , 0 the travelling wave solution for Equation (2) is Equation (30) and it is clarified by Figure 2(j).
• When h ∈ −∞, −c 2 4λ , Equation (2) has the travelling wave solution (21) and it is clarified by Figure 2(b). (2) is regarded as an extension or generalization of the mKdV equation. Consequently, we can construct some exact travelling wave solutions for the mKdV equation (3) as special cases from the new solutions of equation (2). Let us clarify that with some details. To construct a travelling wave solution for the mKdV equation (3), we consider the change of variables (4) with k 1 = 0, k 2 = 1. Inserting those values in Equation (6), we get λ = 2 andλ = 1. Thus, we have

Remark 2: PDE equation
Based on the sign of the parameter c, we study the following three possible cases one by one:

Conclusion
We have studied a (2+1) KdV equation that has been recently introduced in 2019 [19]. This equation has been converted to a dynamical system by using certain transformation. The bifurcation and phase portrait for this system has been investigated. We have introduced some travelling wave solutions for a (2+1) KdV equation and those solutions have been graphically explained. Because a (2+1) KdV equation has been regarded as an extension or generalization for the mKdV equation, we can construct some exact travelling solution for it as special cases from the exact travelling solution for a (2+1) KdV equation for certain values of the included parameters.

Disclosure statement
No potential conflict of interest was reported by the authors.

Funding
This study was funded by the Deanship of Scientific Research at King Faisal University (grant number 17122010).