Abstract
In the present paper, we have obtained an exact biperiodic, one-phase solution of the Kawahara evolution equation. Two classes of real periodic waves generated by the biperiodic solution have been analyzed. A modification of the bilinear-transformation method has been applied allowing to provide a single solution of the residual equation derived from the bidifferential reduction of the considered nonintegrable equation. It is shown that the spatial displacements are individual for each separate harmonic of the real periodic solutions.
1. Introduction
The nonlinear evolution equation was introduced by Kawahara [1] as a model describing one-dimensional waves for which the angle between the front and the gradient of the external field tends to the critical angle: . Here, , , are real parameters , is the elevation of dispersion medium, and , are the ion and electron masses, respectively. Within the zone of the critical angle , the coefficient in front of in (1.1) decreases causing disbalance between the nonlinear and dispersion effects, which is balanced by the member with the highest-order derivative. Therefore, three varieties of the dispersion medium are distinguished: medium with positive dispersion , ; medium with negative dispersion , and medium with mixed dispersion . Yamamoto [2] obtained the same evolution equation (1.1) in describing the dynamics of a nonviscous fluid in the vicinity of the critical depth. In this case, dispersion of the lowest order caused by the gravitation is balanced by the one caused by the surface tension.
The Kawahara evolution equation (1.1) is nonintegrable, which can be easily established by applying the modified criterion of Ablowtz and Segur [3]. Although it does not possess any conservation law and -soliton solutions , (1.1) has been a subject of research by a number of authors employing both analytic and numerical methods.
Kano and Nakayama [4] have found elliptic solution of (1.1), from which by means of appropriate phase modulations they obtained periodic cnoidal solutions, as well as M-type and W-type solitary waves. Yamamoto and Takizawa [5] have found a stable progressive-pulse solution which actually represents one-soliton impulse that can be generated only in a medium with positive or negative dispersion (i.e., for ). In the particular case , (1.1) is called FKdV (Five Korteweg-de Vries) and has been thoroughly studied numerically and analytically by Boyd [6]. Summarizing the numerical analysis accomplished by a number of authors [1, 2, 7], we can say that the Kawahara equation has two types of localized solitary waves: compressed and fictitiously scattered corresponding to a medium with negative and positive dispersion, respectively.
Studies on the existence of periodic solutions of (1.1) are considerably fewer. In 1997 Berloff and Howard [8], by applying the singular manifold method adapted for partial differential equations [9], obtained solitary and approximate periodic solutions of the Kawahara equation.
In the present paper, two families of real periodic solutions of the aforementioned equation are obtained. These exact periodic solutions were proved to be dynamically equivalent. A βspatialβ version [10] of the bilinear-transformation method of Hirota [11] and Matsuno [12] has been applied to derive them. In its classic form, this analytic model is inapplicable to nonintegrable partial differential equations.
2. Preliminaries
We could reduce the Kawahara equation in the so-called normal form if instead of the parameter we put , that is, , where , , and assume that
As will be seen below, the last condition does not lead to loss of generality since after the following rescaling: (1.1) is reduced to its normal form
Hence, it is easy to realize that for , we have a medium with negative dispersion, for along with replacing the variables , , , the medium possesses positive dispersion, while for , we have a medium with mixed dispersion, that is, we have all available variations of the dispersion medium.
The Kawahara equation (2.2) is invariant with respect to the Galilean transformation , , and , where = const, = const, that is, if is a solution of (2.2) with constant phase velocity , then the function is also a solution of the same equation but with increased phase velocity .
3. Biperiodic Solution
We will represent the solution of (2.2) by means of the Hirota-Satsuma [13] transformation where , are unknown parameters so far, and is unknown function, but on the assumption periodic and continuously differentiable to seventh order (with respect to ) in the semi-infinite domain . By substituting from (3.1) into (2.2) and employing the bidifferential identities for the Hirota operator (see Appendix A) we will obtain the following bidifferential form of this equation: where is an integration constant assumed to be nonzero. This constant plays a major role in constructing the periodic solutions though it does not have any dynamic features. If in (3.3) we put , then a sufficient condition for the function to be its solution is to satisfy the following two equations accordingly:
The first of these equations is called bidifferential (since its structure is such), and the second one will be called residual. We will search for the solution of (3.4) in the form where is the Jacobi biperiodic function [14] with phase variable . It is possible that the parameters , , are complex and , and (i.e., ) is the perturbation parameter. The function has a real period and an imaginary one . We substitute from (3.6) into the bidifferential (3.4). We obtain an infinite system of algebraic equations where and having .
The bilinear structure of (3.7) makes it possible to apply the index parity principle therein, which means that if in (3.8) we substitute , we will have the following relations: Summing up separately the even and odd addends in (3.7), we will obtain the following more compact form of this equation: where is the second Jacobi biperiodic function [14] defined by the equality . Accounting for the linear independence of the functions and in (3.10), it is reduced to two equations: and . In the context of equality (3.8) and the functional identities for the -functions, given in Appendix B, the last two equations are reduced to the following algebraic linear system: with respect to and . This system is compatible and definite since , where is the Wronskian determinant, and the wave number is different from zero. We have denoted for convenience , , so that all derivatives of these functions in the system (3.11) are with respect to the parameter . The solution of system (3.11) is as follows:
The residual equation (3.5) does not possess a bilinear structure, which is a typical feature for the nonintegrable partial differential equations. This means that we cannot apply the index parity principle to (3.5), that is, we cannot reduce the infinite system generated by this equation (after substituting with its equal function ) to a system of two equations. Consequently, we will represent the parameter in the formal series by the terms , which are real parameters unknown at this stage. Substituting (3.6) and (3.14) into the residual equation (3.5), we will obtain the infinite system
The infinite series on both sides of the system (3.15) are absolutely convergent for . It is easy to deduce that in summation term by term of the equalities from (3.15) for each whole number , the series (3.14) is an absolutely and uniformly convergent series under the same condition for the perturbation parameter. Thus, we obtain that the terms of the functional series (3.14) are determined for every integer by the formula
The terms of the uniformly convergent series (3.14) will be called spatial displacements. This will be explained further.
And finally in this section, we can make the conclusion that the continuously differentiable function is a family of localized biperiodic solutions of the Kawahara equation (2.2). For convenience, in (3.17), we have denoted
Figure 1 shows the discrete values of the spatial displacements for a mean value of the perturbation parameter , , , and . In addition to the foregoing, the function , expressed by (3.17), is a solution of the considered equation if the phase velocity satisfies the dispersion relation (3.12), and the integration constant is as in (3.13). The wave number is complex in the general case, but we have to exclude those values for which the constant is zero, that is,
4. Analicity Condition and Real Biperiodic Solutions
The obtained exact biperiodic solution (3.17) of the Kawahara evolution equation (2.2) is generally speaking a rational complex function of the phase variable . This function has twofold poles in the lattice of complex points , which are zeros of the function . As we are interested in the physical relevance of the obtained biperiodic solution, we can suitably choose the free parameters , , , so that we can achieve real biperiodic analogues of the solution and at the same time avoid the twofold poles. For this purpose, let us assume that
Under these conditions, to avoid the poles of the solution at the points , it is sufficient to limit the phase variable within the horizontal strip which is actually the analicity condition for the solution. To provide real biperiodic solutions from (3.17), we will consider two possible options for the wave number , under the condition of (3.19).(i)The wave number is real. Under this assumption and provided the condition of (4.1) holds, the phase variable is real for real phase shift . The logarithmic derivative from (3.17) in this case can be expressed in a Fourier series [14] by the equality Thus, from (3.17), we obtain a well-defined, real function which is periodical on the real straight line Figure 2 illustrates the forms of the real periodic solution for a real value of the wave number , (). It can be seen from the last formula that the spatial displacements contribute to each separate harmonic of the biperiodic solution.(ii)The wave number is imaginary. Let , and without loss of generality, we can assume that . Under this hypothesis, it is possible to generate real periodic solutions if the phase velocity is an imaginary number and the phase shift is properly chosen. It is obvious from the dispersion equality (3.12) that the phase velocity is an imaginary number (since contains only odd powers of ). For the phase shift , we choose , and hence . With this transition, the dispersion relation (3.12) remains the same. Applying the quasiperiodic property for , and from formula (3.17), we have for the logarithmic derivative (see [2]) In this case, we also obtain a real periodic function which is well defined in the strip (4.2):
Figure 3 visualizes the periodic solitary-wave profiles of the type for imaginary values of the wave number , and . Despite their structural differences, the two types of real biperiodic solutions of the Kawahara equation obtained in (4.4) and (4.7) are dynamically equivalent. This comes to confirm the nonlinear superposition principle established by Toda [15] for the KdV evolution equation. In conclusion, let us note that the biperiodic real solutions (4.4) and (4.7) actually describe the dynamics of dispersing waves resulting from the dispersion relation (3.12). Indeed, , which means that these waves are propagating in two directions, as periodically repeating groups traveling with group velocity .
5. Conclusions
Regardless of the purely technical difficulties arising in the application of the bilinear transformation method, it turns out to be a versatile theoretical instrument in analyzing nonintegrable, nonlinear partial differential equations, such as the Kawahara equation (1.1). The major intricacy arises from the residual equation after the bilinear reduction of the initial equation. Harmonization of the bilinear and the residual equation practically means to satisfy an algebraic system with infinite number of equations, but with a finite number of unknowns (usually one or two). The solution of a similar infinite system is feasible if a suitable unknown quantity is presented in the form of a convergent numerical or functional series with unknown termsβa procedure often applied in the Fourier method for linear partial differential equations. What is most surprising here is that the terms of these numerical or functional series have unambiguous physical interpretation. They describe spatial displacements individually for each harmonic. For this reason, we called this modification of the bilinear-transformation method, as well in [10], βspatial.β Note that for some nonintegrable nonlinear evolution equations, such as the Kuramoto-Sivashinsky equation spatial displacements are absent. In this case, their role is taken by the βwaveβ displacements, that is, the wave number is presented as a convergent infinite series.
Appendices
A. Logarithmic Derivatives Expressed by the Hirotaβs Bilinear Differential Operators ,
Consider the following:
B. Identities for the Jacobi -Functions
Consider the following:
Acknowledgment
The authors would like to thank the Academic Editors for the useful remarks that helped in improving the presentation of the paper.