Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

On an asymptotically log-periodic solution to the graphical curve shortening flow equation

  • 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 infty(x,t)=α,   lim supty(x,t)=β,   xK

    on any compact subset KR.

    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

    Related Papers:

    [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 infty(x,t)=α,   lim supty(x,t)=β,   xK

    on any compact subset KR.



    We consider, for xRn, t>0, the following system

    {uttθΔ(u+t0ϖ1(ts)u(s)ds)+αut=h1(u,v,w)vttθΔ(v+t0ϖ2(ts)v(s)ds)+αvt=h2(u,v,w)wttθΔ(w+t0ϖ3(ts)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 n3,α>0, the functions hi(.,.,.)(R3,R),i=1,2,3 are given by

    h1(ξ1,ξ2,ξ3)=(q+1)[d|ξ1+ξ2+ξ3|(q1)(ξ1+ξ2+ξ3)+e|ξ1|(q3)/2ξ1|ξ2|(q+1)/2];
    h2(ξ1,ξ2,ξ3)=(q+1)[d|ξ1+ξ2+ξ3|(q1)(ξ1+ξ2+ξ3)+e|ξ2|(q3)/2ξ2|ξ3|(q+1)/2];
    h3(ξ1,ξ2,ξ3)=(q+1)[d|ξ1+ξ2+ξ3|(q1)(ξ1+ξ2+ξ3)+e|ξ3|(q3)/2ξ3|ξ1|(q+1)/2],

    with d,e>0,q>3. The function 1θ(x)ϑ(x)>0 for all xRn is a density such that

    ϑLτ(Rn)withτ=2n2nrn+2rfor2r2nn2. (1.2)

    As in [15], it is not hard to see that there exists a function GC1(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 C0(Rn), as in [18], we have

    H={vL2nn2(Rn)v(L2(Rn))n},

    with respect to the norm vH=(v,v)1/2H for the inner product

    (v,w)H=Rnvwdx,

    and L2ϑ(Rn) as that to the norm vL2ϑ=(v,v)1/2L2ϑ for

    (v,w)L2ϑ=Rnϑvwdx.

    For general r[1,+)

    vLrϑ=(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Ω,t0u(x,0)=u0(x),ut(x,0)=u1(x),xΩ, (1.5)

    where Ω is a bounded domain of Rn, n1 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|q1|v|q+1u,f2=|v|q1|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+ωutt0g(ts)u(s)ds)=u|u|q1, (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(ts)Δu(s)ds+ut=f1(u,v)vttΔv+t0h(ts)Δ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(ts)Δu(s)ds+(|u|k+|v|q)|ut|m1ut=f1(u,v) vttΔv+t0h(ts)Δv(s)ds+(|v|κ+|u|ρ)|vt|r1vt=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=mi,j=1,ij|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

    v2nn2CτvH,

    and

    vL2ϑCPvH,

    for vH.

    Lemma 1.2. [12] Let ϑ satisfy (1.2), then the estimates

    vLrϑCrvH,

    and

    Cr=Cτϑ1rτ,

    hold for vH. Here τ=2n/(2nrn+2r) for 1r2n/(n2).

    We assume that the kernel functions ϖ1,ϖ2,ϖ3C1(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)=maxt0{ϖ1(t),ϖ2(t),ϖ3(t)}, (1.12)

    and

    ϖ0(t)=mint0{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ϑ(u2L(q+1)ϑ+v2L(q+1)ϑ)(q+1)/2(lu2H+mv2H)(q+1)/2, (1.15)

    and

    vw(q+1)/2L(q+1)/2ϑ(mv2H+νw2H)(q+1)/2, (1.16)

    and

    wu(q+1)/2L(q+1)/2ϑ(νw2H+lu2H)(q+1)/2. (1.17)

    Thanks to Minkowski's inequality to give

    u+v+w(q+1)L(q+1)ϑc(u2L(q+1)ϑ+v2L(q+1)ϑ+w2L(q+1)ϑ)(q+1)/2c(u2H+v2H+w2H)(q+1)/2.

    Then there exist η>0 such that

    u+v+w(q+1)L(q+1)ϑ+2uv(q+1)/2L(q+1)/2ϑ+2vw(q+1)/2L(q+1)/2ϑ+2wu(q+1)/2L(q+1)/2ϑη(lu2H+mv2H+νw2H)(q+1)/2. (1.18)

    We need to define positive constants λ0 and E0 by

    λ0η1/(q1)andE0=(121q+1)η2/(q1). (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)}iNR×H of

    θ(x)Δei=λieixRn,

    for any iN, (θ(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 xRn,

    {Rnϑ(x)(utt+αut)φdx+Rnuφdxt0ϖ1(ts)u(s)dsφdx=Rnϑ(x)h1(u,v,w)φdx,Rnϑ(x)(vtt+αvt)ψdx+Rnvψdxt0ϖ2(ts)v(s)dsψdx=Rnϑ(x)h2(u,v,w)ψdx,Rnϑ(x)(wtt+αwt)Ψdx+RnvΨdxt0ϖ3(ts)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<qn+2n2andthatn3. (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, XTC([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:H3R defined by

    J(u,v,w)=(1t0ϖ1(s)ds)u2H+(ϖ1u)+(1t0ϖ2(s)ds)v2H+(ϖ2v)+(1t0ϖ3(s)ds)w2H+(ϖ3w), (2.2)

    where

    (ϖjw)(t)=t0ϖj(ts)w(t)w(s)2Hds,

    for any wL2(Rn),j=1,2,3. The modified energy is defined by

    E(t)=12(ut2L2ϑ+vt2L2ϑ+wt2L2ϑ)+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, XC([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 0t1<t2T, we have

    E(t2)E(t1)=t2t1ddtE(t)dt=12t2t1(ϖ1(t)u2H(ϖ1u))dt12t2t1(ϖ2(t)v2H(ϖ2v))dt12t2t1(ϖ3(t)w2H(ϖ3w))dtα(ut2L2ϑ+vt2L2ϑ+wt2L2ϑ)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(ts)Δu(s)ds+ϑ(x)h1(u,v,w)ϑ(x)(ytt+αyt)Δy=t0ϖ2(ts)Δv(s)ds+ϑ(x)h2(u,v,w)ϑ(x)(ζtt+αζt)Δζ=t0ϖ3(ts)Δ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)dx12J(u,v,w)1q+1u+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[lu2H+mv2H+νw2H](q+1)/212J(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 0E(0)<E0.

    (i) If u02H+v02H+w02H<λ20, then local solution of (1.1) satisfies

    J(u,v,w)<λ20, t[0,T).

    (ii) If u02H+v02H+w02H>λ20, then local solution of (1.1) satisfies

    u2H+v2H+w2H>λ21, t[0,T),λ1>λ0.

    Proof. Since 0E(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

    u02H+v02H+w02Hξ22,

    and

    u2H+v2H+w2Hξ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 0E(0)<E0 and

    u02H+v02H+w02H<λ20.

    By Lemma 2.3 and Lemma 2.4, we have

    12(ut2L2ϑ+vt2L2ϑ+wt2L2ϑ)+lu2H+mv2H+νw2H12(ut2L2ϑ+vt2L2ϑ+wt2L2ϑ)+(1t0ϖ1(s)ds)u2H+(ϖ1u)+(1t0ϖ2(s)ds)u2H+(ϖ2v)+(1t0ϖ3(s)ds)w2H+(ϖ3w)2E(t)+2ηq+1[lu2H+mu2H+νw2H](q+1)/22E(0)+2ηq+1(J(u,v,w))(q+1)/22E0+2ηq+1λq+10=η2/(q1). (2.8)

    This completes the proof.

    Let

    Λ(u,v,w)=12(1t0ϖ1(s)ds)u2H+12(ϖ1u)+12(1t0ϖ2(s)ds)v2H+12(ϖ2v)+12(1t0ϖ3(s)ds)w2H+12(ϖ3w)Rnϑ(x)G(u,v,w)dx, (2.9)
    Π(u,v,w)=(1t0ϖ1(s)ds)u2H+(ϖ1u)+(1t0ϖ2(s)ds)v2H+(ϖ2v)+(1t0ϖ3(s)ds)w2H+(ϖ3w)(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

    u02H+v02H+w02H(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)=q12(q+1)[(1t0ϖ1(s)ds)u2H+(1t0ϖ2(s)ds)v2H+(1t0ϖ3(s)ds)w2H]+q12(q+1)[(ϖ1u)+(ϖ2v)+(ϖ3w)]+1q+1Π(u,v,w)q12(q+1)[lu2H+mv2H+νw2H+(ϖ1u)+(ϖ2v)+(ϖ3w)].

    Owing to (2.3), it follows for (u,v,w)Y

    lu2H+mv2H+νw2H2(q+1)q1Λ(u,v,w)2(q+1)q1E(t)2(q+1)q1E(0). (2.13)

    By (1.18), (3.1) we have

    (q+1)RnG(u(t0),v(t0),w(t0))η(lu(t0)2H+mv(t0)2H+νw(t0)2H)(q+1)/2η(2(q+1)q1E(0))(q1)/2(lu(t0)2H+mv(t0)2H+νw(t0)2H)γ(lu(t0)2H+mv(t0)2H+νw(t0)2H)<(1t00ϖ1(s)ds)u(t0)2H+(1t00ϖ2(s)ds)v(t0)2H+(1t00ϖ3(s)ds)w(t0)2H<(1t00ϖ1(s)ds)u(t0)2H+(1t00ϖ2(s)ds)v(t0)2H+(1t00ϖ3(s)ds)w(t0)2H+(ϖ1u)+(ϖ2v)+(ϖ3w), (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)q1E(0))(q1)/2<1, (3.1)

    there exists t0>0 depending only on ϖ1,ϖ2,ϖ3, λ1 and χ(0) such that

    0E(t)<E(t0)exp(tt0ϖ(s)1ϖ0(t)), (3.2)

    holds for all tt0.

    Proof. (Of Theorem 3.1.) By (1.18) and (2.13), we have for t0

    0<lu2H+mv2H+νw2H2(q+1)q1E(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, t0.

    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(ut2L2ϑ+vt2L2ϑ+wt2L2ϑ)+12[(ϖ1u)+(ϖ2v)+(ϖ3w)]+12(1t0ϖ1(s)ds)u2H+12(1t0ϖ2(s)ds)v2H+12(1t0ϖ3(s)ds)w2H12(ut2L2ϑ+vt2L2ϑ+wt2L2ϑ)+12[(ϖ1u)+(ϖ2v)+(ϖ3w)]+12(1ϖ0(t))[u2H+v2H+w2H].

    Then, by definition of I(t), we have

    I(t)E(t)12I(t)(ut2L2ϑ+vt2L2ϑ+wt2L2ϑ)+12ϖ(t)[u2H+v2H+w2H]+12I(t)[(ϖ1u)+(ϖ2v)+(ϖ3w)], (3.4)

    and Lemma 2.3, we have for all t1,t20

    E(t2)E(t1)12t2t1(ϖ(t)[u2H+v2H+w2H])dt+12t2t1((ϖ1u)+(ϖ2v)+(ϖ3w))dtαt2t1(ut2L2ϑ+vt2L2ϑ+wt2L2ϑ)dt,

    then,

    E(t)12ϖ(t)[u2H+v2H+w2H]+12[(ϖ1u)+(ϖ2v)+(ϖ3w)]α(ut2L2ϑ+vt2L2ϑ+wt2L2ϑ), (3.5)

    Finally, tt0, we have

    E(t)+I(t)E(t)(12I(t)α)(ut2L2ϑ+vt2L2ϑ+wt2L2ϑ)+12[(ϖ1u)+(ϖ2v)+(ϖ3w)]+12I(t)((ϖ1u)+(ϖ2v)+(ϖ3w)),

    and we can choose t0>0 large enough such that

    12I(t)<α,

    then

    E(t)+I(t)E(t)12t0{ϖ1(tτ)+I(t)ϖ2(tτ)}u(t)u(τ)2Hdτ+12t0{ϖ2(tτ)+I(t)ϖ2(tτ)}v(t)v(τ)2Hdτ+12t0{ϖ3(tτ)+I(t)ϖ3(tτ)}w(t)w(τ)2Hdτ12t0{ϖ1(τ)+I(t)ϖ1(τ)}u(t)u(tτ)2Hdτ+12t0{ϖ2(τ)+I(t)ϖ2(τ)}v(t)v(tτ)2Hdτ+12t0{ϖ3(τ)+I(t)ϖ3(τ)}w(t)w(tτ)2Hdτ12t0{χ(ϖ1(τ))+χ(0)ϖ1(τ)}u(t)u(tτ)2Hdτ+12t0{χ(ϖ2(τ))+χ(0)ϖ2(τ)}v(t)v(tτ)2Hdτ+12t0{χ(ϖ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
  • This article has been cited by:

    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
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2554) PDF downloads(116) Cited by(0)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog