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:
where . If we define P as
then Eq. (1.2) can be rewritten as
Moreover, for regular solutions, we have that the total energy
is constant in time. Thus, Eq. (2.1) possesses the -norm conservation law defined as
where . Since , Young’s inequality ensures .
We reformulate system (2.1) into a Lagrangian equivalent semilinear system as follows.
Let denote the solution of system (2.1). For given initial data , we define the corresponding characteristic as the solution of
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
After the change of variables and , we obtain the following expressions for and P, namely
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
Let us introduce another variable such that (it will turn out that ). With these new variables, we now derive an equivalent system of equations (2.1),
where P and are given by (2.8). Differentiating (2.10) w.r.t. ξ yields
which is semilinear w.r.t. the variables , , , and .
System (2.10) can be regarded as an O.D.E. in the Banach space E given by
endowed with the norm
for any . Here W is a Banach space defined as
with the norm . Note that .
3 Global solutions of the equivalent system
In this section, we prove that the equivalent system admits a unique global solution. We first obtain the Lipschitz bounds we need on P and .
Lemma 3.1 (See )
Let and , or be two locally Lipschitz maps. Then the product is also a locally Lipschitz map fromEto , or fromEtoW.
Lemma 3.2For any given , Pand defined by (2.8) are locally Lipschitz continuous fromEto . Moreover, we have
Proof We write
where denotes the indicator function of a given set Ω, and , are the operators which correspond to the two terms of the last identity in (3.2). We rewrite as
where R is the operator from E to given by
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]. □
Theorem 3.1Let any be given. System (2.10) admits a unique local solution defined on some time interval , whereTdepends only on .
Proof To establish the local existence of solutions, one proceeds as in Lemma 3.2, then obtains that , which is defined by
with , is Lipschitz continuous on any bounded set of E. We rewrite the solutions of system (2.10) as
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.
Definition 3.1 The set Γ is composed of all such that
with , where .
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.
Lemma 3.3Given initial data , for some , we consider the local solution of system (2.10) given by Theorem 3.1. We have
(i) for all ,
(ii) for a.e. and a.e. ,
(iii) for all .
Proof (i) For given initial data , to ensure that the solution of (2.10) also belongs to , we have to specify the initial conditions for (2.11). Let Ω be the following set:
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, .
(ii) Define the set . It follows from Fubini’s theorem that
where and . From the above proof, we know that for all , consists of isolated points that are countable. This means that . It follows from (3.6), and since , that
This implies that for almost all t and therefore is strictly increasing and invertible w.r.t. ξ.
(iii) For any given , since and , we know that exist. We have
Let . Since U, P are bounded in and as , it then follows from (3.7) that for all . Since , it follows that for all . □
Theorem 3.2For any initial data , there exists a unique global solution for system (2.10). Moreover, for all , if we equip Γ with the topology endowed with theE-norm, then the map defined as
is a continuous semigroup.
Proof Let be a local solution of (2.10) with initial data . To obtain the global existence of solutions, it suffices to show that
Since is an increasing function w.r.t. ξ for all t and , we have . We consider a fixed and drop it for simplification. Since when , and for a.e. ξ, it follows from (3.5c) that
We can obtain from the governing equation (2.10) that
Thus, . The governing equation (2.10) also implies that and .
From the identity , we can deduce that
which implies that
Therefore, . Similarly, we obtain and the bounds hold for P. Let
After taking the -norms on both sides of (2.10) and (2.11), we obtain
It follows from Gronwall’s lemma that , which implies that is a continuous semigroup by the standard O.D.E. theory. □
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
where Id denotes the identity function. For any , we define the subsets of G as
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
With the above useful property of , it is not hard to prove that the space F is preserved by the governing equation (2.10).
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:
where and μ is a positive finite Radon measure with as its absolute continuous part.
We now establish a bijection between and D to transport the continuous semigroup obtained in the Lagrangian equivalent system (functions in ) into the original system (functions in D).
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
for all Borel set B. Let be defined as
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 are led to the mapping , which conversely transforms the original system into the Lagrangian equivalent system defined as follows.
Definition 4.1 For any , let
where . We define as the equivalence class of .
Remark 4.1 Note that , which satisfies (3.5a)-(3.5c) from the definition of y, U, V, N, H in (4.2a)-(4.2c). Moreover, by the definition (4.2c), we have that . Thus, .
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
and . From (4.1b) and since , for any , we get
Since and , we have
From the definition of , 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, .
Let be given and be the representative of in given by (4.2a)-(4.2c). Then, let . Let g be defined as before by (4.3). The same computation that leads to (4.5) now gives
Given , we consider an increasing sequence converging to , which is guaranteed by (4.2a), and such that . Let i tend to infinity. Since is lower semi-continuous, we have . Take and then we get
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
We obtain that by comparing (4.9) and (4.7). It is clear from the definitions that . Hence, and . □
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.
Let us define the distance on D as
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 denote the corresponding trajectory, then is 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).
Proof To prove that is a weak solution of the original system (2.1), it suffices to show that, for all with compact support,
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
Since and , it then follows from (2.10) that
On the other hand, using the change of variables and , and since y is an increasing function, we have
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).
Since a.e. for almost every , it then follows from (3.5c) that
for any Borel set B. Since y is one-to-one and almost everywhere, then (4.17) implies that
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
where satisfy the explicit system of ordinary differential equations
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.
We consider initial data given by
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
where the variables denote the position of the peaks, and the variables denote the values of u at the peaks. In the following, we will show that this property persists for conservative solutions.
Let us define as
which is a representative of in the Lagrangian equivalent system, that is, . Let . We claim that the functions , , and belong to and even belong to , as the next lemma shows.
Lemma 5.1For given initial data such that , the associated solution of (2.10) belongs to .
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
It suffices to show that P and are Lipschitz from bounded sets of into . Given a bounded set , where is a positive constant, it follows from Lemma 3.2 that
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
System (5.6) is affine w.r.t. , , , , . Hence, on any interval , we have
where C is a constant depending only on , which is bounded. Thus does not blow up from Gronwall’s lemma, and therefore the solution is globally defined in . □
Theorem 5.1Let the initial data be given in (5.3). The solution given by Theorem 4.2 satisfies between the peaks.
Proof Assuming that , we have
For a given , differentiating (5.8) w.r.t. t, we obtain, by using (2.10), (2.11) and (5.6), that
We differentiate (3.5c) w.r.t. ξ and get
After inserting the value of given by (5.10) into (5.9) and multiplying the equation by , we obtain that
It follows from (3.5c) and since that
We claim that is in time. Indeed, we have
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 some constant independent of time. This leads to
which corresponds to the conservation of spatial angular momentum. For the multipeakons at time , we have and for all . Hence,
for all time t and all . Thus, . □
For solutions with multipeakon initial data, we have the following result: If vanishes at some point in the interval , then