With the help of heat equation, we first construct an example of a graphical solution to the curve shortening flow. This solution y(x,t) has the interesting property that it converges to a log-periodic function of the form
Asin(logt)+Bcos(logt)
as t→∞, where A, B are constants. Moreover, for any two numbers α<β, we are also able to construct a solution satisfying the oscillation limits
lim inft→∞y(x,t)=α, lim supt→∞y(x,t)=β, x∈K
on any compact subset K⊂R.
Citation: Dong-Ho Tsai, Xiao-Liu Wang. On an asymptotically log-periodic solution to the graphical curve shortening flow equation[J]. Mathematics in Engineering, 2022, 4(3): 1-14. doi: 10.3934/mine.2022019
[1] | Adel M. Al-Mahdi, Mohammad M. Al-Gharabli, Mohammad Kafini . Asymptotic behavior of the wave equation solution with nonlinear boundary damping and source term of variable exponent-type. AIMS Mathematics, 2024, 9(11): 30638-30654. doi: 10.3934/math.20241479 |
[2] | Khaled zennir, Djamel Ouchenane, Abdelbaki Choucha, Mohamad Biomy . Well-posedness and stability for Bresse-Timoshenko type systems with thermodiffusion effects and nonlinear damping. AIMS Mathematics, 2021, 6(3): 2704-2721. doi: 10.3934/math.2021164 |
[3] | Aihui Sun, Hui Xu . Decay estimate and blow-up for fractional Kirchhoff wave equations involving a logarithmic source. AIMS Mathematics, 2025, 10(6): 14032-14054. doi: 10.3934/math.2025631 |
[4] | Mohammad Kafini, Mohammad M. Al-Gharabli, Adel M. Al-Mahdi . Existence and stability results of nonlinear swelling equations with logarithmic source terms. AIMS Mathematics, 2024, 9(5): 12825-12851. doi: 10.3934/math.2024627 |
[5] | Mohamed Biomy . Decay rate for systems of m-nonlinear wave equations with new viscoelastic structures. AIMS Mathematics, 2021, 6(6): 5502-5517. doi: 10.3934/math.2021326 |
[6] | Keltoum Bouhali, Sulima Ahmed Zubair, Wiem Abedelmonem Salah Ben Khalifa, Najla ELzein AbuKaswi Osman, Khaled Zennir . A new strict decay rate for systems of longitudinal m-nonlinear viscoelastic wave equations. AIMS Mathematics, 2023, 8(1): 962-976. doi: 10.3934/math.2023046 |
[7] | Hicham Saber, Fares Yazid, Fatima Siham Djeradi, Mohamed Bouye, Khaled Zennir . Decay for thermoelastic laminated beam with nonlinear delay and nonlinear structural damping. AIMS Mathematics, 2024, 9(3): 6916-6932. doi: 10.3934/math.2024337 |
[8] | Abdelbaki Choucha, Salah Boulaaras, Asma Alharbi . Global existence and asymptotic behavior for a viscoelastic Kirchhoff equation with a logarithmic nonlinearity, distributed delay and Balakrishnan-Taylor damping terms. AIMS Mathematics, 2022, 7(3): 4517-4539. doi: 10.3934/math.2022252 |
[9] | Adel M. Al-Mahdi, Mohammad M. Al-Gharabli, Maher Nour, Mostafa Zahri . Stabilization of a viscoelastic wave equation with boundary damping and variable exponents: Theoretical and numerical study. AIMS Mathematics, 2022, 7(8): 15370-15401. doi: 10.3934/math.2022842 |
[10] | Qian Li . General and optimal decay rates for a viscoelastic wave equation with strong damping. AIMS Mathematics, 2022, 7(10): 18282-18296. doi: 10.3934/math.20221006 |
With the help of heat equation, we first construct an example of a graphical solution to the curve shortening flow. This solution y(x,t) has the interesting property that it converges to a log-periodic function of the form
Asin(logt)+Bcos(logt)
as t→∞, where A, B are constants. Moreover, for any two numbers α<β, we are also able to construct a solution satisfying the oscillation limits
lim inft→∞y(x,t)=α, lim supt→∞y(x,t)=β, x∈K
on any compact subset K⊂R.
We consider, for x∈Rn, t>0, the following system
{utt−θΔ(u+∫t0ϖ1(t−s)u(s)ds)+αut=h1(u,v,w)vtt−θΔ(v+∫t0ϖ2(t−s)v(s)ds)+αvt=h2(u,v,w)wtt−θΔ(w+∫t0ϖ3(t−s)w(s)ds)+αwt=h3(u,v,w)u(x,0)=u0(x),v(x,0)=v0(x),w(x,0)=w0(x)ut(x,0)=u1(x),vt(x,0)=v1(x),wt(x,0)=w1(x), | (1.1) |
where n≥3,α>0, the functions hi(.,.,.)∈(R3,R),i=1,2,3 are given by
h1(ξ1,ξ2,ξ3)=(q+1)[d|ξ1+ξ2+ξ3|(q−1)(ξ1+ξ2+ξ3)+e|ξ1|(q−3)/2ξ1|ξ2|(q+1)/2]; |
h2(ξ1,ξ2,ξ3)=(q+1)[d|ξ1+ξ2+ξ3|(q−1)(ξ1+ξ2+ξ3)+e|ξ2|(q−3)/2ξ2|ξ3|(q+1)/2]; |
h3(ξ1,ξ2,ξ3)=(q+1)[d|ξ1+ξ2+ξ3|(q−1)(ξ1+ξ2+ξ3)+e|ξ3|(q−3)/2ξ3|ξ1|(q+1)/2], |
with d,e>0,q>3. The function 1θ(x)∼ϑ(x)>0 for all x∈Rn is a density such that
ϑ∈Lτ(Rn)withτ=2n2n−rn+2rfor2≤r≤2nn−2. | (1.2) |
As in [15], it is not hard to see that there exists a function G∈C1(R3,R) such that
uh1(u,v,w)+vh2(u,v,w)+wh3(u,v,w)=(q+1)G(u,v,w), ∀(u,v,w)∈R3. | (1.3) |
satisfies
(q+1)G(u,v,w)=|u+v+w|q+1+2|uv|(q+1)/2+2|vw|(q+1)/2+2|wu|(q+1)/2. | (1.4) |
We define the function spaces H as the closure of C∞0(Rn), as in [18], we have
H={v∈L2nn−2(Rn)∣∇v∈(L2(Rn))n}, |
with respect to the norm ‖v‖H=(v,v)1/2H for the inner product
(v,w)H=∫Rn∇v⋅∇wdx, |
and L2ϑ(Rn) as that to the norm ‖v‖L2ϑ=(v,v)1/2L2ϑ for
(v,w)L2ϑ=∫Rnϑvwdx. |
For general r∈[1,+∞)
‖v‖Lrϑ=(∫Rnϑ|v|rdx)1r. |
is the norm of the weighted space Lrϑ(Rn).
The main aim of this work is to consider an important problem from the point of view of application in sciences and engineering (materials which is something between that of elastic solids and Newtonian fluids), namely, a system of three wave equations having a damping effects in an unbounded domain with strong external forces including damping terms of memory type with past history. Using the Faedo-Galerkin [16] method and some energy estimates, we proved the existence of global solution in Rn owing to the weighted function. By imposing a new appropriate condition, with the help of some special estimates and generalized Poincaré's inequality, we obtained an unusual decay rate for the energy function. For more detail regarding the single equation, we review the following references [7,8]. The paper [7] is one of the pioneer in literature for the single equation, which is the source of inspiration of several researches, while the work [8] is a recent generalization of [7] by introducing less dissipative effects.
To enrich our topic, it is necessary to review previous works regarding the nonlinear coupled system of wave equations, from a qualitative and quantitative study. Let us begin with the single wave equation treated in [13], where the aim goal was mainely on the system
{utt+μut−Δu−ωΔut=uln|u|, (x,t)∈Ω×(0,∞)u(x,t)=0,x∈∂Ω,t≥0u(x,0)=u0(x),ut(x,0)=u1(x),x∈Ω, | (1.5) |
where Ω is a bounded domain of Rn, n≥1 with a smooth boundary ∂Ω. The author firstly constructed a local existence of weak solution by using contraction mapping principle and of course showed the global existence, decay rate and infinite time blow up of the solution with condition on initial energy.
Next, a nonexistence of global solutions for system of three semi-linear hyperbolic equations was introduced in [3]. A coupled system for semi-linear hyperbolic equations was investigated by many authors and a different results were obtained with the nonlinearities in the form f1=|u|q−1|v|q+1u,f2=|v|q−1|u|q+1v. (Please, see [2,5,9,14,24,29]).
In the case of non-bounded domain Rn, we mention the paper recently published by T. Miyasita and Kh. Zennir in [16], where the considered equation as follows
utt+aut−ϕ(x)Δ(u+ωut−∫t0g(t−s)u(s)ds)=u|u|q−1, | (1.6) |
with initial data
{u(x,0)=u0(x),ut(x,0)=u1(x). | (1.7) |
The authors showed the existence of unique local solution and they continued to extend it to be global in time. The rate of the decay for solution was the main result by considering the relaxation function is strictly convex, for more results related to decay rate of solution of this type of problems, please see [6,17,25,26,30,31].
Regarding the study of the coupled system of two nonlinear wave equations, it is worth recalling some of the work recently published. Baowei Feng et al. considered in [10], a coupled system for viscoelastic wave equations with nonlinear sources in bounded domain ((x,t)∈Ω×(0,∞)) with smooth boundary as follows
{utt−Δu+∫t0g(t−s)Δu(s)ds+ut=f1(u,v)vtt−Δv+∫t0h(t−s)Δv(s)ds+vt=f2(u,v). | (1.8) |
Here, the authors concerned with a system in Rn(n=1,2,3). Under appropriate hypotheses, they established a general decay result by multiplication techniques to extends some existing results for a single equation to the case of a coupled system.
It is worth noting here that there are several studies in this field and we particularly refer to the generalization that Shun et al. made in studying a complicate non-linear case with degenerate damping term in [22]. The IBVP for a system of nonlinear viscoelastic wave equations in a bounded domain was considered in the problem
{utt−Δu+∫t0g(t−s)Δu(s)ds+(|u|k+|v|q)|ut|m−1ut=f1(u,v) vtt−Δv+∫t0h(t−s)Δv(s)ds+(|v|κ+|u|ρ)|vt|r−1vt=f2(u,v)u(x,t)=v(x,t)=0,x∈∂Ω,t>0u(x,0)=u0(x),v(x,0)=v0(x)ut(x,0)=u1(x),vt(x,0)=v1(x), | (1.9) |
where Ω is a bounded domain with a smooth boundary. Given certain conditions on the kernel functions, degenerate damping and nonlinear source terms, they got a decay rate of the energy function for some initial data.
The lack of existence (Blow up) is considered one of the most important qualitative studies that must be spoken of, given its importance in terms of application in various applied sciences. Concerning the nonexistence of solution for a more degenerate case for coupled system of wave equations with different damping, we mention the papers [19,20,21,23,27].
In m-equations, paper in [1] considered a system
uitt+γuit−Δui+ui=m∑i,j=1,i≠j|uj|pj|ui|piui, i=1,2,…,m, | (1.10) |
where the absence of global solutions with positive initial energy was investigated.
We introduce a very useful Sobolev embedding and generalized Poincaré inequalities.
Lemma 1.1. [16] Let ϑ satisfy (1.2). For positive constants Cτ>0 and CP>0 depending only on ϑ and n, we have
‖v‖2nn−2≤Cτ‖v‖H, |
and
‖v‖L2ϑ≤CP‖v‖H, |
for v∈H.
Lemma 1.2. [12] Let ϑ satisfy (1.2), then the estimates
‖v‖Lrϑ≤Cr‖v‖H, |
and
Cr=Cτ‖ϑ‖1rτ, |
hold for v∈H. Here τ=2n/(2n−rn+2r) for 1≤r≤2n/(n−2).
We assume that the kernel functions ϖ1,ϖ2,ϖ3∈C1(R+,R+) satisfy
{1−¯ϖ1=l>0for¯ϖ1=∫+∞0ϖ1(s)ds, ϖ′1(t)≤0,1−¯ϖ2=m>0for¯ϖ2=∫+∞0ϖ2(s)ds, ϖ′2(t)≤0,1−¯ϖ3=ν>0for¯ϖ3=∫+∞0ϖ3(s)ds, ϖ′3(t)≤0, | (1.11) |
we mean by R+ the set {τ∣τ≥0}. Noting by
ϖ(t)=maxt≥0{ϖ1(t),ϖ2(t),ϖ3(t)}, | (1.12) |
and
ϖ0(t)=mint≥0{∫t0ϖ1(s)ds,∫t0ϖ2(s)ds,∫t0ϖ3(s)ds}. | (1.13) |
We assume that there is a function χ∈C1(R+,R+) such that
ϖ′i(t)+χ(ϖi(t))≤0,χ(0)=0,χ′(0)>0, i=1,2,3, | (1.14) |
for any ξ≥0.
Hölder and Young's inequalities give
‖uv‖(q+1)/2L(q+1)/2ϑ≤(‖u‖2L(q+1)ϑ+‖v‖2L(q+1)ϑ)(q+1)/2≤(l‖u‖2H+m‖v‖2H)(q+1)/2, | (1.15) |
and
‖vw‖(q+1)/2L(q+1)/2ϑ≤(m‖v‖2H+ν‖w‖2H)(q+1)/2, | (1.16) |
and
‖wu‖(q+1)/2L(q+1)/2ϑ≤(ν‖w‖2H+l‖u‖2H)(q+1)/2. | (1.17) |
Thanks to Minkowski's inequality to give
‖u+v+w‖(q+1)L(q+1)ϑ≤c(‖u‖2L(q+1)ϑ+‖v‖2L(q+1)ϑ+‖w‖2L(q+1)ϑ)(q+1)/2≤c(‖u‖2H+‖v‖2H+‖w‖2H)(q+1)/2. |
Then there exist η>0 such that
‖u+v+w‖(q+1)L(q+1)ϑ+2‖uv‖(q+1)/2L(q+1)/2ϑ+2‖vw‖(q+1)/2L(q+1)/2ϑ+2‖wu‖(q+1)/2L(q+1)/2ϑ≤η(l‖u‖2H+m‖v‖2H+ν‖w‖2H)(q+1)/2. | (1.18) |
We need to define positive constants λ0 and E0 by
λ0≡η−1/(q−1)andE0=(12−1q+1)η−2/(q−1). | (1.19) |
The mainely aim of the present paper is to obtain a novel decay rate of solution from the convexity property of the function χ given in Theorem 3.1.
We denote as in [18,28] an eigenpair {(λi,ei)}i∈N⊂R×H of
−θ(x)Δei=λieix∈Rn, |
for any i∈N, (θ(x))−1≡ϑ(x). Then
0<λ1≤λ2≤⋯≤λi≤⋯↑+∞, |
holds and {ei} is a complete orthonormal system in H.
Definition 1.3. The triplet functions (u,v,w) is said a weak solution to (1.1) on [0,T] if satisfies for x∈Rn,
{∫Rnϑ(x)(utt+αut)φdx+∫Rn∇u∇φdx−∫t0ϖ1(t−s)∇u(s)ds∇φdx=∫Rnϑ(x)h1(u,v,w)φdx,∫Rnϑ(x)(vtt+αvt)ψdx+∫Rn∇v∇ψdx−∫t0ϖ2(t−s)∇v(s)ds∇ψdx=∫Rnϑ(x)h2(u,v,w)ψdx,∫Rnϑ(x)(wtt+αwt)Ψdx+∫Rn∇v∇Ψdx−∫t0ϖ3(t−s)∇w(s)ds∇Ψdx=∫Rnϑ(x)h3(u,v,w)Ψdx, | (1.20) |
for all test functions φ,ψ,Ψ∈H for almost all t∈[0,T].
The next Theorem is concerned on the local solution (in time [0,T]).
Theorem 2.1. (Local existence) Assume that
1<q≤n+2n−2andthatn≥3. | (2.1) |
Let (u0,v0,w0)∈H3 and (u1,v1,w3)∈L2ϑ(Rn)×L2ϑ(Rn)×L2ϑ(Rn). Under the assumptions (1.2)–(1.17) and (1.11)–(1.14). Then (1.1) admits a unique local solution (u,v,w) such that
(u,v,w)∈X3T, XT≡C([0,T];H)∩C1([0,T];L2ϑ(Rn)), |
for sufficiently small T>0.
We prove the existence of global solution in time. Let us introduce the potential energy J:H3→R defined by
J(u,v,w)=(1−∫t0ϖ1(s)ds)‖u‖2H+(ϖ1∘u)+(1−∫t0ϖ2(s)ds)‖v‖2H+(ϖ2∘v)+(1−∫t0ϖ3(s)ds)‖w‖2H+(ϖ3∘w), | (2.2) |
where
(ϖj∘w)(t)=∫t0ϖj(t−s)‖w(t)−w(s)‖2Hds, |
for any w∈L2(Rn),j=1,2,3. The modified energy is defined by
E(t)=12(‖ut‖2L2ϑ+‖vt‖2L2ϑ+‖wt‖2L2ϑ)+12J(u,v,w)−∫Rnϑ(x)G(u,v,w)dx, | (2.3) |
Theorem 2.2. (Global existence) Let (1.2)–(1.17) and (1.11)–(1.14) hold. Under (2.1) and for sufficiently small (u0,u1),(v0,v1),(w0,w1)∈H×L2ϑ(Rn), problem (1.1) admits a unique global solution (u,v,w) such that
(u,v,w)∈X3, X≡C([0,+∞);H)∩C1([0,+∞);L2ϑ(Rn)). | (2.4) |
The next, Lemma will play an important role in the sequel.
Lemma 2.3. For (u,v,w)∈X3T, the functional E(t) associated with problem (1.1) is a decreasing energy.
Proof. For 0≤t1<t2≤T, we have
E(t2)−E(t1)=∫t2t1ddtE(t)dt=−12∫t2t1(ϖ1(t)‖u‖2H−(ϖ′1∘u))dt−12∫t2t1(ϖ2(t)‖v‖2H−(ϖ′2∘v))dt−12∫t2t1(ϖ3(t)‖w‖2H−(ϖ′3∘w))dt−α(‖ut‖2L2ϑ+‖vt‖2L2ϑ+‖wt‖2L2ϑ)≤0, |
owing to (1.11)–(1.14).
We sketch here the outline of the proof for local solution by a standard procedure(See [4,11,31]).
Proof. (Of Theorem 2.1.) Let (u0,u1),(v0,v1),(w0,w1)∈H×L2ϑ(Rn). For any (u,v,w)∈X3T, we can obtain weak solution of the related system
{ϑ(x)(ztt+αzt)−Δz=−∫t0ϖ1(t−s)Δu(s)ds+ϑ(x)h1(u,v,w)ϑ(x)(ytt+αyt)−Δy=−∫t0ϖ2(t−s)Δv(s)ds+ϑ(x)h2(u,v,w)ϑ(x)(ζtt+αζt)−Δζ=−∫t0ϖ3(t−s)Δw(s)ds+ϑ(x)h3(u,v,w)z(x,0)=u0(x),y(x,0)=v0(x),ζ(x,0)=w0(x)zt(x,0)=u1(x),yt(x,0)=v1(x),ζt(x,0)=w1(x). | (2.5) |
We reduces problem (2.5) to Cauchy problem for system of ODE by using the Faedo-Galerkin approximation. We then find a solution map ⊤:(u,v,w)↦(z,y,ζ) from X3T to X3T. We are now ready show that ⊤ is a contraction mapping in an appropriate subset of X3T for a small T>0. Hence ⊤ has a fixed point ⊤(u,v,w)=(u,v,w), which gives a unique solution in X3T.
We will show the global solution. By using conditions on functions ϖ1,ϖ2,ϖ3, we have
E(t)≥12J(u,v,w)−∫Rnϑ(x)G(u,v,w)dx≥12J(u,v,w)−1q+1‖u+v+w‖(q+1)L(q+1)ϑ−2q+1(‖uv‖(q+1)/2L(q+1)/2ϑ+‖vw‖(q+1)/2L(q+1)/2ϑ+‖wu‖(q+1)/2L(q+1)/2ϑ)≥12J(u,v,w)−ηq+1[l‖u‖2H+m‖v‖2H+ν‖w‖2H](q+1)/2≥12J(u,v,w)−ηq+1(J(u,v,w))(q+1)/2=G(ς), | (2.6) |
here ς2=J(u,v,w), for t∈[0,T), where
G(ξ)=12ξ2−ηq+1ξ(q+1). |
Noting that E0=G(λ0), given in (1.19). Then
{G(ξ)≥0in ξ∈[0,λ0]G(ξ)<0in ξ>λ0. | (2.7) |
Moreover, limξ→+∞G(ξ)→−∞. Then, we have the following lemma
Lemma 2.4. Let 0≤E(0)<E0.
(i) If ‖u0‖2H+‖v0‖2H+‖w0‖2H<λ20, then local solution of (1.1) satisfies
J(u,v,w)<λ20, ∀t∈[0,T). |
(ii) If ‖u0‖2H+‖v0‖2H+‖w0‖2H>λ20, then local solution of (1.1) satisfies
‖u‖2H+‖v‖2H+‖w‖2H>λ21, ∀t∈[0,T),λ1>λ0. |
Proof. Since 0≤E(0)<E0=G(λ0), there exist ξ1 and ξ2 such that G(ξ1)=G(ξ2)=E(0) with 0<ξ1<λ0<ξ2.
The case (i). By (2.6), we have
G(J(u0,v0,w0))≤E(0)=G(ξ1), |
which implies that J(u0,v0,w0)≤ξ21. Then we claim that J(u,v,w)≤ξ21, ∀t∈[0,T). Moreover, there exists t0∈(0,T) such that
ξ21<J(u(t0),v(t0),w(t0))<ξ22. |
Then
G(J(u(t0),v(t0),w(t0)))>E(0)≥E(t0), |
by Lemma 2.3, which contradicts (2.6). Hence we have
J(u,v,w)≤ξ21<λ20, ∀t∈[0,T). |
The case (ii). We can now show that
‖u0‖2H+‖v0‖2H+‖w0‖2H≥ξ22, |
and
‖u‖2H+‖v‖2H+‖w‖2H≥ξ22>λ20, |
in the same way as (i).
Proof. (Of Theorem 2.2.) Let (u0,u1),(v0,v1),(w0,w1)∈H×L2ϑ(Rn) satisfy both 0≤E(0)<E0 and
‖u0‖2H+‖v0‖2H+‖w0‖2H<λ20. |
By Lemma 2.3 and Lemma 2.4, we have
12(‖ut‖2L2ϑ+‖vt‖2L2ϑ+‖wt‖2L2ϑ)+l‖u‖2H+m‖v‖2H+ν‖w‖2H≤12(‖ut‖2L2ϑ+‖vt‖2L2ϑ+‖wt‖2L2ϑ)+(1−∫t0ϖ1(s)ds)‖u‖2H+(ϖ1∘u)+(1−∫t0ϖ2(s)ds)‖u‖2H+(ϖ2∘v)+(1−∫t0ϖ3(s)ds)‖w‖2H+(ϖ3∘w)≤2E(t)+2ηq+1[l‖u‖2H+m‖u‖2H+ν‖w‖2H](q+1)/2≤2E(0)+2ηq+1(J(u,v,w))(q+1)/2≤2E0+2ηq+1λq+10=η−2/(q−1). | (2.8) |
This completes the proof.
Let
Λ(u,v,w)=12(1−∫t0ϖ1(s)ds)‖u‖2H+12(ϖ1∘u)+12(1−∫t0ϖ2(s)ds)‖v‖2H+12(ϖ2∘v)+12(1−∫t0ϖ3(s)ds)‖w‖2H+12(ϖ3∘w)−∫Rnϑ(x)G(u,v,w)dx, | (2.9) |
Π(u,v,w)=(1−∫t0ϖ1(s)ds)‖u‖2H+(ϖ1∘u)+(1−∫t0ϖ2(s)ds)‖v‖2H+(ϖ2∘v)+(1−∫t0ϖ3(s)ds)‖w‖2H+(ϖ3∘w)−(q+1)∫Rnϑ(x)G(u,v,w)dx. | (2.10) |
Lemma 2.5. Let (u,v,w) be the solution of problem (1.1). If
‖u0‖2H+‖v0‖2H+‖w0‖2H−(q+1)∫Rnϑ(x)G(u0,v0,w0)dx>0. | (2.11) |
Then under condition (3.1), the functional Π(u,v,w)>0, ∀t>0.
Proof. By (2.11) and continuity, there exists a time t1>0 such that
Π(u,v,w)≥0,∀t<t1. |
Let
Y={(u,v,w)∣Π(u(t0),v(t0),w(t0))=0, Π(u,v,w)>0,∀t∈[0,t0)}. | (2.12) |
Then, by (2.9), (2.10), we have for all (u,v,w)∈Y,
Λ(u,v,w)=q−12(q+1)[(1−∫t0ϖ1(s)ds)‖u‖2H+(1−∫t0ϖ2(s)ds)‖v‖2H+(1−∫t0ϖ3(s)ds)‖w‖2H]+q−12(q+1)[(ϖ1∘u)+(ϖ2∘v)+(ϖ3∘w)]+1q+1Π(u,v,w)≥q−12(q+1)[l‖u‖2H+m‖v‖2H+ν‖w‖2H+(ϖ1∘u)+(ϖ2∘v)+(ϖ3∘w)]. |
Owing to (2.3), it follows for (u,v,w)∈Y
l‖u‖2H+m‖v‖2H+ν‖w‖2H≤2(q+1)q−1Λ(u,v,w)≤2(q+1)q−1E(t)≤2(q+1)q−1E(0). | (2.13) |
By (1.18), (3.1) we have
(q+1)∫RnG(u(t0),v(t0),w(t0))≤η(l‖u(t0)‖2H+m‖v(t0)‖2H+ν‖w(t0)‖2H)(q+1)/2≤η(2(q+1)q−1E(0))(q−1)/2(l‖u(t0)‖2H+m‖v(t0)‖2H+ν‖w(t0)‖2H)≤γ(l‖u(t0)‖2H+m‖v(t0)‖2H+ν‖w(t0)‖2H)<(1−∫t00ϖ1(s)ds)‖u(t0)‖2H+(1−∫t00ϖ2(s)ds)‖v(t0)‖2H+(1−∫t00ϖ3(s)ds)‖w(t0)‖2H<(1−∫t00ϖ1(s)ds)‖u(t0)‖2H+(1−∫t00ϖ2(s)ds)‖v(t0)‖2H+(1−∫t00ϖ3(s)ds)‖w(t0)‖2H+(ϖ1∘u)+(ϖ2∘v)+(ϖ3∘w), | (2.14) |
hence Π(u(t0),v(t0),w(t0))>0 on Y, which contradicts the definition of Y since Π(u(t0),v(t0),w(t0))=0. Thus Π(u,v,w)>0, ∀t>0.
The decay rate for solution is given in the next Theorem
Theorem 3.1. (Decay of solution) Let (1.2)–(1.17) and (1.11)–(1.14) hold. Under condition (2.1) and
γ=η(2(q+1)q−1E(0))(q−1)/2<1, | (3.1) |
there exists t0>0 depending only on ϖ1,ϖ2,ϖ3, λ1 and χ′(0) such that
0≤E(t)<E(t0)exp(−∫tt0ϖ(s)1−ϖ0(t)), | (3.2) |
holds for all t≥t0.
Proof. (Of Theorem 3.1.) By (1.18) and (2.13), we have for t≥0
0<l‖u‖2H+m‖v‖2H+ν‖w‖2H≤2(q+1)q−1E(t). | (3.3) |
Let
I(t)=ϖ(t)1−ϖ0(t), |
where ϖ and ϖ0 defined in (1.12) and (1.13).
Noting that limt→+∞ϖ(t)=0 by (1.11)–(1.13), we have
limt→+∞I(t)=0, I(t)>0, ∀t≥0. |
Then we take t0>0 such that
0<12I(t)<χ′(0), |
with (1.14) for all t>t0. Due to (2.3), we have
E(t)≤12(‖ut‖2L2ϑ+‖vt‖2L2ϑ+‖wt‖2L2ϑ)+12[(ϖ1∘u)+(ϖ2∘v)+(ϖ3∘w)]+12(1−∫t0ϖ1(s)ds)‖u‖2H+12(1−∫t0ϖ2(s)ds)‖v‖2H+12(1−∫t0ϖ3(s)ds)‖w‖2H≤12(‖ut‖2L2ϑ+‖vt‖2L2ϑ+‖wt‖2L2ϑ)+12[(ϖ1∘u)+(ϖ2∘v)+(ϖ3∘w)]+12(1−ϖ0(t))[‖u‖2H+‖v‖2H+‖w‖2H]. |
Then, by definition of I(t), we have
I(t)E(t)≤12I(t)(‖ut‖2L2ϑ+‖vt‖2L2ϑ+‖wt‖2L2ϑ)+12ϖ(t)[‖u‖2H+‖v‖2H+‖w‖2H]+12I(t)[(ϖ1∘u)+(ϖ2∘v)+(ϖ3∘w)], | (3.4) |
and Lemma 2.3, we have for all t1,t2≥0
E(t2)−E(t1)≤−12∫t2t1(ϖ(t)[‖u‖2H+‖v‖2H+‖w‖2H])dt+12∫t2t1((ϖ′1∘u)+(ϖ′2∘v)+(ϖ′3∘w))dt−α∫t2t1(‖ut‖2L2ϑ+‖vt‖2L2ϑ+‖wt‖2L2ϑ)dt, |
then,
E′(t)≤−12ϖ(t)[‖u‖2H+‖v‖2H+‖w‖2H]+12[(ϖ′1∘u)+(ϖ′2∘v)+(ϖ′3∘w)]−α(‖ut‖2L2ϑ+‖vt‖2L2ϑ+‖wt‖2L2ϑ), | (3.5) |
Finally, ∀t≥t0, we have
E′(t)+I(t)E(t)≤(12I(t)−α)(‖ut‖2L2ϑ+‖vt‖2L2ϑ+‖wt‖2L2ϑ)+12[(ϖ′1∘u)+(ϖ′2∘v)+(ϖ′3∘w)]+12I(t)((ϖ1∘u)+(ϖ2∘v)+(ϖ3∘w)), |
and we can choose t0>0 large enough such that
12I(t)<α, |
then
E′(t)+I(t)E(t)≤12∫t0{ϖ′1(t−τ)+I(t)ϖ2(t−τ)}‖u(t)−u(τ)‖2Hdτ+12∫t0{ϖ′2(t−τ)+I(t)ϖ2(t−τ)}‖v(t)−v(τ)‖2Hdτ+12∫t0{ϖ′3(t−τ)+I(t)ϖ3(t−τ)}‖w(t)−w(τ)‖2Hdτ≤12∫t0{ϖ′1(τ)+I(t)ϖ1(τ)}‖u(t)−u(t−τ)‖2Hdτ+12∫t0{ϖ′2(τ)+I(t)ϖ2(τ)}‖v(t)−v(t−τ)‖2Hdτ+12∫t0{ϖ′3(τ)+I(t)ϖ3(τ)}‖w(t)−w(t−τ)‖2Hdτ≤12∫t0{−χ(ϖ1(τ))+χ′(0)ϖ1(τ)}‖u(t)−u(t−τ)‖2Hdτ+12∫t0{−χ(ϖ2(τ))+χ′(0)ϖ2(τ)}‖v(t)−v(t−τ)‖2Hdτ+12∫t0{−χ(ϖ3(τ))+χ′(0)ϖ3(τ)}‖w(t)−w(t−τ)‖2Hdτ≤0, |
by the convexity of χ and (1.14), we have
χ(ξ)≥χ(0)+χ′(0)ξ=χ′(0)ξ. |
Then
E(t)≤E(t0)exp(−∫tt0I(s)ds), |
which completes the proof.
The author would like to thank the anonymous referees and the handling editor for their careful reading and for relevant remarks/suggestions to improve the paper.
The author agrees with the contents of the manuscript, and there is no conflict of interest among the author.
[1] | S. Angenent, On the formation of singularities in the curve shortening flow, J. Diff. Geom., 33 (1991), 601–633. |
[2] | P. Broadbridge, P. J. Vassiliouz, The role of symmetry and separation in surface evolution and curve shortening, SIGMA, 7 (2011), 1–19. |
[3] | K. S. Chou, X. P. Zhu, The curve shortening problem, Chapman: Hall/CRC, 2001. |
[4] | K. S. Chou, X. P. Zhu, Shortening complete plane curves, J. Diff. Geom., 50 (1998), 471–504. |
[5] |
P. W. Doyle, P. J. Vassiliou, Separation of variables for the 1-dimensional non-linear diffusion equation, Int. J. Non-Linear Mech., 33 (1998), 315–326. doi: 10.1016/S0020-7462(97)00013-9
![]() |
[6] | M. A. Grayson, The heat equation shrinks embedded plane curves to round points, J. Diff. Geom., 26 (1987), 285–314. |
[7] | M. E. Gage, R. Hamilton, The heat equation shrinking convex plane curves, J. Diff. Geom., 23 (1986), 69–96. |
[8] | M. E. Gurtin, Thermomechanics of evolving phase boundaries in the plane, New York: The Clarendon Press, 1993. |
[9] |
M. Nara, M. Taniguchi, The condition on the stability of stationary lines in a curvature flow in the whole plane, J. Diff. Eq., 237 (2007), 61–76. doi: 10.1016/j.jde.2007.02.012
![]() |
[10] |
D. H. Tsai, X. L. Wang, On some simple methods to derive the hairclip and paperclip solutions of the curve shortening flow, Acta Math. Sci., 39 (2019), 1674–1694. doi: 10.1007/s10473-019-0616-5
![]() |
1. | Zhehao Huang, Tianpei Jiang, Zhenzhen Wang, On a multiple credit rating migration model with stochastic interest rate, 2020, 43, 0170-4214, 7106, 10.1002/mma.6435 | |
2. | Zhehao Huang, Zhenghui Li, Zhenzhen Wang, Utility Indifference Valuation for Defaultable Corporate Bond with Credit Rating Migration, 2020, 8, 2227-7390, 2033, 10.3390/math8112033 | |
3. | Jie Xiong, Geng Deng, Xindong Wang, Extension of SABR Libor Market Model to handle negative interest rates, 2020, 4, 2573-0134, 148, 10.3934/QFE.2020007 | |
4. | Hiroyuki Kawakatsu, Recovering Yield Curves from Dynamic Term Structure Models with Time-Varying Factors, 2020, 3, 2571-905X, 284, 10.3390/stats3030020 | |
5. | João F. Caldeira, Werley C. Cordeiro, Esther Ruiz, André A.P. Santos, Forecasting the yield curve: the role of additional and time‐varying decay parameters, conditional heteroscedasticity, and macro‐economic factors, 2024, 0143-9782, 10.1111/jtsa.12769 |