Singular solutions of a modified two-component Camassa-Holm equation

The Camassa-Holm equation (CH) is a well known integrable equation describing the velocity dynamics of shallow water waves. This equation exhibits spontaneous emergence of singular solutions (peakons) from smooth initial conditions. The CH equation has been recently extended to a two-component integrable system (CH2), which includes both velocity and density variables in the dynamics. Although possessing peakon solutions in the velocity, the CH2 equation does not admit singular solutions in the density profile. We modify the CH2 system to allow dependence on average density as well as pointwise density. The modified CH2 system (MCH2) does admit peakon solutions in velocity and average density. We analytically identify the steepening mechanism that allows the singular solutions to emerge from smooth spatially-confined initial data. Numerical results for MCH2 are given and compared with the pure CH2 case. These numerics show that the modification in MCH2 to introduce average density has little short-time effect on the emergent dynamical properties. However, an analytical and numerical study of pairwise peakon interactions for MCH2 shows a new asymptotic feature. Namely, besides the expected soliton scattering behavior seen in overtaking and head-on peakon collisions, MCH2 also allows the phase shift of the peakon collision to diverge in certain parameter regimes.Comment: 25 pages, 11 figure

Self-Similar Blowup Solutions to the 2-Component Camassa-Holm Equations

In this article, we study the self-similar solutions of the 2-component Camassa-Holm equations% \begin{equation} \left\{ \begin{array} [c]{c}% \rho_{t}+u\rho_{x}+\rho u_{x}=0 m_{t}+2u_{x}m+um_{x}+\sigma\rho\rho_{x}=0 \end{array} \right. \end{equation} with \begin{equation} m=u-\alpha^{2}u_{xx}. \end{equation} By the separation method, we can obtain a class of blowup or global solutions for $\sigma=1$ or $-1$. In particular, for the integrable system with $\sigma=1$, we have the global solutions:% \begin{equation} \left\{ \begin{array} [c]{c}% \rho(t,x)=\left\{ \begin{array} [c]{c}% \frac{f\left( \eta\right) }{a(3t)^{1/3}},\text{ for }\eta^{2}<\frac {\alpha^{2}}{\xi} 0,\text{ for }\eta^{2}\geq\frac{\alpha^{2}}{\xi}% \end{array} \right. ,u(t,x)=\frac{\overset{\cdot}{a}(3t)}{a(3t)}x \overset{\cdot\cdot}{a}(s)-\frac{\xi}{3a(s)^{1/3}}=0,\text{ }a(0)=a_{0}% >0,\text{ }\overset{\cdot}{a}(0)=a_{1} f(\eta)=\xi\sqrt{-\frac{1}{\xi}\eta^{2}+\left( \frac{\alpha}{\xi}\right) ^{2}}% \end{array} \right. \end{equation} where $\eta=\frac{x}{a(s)^{1/3}}$ with $s=3t;$ $\xi>0$ and $\alpha\geq0$ are arbitrary constants.\newline Our analytical solutions could provide concrete examples for testing the validation and stabilities of numerical methods for the systems.Comment: 5 more figures can be found in the corresponding journal paper (J. Math. Phys. 51, 093524 (2010) ). Key Words: 2-Component Camassa-Holm Equations, Shallow Water System, Analytical Solutions, Blowup, Global, Self-Similar, Separation Method, Construction of Solutions, Moving Boundar

The inverse spectral transform for the conservative Camassa-Holm flow with decaying initial data

We establish the inverse spectral transform for the conservative Camassa-Holm flow with decaying initial data. In particular, it is employed to prove existence of weak solutions for the corresponding Cauchy problem.Comment: 27 page

Two-component {CH} system: Inverse Scattering, Peakons and Geometry

An inverse scattering transform method corresponding to a Riemann-Hilbert problem is formulated for CH2, the two-component generalization of the Camassa-Holm (CH) equation. As an illustration of the method, the multi - soliton solutions corresponding to the reflectionless potentials are constructed in terms of the scattering data for CH2.Comment: 22 pages, 3 figures, draft, please send comment

Self-Similar Blowup Solutions to the 2-Component Degasperis-Procesi Shallow Water System

In this article, we study the self-similar solutions of the 2-component Degasperis-Procesi water system:% [c]{c}% \rho_{t}+k_{2}u\rho_{x}+(k_{1}+k_{2})\rho u_{x}=0 u_{t}-u_{xxt}+4uu_{x}-3u_{x}u_{xx}-uu_{xxx}+k_{3}\rho\rho_{x}=0. By the separation method, we can obtain a class of self-similar solutions,% [c]{c}% \rho(t,x)=\max(\frac{f(\eta)}{a(4t)^{(k_{1}+k_{2})/4}},\text{}0),\text{}u(t,x)=\frac{\overset{\cdot}{a}(4t)}{a(4t)}x \overset{\cdot\cdot}{a}(s)-\frac{\xi}{4a(s)^{\kappa}}=0,\text{}a(0)=a_{0}% \neq0,\text{}\overset{\cdot}{a}(0)=a_{1} f(\eta)=\frac{k_{3}}{\xi}\sqrt{-\frac{\xi}{k_{3}}\eta^{2}+(\frac{\xi}{k_{3}}\alpha) ^{2}}% where $\eta=\frac{x}{a(s)^{1/4}}$ with $s=4t;$ $\kappa=\frac{k_{1}}{2}+k_{2}-1,$ $\alpha\geq0,$ $\xi<0$, $a_{0}$ and $a_{1}$ are constants. which the local or global behavior can be determined by the corresponding Emden equation. The results are very similar to the one obtained for the 2-component Camassa-Holm equations. Our analytical solutions could provide concrete examples for testing the validation and stabilities of numerical methods for the systems. With the characteristic line method, blowup phenomenon for $k_{3}\geq0$ is also studied.Comment: 13 Pages, Key Words: 2-Component Degasperis-Procesi, Shallow Water System, Analytical Solutions, Blowup, Global, Self-Similar, Separation Method, Construction of Solutions, Moving Boundary, 2-Component Camassa-Holm Equation

Integration of the EPDiff equation by particle methods

The purpose of this paper is to apply particle methods to the numerical solution of the EPDiff equation. The weak solutions of EPDiff are contact discontinuities that carry momentum so that wavefront interactions represent collisions in which momentum is exchanged. This behavior allows for the description of many rich physical applications, but also introduces difficult numerical challenges. We present a particle method for the EPDiff equation that is well-suited for this class of solutions and for simulating collisions between wavefronts. Discretization by means of the particle method is shown to preserve the basic Hamiltonian, the weak and variational structure of the original problem, and to respect the conservation laws associated with symmetry under the Euclidean group. Numerical results illustrate that the particle method has superior features in both one and two dimensions, and can also be effectively implemented when the initial data of interest lies on a submanifold

Integrable viscous conservation laws

We propose an extension of the Dubrovin-Zhang perturbative approach to the study of normal forms for non-Hamiltonian integrable scalar conservation laws. The explicit computation of the first few corrections leads to the conjecture that such normal forms are parameterized by one single functional parameter, named viscous central invariant. A constant valued viscous central invariant corresponds to the well-known Burgers hierarchy. The case of a linear viscous central invariant provides a viscous analog of the Camassa-Holm equation, that formerly appeared as a reduction of a two-component Hamiltonian integrable systems. We write explicitly the negative and positive hierarchy associated with this equation and prove the integrability showing that they can be mapped respectively into the heat hierarchy and its negative counterpart, named the Klein-Gordon hierarchy. A local well-posedness theorem for periodic initial data is also proven. We show how transport equations can be used to effectively construct asymptotic solutions via an extension of the quasi-Miura map that preserves the initial datum. The method is alternative to the method of the string equation for Hamiltonian conservation laws and naturally extends to the viscous case. Using these tools we derive the viscous analog of the Painlevé I2 equation that describes the universal behaviour of the solution at the critical point of gradient catastrophe