Journal of Nonlinear Mathematical Physics

Volume 28, Issue 3, September 2021, Pages 321 - 336

Analytical Properties for the Fifth Order Camassa-Holm (FOCH) Model

Authors
Mingxuan Zhu1, Lu Cao2, Zaihong Jiang2, ORCID, Zhijun Qiao3, *, ORCID
1School of Mathematical Sciences, Qufu Normal University, Qufu 273100, P. R. China
2Department of Mathematics, Zhejiang Normal University, Jinhua, 321004, P. R. China
3School of Mathematical and Statistical Science, University of Texas Rio Grande Valley, Edinburg, TX 78539, USA
*Corresponding author. Email: zhijun.qiao@utrgv.edu
Corresponding Author
Zhijun Qiao
Received 21 March 2021, Accepted 6 May 2021, Available Online 16 June 2021.
DOI
10.2991/jnmp.k.210519.001How to use a DOI?
Keywords
The Fifth order Camassa-Holm (FOCH) model; global existence; infinite propagation speed; long time behavior
Abstract

This paper devotes to present analysis work on the fifth order Camassa-Holm (FOCH) model which recently proposed by Liu and Qiao. Firstly, we establish the local and global existence of the solution to the FOCH model. Secondly, we study the property of the infinite propagation speed. Finally, we discuss the long time behavior of the support of momentum density with a compactly supported initial data.

Copyright
© 2021 The Authors. Published by Atlantis Press B.V.
Open Access
This is an open access article distributed under the CC BY-NC 4.0 license (http://creativecommons.org/licenses/by-nc/4.0/).

1. INTRODUCTION

In this paper, we consider the following fifth order Camassa-Holm (FOCH) model [29]:

{mt+umx+buxm=0,t>0,x,m=(1-α2x2)(1-β2x2)u,t>0,x, (1.1)
where b ∈ ℝ is a constant, α, β ∈ ℝ are two parameters, αβ, αβ ≠ 0. Without loss of generality, we only consider the case α > 0, β > 0. When α < 0, β < 0, one can get the similar results by using the corresponding absolute values |α| and |β| instead of α, β.

In what follows, we present some mathematical results related to the topic of this paper. Liu and Qiao [29] obtained some interesting solutions including explicit single pseudo-peakons, two-peakon, and N-peakon solutions. Detailed dynamical interactions for two-pseudo-peakons and three-pseudo-peakons were also investigated in their paper with numerical simulations. There have been extensive studies on high order Camassa-Holm type equations in the mathematics physics fields [5,16,17,20,22,31,3841]. For the case α = β = 1, on the circle, McLachlan and Zhang [31] established the local well-posedness of the solution in Hs with s>72, it was shown that system (1.1) with α = β = 1 does’t admit finite time blow-up solutions. Tang and Liu [38] proved that the Cauchy problem for this equation is locally well-posed in the critical Besov space B2,17/2 or Bp,rs(1p,r+ands>max{3+1p,72}). The peakon-like solution and ill-posedness was also studied in [38]. For the case b = 2, m = uuxx + uxxxx, by using the Kato’s theory, the local well-posedness [39] was studied in the Sobolev space Hs with s>92. Ding [16,17] investigated the stationary solution, generality mild traveling solutions and conservative solution. Coclite, Holden and Karlsen [5] established the existence of global weak solutions. They also presented some invariant spaces under the action of the equation. In [20], the infinite propagation speed was considered for the case m = 4u − 5uxx + uxxxx. They also proved asymptotic behavior of the solution under the condition that the initial data decays exponentially and algebraically.

When β = 0 (or α = 0), it means m = uα2uxx. The Camassa-Holm equation, the Degasperis-Procesi equation, and the Holm-Staley b-family equations are the special cases of equation (1.1) with b = 2, b = 3 and b ∈ ℝ, respectively. These equations arise at various levels of approximation in shallow water theory, and possess a physics background with shallow water propagation, the bi-Hamiltonian structure, Lax pair, and explicit solutions including classical soliton, cuspon, and peakon solutions.

In 1993, Camassa and Holm [3] derived an integrable shallow water equation with peaked solitons, which was called the Camassa-Holm equation. In 1999, Degasperis and Procesi [15] extended the Camassa-Holm equation to a new water wave equation (Degasperis-Procesi equation). Both Camassa-Holm equation and Degasperis-Procesi equation have attracted much attention. They are completely integrable [11,12,14,35]. Infinitely many conservation laws have been shown in [14,27,37]. For the Camassa-Holm equation, in [9,28], They proved the local well-posedness for the initial datum in Hs with s > 3/2. There were many works to study the blow-up phenomenon, such as [810,24,28,32]. McKean [32] (See also [24] for a simple proof) proved that if and only if some portion of the positive part of y0 = u0u0xx lies to the left of some portion of its negative part, then the Camassa-Holm equation blow-up in finite time. The hierarchy properties, related finite-dimensional constrained flows, and algebro-geometric solutions of the Camassa-Holm equation were proposed in [34]. In [1], they studied the global conservative solution for the Camassa-Holm equation. Global dissipative solution have been shown in [2]. Constantin and Strauss [13] studied the orbital stability of the peakons. Himonas and his collaborators [21] obtained the persistence properties and unique continuation of solutions of the Camassa-Holm equation. In [25], the authors deduced the limit of the support of momentum density as t goes to +∞. In [4,6,7,23,26,30,33,35,36,42], they have investigated some mathematical properties for the Degasperis-Procesi equation. For the Holm-Staley b-family equation, mathematical studies have also been presented in [18,19,43].

The paper is organized as follows. In Section 2, we establish the local well-posedness and blow up scenario for the FOCH model. Conditions for global existence are found in Section 3. In Section 4, we establish the property of the infinite propagation speed for the FOCH model. In Section 5, we discuss the long time behavior for the support of momentum density of the FOCH model.

2. LOCAL WELL-POSEDNESS AND BLOW UP SCENARIO

Similar to the Camassa-Holm equation [9], we can establish the following local well-posedness theorem for the FOCH model (1.1).

Theorem 2.1.

Let u0Hs(ℝ) with s>72. Then there exist a T > 0 depending on ǁu0ǁHs, such that the FOCH model (1.1) has a unique solution

uC([0,T);Hs())C1([0,T);Hs-1()).

Morever, the map u0HsuC([0, T); Hs(ℝ)) ∩ C1([0, T); Hs−1(ℝ)) is continuous.

The proof is similar to that of Theorem 2.1 in [9,39]. To make the paper concise, we would like to omit the detail proof here. The maximum value of T in Theorem 2.1 is called the lifespan of the solution, in general. If T < ∞, that is

limtT-uHs=,
we say the solution blows up in finite time.

Before going to the blow up scenario, we have the following Lemma.

Lemma 2.2.

As m=(1-α2x2)(1-β2x2)u, then

u=p*m,p=α2α2-β2p1-β2α2-β2p2,
where p1=12αe-|x|α, p2=12βe-|x|β, αβ, α > 0, β > 0.

Proof. Taking fourier transform to m=(1-α2x2)(1-β2x2)u, we have

m^=(1+α2ξ2)(1+β2ξ2)u^.

Notice that when f(x) = ea|x|, a > 0 then

f^(ξ)=2aξ2+a2.

It follows that

u^(ξ)=11+α2ξ211+β2ξ2m^(ξ)=p^1(ξ)p^2(ξ)m^(ξ),
where p1=12αe-|x|α, p2=12βe-|x|β. Then,
u(x)=-1(p^1(ξ)p^2(ξ)m^(ξ))=p1*p2*m(x).

Let p = p1p2, we have

p(x)=14αβe-|x-y|αe-|y|βdy=12(α2-β2)[αe-|x|α-βe-|x|β]=α2(α2-β2)p1-β2(α2-β2)p2.

By Lemma 2.2, we can rewrite u(x, t) as

u=[α2(α2-β2)p1-β2(α2-β2)p2]*m=α2(α2-β2)[e-xα-xeξαm(ξ,t)dξ+exαx+e-ξαm(ξ,t)dξ]-β2(α2-β2)[e-xβ-xeξβm(ξ,t)dξ+exβx+e-ξβm(ξ,t)dξ]. (2.1)

Differentiating u with respect to x, we have

ux=12(α2-β2)[-e-xα-xeξαm(ξ,t)dξ+exαx+e-ξαm(ξ,t)dξ]+12(α2-β2)[e-xβ-xeξβm(ξ,t)dξ-exβx+e-ξβm(ξ,t)dξ].

Then, we present the precise blow-up scenario.

Theorem 2.3.

Assume that u0H4(ℝ) and let T be the maximal existence time of the solution u(x, t) to equation (1.1), αβ, α > 0, β > 0 with the initial data u0(x).

  1. (1).

    If b>12, then the corresponding solution of the FOCH model (1.1) blows up in finite time if and only if

    limtTinfx{ux(x,t)}=-.

  2. (2).

    If b<12, then the corresponding solution of the FOCH model (1.1) blows up in finite time if and only if

    limtTsupx{ux(x,t)}=+.

Proof. By direct calculation, we have

mL22=[u-(α2+β2)uxx+α2β2uxxxx]2dx=u2+(α2+β2)2uxx2-2(α2+β2)uuxx+α4β4uxxxx2+2α2β2uuxxxx-2(α2+β2)α2β2uxxuxxxxdx=u2+(α2+β2)2uxx2+2(α2+β2)ux2+α4β4uxxxx2+2α2β2uxx2+2(α2+β2)α2β2uxxx2dx.

Hence

cuH42mL22CuH42,
where c and C are positive constants depending on α and β. If b>12, direct calculation we have
ddtm2dx=(1-2b)uxm2dx(1-2b)infx{ux(x,t)}m2dx.

If

infx{ux(x,t)}-M,
then
ddtm2dx-(1-2b)Mm2dx.

By using the Gronwall inequality,

mL22=m2dxe-(1-2b)Mm02dx=e-(1-2b)Mm0L22.

Therefore the H4-norm of the solution is bounded on [0, T).

On the other hand,

u=α2(α2-β2)p1*m-β2(α2-β2)p2*m=α2(α2-β2)p1(x-ξ)m(ξ)dξ-β2(α2-β2)p2(x-ξ)m(ξ)dξ.

By the Sobolev’s embedding ‖ux ≤ ‖uH4, it tells us if H4-norm of the solution is bounded, then the L-norm of the first derivative is bounded.

By the same argument, we can get the similar result for b<12. So, we omit the details and complete the proof of Theorem 2.3.

3. GLOBAL EXISTENCE

In this section, we study the global existence. Before going to our main results, we give the particle trajectory as

{qt=u(q,t),0<t<T,x,q(x,0)=x,x, (3.1)
where T is the lifespan of the solution. Taking derivative (3.1) with respect to x, we obtain
dqtdx=qtx=ux(q,t)qx,t(0,T).

Therefore

{qx=exp{0tux(q,s)ds},0<t<T,x,qx(x,0)=1,x,
which is always positive before the blow-up time. Therefore, the function q(x, t) is an increasing diffeomorphism of the line before blow-up. In fact, direct calculation yields
ddt(m(q)qxb)=[mt(q)+u(q,t)mx(q)+bux(q,t)m(q)]qxb=0.

Hence, we have the following identity

m(q)qxb=m0(x),0<t<T,x. (3.2)

Theorem 3.1.

Assume that u0H4(ℝ), αβ, α > 0, β > 0, if b=12 or b = 2, then the corresponding solution of FOCH model (1.1) will exist globally in time.

Remark 3.1.

If α = 0 or β = 0, system (1.1) reduce to the well-known b-family equation. The global existence for b=12 and Theorem 3.2 can be reduce to the results for b-family equation [18]. The global existence for b = 2 is the new discovery compared to the b-family equation.

Proof. Let

E(t)=u2+(α2+β2)ux2+α2β2uxx2dx.

Differentiating E(t), we have

ddtE(t)=2uut+2(α2+β2)uxuxt+2α2β2uxxuxxtdx=2uut-2(α2+β2)uuxxt+2α2β2uuxxxxtdx=2umtdx=(b-2)u2mxdx.

It yields that E(t) = E(0) when b = 2. By the Sobolev’s imbedding, we have

uxLuH22CE(t)=CE(0).

The global existence for b = 2 is completed by Theorem 2.3. Applying m on (1.1) and integration by parts, we obtain

ddtm2dx=-2buxm2+mmxudx=-2buxm2-m22uxdx=(1-2b)uxm2dx.

If b=12, then ddtm2dx=0. Hence,

uH42mL22=m0L22.

It follows that the corresponding solution of FOCH model (1.1) exists globally when b=12.

Theorem 3.2.

Supposing that u0H4, αβ, α > 0, β > 0, m0=(1-α2x2)(1-β2x2)u0 does not change sign. Then the corresponding solution to (1.1) exists globally.

Proof. We can assume that m0 ≥ 0. It is sufficient to prove ux is bounded for all t. In fact,

ux=12(α2-β2)[exαx+e-ξαm(ξ,t)dξ-e-xα-xeξαm(ξ,t)dξ]+12(α2-β2)[e-xβ-xeξβm(ξ,t)dξ-exβx+e-ξβm(ξ,t)dξ].

If m0 ≥ 0, α > β > 0, then

ux=12(α2-β2)[exαx+e-ξαm(ξ,t)dξ-e-xα-xeξαm(ξ,t)dξ]+12(α2-β2)[e-xβ-xeξβm(ξ,t)dξ-exβx+e-ξβm(ξ,t)dξ]12(α2-β2)[exαxe-ξαm(ξ,t)dξ+e-xβ-xeξβm(ξ,t)dξ]12(α2-β2)[exαe-xαxm(ξ,t)dξ+e-xβexβ-xm(ξ,t)dξ]12(α2-β2)[m(ξ,t)dξ+m(ξ,t)dξ]=1(α2-β2)m0(ξ,t)dξ
and
ux=12(α2-β2)[exαxe-ξαm(ξ,t)dξ-e-xα-xeξαm(ξ,t)dξ]+12(α2-β2)[e-xβ-xeξβm(ξ,t)dξ-exβx+e-ξβm(ξ,t)dξ]12(α2-β2)[-e-xα-xeξαm(ξ,t)dξ-exβx+e-ξβm(ξ,t)dξ]-1(α2-β2)m(ξ,t)dξ=-1(α2-β2)m0(ξ,t)dξ.

That is

|ux|1(α2-β2)m0(ξ,t)dξ.

If m0 ≥ 0, 0 < α < β, then

ux=12(α2-β2)[exαxe-ξαm(ξ,t)dξ-e-xα-xeξαm(ξ,t)dξ]+12(α2-β2)[e-xβ-xeξβm(ξ,t)dξ-exβx+e-ξβm(ξ,t)dξ]12(α2-β2)[-e-xα-xeξαm(ξ,t)dξ-exβx+e-ξβm(ξ,t)dξ]12(α2-β2)[-e-xαexα-xm(ξ,t)dξ-exβe-xβxm(ξ,t)dξ]12(β2-α2)[m(ξ,t)dξ+m(ξ,t)dξ]=1(β2-α2)m0(ξ,t)dξ
and
ux=12(α2-β2)[exαxe-ξαm(ξ,t)dξ-e-xα-xeξαm(ξ,t)dξ]+12(α2-β2)[e-xβ-xeξβm(ξ,t)dξ-exβx+e-ξβm(ξ,t)dξ]12(α2-β2)[exαxe-ξαm(ξ,t)dξ+e-xβ-xeξβm(ξ,t)dξ]-1(β2-α2)m(ξ,t)dξ=-1(β2-α2)m0(ξ,t)dξ.

That is

|ux|1(β2-α2)m0(ξ,t)dξ.

When m0 ≤ 0, via the similar approach that is used above, we could also obtain the global existence result. So, we omit the details and complete the proof of Theorem 3.2.

4. INFINITE PROPAGATION SPEED

The main theorem reads as follows:

Theorem 4.1.

Assume that the initial datum u0(x) ∈ H4(ℝ) is compactly supported in [a, c], then for t ∈ (0, T), the corresponding solution u(x, t) to the FOCH model (1.1) αβ, α > 0, β > 0 has the following property:

u(x,t)={α2(α2-β2)e-xαE1(t)-β2(α2-β2)e-xβE2(t),asx>q(c,t),α2(α2-β2)exαF1(t)-β2(α2-β2)exβF2(t),asx<q(a,t),
where
E1(t)=exαm(x,t)dx,F1(t)=e-xαm(x,t)dx,
and
E2(t)=exβm(x,t)dx,F2(t)=e-xβm(x,t)dx,
denote continuous nonvanishing functions.

Furthermore, if α > 0, 0<β32α, 0bmin{3-2β2α2,53}, E1(t) is strictly increasing function, while F1(t) is strictly decreasing function.

Similarly, if β > 0, 0<α32β, 0bmin{3-2α2β2,53}, E2(t) is strictly increasing function, while F2(t) is strictly decreasing function.

Remark 4.1.

Theorem 4.1 implies that the strong solution u(x, t) doesn’t have compact x-support for any t > 0 in its lifespan, although the corresponding u0(x) is compactly supported.

Proof. Since u0(x) has a compact support in the interval [a, c], so does m0(x)=(1-α2x2)(1-β2x2)u0(x). Equation (3.2) tells us that m(x)=(1-α2x2)(1-β2x2)u(x) is compactly supported in the interval [q(a, t), q(c, t)] well-defined in its lifespan. Hence the following functions are well-defined

E1(t)=exαm(x,t)dx,F1(t)=e-xαm(x,t)dx,E2(t)=exβm(x,t)dx,F2(t)=e-xβm(x,t)dx.

Using (3.2),

m(q(x,t),t)0,x<aorx>c,
we know
u(x,t)=(α2(α2-β2)p1-β2(α2-β2)p2)*m(x,t)=α2(α2-β2)e-|x-ξ|αm(ξ)dξ-β2(α2-β2)e-|x-ξ|βm(ξ)dξ=α2(α2-β2)q(a,t)q(c,t)e-|x-ξ|αm(ξ)dξ-β2(α2-β2)q(a,t)q(c,t)e-|x-ξ|βm(ξ)dξ.

Then, for x > q(c, t), we have

u(x,t)=α2(α2-β2)q(a,t)q(c,t)e-x-ξαm(ξ)dξ-β2(α2-β2)q(a,t)q(c,t)e-x-ξβm(ξ)dξ=α2(α2-β2)e-xαq(a,t)q(c,t)eξαm(ξ)dξ-β2(α2-β2)e-xβq(a,t)q(c,t)eξβm(ξ)dξ=α2(α2-β2)e-xαE1(t)-β2(α2-β2)e-xβE2(t). (4.1)

Similarly, when x < q(a, t), we have

u(x,t)=α2(α2-β2)q(a,t)q(c,t)ex-ξαm(ξ)dξ-β2(α2-β2)q(a,t)q(c,t)ex-ξβm(ξ)dξ=α2(α2-β2)exαq(a,t)q(c,t)e-ξαm(ξ)dξ-β2(α2-β2)exβq(a,t)q(c,t)e-ξβm(ξ)dξ=α2(α2-β2)exαF1(t)-β2(α2-β2)exβF2(t). (4.2)

On the other hand,

dE1(t)dt=eξαmt(ξ,t)dξ.

It is easy to get

mt=-mx-bmux=[(α2+β2)uxxx-α2β2uxxxxx-ux]u-b[u-(α2+β2)uxx+α2β2uxxxx]ux=(α2+β2)uxxxu-α2β2uxxxxxu-(b+1)uux+b(α2+β2)uxxux-bα2β2uxxxxux. (4.3)

Taking (4.3) into dE1(t)dt, we obtain

dE1(t)dt=exαmtdx=exα[(α2+β2)uxxxu-α2β2uxxxxxu-(b+1)uux+b(α2+β2)uxxux-bα2β2uxxxxux]dx=(α2+β2)exαuxxxudx-α2β2exαuxxxxxudx-(b+1)exαuuxdx+b(α2+β2)exαuxxuxdx-bα2β2exαuxxxxuxdx=i=15Ii. (4.4)

I1-I5 can be estimated as follows:

I1=(α2+β2)exαuxxxudx=-α2+β22α3exαu2dx+3(α2+β2)2αexαux2dx, (4.5)
I2=-α2β2exαuxxxxxudx=β22α3exαu2dx-5β22αexαux2dx+5αβ22exαuxx2dx, (4.6)
I3=-(b+1)exαuuxdx=b+12αexαu2dx, (4.7)
I4=b(α2+β2)exαuxxuxdx=-b(α2+β2)2αexαux2dx, (4.8)
I5=-bα2β2exαuxxxxuxdx=bβ22αexαux2dx-3bαβ22exαuxx2dx. (4.9)

Combining (4.5)(4.9) to (4.4), we have

dE1(t)dt=exαmtdx=b2αexαu2dx-(b-3)α2+2β22αexαux2dx+(5-3b)αβ22exαuxx2dx. (4.10)

For α > 0, 0<β32α, 0bmin{3-2β2α2,53}, from (4.10), E1(t) is strictly increasing for nontrivial solution.

Similary,

dF1(t)dt=e-xαmtdx=-b2αe-xαu2dx+(b-3)α2+2β22αe-xαux2dx-(5-3b)αβ22e-xαuxx2dx. (4.11)

For α > 0, 0<β32α, 0bmin{3-2β2α2,53}, from (4.11), F1(t) is strictly decreasing for nontrivial solution.

dE2(t)dt=exβmtdx=b2βexβu2dx-(b-3)β2+2α22βexβux2dx+(5-3b)βα22exβuxx2dx. (4.12)

For β > 0, 0<α32β, 0bmin{3-2α2β2,53}, from (4.12), E2(t) is strictly increasing for nontrivial solution.

dF2(t)dt=e-xβmtdx=-b2βe-xβu2dx+(b-3)β2+2α22βe-xβux2dx-(5-3b)βα22e-xβuxx2dx. (4.13)

For β > 0, 0<α32β, 0bmin{3-2α2β2,53}, from (4.13), F2(t) is strictly decreasing for nontrivial solution.

This complete the proof of Theorem 4.1.

Remark 4.2.

Let

u(x,t)={α2(α2-β2)e-xαE1(t),asx>q(c,t),α2(α2-β2)exαF1(t),asx<q(a,t),u(x,t)={β2(α2-β2)e-xβE2(t),asx>q(c,t),β2(α2-β2)exβF2(t),asx<q(a,t),

We rewrite u = u′u″, as consequences of (4.1) and (4.2), we have

u(x,t)=-αux(x,t)=α2uxx(x,t)=α2(α2-β2)e-xαE1(t),asx>q(c,t),u(x,t)=αux(x,t)=α2uxx(x,t)=α2(α2-β2)exαF1(t),asx<q(a,t).
and
u(x,t)=-βux(x,t)=β2uxx(x,t)=β2(α2-β2)e-xβE2(t),asx>q(c,t),u(x,t)=βux(x,t)=β2uxx(x,t)=β2(α2-β2)exβF2(t),asx<q(a,t).

Theorem 4.2.

Suppose the initial value u0(x) ∈ H4(ℝ), m0=(1-α2x2)(1-β2x2)u0, α > β > 0 m0 doesn’t change sign onand u0 has compact support in the interval [a, c]. Then for t ∈ (0, T), the corresponding solution u(x, t) of equation (1.1) satisfies

12(α+β)e-xα|E1(t)|u(x,t)α2(α2-β2)e-xα|E1(t)|,asx>q(c,t),12(α+β)exα|F1(t)|u(x,t)α2(α2-β2)exα|F1(t)|,asx<q(a,t).
where
E1(t)=eξαm(ξ,t)dξ,F1(t)=e-ξαm(ξ,t)dξ,
denote continuous nonvanishing functions.

Remark 4.3.

We assume α > β > 0 to get the above conclusion in Theorem 4.2, because the position of α, β is symmetric, then β > α > 0, we have results similar to the above conclusions about E2(t)=eξβm(ξ,t)dξ, F2(t)=e-ξβm(ξ,t)dξ,

12(α+β)e-xβ|E2(t)|u(x,t)β2(β2-α2)e-xβ|E2(t)|,asx>q(c,t),12(α+β)exβ|F2(t)|u(x,t)β2(β2-α2)exβ|F2(t)|,asx<q(a,t).

Theorem 4.2 can be seem as a generalization of the result in [20]. Comparing with Theorem 4.1, it show more detailed estimation by adding the additional condition on m0.

Proof. If u0 has a compact support set [a, c], then the corresponding m0 also has a corresponding compact support set [a, c]. It is known from (3.2) that m has the same compact support set [q(a, t), q(c, t)]. We define

u1=(1-βα)α2(α2-β2)e-|x-ξ|αmdξ=12(α+β)e-|x-ξ|αmdξ,u2=α2(α2-β2)e-|x-ξ|αmdξ. (4.14)

According to E1(t)=eξαm(ξ,t)dξ,F1(t)=e-ξαm(ξ,t)dξ, then

u1(x,t)=12(α+β)e-xαE1(t),u2(x,t)=α2(α2-β2)e-xαE1(t),asx>q(c,t),u1(x,t)=12(α+β)exαF1(t),u2(x,t)=α2(α2-β2)exαF1(t),asx<q(a,t).

According to (4.1) and (4.14), we obtain

u2(x,t)-u(x,t)=β2(α2-β2)e-|x-ξ|βmdξ,u(x,t)-u1(x,t)=β2(α2-β2)(e-|x-ξ|α-e-|x-ξ|β)mdξ.

Then, we obtain

{u1(x,t)u(x,t)u2(x,t),m00,u2(x,t)u(x,t)u1(x,t),m00.

If m0 ≥ 0,

{12(α+β)e-xαE1(t)u(x,t)α2(α2-β2)e-xαE1(t),asx>q(c,t),12(α+β)exαF1(t)u(x,t)α2(α2-β2)exαF1(t),asx<q(a,t). (4.15)

If m0 ≤ 0,

{α2(α2-β2)e-xαE1(t)u(x,t)12(α+β)e-xαE1(t),asx>q(c,t),α2(α2-β2)exαF1(t)u(x,t)12(α+β)exαF1(t),asx<q(a,t). (4.16)

The proof of Theorem 4.2 is finished.

5. LONG TIME BEHAVIOR FOR THE SUPPORT OF MOMENTUM DENSITY

After the global existence of solution is established, we will discuss the long time behavior for the support of momentum density of the FOCH model. Now, we give the lemma and main theorem as follows:

Lemma 5.1.

Let α > β > 0, Assume the initial value u0 ≢ 0 has a compact supported set [a, c].

  1. (1).

    If m0(x) ≥ 0(≢ 0), x ∈ [a, c], then we have

    limt+F1(t)=0.

  2. (2).

    If m0(x) ≤ 0(≢ 0), x ∈ [a, c], then we have

    limt+E1(t)=0.

Remark 5.1.

By the same argument, we can get a similar conclusion for β > α > 0. If m0(x) ≥ 0(≢ 0), x ∈ [a, c], then we have

limt+F2(t)=0.

If m0(x) ≤ 0(≢ 0), x ∈ [a, c], then we have

limt+E2(t)=0.

Proof. (1) For m0(x) > 0, from (3.2), we have E1(t) > 0, F1(t) > 0, E2(t) > 0, F2(t) > 0, for all t ≥ 0. As F1(t) > 0, we claim that

limt+F1(t)=0.

Otherwise, there is a constant 0 > 0, for any T > 0, there will exist a t > T, such that F1(t) ≥ 0.

For any d < a, from (4.15) we have

ddtq(d,t)=u(q(d,t),t)12(α+β)eq(d,t)αF1(t)12(α+β)eq(d,t)α0.

It follows that

e-q(d,t)α-02α(α+β)t+e-dαα.

Taking T=2(α+β)0e-dα, however, when t = T + 1,

-02α(α+β)t+e-dαα<0,

This is the contradiction. So our claim is right.

(2). For m0(x) < 0, from (3.2), we have E1(t) < 0, F1(t) < 0, E2(t) < 0, F2(t) < 0, for all t ≥ 0. As F1(t) > 0, we claim that

limt+E1(t)=0.

Otherwise, there is a constant 0 > 0, for any T > 0, for any T > 0, there will exist a t > T, such that E1(t) ≤ −0.

For any h > c, from (4.16) we have

ddtq(h,t)=u(q(h,t),t)12(α+β)e-q(h,t)αE1(t)-02(α+β)e-q(h,t)α.

It follows that

eq(h,t)α-02α(α+β)t+e-dαα.

Taking T=2(α+β)0e-dα, however, when t = T + 1,

-02α(α+β)t+e-dαα<0,

This is the contradiction. So our claim is right.

Theorem 5.2.

If b > 1, α > β > 0, and suppose that m0(x)L1b and u0(x) has a compact supported set [a, c].

  1. (1).

    If m0(x) ≥ 0(≢ 0), x ∈ [a, c], then we have

    eq(c,t)α(b-1)-eq(a,t)α(b-1)+,ast+. (5.1)

  2. (2).

    If m0(x) ≤ 0(≢ 0), x ∈ [a, c], then we have

    e-q(a,t)α(b-1)-e-q(c,t)α(b-1)+,ast+. (5.2)

Remark 5.2.

For the case β > α > 0, by using the properties of E2 and F2 in Remark 5.1, one can replace α with β in (5.1) and (5.2).

Proof. (1) By (3.2) and direct calculation, we have

(ac(m0)1bdx)b=(ac(m(q,t)qxb)1bdx)b=(ac(m(q,t))1bqxdx)b=(q(a,t)q(c,t)(m(ξ,t))1bdξ)b(q(a,t)q(c,t)m(ξ,t)e-ξαdξ)[q(a,t)q(c,t)eξα(b-1)dξ]b-1=F1(t)[α(b-1)(eq(c,t)α(b-1)-eq(a,t)α(b-1))](b-1).

It follows

[α(b-1)(eq(c,t)α(b-1)-eq(a,t)(b-1)α)](b-1)(ac(m0)1bdx)bF1(t).

Using the limit

limt+F1=0,
we can get
eq(c,t)α(b-1)-eq(a,t)α(b-1)+,ast+.

(2). Direct calculation, we have

(ac(-m0)1bdx)b=(ac(-m0)1bdx)b=(ac(-m(q,t)qxb)1bdx)b=(ac(-m(q,t))1bqxdx)b=(q(a,t)q(c,t)(-m(ξ,t))1bdξ)b(q(a,t)q(c,t)(-m(ξ,t)eξα)dξ)[q(a,t)q(c,t)e-ξα(b-1)dξ]b-1=-E1(t)[α(b-1)(e-q(a,t)α(b-1)-e-q(c,t)α(b-1))]b-1.

It follows

[α(b-1)(e-q(a,t)α(b-1)-e-q(c,t)α(b-1))]b-1(ac(-m0)1bdx)b-E1(t).

Using the limit

limt+E1=0,
we can obtain
e-q(a,t)α(b-1)-e-q(c,t)α(b-1)+,ast+.

Theorem 5.3.

If b = 1, suppose that m0(x) ∈ L1 and u0(x) has a compact supported set [a, c].

  1. (1).

    If m0(x) ≥ 0(≢ 0), x ∈ [a, c], then we have

    q(c,t)+,ast+.

  2. (2).

    If m0(x) ≤ 0(≢ 0), x ∈ [a, c], then we have

    q(a,t)-,ast+.

Proof. We only present the proof for α > β > 0. The case β > α > 0 can be proved by the same argument. (1) As m0(x) ≥ 0, for any t ≥ 0, we have F1(t)> 0. According to Lemma 5.1, we know

limt+F1(t)=0.

Direct calculation, we have

acm0dx=acm(q,t)qxdxeq(c,t)αq(a,t)q(c,t)m(ξ,t)e-ξαdξ=eq(c,t)αF1(t).

It follows

eq(c,t)αacm0dxF1(t)+,ast+,
then we can get
q(c,t)+,ast+.

(2). As m0(x) ≤ 0, for any t ≥ 0, we have E1(t) < 0. According to Lemma 5.1, we know

limt+E1(t)=0.

Direct calculation, we have

ac(-m0)dx=ac(-m(q,t)qx)dxe-q(a,t)αq(a,t)q(c,t)(-m(ξ,t))eξαdξ=-e-q(a,t)αE1(t).

It follows

e-q(a,t)αac(-m0)dx-E1(t)+,ast+,
then we can get
q(a,t)-,ast+.

Theorem 5.4.

If 0 < b < 1, α > β > 0, m0(x)L1b or b = 0, m0L. Suppose that u0(x) has a compact supported set [a, c].

  1. (1).

    If m0(x) ≥ 0(≢ 0) for x ∈ [a, c], then we have

    eq(c,t)α-eq(a,t)α+,ast+. (5.3)

  2. (2).

    If m0(x) ≤ 0(≢ 0) for x ∈ [a, c], then we have

    e-q(a,t)α-e-q(c,t)α+,ast+. (5.4)

Remark 5.3.

For the case β > α > 0, by using the properties of E2 and F2 in Remark 5.1, one can replace α with β in (5.3) and (5.4).

Proof. (1). For m0(x) ≥ 0, we have F1(t) > 0 for all t ≥ 0. From Lemma 5.1, we know

limt+F1(t)=0.

According to the conservation law

mdx=m0dx,m1bdx=m01bdx.If0<b<1and{γ+ηb=1,2γ+η=1,0<γ,η<1.{0<η=22-b-1<1,0<γ=1+1b-2<1.

By direct calculation, we obtain

m0dx=mdx=m(q,t)qxdx=[ac(me-qαqx)γ(m1bqx)η(eqαqx)γdx](acme-qαqxdx)γ(acm1bqxdx)η(aceqαqxdx)γ=(q(a,t)q(c,t)me-ξαdξ)γ(m1bdξ)η(q(a,t)q(c,t)eξαdξ)γ=(F1(t))γ(m1bdξ)η(q(a,t)q(c,t)eξαdξ)γ=(F1(t))γ(m01bdξ)η(αeq(c,t)α-αeq(a,t)α)γ.

It follows

(αeq(c,t)α-αeq(a,t)α)γm0dx(F1(t))γ(m01bdξ)η+,
then we can obtain
eq(c,t)α-eq(a,t)α+,ast+.

If b = 0, we can obtain

m0dx=limb0mdx=m(q,t)qxdx=limb0[ac(me-qαqx)γ(m1bqx)η(eqαqx)γdx]limb0(acme-qαqxdx)γ(acm1bqxdx)η(aceqαqxdx)γ=limb0(q(a,t)q(c,t)me-ξαdξ)γ(m1bdξ)η(q(a,t)q(c,t)eξαdξ)γ=limb0(F1(t))γ(m1bdξ)η(q(a,t)q(c,t)eξαdξ)γ=limb0(F1(t))γ(m01bdξ)η(αeq(c,t)α-αeq(a,t)α)γ.

It follows

(αeq(c,t)α-αeq(a,t)α)γm0dx(F1(t))γ(limb0m01bdξ)η+,
then we can obtain
eq(c,t)α-eq(a,t)α+,ast+.

(2). For m0(x) ≤ 0, we have E1(t) < 0 for all t ≥ 0. From Lemma 5.1, we know

limt+E1(t)=0.

Similarly, according to the conservation law

mdx=m0dx,m1bdx=m01bdx.If0<b<1and{γ+ηb=1,2γ+η=1,0<γ,η<1.{0<η=22-b-1<1,0<γ=1+1b-2<1.

By direct calculation, we obtain

-m0dx=-mdx=-m(q,t)qxdx=[ac(-meqαqx)γ((-m)1bqx)η(e-qαqx)γdx](ac-meqαqxdx)γ(ac(-m)1bqxdx)η(ace-qαqxdx)γ=(q(a,t)q(c,t)-meξαdξ)γ((-m)1bdξ)η(q(a,t)q(c,t)e-ξαdξ)γ=(-E1(t))γ(m1bdξ)η(q(a,t)q(c,t)e-ξαdξ)γ=(-E1(t))γ(m01bdξ)η(αe-q(a,t)α-αe-q(c,t)α)γ.

It follows

(αe-q(a,t)α-αe-q(c,t)α)γ-m0dx(-E1(t))γ((-m0)1bdξ)η+,
then we can obtain
e-q(a,t)α-e-q(c,t)α+,ast+.

If b = 0, we can obtain

-m0dx=-limb0mdx=-m(q,t)qxdx=limb0[ac(-meqαqx)γ((-m)1bqx)η(e-qαqx)γdx]limb0(ac-meqαqxdx)γ(ac(-m)1bqxdx)η(ace-qαqxdx)γ=limb0(q(a,t)q(c,t)-meξαdξ)γ((-m)1bdξ)η(q(a,t)q(c,t)e-ξαdξ)γ=limb0(-E1(t))γ(m1bdξ)η(q(a,t)q(c,t)e-ξαdξ)γ=limb0(-E1(t))γ(m01bdξ)η(αe-q(a,t)α-αe-q(c,t)α)γ.

It follows

(αe-q(a,t)α-αe-q(c,t)α)γ-m0dx(-E1(t))γ(limb0(-m0)1bdξ)η+,
then we can obtain
e-q(a,t)α-e-q(c,t)α+,ast+.

6. CONCLUSION

We have considered the FOCH model αβ, α > 0, β > 0. When α, β is negative, one can get the same results by taking absolute value |α| and |β|. This model is highly related to the classical Camassa-Holm equation, the Degasperis-Procesi equation and the Holm-Staley b-family equation. We have studied some mathematical property, such as global existence, infinite propagation speed and long time behavior of the support of momentum density. Another highly related equation is (1.1) with α = β. Due to (1.1) with α = β doesn’t have the structure in Lemma 2.2 and (2.1), some results in this manuscript may can’t been realized for α = β.

CONFLICTS OF INTEREST

The authors declare they have no conflicts of interest.

ACKNOWLEDGEMENTS

This work has partially been supported by the National Natural Science Foundation of China (No. 12071439 and No. 11971475) and the 2019 - 2020 Hunan overseas distinguished professorship project (No. 2019014). The third author Jiang was partially supported by NSFC (grant No. 11101376, 11671309, 11671364) and ZJNSF (grant No. LY19A010016). The fourth author Qiao thanks the UT President’s Endowed Professorship (Project # 450000123) for its partial support.

REFERENCES

[14]A Degasperis, DD Holm, and ANI Khon, A new integrable equation with peakon solutions, Teoret. Mat. Fiz., Vol. 133, 2002, pp. 170-183.
[15]A Degasperis and M Procesi, Asymptotic integrability, A Degasperis and G Gaeta (editors), Symmetry and Perturbation Theory, World Scientific, Singapore, 1999, pp. 23-37.
[22]DD Holm and ANW Hone, Nonintegrability of a fifth-order equation with integrable two-body dynamics, Theoret. Math. Phys., Vol. 137, 2003, pp. 1459-1471.
[32]HP McKean, Breakdown of a shallow water equation, Asian J. Math., Vol. 2, 1998, pp. 867-874.
[33]OG Mustafa, A note on the Degasperis-Procesi equation, J. Nonlinear Math. Phys., Vol. 12, 2005, pp. 10-14.
[41]F Wang, F Li, and Z Qiao, On the Cauchy problem for a higher-order μ-Camassa–Holm equation, Discrete Contin. Dyn. Syst., Vol. 38, 2018, pp. 4163-4187.
[43]Y Zhou, On solutions to the Holm Staley b-family of equations, Nonlinearity, Vol. 23, 2010, pp. 369-38.
Journal
Journal of Nonlinear Mathematical Physics
Volume-Issue
28 - 3
Pages
321 - 336
Publication Date
2021/06/16
ISSN (Online)
1776-0852
ISSN (Print)
1402-9251
DOI
10.2991/jnmp.k.210519.001How to use a DOI?
Copyright
© 2021 The Authors. Published by Atlantis Press B.V.
Open Access
This is an open access article distributed under the CC BY-NC 4.0 license (http://creativecommons.org/licenses/by-nc/4.0/).

Cite this article

TY  - JOUR
AU  - Mingxuan Zhu
AU  - Lu Cao
AU  - Zaihong Jiang
AU  - Zhijun Qiao
PY  - 2021
DA  - 2021/06/16
TI  - Analytical Properties for the Fifth Order Camassa-Holm (FOCH) Model
JO  - Journal of Nonlinear Mathematical Physics
SP  - 321
EP  - 336
VL  - 28
IS  - 3
SN  - 1776-0852
UR  - https://doi.org/10.2991/jnmp.k.210519.001
DO  - 10.2991/jnmp.k.210519.001
ID  - Zhu2021
ER  -