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
[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 |
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), s∈R/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
κ=κ⋅ν,˜κ=κ−1L∫L0κ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) (s∈R/λ−1LZ), then,
g(fλ)=g(f). |
Here we study three cases of g:
(AP) If we set g≡0, then our equation represents the area-preserving flow. In fact, we set A as
A=−12∫L0f⋅νds |
which is the enclosed area when Imf is a simple curve. Consequently, it holds that
dAdt=0. |
(LP) Let g=L(∫L0κds)−1∫L0˜κ2ds. Here, the equation represents the length-preserving flow:
dLdt=0. |
(JP) Jiang-Pan considered an equation with g=L22A−∫L0κds in [5]. Here, the isoperimetric ratio does not increase along with the flow:
ddtL2A=−2LA∫L0‖∂tf‖2ds. |
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
I−1=1−4πAL2 |
played an important role. First, we show the decay of I−1. Set
Iℓ=L2ℓ+1∫L0|˜κ(ℓ)|2ds for ℓ∈{0}∪N. |
In [7], we showed the inequality
Ij≦C(Iℓ−j2−1Iℓ+Iℓ−jℓ+1−1Ij+1ℓ+1ℓ) | (1.1) |
for an integer j∈[0,ℓ] with a positive constant C=C(j,ℓ) independent of the total length of curve. Since I−1 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
I−1=1−4nπAL2. |
The isoperimetric inequality shows I−1≧0 when n=1. However, I−1 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, I−1 gives us some useful information. For example, we can show that if I−1 is negative for t=0, then the solution blows up in finite time. See our first main result, Theorem 3.1. This implies I−1≧0 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:
˜I−1=1L‖2πnL(f−1L∫L0fds)+ν‖2L2. |
Then we can show
Ij≦C(˜Iℓ−j2−1Iℓ+˜Iℓ−jℓ+1−1Ij+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 n≧1. Set
f=f1+if2,ν=ν1+iν2=−f′2+if′1=if′. |
The functions φk=1√Lexp(2πiksL) for k∈Z 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.
∑k∈Zk|ˆf(k)|2= LAπ, | (2.1) |
∑k∈Zk2|ˆf(k)|2= (L2π)2∫L0κ0ds=L34π2, | (2.2) |
∑k∈Zk3|ˆf(k)|2= (L2π)3∫L0κds=nL34π2, | (2.3) |
∑k∈Zk4|ˆf(k)|2= (L2π)4∫L0κ2ds, | (2.4) |
∑k∈Zk5|ˆf(k)|2= (L2π)5∫L0κ3ds, | (2.5) |
∑k∈Zk6|ˆf(k)|2= (L2π)6∫L0{κ4+(κ′)2}ds. | (2.6) |
Note that we have
∑k∈Zk2(k−n)|ˆ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π4L3∑k∈Zk3(k−n)|ˆf(k)|2 | (2.8) |
= 16π4L3∑k∈Zk2(k−n)2|ˆf(k)|2. | (2.9) |
Proof. We obtain (2.8) as
I0= L∫L0˜κ2ds=L∫L0˜κκds=L(∫L0κ2ds−2πnL∫L0κds)= 16π4L3∑k∈Zk3(k−n)|ˆ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 L2−4πA≧0.
Proof. It follows from (2.2) and (2.1) that
L2−4πA=4π2L(L34π2−LAπ)=4π2L∑k∈Zk(k−1)|ˆf(k)|2≧0. |
Similarly, I−1 has two expressions.
Lemma 2.4. We have
I−1=4π2L3∑k∈Zk(k−n)|ˆf(k)|2=−4π2nL3∑k∈Z∖{0}k(k−n)2|ˆf(k)|2. |
Proof. It follows from (2.2) and (2.1) that
I−1=1−4πnAL2=4π2L3(L34π2−nLAπ)=4π2L3∑k∈Zk(k−n)|ˆf(k)|2. |
The second expression of I−1 is obtained from the above and (2.7).
Since k(k−n) is not necessarily non-negative when n>1, we know the same holds for I−1. However, the modulus of I−1 can be estimated by I0 for n≧1 as follows. This is Wirtinger's inequality when n=1.
Lemma 2.5. It holds that 4π2n|I−1|≦I0.
Proof. From Lemmas 2.2–2.4 we obtain
I0±4π2nI−1= 16π4L3∑k∈Z{k2(k−n)2∓k(k−n)2}|ˆf(k)|2= 16π4L3∑k∈Zk(k∓1)(k−n)2|ˆf(k)|2≧0. |
Here, we use k(k∓1)≧0 for k∈Z.
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=−∫L0∂tf⋅κds, we have
dL2dt=−2L∫L0(˜κ−gL)κds=−2L∫L0˜κ2ds+4πng, |
that is,
dL2dt+2I0=4πng. | (2.10) |
Similarly, we have
dAdt=−∫L0∂tf⋅νds=−∫L0(˜κ−gL)ds=g. | (2.11) |
It follows from the above that
ddt(L2I−1)+2I0=ddt(L2−4π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,A≡A(0)>0,dL2dt≦0,dI−1dt≦0. |
2). For solutions of (LP),
dAdt≧0,A≧A(0)>0,dL2dt=0,dI−1dt≦0. |
3). For solutions of (JP),
A>0,dI−1dt≦0. |
4). For solutions of (AP), (LP), (JP),
1−n≦I−1≦I−1(0). |
In other words,
4π≦L2A≦L(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
dI−1dt=−ddt4πnAL2=−4πnL2dAdt+4πnAL4dL2dt≦0. |
In the case of (JP), we prove the positivity of A by applying the contradiction argument. In this case,
g=L2I−12A. | (2.13) |
It follows from (2.11) that
dA2dt=2Ag=L2I−1. | (2.14) |
Assume that A(t0)2=0 for some first time t0∈(0,T). Since A2≧0, 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(L2I−1)=−2I0≦0. |
Therefore, by (2.16)
dA2dt(t0)=dA2dt(t1)+∫t0t1d2A2dt2dt≦dA2dt(t1)<0. |
This contradicts (2.15). Hence, A>0 on (0,T). Using (2.12), (2.10), I−1−1=−4πnAL2, and (2.13), we have
L2dI−1dt= −I−1dL2dt−2I0= −I−1(4πng−2I0)−2I0=−4πngI−1+2(I−1−1)I0= −4πnL2(L2gI−1+2AI0)=−4πnL2(L4I2−12A+2AI0)≦0. |
Since I−1 is non-increasing, we have I−1≦I−1(0). Lemma 2.3 gives us
I−1=1−4πnAL2=1−n+n(1−4πAL2)≧1−n. |
The non-positivity of I−1(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, I−1(0)<0. Then, the solution blows up in finite time. The blow-up time T is estimated from above as follows:
(AP)T≦L(0)2−4πA(0)−8π2nI−1(0),
(LP)T≦L(0)2−4πA(0)−8π2I−1(0),
(JP)T≦L(0)2−8π2nI−1(0).
Proof. In the case of (AP), g≡0. It follows from Proposition 2.1 that I−1(t)≦I−1(0)<0. By (2.10) and Lemma 2.5, we have
dL2dt=−2I0(t)≦8π2nI−1(t)≦8π2nI−1(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π4nI−1(0)t. |
Since the first side is non-positive by the isoperimetric inequality (Lemma 2.3), t must satisfy
t≦L(0)2−4πA(0)−8π2nI−1(0). |
In the case of (LP), g=I02πn≧0. Proposition 2.1 shows I−1(t)≦I−1(0)<0. From (2.11) and Lemma 2.5, we have
−dAdt=−12πnI0(t)≦2πI−1(t)≦2πI−1(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)2≦4π(A(0)−A(t))≦8πI−1(0)t. |
Consequently, t must satisfy
t≦L(0)2−4πA(0)−8π2I−1(0). |
In the case of (JP), g=L2I−12A. It follows from (2.10), Proposition 2.1, and Lemma 2.5 that
dL2dt=−2I0(t)+2πnL(t)2A(t)I−1(t)≦−2I0(t)≦8π2nI−1(t)≦8π2nI−1(0). |
We integrate this from 0 to t∈(0,T). Using Lemma 2.3, we obtain
−L(0)2≦L(t)2−L(0)2≦8π2nI−1(0)t. |
Consequently t must satisfy
t≦L(0)2−8π2nI−1(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 I−1(0)=0, but it is not an n-fold circle. Then, T<∞.
Proof. Assume T=∞. Then, Theorem 3.1 implies that I−1(t)≧0 for all t∈[0,∞). On the other hand, (2.12) with I−1(0)=0 shows that I−1(t)≦0. Hence, I−1(t)≡0. When t>0,
∫L0˜κ2ds=I0L=−12Lddt(L2I−1)=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=I−1=0 by Lemmas 2.2 and 2.4. Hence, ˜κ−gL≡0 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 t↗T<∞. Then, we have
lim supt↗TI0(t)=∞. |
Indeed, if lim supt↗TI0(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 t↗T. 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 supt↗TW(t)=∞.
Proof. Set
R=∫L0κds, |
and we have
LW=L∫L0{˜κ2+(RL)2}ds=I0+R2. |
Hence,
lim supt↗TL(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 L≧√4πA=√4πA0. Consequently, L(t) converges to a positive constant as t↗T, 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 inft↗TA(t)=0. I−1 is monotone by Proposition 2.1. Therefore, it follows from
dAdt=L22AI−1 |
that A does not oscillate near t=T. Hence, we may assume limt↗TA(t)=0. From the above relation and Proposition 2.1, we find that dAdt is bounded. Consequently, the estimate
0<A(t)≦C(T−t) |
holds. Thus, we have
0≦A(t)2T−t≦C(T−t)2T−t→0 as t↗T, |
and therefore,
limt↗TA(T−0)2−A(t)2T−t=0. |
This implies that the left derivative of A2 at T vanishes:
dA2dt(T−0)=0. | (3.1) |
However, A(0)2>0 and A(T−0)2=0 show the existence of t∗∈(0,T) such that
dA2dt(t∗)<0. |
Since
d2A2dt2=−2I0≦0, |
we have
dA2dt(t)<dA2dt(t∗)<0 |
for t∈(t∗,T). This contradicts (3.1). Now, we prove lim inft↗TA(t)>0. Since
dAdt=L22AI−1 |
has a constant sign near T, we conclude that limt↗TA(t)>0. The convergence of limt↗TL(t) follows from the convergence of A, and the monotonicity and boundedness of I−1. 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=Lp−1∫L0˜κpds(p∈N∖{1}), |
which are scale-invariant quantities. Note that I0=J2.
Lemma 3.2. It holds that
dWdt=1L3{−2I1+J4+(3R−g)J3+3R(R−g)J2−R3g}. |
Proof. The proof is a direct calculation:
dWdt= ∫L0∂tf⋅(2∂2sκ+κ3)ds=∫L0(˜κ−gL)(2∂2sκ+κ3)ds= −2∫L0(∂s˜κ)2ds+∫L0(˜κ−gL)(˜κ+RL)3ds= −2I1L3+∫L0(˜κ3+3R˜κ2L+3R2˜κL2+R3L3)(˜κ−gL)ds= −2I1L3+∫L0{˜κ4+(3RL−gL)˜κ3+(3R2L2−3RgL2)˜κ2−R3gL4}ds= 1L3{−2I1+J4+(3R−g)J3+3R(R−g)J2−R3g}. |
Thirdly, we estimate dWdt from above.
Lemma 3.3. We have
dWdt≦W32M2. |
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 θ=12−1p. Then, Gagliardo-Nirenberg's inequality yields
|Jp|≦C(I1−θ0Iθ1)p2=CIp4+120Ip4−121. |
Hence,
dWdt+2I1L3≦ CL3(I320I121+RI540I141+R2I0)≦ I1L3+CL3(I30+R43I530+R2I0). |
Since 0≦I0≦LW and R2≦LW, we obtain
I30≦L3W3,I530≦L53W53=(LW)−43L3W3≦R−83L3W3,I0≦LW=(LW)−4L3W3≦R−8L3W3. |
Furthermore,
R=2πn≧2π. |
Consequently, we conclude that
dWdt≦C(1+R−43+R−6)W3≦CW3. |
In the case of (LP), since g=I0R, we have
dWdt+1L3(2I1+3I20+R2I0)= 1L3{J4+(3R−I0R)J3+3R2I0}≦ CL3(I320I121+RI540I141+R−1I940I141+R2I0)≦ I1L3+CL3(I30+R43I530+R−43I30+R2I0)≦ I1L3+C(1+R−43+R−6)W3≦ I1L3+CW3. |
In the case of (JP), since g=L22A−R, we have
dWdt+1L3(2I1+3RL22AI0+R3L22A)=1L3{J4+(3R−L22A+R)J3+6R2J2+R4}≦CL3{I320I121+RI540I141+L2AI540I141+R2I0+R−2(LW)3}≦I1L3+CL3[I30+{R+(L2A)}43I530+R2I0+R−2L3W3]≦I1L3+C[1+{1R+(L2R2A)}43+R−6+R−2]W3≦I1L3+C{1+(L2A)43}W3. |
By Proposition 2.1, we have
(L2A)43≦(L20A0)43. |
Consequently, we can conclude that
dWdt≦C{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)≧M√T−t, |
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
ddtW−2≧−M−2. |
Due to Lemma 3.1, there exists a sequence {tn} such that tnt↗T and W(tn)−2→0 as n→∞. Integrating the differential inequality from t to tn, we have
W(t)−2−W(tn)−2≦M−2(tn−t). |
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
limt↗Tmaxs∈R/L(t)Zκ(s,t)=∞ |
or
limt↗Tmins∈R/L(t)Zκ(s,t)=−∞. |
If a positive/negative loop of Imf shrinks as t↗T, 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 t↗T. For example, if a negative loop shrinks as t↗T 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 supt↗Tmaxs∈R/L(t)Zκ(s,t)=∞, |
then it satisfies
maxs∈R/L(t)Zκ(s,t)≧1√2(T−t). |
Proof. Set
K(t)= maxs∈R/L(t)Zκ(s,t),d+Kdt(t)= lim suph→+0K(t+h)−K(t)h. |
Define the set St by St={s∈R/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)=maxs∈St∂tκ(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 g≧0. In the case of (JP),
R+g=L2A≧0. |
∂2sκ≦0 holds for s∈St. Hence, we have
∂2sκ+κ2(κ−R+gL)≦κ3=K3 |
for s∈St, and
d+Kdt(t)≦maxs∈St∂tκ≦K3(t). |
We calculate Dini's derivative of K−2 as
d+dtK−2(t)= lim suph→+0K−2(t+h)−K−2(t)h= lim suph→+0(K(t)+K(t+h))(K(t)−K(t+h))K2(t+h)K2(t)h= −2K−3(t)lim infh→+0K(t+h)−K(t)h≧ −2K−3(t)lim suph→+0K(t+h)−K(t)h= −2K−3(t)d+Kdt(t)≧−2. |
According to the assumption of the theorem, there exists a sequence {tk}k∈N such that tk↗T and K(tk)−2→0 as k→∞. Using [4,Theorem 3], we have
K−2(tk)−K−2(t)≧∫_ tktd+dtK−2(t)dt≧−2(tk−t) |
for tk∈(t,T). Therefore, we can conclude that
K−2(t)≦2(T−t) |
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)maxs∈R/L(t)Zκ(s,t)<∞. |
For the solution of (AP),
mins∈R/L(t)Zκ(s,t)≦−1√4(T−t) |
holds.
For the solution of (LP),
mins∈R/L(t)Zκ(s,t)≦−{2πn9L(0)(T−t)}13 |
holds.
For the solution of (JP), there exists a time T∗∈[0,T) such that
−mins∈R/L(t)Zκ(s,t)≧maxs∈R/L(t)Zκ(s,t) |
holds for t∈[T∗,T). Additionally, it holds that
mins∈R/L(t)Zκ(s,t)≦−1√2C∗(T−t), |
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)= −mins∈R/L(t)Zκ(s,t),d+Kdt(t)= lim suph→+0K(t+h)−K(t)h. |
Define the set St by St={s∈R/L(t)Z|−κ(s,t)=K(t)}. As shown before, it holds that
d+Kdt(t)=maxs∈St∂t(−κ). |
−κ satisfies
∂t(−κ)=∂2s(−κ)+(−κ)2{(−κ)+R+gL}. |
Since ∂2s(−κ)≦0 and −κ=K for s∈St,
∂t(−κ)≦K3+(R+g)K2L. |
If κ≦C<∞ holds on [0,T), then,
Lmax{C2+K2}≧∫L0κ2ds=W→∞ as t↗T |
by Theorem 3.2. Since L is bounded, we conclude that K→∞ as t↗T. Therefore, |κ|≦max{C,K}≦K near T. Hence, there exists T∗∈[0,T) as mentioned in the statement. Considering t≧T∗, we may assume that |κ|≦K.
In the case of (AP), since g=0,
(R+g)K2L=RK2L. |
Using this and
R=∫L0κds≦∫L0|κ|ds≦LK, |
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=K2R∫L0˜κ2ds≦K2R∫L0κ2ds≦LK4R. |
The estimate RL≦K holds for all cases. Hence,
K3=LR⋅RL⋅K3≦LK4R,K2RL=(RL)2LK2R≦LK4R. |
Consequently, we have
d+Kdt(t)≦3LK4R=3L(0)K42πn. |
Here, we use L≡L(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=L22A⋅RL⋅K2R≦L(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
mins∈R/L(t)Zκ(s,t)≦−{2πn9L(0)(T−t)}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
mins∈R/Lλ(t)Zκλ(s,t)=λminλs∈R/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=−λ1−3p{2πn9Lλ(0)(Tλ−t)}p. |
Hence, p must be 13. The L(0) in braces comes from the estimate K2gL≦LK4R in the proof. If we can improve this as K2gL≦CK3, 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 I−1.
Lemma 4.1. For the global solution above, I−1(t) fulfills
0≦I−1(t)≦L(0)2I−1(0)L(t)2exp(−∫t08π2nL(τ)2dτ). |
In particular, the estimate
0≦I−1(t)≦L(0)2I−1(0)4πnA(0)exp(−8π2nL(0)2t) |
is satisfied with respect to the global solution for (AP); the estimate
0≦I−1(t)≦I−1(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
0≦I−1(t)≦L(0)2I−1(0)4πnA(0)exp(−8π2nˉL2t). |
Proof. For global solutions, we know, from Theorem 3.1, that I−1(t)≧0. Hence, we have
4π2nI−1(t)≦I0(t) | (4.1) |
by Lemma 2.5. Consequently, (2.12) becomes
ddt(L2I−1)+8π2nL2(L2I−1)≦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
L2I−1+2∫t0I0dτ=L20I−1(0). |
L2A is uniformly positive and bounded by Proposition 2.1. From this, (2.10) with g=L2I−12A and (4.1), we have
dL2dt+2I0=2πnL2AI−1≦L22πAI0≦CI0. |
Integrating this, we have
L2+2∫t0I0(τ)dτ≦L20+C∫t0I0(τ)dτ≦L20(1+CI−1(0)). |
Hence, ˉL<∞. It follows from (2.11) and g=L2I−12A≧0 that
dA2dt=L2I−1≧0. |
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
I−1(t)≦Ce−λ−1t. |
Corollary 4.1. For the global solution above, there exists L∞>0 and A∞>0 such that
|L−L∞|+|A−A∞|≦Ce−λ−1t. |
Proof. In the case of (AP), by Proposition 2.1, we have dLdt≦0. Hence, we conclude the convergence of limt→∞L(t). Set the limit value as L∞. Since A(t)≡A(0), and since limt→∞I−1(t)=0, it holds that
L2∞=limt→∞4πnA(t)=4πnA(t)=4πnA(0)>0 |
and L∞≦L≦L(0). Therefore,
0≦ L−L∞=L2−L2∞L+L∞=L2−4πnAL+L∞=L2I−1L+L∞≦ L(0)2I−12L∞=L(0)2I−14√πnA(0)≦Ce−λ−1t. |
In the case of (LP), since dAdt≧0 and since 4πA≦L2=L(0)2, we conclude the convergence of limt→∞A(t). Set the limit value as A∞. Since L(t)≡L(0), and limt→∞I−1(t)=0, it holds that 4πnA∞=L(0)2. Consequently, (2.11) with g=I02πn yields
0≦A∞−A=∫∞tI02πndt=L204πnI−1(t)≦Ce−λ−1t. |
Here, we use (2.12) and Lemma 4.1.
In the case of (JP), dAdt=L2I−12A≧0. By Proposition 2.1, AL2 is uniformly positive and bounded. Combining the above two statements with Lemma 4.1, we conclude
0≦A∞−A=∫∞tL2I−12Adt≦C∫∞tI−1dt≦Ce−λ−1t. |
Furthermore, we estimate that
|L−L∞|= |L2−L2∞|L+L∞=|L2I−1+4πnA−4πnA∞|L+L∞≦ L2I−1+4πn|A−A∞|L∞≦Ce−λ−1t. |
Corollary 4.2. For the global solution above, it holds that
∫∞tI0dt≦Ce−λ−1t. |
Proof. We know that L is uniformly bounded for all cases. Therefore, (2.12) implies that
∫∞tI0dt=L2I−12≦Ce−λ−1t. |
Lemma 4.2. For the global solution above, there exists λ0>0 such that
I0≦Ce−λ0t. |
Proof. As in Section 3, we set
W=∫L0κ2ds,R=∫L0κds,Jp=Lp−1∫L0˜κpds. |
As we know that L→L∞>0 as t→∞, it is enough to show that
L2I0≦Ce−λ0t. |
Since I0=J2=LW−R2, we have from (2.10) and Lemma 3.2
ddt(L2I0)= ddt(L3W−R2L2)=L3dWdt+(32LW−R2)dL2dt= −2I1+J4+(3R−g)J3+3R(R−g)J2−R3g+(32I0+12R2)(−2I0+2Rg)= −2I1−3I20+J4+(3R−g)J3+2R2J2. |
We obtain
ddt(L2I0)+I1+3I20≦C(I30+I0+I530+|g|43I530) |
in a manner similar to the proof of Lemma 3.3.
Since g=0 in (AP), and g=L2I−12A in (JP), |g| is uniformly bounded for these cases. In (LP), g=R−1I0. Hence, it holds for every case that
ddt(L2I0)+I1+3I20≦C(I0+I30). |
This can be presented as
ddt(L2I0)+I1+I20(3−CI0)≦CI0. |
By Corollary 4.2, there exists t0>0 such that
I0(t0)≦1C,∫∞t0I0dt≦L2C. |
Set
t1=sup{t∈[t0,∞)|I0(t)<3C(t∈[t0,∞))}. |
If t1<∞, then,
lim supt→t1−0I0(t)=3C<∞. |
For t∈(t0,t1), we have
ddt(L2I0)≦CI0, |
and therefore,
I0(t)≦I0(t0)+1L2∫tt0I0dt≦2C=23lim supt→t1I0(t). |
Letting t↗t1, 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+3I20≦CI0. |
It follows from Wirtinger's inequality and the uniform estimate of L2 that
ddt(L2I0)+2λL2I0≦CI0 |
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)+C∫tt2e2λτI0(τ)dτ≦ Ceλt+Ce2λt∫∞t2I0(τ)dτ. |
That is, we have
L(t)2I0(t)≦Ce−λt+C∫∞t2I0(τ)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 ˜I−1 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∞=L∞2π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 c∞∈R2, r∞=L∞2πn>0, σ∞∈R/L∞Z, λ>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 n≧1. 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
˜I−1=4π2L3∑k∈Z∖{0}(k−n)2|ˆf(k)|2. |
Proposition 5.1. We have
˜I−1=1L‖2πnL(f−1L∫L0fds)+ν‖2L2. |
˜I−1 vanishes if and only if Imf is an n-fold circle.
Proof. Setting
˜f=f−1L∫L0fds, |
we have
‖˜f‖2L2=∑k∈Z∖{0}|ˆf(k)|2. |
The squared L2-norm of ν is
‖ν‖2L2=‖f′‖2L2=∑k∈Z(2πkL)2|ˆf(k)|2=4π2L2∑k∈Z∖{0}k2|ˆf(k)|2. |
On the other hand, we have
⟨˜f,ν⟩L2=⟨˜f,if′⟩=−∑k∈Z∖{0}2πkL|ˆf(k)|2=−4π2L2∑k∈Z∖{0}kL2π|ˆf(k)|2. |
Since the last right-hand side expression is a real number, it holds that
4π2L2∑k∈Z∖{0}k|ˆf(k)|2=−2πLℜ⟨˜f,ν⟩L2. |
Consequently, we obtain
4π2L2∑k∈Z∖{0}(k−n)2|ˆf(k)|2= ‖ν‖2L2+4nπLℜ⟨˜f,f′⟩L2+(2πnL)2‖˜f‖2L2= ‖2πnL˜f+ν‖2L2= ‖2πnL(f−1L∫L0fds)+ν‖2L2. |
˜I−1 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 ˜I−1 as well.
Lemma 5.1. It holds that 4π2˜I−1≦I0.
Proof. Since k2−1≧0 for k∈Z∖{0}, we have
I0−4π2˜I−1=16π4L3∑k∈Z∖{0}(k2−1)(k−n)2|ˆf(k)|2≧0. |
The next proposition corresponds to [7,Theorem 2.2].
Proposition 5.2. It holds that
I0≦˜I12−1[∫L0L3{κ4+(κ′)2}ds.] |
Proof. It follows from Lemma 2.2, Schwarz' inequality, and (2.6) that
I0= 16π4L3∑k∈Z∖{0}k3(k−n)|ˆf(k)|2≦ 8π3L32{4π2L3∑k∈Z∖{0}(k−n)2|ˆf(k)|2}12{∑k∈Z∖{0}k6|ˆf(k)|2}12= 8π3L32˜I12−1{∑k∈Z∖{0}k6|ˆf(k)|2}12= ˜I12−1[∫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
Ij≦C(˜Iℓ−j2−1Iℓ+˜Iℓ−jℓ+1−1Ij+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
I0≦C˜I12−1(I1+˜I121) | (5.1) |
from Proposition 5.2 and Gagliardo-Nirenberg's inequality
(L(j+1)p−1∫L0|˜κ(j)|pds)1p≦C(j,m,p)I12m(j−1p+12)mI12{1−1m(j−1p+12)}0 | (5.2) |
for p≧2 and j≦m. Here C(j,m,p) is independent of L. It follows from (5.2) that
Ij≦C(j,n)IjmnI1−jm0. | (5.3) |
Combining this together with (5.1), we obtain the assertion for j=0. It gives also the assertion for j≧1 with help of (5.3).
For the proof of convergence of global flow to a circle, we use in [7] the following properties of I−1:
(ⅰ) I−1≧0,
(ⅱ) I−1=0 holds if and if the image of f is a circle,
(ⅲ) C−1I−1≦I0 (an inequality of Wirtinger's type).
These are satisfied when n=1, but not when n>1. The quantity ˜I−1 satisfies
(ⅰ) ˜I−1≧0,
(ⅱ) ˜I−1=0 holds if and if the image of f is an n-fold circle,
(ⅲ) C−1˜I−1≦I0 (an inequality of Wirtinger's type).
Hence, it is an alternative quantity to I−1.
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. |
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 |