Research Article | Open Access
A Generalized KdV Equation of Neglecting the Highest-Order Infinitesimal Term and Its Exact Traveling Wave Solutions
We study a generalized KdV equation of neglecting the highest order infinitesimal term, which is an important water wave model. Some exact traveling wave solutions such as singular solitary wave solutions, semiloop soliton solutions, dark soliton solutions, dark peakon solutions, dark loop-soliton solutions, broken loop-soliton solutions, broken wave solutions of U-form and C-form, periodic wave solutions of singular type, and broken wave solution of semiparabola form are obtained. By using mathematical software Maple, we show their profiles and discuss their dynamic properties. Investigating these properties, we find that the waveforms of some traveling wave solutions vary with changes of certain parameters.
In 1995, based on the physical and asymptotic considerations, Fokas  derived the following generalized KdV equation: which is an important water wave model, where , , , , , , , , . Regarding the as free parameters and using the to replace the (i.e., ), (1) becomes the following PDE: which is given by Tzirtzilakis et al. in ; they called it high-order wave equation of KdV type. Just as Tzirtzilakis  said these two equations are both water wave equations of KdV type, which are more physically and practically meaningful. The motion described by the model (1) or (2) is a 2-dimensional, inviscid, incompressible fluid (water) lying above a horizontal fiat bottom located at ( is a constant), and letting the air above the water. It turns out that, for such a system if the vorticity is zero initially, it remains zero. The fluids (water) analyzed by Fokas are only irrotational flows. In addition, this system is characterized by two parameters , with , , where and are two typical values of the amplitude and of the wavelength of the waves. The parameters , , satisfy the condition because the system is a model of short amplitude and long wavelength. Therefore the parameters and satisfy the condition , .
Assuming that the waves are unidirectional and neglecting terms of , (1) can be reduced to the classical KdV equation: In , Fokas assumed that is less than ; this implies . According to this assumption, we easily know that and . Neglecting two high-order infinitesimal terms of , (1) can be reduced to another high-order wave equation of KdV type [2–5] as follows: Equation (4) is a special case of (1) for . In , it was observed that (4) can be reduced by the local transformation of coordinates to a completely integrable PDE as follows: Equation (6) was first derived in  by using the method of bi-Hamiltonian systems, which based on the physical considerations, its Lax Pair was also given in .
If only neglecting the highest order infinitesimal term of , then (1) can be reduced to a new generalized KdV equation as follows: We call it a generalized KdV equation of neglecting the highest order infinitesimal term. In fact, (7) is another special case of (1) for ; it is also third-order approximate equation of KdV type. Of course, on describing dynamical behaviors of water waves, (4) is only a rough approximative model of (1) compared with (7); that is, the precision of model (7) is better than that of model (4) on describing dynamical behaviors of water waves. In other words, model (7) exhibits much richer phenomenology than model (4). Therefore, the investigation of exact traveling wave solutions for (7) are more practically meaningful than that of (4).
On the other hand, under the local transformation (1) can be reduced to the following generalized Gardner equation : where are certain expressions of . Clearly, (9) is equivalent to the following generalization form of the modified KdV equation [1, 8] under the transformation ,
From the above references and the references cited therein, we can know that (1) and (2) are very important water wave models. However, (1) and (2) are too complex to obtain their exact solution under universal conditions. Only under some special parametric conditions, their exact solutions were obtained in the existing literature. Next, let us briefly review the research backgrounds for the above equations.
In , Tzirtzilakis et al. only obtained two soliton-like solutions of (2) in the two groups of special conditions , , , , and , , , , . In , by using the planar bifurcation method of dynamical systems, under the four groups of special conditions , , , ; , , ; , , ; , , , , Li et al. studied (2); the existence of all kinds of traveling wave solutions were discussed completely, but its exact solutions were not investigated although some results of numerical simulation were obtained in this literature. In , in order to answer what is the dynamical behavior of one-loop soliton solution, Li studied the special case of (1) for . In , under different kinds of parametric conditions, Marinakis discussed two integrable cases for the third-order approximation model (1). In , Marinakis proved that (1) and some its of special cases are integrable. In , Gandarias and Bruzon proved that (1) is self-adjoint if and only if . In , by using the method as in , Li et al. studied (10); the existence of solitary wave, kink, and antikink wave solutions and uncountably infinite many smooth and nonsmooth periodic wave solutions was discussed except exact solutions. In , Bi also obtained some results of numerical simulation of (10); it is a pity that the exact traveling wave solutions of (10) were not obtained yet. In [4, 5], we obtained some exact traveling wave solutions of (4) under the parametric conditions or . In , under a new ansätze, Khuri studied (4); some exact solitary wave solutions and periodic wave solutions were obtained. In , by using method of planar dynamical system, Long et al. studied (6); the existence of smooth solitary wave and uncountably infinite many smooth and nonsmooth periodic wave solutions was proved in this literature.
From the above research backgrounds of (1) and (2), we can see that their exact solutions in universal conditions are hard to obtain because they are highly nonlinear equations and most probably they are not integrable equations in general. Thus, large numbers of research results are still concentrated in the classical KdV equation and some other high-order equations with KdV type, such as KdV-Burgers equation [17, 18] and KdV-Burgers-Kuramoto equation , at present. Therefore, the investigation of the more exact solutions for (1) is very important and necessary. However, by using the current methods, we can not obtain exact solutions of (1) in universal conditions; the next best thing is the investigation of exact solutions of (7). In this paper, still regarding the as free parameters and by using the integral bifurcation method [20, 21], we will investigate exact traveling wave solutions and their properties of (7).
The rest of this paper is organized as follows. In Section 2, we will derive two-dimensional planar system which is equivalent to (7) and give its first integral equation. In Section 3, by using the integral bifurcation method, we will obtain some new traveling wave solutions and discuss their dynamic properties.
2. Two-Dimensional Planar Dynamical System of (7) and Its First Integral and Conservation of Energy
Making a transformation with , (7) can be reduced to the following ODE: where is wave velocity which moves along the direction of -axis and . Integrating (11) once and setting the integral constant as zero yields Let . Thus (12) can be reduced to a planar system and a linear equation Obviously, the solutions of (12) cover the solutions of (13) and (14). We notice that the second equation in (13) is not continuous when ; that is, the function is not defined by . So (13) is a singular system; the line is called the singular line. Thus, we make the following transformation: where is a free parameter. Under the transformation (15), (13) can be rewritten as the following system Clearly, (16) is equivalent to (12). Except for the singular line , (13) and (16) have the same first integral.
First, from the parametric conditions in the model (1), we easily know that . Second, from [9–11, 13], we know that (1) contains many good properties including loop soliton, integrability, and self-adjoint property when . In addition, in the condition , it becomes easy when we solve (7). Therefore, we only consider the special case of (7) for in this paper. Thus, by using this assumption, we obtain the first integral of (13) and (16) as follows: where is an integral constant. For the convenience of discussion, taking the integral constant in (17) yields where denotes the particle's mass of system. Let , , . Thus, (18) can be rewritten as the following equation of conservation of energy: where denotes kinetic energy, denotes external energy, and denotes potential energy. In (19), the kinetic energy and potential energy are not conserved because the external energy , but the global energy () are subject to the energy conservation. if only if . Obviously, once ; under this case, (7) becomes the following equation: Similarly, (20) can be reduced to the following planar system: Obviously, the system (21) is a regular system; it has no singular line. Taking the integral constant as zero, we obtain (21)'s first integral as follows: Equation (22) can be rewritten as the following equation of conservation of energy: where are given above. The kinetic energy and potential energy are conserved in (23); in other words, the particle motion satisfies the conservation of kinetic energy and potential energy, which converts the kinetic energy from the potential energy then it converts the potential energy from the kinetic energy and go round and round. Therefore, (20) only has nonsingular traveling wave solutions, and all solutions are smooth; this is very different from (7).
In order to discuss singular or nonsingular traveling wave solutions of (7), we will consider (17). When , (17) can be reduced to Equation (24) can be further simplified under the following parametric conditions.
Case 1. One has .
Case 2. One has .
Case 3. One has .
Case 4. One has .
Case 5. One has or or .
Under the parametric conditions of Case 1, (24) can be reduced to Under the parametric conditions of Case 2, (24) can be reduced to Under the parametric conditions of Case 3, (24) can be reduced to Under the parametric conditions of Case 4, (24) can be reduced to Under the the parametric conditions of Case 5, (24) can be reduced to
3. Exact Traveling Wave Solutions of (7) and Their Dynamic Properties
3.1. The Exact Solutions under the Parametric Conditions of Case 1
Substituting (25) into the left expression (i.e., ) of (16) yields For the sake of convenience, in all the following discussions, we only discuss the case where the right of equation is “+” sign; the case of “−” sign can be similarly discussed. By the way, the solutions obtained by taking “+” sign and the solutions obtained by taking “−” sign are same when the is even function.
(i) When , completing the integrals (33) yields where , and with .
(ii) When , completing the integrals (33) yields Similarly, by using (41) and the transformation , we obtain four periodic solutions of (31) as follows: where ; that is, . Respectively, substituting (42) and (43) into (15) yields
Among the above traveling wave solutions, it is very worthy to mention solutions (40), (46), and (47); they have some peculiar dynamical properties and their traveling wave phenomena are very interesting. Solution (40) denotes a singular solitary wave; its shape is very similar to bright soliton, but it is not a normal soliton; the left parts of waveform degenerate into a crook. Solution (46) denotes a bounded wave with level asymptote; its shape is a half of a whole loop soliton; we call it semiloop soliton. In order to describe the dynamic properties of these two traveling wave solutions intuitively, we, respectively, plot profile of solutions (40) and (46) by using software Maple, which are shown in Figure 1.
(a) Singular solitary wave
(b) Semiloop soliton
It is very interesting that the solution (47) has six kinds of waveforms. Solution (47) with “+”, respectively, denotes a dark peakon, a dark loop soliton, and a broken loop soliton when the parameter decreases from to , which are shown in Figures 2(a)–2(d). Solution (47) with “−”, respectively, denotes a bright soliton and a singular compacton when the parameter decreases from to , which are shown in Figures 2(e) and 2(f). In other words, its waveforms depend on parameter extremity. Similar traveling wave phenomena that one solution contains multiwaveform first appeared in the investigation of the Degasperis-Procesi equation .
(a) Dark soliton
(b) Dark peakon
(c) Dark loop soliton
(d) Broken loop soliton (broken wave of C-form)
(e) Bright soliton
(f) Singular compacton
Figure 2(a) shows a shape of dark soliton under the fixed parameters , , , and . Figure 2(b) shows a shape of dark peakon under the fixed parameters , , , and . Figure 2(c) shows a shape of dark loop soliton under the fixed parameters , , , and . Similar traveling wave phenomena also appear in , in “Concluding remarks” of this literature; by using these phenomena, we successfully explain the movement of water waves. Figure 2(d) shows a shape of broken wave of C-form (upward hatch C) under the fixed parameters , , , and ; we also call it broken loop soliton or open-orbicular wave. Figure 2(e) shows a shape of bright soliton under the fixed parameters , , , and . Figure 2(f) shows a shape of singular compacton under the fixed parameters , , , and ; we call it singular compacton because it is not a normal compacton. Specifically, both sides of compacton wave become upward crook; it very much likes solitary waves but it is not a solitary wave after all because the left and right parts of waveform do not extend as .
3.2. The Exact Solutions under the Parametric Conditions of Case 2
(i) When , (49) can be reduced to By the Ferrari method or the Descartes method, it is easy to know that the equation has four roots as follows: where with . Especially, when , these four roots can be reduced to . For the convenience of discussion, we denote the above two real roots by the signs and always assume that the biggest root is denoted by . We also denote the above two complex roots by the signs . Thus (50) can be rewritten as where . Taking as the initial values and integrating (53), we obtain implicit solution of (7) as follows: where , , , , (ii) When , (49) can be reduced to Similarly, (56) can be rewritten as where are four roots of (51) and . Taking as the initial values then integrating (56), we obtain another implicit solution of (7) as follows: where the constants are given above, the , and
3.3. The Exact Solutions under the Parametric Conditions of Case 3
where with . When , (63) has another four roots .
(i) Particularly, when and , (62) can be reduced to
According to the above cases, we find that (63) still has two real roots and two complex roots; this case is very similar to that in (51) except their especial cases are different. Therefore, the types of their solutions are the same when and . Thus, we omit these parts of discussions. The especial case of (62) (i.e., (65)) can be reduced to where . Taking the as initial values and integrating (66), we obtain (7)'s periodic wave solution of implicit function as follows:
3.4. The Exact Solutions under the Parametric Conditions of Case 4
When , (71) can be reduced to Next, we discuss the roots of the following equation:
When , (73) has four real roots as follows:
When , (73) has still four real roots as follows: where is a double root.
When , (73) has two real roots and two complex roots as follows: Thus, by using the above parametric conditions, we can obtain three kinds of exact solutions of implicit function type; see the following discussions.
Under the conditions and , (72) can be reduced to Equation (79) can be rewritten as Taking as initial values and integrating (80), we obtain (7)'s solitary wave solution of implicit function type as follows: