355
Views
1
CrossRef citations to date
0
Altmetric
Research Article

A study on the stability of a modified Degasperis–Procesi equation

| (Reviewing Editor)
Article: 1251875 | Received 12 Jun 2016, Accepted 12 Oct 2016, Published online: 07 Nov 2016

Abstract

A modified Degasperis–Procesi equation is investigated. The local existence and uniqueness of the strong solution for the equation are established in the Sobolev space Hs(R) with s>32. The L1(R) local stability for the strong solution is obtained under certain assumptions on the initial data.

AMS Subject Classifications:

Public Interest Statement

In this paper, we have studied the generalized Degasperis–Procesi equation. Using the partial differential operator, the equivalent form Equation (6) has been derived. In the previous sections, we have used the Kato Theorem and double variables method to establish existence and stability of the solution for Equation (6). The approaches presented in this paper can be summarized to discuss other partial differential equations with initial value.

1. Introduction and main results

The Degasperis–Procesi (DP) equation of the form(1) ϕt-ϕtxx+4ϕϕx=3ϕxϕxx+ϕϕxxx,t>0,xR,(1)

which represents a model for shallow water dynamics, has been investigated by many scholars (see Coclite & Karlsen, Citation2006,Citation2007; Degasperis, Holm, & Hone, Citation2002; Escher, Liu, & Yin, Citation2006; Lai, Yan, Chen, & Wang, Citation2014; Lai & Wu, Citation2010; Lin & Liu, Citation2009; Lenells, Citation2005; Yin, Citation2003). Coclite and Karlsen (Citation2006) established existence and L1 stability results for entropy weak solutions of Equation (3) in the space L1BV and extended the results to a kind of generalized Degasperis–Procesi equations. Escher et al. (Citation2006) discussed several qualitative properties of the Degasperis–Procesi equation. The existence and uniqueness of global weak solutions for Equation (1) have been established, provided that the initial data satisfy appropriate conditions in Escher (Citation2006). Lenells (Citation2005) dealt with the travelling wave solutions of Equation (1) and classified all weak travelling wave solutions of the Degasperis–Procesi equation. Recently, Lai et al. (Citation2014) studied the generalized Degasperis–Procesi equation(2) ϕt-ϕtxx+kϕx+mϕϕx=3ϕxϕxx+ϕϕxxx,(2)

where k0 and m>0 are constants. Lai et al. (Citation2014) derived the L2(R) conservation law and established the L1(R) stability of local strong solutions to Equation (2) by assuming that its initial value belongs to the space Hs(R) with s>32. For other approaches to study the DP equation and related partial differential equations, the reader is referred to Kato (Citation1975), Kruzkov (Citation1970), Rodriguez-Blanco (Citation2001) and the references therein.

The objective of this work is to investigate the modified Degasperis–Procesi equation in the form(3) ϕt-ϕtxx+βf(ϕ)ϕx=f(ϕ)ϕx3+3f(ϕ)ϕxϕxx+f(ϕ)ϕxxx+αm[(m-1)ϕm-2ϕx2+ϕm-1ϕxx]-αϕm,(3)

where α,βR,m2, function f(ϕ) is a polynomial of order n(n2). Letting α=0,β=4,f(ϕ)=ϕ22, Equation (3) reduces to the Degasperis–Procesi Equation (1). Applying the operator 1-x2-1 to Equation (3), we obtain its equivalent form(4) ϕt+f(ϕ)ϕx+αϕm+(β-1)(1-x2)-1xf(ϕ)=0,(4)

where (1-x2)-1I(t,x)=12Re-|x-y|I(t,y)dy. Assuming that the initial value ϕ(0,x) of Equation (4) belongs to Hs(R)s>32, we will prove the existence and uniqueness of local solution for Equation (4) using the Kato Theorem (see Kato, Citation1975) in the space C[0,);Hs(R)C1[0,);Hs-1(R). Furthermore, we will use the approaches presented in Kruzkov (Citation1970) to establish the L1(R) local stability of the solution for the modified Degasperis–Procesi Equation (4). The results obtained in this paper extend parts of results presented in Lai et al. (Citation2014).

We let Hs(R) (where s is a real number) denote the Sobolev space with the norm defined byhHs=R1+|ξ|2sh^(t,ξ)212<,

where h^(t,ξ)=Re-ixξh(t,x)dx. For T>0 and s0, we let C[0,);Hs(R) denote the Fréchet space of all continuous Hs-valued functions on [0, T). Set Λ=(1-x2)12. For simplicity, we let c denote any positive constants.

We consider the Cauchy problem of Equation (3)(5) ϕt-ϕtxx+βf(ϕ)ϕx=f(ϕ)ϕx3+3f(ϕ)ϕxϕxx+f(ϕ)ϕxxx+αm[(m-1)ϕm-2ϕx2+ϕm-1ϕxx]-αϕm,ϕ(0,x)=ϕ0(x),(5)

which is equivalent to the problem(6) ϕt+f(ϕ)ϕx+αϕm+(β-1)Λ-2xf(ϕ)=0,ϕ(0,x)=ϕ0(x).(6)

Now we state the main results of our work.

Theorem 1.1

Let ϕ0(x)Hs(R) with s>32. There exists a T>0 depending on ϕ0Hs(R) such that problem (5) or (6) has a unique solution ϕ(t,x)C[0,T);Hs(R)C1[0,T);Hs-1(R).

Theorem 1.2

Assume that ϕ(t,x) and ψ(t,x) are two local strong solutions of problem (5) or (6) with initial data ϕ0(x),ψ0(x)L1(R)Hs(R)s>32, respectively. Let T>0 be the maximum existence time of ϕ(t,x) and ψ(t,x). For any t[0,T), it holds thatϕ(t,x)-ψ(t,x)L1ect-+|ϕ0(x)-ψ0(x)|dx,

where c is a positive constant depending on ϕ0L(R) and ψ0L(R).

This paper is organized as follows. Section 2 gives the proof of Theorem 1.1. The proof of Theorem 1.2 is given in Section 3.

2. Proof of Theorem 1.1

Firstly, we introduce the abstract quasi-linear evolution equation(7) dudt+A(u)u=W(u),t0,u(0)=u0.(7)

Let X and Y be Hilbert spaces where Y is continuously and densely embedded in X, and Q:YX be a topological isomorphism. We define L(YX) as the space of all bounded linear operators from Y to X. We denote L(XX) by L(X). Note that ρ1,ρ2,ρ3 and ρ4 in the following are constants and depend on max{yY,zY}.

(I)

A(y)L(Y,X) for yX with A(y)-A(z)wXρ1y-zXwY,y,z,wY,A(y)G(X,1,ξ)(ξ>0) and uniform on bounded sets in Y.

(II)

QA(y)Q-1=A(y)+B(y), in which B(y)L(X) is bounded and uniform on bounded sets in Y and B(y)-B(z)wXρ2y-zXwY,y,zY,wX.

(III)

W:YY extends to a map from X to X, is bounded on bounded sets in Y and satisfies W(y)-W(z)Yρ3y-zY,y,zY,W(y)-W(z)Xρ4y-zX,y,zY.

Kato Theorem (Citation1975) . Assume that (I), (II) and (III) hold. If u0Y, there is a maximal T>0 depending only on u0Y and a unique solution u to problem (7) such thatu=u(.,u0)C[0,T);YC1[0,T);X.

Moreover, the map u0u(.,u0) is a continuous map from Y to the space C[0,T);YC1[0,T);X.

For problem (6), we set A(ϕ)=f(ϕ)x,Y=Hs(R), X=Hs-1(R), W(ϕ)=-αϕm-(β-1)Λ-2xf(ϕ), and Q=Λ. Then, we will verify that A(ϕ) and W(ϕ) satisfy conditions (I)–(III). We cite several conclusions presented in Rodriguez–Blanco (Citation2001).

Lemma 2.1

The operator A(ϕ)=f(ϕ)x with ϕHs(R)s>32, belongs to G(Hs-1,1,ξ).

Lemma 2.2

For ϕ,z,wHs(R) with s>32, A(ϕ)L(Hs,Hs-1), it holds thatA(ϕ)-A(z)wHs-1ρ1ϕ-zHs-1wHs.

Lemma 2.3

For ϕ,zHs(R) and wHs-1(R)s>32, it holds that B(ϕ)=Λ,f(ϕ)xΛ-1L(Hs-1) andB(ϕ)-B(z)wHs-1ρ2ϕ-zHswHs-1.

Lemma 2.4

(Kato, Citation1975).   Let r and q be real numbers such that -r<qr. Then,ϕψHqleqcϕHrψHq,ifr>12,ϕψHr+q-12cϕHrψHq,ifr<12.

Lemma 2.5

Let ϕ,zHs(R) with s>32 and W(ϕ)=-αϕm-(β-1)Λ-2xf(ϕ). Then, W is bounded on bounded sets in Hs and satisfies(8) W(ϕ)-W(z)Hsρ3ϕ-zHs,(8) (9) W(ϕ)-W(z)Hs-1ρ4ϕ-zHs-1.(9)

Proof

For s>32, we have ϕLcϕHs and ϕHs-1cϕHs. Applying the algebra property of Hs(R) and Lemma 2.4, we getW(ϕ)-W(z)Hs|α|ϕm-zmHs+|β-1|Λ-2xf(ϕ)-f(z)Hscϕm-zmHs+cf(ϕ)-f(z)Hs-1cϕ-zHsϕm-1+ϕm-2z++zm-1Hs+cf(ϕ)-f(z)Hsρ3ϕ-zHs,

which completes the proof of (8). Similarly, we getW(ϕ)-W(z)Hs-1αϕm-αzmHs-1+|β-1|Λ-2xf(ϕ)-f(z)Hs-1cϕm-zmHs-1+cf(ϕ)-f(z)Hs-2cϕ-zHs-1ϕm-1+ϕm-2z++zm-1Hs-1+cf(ϕ)-f(z)Hs-1ρ4ϕ-zHs-1,

which completes the proof of (9).

Proof of Theroem 1.1 Using Lemmas 2.1–2.3 and Lemma 2.5, we know that the conditions (I), (II) and (III) hold. Applying the Kato Theorem, we find that problem (5) or (6) has a unique local solutionϕ=ϕ(t,x)C[0,T);Hs(R)C1[0,T);Hs-1(R),

where T>0 depends on ϕ0Hs.

Remark 2.6

Let T>0 be described in Theorem 1.1. Using the Sobolev embedding Theorem, we ensure the boundedness of solution ϕ(t,x) to problem (6) in the domain [0,T)×R. Namely, provided that ϕ0Hs(R) with s>32, we have ϕL(R)MT, where MT is a positive constant.

3. Proof of Theorem 1.2

Let P(t,x,ϕ)=αϕm+(β-1)Λ-2xf(ϕ) in the first equation of (6); we get,(10) ϕt+f(ϕ)ϕx+P(t,x,ϕ)=0,ϕ(0,x)=ϕ0(x).(10)

Assume that ϕ(t,x) and ψ(t,x) are solutions of problem (10) in the domain [0,T)×R with initial functions ϕ0(x) and ψ0(x)Hs(R)s>32. Now we give several lemmas.

Lemma 3.1

Let ϕ(t,x) be the solution of problem (10) and ϕ0(x)Hs(R) with s>32. Then, ϕL(R)MT andP(t,x,ϕ)L(R)cMTk,

where positive constant c depends on α,β,m,n, ϕ0L and k=max{m,n}.

Proof

We haveP(t,x,ϕ)=αϕm(t,x)+(β-1)xΛ-2f(ϕ)dy=αϕm(t,x)+β-12Re-|x-y|sign(y-x)f(ϕ)dy|α|ϕm(t,x)+|β-1|2Re-|x-y|f(ϕ)dy.

Applying Remark 2.6 and the integral Re-|x-y|dy=2, we complete the proof.

Lemma 3.2

Assume that ϕ(t,x) and ψ(t,x) are solutions of problem (10) in the domain [0,T)×R with initial functions ϕ0(x) and ψ0(x)Hs(R)(s>32), respectively. Then,(11) -+|P(t,x,ϕ)-P(t,x,ψ)|dxc-+|ϕ-ψ|dx,(11)

where c>0 depends on α,β,m,n,ϕ0L(R),ψ0L(R) and T.

Proof

Using the property of the operator Λ-2 and Remark 2.6, we get-+|P(t,x,ϕ)-P(t,x,ψ)|dx|α|-+|ϕm-ψm|dx+|β-1|2-+xΛ-2(f(ϕ)-f(ψ))dxc-+|ϕ-ψ|dx+c-+dx-+e-|x-y|sign(y-x)(f(ϕ)-f(ψ))dyc-+|ϕ-ψ|dx+c-+|f(ϕ)-f(ψ)|dx-+e-|x-y|dyc-+|ϕ-ψ|dx,

in which we apply the Tonelli Theorem to complete the proof.

We introduce a function δ(σ) which is infinitely differential on (-,+). Note that δ(σ)0, δ(σ)0 for |σ|1, -+δ(σ)dσ=1. Let δε(σ)=δ(ε-1σ)ε, where ε is an arbitrary positive constant. It is found that δε(σ)C0(-,+) and(12) δε(σ)0,δε(σ)=0for|σ|ε,|δε(σ)|cε,-+δε(σ)dσ=1.(12)

Let the function v(x) be defined and locally integrable on (-,+). Let vε(x) denote the approximation function of v(x) as(13) vε(x)=1ε-+δx-yεv(y)dy.(13)

We call x0 a Lebesgue point of the function v(x) iflimε01εx0-εx0+ε|v(x)-v(x0)|dx=0.

At any Lebesgue point x0, we getlimε0vε(x0)=v(x0).

Since the set of points which are not Lebesgue points of v(x) has measure zero, we have vε(x)v(x) as ε0 almost everywhere.

For any T1[0,T), we denote the band {(t,x):[0,T1]×R} by πT1. Let Kr={x:|x|r} andΠ=(t,x,τ,y):t-τ2ε,ρt+τ2T1-ρ,x-y2ε,x+y2r-ρ,

where r>0, ρ>0.

Lemma 3.3

(Kruzkov, Citation1970).    Let the function v(tx) be bounded and measurable in cylinder [0,T1]×Kr. If for any ρ(0,min[r,T1]) and any ε(0,ρ), the functionVε=1ε2Π|v(t,x)-v(τ,y)|dtdxdτdy

satisfieslimε0Vε=0.

Lemma 3.4

(Kruzkov, Citation1970).   If F(u)u is bounded, then the function H(u,v)=sign(u-v)(F(u)-F(v)) satisfies the Lipschitz condition in u and v.

We state the concept of a characteristic cone. Let T be described in Theorem 1.1 and ϕL(R)MT. For any T1[0,T) and R1>0, we defineN>max(t,x)[0,T1]×KR1f(ϕ).

Let Ω represent the cone (t,x):|x|R1-Nt,0tT0=min(T1,R1N-1) and Sτ designate the cross section of the cone Ω by the plane t=τ,τ[0,T0].

Lemma 3.5

Let ϕ(t,x) be the solution of problem (10) on πT1, g(t,x)C0(πT1); it holds that(14) πT1|ϕ-k|gt+sign(ϕ-k)[f(ϕ)-f(k)]gx-sign(ϕ-k)g(t,x)P(t,x,ϕ)dtdx=0,(14)

where k is an arbitrary constant.

Proof

Suppose that Φ(ϕ) is a twice differential function. Multiplying the first equation of problem (10) by Φ(ϕ)g(t,x) and integrating over πT1, we get(15) πT1Φ(ϕ)gϕt+Φ(ϕ)gf(ϕ)ϕx+Φ(ϕ)gP(t,x,ϕ)dtdx=0.(15)

Using the method of integration by parts, we get(16) πT1Φ(ϕ)gϕtdtdx=-πT1Φ(ϕ)gtdtdx.(16)

Notice thatkϕΦ(z)f(z)dz=Φ(ϕ)f(ϕ)ϕx.

Thus,-+kϕΦ(z)f(z)dzgxdx=--+Φ(ϕ)f(ϕ)ϕxgdx.

Then, we have(17) πT1Φ(ϕ)f(ϕ)ϕxgdtdx=-πT1kϕΦ(z)f(z)dzgxdtdx=-πT1Φ(ϕ)[f(ϕ)-f(k)]-kϕΦ(z)f(z)-f(k)dzgxdtdx.(17)

Substitute Equations (16) and (17) into Equation (15). Let Φε(ϕ) be an approximation of the function Φ(ϕ)=|ϕ-k|. When ε0, Φε(ϕ)Φ(ϕ), we obtain Equation (14).

We will give the proof of Theorem 1.2. Set function g(t,x)C0(πT1),g(t,x)0 outside the cylinder ϖ={(t,x)}=[ρ,T1-2ρ]×Kr-2ρ, where Kr-2ρ={|x|:|x|r-2ρ},r>0,0<2ρ<min(T1,r).

Proof of Theorem 1.2.

We define(18) F(t,x,τ,y)=gt+τ2,x+y2δεt-τ2,x-y2=g()λε(),(18)

in which ()=t+τ2,x+y2,=t-τ2,x-y2. Thus, we obtainFt+Fτ=gt()λε(),Fx+Fy=gx()λε().

Using Lemma 3.5 and setting k=ψ(τ,y),g(t,x)=F(t,x,τ,y), we get(19) πT1×πT1|ϕ(t,x)-ψ(τ,y)|Ft+sign(ϕ(t,x)-ψ(τ,y))×f(ϕ(t,x))-f(ψ(τ,y))Fx-sign(ϕ(t,x)-ψ(τ,y))×P(t,x,ϕ(t,x))Fdtdxdτdy=0.(19)

Similarly, we have(20) πT1×πT1|ψ(τ,y)-ϕ(t,x)|Fτ+sign(ψ(τ,y)-ϕ(t,x))×f(ψ(τ,y))-f(ϕ(t,x))Fy-sign(ψ(τ,y)-ϕ(t,x))×P(τ,y,ψ(τ,y))Fdtdxdτdy=0.(20)

Adding (19) and (20), we obtain(21) 0πT1×πT1ϕ(t,x)-ψ(τ,y)gtλε+sign(ϕ(t,x)-ψ(τ,y))×f(ϕ(t,x))-f(ψ(τ,y))gxλεdtdxdτdy+πT1×πT1sign(ϕ(t,x)-ψ(τ,y))×P(t,x,ϕ(t,x))-P(τ,y,ψ(τ,y))gλεdtdxdτdy.(21)

We note that the first two terms of the integrand of (21) have the form(22) Kε=Kt,x,τ,y,ϕ(t,x),ψ(τ,y)λε(),(22)

