![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
Abstract
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 clarified graphically.
1. 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 [Citation1–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 [Citation5], Darboux transformations [Citation6,Citation7], extended tanh method [Citation8,Citation9], inverse scattering transform [Citation10] and so on. It is well known that the traditional equation
(1)
(1) is a famous one that is utilized to characterize the waves on shallow water surfaces. Several authors have studied numerous versions of the
equation with distinct procedures and techniques from diverse points of view (see, e.g. [Citation11–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 [Citation19]. It takes the following form:
(2)
(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 does not depend on
independently, i.e.
, Equation (2) reduces immediately
(3)
(3)
Re-expressing the variable as
we get the traditional
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.
2. Travelling wave transformation
In this subsection, we apply the transformation:
(4)
(4) to Equation (2) to be able to construct a travelling wave solution, we obtain
(5)
(5) where ′ denotes derivatives with respect to
while
(6)
(6) where
and
are arbitrary constants and
is the speed of the wave. Integrating Equation (5) with respect to
and equating the integration constant to zero, we obtain
(7)
(7) Equation (7) can be expressed as a system of first-order differential equations in the following form
(8)
(8) Preforming the transformation
(9)
(9) to the dynamical system (8), we obtain
(10)
(10) 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
(11)
(11) where
is an arbitrary constant denoting the value of the energy. The Hamiltonian system (10) is one-dimensional 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. [Citation20–23]).
3. 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 are the zeros of
Thus, the Hamiltonian system (10) has
as a unique equilibrium point if
, while if
, it has three equilibrium points that are
and
. The determinant of the Jacobi matrix corresponding to the Hamiltonian system (10) admits the following form
(12)
(12) where
is any equilibrium point for the system (10). According to the theory of planar dynamical system, the equilibrium point
of the planar integrable dynamical system (10) is either saddle point if
or centre point if
or cusp point if
besides the poincare index of this equilibrium point is zero. The values of the energy
at those equilibrium points are
(13)
(13) The determinants of the Jacobi matrix (12) calculated at the equilibrium points
are
(14)
(14)
(15)
(15) We have three bifurcation curves:
(16)
(16) Those curves divide the space of the parameters into 10 regions that are defined as
(17)
(17)
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
Case I: When , the dynamical system (10) has a unique equilibrium point
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 (a).
Figure 1. The level curves defined by for different values of
and the parameters
The solid black box indicates the equilibrium points.
![Figure 1. The level curves defined by Hu,dudξ=h for different values of h and the parameters c,λ,λ¯. The solid black box indicates the equilibrium points.](/cms/asset/a3c10ffc-12e6-49c1-b4a7-16421325754b/tusc_a_1709271_f0001_oc.jpg)
Case II: When ,
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 (a).
Case III: When , the point,
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 (b).
Case IV: When , the equilibrium point
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 (b)
Case V: Assuming , the system (10) has a unique equilibrium point
which is a centre point. The phase portrait for this case is clarified in Figure (b)
Case VI: For, the dynamical system (10) has a unique equilibrium point
that is a saddle point. The phase portrait for the present case is outlined in Figure (a).
Case VII: For , the dynamical system (10) has three equilibrium points
and
. As a result of that, the determinant of the Jacobi matrix (14) is positive,
it is a centre point. While the determinant of the Jacobi matrix (15) is negative, the equilibrium points
are saddle points. The phase portrait is clarified in Figure (c).
Case VIII: If , the dynamical system (10) has three equilibrium points
and
. The equilibrium point
is a saddle point, while the other two are centres as it is outlined in Figure (d).
Case IX: When the dynamical system has three equilibrium points
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 (d).
Case X: If , the dynamical system (10) has three equilibrium points,
and
and is a center point, while the other two equilibrium points are saddle points. The phase portrait for this case is outlined in Figure (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 , while the similar curve in case X represents the orbits for the dynamical system (10) on the level
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
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.
4. Exact travelling wave solution
Assume that is a continuous solution for the
equation for
and
When
, the solution
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
is equivalent to construct exact travelling wave solutions for Equation (2). We utilize the energy integral (11) to find such solutions. It gives
(18)
(18) where
is a polynomial of degree four in
and it takes the following form
(19)
(19) Depending on the values of
and on a certain level of energy
, we can evaluate the integral in Equation (18) and construct a travelling wave solution for Equation (2).
4.1. Suppose (![](//:0)
![](//:0)
As it is outlined in Figure (a), there are three families of orbits for different values of the parameters . The explicit representation for the travelling wave solution of Equation (2) is introduced by using Equation (18) for different values of the parameter
:
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 (a). Thus, the polynomial (19) is
The expression (18) gives
(20)
(20)
On a positive level of the energy h > 0, there is a family of orbits with blue colour as it is outlined in Figure (a). These orbits do not intersect the
-axis (
) and so the polynomial (19) does not have real zeroes, i.e.
The expression (18) gives
(21)
(21) where
and
are the Jacobi elliptic functions [Citation24]. The travelling wave solution (21) is illustrated in Figure (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 (a). These orbits intersect the u-axis in two points, so the polynomial (19) has two real zeroes, i.e.
The expression (18) implies the following.
(22)
(22)
The travelling wave solution (22) is illustrated in Figure (c).
4.2 . Suppose ![](//:0)
![](//:0)
There are three families of orbits for the Hamiltonian system (10). These families are defined by for different values of
. In Figure (a), the family of orbits with red colour is defined on a zero level of the constant
see, Figure (a) and the parametric representation of the travelling wave solution is given by equation (20) and it is clarified in Figure (a). On a negative level of h, the orbits with blue colour are defined by
, see, Figure (b). The parametric representation of the travelling wave solution is obtained by Equation (21) and it is outlined in Figure (b). The orbits with green colour correspond to
, where
, see, Figure (a). The parametric representation of the travelling wave solution reads as Equation (22) and it is illustrated by Figure (c).
4.3. Suppose ![](//:0)
![](//:0)
There are a family of periodic orbits around the origin that is defined by , where as it is outlined in Figure (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 (
), i.e. the polynomial (19) has two real zeroes and so
where
Consequently, the expression (18) leads to
(23)
(23)
It is clear that the travelling wave solution (23) is periodic with period
(24)
(24) where
is a complete elliptic integral of the first type. Figure (d) illustrates the solution (23).
4.4. Suppose ![](//:0)
![](//:0)
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 (d).
4.5. Suppose ![](//:0)
![](//:0)
There is a family of periodic orbits around the origin as it is outlined in Figure (b). This family corresponds to the curves defined by , where
. Following Figure (b), any orbit of this family intersects the u-axis (
) in two points, so the polynomial
has two real zeroes, i.e.
where
Consequently, the expression (18) leads to
(25)
(25) Notice that the solution (25) does not contain t explicitly due to
in the region
. This type of solution is named static solution and it is is clarified by Figure (e).
4.6. Suppose ![](//:0)
![](//:0)
There is a family of orbits as outlined in Figure (a) with green colour. This family is defined by , where
. For this case, any orbit of this family intersects the u-axis (
) in two points and having a green colour as outlined in Figure (a). Thus, the polynomial (19) has two real zeroes and it takes the following form
The expression (18) implies to
(26)
(26)
The solution (26) represents a static solution for Equation (2). Figure (f) illustrates solution (26).
4.7. Suppose ![](//:0)
![](//:0)
For this case, there are four families of orbits that are defined by the curves , where
is either
or
or
or
, as outlined by Figure (c). Let us study each case individually:
When
there is a family of orbits connecting the two saddle points and they are defined by
These orbits are characterized by red colour as it is outlined by Figure (c). Any orbits of this family intersect the u-axis
in two points that are saddle points, so the polynomial (19) has two repeated real zeroes and it takes the following form
Thanks to the expression (18), we get the explicit representation of the travelling wave solution in the following form
(27)
(27)
Figure (g) illustrates the travelling wave solution (27).
When
, there is a family of orbits that are outlined in Figure (c) with a blue colour. Any orbit of this family does not intersect the
axis (
). The polynomial (19) reads as
and under the expression (18), the travelling wave solution for equation (2) admits the following form
(28)
(28)
On the level,
there is a family of orbits having a green colour in Figure (c). These orbits do not intersect the u-axis (
and the polynomial (19) has complex roots. The polynomial (19) has the form
where
Using the expression (18), we obtain a travelling wave solution in the form
(29)
(29)
On the level,
there is a family of periodic orbits around the origin. This family is characterized by
and it is outlined by cyan colour in Figure (c). Any orbit of this family cuts the u-axis (
) 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
(30)
(30)
4.8. Suppose ![](//:0)
![](//:0)
For this case there are different families of orbits for the Hamiltonian system (10) and these orbits are defined by for different values of the constant
, where
,
Let us study and examine each case individually:
On the level h > 0, there is a family of orbits which intersect the u-axis (
) in two points, as outlined by red colour in Figure (d). Thus, the polynomial (19) has two real zeroes, so it becomes
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 (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 (d). This type of orbits is named homoclinic orbit, which refers usually to the existence of solitary wave solution. It is clear from Figure (d), such orbit intersects the u-axis (
) in three points and so, the polynomial (19) becomes
Thanks the expression (18), we obtain
(31)
(31)
For
equation
defines periodic orbits around the two equilibrium points
. Any orbit of this family intersects the
axis (
) in four points, so the polynomial (19) has four real zeroes, so it can be written as
The periodic travelling wave solution for this case can be obtained from Equation (30) by replacing
by
and it is illustrated by Figure (j).
4.9. Suppose ![](//:0)
![](//:0)
There are different families of orbits for the Hamiltonian system (10) that are defined by different values of the constant
They are outlined in Figure (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
where
, and they are clarified by red colour in Figure (d). Any orbit of this family intersects the u-axis (
) in two points. This clarifies the polynomial (18) has two real zeroes, so it reads as
Thanks to expression (18), the travelling wave solution for Equation (2) can be obtained from Equation (22) by replacing
by
. This solution is illustrated in Figure (c).
On a zero level of the energy h, there is a heteroclinic orbit that is defined by
and it is outlined by blue colour in Figure (d). The travelling wave solution for Equation (2) is given by Equation (28) and it is illustrated by Figure (h).
Equation
where
defines two families of periodic orbits around the two centre points
and they are illustrated by Figure (d) by green colour. Figure (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
(32)
(32)
4.10. Suppose ![](//:0)
![](//:0)
There are different families of orbits for the Hamiltonian system (10) and they are defined by for different values of h as it is outlined in Figure (c) with different colours.
On the level
the corresponding family of heteroclinic orbits is illustrated by Figure (c) with red colour. The travelling wave solution for equation (2) can be obtained by equation (27) and it is clarified by Figure (g).
When
, the travelling wave solution for Equation (2) can be obtained by Equation (22) and it is illustrated in Figure (c).
When
the travelling wave solution for Equation (2) is Equation (30) and it is clarified by Figure (j).
When
, Equation (2) has the travelling wave solution (21) and it is clarified by Figure (b).
Remark 2: PDE equation (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 Inserting those values in Equation (6), we get
and
Thus, we have
(33)
(33) Based on the sign of the parameter
, we study the following three possible cases one by one:
Case (I): For the mKdV equation (3) has a periodic travelling wave solution in the form (23) with
and
Case (II): When
we can construct a static solution for the mKdV equation in the form (25) due to
in this case, so the solution does not contain the time explicitly.
Case (III): If there are three possible cases depending on the parameter
is either positive or equals zero or lies between
and zero. For positive values of the parameter
the mKdV equation (3) has a solution as in case (I), while when
the travelling wave solution for the mKdV equation is given by (31) which is a sech-solution. If
the mKdV equation (3) possesses a solution in the form (30) with
and
is replaced by
These solutions have been shown in [Citation25].
5. Conclusion
We have studied a (2+1) KdV equation that has been recently introduced in 2019 [Citation19]. 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.
Additional information
Funding
References
- Johnson RS. Anon-linear equation incorporating damping and dispersion. J Fluid Mech. 1970;42:49–60. doi: 10.1017/S0022112070001064
- Younis M, Ali S. Solitary wave and shock wave solutions to the transmission line model for nano-ionic currents along microtubules. Appl Math Comput. 2014;246:460–463.
- Younis M, Rizvi STR, Ali S. Analytical and soliton solutions: nonlinear model of nanobioelectronics transmission lines. Appl Math Comput. 2015;265:994–1002.
- Razborova P, Moraru L, Biswas A. Perturbation of dispersive shallow water waves with Rosenau- KdV RLW equation and power law nonlinearity. Rom J Phys. 2014;59:7–8.
- Fan EG, Zhang HQ. A note on the homogeneous balance method. Phys Lett A. 1998;246:403–406. doi: 10.1016/S0375-9601(98)00547-7
- Wadati M, Sanuki H, Konno K. Relationships among inverse method. Bäcklund transformation and an infinite number of conservation laws. Prog Theor Phys. 1975;53:419–436. doi: 10.1143/PTP.53.419
- Matveev VB, Salle MA. Darboux transformation and Soliton. Berlin: Springer; 1991.
- Wazwaz AM. The tanh method for travelling wave solutions of nonlinear equations. Appl Math Comput. 2004;154:713–723.
- Wazwaz AM. The extended tanh method for abundant solitary wave solutions of nonlinear wave equations. Appl Math Comput. 2004;154:1131–1142.
- Garder CS, Green JM, Kruskal MD, et al. Method for solving the KdV equation. Phys Rev Lett. 1967;19:1095–1097. doi: 10.1103/PhysRevLett.19.1095
- Ablowitz MJ, Segur H. Solitons and the inverse scattering transform. Philadelphia: SIAM; 1981.
- Hirota R. The Direct method in Soliton theory. Japan: Osaka City University; 2004.
- Olver PJ. Application of Lie Group to differential equation. New York: Springer; 1986.
- Wazwaz AM. A KdV6 hierarchy: integrable members with distinct dispersion relations. Appl Math Lett. 2015;45:86–92. doi: 10.1016/j.aml.2015.01.014
- Geng X, Xue B. N-soliton and quasi-periodic solutions of the KdV6 equations. Appl Math Comp. 2012;219:3504–3510. doi: 10.1016/j.amc.2012.09.025
- Wazwaz AM, Xu GQ. An extended modified KdV equation and its Painlevé integrability. Nonlin Dyn. 2016;86:1455–1460. doi: 10.1007/s11071-016-2971-2
- Zhang Y, Dang XL, Xu HX. Bäcklund transformations and soliton solutions for the KdV6 equation. Appl Math Comp. 2011;217:6230–6236. doi: 10.1016/j.amc.2010.12.108
- Wen XY, Gao YT, Wang L. Darboux transformation and explicit solutions for the integrable sixth-order KdV equation for nonlinear waves. Appl Math Comp. 2011;218:55–60. doi: 10.1016/j.amc.2011.05.045
- Wang G, Kara AH. A (2+ 1)-dimensional KdV equation and mKdV equation: Symmetries, group invariant solutions and conservation laws. Phys Lett A. 2019;383:728–731. doi: 10.1016/j.physleta.2018.11.040
- Elmandouh AA. On the integrability of 2D Hamiltonian systems with variable Gaussian curvature. Nonlinear Dyn. 2018;93:933–943. doi: 10.1007/s11071-018-4237-7
- Elmandouh AA. New integrable problems in a rigid body dynamics with cubic integral in velocities. Results Phys. 2018;8:559–568. doi: 10.1016/j.rinp.2017.12.050
- Yehia HM, Elmandouh AA. A new conditional integrable case in the dynamics of a rigid body gyrostat. Mech Res Commun. 2016;78:25–27. doi: 10.1016/j.mechrescom.2016.09.007
- Elmandouh AA. On the integrability of the motion of 3D-Swinging Atwood machine and related problems. Phys Lett A. 2016;380:989–991. doi: 10.1016/j.physleta.2016.01.021
- Byrd PF, Fridman MD. Handbook of elliptic integrals for engineers and scientists. Berlin: Springer; 1971.
- Lv X, Shao T, Chen J. The study of the solution to a generalized KdV-mKdV equation. In Abst Appl Anal. 2013;2013:ID 249043.