The continuation of solutions for the two-component Camassa-Holm system after wave breaking is studied in this paper. The global conservative solution is derived first, from which a semigroup and a multipeakon conservative solution are established. In developing the solution, a system transformation based on a skillfully defined characteristic and a set of newly introduced variables is used. It is the transformation, together with the associated properties, that allows for the establishment of the results for continuity of the solution beyond collision time.
Keywords:two-component Camassa-Holm system; Lagrangian system; global conservative solutions; conservative multipeakon solutions
Because of its capabilities of describing the dynamic behavior of water wave, the following Camassa-Holm (CH) equation
modeling the unidirectional propagation of shallow water waves in irrotational flow over a flat bottom, with representing the fluid velocity at time t in the horizontal direction, has attracted considerable attention [1-10]. The CH equation is a quadratic order water wave equation in an asymptotic expansion for unidirectional shallow water waves described by the incompressible Euler equations, which was found earlier by Fuchssteiner and Fokas  as a bi-Hamiltonian generalization of the KdV equation. It is completely integrable [2,3] and possesses an infinite number of conservation laws. A remarkable property of the CH equation is the existence of the non-smooth solitary wave solutions called peakons [2,11]. The peakon , , is smooth except at its crest and the tallest among all waves of fixed energy. Another remarkable fact for the CH equation is that it can model wave breaking [2,12], which means that the solution remains bounded while its slope becomes unbounded in finite time [12,13], setting it apart from the classical soliton equations such as KdV. After wave breaking, the solutions of the CH equation can be continued uniquely as either global conservative [4-6] or global dissipative solutions .
with and , (or in the ‘short wave’ limit, ), which is an extension of the CH equation by combining its integrability property with compressibility or free-surface elevation dynamics in its shallow water interpretation [11,17]. This system appeared originally in  as could be identified with the first negative flow of AKNS hierarchy, and then it was derived by Constantin and Ivanov  in the context of shallow water theory, with representing the horizontal velocity of the fluid and in connection with the free-surface elevation from equilibrium with the boundary assumptions and as . It is formally integrable [14-16] in the sense that it can be written as a compatibility condition of two linear systems (Lax pair) with a spectral parameter ζ:
It also has a bi-Hamiltonian structure corresponding to the Hamiltonian
and the Hamiltonian
The Cauchy problem for the two-component Camassa-Holm system has been studied extensively [18-24]. It was shown that the CH2 system is locally well posed with initial data , . The system also has global strong solutions which blow-up in finite time [19,21,22] and a global weak solution . However, the problem about continuation of the solutions beyond wave breaking, although interesting and important, has not been explicitly addressed yet. In our recent work , we studied the continuation beyond wave breaking by applying an approach that reformulated system (1.2) as a semilinear system of O.D.E. taking values in a Banach space. Such treatment makes it possible to investigate the continuity of the solution beyond collision time, leading to a global conservative solution where the energy is conserved for almost all times.
It should be stressed that both global conservation and multipeakon conservation are two important aspects worthy of investigation. To our best knowledge, however, little effort has been made in studying the multipeakon conservation associated with the CH2 system in the literature. As a compliment and extension to the previous work , we develop a novel approach in this work to construct the multipeakon conservative solution for the CH2 system. Different from the work , we reformulate the problem by utilizing a skillfully defined characteristic and a new set of variables, of which the associated energy serves as an additional variable to be introduced such that a well-posed initial-value problem can be obtained, making it convenient to study the dynamic behavior of wave breaking. Because of the introduction of the new variables, we are able to establish the multipeakon conservative solution from the global conservative solution for the CH2 system.
Some related earlier works [4,5] studied the global existence of solutions to the CH equation. However, the system considered in this work is a heavily coupled one, in which the mutual effect between the two components makes the analysis quite complicated and involved as compared with the system with a single component as studied in [4,5]. The key and novel effort made in this work to circumvent the difficulty is the utilization of the skillfully defined characteristic and the new set of variables, as well as careful estimates for each iterative approximate component of the solutions, which allows us to establish the global conservative solutions of system (1.2). It is shown that the multipeakon structure is preserved by the semigroup of a global conservative solution and the multipeakon solution is obtained by carefully computing the convolution equations and (), where, in contrast to the existing works, the inherent mutual effect between the two components is well reflected.
The remainder of this paper is organized as follows. Section 2 presents the transformation from the original system to a Lagrangian semilinear system. The global solutions of the equivalent semilinear system are obtained in Section 3, which are transformed into the global conservative solutions of the original system in Section 4. Finally, we establish the multipeakon conservative solutions for the original system in Section 5.
2 The original system and the equivalent Lagrangian system
We first present the original system. For simplicity, we consider here the associated evolution for positive times (of course, one would get similar results for negative times just by changing the initial condition into ). Let us introduce an operator , which can be expressed by its associated Green’s function such as for all . Thus, we can rewrite Eq. (1.2) as a form of a quasi-linear evolution equation:
then Eq. (1.2) can be rewritten as
Moreover, for regular solutions, we have that the total energy
We reformulate system (2.1) into a Lagrangian equivalent semilinear system as follows.
and define the Lagrangian cumulative energy distribution H as
It is not hard to check that
Then it follows from (2.3) and (2.5) that
Throughout the following, we use the notation
where we have dropped the variable t for simplicity and taken that y is an increasing function for any fixed time t for granted (the validity will be proved later). Using , we can rewrite and P in (2.7) as
From the definition of the characteristic, it follows that
System (2.10) can be regarded as an O.D.E. in the Banach space E given by
endowed with the norm
3 Global solutions of the equivalent system
Lemma 3.1 (See )
Proof We write
Since the operator Λ (defined as in Section 2) is linear and continuous from to , and is continuously embedded in , we have . It is not hard to know that R is locally Lipschitz from E into and therefore from E into . Thus, is locally Lipschitz from E to . Since the mapping is locally Lipschitz from E to W, it then follows from Lemma 3.1 that is locally Lipschitz from E to . Similarly, is also locally Lipschitz and therefore is locally Lipschitz from E to . We can obtain that P defined by (2.8) is locally Lipschitz continuous from E to in the same way. By using the chain rule, the formulas in (3.1) are obtained by direct computation, see [, p.129]. □
Then the theorem follows from the standard contraction argument on Banach spaces. □
Theorem 3.1 gives us the existence of local solutions to (2.10) for initial data in E. It remains to prove that the local solutions can be extended to global solutions. Note that the global solutions of (2.10) may not exist for all initial data in E. However, they exist when the initial data belongs to the set Γ, which is defined as follows.
The global existence of the solution for initial data in Γ relies essentially on the fact that the set Γ is preserved by the flow as the next lemma shows.
Note that . For , we take . For , we define . We consider U, P and as given functions in , which is guaranteed by Lemma 3.2 and V, N in . Thus system (2.11) is affine (it consists of a sum of a linear transformation and a constant) and, therefore, by using a contraction argument, it admits a unique local solution defined on some time interval . Thus, for the given initial condition , the solution of (2.10) given by Theorem 3.1 also belongs to , which implies that satisfies (3.5a) for all . We claim that (3.5c) holds for any and therefore almost everywhere. Consider a fixed and drop it in the notation. On the one hand, it follows from (2.11) that
and on the other hand,
Thus, . Notice that , which implies that for all . Thus, (3.5c) holds. It remains to prove that the inequalities in (3.5b) hold. Set . Assume that . Since is continuous w.r.t. t, we have . It follows from (3.5c) that . Furthermore, (2.11) implies and . If , then , which implies for all by the uniqueness of the solution of system (2.11). This contradicts the fact that for all . If , then . Since , there exists a neighborhood ϖ of such that for all . This contradicts the definition of . Hence, . We now have , which conversely implies for all , which contradicts the fact that . Thus, we have proved for all . We now prove that for all . This follows from (3.5c) when . If , then from (3.5c). As we have seen, would imply that for some in a punctured neighborhood of t, which is impossible. Hence, for all . Now we have for all . If for some , it then follows that , which implies for all , which contradicts the fact that for all . Hence, .
is a continuous semigroup.
We can obtain from the governing equation (2.10) that
which implies that
4 Global solutions for the original system
We transform the global solution of the equivalent system (2.10) into the global conservative solution of the original system (2.1) in this section. It suffices to establish the correspondence between the Lagrangian equivalent system and the original system.
We first introduce a set G as the set of relabeling functions defined by
with a useful property: If (), then almost everywhere. Conversely, if f is absolutely continuous, and there exists such that almost everywhere, then for some α depending only on c and . We now define the subsets F and of Γ such that
Notice that the map given by defines a group action of G on F, we then consider the quotient space of F w.r.t. the group action. The equivalence relation on F is defined as: for any , if there exists such that , we claim that X and are equivalent. We denote the projection by . For any , we introduce the map given by . It is not hard to prove that when , and for any and . Hence, we can define the map as for any representative of . For any , we have . Hence, . Note that any topology defined on is naturally transported into by this isomorphism, that is, if we equip with the metric induced by the E-norm, i.e., for all , which is complete, then for any , the topology on is defined by a complete metric given by .
For any initial data , we denote the continuous semigroup with the solution of system (2.10) by . As we indicated earlier, Eq. (2.1) is invariant w.r.t. relabeling. That is, , for any and . Thus, the map defined by is valid, which generates a continuous semigroup.
To derive the correspondence between the Lagrangian equivalent system and the original system, we have to consider the space D, which characterizes the solutions in the original system:
We first introduce the mapping M, which corresponds to the transformation from the Lagrangian equivalent system into the original system. In the other direction, we obtain the energy density μ in the original system, by pushing forward by y the energy density in the Lagrangian equivalent system, where the push-forward of a measure ν by a measurable function f is defined by
where , . We have that , which does not depend on the representative of we choose. We denote by the mapping to any and given by (4.1a) and (4.1b), which transforms the Lagrangian equivalent system into the original system.
We claim that the transformation from the original system into the Lagrangian equivalent system is a bijection.
Theorem 4.1The mapsMandLare invertible, that is,
Proof Let be given. We consider for a representative of and given by (4.1a) and (4.1b) for this particular X. From the definition of , we have . Let be the representative of in given by (4.2a)-(4.2c). To derive , it suffices to show that . Let
Using the fact that y is increasing and continuous, it follows that
For any given , using the fact that y is increasing and (4.4), it follows that . If , there then exists x such that and (4.6) implies that . Conversely, since y is increasing, we have , which implies that . This is a contradiction. Hence, we have that . Since , it follows directly from the definitions that , , and . Hence, .
By the definition of g, there exists an increasing sequence converging to such that . It follows from the definition of y in (4.2a) that . Passing to the limit, we obtain which, together with (4.8), yields
Our next task is to transport the topology defined in into D, which is guaranteed by the fact that we have established a bijection between the two equivalent systems and then obtained a continuous semigroup of solutions for the original system.
which makes the bijection L between D and into an isometry. Since equipped with is a complete metric space, it is not hard to know that D equipped with the metric is also a complete metric space. For each , we define the mapping as
Theorem 4.2Given, if we denotethe corresponding trajectory, thenis a weak solution of the two-component Camassa-Holm equations (2.1), which constructs a continuous semigroup. Moreover, μis a weak solution of the following transport equation:
Furthermore, we have
Thus, the unique solution described here is a conservative weak solution of system (2.1).
where is given by (2.1). We denote by the solution of (2.10) a representative of . On the one hand, since is Lipschitz and invertible w.r.t. ξ for almost all t, we can use the change of variables , then we get
It follows from the identity (3.5c) that
By comparing (4.15) and (4.16), we know that
Thus, the first identity in (4.13) follows directly from (4.14) and the second identity in (4.13) follows in the same way. It is not hard to check that is the solution of (4.10). From the definition μ in (4.1b), we can get that
which is constant in time from Lemma 3.3(iii). Thus, we have proved (4.11).
Thus, (4.12) holds and the proof is completed. □
5 Multipeakon solutions of the original system
We derive a new system of ordinary differential equations for the multipeakon solutions which is well posed even when collisions occur in this section, and the variables are used to characterize multipeakons in a way that avoids the problems related to blowing up.
Solutions of the two-component Camassa-Holm system may experience wave breaking in the sense that the solution develops singularities in finite time, while keeping the norm finite. Extending the solution beyond wave breaking imposes significant challenge as can be illustrated in the case of multipeakons given by
Peakons interact in a way similar to that of solitons of the CH equation, and wave breaking may appear when at least two of the coincide. Clearly, if the remain distinct, system (5.2) allows for a global smooth solution. In the case where has the same sign for all , the remain distinct, and (5.2) admits a unique global solution. In this case, the peakons are traveling in the same direction. However, when two peakons have opposite signs, collisions may occur, and if so, system (5.2) blows up.
Without loss of generality, we assume that the and are all nonzero, and that the are all distinct. From Theorem 4.2 we know that there exists a unique and global weak solution with initial data (5.3), and the aim is to characterize this solution explicitly. We consider the following characterization of multipeakons. The multipeakons are given as continuous solutions u defined on intervals as the solutions of the Dirichlet problem
Proof To prove this lemma, one proceeds as in Theorem 3.1 by using the contraction argument. The Banach space E is replaced by
endowed with the norm
for a constant C depending only on . From the derivative of given by (3.1) and Lemma 3.1, we have that is locally Lipschitz from into . Similarly, we obtain the same result for P. We compute the derivative of and on A as follows:
Since and are locally Lipschitz maps from into , we have that P and are locally Lipschitz from into . The local solution of (2.10) in then can be obtained by the standard contraction argument. As we know, as far as global existence is concerned, does not blow up with initial data in (see Lemma 3.3(i)). For any , we have that
We differentiate (3.5c) w.r.t. ξ and get
for some polynomial J. Since , we have X, and are in time. Since remains in Γ for all t, from (3.5b), we have and therefore is in time, which implies that is in time. For any time t such that , we have
For solutions with multipeakon initial data, we have the following result: If vanishes at some point in the interval , then vanishes everywhere in . Furthermore, for the given initial multipeakon solution , let be the solution of system (2.10) with initial data given by (5.4a)-(5.4e), then between adjacent peaks, if , the solution is twice differentiable with respect to the space variable, and for , we have that .
We now start the derivation of a system of ordinary differential equations for multipeakons.
At this point, we can get some more understanding of what is happening at the time of collision. Let be the time when the two peaks located at and collide, i.e., such that . The function remains continuous because the solution remains in for all time, thus we have . Still, , , and may have a finite limit when t tends to . However, the first derivative blows up, which implies and . Thus and tend to zero but slower than . In fact, if we let t tend to in (5.18), to first order in , we obtain
Inserting (5.19) into (5.14), we obtain
From (5.16) and (5.18), we get
It then follows from (5.21)-(5.26) that
Therefore, the above formulas (5.27)-(5.32) imply that
which can also be written in the following form
Now we can summarize the result as follows.
Theorem 5.2Let, andforwith a multipeakon initial datagiven by (5.3). Then with initial datathere exists a global solutionof (5.13), (5.33), (5.34). For each timet, is defined as the solution of the Dirichlet problemwith boundary conditions, on each interval. Thus, is a conservative solution of the two-component Camassa-Holm system, which is the multipeakon conservative solution.
The authors declare that they have no competing interests.
YS conceived of the multipeakon conservative solution of the two-component Camassa-Holm system and participated in its coordination and helped to draft the manuscript. YW participated in the design and coordination of the study, carried out the global conservative solution and the multipeakon conservative solution for the two-component Camassa-Holm system and drafted the manuscript. All authors read and approved the final manuscript.
The paper is supported by the Major State Basic Research Development Program 973 (No. 2012CB215202), the National Natural Science Foundation of China (No. 61134001) and the Fundamental Research Funds for the Central Universities (No. CDJXS12170003). The authors would like to thank the referees for constructive suggestions and comments.
Fuchssteiner, B, Fokas, AS: Symplectic structures, their Bäcklund transformations and hereditary symmetries. Physica D. 4, 47–66 (1981/1982). Publisher Full Text
Constantin, A: On the scattering problem for the Camassa-Holm equation. Proc. R. Soc. Lond., Ser. A, Math. Phys. Eng. Sci.. 457, 953–970 (2001). Publisher Full Text
Bressan, A, Constantin, A: Global conservative solutions of the Camassa-Holm equation. Arch. Ration. Mech. Anal.. 183, 215–239 (2007). Publisher Full Text
Holden, H, Raynaud, X: Global conservative solutions of the Camassa-Holm equation - a Lagrangian point of view. Commun. Partial Differ. Equ.. 32, 1511–1549 (2007). Publisher Full Text
Holden, H, Raynaud, X: Global conservative multipeakon solutions of the Camassa-Holm equation. J. Hyperbolic Differ. Equ.. 4, 39–64 (2007). Publisher Full Text
Bressan, A, Constantin, A: Global dissipative solutions of the Camassa-Holm equation. Appl. Anal.. 5, 1–27 (2007). Publisher Full Text
Pan, C, Jia, L, Cai, R, Ding, Y: Modeling and simulation of channel for underwater communication network. Int. J. Innov. Comput. Inf. Control. 8, 2149–2156 (2012). PubMed Abstract
Constantin, A, Lannes, D: The hydrodynamical relevance of the Camassa-Holm and Degasperis-Procesi equations. Arch. Ration. Mech. Anal.. 192, 165–186 (2009). Publisher Full Text
Constantin, A, Escher, J: Wave breaking for nonlinear nonlocal shallow water equations. Acta Math.. 181, 229–243 (1998). Publisher Full Text
Constantin, A: Global existence of solutions and breaking waves for a shallow water equation: a geometric approach. Ann. Inst. Fourier (Grenoble). 50, 321–362 (2000). Publisher Full Text
Chen, M, Liu, SQ, Zhang, Y: A 2-component generalization of the Camassa-Holm equation and its solutions. Lett. Math. Phys.. 75, 1–15 (2006). Publisher Full Text
Falqui, G: On a Camassa-Holm type equation with two dependent variables. J. Phys. A. 39, 327–342 (2006). Publisher Full Text
Constantin, A, Ivanov, R: On an integrable two-component Camassa-Holm shallow water system. Phys. Lett. A. 372, 7129–7132 (2008). Publisher Full Text
Gui, G, Liu, Y: On the Cauchy problem for the two-component Camassa-Holm system. Math. Z.. 268, 45–66 (2011). Publisher Full Text
Wang, Y, Huang, J: Global conservative solutions of the two-component Camassa-Holm shallow water system. Int. J. Nonlinear Sci.. 9, 379–384 (2010). PubMed Abstract
Guan, C, Yin, ZY: Global existence and blow-up phenomena for an integrable two-component Camassa-Holm shallow water system. J. Differ. Equ.. 248, 2003–2014 (2010). Publisher Full Text
Guan, C, Yin, Z: Global weak solutions for a two-component Camassa-Holm shallow water system. J. Funct. Anal.. 260, 1132–1154 (2011). Publisher Full Text