On a degenerate parabolic system describing the mean curvature flow of rotationally symmetric closed surfaces
Harald Garcke, Bogdan-Vasile Matioc

TL;DR
This paper reformulates the mean curvature flow of rotationally symmetric closed surfaces as a coupled system of an evolution PDE and ODEs, establishing well-posedness and smoothing effects for the degenerate parabolic problem.
Contribution
It introduces a new formulation of the mean curvature flow for symmetric surfaces and proves well-posedness and smoothing properties for the resulting nonlinear degenerate parabolic system.
Findings
Well-posedness of the formulated evolution problem
Parabolic smoothing effect demonstrated
Coupled PDE-ODE system accurately describes the flow
Abstract
We show that the mean curvature flow for a closed and rotationally symmetric surface can be formulated as an evolution problem consisting of an evolution equation for the square of the function whose graph is rotated and two ODEs describing the evolution of the points of the evolving surface that lie on the rotation axis. For the fully nonlinear and degenerate parabolic problem we establish the well-posedness property in the setting of classical solutions. Besides we prove that the problem features the effect of parabolic smoothing.
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.
On a degenerate parabolic system describing the mean curvature flow of rotationally symmetric closed surfaces
Harald Garcke
Fakultät für Mathematik, Universität Regensburg, 93040 Regensburg, Deutschland.
and
Bogdan–Vasile Matioc
Fakultät für Mathematik, Universität Regensburg, 93040 Regensburg, Deutschland.
Abstract.
We show that the mean curvature flow for a closed and rotationally symmetric surface can be formulated as an evolution problem consisting of an evolution equation for the square of the function whose graph is rotated and two ODEs describing the evolution of the points of the evolving surface that lie on the rotation axis. For the fully nonlinear and degenerate parabolic problem we establish the well-posedness property in the setting of classical solutions. Besides we prove that the problem features the effect of parabolic smoothing.
Key words and phrases:
Mean curvature flow, degenerate parabolic equation, maximal regularity, parabolic smoothing
2010 Mathematics Subject Classification:
35K55, 53C44, 35R35, 35K93
1. Introduction
Mean curvature flow is the most efficient way to decrease the surface area of a hypersurface. It hence has been of great interest in geometry as well as in materials science and image analysis, see [29, 8, 20, 9, 14, 18, 24]. Since the pioneering work of Brakke [8] and Huisken [20] many results have been shown for mean curvature flow and we refer to [24] and the references therein for more information about the subject. The case of rotationally symmetric evolutions lead to spatially one-dimensional problems and due to the reduced complexity this situation has been studied by several authors analytically [13, 21, 1, 25, 22] as well as numerically [27, 7].
In particular, rotationally symmetric mean curvature flow has been helpful to understand singularity formation in curvature flows, see [21, 13, 1, 25, 15]. Most of the analytical results have been restricted to the case of surfaces with boundary or periodic unbounded situations. The situation becomes analytically far more involved if one considers closed surfaces, i.e., compact surfaces without boundaries. In this context the governing equation can be recast, provided that the points on the rotation axis have positive curvature, as a free boundary problem which involves both degenerate and singular terms. This paper gives first well-posedness and parabolic smoothing results for the free boundary problem describing compact rotationally symmetric surfaces evolving by mean curvature derived herein.
Let us now precisely formulate the analytic problem. We study the evolution of a family of rotationally symmetric surfaces by the mean curvature flow. Given , we assume that
[TABLE]
is the surface obtained by rotating the graph of the unknown function around the -axis. Moreover, we consider herein the case when the surfaces are closed, meaning in particular that also the domain of definition of is unknown. Since the motion of the surfaces is governed by the equation
[TABLE]
where is the normal velocity of and the mean curvature of , with denoting the principle curvatures of , we obtain the following evolution equation for the unknown function :
[TABLE]
The evolution of the boundaries: The first approach. As the functions and are unknown, we have to derive equations describing the evolution of these two boundaries. If we want to evaluate the normal velocity at , it follows from (1.2e) that The next goal is to express at in terms of To this end we assume, for some , that
[TABLE]
is invertible with the inverse function , so that in particular is a -surface close to . Then due to the fact that for we have and
[TABLE]
hence Noticing that
[TABLE]
in the case when we obtain the following relation
[TABLE]
Similarly, assuming that, for some ,
[TABLE]
is invertible and , we find for the evolution equation
[TABLE]
The evolution of the boundaries: The second approach. A major drawback of the (formally) quasilinear parabolic equation (1.2a) is that the boundary conditions (1.2e) make the equation highly degenerate as:
- (I)
The diffusion coefficient vanishes in the limit and ;
- (II)
The term becomes unbounded for and .
In order to overcome (II) we introduce, motivated also by (1.2f)-(1.2g), a new unknown via
[TABLE]
Then also vanishes at the boundary points and and (1.2a) can be expressed as
[TABLE]
This equation is also (formally) quasilinear parabolic and also degenerate – as the diffusion coefficient vanishes for and , cf. (1.4) below – but now none of the coefficients is singular. In order to obtain an evolution equation also for the functions describing the boundaries, we assume that
[TABLE]
and
[TABLE]
Note that (1.4) implies in particular that the corresponding function satisfies Furthermore, (1.5) is a nonlinear boundary condition for which is equivalent to our former assumption that is a -surface, cf. Lemma A.1. Differentiating now the relation , , with respect to time, it follows in virtue of (1.3) and (1.5), that
[TABLE]
These are the very same relations as in (1.2f)-(1.2g). It is not difficult to see, cf. Lemma A.1, that the two approaches are equivalent.
Summarizing, we may formulate the problem by using as an unknown and we arrive at the evolution problem
[TABLE]
In the following we use the formulation (1.6) in order to investigate the mean curvature flow (1.1). We are interested here to prove the existence and uniqueness of solutions which satisfy the equations in a classical sense (a weak formulation of (1.6) is not available yet). The formulation (1.6) has two advantages compared to the classical approach followed in [17, 24, 12] for example. Firstly, the equations are explicit (we do not need to work with local charts) and, secondly because the maximal solutions to (1.6) are defined in general on a larger time interval compared to the ones in [17, 24, 12] (the solutions in [17, 24, 12] exist only in a small neighborhood of a fixed reference manifold over which they are parameterized). A disadvantage of our approach is that herein the initially surfaces are necessary of class while in [17] only -regularity, for some fixed , is required. An interesting research topic which we next plan to follow is to determine initial data for which the closed rotationally symmetric surface evolves such that neck pinching at the origin occurs in finite time. This topic has been already studied in the context of (1.2a), but in the special setting when and are kept fixed, the function is strictly positive, and suitable boundary conditions (either of Neumann or Dirichlet type) are imposed at these two fixed boundary points, cf. [25, 15, 13, 21, 28]. In the context of closed surfaces without boundary considered herein there are several results establishing the convergence of initially convex surfaces towards a round point in finite time, cf. e.g. [28, 3, 20], but to the best of our knowledge no result establishing neck pinching at the origin is available. It is worth to emphasize that in this context however, by using maximum principles and some explicit solutions to the mean curvature flow, such as spheres, hyperboloids, or shrinking donuts, there are several examples of dumbbell shaped surfaces which develop singularities in finite time, cf. [6, 19, 14].
Remark 1.1*.*
- (i)
If , , then the surfaces under consideration are spheres and the radius solves the ODE
[TABLE]
- (ii)
The conditions (1.4)-(1.5) impose some restrictions on the initial data. Lemma A.1 shows that any rotationally symmetric surface of class with mean curvature that does not vanish at the points on the rotation axis satisfies (1.4)-(1.5). These properties are then preserved by the flow.
We will solve the degenerate parabolic system in the setting of small Hölder spaces. The small Hölder space , , is defined as the closure of the smooth periodic functions (or equivalently of , ) in the classical Hölder space of -periodic functions on the line with -Hölder continuous -th derivatives. Besides, , denotes the subspace of consisting only of even functions. By definition, the embedding , , is dense and moreover it holds
[TABLE]
Here denotes the continuous interpolation functor introduced by Da Prato and Grisvard [11].
The main result of this paper is the following theorem.
Theorem 1.2**.**
Let be a fixed Hölder exponent, , and let be positive in such that and
[TABLE]
Then the evolution problem (1.6) has a unique maximal solution such that
[TABLE]
where
[TABLE]
and . Moreover, it holds that
[TABLE]
Remark 1.3*.*
- (i)
The choice of the small Hölder spaces is essential. Indeed, using a singular transformation from [4], we may recast the evolution problem (1.6) as a fully nonlinear evolution equation with the leading order term in having in the linearisation - when working within this class of functions - a positive and bounded coefficient. Besides, the setting of small Hölder spaces is a smart choice when dealing with fully nonlinear parabolic equations, cf. e.g. [11, 23]. A further departure of these spaces from the classical Hölder spaces is illustrated in Lemma A.2.
- (ii)
The problem considered in [4] is general enough to include also (1.6). However, the technical details, see Section 3, are different from those in [4] and also simpler. Besides, the parabolic smoothing property for in Theorem 1.2 is a new result in this degenerate parabolic setting and it extends also to the general problem considered in [4]. In particular, this proves in the context of the porous medium equation (which is the equation that motivates the analysis in [4]) that the interface separating a fluid blob, that expends under the effect of gravity, from air is real-analytic in the positivity set, see [31] for more references on this topic.
- (iii)
If satisfies in , and
[TABLE]
then can be chosen as an initial condition in (1.6). However, the initial data in Theorem 1.2 are not required to be twice differentiable at and . For example if , then
[TABLE]
can be chosen as an initial condition in (1.6) as , but .
- (iv)
The assumption that guarantees that the nonlinear boundary condition holds at . Indeed, since
[TABLE]
l’Hospital’s rule shows that .
2. The transformed problem
In order to study (1.6) we use an idea from [4] and transform the evolution problem (1.6) into a system defined in the setting of periodic functions by using a diffeomorphism that has a first derivative which is singular at the points and . More precisely, we introduce the new unknown
[TABLE]
where
[TABLE]
Given , is a -periodic function on which is even and merely the continuous differentiability of implies that
[TABLE]
In terms of the new variable the problem (1.6) can be recast as follows
[TABLE]
We point out that nonlinear boundary condition has not been taken into account in the transformed system (2.1). This is due to the choice of the function spaces below as, similarly as in Remark 1.3 (iv), requiring that ensures that holds at time .
In order to study (2.1) we choose as an appropriate framework the setting of periodic small Hölder spaces. For a fixed we define the Banach spaces
[TABLE]
with the corresponding norms It is important to point out that the embedding is dense. Though at formal level the equation has a quasilinear structure, our analysis below shows that the problem (2.1) is actually (as a result of the boundary conditions) fully nonlinear (see Lemma 2.1 and the subsequent discussion). This loss of linearity is however compensated by the fact that none of the terms on the right hand side of is singular when choosing . Moreover the function multiplying in is now -Hölder continuous and positive.
Lemma 2.1**.**
The operators
[TABLE]
are bounded.
Proof.
See [4, Lemma 2.1]. ∎
We emphasize that it is not possible to choose in Lemma 2.1 as target space a small Hölder space with . In particular, the terms and on the right-hand side of have the same importance as when linearizing this expression.
We now set
[TABLE]
Then, is an open subset of . Let further
[TABLE]
be the operator defined by
[TABLE]
It is not difficult to check that so that is well-defined. In virtue of Lemma 2.1 it further holds that
[TABLE]
Hence, we are led to the fully nonlinear evolution problem
[TABLE]
with . We shall establish the existence and uniqueness of strict solutions (in the sense of [23]) to (2.3) by using the fully nonlinear parabolic theory presented in the monograph [23]. To this end we next identify the Fréchet derivative and we prove that it generates, for each , a strongly continuous and analytic semigroup. In the notation of Amann [2] this means by definition
[TABLE]
In fact, in view of [2, Corollary I.1.6.3], we only need to show that the partial derivative generates a strongly continuous analytic semigroup in Given , it holds that
[TABLE]
where
[TABLE]
We note that , , with being positive. Moreover, it holds that .
Since
[TABLE]
the operator may be viewed as being a lower order perturbation of , cf. [2, Theorem I.1.3.1 (ii)]. The following result enables us to regard also other terms of as being lower order perturbations.
Lemma 2.2**.**
Let . Then, given , there exists a constant such that
[TABLE]
Proof.
Letting , it is not difficult to verify that
[TABLE]
In view of this equivalence, the claim for follows from the observation that
[TABLE]
where the functions and belong to and respectively. The proof of the second claim follows by similar arguments. ∎
Recalling that , Lemma 2.2 implies that also can be viewed as being a perturbation. Let us now notice that
[TABLE]
Observing that
[TABLE]
where , we may regard in view of Lemma 2.2 also the operator
[TABLE]
as being a perturbation and we are left to prove the generator property for
[TABLE]
In fact, it suffices to establish the generator property for the operator
[TABLE]
where we have dropped the lower order term Indeed, assuming that , it follows . This latter property is equivalent to the existence of constants and such that
[TABLE]
cf. [2, Chapter I]. The relation holds in particular for . In order to conclude that we are thus left to show that is an isomorphism too. Hence, given , for with we set
[TABLE]
Taking into account that it follows that A simple computation shows that , so that also . We may thus conclude that , so that holds also when replacing with , . The nontrivial property is established in detail in Section 3 below, cf. Theorem 3.1.
3. The generator property
The first goal of this section is to establish Theorem 3.1, which is a main ingredient in the proof of the main result.
Theorem 3.1**.**
Given , it holds that .
We consider for , with sufficiently small, partitions of the interval and corresponding families with the following properties
- •
in ;
- •
, , {\rm supp\,}(\pi_{3}^{\varepsilon})=I\setminus\big{(}[-2\varepsilon,2\varepsilon]\cup[\pi-2\varepsilon,\pi+2\varepsilon]);
- •
on , ;
- •
, , {\rm supp\,}(\chi_{3}^{\varepsilon})=I\setminus\big{(}[-\varepsilon,\varepsilon]\cup[\pi-\varepsilon,\pi+\varepsilon]);
- •
and are even on ;
- •
has an even and periodic extension in .
Extending and by zero in , we may view these functions as being smooth and even functions on .
As a first step towards proving Theorem 3.1 we approximate locally by certain operators which are simpler to analyze.
Lemma 3.2**.**
Let be given. Then, there exists , a constant , and a partition such that the operator introduced in (2.6) satisfies
[TABLE]
for and , where
[TABLE]
Proof.
Observing that , it follows that
[TABLE]
which proves (3.1) for .
Furthermore, it holds that
[TABLE]
where
[TABLE]
Using , we now obtain
[TABLE]
provided that is sufficiently small.
Concerning the second term we write
[TABLE]
where
[TABLE]
The arguments in the proof of Lemma 2.2 yield
[TABLE]
Besides, since , we get
[TABLE]
Finally, it is not difficult to see that the function
[TABLE]
satisfies . Therewith we have
[TABLE]
and we conclude that
[TABLE]
provided that is sufficiently small. This proves (3.1) for . The proof of the claim for is similar and we therefore omit it. ∎
We now consider the operators , , found in Lemma 3.2 in suitable functional analytic settings. Regarding as an element of , it is well-known that generates an analytic semigroup in . In particular, there exist constants and such that
[TABLE]
cf. [2, Theorem I.1.2.2]. The operator can be viewed as an element of 111For a definition of , , see [23]. Again, , , denotes the closed subspace of consisting of even functions.. Furthermore, in this context appears as the restriction of to the subset of rotationally symmetric functions. Indeed, given , let
[TABLE]
One can show that the radially symmetric function belongs to and that
[TABLE]
with a constant independent of . Recalling that , cf. [23, Theorem 3.1.14 and Corollary 3.1.16], there exist constants and such that
[TABLE]
In particular it holds that
[TABLE]
for and Moreover, in virtue of
[TABLE]
we conclude that
[TABLE]
The constants and can be chosen such that (3.6) holds true also when replacing by where , denotes the right translation by .
In particular (2.5), (3.5), and (3.6) ensure there exists and such that
[TABLE]
for all and all . The estimate (3.7) together with the observation that the map
[TABLE]
defines a norm on which is equivalent to the standard Hölder norm are essential for establishing the following result.
Lemma 3.3**.**
There exist and such that
[TABLE]
for all and all .
Proof.
Letting and denote the constants in (3.7), we chose in Lemma 3.2. Lemma 3.2 together with (3.7) yields
[TABLE]
for , , and In virtue of (3.8) and of (2.5) it now follows that there exists a constant such that
[TABLE]
for and Finally, the interpolation property (1.7), the latter estimate, and Young’s inequality ensure that there exist constants and such that (3.9) is satisfied. ∎
To derive the desired generation result we are left to show that To this end we infer from (3.9) that is one-to-one. Having shown that is a Fredholm operator of index zero, the isomorphism property follows then in view of the compactness of the embedding
Lemma 3.4**.**
* is a Fredholm operator of index zero.*
Proof.
Since , the equation is equivalent to
[TABLE]
hence . The kernel of consists thus only of constant functions.
It is easy to see that the range of is contained in
[TABLE]
which is a closed subspace of of codimension . To show that the range of coincides with we associate to the function
[TABLE]
Using the property defining , it is not difficult to check that is twice continuously differentiable with
[TABLE]
The second last identity above follows by using appropriate substitutions in the second integral. Moreover, it holds that , , and
[TABLE]
as we may extend by periodicity to . Some standard (but lengthy) arguments show that lies in , which implies that . Thus, belongs to the range of and the claim follows. ∎
Proof of Theorem 3.1.
In view of Lemma 3.3 it remains to show that is an isomorphism. This property is an immediate consequence of the estimate (3.9), which implies in particular that is injective, and of the fact that is a Fredholm operator of index zero, cf. Lemma 3.4 (we recall at this point that the embedding is compact). ∎
We conclude this section with the proof of the well-posedness result stated in Theorem 1.2. The proof of the parabolic smoothing property for the function is postponed to Section 4.
Proof of Theorem 1.2.
We first address the solvability of (2.1). As a direct consequence of Theorem 3.1 we have that
[TABLE]
for all . Recalling also (2.2) and the interpolation property of the small Hölder spaces (1.7), the assumptions of [23, Theorem 8.4.1] are all satisfied in the context of (2.1). Hence, for each , (2.1) possesses a unique maximal strict solution
[TABLE]
such that
[TABLE]
where . Since by assumption , the existence and uniqueness claim in Theorem 1.2 follows. That is a straight forward consequence of [23, Corollary 8.4.6]. The real-analyticity property for (or ) is however more subtle and is established in Section 4 below. ∎
4. Parabolic smoothing
In the following we consider a solution to (1.6) with maximal existence time as found in Theorem 1.2. and we prove that the associated function
[TABLE]
is real-analytic. In this way we establish the parabolic smoothing property for the function as stated in Theorem 1.2. The proof below exploits a parameter trick which has been used in other variants also in [5, 10, 16, 30, 26] to improve the regularity of solutions to parabolic or elliptic equations. The degenerate parabolic setting considered herein raises new difficulties, in particular due to the fact that the solutions vanish at [math] and , which hinder us to establish real-analyticity of in a neighborhood of these points.
To start, we fix an arbitrary constant such that . Given with
[TABLE]
and , we introduce the function with
[TABLE]
The smallness condition (4.1) ensures that is a real-analytic diffeomorphism. We associate to the function , , , Let further . Since for and , Theorem 1.2 yields
[TABLE]
Clearly, is -periodic and even with respect to . Observing that
[TABLE]
tedious computations show that
[TABLE]
We emphasize that Lemma A.2 (ii) plays a key role in the proof of (4.2). Moreover, given , it holds that
[TABLE]
together with
[TABLE]
Furthermore, the pair solves the parameter dependent evolution problem
[TABLE]
where is defined by
[TABLE]
and
[TABLE]
Recalling (4.1), it then follows that
[TABLE]
Observing that is a bounded operator which can be estimated in a similar way as the operators in Lemma 2.2, we may repeat the arguments in Sections 2-3 to conclude that
[TABLE]
for all Applying [23, Theorem 8.4.1], it follows that (4.3) possesses for each a unique maximal strict solution with
[TABLE]
where is the maximal existence time. In view of [23, Corollary 8.4.6] we may conclude that the mapping
[TABLE]
where
[TABLE]
is real-analytic. Let now be fixed. Since is a real-analytic map, we obtain for the function determined by the solution considered above, in particular that
[TABLE]
is real-analytic too. Additionally, given , for sufficiently small it holds that
[TABLE]
is well-defined and real-analytic. Here we use the real-analyticity of in which we already established. Composing the mappings (4.4) and (4.5), it follows in view of the fact that is arbitrary that
[TABLE]
is real-analytic. Recalling that the property
[TABLE]
follows at once.
Appendix A
The next result shows that the two approaches used in the Introduction to derive evolution equations for the functions and require the same assumptions. In particular, it shows that the solutions to the problem (1.6) describe closed -surfaces without boundary and with positive curvature at the points on the rotation axis.
Lemma A.1**.**
Let and let satisfy for all and . Then, the following are equivalent:
- (i)
* and*
[TABLE]
- (ii)
There exists such that is invertible and the inverse satisfies , , and .
**
- (iii)
The function satisfies , , and
Proof.
We first prove the implication (i)(ii). It is obvious that if is sufficiently small, then has an inverse function that satisfies Furthermore, it holds that and
[TABLE]
Hence, is twice differentiable in [math] and . Furthermore, the mean value theorem yields the existence of a sequence such that . Since
[TABLE]
we obtain the following relation
[TABLE]
and therewith we get that
We now prove the implication (ii)(iii). We may assume that for . Invoking (A.1) we get that and . Moreover it holds that
[TABLE]
and this proves (iii).
We conclude with the proof of (iii)(i). The relations and are immediate and together with
[TABLE]
we have completed the proof. ∎
Lemma A.2 provides a continuity result which is used to establish (4.2). This lemma also exemplifies why the small Hölder spaces are to be preferred in certain applications to the classical ones.
Lemma A.2**.**
Let .
- (i)
Given the mapping
[TABLE]
is in general not continuous.
**
- (ii)
Given the mapping
[TABLE]
is continuous.
Proof.
It is easy to verify that is well-defined. The following example shows that this (nonlinear) mapping is in general not continuous. Indeed, let be a function which satisfies on and in . The -periodic extension of
[TABLE]
satisfies Given , let and denote the -periodic extensions of
[TABLE]
It then holds and in . Since
[TABLE]
it follows that
[TABLE]
which proves (i).
We now prove (ii). Let thus with as in and be a sequence with in It follows that for . Moreover, it holds
[TABLE]
where
[TABLE]
These estimates show that in and that each ball in centered in contains a function with suitably large. Since , it follows that also and this completes the proof. ∎
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] S. Altschuler, S. B. Angenent, and Y. Giga. Mean curvature flow through singularities for surfaces of rotation. J. Geom. Anal. , 5(3):293–358, 1995.
- 2[2] H. Amann. Linear and Quasilinear Parabolic Problems. Vol. I , volume 89 of Monographs in Mathematics . Birkhäuser Boston, Inc., Boston, MA, 1995. Abstract linear theory.
- 3[3] B. Andrews and J. Mc Coy. Convex hypersurfaces with pinched principal curvatures and flow of convex hypersurfaces by high powers of curvature. Trans. Amer. Math. Soc. , 364(7):3427–3447, 2012.
- 4[4] S. Angenent. Local existence and regularity for a class of degenerate parabolic equations. Math. Ann. , 280(3):465–482, 1988.
- 5[5] S. B. Angenent. Nonlinear analytic semiflows. Proc. Roy. Soc. Edinburgh Sect. A , 115(1-2):91–107, 1990.
- 6[6] S. B. Angenent. Shrinking doughnuts. In Nonlinear diffusion equations and their equilibrium states, 3 (Gregynog, 1989) , volume 7 of Progr. Nonlinear Differential Equations Appl. , pages 21–38. Birkhäuser Boston, Boston, MA, 1992.
- 7[7] J. W. Barrett, H. Garcke, and R. Nürnberg. Variational discretization of axisymmetric curvature flows. Numer. Math. , 141(3):791–837, 2019.
- 8[8] K. A. Brakke. The motion of a surface by its mean curvature , volume 20 of Mathematical Notes . Princeton University Press, Princeton, N.J., 1978.
