On the rate of convergence of the $p$-curve shortening flow
Jean C. Cortissoz, Andr\'es Galindo, Alexander Murcia

TL;DR
This paper establishes improved and likely sharp rates of convergence for the p-curve shortening flow when p is an integer greater than or equal to 1, advancing understanding of the flow's behavior.
Contribution
It provides new, sharper convergence rate estimates for the p-curve shortening flow, enhancing previous results and likely representing the optimal rates.
Findings
Improved convergence rate estimates for p-curve shortening flow.
Results are probably sharp, indicating optimality.
Advances understanding of flow dynamics for p ≥ 1.
Abstract
In this paper we give rates of convergence for the -curve shortening flow for an integer, which improves on the known estimates and which are probably sharp.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsGeometric Analysis and Curvature Flows · Geometry and complex manifolds · Point processes and geometric inequalities
On the rate of convergence of the -curve shortening flow
Jean C. Cortissoz, Andrés Galindo and Alexander Murcia
Department of Mathematics, Universidad de los Andes, Bogotá DC, Colombia.
Abstract.
In this paper we give rates of convergence for the -curve shortening flow for an integer, which improves on the known estimates and which are probably sharp.
1. Introduction
Let us first introduce the main character of this story, the -curve shortening flow, with a positive integer. So, we let
[TABLE]
be a family of smooth convex embeddings of , the unit circle, into . We say that satisfies the -curve shortening flow, , if satisfies
[TABLE]
where is the curvature of the embedding and is the normal vector pointing outwards the region bounded by .
This is just a natural generalisation of the well known and well studied curve shortening flow. A solution to 1 starting from an embedded convex simple curve will contract, via embedded convex curves, towards a round point in finite time: this means that if we start with a simple convex curve, via the -curve shortening, after a convenient normalisation, which includes a time reparametrisation, the embedded curves converge smoothly to a circle (see [1]). It is also known that this convergence is exponential in the following sense (here denotes the curvature of the embedded curves after normalisation)
[TABLE]
with where represents the -derivative of with respect to the arclength parameter in , , and . For the curve shortening flow, we can use as any for , this was proved by Gage and Hamilton in their by now famous (by mathematical standards) paper [6]. For the -curve shortening, Huang in [7] showed that can be taken as , with the same restrictions on . Interestingly enough, with the exception of the curve shortening flow (), it has not been showed that exponentially! For the curve shortening flow (), in the book [3] exponential convergence of the curvature towards 1 is shown, and Andrews and Bryan showed in [2] (although they did not stated explicitly) that as fast as .
Related to this problem is the mean curvature flow, and Sesum in [9], using Huisken’s work as a departure point, has given sharp rates of convergence for this flow.
The main goal of this paper to give better rates of convergence for the -curve shortening flow, an integer, than the ones previously known. Our main result, from which the said rates of convergence can be deduced, is the following.
Theorem 1**.**
Let be the curvature of the initial condition to (1). Then there exists a constant such that if
[TABLE]
then the solution to the normalised -curve shortening flow ( a positive integer), that is for the curvature of the curves given by the rescaled embedding , with rescaled time parameter , where is the maximum time of existence for (1), it holds that
[TABLE]
*where is a constant that only depends on and . *
Together with Theorem I1.1 from [1], this gives the following
Theorem 2**.**
For any simple convex curve as initial data, the normalised version of (1), converges towards a circle smoothly and the curvature of the normalised embeddings satisfy
[TABLE]
where is a constant that only depends on and the curvature of the initial condition.
Indeed, by the theorem of Andrews referred to above, (2) eventually holds if we start (1) with a given convex simple curve as initial data. The rate of convergence given by our main result seems to be the sharpest possible rate of convergence for the -curve shortening flow (see the remark at the end of this paper).
A naive idea for proving Theorem 1 would be to use the Parabolic PDE to which the normalised version of -curve shortening flow is equivalent to (see equation (14) in Section 3.2), and then linearise around the steady solution to obtain exponential convergence. However, if we linearise around the steady solution, the elliptic part of the parabolic operator corresponding to the -curveshortening flow has a negative eigenvalue, so no exponential convergence should be expected (see the discussion in [5] right after Theorem 2.2, and notice that when , a negative eigenvalue occurs). The good news here is that, being a curvature, it satisfies an important identity, which is responsible for us being able to get this exponential convergence.
Our methods are based on the techniques employed in [5], that is to say on the Fourier method. Hence, we will transform our problem into (finite dimensional) approximations of an infinite dimensional dynamical system, for which appropriate estimates will be proved, and which will finally lead to a proof of Theorem 1, proof which is given in the final section of this paper. The intermediate sections are devoted to show these appropriate estimates, which, in short, amount to controlling the Fourier wavenumbers of a solution to (1) in terms of the average of the curvature; from this we will be able to show a time decay for the Fourier wavenumbers different from the average (which in fact blows-up), and which, as we said before, will lead to a proof of Theorem 1.
2. Basic definitions and notation
When the initial curve is convex, the -curve shortening flow is equivalent to the following Boundary Value Problem:
[TABLE]
, with periodic boundary conditions, and a strictly positive function. Notice that the Maximum Principle implies that must remain positive for all times (i.e. a convex curve remains convex). We will need to compute finite dimensional approximations of the previous partial differential equation in Fourier space, so we must establish some definitions and notation. Recall that For , its Fourier expansion is given by:
[TABLE]
where,
[TABLE]
We shall refer to as the Fourier wavenumbers of .
We will also adopt the notation
[TABLE]
[TABLE]
and define the following sets
[TABLE]
[TABLE]
and,
[TABLE]
From now on will denote a finite set of integers which contains 0 (i.e, ), and which is symmetric around 0 (i.e., if then ).
Using this notation, in Fourier space, the -curve shortening flow can be approximated by the following finite dimensional dynamical system:
[TABLE]
with initial condition
[TABLE]
Formally, the in the system right above should bear, for instance, a subindex which makes its dependence on explicit, but as this is understood from now on, we will suppress it in what follows (and as our estimates will not depend on , this should be of no importance).
Notice also that (4) is an autonomous system, so there is a unique and smooth solution for a short time (see [4]). We will also make use of the seminorms , which are defined as in [5] as follows:
[TABLE]
As usual, we define , as the space of functions with continuous derivatives of order , equipped with the norm
[TABLE]
3. Technical Lemmas and intermediate results
We shall follow closely the arguments presented in [8]. Therefore we must show that for given a solution to (4), we can control the Fourier wavenumbers different from [math] in terms of the [math]-th wavenumber. The fact that the first eigenvalue is the main difficulty we must face, as this makes difficult to control the - wave number in terms of the [math] th wave number. Once we have done this, all that is left is to follow the arguments presented in [5, 8]. The key to our proofs is that a curvature function of a locally convex curve satisfies
[TABLE]
since this identity, once we have control over the higher Fourier wavenumbers (those with ) assuming control over the wavenumbers, allows us to control the Fourier wavenumbers.
The careful reader must notice that the proofs given in this paper, our estimates are given for system (4), and that this estimates are independent of , this allows us to take a limit so the results are valid for the full system (3).
3.1. Controlling the Fourier wavenumbers
We start with a technical lemma.
Lemma 1**.**
There is a such that if the initial condition of (3) satisfies:
[TABLE]
and
[TABLE]
for , then is non decreasing.
Proof.
From the hypothesis of the lemma,
[TABLE]
and the implicit constnat in the big notation does not depend on . Hence, for small enough, the term dominates the the term in differential equation for . This implies that , and the conclusion of the lemma follows.
∎
Now we show some control estimates for the Fourier modes,
Lemma 2**.**
There is a such that if the initial condition of (3) satisfies:
[TABLE]
holds, and for
[TABLE]
and
[TABLE]
Then
[TABLE]
Proof.
Let us consider the quantity and we prove that is nonincreasing for fixed.
We compute:
[TABLE]
where the terms are given by:
[TABLE]
We bound ,
[TABLE]
Now if , then
[TABLE]
and hence
[TABLE]
with independent of . Since , we get
[TABLE]
Splitting the sums and using similar calculations as in (5), we obtain
[TABLE]
and again is independent of .
Since the sum of the absolute value of all these terms can be made smaller than
[TABLE]
for by taking small enough, then the term is non increasing for small enough. From this the conclusion of the lemma follows. ∎
As we have been doing so far, in what follows the wave numbers are restricted to a fixed but arbitrary set , so keep this in mind. And, as announced at be beginning of this section, since the estimates are independent of , a limiting procedure will give the result for when we take as the whole set of integers.
Lemma 3**.**
Let be such that . There is a such that if and then whenever for all , for all , then we also have for on the same time interval.
Proof.
Since we have that , we can choose a such that on (remember we are working with an arbitrary but final dimensional approximation of the -curve shortening flow). We have the following identity
[TABLE]
where . It can be easily seen that for the Fourier modes of we have:
[TABLE]
Now taking the Fourier transform, this implies
[TABLE]
Since is the curvature of a convex curve, we have , so
[TABLE]
then
[TABLE]
In order to estimate the sum in the right side, let us notice that
[TABLE]
Now we proceed to estimate ,
[TABLE]
where is a constant independent of and .
As we have that , we get
[TABLE]
Therefore
[TABLE]
as long as . ∎
Lemma 4**.**
There is a , independent of , such that if the initial condition of (3) satisfies and
[TABLE]
then for all times (as long as the solution to (4) exists),
[TABLE]
Proof.
There exists a be such that the interval is maximal with respect to the following property: for all we have
[TABLE]
and
[TABLE]
Now applying Lemmas 1 and 2, we obtain that
[TABLE]
whenever . Now Applying Lemma 3, we get that
[TABLE]
when . Hence we have that and if we apply the same arguments as before we can show that there is a such that if then , contradicting the maximality of . ∎
For small enough, assuming that for the initial condition we have , we have now control over all the Fourier wavenumbers of the solution. The arguments in [5] now apply almost verbatim: see the upcoming sections.
3.2. Decay of the Fourier wavenumbers
Again, all the estimates proved in this section are valid for any choice of , and are also independent of the choice. Our main purpose is to show that the Fourier wavenumbers , , go to 0 as . To begin, we have, as in [5], the Trapping Lemma (Lemma 3.2 in [5]). Keep in mind that we are always under the assumption that is the curvature function of a simple convex closed curve (or equivalently, the identity holds).
Theorem 3** (Trapping Lemma).**
There exists a constant independent of the choice of such that if the initial datum satisfies the following inequality:
[TABLE]
then there exists a that depends on such that the solution to (4) satisifies:
[TABLE]
Also, in the same way as Lemmas 3.3 and 3.4 are obtained in [5], we have a Blow-up Lemma.
Lemma 5** (Blow-up).**
There is a (the same as in the Trapping Lemma) such that if the initial condition of (3) satisfies
[TABLE]
then there are constants a number such that:
[TABLE]
From now on, we assume that satisfies
[TABLE]
where is such that the Trapping Lemma holds.
We also have a few important observations. First, integrating the ODE for , we obtain
[TABLE]
where is given by:
[TABLE]
Applying the Trapping Lemma, we get
[TABLE]
Also, from the Trapping Lemma, there exists such that
[TABLE]
We shall use these observations in proving the following decay (in time) estimates for the Fourier wave numbers of .
Proposition 6**.**
There exists which depends on , and a that depends also on and on , such that if then there is a constant such that for any , for , the following estimate holds for the solution of (3),
[TABLE]
Proof.
(See also the proof of Lemma 3.6 in [5]) First we have
[TABLE]
where
[TABLE]
We are going to estimate the term inside the integral in the last inequality. As before we split into sums of the form
[TABLE]
using this, the Trapping Lemma and the observations after its statement, we get
[TABLE]
and since , we finally obtain
[TABLE]
Then we have
[TABLE]
Without lost of generality we can assume . Now using again (10) and the fact that , obtain
[TABLE]
then we have,
[TABLE]
for any . ∎
Next we are going to improve on the decay estimates of the Fourier coefficients. To be able to do this, we will need the following lemma.
Lemma 7**.**
There is a such that if , then we have the estimates
[TABLE]
and
[TABLE]
Proof.
We have that the following differential inequality
[TABLE]
holds for a constant independent of . This is equivalent to
[TABLE]
Using Lemma 5, from the previous differential inequality we obtain
[TABLE]
The result follows by integration. For the other inequality, notice that there exists a constant so we also have the following differential inequality
[TABLE]
∎
In order to proceed we use previous proposition to estimate the integral
[TABLE]
from below, from previous lemma, since is small, using Taylor’s Theorem, this shows
[TABLE]
and using this and (9), we get
[TABLE]
here
[TABLE]
Using the estimate of the proposition and the fact that if then has at least two entries different from [math],we get
[TABLE]
If we introduce the bound from Proposition 6 in (3.2) we get,
[TABLE]
for a well chosen , and from which we obtain the estimate (we use again the fact )
[TABLE]
If we assume that using this new bound an pluggin it into (3.2), we improve again our estimate on :
[TABLE]
,
Finally , if we repeat this procedure a finite number of times we arrive at
[TABLE]
where is a constant independent of , and is defined by (12) (so is value its ).
Now, we must show now that the wavenumbers satisfy the same estimate. In this case we write
[TABLE]
And we have an identity which follows from (here we use that is the conjugate of , as is real valued)
[TABLE]
where the prime (′) in the inner sum of the righthand side indicates that at least one of the . Using similar computations as in the proof of Lemma 3, together with (13), we can conclude that
[TABLE]
hence, if is small enough, we can deduce that
[TABLE]
which is just that . So we have proved
Proposition 8**.**
Let is a smooth -periodic function which satisfies that . There exists a positive constant such that if
[TABLE]
then a solution to (3) satisfies
[TABLE]
where is the blow-up time and is a constant that depends only on and .
We normalize the solution of (3) by means of the following transformation:
[TABLE]
Applying chain rule, we obtain the following normalized version of (3)
[TABLE]
Using this normalisation, Proposition 8 translates into:
Corollary 9**.**
Let is a smooth -periodic function which satisfies that . There exists a positive constant such that if
[TABLE]
then the normalization of satisfies:
[TABLE]
where is a positive constant that depends only on and .
4. Exponential convergence of the normalised curvature towards 1: Proof of the main result
We will need the following improvement over Lemma 7.
Lemma 10**.**
The following estimates hold
[TABLE]
and
[TABLE]
Proof.
Notice that using (13) and the equation satisfied by , we have the differential inequality, which is valid for a constant
[TABLE]
Integrating, from to , we obtain the second inequality. Analogously for a constant , we have the differential inequality
[TABLE]
which by integration gives the first inequality. ∎
Finally we have our main result.
Theorem 4**.**
Let be the initial condition of (3) (so it is the curvature function of a convex simple curve). Then there exists a constant such that if
[TABLE]
then the solution to (14) satisfies
[TABLE]
where is a constant that only depends on the initial condition and and .
Proof.
Let
[TABLE]
now we compute
[TABLE]
Analogously,
[TABLE]
In this case, , then
[TABLE]
Applying the triangular inequality and Corollary 9 we can conclude that
[TABLE]
for some constant that depends on and .
∎
4.1. Final Remarks
The rate of convergence obatined in Theorem, seems to be the best possible in general. We have not been able to produce an example where the rate given in the theorem is met; however, to justify our claim, we refer to the comments after the statement of Theorem 2.2 in [5]: The first positive eigenvalue of the elliptic part of (14), i.e., the left hand side of the equation, when linearised around the steady solution is precisely .
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] B. Andrews, Evolving convex curves, Calc. Var. Partial Differential Equations 7 (1998) 315–371.
- 2[2] B. Andrews and P. Bryan, Curvature bound for curve shortening flow via distance comparison and a direct proof of Grayson’s theorem. J. Reine Angew. Math. 653 (2011), 179–187.
- 3[3] K.-S. Chou and X.-P. Zhu, The curve shortening problem. Chapman and Hall/CRC, Boca Raton, FL, 2001. x+255 pp.
- 4[4] E. A. Coddington and N. Levinson, Theory of ordinary differential equations. Mc Graw-Hill Book Company, Inc., New York-Toronto-London, 1955. xii+429 pp.
- 5[5] J. C. Cortissoz and A. Murcia, On the stability of m-fold circles and the dynamics of generalized curve shortening flows. J. Math. Anal. Appl. 402 (2013), no. 1, 57–70.
- 6[6] M. Gage and R.S. Hamilton, The heat equation shrinking convex plane curves. J. Differential Geom. 23 (1986), no. 1, 69–96.
- 7[7] R. L. Huang, Blow-up rates for the general curve shortening flow. J. Math. Anal. Appl. 383 (2011), no 2, 482–487.
- 8[8] A. Murcia, The Ricci Flow On Surfaces with boundary and Quasilinear Evolution Equations in 𝕊 1 superscript 𝕊 1 \mathbb{S}^{1} . Ph D thesis, Universidad de los Andes, 2015.