where K satisfies the Lipschitz condition in all its variables. Then,πT1×πT1Kεdtdxdτdy=πT1×πT1Kt,x,τ,y,ϕ(t,x),ψ(τ,y)λεdtdxdτdy=πT1×πT1Kt,x,τ,y,ϕ(t,x),ψ(τ,y)-Kt,x,t,x,ϕ(t,x),ψ(t,x)λεdtdxdτdy+πT1×πT1Kt,x,t,x,ϕ(t,x),ψ(t,x)λεdtdxdτdy=J11(ε)+J12

Note that Kε=0 outside the region Π. Applying the estimate λε()|cε2 and Lemma 3.4, we get J11(ε)cε+1ε2Πψ(t,x)-ψ(τ,y)dtdxdτdy,

where c is a positive constant independent of ε. Using Lemma 3.3, we know J11(ε)0 as ε0.

For the term J12, we substitute t=α1,t-τ2=β1,x=η,x-y2=γ. Combining with the identity-εε-+λε(β1,γ)dβ1dγ=1,

we obtainJ12=22πT1×πT1Kα1,η,α1,η,ϕ(α1,η),ψ(α1,η)-εε-+λε(β1,γ)dβ1dγdα1dη=4πT1Kt,x,t,x,ϕ(t,x),ψ(t,x)dtdx.

Thus, we have(23) limε0πT1×πT1Kε=4πT1Kt,x,t,x,ϕ(t,x),ψ(t,x)dtdx.(23)

Similarly, the integrand of the third term in (21) can be represented as(24) K¯ε=signϕ(t,x)-ψ(τ,y)P(t,x,ϕ(t,x))-P(τ,y,ψ(τ,y))gλε=K¯t,x,τ,y,ϕ(t,x),ψ(τ,y)λε().(24)

Then,πT1×πT1K¯εdtdxdτdy=πT1×πT1K¯t,x,τ,y,ϕ(t,x),ψ(τ,y)-K¯t,x,t,x,ϕ(t,x),ψ(t,x)λεdtdxdτdy+πT1×πT1K¯t,x,t,x,ϕ(t,x),ψ(t,x)λεdtdxdτdy=J21(ε)+J22.

Using Lemma 3.4, we have J21(ε)cε+1ε2Πψ(t,x)-ψ(τ,y)dtdxdτdy.

Using Lemma 3.3, it yields J21(ε)0 as ε0. Repeating the steps as before, we have(25) J22=4πT1K¯t,x,t,x,ϕ(t,x),ψ(t,x)dtdx.(25)

From (21) to (25), we get(26) πT1ϕ(t,x)-ψ(t,x)gt+sign(ϕ(t,x)-ψ(t,x))f(ϕ)-f(ψ)gxdtdx+πT1sign(ϕ(t,x)-ψ(t,x))P(t,x,ϕ)-P(t,x,ψ)gdtdx0.(26)

We set(27) h(t)=-+ϕ(t,x)-ψ(t,x)dx(27)

and(28) με(σ)=-σδε(σ)dσ.(28)

Take two numbers ρ, τ(0,T0) and ρ<τ. In (26), we let(29) g(t,x)=με(t-ρ)-με(t-τ)χ(t,x),ε<min(ρ,T0-τ),(29)

in which(30) χ(t,x)=χθ(t,x)=1-μθ(|x|+Nt-R1+θ),(30)

where θ is a small positive constant and χ(t,x)=0 outside the cone Ω. When θ0,R1+, we observe that χθ1. By the definition of the number N, we have0=χt+N|χx|χt+Nχx,(t,x)Ω.

Applying (26)–(30), we get(31) πT0ϕ(t,x)-ψ(t,x)δε(t-ρ)-δε(t-τ)χθ(t,x)dtdx+0T0dt-+P(t,x,ϕ)-P(t,x,ψ)×με(t-ρ)-με(t-τ)χθ(t,x)dx0.(31)

In (31), sending θ0,R1+ and using Lemma 3.2, we obtain(32) 0T0δε(t-ρ)-δε(t-τ)h(t)dt+c0T0με(t-ρ)-με(t-τ)h(t)dt0,(32)

where c is independent of ε.

Applying the properties of the function δε for εmin(ρ,T0-ρ), we get0T0δε(t-ρ)h(t)-h(ρ)dt=0T0δε(t-ρ)[h(t)-h(ρ)]dtcερ-ερ+εh(t)-h(ρ)dt.

Then,(33) 0T0δε(t-ρ)h(t)dth(ρ)asε0.(33)

LetG(ρ)=0T0με(t-ρ)h(t)dt=0T0dt-t-ρδε(σ)h(t)dσ.

We observe thatG(ρ)=-0T0δε(t-ρ)h(t)dt.

Let ε0; it derives thatG(ρ)-h(ρ)

andG(ρ)G(0)-0ρh(t)dt,G(τ)G(0)-0τh(t)dt.

Therefore, we have(34) G(ρ)-G(τ)ρτh(t)dtasε0.(34)

From (32)–(34), we obtain inequalityh(ρ)+cρτh(t)dth(τ).

Let ρ0,τt; we get-+|ϕ(0,x)-ψ(0,x)|dx+c0tdt-+|ϕ(t,x)-ψ(t,x)|dx-+|ϕ(t,x)-ψ(t,x)|dx.

Applying the Gronwall inequality, we complete the proof of Theorem 1.2.

Additional information

Funding

This work was supported by National Natural Science Foundation of China (CN) [grant number 11471263].

Notes on contributors

Jing Chen

Jing Chen was born in Sichuan, China, on July 1978. She received her MS degree in mathematics from Sichuan Normal University, in 2007.She has been working on the research of the existence and stability of solutions of partial differential equations and has published many articles in international journals. At present, she is teaching at Southwest University of Science and Technology.

References

  • Coclite, G. M., & Karlsen, K. H. (2006). On the well-posedness of the Degasperis-Procesi equation. Journal of Functional Analysis, 233, 60–91.
  • Coclite, G. M., & Karlsen, K. H. (2007). On the uniqueness of discontinuous solutions to the Degasperis-Procesi equation. Journal of Differential Equations, 234, 142–160.
  • Degasperis, A., Holm, D., & Hone, A. (2002). A new integral equation with peakon solutions. Theoretical and Mathematics Physis, 133, 1461–1472.
  • Escher, J., Liu, Y., & Yin, Z. Y. (2006). Global weak solutions and blow-up structre for the Degasperis-Procesi equation. Journal of Functional Analysis, 241, 457–485.
  • Kato, T. (1975). Quasi-linear equations of evalution with applications to partial differential equations. Spectral Theory and Differential Equations. In Lecture Notes in Math (Vol. 448, pp. 25–70). Berlin: Springer-Verlag.
  • Kruzkov, S. N. (1970). First order quasi-linear equations in several independent variables. Math USSR-Sb, 10, 217–243.
  • Lai, S. Y., & Wu, Y. H. (2010). Global solutions and blow-up phenomena to a shallow water equation. Journal of Differential Equations, 249, 693–706.
  • Lai, S. Y., Yan, H. B., Chen, H. J., & Wang, Y. (2014). The stability of local strong solutions for a shallow water equation. Journal of Inequalities and Applications, 410, 410–423.
  • Lenells, J. (2005). Travelling wave solutions of the Degasperis-Procesi equation. Journal of Mathematics Analysis and Applications, 306, 72–82.
  • Lin, Z., & Liu, Y. (2009). Stability of peakons for the Degasperis-Procesi equation. Communications on Pure and Applied Mathematics, 62, 125–146.
  • Rodriguez-Blanco, G. (2001). On the Cauthy problem for the Camassa-Holm equation. Nonlinear Analysis, 46, 309–327.
  • Yin, Z. Y. (2003). On the Cauthy problem for an integrable equation with peakon solutions. Illinois Journal of Mathematics, 47, 649–666.