Processing math: 100%
Research article Special Issues

Asymptotic analysis for non-local curvature flows for plane curves with a general rotation number

  • Several non-local curvature flows for plane curves with a general rotation number are discussed in this work. The types of flows include the area-preserving flow and the length-preserving flow. We have a relatively good understanding of these flows for plane curves with the rotation number one. In particular, when the initial curve is strictly convex, the flow exists globally in time, and converges to a circle as time tends to infinity. Even if the initial curve is not strictly convex, a global solution, if it exists, converges to a circle. Here, we deal with curves with a general rotation number, and show, not only a similar result for global solutions, but also a blow-up criterion, upper estimates of the blow-up time, and blow-up rate from below. For this purpose, we use a geometric quantity which has never been considered before.

    Citation: Takeyuki Nagasawa, Kohei Nakamura. Asymptotic analysis for non-local curvature flows for plane curves with a general rotation number[J]. Mathematics in Engineering, 2021, 3(6): 1-26. doi: 10.3934/mine.2021047

    Related Papers:

    [1] Daniel N. Riahi, Saulo Orizaga . On rotationally driven nonlinear inclined polymeric jet with gravity effect. Mathematics in Engineering, 2022, 4(2): 1-18. doi: 10.3934/mine.2022014
    [2] Tatsuya Miura . Polar tangential angles and free elasticae. Mathematics in Engineering, 2021, 3(4): 1-12. doi: 10.3934/mine.2021034
    [3] Bruno Bianchini, Giulio Colombo, Marco Magliaro, Luciano Mari, Patrizia Pucci, Marco Rigoli . Recent rigidity results for graphs with prescribed mean curvature. Mathematics in Engineering, 2021, 3(5): 1-48. doi: 10.3934/mine.2021039
    [4] Qiang Guang, Qi-Rui Li, Xu-Jia Wang . Flow by Gauss curvature to the $ L_p $ dual Minkowski problem. Mathematics in Engineering, 2023, 5(3): 1-19. doi: 10.3934/mine.2023049
    [5] Francesca Tedeschi, Giulio G. Giusteri, Leonid Yelash, Mária Lukáčová-Medvid'ová . A multi-scale method for complex flows of non-Newtonian fluids. Mathematics in Engineering, 2022, 4(6): 1-22. doi: 10.3934/mine.2022050
    [6] Eleonora Cinti, Roberto Ognibene, Berardo Ruffini . A quantitative stability inequality for fractional capacities. Mathematics in Engineering, 2022, 4(5): 1-28. doi: 10.3934/mine.2022044
    [7] Jerome D. Wettstein . Half-harmonic gradient flow: aspects of a non-local geometric PDE. Mathematics in Engineering, 2023, 5(3): 1-38. doi: 10.3934/mine.2023058
    [8] Virginia Agostiniani, Lorenzo Mazzieri, Francesca Oronzio . A geometric capacitary inequality for sub-static manifolds with harmonic potentials. Mathematics in Engineering, 2022, 4(2): 1-40. doi: 10.3934/mine.2022013
    [9] Alessandra De Luca, Veronica Felli . Unique continuation from the edge of a crack. Mathematics in Engineering, 2021, 3(3): 1-40. doi: 10.3934/mine.2021023
    [10] Anoumou Attiogbe, Mouhamed Moustapha Fall, El Hadji Abdoulaye Thiam . Nonlocal diffusion of smooth sets. Mathematics in Engineering, 2022, 4(2): 1-22. doi: 10.3934/mine.2022009
  • Several non-local curvature flows for plane curves with a general rotation number are discussed in this work. The types of flows include the area-preserving flow and the length-preserving flow. We have a relatively good understanding of these flows for plane curves with the rotation number one. In particular, when the initial curve is strictly convex, the flow exists globally in time, and converges to a circle as time tends to infinity. Even if the initial curve is not strictly convex, a global solution, if it exists, converges to a circle. Here, we deal with curves with a general rotation number, and show, not only a similar result for global solutions, but also a blow-up criterion, upper estimates of the blow-up time, and blow-up rate from below. For this purpose, we use a geometric quantity which has never been considered before.


    In this paper, we deal with curvature flows comprising non-local terms for plane curves with a general rotation number. Let f be an R2-valued function on R/L(t)Z×[0,T) such that for a fixed t[0,T), it is an arc-length parametrization of a closed plane curve with total length L(t). In the following text, we simply denote L(t) as L in many cases. To explain the curvature flow that is considering in this work, we introduce a certain geometric quantity. For a fixed t[0,T), sR/LZ is an arc-length parameter. Then, τ=sf and κ=2sf are the unit tangent vector and the curvature vector respectively. The vector ν is a unit normal vector given by rotating τ counter-clockwise by π2. The curvature κ and its deviation ˜κ are given by

    κ=κν,˜κ=κ1LL0κds.

    Here, ˜κ is a non-local quantity. The equation we consider is of the following form:

    tf=(˜κgL)ν.

    Here, we assume that the function g is a scale-invariant non-local quantity determined by f. That is, set fλ(s)=1λf(λs) (sR/λ1LZ), then,

    g(fλ)=g(f).

    Here we study three cases of g:

    (AP) If we set g0, then our equation represents the area-preserving flow. In fact, we set A as

    A=12L0fνds

    which is the enclosed area when Imf is a simple curve. Consequently, it holds that

    dAdt=0.

    (LP) Let g=L(L0κds)1L0˜κ2ds. Here, the equation represents the length-preserving flow:

    dLdt=0.

    (JP) Jiang-Pan considered an equation with g=L22AL0κds in [5]. Here, the isoperimetric ratio does not increase along with the flow:

    ddtL2A=2LAL0tf2ds.

    Let

    n=12πL0κds

    be the rotation number. For classical solutions, the rotation number n is independent of t. There are a multitude of literature available considering the case when n=1 in the above equations. First of all, we should mention Gage's result [3]. Assume that Imf(0) is a strictly convex, closed curve with a rotation number equal to 1 in the class of C2. Then, the solution f with the initial data f(0) exists globally in time, and Imf(t) converges to a circle with a surrounding area A(0) as t. Similar results for (LP) and (JP) were proved by [6] and [5] respectively under the convexity condition. The authors considered flows without the convexity condition in [7,8]. Instead of convexity, we assume the global existence of the solution. Then the solution of (AP), (LP), or (JP) converges to a circle as t exponentially. As a result, the curvature uniformly converges to a positive constant, and thus, the curve becomes convex in finite time. In our previous works, the isoperimetric deficit

    I1=14πAL2

    played an important role. First, we show the decay of I1. Set

    I=L2+1L0|˜κ()|2ds for {0}N.

    In [7], we showed the inequality

    IjC(Ij21I+Ij+11Ij+1+1) (1.1)

    for an integer j[0,] with a positive constant C=C(j,) independent of the total length of curve. Since I1 is small for a sufficiently large t, we can regard this inequality as an embedding with a small embedding constant. We showed the exponential decay of I using the standard energy method, combining the above inequality. Finally, using the decay of I, we showed the convergence of Imf to a circle.

    In this paper we study the case of n>1, when the isomerimetric deficit is

    I1=14nπAL2.

    The isoperimetric inequality shows I10 when n=1. However, I1 is not necessarily non-negative for n>1. This implies the technique used in [7,8] is not applicable for n>1. In spite of this, I1 gives us some useful information. For example, we can show that if I1 is negative for t=0, then the solution blows up in finite time. See our first main result, Theorem 3.1. This implies I10 for global solutions, and that sounds a good information. However, the inequality (1.1) does not hold for n>1. There are at least two approaches for dealing with this difficulty. One is to give a proof without using (1.1), and another is to show an alternative inequality to (1.1). In this paper, we show that both are in success. For the second approach, we use a geometric quantity which has never been considered before, given as follows:

    ˜I1=1L2πnL(f1LL0fds)+ν2L2.

    Then we can show

    IjC(˜Ij21I+˜Ij+11Ij+1+1). (1.2)

    We prepare several inequalities and estimates for closed curves with a rotation number n, in § 2.1. And we describe some basic properties of the flows (AP), (LP) and (JP), in § 2.2. Using these, in § 3, we discuss blow-up solutions with blow-up time estimates, blow-up quantities, and blow-up rates. In § 4, the convergence to an n-fold circle of global solutions is proved without using (1.2). Finally, we show (1.2) in the final section.

    In this section, we provide several estimates and inequalities for plane curves. Those in § 2.1 hold for curves which are not necessarily solutions of the flows. We derive the basic properties of flows in § 2.2.

    Let f=(f1,f2) be an arc-length parametrization of a plane curve with the rotation number n1. Set

    f=f1+if2,ν=ν1+iν2=f2+if1=if.

    The functions φk=1Lexp(2πiksL) for kZ generate the standard complete orthogonal system of L2(R/LZ). Let ˆf(k) be the Fourier coefficient of f. Subsequently, we can derive the following relations in a manner similar to [7,Corollary 2.1], where we dealt with the case of n=1. The difference is just "n" in (2.3) which comes exactly from the definition of the rotation number. We can find similar argument in [1,10]

    Lemma 2.1.

    kZk|ˆf(k)|2= LAπ, (2.1)
    kZk2|ˆf(k)|2= (L2π)2L0κ0ds=L34π2, (2.2)
    kZk3|ˆf(k)|2= (L2π)3L0κds=nL34π2, (2.3)
    kZk4|ˆf(k)|2= (L2π)4L0κ2ds, (2.4)
    kZk5|ˆf(k)|2= (L2π)5L0κ3ds, (2.5)
    kZk6|ˆf(k)|2= (L2π)6L0{κ4+(κ)2}ds. (2.6)

    Note that we have

    kZk2(kn)|ˆf(k)|2=0 (2.7)

    from (2.2) and (2.3). The above is very useful for our analysis.

    Lemma 2.2. We have

    I0= 16π4L3kZk3(kn)|ˆf(k)|2 (2.8)
    = 16π4L3kZk2(kn)2|ˆf(k)|2. (2.9)

    Proof. We obtain (2.8) as

    I0= LL0˜κ2ds=LL0˜κκds=L(L0κ2ds2πnLL0κds)= 16π4L3kZk3(kn)|ˆf(k)|2

    from (2.4) and (2.3). Combining this with (2.7), we obtain (2.9).

    Though I0 must be non-negative by the definition, it is not obvious to see that from the first expression (2.8). However, it can be seen from the second one (2.9). Furthermore, we see from (2.9) that I0=0 if and only if Imf is an n-fold circle.

    The isoperimetric inequality holds even if n is not 1.

    Lemma 2.3. We have L24πA0.

    Proof. It follows from (2.2) and (2.1) that

    L24πA=4π2L(L34π2LAπ)=4π2LkZk(k1)|ˆf(k)|20.

    Similarly, I1 has two expressions.

    Lemma 2.4. We have

    I1=4π2L3kZk(kn)|ˆf(k)|2=4π2nL3kZ{0}k(kn)2|ˆf(k)|2.

    Proof. It follows from (2.2) and (2.1) that

    I1=14πnAL2=4π2L3(L34π2nLAπ)=4π2L3kZk(kn)|ˆf(k)|2.

    The second expression of I1 is obtained from the above and (2.7).

    Since k(kn) is not necessarily non-negative when n>1, we know the same holds for I1. However, the modulus of I1 can be estimated by I0 for n1 as follows. This is Wirtinger's inequality when n=1.

    Lemma 2.5. It holds that 4π2n|I1|I0.

    Proof. From Lemmas 2.2–2.4 we obtain

    I0±4π2nI1= 16π4L3kZ{k2(kn)2k(kn)2}|ˆf(k)|2= 16π4L3kZk(k1)(kn)2|ˆf(k)|20.

    Here, we use k(k1)0 for kZ.

    In this subsection, we derive the basic properties of the flows, which we use in following sections. Let f be a classical solution of one of (AP), (LP), or (JP) on [0,T), and let T be the maximum existence time. Since dLdt=L0tfκds, we have

    dL2dt=2LL0(˜κgL)κds=2LL0˜κ2ds+4πng,

    that is,

    dL2dt+2I0=4πng. (2.10)

    Similarly, we have

    dAdt=L0tfνds=L0(˜κgL)ds=g. (2.11)

    It follows from the above that

    ddt(L2I1)+2I0=ddt(L24πnA)+2I0=0. (2.12)

    From these, we summarize the basic properties of each solution as follows.

    Proposition 2.1. Assume that the initial curve is smooth, and that A(0) is positive. Let f be a classical solution of one of (AP), (LP), or (JP) on [0,T) and let T be the maximum existence time. Then, the following holds for t(0,T).

    1). For solutions of (AP),

    dAdt=0,AA(0)>0,dL2dt0,dI1dt0.

    2). For solutions of (LP),

    dAdt0,AA(0)>0,dL2dt=0,dI1dt0.

    3). For solutions of (JP),

    A>0,dI1dt0.

    4). For solutions of (AP), (LP), (JP),

    1nI1I1(0).

    In other words,

    4πL2AL(0)2A(0).

    Proof. In the cases of (AP) and (LP), the signs of dAdt and dL2dt immediately follow from (2.11) and (2.10). Therefore, A>0 and

    dI1dt=ddt4πnAL2=4πnL2dAdt+4πnAL4dL2dt0.

    In the case of (JP), we prove the positivity of A by applying the contradiction argument. In this case,

    g=L2I12A. (2.13)

    It follows from (2.11) that

    dA2dt=2Ag=L2I1. (2.14)

    Assume that A(t0)2=0 for some first time t0(0,T). Since A20, we have

    dA2dt(t0)=0. (2.15)

    Since A(0)2>0, there exists t1(0,t0) such that

    dA2dt(t1)<0. (2.16)

    It follows from (2.14) and (2.12) that

    d2A2dt2=ddt(L2I1)=2I00.

    Therefore, by (2.16)

    dA2dt(t0)=dA2dt(t1)+t0t1d2A2dt2dtdA2dt(t1)<0.

    This contradicts (2.15). Hence, A>0 on (0,T). Using (2.12), (2.10), I11=4πnAL2, and (2.13), we have

    L2dI1dt= I1dL2dt2I0= I1(4πng2I0)2I0=4πngI1+2(I11)I0= 4πnL2(L2gI1+2AI0)=4πnL2(L4I212A+2AI0)0.

    Since I1 is non-increasing, we have I1I1(0). Lemma 2.3 gives us

    I1=14πnAL2=1n+n(14πAL2)1n.

    The non-positivity of I1(0) implies that the blow-up phenomena occurs in finite time.

    Theorem 3.1. Let f be a classical solution of one of (AP), (LP), or (JP) on [0,T) and let T be the maximum existence time. Assume that the initial curve is smooth, and satisfies A(0)>0, I1(0)<0. Then, the solution blows up in finite time. The blow-up time T is estimated from above as follows:

    (AP)TL(0)24πA(0)8π2nI1(0),

    (LP)TL(0)24πA(0)8π2I1(0),

    (JP)TL(0)28π2nI1(0).

    Proof. In the case of (AP), g0. It follows from Proposition 2.1 that I1(t)I1(0)<0. By (2.10) and Lemma 2.5, we have

    dL2dt=2I0(t)8π2nI1(t)8π2nI1(0).

    Integrating this from 0 to t(0,T), and using Lemma 2.3, we obtain

    4πA(0)L2(0)=4πA(t)L2(0)L2(t)L2(0)8π4nI1(0)t.

    Since the first side is non-positive by the isoperimetric inequality (Lemma 2.3), t must satisfy

    tL(0)24πA(0)8π2nI1(0).

    In the case of (LP), g=I02πn0. Proposition 2.1 shows I1(t)I1(0)<0. From (2.11) and Lemma 2.5, we have

    dAdt=12πnI0(t)2πI1(t)2πI1(0).

    We integrate this from 0 to t(0,T). Using Lemma 2.3, we obtain

    4πA(0)L(0)2=4πA(0)L(t)24π(A(0)A(t))8πI1(0)t.

    Consequently, t must satisfy

    tL(0)24πA(0)8π2I1(0).

    In the case of (JP), g=L2I12A. It follows from (2.10), Proposition 2.1, and Lemma 2.5 that

    dL2dt=2I0(t)+2πnL(t)2A(t)I1(t)2I0(t)8π2nI1(t)8π2nI1(0).

    We integrate this from 0 to t(0,T). Using Lemma 2.3, we obtain

    L(0)2L(t)2L(0)28π2nI1(0)t.

    Consequently t must satisfy

    tL(0)28π2nI1(0).

    Corollary 3.1. Let f be a classical solution of one of (AP), (LP), or (JP) on [0,T) and let T be the maximum existence time. Assume that the initial curve is smooth, and that satisfies A(0)>0, and I1(0)=0, but it is not an n-fold circle. Then, T<.

    Proof. Assume T=. Then, Theorem 3.1 implies that I1(t)0 for all t[0,). On the other hand, (2.12) with I1(0)=0 shows that I1(t)0. Hence, I1(t)0. When t>0,

    L0˜κ2ds=I0L=12Lddt(L2I1)=0.

    Combining this with the rotation number n, we find that Imf(t) is an n-fold circle. However, this does not satisfy the initial condition.

    Corollary 3.2. f is a classical stationary solution of one of (AP), (LP), or (JP), if and only of it is an n-fold circle.

    Proof. Assume that Imf is an n-fold circle. Then, ˜κ0. Since f=ˆf(0)φ0+ˆf(n)φn, we see I0=I1=0 by Lemmas 2.2 and 2.4. Hence, ˜κgL0 for each case. Consequently, it is a stationary solution.

    Conversely, assume that f is a stationary solution. It follows from (2.12) that I0(t)0. Hence, we can conclude that Imf(t) is an n-fold circle in a manner similar to the proof of the previous corollary.

    Suppose now f blows up as tT<. Then, we have

    lim suptTI0(t)=.

    Indeed, if lim suptTI0(t)<, then supt(0,T)I0(t) is bounded. We can show the boundedness of supt(0,T)I(t) by the standard energy method. Using this and the equation of the flow, we can see that f(t) converges to a smooth function as tT. Consequently, the solution can be expanded beyond T. This is a contradiction.

    Set

    W=L0κ2ds.

    We will show the blow-up of W and its blow-up rate. Firstly, we consider the limit supremum of W.

    Lemma 3.1. It holds that lim suptTW(t)=.

    Proof. Set

    R=L0κds,

    and we have

    LW=LL0{˜κ2+(RL)2}ds=I0+R2.

    Hence,

    lim suptTL(t)W(t)=.

    Therefore, the assertion immediately follows in the case of (LP).

    In the case of (AP), L is non-increasing by Proposition 2.1. Lemma 2.3 implies that L4πA=4πA0. Consequently, L(t) converges to a positive constant as tT, and the assertion follows.

    We show that L(t) converges to a positive constant in the case of (JP) as well. We assume that lim inftTA(t)=0. I1 is monotone by Proposition 2.1. Therefore, it follows from

    dAdt=L22AI1

    that A does not oscillate near t=T. Hence, we may assume limtTA(t)=0. From the above relation and Proposition 2.1, we find that dAdt is bounded. Consequently, the estimate

    0<A(t)C(Tt)

    holds. Thus, we have

    0A(t)2TtC(Tt)2Tt0 as tT,

    and therefore,

    limtTA(T0)2A(t)2Tt=0.

    This implies that the left derivative of A2 at T vanishes:

    dA2dt(T0)=0. (3.1)

    However, A(0)2>0 and A(T0)2=0 show the existence of t(0,T) such that

    dA2dt(t)<0.

    Since

    d2A2dt2=2I00,

    we have

    dA2dt(t)<dA2dt(t)<0

    for t(t,T). This contradicts (3.1). Now, we prove lim inftTA(t)>0. Since

    dAdt=L22AI1

    has a constant sign near T, we conclude that limtTA(t)>0. The convergence of limtTL(t) follows from the convergence of A, and the monotonicity and boundedness of I1. Since L2A is strictly positive by Proposition 2.1, the limit of L is positive.

    Next, we derive the time derivative of W. Set

    Jp=Lp1L0˜κpds(pN{1}),

    which are scale-invariant quantities. Note that I0=J2.

    Lemma 3.2. It holds that

    dWdt=1L3{2I1+J4+(3Rg)J3+3R(Rg)J2R3g}.

    Proof. The proof is a direct calculation:

    dWdt= L0tf(22sκ+κ3)ds=L0(˜κgL)(22sκ+κ3)ds= 2L0(s˜κ)2ds+L0(˜κgL)(˜κ+RL)3ds= 2I1L3+L0(˜κ3+3R˜κ2L+3R2˜κL2+R3L3)(˜κgL)ds= 2I1L3+L0{˜κ4+(3RLgL)˜κ3+(3R2L23RgL2)˜κ2R3gL4}ds= 1L3{2I1+J4+(3Rg)J3+3R(Rg)J2R3g}.

    Thirdly, we estimate dWdt from above.

    Lemma 3.3. We have

    dWdtW32M2.

    Here,

    M={Cfor (AP) and (LP),C{1+(L20A0)43}12for (JP)

    with the constant C being independent of the initial curve and the rotation number.

    Proof. Here, we use Lemma 3.2. In the case of (AP), since g=0, we have

    dWdt+2I1L3=1L3(J4+3RJ3+3R2J2).

    Set θ=121p. Then, Gagliardo-Nirenberg's inequality yields

    |Jp|C(I1θ0Iθ1)p2=CIp4+120Ip4121.

    Hence,

    dWdt+2I1L3 CL3(I320I121+RI540I141+R2I0) I1L3+CL3(I30+R43I530+R2I0).

    Since 0I0LW and R2LW, we obtain

    I30L3W3,I530L53W53=(LW)43L3W3R83L3W3,I0LW=(LW)4L3W3R8L3W3.

    Furthermore,

    R=2πn2π.

    Consequently, we conclude that

    dWdtC(1+R43+R6)W3CW3.

    In the case of (LP), since g=I0R, we have

    dWdt+1L3(2I1+3I20+R2I0)= 1L3{J4+(3RI0R)J3+3R2I0} CL3(I320I121+RI540I141+R1I940I141+R2I0) I1L3+CL3(I30+R43I530+R43I30+R2I0) I1L3+C(1+R43+R6)W3 I1L3+CW3.

    In the case of (JP), since g=L22AR, we have

    dWdt+1L3(2I1+3RL22AI0+R3L22A)=1L3{J4+(3RL22A+R)J3+6R2J2+R4}CL3{I320I121+RI540I141+L2AI540I141+R2I0+R2(LW)3}I1L3+CL3[I30+{R+(L2A)}43I530+R2I0+R2L3W3]I1L3+C[1+{1R+(L2R2A)}43+R6+R2]W3I1L3+C{1+(L2A)43}W3.

    By Proposition 2.1, we have

    (L2A)43(L20A0)43.

    Consequently, we can conclude that

    dWdtC{1+(L20A0)43}W3.

    Now, we prove the following theorem.

    Theorem 3.2. Let T(0,T) be the blow-up time for a solution of one of (AP), (LP), or (JP). Then, W(t) blows up as

    W(t)MTt,

    where

    M={Cfor (AP) and (LP),C{1+(L20A0)43}12for (JP)

    with a constant C that is independent of the initial curve and the rotation number.

    Proof. It follows from Lemma 3.3 that

    ddtW2M2.

    Due to Lemma 3.1, there exists a sequence {tn} such that tntT and W(tn)20 as n. Integrating the differential inequality from t to tn, we have

    W(t)2W(tn)2M2(tnt).

    Therefore, we obtain the theorem as n.

    The curve Imf may have several loops. When the orientation of a loop is counter-clockwise as s increases, it is called a positive loop. Otherwise, it is called a negative loop. It has already been shown that L(t) converges to a positive constant as t. Therefore, from the above theorem we know that

    limtTmaxsR/L(t)Zκ(s,t)=

    or

    limtTminsR/L(t)Zκ(s,t)=.

    If a positive/negative loop of Imf shrinks as tT, the maximum/minimum value of the curvature may not remain bounded. On the other hand, there is a possibility of the maximum or minimum remaining bounded as tT. For example, if a negative loop shrinks as tT before the positive loops shrink, the minimum value of the curvature goes to , but the maximum remains bounded. In the last part of this section, we discuss the blow-up of the maximum and minimum.

    Theorem 3.3. Let T(0,) be the blow-up time for a solution of one of (AP), (LP), or (JP). Assume that

    lim suptTmaxsR/L(t)Zκ(s,t)=,

    then it satisfies

    maxsR/L(t)Zκ(s,t)12(Tt).

    Proof. Set

    K(t)= maxsR/L(t)Zκ(s,t),d+Kdt(t)= lim suph+0K(t+h)K(t)h.

    Define the set St by St={sR/L(t)Z|κ(s,t)=K(t)}. After re-parametrizing f(,t) by a new parameter that is independent of t, we apply [2,Lemma B.40]. Consequently, we can conclude that K is a continuous function of t, and that

    d+Kdt(t)=maxsSttκ(s,t).

    κ satisfies the equation

    tκ=2sκ+κ2(˜κgL)=2sκ+κ2(κR+gL).

    For the cases of (AP) and (LP), R+g>0 as R>0 and g0. In the case of (JP),

    R+g=L2A0.

    2sκ0 holds for sSt. Hence, we have

    2sκ+κ2(κR+gL)κ3=K3

    for sSt, and

    d+Kdt(t)maxsSttκK3(t).

    We calculate Dini's derivative of K2 as

    d+dtK2(t)= lim suph+0K2(t+h)K2(t)h= lim suph+0(K(t)+K(t+h))(K(t)K(t+h))K2(t+h)K2(t)h= 2K3(t)lim infh+0K(t+h)K(t)h 2K3(t)lim suph+0K(t+h)K(t)h= 2K3(t)d+Kdt(t)2.

    According to the assumption of the theorem, there exists a sequence {tk}kN such that tkT and K(tk)20 as k. Using [4,Theorem 3], we have

    K2(tk)K2(t)_  tktd+dtK2(t)dt2(tkt)

    for tk(t,T). Therefore, we can conclude that

    K2(t)2(Tt)

    by k

    Theorem 3.4. Let T(0,) be the blow-up time for a solution of one of (AP), (LP), or (JP). Assume that

    supt[0,T)maxsR/L(t)Zκ(s,t)<.

    For the solution of (AP),

    minsR/L(t)Zκ(s,t)14(Tt)

    holds.

    For the solution of (LP),

    minsR/L(t)Zκ(s,t){2πn9L(0)(Tt)}13

    holds.

    For the solution of (JP), there exists a time T[0,T) such that

    minsR/L(t)Zκ(s,t)maxsR/L(t)Zκ(s,t)

    holds for t[T,T). Additionally, it holds that

    minsR/L(t)Zκ(s,t)12C(Tt),

    where

    C=1+L(T)24πnA(T).

    Remark 3.1. The time T above exists for all cases. And for the proof, it does not need to be the first or last such time.

    Proof. Here, we set

    K(t)= minsR/L(t)Zκ(s,t),d+Kdt(t)= lim suph+0K(t+h)K(t)h.

    Define the set St by St={sR/L(t)Z|κ(s,t)=K(t)}. As shown before, it holds that

    d+Kdt(t)=maxsStt(κ).

    κ satisfies

    t(κ)=2s(κ)+(κ)2{(κ)+R+gL}.

    Since 2s(κ)0 and κ=K for sSt,

    t(κ)K3+(R+g)K2L.

    If κC< holds on [0,T), then,

    Lmax{C2+K2}L0κ2ds=W as tT

    by Theorem 3.2. Since L is bounded, we conclude that K as tT. Therefore, |κ|max{C,K}K near T. Hence, there exists T[0,T) as mentioned in the statement. Considering tT, we may assume that |κ|K.

    In the case of (AP), since g=0,

    (R+g)K2L=RK2L.

    Using this and

    R=L0κdsL0|κ|dsLK,

    we have t(κ)2K3 on St, i.e.,

    d+Kdt(t)2K3.

    Consequently, we obtain the assertion as before.

    In the case of (LP),

    K2gL=K2I0RL=K2RL0˜κ2dsK2RL0κ2dsLK4R.

    The estimate RLK holds for all cases. Hence,

    K3=LRRLK3LK4R,K2RL=(RL)2LK2RLK4R.

    Consequently, we have

    d+Kdt(t)3LK4R=3L(0)K42πn.

    Here, we use LL(0). The statement follows from the above, as shown before.

    In the case of (JP), using R+g=L22A and Lemma 2.1, we have

    K2(R+g)L=K2L2A=L22ARLK2RL(T)22A(T)K3R=L(T)2K34πnA(T).

    Hence, it holds that

    d+Kdt(t)(1+L(T)24πnA(T))K3,

    which leads to the required conclusion and ends the proof.

    Remark 3.2. At a glance, the power 13 of blow-up rate in (LP) seems to be curious. The difference with other cases is that there is the length L(0) in the braces. If an estimate

    minsR/L(t)Zκ(s,t){2πn9L(0)(Tt)}p

    holds, then the power p must be 13. To see this, assume that f is a solution of (LP) which blows up at T<. For a positive constant λ, set

    fλ(s,t)=λ1f(λs,λ2t).

    We denote quantities of fλ the notation with the suffix λ; for example κλ is its curvature. Then, fλ satisfies (LP) with the length Lλ=λ1L(0), and blows up at Tλ=λ2T. The minimum of curvature is

    minsR/Lλ(t)Zκλ(s,t)=λminλsR/L(λ2t)Zκ(λs,λ2t)λ{2πn9L(0)(Tλ2t)}p.

    Using L(0)=λLλ(0) and T=λ2Tλ, we have

    λ{2πn9L(0)(Tλ2t)}p=λ13p{2πn9Lλ(0)(Tλt)}p.

    Hence, p must be 13. The L(0) in braces comes from the estimate K2gLLK4R in the proof. If we can improve this as K2gLCK3, then the blow-up rate coincides with other cases.

    In this section, we assume that f is a classical global solution of one of (AP), (LP), or (JP), and that the initial curve satisfies A(0)>0. We prove that Imf converges to an n-fold circle exponentially as t.

    Remark 4.1. However, this conclusion is meaningless if n-fold circles are only global solutions. At least, in the case of (AP), under suitable assumptions on the initial curve, regarding symmetry and convexity, solutions exist globally in time even if n>1. See [9].

    Firstly we prove the decay of I1.

    Lemma 4.1. For the global solution above, I1(t) fulfills

    0I1(t)L(0)2I1(0)L(t)2exp(t08π2nL(τ)2dτ).

    In particular, the estimate

    0I1(t)L(0)2I1(0)4πnA(0)exp(8π2nL(0)2t)

    is satisfied with respect to the global solution for (AP); the estimate

    0I1(t)I1(0)exp(8π2nL(0)2t)

    for the global solution of (LP). In the case of (JP), setting ˉL=supt[0,)L(t), we have ˉL<, and

    0I1(t)L(0)2I1(0)4πnA(0)exp(8π2nˉL2t).

    Proof. For global solutions, we know, from Theorem 3.1, that I1(t)0. Hence, we have

    4π2nI1(t)I0(t) (4.1)

    by Lemma 2.5. Consequently, (2.12) becomes

    ddt(L2I1)+8π2nL2(L2I1)0.

    Solving this differential inequality, we obtain the first assertion.

    We use 4nπA(0)L(t)L(0) for (AP), and L(t)L(0) for (LP). Then, the second assertion follows for these two cases.

    Now, we consider the case of (JP). Integrating (2.12), we have

    L2I1+2t0I0dτ=L20I1(0).

    L2A is uniformly positive and bounded by Proposition 2.1. From this, (2.10) with g=L2I12A and (4.1), we have

    dL2dt+2I0=2πnL2AI1L22πAI0CI0.

    Integrating this, we have

    L2+2t0I0(τ)dτL20+Ct0I0(τ)dτL20(1+CI1(0)).

    Hence, ˉL<. It follows from (2.11) and g=L2I12A0 that

    dA2dt=L2I10.

    Therefore, the lower bound L follows from L(t)4(4πnA(t))2(4πnA(0))2. Consequently, we obtain the second assertion for (JP).

    We denote the relevant statement of Lemma 4.1 as

    I1(t)Ceλ1t.

    Corollary 4.1. For the global solution above, there exists L>0 and A>0 such that

    |LL|+|AA|Ceλ1t.

    Proof. In the case of (AP), by Proposition 2.1, we have dLdt0. Hence, we conclude the convergence of limtL(t). Set the limit value as L. Since A(t)A(0), and since limtI1(t)=0, it holds that

    L2=limt4πnA(t)=4πnA(t)=4πnA(0)>0

    and LLL(0). Therefore,

    0 LL=L2L2L+L=L24πnAL+L=L2I1L+L L(0)2I12L=L(0)2I14πnA(0)Ceλ1t.

    In the case of (LP), since dAdt0 and since 4πAL2=L(0)2, we conclude the convergence of limtA(t). Set the limit value as A. Since L(t)L(0), and limtI1(t)=0, it holds that 4πnA=L(0)2. Consequently, (2.11) with g=I02πn yields

    0AA=tI02πndt=L204πnI1(t)Ceλ1t.

    Here, we use (2.12) and Lemma 4.1.

    In the case of (JP), dAdt=L2I12A0. By Proposition 2.1, AL2 is uniformly positive and bounded. Combining the above two statements with Lemma 4.1, we conclude

    0AA=tL2I12AdtCtI1dtCeλ1t.

    Furthermore, we estimate that

    |LL|= |L2L2|L+L=|L2I1+4πnA4πnA|L+L L2I1+4πn|AA|LCeλ1t.

    Corollary 4.2. For the global solution above, it holds that

    tI0dtCeλ1t.

    Proof. We know that L is uniformly bounded for all cases. Therefore, (2.12) implies that

    tI0dt=L2I12Ceλ1t.

    Lemma 4.2. For the global solution above, there exists λ0>0 such that

    I0Ceλ0t.

    Proof. As in Section 3, we set

    W=L0κ2ds,R=L0κds,Jp=Lp1L0˜κpds.

    As we know that LL>0 as t, it is enough to show that

    L2I0Ceλ0t.

    Since I0=J2=LWR2, we have from (2.10) and Lemma 3.2

    ddt(L2I0)= ddt(L3WR2L2)=L3dWdt+(32LWR2)dL2dt= 2I1+J4+(3Rg)J3+3R(Rg)J2R3g+(32I0+12R2)(2I0+2Rg)= 2I13I20+J4+(3Rg)J3+2R2J2.

    We obtain

    ddt(L2I0)+I1+3I20C(I30+I0+I530+|g|43I530)

    in a manner similar to the proof of Lemma 3.3.

    Since g=0 in (AP), and g=L2I12A in (JP), |g| is uniformly bounded for these cases. In (LP), g=R1I0. Hence, it holds for every case that

    ddt(L2I0)+I1+3I20C(I0+I30).

    This can be presented as

    ddt(L2I0)+I1+I20(3CI0)CI0.

    By Corollary 4.2, there exists t0>0 such that

    I0(t0)1C,t0I0dtL2C.

    Set

    t1=sup{t[t0,)|I0(t)<3C(t[t0,))}.

    If t1<, then,

    lim suptt10I0(t)=3C<.

    For t(t0,t1), we have

    ddt(L2I0)CI0,

    and therefore,

    I0(t)I0(t0)+1L2tt0I0dt2C=23lim suptt1I0(t).

    Letting tt1, we obtain a contradiction. Consequently, t1=, that is, I0(t)<3C for t[t0,). Since we know that I0 is uniformly bounded, we obtain

    ddt(L2I0)+I1+3I20CI0.

    It follows from Wirtinger's inequality and the uniform estimate of L2 that

    ddt(L2I0)+2λL2I0CI0

    for some constant λ>0. Multiplying both sides by e2λt, and integrating from t2 to t, we have

    e2λtL(t)2I0(t) CeλtL(t2)2I0(t2)+Ctt2e2λτI0(τ)dτ Ceλt+Ce2λtt2I0(τ)dτ.

    That is, we have

    L(t)2I0(t)Ceλt+Ct2I0(τ)dτ.

    Using the uniform estimate of L and the exponential decay of t2I0dt, we finally obtain the exponential decay of I0.

    Once we obtain the exponential decay of ˜I1 and I0, we can obtain the convergence of Imf to an n-fold circle as t.

    Theorem 4.1. Let f be a classical global solution of one of (AP), (LP), or (JP), with the smooth initial curve satisfying A(0)>0. Then, Imf converges to an n-fold circle with centre c, and radius r=L2πn in the following sense. Set

    f(s,t)=c(t)+r(t)(cos2πn(s+σ(t))L(t),sin2πn(s+σ(t))L(t))+ρ(s,t),c(t)=1L(t)L(t)0f(s,t)ds,r(t)=L(t)2πn,

    with the R/L(t)Z-valued function σ defined by

    ˆf(n)(t)=L(t)r(t)exp(2πinσ(t)L(t)).

    Then, there exist cR2, r=L2πn>0, σR/LZ, λ>0, and C>0 such that

    c(t)c+|r(t)r|+|σ(t)L(t)σL|Ceλt.

    Furthermore, for k{0}N, there exist γk>0 and Ck>0 such that

    ρ(,t)Ck(R/L(t)Z)Ckeγkt.

    When n=1, we used (1.1) for the proof of this theorem in [7,§ 4], and [8,§ 2.2]. The most crucial part is to show the decay of I0. As above, we have already obtained a decay estimate of I0 without using (1.1) for n1. Once we obtain it, to show the theorem, we can perform the standard energy method with help of usual Gagliardo-Nirenberg's inequality rather than (1.1) as the previous papers. In this sense, (1.1) is not absolutely necessary, however, we need several modification of argument. Using (1.2) which is an alternative inequality to (1.1), we can develop the argument almost word to word as the previous papers. Thus, we deal with (1.2) in the next section.

    We discuss (1.2) in this section. Set

    ˜I1=4π2L3kZ{0}(kn)2|ˆf(k)|2.

    Proposition 5.1. We have

    ˜I1=1L2πnL(f1LL0fds)+ν2L2.

    ˜I1 vanishes if and only if Imf is an n-fold circle.

    Proof. Setting

    ˜f=f1LL0fds,

    we have

    ˜f2L2=kZ{0}|ˆf(k)|2.

    The squared L2-norm of ν is

    ν2L2=f2L2=kZ(2πkL)2|ˆf(k)|2=4π2L2kZ{0}k2|ˆf(k)|2.

    On the other hand, we have

    ˜f,νL2=˜f,if=kZ{0}2πkL|ˆf(k)|2=4π2L2kZ{0}kL2π|ˆf(k)|2.

    Since the last right-hand side expression is a real number, it holds that

    4π2L2kZ{0}k|ˆf(k)|2=2πL˜f,νL2.

    Consequently, we obtain

    4π2L2kZ{0}(kn)2|ˆf(k)|2= ν2L2+4nπL˜f,fL2+(2πnL)2˜f2L2= 2πnL˜f+ν2L2= 2πnL(f1LL0fds)+ν2L2.

    ˜I1 vanishes if and only if

    f=ˆf(0)φ0+ˆf(n)φn.

    Hence, Imf is an n-fold circle.

    An estimate similar to Lemma 2.5 holds for ˜I1 as well.

    Lemma 5.1. It holds that 4π2˜I1I0.

    Proof. Since k210 for kZ{0}, we have

    I04π2˜I1=16π4L3kZ{0}(k21)(kn)2|ˆf(k)|20.

    The next proposition corresponds to [7,Theorem 2.2].

    Proposition 5.2. It holds that

    I0˜I121[L0L3{κ4+(κ)2}ds.]

    Proof. It follows from Lemma 2.2, Schwarz' inequality, and (2.6) that

    I0= 16π4L3kZ{0}k3(kn)|ˆf(k)|2 8π3L32{4π2L3kZ{0}(kn)2|ˆf(k)|2}12{kZ{0}k6|ˆf(k)|2}12= 8π3L32˜I121{kZ{0}k6|ˆf(k)|2}12= ˜I121[L0L3{κ4+(κ)2}ds].

    Using this proposition, we can prove the following estimates.

    Theorem 5.1. Let j[0,] be an integer. Then, there exists a positive constant C=C(j,) independent of L such that

    IjC(˜Ij21I+˜Ij+11Ij+1+1).

    Proof. Since the assertion can be proven in a manner similar to the proof of [7,Theorem 3.1], we give only the sketch. Firstly, we derive

    I0C˜I121(I1+˜I121) (5.1)

    from Proposition 5.2 and Gagliardo-Nirenberg's inequality

    (L(j+1)p1L0|˜κ(j)|pds)1pC(j,m,p)I12m(j1p+12)mI12{11m(j1p+12)}0 (5.2)

    for p2 and jm. Here C(j,m,p) is independent of L. It follows from (5.2) that

    IjC(j,n)IjmnI1jm0. (5.3)

    Combining this together with (5.1), we obtain the assertion for j=0. It gives also the assertion for j1 with help of (5.3).

    For the proof of convergence of global flow to a circle, we use in [7] the following properties of I1:

    (ⅰ) I10,

    (ⅱ) I1=0 holds if and if the image of f is a circle,

    (ⅲ) C1I1I0 (an inequality of Wirtinger's type).

    These are satisfied when n=1, but not when n>1. The quantity ˜I1 satisfies

    (ⅰ) ˜I10,

    (ⅱ) ˜I1=0 holds if and if the image of f is an n-fold circle,

    (ⅲ) C1˜I1I0 (an inequality of Wirtinger's type).

    Hence, it is an alternative quantity to I1.

    The first author is partly supported by Grant-in-Aid for Scientific Research (C) (17K05310), and (B) (20H01813), Japan Society for the Promotion Science. The authors express their appreciation to the anonymous referee for his/her suggestive comments and information of related articles [1,10].

    The authors declare no conflict of interest.



    [1] B. Andrews, J. McCoy, G. Wheeler, V. M. Wheeler, Closed ideal planar curves, arXiv: 1810.06154.
    [2] B. Chow, P. Lu, L. Ni, Hamilton's Ricci flow, New York: American Mathematical Society, 2006.
    [3] M. Gage, On an area-preserving evolution equation for plane curves, In: Nonlinear problems in geometry, Providence: American Mathematical Society, 1986, 51–62.
    [4] J. W. Hagood, B. S. Thomson, Recovering a function from a Dini derivative, Am. Math. Mon., 113 (2006), 34–46. doi: 10.1080/00029890.2006.11920276
    [5] L. Jiang, S. Pan, On a non-local curve evolution problem in the plane, Commun. Anal. Geom., 16 (2008), 1–26. doi: 10.4310/CAG.2008.v16.n1.a1
    [6] L. Ma, A. Zhu, On a length preserving curve flow, Monatsh. Math., 165 (2012), 57–78. doi: 10.1007/s00605-011-0302-8
    [7] T. Nagasawa, K. Nakamura, Interpolation inequalities between the deviation of curvature and the isoperimetric ratio with applications to geometric flows, Adv. Differential Equ., 24 (2019), 581–608.
    [8] K. Nakamura, An application of interpolation inequalities between the deviation of curvature and the isoperimetric ratio to the length-preserving flow, Discrete Contin. Dyn. S, doi: 10.3934/dcdss.2020385.
    [9] X. L. Wang, L. H. Kong, Area-preserving evolution of nonsimple symmetric plane curves, J. Evol. Equ., 14 (2014), 387–401. doi: 10.1007/s00028-014-0219-5
    [10] G. Wheeler, Convergence for global curve diffusion flows, arXiv: 2004.08494.
  • This article has been cited by:

    1. Xiao-Liu Wang, The evolution of area-preserving and length-preserving inverse curvature flows for immersed locally convex closed plane curves, 2023, 284, 00221236, 109744, 10.1016/j.jfa.2022.109744
    2. Takeyuki Nagasawa, Lower estimate of blow-up time and blow-up rate for non-local curvature flows for plane curves, 2025, 0916-7005, 10.1007/s13160-025-00690-w
  • Reader Comments
  • © 2021 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(2898) PDF downloads(624) Cited by(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog