Home Flat flow solution to the mean curvature flow with volume constraint
Article Open Access

Flat flow solution to the mean curvature flow with volume constraint

  • Vesa Julin ORCID logo EMAIL logo
Published/Copyright: April 24, 2024

Abstract

In this paper I will revisit the construction of a global weak solution to the volume preserving mean curvature flow via discrete minimizing movement scheme by Mugnai, Seis and Spadaro [L. Mugnai, C. Seis and E. Spadaro, Global solutions to the volume-preserving mean-curvature flow, Calc. Var. Partial Differential Equations 55 2016, 1, Article ID 18]. This method is based on the gradient flow approach due to Almgren, Taylor and Wang [F. Almgren, J. E. Taylor and L. Wang, Curvature-driven flows: a variational approach, SIAM J. Control Optim. 31 1993, 2, 387–438] and Luckhaus and Sturzenhecker [S. Luckhaus and T. Sturzenhecker, Implicit time discretization for the mean curvature flow equation, Calc. Var. Partial Differential Equations 3 1995, 2, 253–271] and my aim is to replace the volume penalization with the volume constraint directly in the discrete scheme, which from practical point of view is perhaps more natural. A technical novelty is the proof of the density estimate which is based on second variation argument.

MSC 2020: 35K93; 53C42

1 Introduction

A smooth family of set ( E t ) t 0 is said to evolve according to volume preserving mean curvature flow if the normal velocity V t is proportional to the mean curvature H E t as

(1.1) V t = - ( H E t - H ¯ E t ) on  E t ,

where H ¯ E t = E t H E t d n . Such a geometric equation has been proposed in the physical literature to model coarsening phenomena, where the system consisting of several subdomains evolves such that it decreases the interfacial area while keeping the total volume unchanged [7, 18]. From purely mathematical point of view the equation (1.1) can be seen as the L 2 -gradient flow of the surface area under the volume constraint [18]. One has to be careful in this interpretation as the Riemannian distance between two sets is in general degenerate [16]. In order to overcome this one may use the idea due to Almgren, Taylor and Wang [2] and Luckhaus and Sturzenhecker [14] and to view (1.1) as the gradient flow of the surface area with respect to a different, non-degenerate, distance. Using the gradient flow structure, one may then construct a discrete-in-time approximation to the solution of (1.1) via the Euler implicit method, also known as the minimizing movements scheme. By letting the time step to zero, one then obtains a candidate for a weak solution of (1.1) called flat flow, as the convergence is measured in terms of the “flat norm”. This method is implemented to the volume preserving setting in [19].

In [19] the authors observe that from technical point of view it is easier to replace the volume constraint of the problem with volume penalization, as this simplifies certain regularity issues at the level of the discrete approximation. My aim here is to show that one may construct the flat flow solution to (1.1) by implementing the volume constraint in the minimizing movements scheme directly and thus avoid the volume penalization.

Let me quickly recall the discrete minimizing movements scheme for (1.1). One defines a sequence of sets ( E k h ) k , with fixed time step h > 0 , iteratively such that E 0 h = E 0 , where E 0 is the given initial set, and E k + 1 h is a minimizer of the functional

P ( E ) + 1 h E d ¯ E k h d x under the constraint  | E | = | E k h | .

Here P ( E ) denotes the perimeter (generalized surface area) of the set E and d ¯ F is the signed distance function of the set F (see next section). One then defines an approximative flat flow solution to (1.1) ( E t h ) t 0 from the previous sequence by E t h = E k h for t [ k h , ( k + 1 ) h ) . Any cluster point of ( E t h ) t 0 is then defined as flat flow solution to (1.1). The advantage is that such a solution is defined for all times and for rough initial data. The main result in the paper is the existence of a flat flow solution.

Theorem 1.

Assume that E 0 R n + 1 is an open and bounded set with finite perimeter and let ( E t h ) t 0 be an approximative flat flow solution to (1.1) stating from E 0 (see Definition 2.1). Then there exists a family of bounded sets of finite perimeter ( E t ) t 0 and a subsequence h k 0 such that

lim h k 0 | Δ E t h k E t | = 0 for a.e.  t 0

and for every 0 < t < s it holds | E t | = | E 0 | , P ( E t ) P ( E 0 ) and

| Δ E t E s | C s - t ,

where C depends on the dimension and on E 0 . Moreover, if the initial set E 0 is C 1 , 1 -regular, then any such limit flow ( E t ) t 0 agrees with the unique classical solution of (1.1) as long as the latter exists.

The above theorem thus provides the existence of a flat flow solution and guarantees that this notion is consistent with the classical solution when the initial set is regular enough. The disadvantage of the flat flow is that it is not clear if it provides a solution to the original equation (1.1) in any weak sense after the first singular time. However, the conditional result in the spirit of Luckhaus and Sturzenhecker [14] holds also in this case.

Theorem 2.

Let ( E t h ) t 0 be an approximative flat flow solution to (1.1) and let ( E t h k ) t 0 be the converging subsequence in Theorem 1. Assume further that it holds

lim h k 0 P ( E t h k ) = P ( E t ) for a.e.  t 0 .

Then for n 6 the flat flow ( E t ) t 0 is a distributional solution to (1.1) (see Definition 4.5).

One may also try to view equation (1.1) as a mean curvature flow with forcing, where the forcing term depends on the flow itself. In this way one may try to use different methods to construct a solution to the equation see, e.g., [5, 6]. I also refer the recent work [13] for a weak-strong uniqueness result related to (1.1).

As I already mentioned, the flat flow is defined for all times and one may study its asymptotical behavior. Indeed, by using the methods from [10, 11] one may deduce the convergence of the flow in low dimensions. I will state this merely as a remark as it follows from the above methods without any modifications.

Remark 1.1.

Assume that E 0 n + 1 , with n 2 , is as in Theorem 1 and let ( E t ) t 0 be a limit flat flow. When n = 1 , the flow E t converges to a union of disjoint balls exponentially fast and when n = 2 the flow converges to a union of disjoint balls up to a possible translation of the components.

The main technical challenge in proving Theorem 1 is to obtain the sharp density estimate for the discrete flow. This is also the main technical novelty of this paper. There are several techniques to deal with the volume constraint in variational problems, e.g., by using the argument from [3] (see also [15, Lemma 17.21]) or from [8] (see also [4, 9]). However, due to the presence of the dissipation term in the energy it is not obvious how to apply these arguments in order to obtain sharp density estimates in terms of the time step h. I will use an argument which is based on the second variation condition of the energy to prove the density estimate in Proposition 3.1. After this the proof of Theorem 1 follows exactly as in [14, 19] and the consistency follows almost directly using the argument in [12]. The proof also provides the dissipation inequality and therefore the results in [10, 11] hold and one obtains the result stated in Remark 1.1. Finally, I would like to point out that this article is not self-consistent as I will take several well-known arguments for granted, in particular, in Section 4.

2 Preliminaries

In this section I will briefly introduce the notation, the definition of the flat flow solution and recall some of its basic properties.

Given a set E n + 1 , the distance function dist ( , E ) : n + 1 [ 0 , ) is defined, as usual, as

dist ( x , E ) := inf y E | x - y |

and denote the signed distance function by d ¯ E : n + 1 ,

d ¯ E ( x ) := { - dist ( x , E ) for  x E , dist ( x , E ) for  x n + 1 E .

Then clearly it holds dist ( , E ) = | d ¯ E | . I denote the ball with radius r centered at x by B r ( x ) and by B r if it is centered at the origin.

For a measurable set E n + 1 the perimeter in an open set U n + 1 is defined as

P ( E , U ) := sup { E div X d x : X C 0 1 ( U , n + 1 ) , X L 1 }

and write P ( E ) = P ( E , n + 1 ) . If P ( E ) < , then E is called a set of finite perimeter. For an introduction to the topic I refer to [15]. The reduced boundary of a set of finite perimeter E is denoted by * E and the generalized unit outer normal by ν E . Note that it holds P ( E , U ) = n ( * E U ) for open sets U. Recall also that if E is regular enough, say with Lipschitz boundary, then P ( E ) = n ( E ) . For a given vector field X C 1 ( n + 1 , n + 1 ) and a set of finite perimeter E denote the tangential divergence on E * as div τ X = div X - D X ν E , ν E . The distributional mean curvature H E L 1 ( * E , ) is defined via the divergence theorem such that for every test vector field X C 0 1 ( n + 1 , n + 1 ) it holds

* E div τ X d n = * E H E X , ν E d n .

I will consider a flat flow solution to (1.1) in the spirit of Almgren, Taylor and Wang [2] and Luckhaus and Sturzenhecker [14]. To this end, for a fixed h ( 0 , 1 ) and a given (open) set F n + 1 , I define the functional

(2.1) h ( E , F ) = P ( E ) + 1 h E d ¯ F d x .

The flat flow solution is defined analogously as in [19].

Definition 2.1.

Let E 0 n + 1 be an open and bounded set of finite perimeter and fix h ( 0 , 1 ) . Define the sequence of sets ( E k h ) k = 0 iteratively as E 0 h = E 0 and E k + 1 h is a minimizer of the problem

min { h ( E , E k h ) : | E | = | E 0 | } .

Moreover, define an approximative flat flow ( E t h ) t 0 for (1.1) starting from E 0 as

E t h = E k h for  t [ k h , ( k + 1 ) h ) .

One has to be careful in the definition of the functional (2.1) if the set F is merely a set of finite perimeter as its value depends on the choice of the representative of F. One may overcome this by choosing a proper representative of the set F. However, this is not necessary as the regularity theorem below implies that one may in fact assume the sets E k h to be open. The difference in Definition 2.1 to the scheme in [19] is that here the minimizing problem is under volume constraint. On one hand this makes the minimization problem more natural, but on the other hand, it makes the quantitative density estimates more difficult to prove.

For a given open and bounded set F n + 1 consider the minimization problem

(2.2) min { h ( E , F ) : | E | = | F | } ,

where h ( , F ) is defined in (2.1). One may use an argument similar to [9] or [15, Lemma 17.21] to remove the volume constraint in (2.2) and deduce that a minimizer of (2.2) is a minimizer also for

(2.3) min { h ( E , F ) + Λ ~ | | E | - | F | | } ,

when Λ ~ is chosen large. Note that the constant Λ ~ may have nonoptimal dependence on E and on h. However, the property (2.3) is enough to deduce qualitative regularity properties since it implies that the minimizer inherits the regularity from the theory of the perimeter minimizers [15]. One may also write the Euler–Lagrange equation and by standard calculations (see, e.g., [1]) we have the second variation condition. We state this in the following proposition.

Proposition 2.2.

Let F R n + 1 be an open and bounded set, fix h ( 0 , 1 ) and let E be a minimizer of (2.2). Then E can be chosen to be open, which topological boundary is C 2 , α -regular up to a relatively closed singular set which Hausdorff dimension is at most n - 7 . The regular part is exactly the reduced boundary * E .

The Euler–Lagrange equation

(2.4) d F h = - H E + λ ,

where λ R is the Lagrange-multiplier, holds point-wise on * E and in a distributional sense on E . The quadratic form associated with the second variation of the energy is non-negative, i.e., for all φ H 1 ( * E ) with * E φ d H n = 0 it holds

(2.5) * E | τ φ | 2 - | B E | 2 φ 2 d n + 1 h * E d ¯ F , ν E φ 2 d n 0 ,

where B E ( x ) denotes the second fundamental form of E at x * E .

Proof.

Since the argument is standard, I will only give the outline. As I already mentioned, the minimizer E is also a minimizer of problem (2.3) for some large constant Λ ~ , which depends on h and on E itself. This implies that the set E is a Λ-minimizer of the perimeter and thus the reduced boundary * E is relatively open, C 1 , α -regular hypersurface and the singular set E * E has dimension at most n - 7 (see [15]). The C 2 , α -regularity then follows from the Euler–Lagrange equation and from standard Schauder-estimates for elliptic PDEs.

One may obtain the second variation condition (2.5) by using the argument from [1]. Indeed, given a function φ C 0 2 ( * E ) with * E φ d n = 0 , we may construct a family of diffeomorphisms Φ t such that Φ 0 = id , | Φ t ( E ) | = | E | and t | t = 0 Φ t ( x ) ν E = φ . Then the inequality follows from the minimality of E as

2 t 2 | t = 0 h ( Φ t ( E ) , F ) 0

and following the standard calculation of the second variation (see, e.g., [1]). Finally, one obtains (2.5) for all φ H 1 ( * E ) by approximation argument and by the fact that the singular set has zero capacity. ∎

3 Density estimates

This section is the theoretical core of the paper. The aim is to prove the following density estimate.

Proposition 3.1.

Let F R n + 1 be an open and bounded set of finite perimeter, fix h ( 0 , 1 ) and let E be a minimizer of (2.2). Then there is a constant c > 0 , which depends on the dimension n, | F | and on P ( F ) such that for all r h and all x E it holds

min { | E B r ( x ) | , | B r ( x ) E | } c r n + 1

and for all r C 0 h , where C 0 1 , it holds

c r n P ( E , B r ( x ) ) C 1 r n ,

where C 1 depends also on C 0 . Moreover, the following estimates hold

H E L ( * E ) 1 c h 𝑎𝑛𝑑 d ¯ F L ( E ) c - 1 h .

It is interesting that in [19, Corollary 3.3] the authors obtain similar result for their scheme for a constant which is independent of P ( F ) .

I need several lemmas in order to prove Proposition 3.1 and therefore I postpone its proof to the end of the section. Before proceeding to technical details, I state a useful consequence of Proposition 3.1.

Proposition 3.2.

Let F , E R n + 1 be as in Proposition 3.1. Then there are constants C 1 , c > 0 and h 0 > 0 , depending on the dimension, | F | and P ( F ) such that E is the ( Λ , r ) -minimizer of the perimeter for Λ = C h , r = c h and for h < h 0 . To be more precise, for sets G R n + 1 with Δ E G B c h ( x 0 ) it holds

P ( E ) P ( G ) + C h | Δ E G | .

Proof.

The argument is standard but I recall it for the reader’s convenience. Let me first show that there is x E and c ~ > 0 such that for ρ = c ~ h it holds B ρ ( x ) E . Fix ρ and apply the Besicovitch covering theorem to find disjoint balls { B ρ ( x i ) } i = 1 N such that x i E and

(3.1) i = 1 N | B ρ ( x i ) | = N | B 1 | ρ n + 1 c | E | .

I claim that for some i = 1 , 2 , , N it holds B ρ 2 ( x i ) E . Indeed, if this is not the case, then Proposition 3.1 implies

P ( E , B ρ ( x i ) ) c ρ n

for all i. Since the balls are disjoint, one has by the above and by (3.1) that

P ( E ) i = 1 N P ( E , B ρ ( x i ) ) c N ρ n c | E | ρ c h .

This is a contradiction when h is small enough.

Fix x 0 and G as in the claim. Note that in general the set G does not have the same measure as E and one needs to modify it to G ~ with | G ~ | = | E | , e.g., by using the argument from [9] as follows. Assume that | G | < | E | (the case | G | > | E | follows from similar argument). Since B ρ ( x ) E , by decreasing ρ and r if needed, it holds B ρ ( x ) G . By continuity there is z n + 1 such that | z - x 0 | 2 ρ and | G B ρ ( z ) | = | E | . Define G ~ = G B ρ ( z ) . Then by the minimality of E and Proposition 3.1 it holds

P ( E ) P ( G ~ ) + C h | Δ G ~ E | .

Arguing as in [9], one then deduces

P ( G ~ ) - P ( G ) n ( B ρ ( z ) G ) - n ( G B ρ ( z ) ) C ρ | B ρ ( z ) G | C h | Δ G ~ E |

and the claim follows as | Δ G ~ E | 2 | Δ G E | . ∎

The first technical result which I need is the classical density estimate which can be found, e.g., in [20].

Lemma 3.3.

Assume E R n + 1 is a set of finite perimeter with distributional mean curvature H E which satisfies H E L ( B 2 R ( x 0 ) ) Λ . Then for all x B R ( x 0 ) which are on the boundary of E and r min { R , Λ - 1 } it holds

P ( E , B r ( x ) ) c n r n

for a dimensional constant c n > 0 .

For a minimizer of (2.2) it holds the inverse of the isoperimetric inequality.

Lemma 3.4.

Let F R n + 1 be an open and bounded set of finite perimeter, fix h ( 0 , 1 ) and let E be a minimizer of (2.2). Then for all x E and r C 0 h it holds

P ( E , B r ( x ) ) C r min { | E B 2 r ( x ) | , | B 2 r ( x ) E | }

for a constant which depends on the dimension and on C 0 > 0 . In, particular it holds P ( E , B r ( x ) ) C r n .

Proof.

Fix h ( 0 , 1 ) , x and r > 0 as in the claim and without loss of generality assume that x = 0 . One may also consider only the case | E B 2 r | | B 2 r E | as the other case is similar. In particular, it holds | E B 2 r | 1 2 | B 2 r | . Since

0 2 r n ( B ρ E ) = | E B 2 r | ,

there is ρ ( r , 2 r ) such that

(3.2) n ( B ρ E ) C n | E B 2 r | r and | B ρ | 2 3 | B 2 r | .

Consider first the set E 1 = E B ¯ ρ . In order to have a competing set with the volume of E, define ρ ~ ρ to be a radius such that | B ρ ~ | = | E B ρ | and define E 2 = E 1 B ρ ~ . Then it holds by construction that | E 2 | = | E | , ρ ~ < ρ and

(3.3) n ( B ρ ~ ) = c n | B ρ ~ | n n + 1 = c n | E B ρ | n n + 1 c n | E B 2 r | n n + 1 C n | E B 2 r | r .

By the minimality of E we have

P ( E ) + 1 h E d ¯ F d x P ( E 2 ) + 1 h E 2 d ¯ F d x .

Estimate the perimeter of E 2 using (3.2) and (3.3) as

P ( E 2 ) P ( E , n + 1 B ρ ) + n ( B ρ E ) + n ( B ρ ~ ) P ( E , n + 1 B ρ ) + C n | E B 2 r | r .

Use then Δ E E 2 B 2 r , | E | = | E 2 | and the fact that the signed distance function is 1-Lipschitz to estimate

| E d ¯ F d x - E 2 d ¯ F | 4 r | E B ρ | 4 r | E B 2 r | .

Therefore one obtains by combining the three above inequalities and r C 0 h

P ( E , B r ) P ( E , B ρ ) C n | E B 2 r | r + 4 r h | E B 2 r | C | E B 2 r | r .

By Lemma 3.3 and Lemma 3.4 it is clear that for Proposition 3.1 it is crucial to prove the curvature estimate H E L C h . The next lemma is a step towards this.

Lemma 3.5.

Let F , E and h be as in Proposition 3.1. Then it holds

H E L 2 ( * E ) C 1 h ,

where the constant C 1 depends on the dimension and on | F | and P ( F ) .

Proof.

The proof relies on the second variation inequality in Proposition 2.2. I would like to point out that in the case of the mean curvature flow, when there is no volume constraint, the proof is considerable easier as one could choose constant function in (2.5). In the volume preserving case I will choose a cut-off function for a test function.

To this end, use first [11, Proposition 2.3] (see also [17, Lemma 2.1]) to find a point x 0 n + 1 and a radius r ( c , 1 ) , where c = c ( n , | F | , P ( F ) ) , such that

| E B r ( x 0 ) | = 1 2 | B r | .

Note that the minimality of E yields P ( E ) P ( F ) C . Moreover, by the isoperimetric inequality it holds P ( E ) c n | E | n n + 1 = c n | F | n n + 1 c . These estimates are used repeatedly from now on without mentioning. Without loss of generality assume that x 0 = 0 . Choose ρ < r such that | B r B ρ | = 1 4 | B r | . Note that then r - ρ c n > 0 and

(3.4) 3 4 | B ρ | | E B ρ | 1 4 | B ρ | .

Define first a cut-off function ζ C 0 1 ( n + 1 ) such that 0 ζ 1 , ζ = 1 in B ρ , ζ = 0 outside B r and | ζ | C n . Choose then φ = ζ - ζ ¯ , where ζ ¯ = * E ζ d n , as a test function in (2.5), use | d ¯ F , ν E | 1 and | φ | 1 , and obtain

(3.5) * E | B E | 2 ( ζ - ζ ¯ ) 2 d n * E | τ ζ | 2 d n + P ( E ) h .

Since H E = Trace ( B E ) , it holds point-wise on * E

(3.6) | B E | 2 H E 2 n .

Recall that 0 ζ 1 . Moreover, by the isoperimetric inequality and by (3.4) it holds P ( E , B ρ ) c n | E B ρ | n n + 1 c . Therefore ζ ¯ c for c = c ( n , | F | , P ( F ) ) . In particular, it holds | ζ ( x ) - ζ ¯ | c for x * E B r . Hence, we have by (3.5),

* E B r | H E | 2 d n C h P ( F ) .

We repeat the same argument by defining a cut-off function ζ C 1 ( n + 1 ) as ζ = 0 in B r , ζ = 1 outside B R and | ζ | C n , where R > r is such that | B R B r | = 1 4 | B r | . Using φ = ζ - ζ ¯ in (2.5) and arguing as above yields

* E B r | H E | 2 d n C h P ( F )

and the claim follows. ∎

The last lemma I need is a bound on the Lagrange multiplier in the Euler–Lagrange equation (2.4).

Lemma 3.6.

Let F , E and h be as in Proposition 3.1. Then for the Lagrange multiplier in (2.4), i.e.,

d ¯ F h = - H E + λ on  * E

it holds

| λ | C 2 h ,

where the constant C 2 depends on the dimension, on | F | and on P ( F ) .

Proof.

Let Λ 0 be such that | λ | = Λ h . Below all the constants depend on n , | F | and P ( F ) . I only treat the case when λ is positive as in the negative case the proof is similar. Define the set

Σ = { x * E : | H E ( x ) | < C ^ h } .

I claim that we may choose C ^ > 2 such that it depends on n , | F | , P ( F ) and on C 1 from Lemma 3.5 and it holds

(3.7) n ( Σ ) P ( E ) 2 .

Indeed, by Lemma 3.5 and by C ^ > 2 it holds

C ^ 2 h n ( * E Σ ) * E Σ H E 2 d n * E H E 2 d n C 1 2 h .

By choosing C ^ large enough one then obtains n ( * E Σ ) < P ( E ) 2 and (3.7) follows.

By the Besicovitch covering theorem one finds disjoint balls of radius h , denote them by { B h ( x i ) } i = 1 N , with x i Σ such that

(3.8) i = 1 N P ( E , B h ( x i ) ) c n n ( Σ ) c P ( E ) .

By the Euler–Lagrange equation (2.4) and by the definition of the set Σ it holds for all x * E B h ( x i ) , with x i Σ , that

(3.9) | H E ( x ) | | d ¯ F ( x ) h - λ | | d ¯ F ( x ) - d ¯ F ( x i ) | h + | d ¯ F ( x i ) h - λ | 1 h + | H E ( x i ) | 2 C ^ h .

Therefore by Lemma 3.3 it holds

P ( E , B h 2 ( x i ) ) P ( E , B h 2 C ^ ( x i ) ) c h n 2 .

On the other hand, applying Lemma 3.4 first with r = h 2 yields

| E B h ( x i ) | c h P ( E , B h 2 ( x i ) )

and then with r = h yields h n 2 c P ( E , B h ( x i ) ) . In conclusion, it holds

(3.10) | E B h ( x i ) | c h P ( E , B h ( x i ) )

for all balls in the cover.

The minimality of E implies

(3.11) 1 h E F d ¯ F d x P ( E ) + 1 h Δ E F | d ¯ F | d x P ( F ) .

Note that by (3.9), by λ = Λ h and by the Euler–Lagrange equation (2.4) we have for all x B h ( x i ) that

d ¯ F ( x ) λ h - | H ( x ) | h ( Λ - 2 C ^ ) h .

Therefore either Λ 4 C ^ , in which case the claim follows trivially, or

d ¯ F ( x ) Λ 2 h .

I assume the latter and show that also in this case Λ is bounded. Indeed, by the above discussion the balls B h ( x i ) are in the exterior of F. Therefore we estimate by (3.10) and by (3.8) that

1 h E F d ¯ F d x 1 h i = 1 N E B h ( x i ) d ¯ F d x 1 h i = 1 N ( Λ 2 h | E B h ( x i ) | ) c Λ i = 1 N P ( E , B h ( x i ) ) c Λ n ( Σ ) c Λ P ( E ) .

Since P ( E ) c n | E | n n + 1 = c n | F | n n + 1 , the above and (3.11) gives a bound for Λ and the claim follows. ∎

Here is the proof of the density estimate.

Proof of Proposition 3.1.

By Lemma 3.3, Lemma 3.4, Lemma 3.6 and by the Euler–Lagrange equation (2.4) it is enough to prove

(3.12) d ¯ F L ( * E ) C h .

Argue by contradiction and assume that there is x 0 * E such that

| d ¯ F ( x 0 ) | = Λ h

for large Λ 1 . Without loss of generality assume that x 0 = 0 and consider only the case d ¯ F ( x 0 ) > 0 as the case d ¯ F ( x 0 ) < 0 is similar.

By the Euler–Lagrange equation (2.4), by Lemma 3.3 and Lemma 3.6 it holds for r 0 = h 2 Λ that

(3.13) P ( E , B r 0 ) c r 0 n .

Define radii r k = k h + r 0 for k = 0 , 1 , , n + 1 , where n + 1 is the dimension of the ambient space. For every k = 0 , 1 , , n + 1 choose a cut-off function ζ k C 0 1 ( n + 1 ) such that 0 ζ 1 , ζ k = 1 in B r k , ζ k + 1 = 0 in n + 1 B r k + 1 and | ζ k | 2 h . Choose

φ = ζ k - ζ ¯ k ,

where ζ ¯ k = * E ζ k d n , as a test function in the second variation condition (2.5), use | d ¯ F , ν E | 1 and (3.6), and obtain

(3.14) 1 n * E | H E | 2 ( ζ k - ζ ¯ k ) 2 d n 1 h * E ( ζ k - ζ ¯ k ) 2 d n + * E | τ ζ k | 2 d n .

Since ζ k = 0 outside B r k + 1 , one may estimate

* E | τ ζ k | 2 d n 4 h P ( E , B r k + 1 ) .

Moreover, since 0 ζ k 1 , it holds

* E ( ζ k - ζ ¯ k ) 2 d n = * E ζ k 2 - ζ ¯ k 2 d n * E ζ k 2 d n P ( E , B r k + 1 ) .

When Λ is large enough, for all x * E B r k + 1 and all k n + 1 it holds d ¯ F ( x ) Λ 2 h . By the Euler–Lagrange equation (2.4), by d ¯ F ( x ) Λ 2 h , and by Lemma 3.6 it holds

| H E | Λ 4 h on  * E B r k

when Λ is large. Therefore, by using the notation from Lemma 3.6, it holds

* E B r k + 1 { x * E : | H E ( x ) | C ^ h } = * E Σ

when Λ is large. Then one may deduce from (3.7) that P ( E , B r k + 1 ) 1 2 P ( E ) . This yields

0 ζ ¯ k 1 2 .

Therefore it holds

1 n * E | H E | 2 ( ζ k - ζ ¯ k ) 2 d n 1 n * E B r k | H E | 2 ( 1 - ζ ¯ k ) 2 d n c n Λ 2 h P ( E , B r k ) .

Combining the three above estimates with (3.14) yields

c n Λ 2 h P ( E , B r k ) 5 h P ( E , B r k + 1 ) .

For Λ large enough this implies

(3.15) Λ P ( E , B r k ) P ( E , B r k + 1 ) .

Use (3.15) ( n + 1 ) -times from k = 0 to k = n , use then (3.13) and recall that r 0 = h 2 Λ and obtain finally that

(3.16) P ( E , B r n + 1 ) Λ n + 1 P ( E , B r 0 ) c Λ n + 1 r 0 n = c Λ n + 1 ( h 2 Λ ) n = c Λ h n 2 .

But now since r n + 1 = ( n + 1 ) h + r 0 2 ( n + 1 ) h , one obtains from Lemma 3.4 with r = 2 ( n + 1 ) h that

P ( E , B r n + 1 ) P ( E , B r ) C h n 2 ,

which contradicts (3.16) when Λ is large. ∎

4 Existence of the flat flow

Now that the density estimates are proven, the proof of Theorem 1 follows from the arguments from [14, 19] without major changes. In this section I consider the approximative flat flow ( E t h ) t 0 and the associated sequence ( E k h ) k 0 as in Definition 2.1, starting from an open and bounded set of finite perimeter E 0 . The proof for the following “interpolation” result can be found in [14, Lemma 1.5].

Lemma 4.1.

Let ( E t h ) t 0 be an approximative flat flow starting from E 0 and fix h ( 0 , 1 ) and t > h . Then for all l h it holds

| Δ E t h E t - h h | C ( l P ( E t - h h ) + 1 l Δ E t h E t - h h | d ¯ E t - h h | d x ) .

By the regularity result stated in Proposition 2.2, the Euler–Lagrange equation

(4.1) d ¯ E t - h h h = - H E t h + λ t , h

holds point-wise on * E t h and in a distributional sense on E t h . Here λ t , h is the Lagrange multiplier. Using the minimality of E t h against the previous set E t - h h , one obtains the important inequality

(4.2) P ( E t h ) + 1 h Δ E t h E t - h h | d ¯ E t - h h | d x P ( E t - h h ) .

Using (4.2) and the argument [14, Lemma 2.1] or [19, Lemma 3.6], one obtains the following dissipation inequality.

Lemma 4.2.

Let ( E t h ) t 0 be an approximative flat flow starting from E 0 and fix h ( 0 , 1 ) . Then for all T 2 > T 1 h it holds

T 1 T 2 H E t h - λ t , h L 2 ( * E t h ) 2 d t C ( P ( E T 1 - h ) - P ( E T 2 ) ) .

Moreover, it holds

T 1 T 2 ( H E t h L 2 ( * E t h ) 2 + λ t , h 2 ) d t C ( 1 + T 2 - T 1 ) .

The constant depends on the dimension, | E 0 | and P ( E 0 ) .

Proof.

I will only sketch the proof. Let ( E k h ) be the sequence of sets associated with ( E t h ) t 0 . For l with 2 l 2 C h - 1 2 set

K ( l ) = { x n + 1 : 2 l h < | d ¯ E t - h h ( x ) | 2 l + 1 h } .

Here C is such that | d ¯ E t - h h | C h on E t h . Proposition 3.1 yields that for every x E t h ,

| E t h B 2 l h ( x ) | c ( 2 l h ) n + 1 and n ( E t h B 2 l h ( x ) ) C ( 2 l h ) n .

Therefore for all x E t h K ( l ) it holds

B 2 l h ( x ) Δ E t h E t - h h | d ¯ E t - h h | d x c ( 2 l h ) n + 2 and B 2 l h ( x ) E t h d ¯ E t - h h 2 d n C ( 2 l h ) n + 2 .

Combining these two yields

B 2 l h ( x ) E t h d ¯ E t - h h 2 d n C B 2 l h ( x ) Δ E t h E t - h h | d ¯ E t - h h | d x .

By applying the Besicovitch covering theorem and summing over l (see [19, Lemma 3.6] for details) yields

E t h d ¯ E t - h h 2 d n Δ E t h E t - h h | d ¯ E t - h h | d x .

This together with (4.1) and (4.2) implies

h E t h ( H E t h - λ t , h ) 2 d n C ( P ( E t - h h ) - P ( E t h ) ) .

The first inequality then follows by iterating the above.

By [11, Lemma 2.4] it holds

| λ t , h | C ( 1 + H E t h - λ t , h L 1 ( * E t h ) )

for a constant that depends on the dimension, on | E 0 | and P ( E 0 ) . Note that then

H E t h L 2 2 + λ t , h 2 C ( 1 + H E t h - λ t , h L 2 2 ) .

Therefore by the first inequality one obtains

T 1 T 2 ( H E t h L 2 2 + λ t , h 2 ) d t C T 1 T 2 ( 1 + H E t h - λ t , h L 2 2 ) d t C ( 1 + T 2 - T 1 ) .

The third lemma we need is a quantitative bound on the diameter of the sets ( E t h ) , which is essentially the same as [19, Lemma 3.8].

Lemma 4.3.

Let ( E t h ) t 0 be an approximative flat flow starting from E 0 for h ( 0 , 1 ) . Then for all T > 0 there is R T , which depends on T, on the dimension and on the diameter of the initial set E 0 , such that E t h B R T for all t T .

Proof.

As in [19, Lemma 3.8] define r t for all t T as

r t := inf { r > 0 : E t h B r } .

Arguing as in [19, Lemma 3.8], one deduces that at the point y B r t E t h it holds H E t h ( y ) 0 and therefore by (4.1)

r t r t - h + h | λ t , h | .

Iterating this and using Lemma 4.2 yields

R T - R 0 0 T | λ t , h | d t 0 T ( 1 + λ t , h 2 ) d t C ( 1 + T ) .

Proof of Theorem 1.

Let ( E t h ) t 0 be an approximative flat flow starting from E 0 for h ( 0 , 1 ) and fix T 1 . Then by (4.2) it holds P ( E t h ) P ( E 0 ) and by Lemma 4.3 it holds E t h B R T for all t T . I claim that for 0 < t < s with s - t h it holds

(4.3) | Δ E t h E s h | C t - s .

Once (4.3) is obtained, then the convergence of a subsequence E t h k E t in measure follows as in [14, 19].

Let j , k be such that s [ j h , ( j + 1 ) h ) and t [ ( j + k ) h , ( j + k + 1 ) h ) . Then by applying Lemma 4.1 for l = h s - t and by (4.2) one obtains

| Δ E t h E s h | m = j j + k | Δ E m h h E ( m + 1 ) h h | C m = j j + k ( h s - t P ( E m h h ) + s - t h Δ E ( m + 1 ) h h E m h h | d ¯ E m h h | d x ) C m = j j + k h s - t P ( E 0 ) + C s - t m = j j + k ( P ( E m h h ) - P ( E ( m + 1 ) h h ) ) C k h s - t P ( E 0 ) + C s - t P ( E 0 ) .

Since k h 2 ( s - t ) , one obtains (4.3).

The proof of the consistency principle for C 1 , 1 -regular initial sets follows by using the arguments in [12]. The volume penalization is used only in [12, Lemma 3.2], but one may overcome this by using the lemma below. ∎

Lemma 4.4.

Let F R n + 1 be an open and bounded set which satisfies interior and exterior ball condition with radius r 0 > 0 and let E be a minimizer of (2.2). There are ρ 0 and h 0 with the property that if G is a set of finite perimeter such that

Δ G E B ρ ( x ) 𝒩 C 0 h ( F )

for ρ ρ 0 and h h 0 , where N C 0 h ( F ) = { x : dist ( x , F ) < C 0 h } , then it holds

P ( E ) P ( G ) + C ρ n + 1 .

Above the constant depends on the dimension, on r 0 , C 0 , | F | and P ( F ) .

Proof.

By approximation one may assume G to be smooth. Since F satisfies interior and exterior ball condition, by [12, Lemma 3.1] it holds

(4.4) max x Δ E F dist ( x , F ) C h

when h h 0 . As in the proof of Proposition 3.2 the set G may not have the same measure as E and one has to modify it to G ~ with | G ~ | = | E | . Assume that | G | < | E | . Since F satisfies interior ball condition with radius r 0 , there is y G such that B r 0 2 ( y ) G . By continuity there is z n + 1 such that | z - x | 2 ρ 0 and | G B r 0 2 ( z ) | = | E | when ρ 0 is small. Define G ~ = G B r 0 2 ( z ) . By the minimality of E, by (4.4) and by the assumption Δ G E B ρ ( x ) 𝒩 C 0 h ( F ) it holds

P ( E ) P ( G ~ ) + C | Δ G ~ E | P ( G ~ ) + C ρ n + 1 .

Arguing as in [9], one then deduces

P ( G ~ ) - P ( G ) n ( B r 0 2 ( z ) G ) - n ( G B r 0 2 ( z ) ) 2 ( n + 1 ) | B 1 | r 0 | B r 0 2 ( z ) G | C | Δ G ~ E | C ρ n + 1

and the claim follows. ∎

The paper concludes with Theorem 2. To this end, I recall the definition of a distributional solution of (1.1) from [14].

Definition 4.5.

A family of sets of finite perimeter ( E t ) t 0 is a distributional solution to (1.1) starting from E 0 n + 1 if the following holds:

  1. For almost every t > 0 the set E t has mean curvature H E t L 2 ( * E t ) in a distributional sense and for every T > 0 ,

    0 T H E t L 2 ( * E t ) 2 d t < .

  2. There exists v : n + 1 × ( 0 , ) with v L 2 ( 0 , T ; L 2 ( * E t ) ) such that for every ϕ C 0 1 ( n + 1 × [ 0 , ) ) it holds

    - 0 T * E t v ϕ d n d t = 0 T * E t ( H E t - H ¯ E t ) ϕ d n d t , 0 T E t t ϕ d x d t + E 0 ϕ ( , 0 ) d x = - 0 T * E t v ϕ d n d t .

Proof of Theorem 2.

The proof is exactly the same as [19, Theorem 2.3]. Note that Proposition 3.2 implies that the sets E t h are ( C h - 1 2 , c h ) -minimizers of the perimeter, i.e., for every F with Δ F E t h B c h ( x 0 ) it holds

P ( E t h ) P ( F ) + C h | Δ E t h F | .


Communicated by Frank Duzaar


Funding source: Academy of Finland

Award Identifier / Grant number: 314227

Award Identifier / Grant number: 347550

Funding statement: The author is supported by the Academy of Finland, grants no. 314227 and no. 347550.

References

[1] E. Acerbi, N. Fusco and M. Morini, Minimality via second variation for a nonlocal isoperimetric problem, Comm. Math. Phys. 322 (2013), no. 2, 515–557. 10.1007/s00220-013-1733-ySearch in Google Scholar

[2] F. Almgren, J. E. Taylor and L. Wang, Curvature-driven flows: a variational approach, SIAM J. Control Optim. 31 (1993), no. 2, 387–438. 10.1137/0331020Search in Google Scholar

[3] F. J. Almgren, Jr., Existence and regularity almost everywhere of solutions to elliptic variational problems with constraints, Mem. Amer. Math. Soc. 4 (1976), no. 165, 1–199. 10.1090/memo/0165Search in Google Scholar

[4] G. Antonelli, E. Pasqualetto and M. Pozzetta, Isoperimetric sets in spaces with lower bounds on the Ricci curvature, Nonlinear Anal. 220 (2022), Article ID 112839. 10.1016/j.na.2022.112839Search in Google Scholar

[5] G. Bellettini, V. Caselles, A. Chambolle and M. Novaga, The volume preserving crystalline mean curvature flow of convex sets in N , J. Math. Pures Appl. (9) 92 (2009), no. 5, 499–527. 10.1016/j.matpur.2009.05.016Search in Google Scholar

[6] L. Bronsard and B. Stoth, Volume-preserving mean curvature flow as a limit of a nonlocal Ginzburg–Landau equation, SIAM J. Math. Anal. 28 (1997), no. 4, 769–807. 10.1137/S0036141094279279Search in Google Scholar

[7] W. Carter, A. Roosen, J. Cahn and J. Taylor, Shape evolution by surface diffusion and surface attachment limited kinetics on completely faceted surfaces., Acta Metall. Materialia 43 (1995), 4309–4323. 10.1016/0956-7151(95)00134-HSearch in Google Scholar

[8] E. Gonzalez, U. Massari and I. Tamanini, On the regularity of boundaries of sets minimizing perimeter with a volume constraint, Indiana Univ. Math. J. 32 (1983), no. 1, 25–37. 10.1512/iumj.1983.32.32003Search in Google Scholar

[9] M. Grüter, Boundary regularity for solutions of a partitioning problem, Arch. Ration. Mech. Anal. 97 (1987), no. 3, 261–270. 10.1007/BF00250810Search in Google Scholar

[10] V. Julin, M. Morini, M. Ponsiglione and E. Spadaro, The asymptotics of the area-preserving mean curvature and the Mullins–Sekerka flow in two dimensions, Math. Ann. 387 (2023), no. 3–4, 1969–1999. 10.1007/s00208-022-02497-3Search in Google Scholar

[11] V. Julin and J. Niinikoski, Quantitative Alexandrov theorem and asymptotic behavior of the volume preserving mean curvature flow, Anal. PDE 16 (2023), no. 3, 679–710. 10.2140/apde.2023.16.679Search in Google Scholar

[12] V. Julin and J. Niinikoski, Consistency of the flat flow solution to the volume preserving mean curvature flow, Arch. Ration. Mech. Anal. 248 (2024), no. 1, Paper No. 1. 10.1007/s00205-023-01944-ySearch in Google Scholar

[13] T. Laux, Weak-strong uniqueness for volume-preserving mean curvature flow, Rev. Mat. Iberoam. 40 (2024), no. 1, 93–110. 10.4171/rmi/1395Search in Google Scholar

[14] S. Luckhaus and T. Sturzenhecker, Implicit time discretization for the mean curvature flow equation, Calc. Var. Partial Differential Equations 3 (1995), no. 2, 253–271. 10.1007/BF01205007Search in Google Scholar

[15] F. Maggi, Sets of Finite Perimeter and Geometric Variational Problems. An Introduction to Geometric Measure Theory, Cambridge Stud. Adv. Math. 135, Cambridge University, Cambridge, 2012. 10.1017/CBO9781139108133Search in Google Scholar

[16] P. W. Michor and D. Mumford, Riemannian geometries on spaces of plane curves, J. Eur. Math. Soc. (JEMS) 8 (2006), no. 1, 1–48. 10.4171/jems/37Search in Google Scholar

[17] M. Morini, M. Ponsiglione and E. Spadaro, Long time behavior of discrete volume preserving mean curvature flows, J. Reine Angew. Math. 784 (2022), 27–51. 10.1515/crelle-2021-0076Search in Google Scholar

[18] L. Mugnai and C. Seis, On the coarsening rates for attachment-limited kinetics, SIAM J. Math. Anal. 45 (2013), no. 1, 324–344. 10.1137/120865197Search in Google Scholar

[19] L. Mugnai, C. Seis and E. Spadaro, Global solutions to the volume-preserving mean-curvature flow, Calc. Var. Partial Differential Equations 55 (2016), no. 1, Article ID 18. 10.1007/s00526-015-0943-xSearch in Google Scholar

[20] L. Simon, Introduction to geometric measure theory, Tsinghua Lectures, 2014. Search in Google Scholar

Received: 2023-04-26
Accepted: 2024-03-10
Published Online: 2024-04-24
Published in Print: 2024-10-01

© 2024 Walter de Gruyter GmbH, Berlin/Boston

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Frontmatter
  2. Another proof of the existence of homothetic solitons of the inverse mean curvature flow
  3. A Weierstrass extremal field theory for the fractional Laplacian
  4. Minimizing movements for anisotropic and inhomogeneous mean curvature flows
  5. A singular Yamabe problem on manifolds with solid cones
  6. Uniqueness for volume-constraint local energy-minimizing sets in a half-space or a ball
  7. Limit of solutions for semilinear Hamilton–Jacobi equations with degenerate viscosity
  8. Monotonicity of entire solutions to reaction-diffusion equations involving fractional p-Laplacian
  9. Hierarchy structures in finite index CMC surfaces
  10. No breathers theorem for noncompact harmonic Ricci flows with asymptotically nonnegative Ricci curvature
  11. Wolff potentials and local behavior of solutions to elliptic problems with Orlicz growth and measure data
  12. Asymptotic analysis of single-slip crystal plasticity in the limit of vanishing thickness and rigid elasticity
  13. On the regularity of optimal potentials in control problems governed by elliptic equations
  14. Sobolev embeddings and distance functions
  15. Effective quasistatic evolution models for perfectly plastic plates with periodic microstructure: The limiting regimes
  16. On the area-preserving Willmore flow of small bubbles sliding on a domain’s boundary
  17. Sobolev contractivity of gradient flow maximal functions
  18. Discrete approximation of nonlocal-gradient energies
  19. Symmetry and monotonicity of singular solutions to p-Laplacian systems involving a first order term
  20. Flat flow solution to the mean curvature flow with volume constraint
Downloaded on 14.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/acv-2023-0047/html
Scroll to top button