Home On the existence of canonical multi-phase Brakke flows
Article Open Access

On the existence of canonical multi-phase Brakke flows

  • Salvatore Stuvard ORCID logo EMAIL logo and Yoshihiro Tonegawa
Published/Copyright: July 23, 2022

Abstract

This paper establishes the global-in-time existence of a multi-phase mean curvature flow, evolving from an arbitrary closed rectifiable initial datum, which is a Brakke flow and a BV solution at the same time. In particular, we prove the validity of an explicit identity concerning the change of volume of the evolving grains, showing that their boundaries move according to the generalized mean curvature vector of the Brakke flow. As a consequence of the results recently established in [J. Fischer, S. Hensel, T. Laux and T. M. Simon, The local structure of the energy landscape in multiphase mean curvature flow: Weak-strong uniqueness and stability of evolutions, preprint 2020, https://arxiv.org/abs/2003.05478], under suitable assumptions on the initial datum, such additional property resolves the non-uniqueness issue of Brakke flows.

MSC 2010: 53E10; 49Q15

1 Introduction

Arising as the gradient flow of the area functional, the mean curvature flow (henceforth abbreviated as MCF) is arguably the most fundamental geometric flow involving extrinsic curvatures. The unknown of MCF is a one-parameter family { Γ ( t ) } t 0 of surfaces in the Euclidean space (or in an ambient Riemannian manifold) such that the normal velocity of the flow equals the mean curvature vector at each point for every time t. The initial value problem for MCF starting with a smooth closed surface Γ ( 0 ) = Γ 0 is locally well-posed in time, until the appearance of singularities such as shrinking or neck pinching. Numerous frameworks of generalized solutions past singularities have been proposed: we mention, among others, the Brakke flows [6], level set flows [8, 12], BV solutions [25, 23] and L 2 flows [29, 5]. Existence of these possibly different generalized solutions to the MCF as well as their relations have been studied intensively in the past 40 years or so.

The aim of the present paper is to establish the global-in-time existence of a canonical multi-phase Brakke flow evolving from an arbitrary rectifiable initial datum. The attribute multi-phase here refers to the fact that the evolving surfaces are, in fact, boundaries of finitely many, but at least two, labelled open subsets of n + 1 (henceforth referred to also as grains or phases). The MCF evolution of such objects is strongly motivated by materials science, as it describes the motion and growth of crystallites in polycrystalline materials; see, e.g., [30]. While the literature concerning the two-phase MCF is rich, fewer works have been dedicated to the more general case of multi-phase MCF with at least three grains, despite its relevance in the modeling of physical processes governed by surface tension type energies.

In the present paper, we work with an arbitrary (but finite) number of grains. The solution we construct consists of two objects: the flow of the evolving grains and a Brakke flow, intertwined as follows. The Brakke flow – a measure-theoretic generalization of MCF, particularly suited to describe the evolution of surfaces through singularities (see Definition 2.1) – is essentially supported on the topological boundary of the grains, and it keeps track of multiplicities. Additionally, the mean curvature of the Brakke flow determines the distributional velocity at which the reduced boundary of each grain moves. As a result of the latter property, the change of volume of each grain between two instants of time can be recovered by integrating the mean curvature over the reduced boundary, a property certainly expected for a smooth MCF but quite non-trivial in a context where singularities and multiplicities occur. The attribute canonical refers to this very precise interplay between the Brakke flow and the evolution of the grains. Note that Brakke flows are non-unique in general due to the nature of the formulation, but the existence of the grains prevents redundant non-uniqueness such as sudden vanishing, for example. In the absence of higher multiplicities of the Brakke flow, we show that the collection of the grains (or, more precisely, the collection of their indicator functions) constitutes a BV solution of the MCF (see Definition 2.4). In this case, the grain boundaries are, in fact, a smooth MCF almost everywhere in space and time by Brakke’s partial regularity theory [6, 20, 38].

In certain instances, the additional BV characterization may lead to the uniqueness of the canonical Brakke flow. For instance, in ambient dimension n + 1 = 2 , the recent work of Fischer, Hensel, Laux and Simon [13] shows that, as soon as a strong solution of the network flow exists (see [13, Definition 16] for the definition of strong solution), any BV solution must coincide with it at least for all times until the first topology changes occur. One conclusion derived from [13] and the present paper is then that, if n = 1 and if a regular network flow starting from the given initial datum exists, then the canonical Brakke flow constructed in the present paper from the same initial datum necessarily coincides with that regular flow until the first topology changes; see [13, Theorem 19] for the precise statement (also see [17] for a similar uniqueness result in n = 2 ). In fact, the work [13] inspired the study carried out in the present paper.

In more precise and technical terms, the highlight of the main results of the present paper may be stated as follows (see the complete statements in Section 2.3).

Theorem A.

Let E 0 , 1 , , E 0 , N R n + 1 be mutually disjoint non-empty open sets with N 2 such that Γ 0 := R n + 1 i = 1 N E 0 , i is countably n-rectifiable. Assume that the n-dimensional Hausdorff measure of Γ 0 is finite or grows at most exponentially fast at infinity. Then there exist a Brakke flow { V t } t 0 (see Definition 2.1 below) as well as one-parameter families { E i ( t ) } t 0 ( i { 1 , , N } ) of open sets, with { E 1 ( t ) , , E N ( t ) } mutually disjoint for each t 0 , such that

V 0 = n Γ 0 , E i ( 0 ) = E 0 , i for every  i = 1 , , N ,

and satisfying the following properties, writing Γ ( t ) := R n + 1 i = 1 N E i ( t ) :

  1. n - 1 + δ ( Γ ( t ) spt V t ) = 0 for any δ > 0 and for a.e. t 0 .

  2. For each i = 1 , , N and for any arbitrary test function ϕ = ϕ ( x , t ) , we have, in the sense of distributions on [ 0 , ) , that

    (1.1) d d t E i ( t ) ϕ 𝑑 x = * E i ( t ) ϕ h ν i 𝑑 n + E i ( t ) ϕ t 𝑑 x .

  3. If the Brakke flow is locally a unit density flow, then, locally, we have

    n ( Γ ( t ) i = 1 N * E i ( t ) ) = 0 𝑎𝑛𝑑 V t = n i = 1 N * E i ( t )

    for a.e. t.

In the above statement, A B denotes the symmetric difference of two sets A and B, and spt V t is the support of the weight measure V t of the varifold V t . The symbol h denotes the generalized mean curvature vector of V t , and ν i is the outer unit normal vector field to the reduced boundary * E i ( t ) of E i ( t ) . Since Γ ( t ) = i = 1 N E i ( t ) for all t > 0 (see Theorem 2.11 (iii)), claim (i) shows that the support of the Brakke flow coincides – up to a lower-dimensional set – with the union of the topological boundaries of the grains for a.e. t > 0 . Claim (ii) states that each reduced boundary * E i ( t ) is a solution to the MCF in the integral sense specified in (1.1): that is, the generalized velocity of * E i ( t ) , defined as the distributional time derivative of the indicator function of E i ( t ) , is precisely h ν i n * E i ( t ) .

When integrated, formula (1.1) provides, as a byproduct, the change of the ( n + 1 ) -dimensional volume of each grain E i ( t ) in any bounded open set U:

(1.2) n + 1 ( U E i ( t 2 ) ) - n + 1 ( U E i ( t 1 ) ) = t 1 t 2 U * E i ( t ) h ν i 𝑑 n 𝑑 t .

We emphasize the following point in particular: formulae (1.1) and (1.2) hold true even if there is a “higher multiplicity portion” of V t on * E i ( t ) , or some “interior boundary” E i ( t ) * E i ( t ) . As far as the authors are aware of, for generalized MCF, (1.2) had never been established prior to the present paper, even for the two-phase MCF. It is worth mentioning that, in the latter context, the existence of a square-integrable generalized velocity advecting the common boundary of the two phases had been proved before, in the setting of sharp interface limits of minimizers of the Allen–Cahn action functional in [29], and of solutions to the Allen–Cahn equation in [16]. A fundamental new achievement of the present paper is, besides the multi-phase setting, the fact that the generalized velocity advecting each phase is precisely the dot product of the generalized mean curvature of the underlying varifold with the appropriate unit normal. Even though we have V 0 = n Γ 0 , and thus the unit density condition is satisfied at the initial time of the Brakke flow, it is not possible to exclude, in the very general framework under consideration, the occurrence of higher multiplicities at a later time. Despite all those possible singular behaviors, (1.1) and (1.2) are guaranteed time-globally. Claim (iii) states that, if the higher multiplicity of V t does not occur for n -a.e. x and a.e. t locally in space-time, then the reduced boundary measure and the Brakke flow may be identified with one another in that region, and we may say that the N-tuple χ = ( χ E 1 , , χ E N ) is a BV solution to MCF (see Definition 2.4). We can guarantee the existence of a unit density flow for some initial time interval [ 0 , T 0 ) if we additionally assume a suitable density ratio upper bound on Γ 0 (see Theorem 2.13). Such assumption would still allow for an initial datum Γ 0 which consists of a union of Lipschitz curves joined by triple junctions in n = 1 and Lipschitz bubble clusters with tetrahedral singularities in n = 2 , for example.

For general Brakke flows, there is no clear pathway leading from the characterization of Brakke flow, which consists of a variational inequality dictating an upper bound on the rate of change of the mass of the evolving surfaces, to formula (1.1), even under the unit density assumption. In fact, as mentioned already, by the partial regularity theorem for unit density Brakke flows [6, 20, 38], in this latter case it is known that Γ ( t ) is a C MCF in a space-time neighborhood of ( x 0 , t 0 ) for a.e. t 0 and for n -a.e. x 0 Γ ( t 0 ) : nonetheless, this alone is not sufficient to guarantee that χ = ( χ E 1 , , χ E N ) is a BV solution. The same remark goes for the opposite implication, i.e., from BV solution to Brakke flow. Since there is no known partial regularity theory for general BV solutions, these two notions appear far from being equivalent in any case. In the present paper, instead, the Brakke flow arises as the limit of a suitable time-discrete approximation scheme, analogous to that introduced by Kim and the second-named author in [21], and we prove (1.1) by showing that an analogous identity holds approximately true for the approximating flows, with vanishing errors in the limit. In order to gain enough control on the change of volumes in the approximation scheme and consequently obtain good estimates on the error terms, we will need to implement an appropriate modification to the construction of the time-discrete approximate flows devised in [21]: the details of such modification will be explained thoroughly later; see Section 3.1 and Section A.

Next, we discuss closely related works, particularly on the aspect of existence of generalized MCF. For two-phase MCF, the level-set method [8, 12] provides a general existence and uniqueness result even past the time after singularities appear. On the other hand, the level-set may develop a non-trivial interior, a phenomenon called “fattening”, due to the singular behavior of the MCF. Also, the uniqueness of the level-set solutions depends essentially on the maximum principle and it cannot handle general multi-phase MCF of more than two phases.

For the general multi-phase problem, it is natural to consider an initial datum Γ 0 with singularities to start with. For example, in dimension n = 1 , a typical Γ 0 in a three-phase problem is a union of curves meeting at triple junctions. In the parametric setting, Bronsard and Reitich [7] first showed the short-time existence of a unique solution for C 2 , α initial datum. Since then, there have been numerous studies (mostly for n = 1 , but also for higher dimensions [14, 9]), and we refer the reader to the survey [27] for the references on the parametric approaches. Due to the nature of the solutions and the need to heavily employ PDE techniques, these existence results do not extend beyond the time of topological changes. With a non-parametric approach and for the existence of MCF with regular triple junctions, one can adopt the elliptic regularization [19] for the class of flat chains with coefficients in a finite group; see [32].

Luckhaus and Sturzenhecker [25] introduced the formulation of BV solution of the two-phase MCF, which can be extended naturally to the multi-phase MCF. Their existence result is conditional, in the sense that a BV solution is shown to exist under the assumption that the time-discrete approximate solutions converge to their limit without loss of surface energy. Laux and Otto [23] proved that a sequence arising from the thresholding scheme of Merriman, Bence and Osher converges conditionally to a BV solution using the interpretation in terms of minimizing movements due to Esedoḡlu and Otto [10], again under an assumption similar to [25] (see also [24] for a similar convergence result of the parabolic Allen–Cahn system). The BV solution of [25] was partly motivated by the minimizing movements scheme of Almgren, Taylor and Wang [2], and the multi-phase version has been studied recently by Bellettini and Kholmatov [4].

On the side of Brakke flows, Ilmanen [18] proved the existence of a rectifiable Brakke flow arising as a limit of solutions to a parabolic Allen–Cahn equation, and the second-named author proved the integrality of the Brakke flow [37]. These results are for two-phase MCF, but the relation to the BV solution as formulated in [25] remained obscure. In a different but related problem – the Allen–Cahn action functional –, Mugnai and Röger [29] introduced a notion of L 2 flow to describe a weak formulation of MCF with additional L 2 forcing term for n = 1 , 2 . The existence of the L 2 flow depends on the result of Röger and Schätzle [31] which solved one of De Giorgi’s conjectures. The work [29] essentially contains the result that the limit phase boundary of the parabolic Allen–Cahn equation satisfies an analogous equation to (1.1) (see [29, Proposition 4.5]). The solution constructed in the present paper is, in fact, also an L 2 flow in the sense of Mugnai and Röger, and, even though the approach leading to (1.1) is different, some properties of generalized velocities of L 2 flows established in [29] are used in the present paper as well.

Finally, we mention again that a global-in-time existence result of a multi-phase Brakke flow which is equipped with moving grain boundaries was given by Kim and the second-named author in [21], reworking the pioneering paper by Brakke [6] within a different formulation. The grains in [21] move continuously with respect to the Lebesgue measure, but the problem concerning the validity of an exact identity involving the volume change was not addressed in there. Previous works by the authors of the present paper (see [35, 34]), in which certain Brakke flows are constructed with an approximation scheme analogous to that introduced in [21], could be reworked so that the additional conclusions concerning the interplay between the flow of the grains and the Brakke flow can be drawn in those contexts as well: in particular, it is possible to have the Brakke flow with prescribed boundary constructed in [35] satisfy formulae (1.1) and (1.2) (see Section 7.2).

2 Notation and main results

2.1 Basic notation

We shall use the same notation adopted in [21, Section 2]. In particular, the ambient space we will be working in is the Euclidean space n + 1 , and + will denote the interval [ 0 , ) . Coordinates ( x , t ) are set in the product space n + 1 × + , and t will be thought of and referred to as “time”. For a subset A of the Euclidean space, clos A , int A , A and conv A will denote the closure, interior, boundary and convex hull of A, respectively. If A n + 1 is (Borel) measurable, n + 1 ( A ) or | A | will denote the Lebesgue measure of A, whereas k ( A ) denotes the k-dimensional Hausdorff measure of A. When x n + 1 and r > 0 , then U r ( x ) and B r ( x ) denote the open ball and the closed ball centered at x with radius r, respectively. More generally, if k is an integer, then U r k ( x ) and B r k ( x ) will denote open and closed balls in k , respectively, and ω k := k ( U 1 k ( 0 ) ) .

A positive Radon measure μ on n + 1 (or in “space-time” n + 1 × + ) is always also regarded as a positive linear functional on the space C c ( n + 1 ) of continuous and compactly supported functions on n + 1 , with the pairing denoted by μ ( ϕ ) for ϕ C c ( n + 1 ) . The restriction of μ to a Borel set A is denoted by μ A , so that ( μ A ) ( E ) := μ ( A E ) for any E n + 1 . The support of μ is denoted by spt μ , and it is the closed set defined by

spt μ := { x n + 1 : μ ( B r ( x ) ) > 0  for every  r > 0 } .

The upper and lower k-dimensional densities of a Radon measure μ at x n + 1 are

Θ * k ( μ , x ) := lim sup r 0 + μ ( B r ( x ) ) ω k r k , Θ * k ( μ , x ) := lim inf r 0 + μ ( B r ( x ) ) ω k r k ,

respectively. If Θ * k ( μ , x ) = Θ * k ( μ , x ) , then the common value is denoted by Θ k ( μ , x ) and is called the k-dimensional density of μ at x. For 1 p , the space of p-integrable (resp. locally p-integrable) functions with respect to μ is denoted by L p ( μ ) (resp. L loc p ( μ ) ). If U n + 1 is an open set, L p ( n + 1 U ) and L loc p ( n + 1 U ) are simply written L p ( U ) and L loc p ( U ) . For a signed or vector-valued measure μ, | μ | denotes its total variation.

Given an open set U n + 1 , we say that a function f L 1 ( U ) has bounded variation in U, written f BV ( U ) , if

sup { U f div g d x : g C c 1 ( U ; n + 1 )  with  g C 0 1 } < .

If f BV ( U ) , then there exists an n + 1 -valued Radon measure on U, which we will call the measure derivative of f and denote by f , such that

U f div g d x = - U g d f for all  g C c 1 ( U ; n + 1 ) .

We say that f BV loc ( U ) if f BV ( U ) for all U U .

For a set E n + 1 , χ E is the characteristic (or indicator) function of E, defined by χ E ( x ) = 1 if x E , and χ E ( x ) = 0 otherwise. We say that E has locally finite perimeter in n + 1 if χ E BV loc ( n + 1 ) . When E is a set of locally finite perimeter, then the measure derivative χ E is the associated Gauss–Green measure, and its total variation | χ E | is the perimeter measure; by De Giorgi’s structure theorem, | χ E | = n * E , where * E is the reduced boundary of E, and χ E = - ν E | χ E | = - ν E n * E , where ν E is the outer pointing unit normal vector field to * E .

A subset Γ n + 1 is countably k-rectifiable if it is k -measurable and it admits a covering

Γ Z h f h ( k ) ,

where k ( Z ) = 0 and f h : k n + 1 is Lipschitz. A countably k-rectifiable set Γ is (locally) k -rectifiable if, moreover, k ( Γ ) is (locally) finite. Countably k-rectifiable sets Γ are characterized by the existence of approximate tangent planes k -almost everywhere (see [33, Theorem 11.6]). In other words, Γ is countably k-rectifiable if and only if there exists a positive function g L loc 1 ( k Γ ) such that the following holds: for k -a.e. x Γ there exists a k-dimensional linear subspace of n + 1 , denoted by T x Γ or Tan ( Γ , x ) and referred to as the (approximate) tangent plane to Γ at x, such that

g ( x + r ) k Γ - { x } r * g ( x ) k T x Γ in the sense of measures, as  r 0 + .

A (positive) measure μ on n + 1 is said to be k-rectifiable if there are a countably k-rectifiable set Γ and a positive function g L loc 1 ( k Γ ) such that μ = g k Γ . If μ is k-rectifiable, i.e., μ = g k Γ , then for μ-a.e. x the measures μ x , r defined, for r > 0 , by μ x , r ( A ) := r - k μ ( x + r A ) satisfy

μ x , r * g ( x ) k T x Γ as  r 0 + ,

that is, g ( x ) k T x Γ is the tangent measure of μ at x. In this case, the notation T x μ may be used interchangeably with T x Γ to denote the approximate tangent plane to the k-rectifiable measure μ at x for μ-a.e. x. More generally, given a Radon measure μ on n + 1 and x n + 1 , we say that μ has an approximate tangent k-plane at x if there exist g ( x ) ( 0 , ) and a k-dimensional linear subspace π n + 1 such that

μ x , r * g ( x ) k π as  r 0 + .

When this happens, the plane π is unique, and it is denoted by T x μ .

A k-dimensional varifold in n + 1 (see [1, 33]) is defined as a positive Radon measure V on the space n + 1 × 𝐆 ( n + 1 , k ) , where 𝐆 ( n + 1 , k ) is the Grassmannian of (unoriented) k-dimensional linear subspaces of n + 1 . If V is a k-varifold in n + 1 , we write V 𝐕 k ( n + 1 ) , and we let V and δ V denote its weight and first variation, respectively. When δ V is locally bounded and absolutely continuous with respect to V , we let h ( , V ) L loc 1 ( V ; n + 1 ) denote the generalized mean curvature vector of V, so that δ V = - h ( , V ) V in the sense of n + 1 -valued measures on n + 1 . If Γ is countably k-rectifiable, and θ L loc 1 ( k Γ ) is positive and integer-valued, we let 𝐯𝐚𝐫 ( Γ , θ ) denote the varifold 𝐯𝐚𝐫 ( Γ , θ ) = θ k Γ δ T Γ . When V admits a representation V = 𝐯𝐚𝐫 ( Γ , θ ) as above, we say that V is an integral k-varifold, and we write V 𝐈𝐕 k ( n + 1 ) . All above notions concerning measures (and varifolds) in the Euclidean space n + 1 can be immediately localized to open sets U n + 1 .

2.2 Three weak notions of MCF

As anticipated in Section 1, in the last few decades several alternative notions of weak solution to the MCF have been proposed. In this subsection, we briefly define and comment upon the three of interest in the present paper: Brakke flows, L 2 flows and BV flows. We begin with the notion of Brakke flow, introduced by Brakke in [6].

Definition 2.1 (Brakke flow).

Let 0 < T , and let U n + 1 be an open set. A k-dimensional (integral) Brakke flow in U is a one-parameter family of varifolds { V t } t [ 0 , T ) in U such that all of the following conditions hold:

  1. For a.e. t [ 0 , T ) , V t 𝐈𝐕 k ( U ) .

  2. For a.e. t [ 0 , T ) , δ V t is locally bounded and absolutely continuous with respect to V t .

  3. The generalized mean curvature h ( , V t ) (which exists for a.e. t by (ii)) satisfies h ( , V t ) L loc 2 ( V t ; n + 1 ) , and for every compact set K U and for every t < T it holds sup s [ 0 , t ] V s ( K ) < .

  4. For all 0 t 1 < t 2 < T and ϕ C c 1 ( U × [ 0 , T ) ; + ) , it holds

    V t 2 ( ϕ ( , t 2 ) ) - V t 1 ( ϕ ( , t 1 ) )
    (2.1) t 1 t 2 U { - ϕ ( x , t ) | h ( x , V t ) | 2 + ϕ ( x , t ) h ( x , V t ) + ϕ t ( x , t ) } d V t ( x ) 𝑑 t .

The inequality in (2.1) is typically referred to as Brakke’s inequality. It is not difficult to show that if { Γ ( t ) } t [ 0 , T ) is a flow of smooth submanifolds with mean curvature h ( x , t ) = h ( x , Γ ( t ) ) and normal velocity v ( x , t ) , then one has, for any 0 t 1 < t 2 < T and for every ϕ C c 1 ( U × [ t 1 , t 2 ] ) , the identity

(2.2) Γ ( t ) ϕ ( x , t ) d k | t = t 1 t 2 = t 1 t 2 Γ ( t ) { - ϕ v h + ϕ v + ϕ t } d k d t .

In particular, if { Γ ( t ) } t [ 0 , T ) is a smooth MCF, setting V t := 𝐯𝐚𝐫 ( Γ ( t ) , 1 ) defines a Brakke flow for which (2.1) is satisfied with an equality. Conversely, if V t = 𝐯𝐚𝐫 ( Γ ( t ) , 1 ) with Γ ( t ) smooth, then inequality (2.1) is sufficient to guarantee that { Γ ( t ) } t [ 0 , T ) is a classical solution to the MCF. For further details on the definition, the reader can consult the original work by Brakke in [6] or the more recent monograph [39]. Next, we define the notion of unit density flow.

Definition 2.2.

A Brakke flow is said to be a unit density flow in U × [ t 1 , t 2 ] if Θ k ( V t , x ) = 1 for V t -a.e. x U and 1 -a.e. t [ t 1 , t 2 ] . If U = n + 1 , we may simply say that { V t } t [ t 1 , t 2 ] is a unit density Brakke flow.

The following definition of L 2 flow has been given by Mugnai and Röger in [29].

Definition 2.3 ( L 2 flow).

Let 0 < T < , and let U n + 1 be an open and bounded set. A one-parameter family { V t } t [ 0 , T ) of varifolds in U is a k-dimensional L 2 flow if it satisfies (i)–(ii) in Definition 2.1 as well as the following conditions:

  1. The generalized mean curvature h ( , V t ) (which exists for a.e. t [ 0 , T ) by (ii)) satisfies

    h ( , V t ) L 2 ( V t ; n + 1 ) ,

    and d μ := d V t d t is a Radon measure on U × ( 0 , T ) .

  2. There exist a vector field v L 2 ( μ ; n + 1 ) and a positive constant C such that

    1. v ( x , t ) T x V t for μ-a.e. ( x , t ) U × ( 0 , T ) .

    2. For every ϕ C c 1 ( U × ( 0 , T ) ) , it holds

      (2.3) | 0 T U ϕ t ( x , t ) + ϕ ( x , t ) v ( x , t ) d V t ( x ) d t | C ϕ C 0 .

Any function v L 2 ( μ ; n + 1 ) satisfying (iv’) above is called a generalized velocity vector for the flow.

The definition can be easily motivated by considering (2.2) once again: the latter, indeed, implies (2.3) if both h and v are in L 2 ( μ ; n + 1 ) . Then L 2 flows can be described as flows of generalized surfaces with generalized mean curvature and velocity vectors in L 2 .

Finally, the following definition of BV flow is related to a motion of hypersurfaces which are the boundaries of a finite family of sets of locally finite perimeter; see [25, 23]. Given N 2 , we say that { E 1 , , E N } is an n + 1 -partition of n + 1 if E i n + 1 for every i, they are pairwise disjoint, and n + 1 ( n + 1 i = 1 N E i ) = 0 .

Definition 2.4 ( BV flow).

Suppose N 2 is an integer, and let 0 < T < . Then N one-parameter families { E i ( t ) } t [ 0 , T ) ( i = 1 , , N ) identify a BV solution for multi-phase MCF in n + 1 if all of the following assertions hold:

  1. For a.e. t [ 0 , T ) , { E 1 ( t ) , , E N ( t ) } is an n + 1 -partition of n + 1 , E i ( t ) is a set of locally finite perimeter, and, setting I i , j ( t ) := * E i ( t ) * E j ( t ) for i j ,

    (2.4) ess sup t [ 0 , T ) i , j = 1 , i j N n ( I i , j ( t ) ) < .

  2. There exist scalar functions v 1 , , v N such that

    (2.5) 0 T * E i ( t ) | v i ( x , t ) | 2 𝑑 n ( x ) 𝑑 t < for every  i ,

    and with the property that

    (2.6) E i ( t ) ϕ ( x , t ) d x | t = t 1 t 2 = t 1 t 2 E i ( t ) ϕ t ( x , t ) d x d t + t 1 t 2 * E i ( t ) ϕ ( x , t ) v i ( x , t ) d n ( x ) d t

    for a.e. 0 t 1 < t 2 < T and for all ϕ C c 1 ( n + 1 × [ 0 , T ) ) .

  3. Setting ν i ( x , t ) = ν E i ( t ) ( x ) for the outer unit normal to the reduced boundary of E i ( t ) at x, it holds

    v i ( , t ) ν i ( , t ) = v j ( , t ) ν j ( , t ) n -a.e. on  I i , j ( t ) ,  for a.e.  0 t < T .

  4. The functions v i further satisfy

    (2.7) i j 0 T I i , j ( t ) div g - ( ν i ν i ) g d n d t = - i j 0 T I i , j ( t ) v i ν i g 𝑑 n 𝑑 t

    for all vector fields g C c 1 ( n + 1 × [ 0 , T ] ; n + 1 ) .

  1. The following inequality holds for a.e. 0 t < T :

    (2.8) i , j = 1 , i j N n ( I i , j ( t ) ) + i , j = 1 , i j N 0 t I i , j ( s ) | v i ( x , s ) | 2 𝑑 n ( x ) 𝑑 s i , j = 1 , i j N n ( I i , j ( 0 ) ) .

Equation (2.6) characterizes v i ( , t ) as the velocity of * E i ( t ) , as anticipated in Section 1. By [26, Proposition 29.4], for each i = 1 , , N , we have

(2.9) n ( * E i ( t ) j = 1 , j i N I i , j ( t ) ) = 0 and n ( I i , j ( t ) I i , j ( t ) ) = 0    for  j j ,

and thus (2.7) formally characterizes v i ν i as the mean curvature vector of i = 1 N * E i ( t ) on I i , j ( t ) . Since

(2.10) n ( i = 1 N * E i ( t ) ) = 1 2 i , j = 1 , i j N n ( I i , j ( t ) ) ,

(v) corresponds precisely to the energy dissipation inequality for the hypersurface measures of the MCF.

2.3 Main results

Definition 2.5.

We will denote by Ω a function in C 2 ( n + 1 ) such that

0 < Ω ( x ) 1 , | Ω ( x ) | c 1 Ω ( x ) , 2 Ω ( x ) c 1 Ω ( x )

for all x n + 1 , where c 1 0 is a constant and 2 Ω ( x ) is the Hilbert–Schmidt norm of the Hessian matrix 2 Ω ( x ) .

The function Ω is introduced as a weight function to treat unbounded MCF, which may have infinite n-dimensional Hausdorff measure in n + 1 . A typical choice of Ω in this case can be Ω ( x ) = exp ( - 1 + | x | 2 ) with a suitable choice of c 1 . If one is interested in MCF with finite measure, one can choose Ω 1 , with c 1 = 0 in this case. With a function Ω as specified above, we consider an initial datum Γ 0 satisfying the following assumption.

Assumption 2.6.

Suppose that Γ 0 n + 1 is a closed, countably n-rectifiable set such that

(2.11) n Ω ( Γ 0 ) := Γ 0 Ω ( x ) 𝑑 n ( x ) < .

Moreover, assume that there are N 2 mutually disjoint non-empty open sets E 0 , 1 , , E 0 , N such that Γ 0 = n + 1 i = 1 N E 0 , i . In particular, { E 0 , 1 , , E 0 , N } is an n + 1 -partition of n + 1 .

By (2.11) and Ω > 0 , Γ 0 has locally finite n -measure, and thus Γ 0 has no interior points. In particular, we have

Γ 0 = i = 1 N E 0 , i .

For a given Γ 0 as in Assumption 2.11, the assignment of the partition n + 1 Γ 0 = i = 1 N E 0 , i is certainly non-unique. Each E 0 , i may not be connected, for example, and the different choice may result in different MCF. In general, E 0 , i has locally finite perimeter, and * E 0 , i spt | χ E 0 , i | E 0 , i for each i = 1 , , N .

The following Theorems 2.72.12 all hold under Assumption 2.11. First, we claim the existence of a Brakke flow starting with Γ 0 which is also an L 2 flow in any bounded domain of n + 1 , and whose generalized velocity vector coincides with the generalized mean curvature vector of the evolving varifolds.

Theorem 2.7.

There exists an n-dimensional Brakke flow { V t } t 0 in R n + 1 such that the following assertions hold:

  1. V 0 = n Γ 0 .

  2. If n ( i = 1 N ( E 0 , i * E 0 , i ) ) = 0 , then lim t 0 + V t = V 0 .

  3. V t ( Ω ) n Ω ( Γ 0 ) exp ( c 1 2 t / 2 ) and

    0 t n + 1 | h ( , V s ) | 2 Ω d V s 𝑑 s <

    for all t > 0 .

  4. If n ( Γ 0 ) < , and thus if one can choose Ω 1 in Assumption 2.11 , then

    V t ( n + 1 ) + 0 t n + 1 | h ( , V s ) | 2 d V s 𝑑 s n ( Γ 0 ) for all  t > 0 .

  5. For any 0 < T < and for any open and bounded set U n + 1 , the one-parameter family { V t U } t [ 0 , T ) (where we regard V t U as a varifold in U ) is an n -dimensional L 2 flow with generalized velocity vector v ( x , t ) = h ( x , V t ) on U × ( 0 , T ) .

Remark 2.8.

If n = 1 , the above Brakke flow satisfies the additional regularity properties obtained in [22]: for 1 -a.e. t > 0 , spt V t is locally the union of a finite number of W 2 , 2 curves meeting at junctions with angles of either 0, 60 or 120 degrees for N 3 , and only 0 degree (no transverse crossing) for N = 2 ; see [22, Theorems 2.2 and 2.3] for further details.

Definition 2.9.

With { V t } t 0 as in Theorem 2.7, suppose μ is the Radon measure on n + 1 × + given by d μ = d V t d t , and for t + define the closed set

( spt μ ) t := { x n + 1 : ( x , t ) spt μ } .

The following Theorem 2.10 is satisfied in general for Brakke flows, and thus it holds, in particular, for the flow produced in Theorem 2.7.

Theorem 2.10.

The Radon measure μ and the Brakke flow { V t } t 0 are related as follows:

  1. spt V t ( spt μ ) t and n ( B r ( spt μ ) t ) < for all t > 0 and r > 0 .

  2. n - 1 + δ ( ( spt μ ) t spt V t ) = 0 for every δ > 0 and for 1 -a.e. t + .

  3. spt V t is countably n -rectifiable for 1 -a.e. t + .

  4. V t = 𝐯𝐚𝐫 ( spt V t , θ t ) with θ t ( x ) = Θ n ( V t , x ) ( V t -a.e. x ) for 1 -a.e. t + .

The following theorem shows that, in addition to the Brakke flow of Theorem 2.7, there are evolving domains { E i ( t ) } t 0 starting from E 0 , i defining a (generalized) BV solution for multi-phase MCF in n + 1 .

Theorem 2.11.

For each i = 1 , , N , there exists a family of open sets { E i ( t ) } t 0 such that, setting

Γ ( t ) := n + 1 i = 1 N E i ( t ) ,

the following assertions hold:

  1. E i ( 0 ) = E 0 , i for i = 1 , , N .

  2. E 1 ( t ) , , E N ( t ) are pairwise disjoint for t + .

  3. ( spt μ ) t = Γ ( t ) = i = 1 N E i ( t ) for all t > 0 .

  4. For all t + , V t | χ E i ( t ) | for every i = 1 , , N , and 2 V t i = 1 N | χ E i ( t ) | .

  5. S ( i ) := { ( x , t ) : x E i ( t ) , t + } is open in n + 1 × + for i = 1 , , N .

  6. For every 0 < T < , the families { E i ( t ) } t [ 0 , T ) define a generalized BV solution for multi-phase MCF in the following sense: referring to Definition 2.4, (i”) holds when n Ω replaces n in ( 2.4 ); (ii”) holds when n Ω replaces n in ( 2.5 ); (iii”) holds. In fact, the scalar velocities v i are precisely

    v i ( x , t ) = h ( x , V t ) ν i ( x , t ) ,

    so that ( 2.6 ) reads as follows: for every i { 1 , , N } ,

    (2.12) E i ( t ) ϕ ( x , t ) d x | t = t 1 t 2 = t 1 t 2 ( * E i ( t ) ϕ h ν i d n + E i ( t ) ϕ t d x ) d t

    for any 0 t 1 < t 2 < and ϕ C c 1 ( n + 1 × + ) .

  7. If N 3 , then, for 1 -a.e. t + and n -a.e. x Γ ( t ) , θ t ( x ) = 1 implies x i = 1 N * E i ( t ) .

  8. If N = 2 , then, for 1 -a.e. t + and n -a.e. x Γ ( t ) ,

    θ t ( x ) = { odd integer for  x * E 1 ( t ) ( = * E 2 ( t ) ) , even integer for  x Γ ( t ) * E 1 ( t ) .

Notice that in (vi) we speak about generalized BV flow because we cannot guarantee, in general, the validity of (2.7). We will comment further on this point in Section 7.1. We anticipate, nonetheless, that further information on the flow, including the validity of (2.7), can be inferred when the Brakke flow is unit density. Even though the result follows immediately from Theorem 2.11 (vii), we state it separately.

Theorem 2.12.

Suppose that, for 0 t 1 < t 2 < and for an open set U, the Brakke flow in Theorem 2.7 is a unit density flow in U × ( t 1 , t 2 ) . Then, for L 1 -a.e. t ( t 1 , t 2 ) ,

(2.13) V t U = n U i = 1 N * E i ( t ) .

Furthermore, (2.7) holds true with v i = h ν i for all vector fields g C c 1 ( U × [ t 1 , t 2 ] ; R n + 1 ) .

If we assume the following additional conditions on Γ 0 , we can guarantee that the resulting Brakke flow { V t } t 0 as above is unit density for short time; see also [36, Theorem 2.2 (4)] for a similar statement valid in the context of mean curvature flows with a transport term.

Theorem 2.13.

Suppose that H n ( Γ 0 ) < and there exist r 0 > 0 and δ 0 > 0 such that

sup x n + 1 ,  0 < r < r 0 n ( Γ 0 B r ( x ) ) ω n r n < 2 - δ 0 .

Then there exists T 0 = T 0 ( n , r 0 , δ 0 , H n ( Γ 0 ) ) ( 0 , ) such that { V t } t [ 0 , T 0 ) in Theorem 2.7 is a unit density Brakke flow and { ( χ E 1 ( t ) , , χ E N ( t ) ) } t [ 0 , T 0 ) in Theorem 2.11 is a BV solution of MCF.

If we set “the maximal existence time” of unit density Brakke flow as

T ^ := sup { t 0 : 0 t V s ( { x : θ s ( x ) 2 } ) 𝑑 s = 0  and  V t 0 } ,

under the assumption of Theorem 2.13, either V t = 0 before T 0 , or we have T ^ T 0 and V t is a non-trivial unit density flow on [ 0 , T ^ ) and

{ ( χ E 1 ( t ) , , χ E N ( t ) ) } t [ 0 , T ^ )

is also a BV solution of MCF.

Finally, the validity of identity (2.12) allows to employ a simple argument to deduce a lower bound on the extinction time for V t , namely on the quantity

T * := sup { t + : V t 0 } .

Before stating the result, let us notice that if n ( Γ 0 ) < , then, by an argument using the relative isoperimetric inequality, there is only one i 0 { 1 , , N } such that | E 0 , i 0 | = , and we can let i 0 = N without loss of generality. We then set E ( t ) := i = 1 N - 1 E i ( t ) .

Theorem 2.14.

If H n ( Γ 0 ) < , then the extinction time T * satisfies

(2.14) T * 2 ( | E ( 0 ) | n ( Γ 0 ) ) 2 .

Observe that if Γ 0 is bounded, then T * must necessarily be finite, as one can see using the avoidance lemma of Brakke flows (see, for example, [19, 10.6 and 10.7]) and a comparison with a shrinking sphere. Estimate (2.14) was first proved by Giga and Yama-uchi in [15] for two-phase MCF for the solution by level-set method, and it is sharp. In contrast, as far as the authors are aware of, such sharp lower bound had never been established before in the context of multi-phase MCF. Theorem 2.14 will be proved in Section 7.3.

2.4 Outline of the proofs and plan of the paper

In what follows, we assume to have fixed the following:

  1. A function Ω as in Definition 2.5.

  2. A set Γ 0 and domains { E 0 , i } i = 1 N as in Assumption 2.11.

The strategy towards the proof of the existence of the Brakke flow { V t } t 0 and of the one-parameter families { E i ( t ) } t 0 of open sets for i { 1 , , N } from Theorems 2.7 and 2.11, respectively, is analogous to that employed by Kim and the second-named author in [21], with an important technical modification which is crucial to gain enough control on the change of volume of the grains and conclude identity (2.12). Let us explain this point in further detail. The scheme introduced in [21] can be roughly summarized as follows. One constructs, starting from { E 0 , i } i = 1 N , a sequence (indexed by j ) of piecewise constant-in-time flows of open partitions of n + 1 of N elements (see Section 2.5): more precisely, for each j the corresponding flow { j ( t ) } t [ 0 , j ] consists of constant open partitions j ( t ) = j , k = { E j , k , i } i = 1 N for t in intervals (epochs) ( ( k - 1 ) Δ t j , k Δ t j ] of length Δ t j 0 + . The goal is then to show that if the partition j , k + 1 is constructed from the partition j , k appropriately, then, along a suitable subsequence j , the (varifolds associated to) the boundaries j ( t ) converge to the desired Brakke flow V t .

The open partition j , k + 1 at a given epoch is constructed from the open partition j , k at the previous epoch by applying two operations, which we call steps. The first step is a small Lipschitz deformation of partitions with the effect of regularizing singularities by locally minimizing the area of the boundary of partitions at a suitably small scale; the second step consists of flowing the boundary of partitions by an appropriately defined approximate mean curvature vector, obtained by smoothing the surfaces’ first variation via convolution with a localized heat kernel.

The only difference between the scheme employed in [21] and the one devised in the present paper is in the choice of the Lipschitz deformations in the first step. While in [21] one only requires that the change of volumes of the grains due to Lipschitz deformation is small (for a certain smallness scale), in the present paper we shall require that the change of volume is small compared to the reduction in surface measure. We define such new class of volume-controlled Lipschitz deformations in Section 3.1, and then we claim that the construction of [21] can be carried over using volume-controlled deformations, as outlined in Section 3.2: by the results of [21], this completes the proof of Theorem 2.7 (i)–(iv) (with some additional remarks), Theorem 2.10 (i) and Theorem 2.11 (i)–(v). There are technical details for verifying that the modifications to the volume-controlled deformations do not pose any difficulties, and we defer it to Section A.

In Section 4, the main result is Theorem 2.11 (vi), namely identity (2.12). For that, we first prove that any Brakke flow is L 2 flow in Section 4.1, which is Theorem 2.7 (v). This fact seems to be first noted in [5], but we include the proof for completeness. Next, Section 4.2 is the main part of the proof, showing the existence of measure-theoretic velocities. The proof consists of calculating explicitly the rate of change of integrals of test functions across different epochs in the approximating flows. We will prove the validity of an approximate identity for the sets E j , i ( t ) , and the use of volume-controlled Lipschitz deformations will be crucial to conclude that the errors vanish in the limit as . Section 4.3 shows the existence of tangent space of μ on the reduced boundary of grains (seen as a set in space-time) and concludes the proof of Theorem 2.11 (vi).

In Section 5, we prove Theorem 2.10 (ii)–(iv), which, in particular, implies that the support spt V t of the evolving varifolds is n -equivalent to the boundary of partition Γ ( t ) = i = 1 N E i ( t ) . This result is used in Section 6. The argument we follow was suggested by Ilmanen in [18, Section 7.1], and we include it for the sake of completeness. It is based on the so-called clearing out lemma (see Lemma 5.1), a very robust result, which in turn follows from Huisken’s monotonicity formula, stating roughly that if the localized density of V t at some point y is too small at a scale r > 0 , then necessarily the point ( y , t + r 2 ) does not belong to the support of the space-time measure μ.

Section 6 contains the proofs of Theorem 2.11 (vii) and (viii). They are improved versions of the integrality theorem [21, Theorem 8.6] in that the behaviors of approximating flows and their grains are tracked more in detail. A characterization in [22, Section 4] of limiting behavior within a length scale of o ( 1 / j 2 ) , where measure minimizing property dominates, is essentially used. At the end, we describe the proofs of Theorems 2.12 and 2.13.

Finally, Section 7 contains some final remarks on the notion of generalized BV solutions introduced in Theorem 2.11 (vi) and on the construction of canonical multi-phase Brakke flows with fixed boundary conditions in the spirit of [35], as well as the proof of Theorem 2.14.

2.5 Further notation

We collect here some further notation and results which will be extensively used throughout the paper. We begin with the notion of open partitions of n + 1 of N elements and corresponding admissible maps.

Definition 2.15.

Let N 2 be an integer, and let Ω be as in Definition 2.5. A collection = { E 1 , , E N } of subsets E i n + 1 is called an Ω-finite open partition of N elements if the following conditions hold:

  1. E 1 , , E N are open and mutually disjoint.

  2. It holds

    n + 1 i = 1 N E i Ω ( x ) 𝑑 n ( x ) < .

  3. The set i = 1 N E i is countably n-rectifiable.

The set of all Ω-finite open partitions of N elements is denoted by 𝒪 𝒫 Ω N .

Note that it is allowed for some E i to be the empty set . Since Ω > 0 everywhere, property (ii) implies that the closed set n + 1 i = 1 N E i has locally finite n -measure, and thus no interior point. In particular, it holds

n + 1 i = 1 N E i = i = 1 N E i .

We define

:= i = 1 N E i ,

and with a slight abuse of notation we use the same symbol to denote the varifold

𝐯𝐚𝐫 ( i = 1 N E i , 1 ) 𝐈𝐕 n ( n + 1 ) ,

namely the unit density varifold induced from the countably n-rectifiable set .

Definition 2.16.

Given = { E i } i = 1 N 𝒪 𝒫 Ω N , a function f : n + 1 n + 1 is -admissible if it is Lipschitz and if it satisfies the following conditions, where we set E ~ i := int ( f ( E i ) ) for every i:

  1. E ~ 1 , , E ~ N are pairwise disjoint.

  2. It holds

    n + 1 i = 1 N E ~ i f ( i = 1 N E i ) .

  3. sup x n + 1 | f ( x ) - x | < .

The following lemma is [21, Lemma 4.4].

Lemma 2.17.

Let E = { E i } i = 1 N O P Ω N , and let f be E -admissible. Set E ~ := { E ~ i } i = 1 N , where E ~ i = int ( f ( E i ) ) . Then E ~ O P Ω N . We will call E ~ the push-forward of E through f, denoted by E ~ = : f E .

Next, we define a class of test functions with good properties. For j , we set

𝒜 j := { ϕ C 2 ( n + 1 ; + ) : ϕ ( x ) Ω ( x ) , | ϕ ( x ) | j ϕ ( x ) , 2 ϕ ( x ) j ϕ ( x )  for every  x n + 1 } .

Note that Ω 𝒜 j for all j c 1 .

Finally, we give the notion of smoothed mean curvature vector of a varifold. Let ψ C ( n + 1 ) be a radially symmetric cut-off function such that

ψ ( x ) = 1 for  | x | 1 2 ,
ψ ( x ) = 0 for  | x | 1 ,
0 ψ ( x ) 1 , | ψ ( x ) | 3 , 2 ψ ( x ) 9 for all  x n + 1 .

Then, for every ε ( 0 , 1 ) , we define

Φ ^ ε ( x ) := 1 ( 2 π ε 2 ) n + 1 2 exp ( - | x | 2 2 ε 2 ) , Φ ε ( x ) := c ( ε ) ψ ( x ) Φ ^ ε ( x ) ,

where 1 < c ( ε ) c ( n ) is a normalization constant chosen in such a way that

n + 1 Φ ε ( x ) 𝑑 x = 1 .

We shall call Φ ε the smoothing kernel at scale ε. For a varifold V 𝐕 n ( n + 1 ) , we then define the ε-smoothed mean curvature vector of V to be the vector field h ε ( , V ) C ( n + 1 ; n + 1 ) defined by

h ε ( , V ) := - Φ ε ( Φ ε δ V Φ ε V + ε Ω - 1 ) .

In the above formula, Φ ε V is the measure on n + 1 defined by

( Φ ε V ) ( ϕ ) := V ( Φ ε ϕ ) = n + 1 n + 1 Φ ε ( x - y ) ϕ ( y ) 𝑑 y d V ( x )

for all ϕ C c ( n + 1 ) , identified with the smooth function

( Φ ε V ) ( x ) := n + 1 Φ ε ( y - x ) d V ( y )

by means of the identity

( Φ ε V ) ( ϕ ) = Φ ε V , ϕ L 2 ( n + 1 ) .

Analogously, Φ ε δ V is the n + 1 -valued measure on n + 1 defined by

( Φ ε δ V ) ( g ) := δ V ( Φ ε g ) = n + 1 g ( y ) n + 1 × 𝐆 ( n + 1 , n ) S ( Φ ε ( x - y ) ) 𝑑 V ( x , S ) 𝑑 y

for all g C c ( n + 1 ; n + 1 ) , identified with the smooth vector field

( Φ ε δ V ) ( x ) := n + 1 × 𝐆 ( n + 1 , n ) S ( Φ ε ( y - x ) ) 𝑑 V ( y , S ) .

We state the following lemma concerning the smoothed mean curvature vector to be used in the subsequent sections. For the proof, the reader can consult [21, Lemma 5.1].

Lemma 2.18.

For every M > 0 , there exists a constant ε 1 ( 0 , 1 ) , depending only on n, c 1 and M, such that the following holds. Let V V n ( R n + 1 ) be an n-dimensional varifold in R n + 1 such that V ( Ω ) M , and for every ε ( 0 , ε 1 ) let h ε ( , V ) be its smoothed mean curvature vector. Then

(2.15) | h ε ( x , V ) | 2 ε - 2 ,
(2.16) h ε ( x , V ) 2 ε - 4 .

3 Existence of multi-phase Brakke flow

3.1 Volume-controlled Lipschitz deformations

In this subsection, we introduce the modified class of Lipschitz deformations used in the present paper to gain improved control on the volume change of partitions. For 𝒪 𝒫 Ω N and j , the class is denoted by 𝐄 vc ( , j ) (where vc indicates “volume controlled”), and it replaces the class 𝐄 ( , j ) defined in [21, Definition 4.8].

Definition 3.1.

For = { E i } i = 1 N 𝒪 𝒫 Ω N and c 1 j , define 𝐄 vc ( , j ) to be the set of all -admissible functions f : n + 1 n + 1 such that, writing { E ~ i } i = 1 N := f , the following conditions hold:

  1. | f ( x ) - x | 1 / j 2 for all x n + 1 .

  2. It holds

    n + 1 ( E ~ i E i ) { ( Ω ) - f ( Ω ) } j for all  i = 1 , , N .

  3. f ( ϕ ) ( ϕ ) for all ϕ 𝒜 j .

The difference between the above class 𝐄 vc ( , j ) and the class 𝐄 ( , j ) in [21, Definition 4.8] lies in condition (ii), which in [21] was simply n + 1 ( E ~ i E i ) 1 / j . Since Ω 𝒜 j for all j c 1 , (iii) implies

( Ω ) f ( Ω ) ,

and thus the right-hand side of (ii) is non-negative. In particular, the identity map f ( x ) = x belongs to 𝐄 vc ( , j ) for j c 1 . We similarly modify the definitions of Δ j ( Ω ) , 𝐄 ( , C , j ) and Δ j ( C ) in [21, (4.11)–(4.13)] by introducing their volume-controlled counterparts as follows.

Definition 3.2.

For 𝒪 𝒫 Ω N and c 1 j , define

(3.1) Δ j vc ( Ω ) := inf f 𝐄 vc ( , j ) { f ( Ω ) - ( Ω ) } 0 ,

and for a compact set C n + 1 ,

𝐄 vc ( , C , j ) := { f 𝐄 vc ( , j ) : { x : f ( x ) x } { f ( x ) : f ( x ) x } C } ,
(3.2) Δ j vc ( C ) := inf f 𝐄 vc ( , C , j ) { f ( C ) - ( C ) } .

Even with this modification, the claims in [21, Section 4] remain essentially the same.

Lemma 3.3.

For compact sets C C ~ , we have

Δ j vc ( C ~ ) Δ j vc ( C )

and

Δ j vc ( Ω ) ( max C Ω ) { Δ j vc ( C ) + ( 1 - exp ( - c 1 diam C ) ) ( C ) } .

Lemma 3.4.

Suppose that { C l } l = 1 is a sequence of compact sets which are mutually disjoint and suppose that C is a compact set with l = 1 C l C . Then

(3.3) Δ j vc ( C ) l = 1 Δ j vc ( C l ) .

We remark that the corresponding statement in [21, Lemma 4.11] contains the additional assumption n + 1 ( C ) < 1 / j , which we do not need to assume in the present setting. For completeness, we provide the proof.

Proof.

By (3.2),

Δ j vc ( C ) - ( C ) > -

for any compact set C. Let m and ε ( 0 , 1 ) be arbitrary. For all l m , choose f l 𝐄 vc ( , C l , j ) such that

Δ j vc ( C l ) + ε ( f l ) ( C l ) - ( C l ) .

We define a map f : n + 1 n + 1 by setting f C l ( x ) = ( f l ) C l ( x ) for l = 1 , , m , and f n + 1 l = 1 m C l ( x ) = x . The -admissibility of f follows from that of f l and from the fact that { C l } are mutually disjoint. To prove f 𝐄 vc ( , j ) , we need to check Definition 3.1 (i)–(iii) and (i) follows immediately. Writing { E ~ i } i = 1 N := f , we have E ~ i E i = l = 1 m C l ( E ~ i E i ) , and

n + 1 ( E ~ i E i ) = l = 1 m n + 1 ( C l ( E ~ i E i ) )
l = 1 m { ( Ω ) - ( f l ) ( Ω ) } j
= l = 1 m { C l ( Ω ) - ( f l ) C l ( Ω ) } j
= { ( Ω ) - f ( Ω ) } j ,

where we used (ii) for f l 𝐄 vc ( , C l , j ) to conclude (ii) for f. Condition (iii) can be checked similarly. These show that f 𝐄 vc ( , j ) , and in fact f 𝐄 vc ( , C , j ) by construction. By the definition of Δ j vc ( C ) , we have

Δ j vc ( C ) f ( C ) - ( C )
= l = 1 m ( ( f l ) ( C l ) - ( C l ) )
m ε + l = 1 m Δ j vc ( C l ) .

By letting ε 0 + first and then m , we obtain (3.3). ∎

The following lemma is similar to [21, Lemma 4.12] with a minor change. The proof is identical.

Lemma 3.5.

Suppose that E = { E i } i = 1 N O P Ω N , c 1 j N , C is a compact subset of R n + 1 , f : R n + 1 R n + 1 is an E -admissible function such that, writing { E ~ i } i = 1 N := f E , the following conditions hold:

  1. { x : f ( x ) x } { f ( x ) : f ( x ) x } C .

  2. | f ( x ) - x | 1 / j 2 for all x n + 1 .

  3. n + 1 ( E ~ i E i ) { ( Ω ) - f ( Ω ) } / j for all i = 1 , , N .

  4. f ( C ) exp ( - j diam C ) ( C ) .

Then we have f E vc ( E , C , j ) .

3.2 The constructive scheme

In this subsection, we provide the detailed construction of the sequence of piecewise constant-in-time approximating flows of open partitions leading to the existence result of a multi-phase Brakke flow. We let the weight function Ω be as in Definition 2.5, and we consider an initial rectifiable set Γ 0 with a corresponding Ω-finite open partition of N elements 0 as in Assumption 2.11. For every natural number j c 1 , and for times t [ 0 , j ] , we define open partitions j ( t ) = { E j , 1 ( t ) , , E j , N ( t ) } according to the following rule:

(3.4) j ( 0 ) = 0 ,
(3.5) j ( t ) = j , k for all  t ( ( k - 1 ) Δ t j , k Δ t j ] .

In (3.5), the epoch length is Δ t j = 2 - p j for some p j , and k { 1 , , j 2 p j } . For each k, the open partition j , k is obtained from the open partition j , k - 1 (with the convention j , 0 = 0 ) through successive modifications, encoded in the following two-step algorithm:

  1. First, one chooses f 1 𝐄 vc ( j , k - 1 , j ) with the property that

    ( f 1 ) j , k - 1 ( Ω ) - j , k - 1 ( Ω ) ( 1 - j - 5 ) Δ j vc j , k - 1 ( Ω ) ,

    and sets

    j , k * := ( f 1 ) ( j , k - 1 ) .

    Thus, in particular,

    E j , k , i * := int ( f 1 ( E j , k - 1 , i ) ) for every  i { 1 , , N } .

  2. Next, one defines the map

    f 2 ( x ) := x + Δ t j h ε j ( x , j , k * ) ,

    where ε j ( 0 , 1 ) and h ε j ( , j , k * ) is the ε j -smoothed mean curvature vector of the multiplicity one varifold on j , k * . Notice that f 2 is a diffeomorphism of n + 1 due to Lemma 2.18 as soon as Δ t j ε j 4 . We set

    j , k := ( f 2 ) j . k * ,

    and thus

    E j , k , i := f 2 ( E j , k , i * ) for every  i { 1 , , N } .

Notice that the scheme just defined differs from that adopted in [21] only in step (1) of the algorithm, where volume-controlled Lipschitz deformations are used. In spite of such modification, we claim that the proof from [21] can be essentially carried over to this framework, and thus leading to the following theorem.

Theorem 3.6.

There is a constant c 2 = c 2 ( n ) 1 with the following property. Let Ω, Γ 0 and E 0 O P Ω N be as in Definition 2.5 and Assumption 2.11. Then there exist

  1. a subsequence j of ,

  2. reals ε j ( 0 , j - 6 ) with lim ε j = 0 ,

  3. integers p j with Δ t j = 2 - p j ( 2 - 1 ε j c 2 , ε j c 2 ] ,

  4. a family { μ t } t + of Radon measures on n + 1 ,

  5. a family ( t ) = { E 1 ( t ) , , E N ( t ) } t 0 of open sets

such that the approximating flow of open partitions E j ( t ) defined by (3.4) and (3.5) satisfies for all T < ,

(3.6) lim sup sup t [ 0 , T ] j ( t ) ( Ω ) 0 ( Ω ) exp ( c 1 2 T 2 ) ,
(3.7) lim sup 0 T ( n + 1 | Φ ε j δ ( j ( t ) ) | 2 Ω Φ ε j j ( t ) + ε j Ω - 1 𝑑 x - 1 Δ t j Δ j vc j ( t ) ( Ω ) ) 𝑑 t < ,
(3.8) lim j 2 ( n + 1 ) Δ j vc j ( t ) ( Ω ) = 0 for a.e.  t + ,
(3.9) lim j ( t ) ( ϕ ) = μ t ( ϕ ) for all  ϕ C c ( n + 1 ) and any  t + ,
(3.10) χ E j , i ( t ) χ E i ( t ) in  L loc 1 ( n + 1 ) as  for every  i { 1 , , N } and for every  t + .

Furthermore, the following assertions hold:

  1. There exists a subset Z + with 1 ( Z ) = 0 such that, for every t + Z , μ t is integral: that is, there exists V t 𝐈𝐕 n ( n + 1 ) such that V t = μ t .

  2. If V t is defined to be an arbitrary varifold in 𝐕 n ( n + 1 ) with V t = μ t also for t Z , then the family { V t } t + is an n -dimensional Brakke flow in n + 1 satisfying the conclusions of Theorem 2.7 (i)(iv) and Theorem 2.10 (i).

  3. The flow of grains E i ( t ) satisfies the conclusions of Theorem 2.11 (i)(v).

The proofs of the claims contained in the statement of Theorem 3.6 are analogous to the corresponding ones outlined in [21], modulo a few technical modifications in consequence of the use of volume-controlled Lipschitz deformations. More precisely, the following assertions hold:

  1. The conclusions (3.6) to (3.9) are contained in [21, Proposition 6.1, Proposition 6.4].

  2. The existence of the flow of grains E i ( t ) , the convergence in (3.10), and the conclusion (iii) is [21, Theorem 3.5].

  3. The conclusion (i) is [21, Lemma 9.1].

  4. The conclusion (ii) is [21, Theorem 3.2, Proposition 3.4], except that Theorem 2.7 (ii) follows from the argument in [35, Proposition 6.10] and Theorem 2.7 (iv) follows from Brakke’s inequality and Theorem 2.7 (iii) with Ω = 1 : Fix ϕ C c ( n + 1 ; [ 0 , 1 ] ) with ϕ ( x ) = 1 on B 1 and ϕ ( x ) = 0 on n + 1 B 2 , and for k , set ϕ k ( x ) := ϕ ( x / k ) . Use ϕ k in (2.1) with t 1 = 0 and arbitrary t = t 2 > 0 and let k . In doing so,

    0 t | ϕ k | | h | d V s 𝑑 s ( 0 t | ϕ k | 2 d V s 𝑑 s ) 1 2 ( 0 t | h | 2 d V s 𝑑 s ) 1 2 0

    as k , due to | ϕ k | = | ϕ | / k and uniform bounds from Theorem 2.7 (iii). Then one can see that Theorem 2.7 (iv) holds true.

We will not repeat the proofs of the above results, but we will detail the aforementioned changes due to volume-controlled Lipschitz deformations in Section A. For the time being, we will assume the validity of Theorem 3.6, and we will proceed with the derivation of all remaining results claimed in Section 2.

4 BV flow: Proof of Theorem 2.11 (vi)

The main result of this section is the proof of Theorem 2.11 (vi), which we isolate as the following theorem.

Theorem 4.1.

The one-parameter family { E i ( t ) } t R + of open partitions defined in Theorem 3.6 is a generalized BV solution to multi-phase MCF with scalar velocities v i = h ν i . More precisely, for every i { 1 , , N } it holds

(4.1) E i ( t ) ϕ ( x , t ) d x | t = t 1 t 2 = t 1 t 2 E i ( t ) ϕ t ( x , t ) d x d t + t 1 t 2 * E i ( t ) ϕ ( x , t ) h ( x , V t ) ν E i ( t ) ( x ) d n ( x ) d t

for every 0 t 1 < t 2 < and all test functions ϕ C c 1 ( R n + 1 × R + ) .

Remark 4.2.

Introducing the notation χ i for the indicator function χ i ( x , t ) := χ E i ( t ) ( x ) , and recalling that χ E i ( t ) = - ν E i ( t ) n * E i ( t ) as n + 1 -valued Radon measures on n + 1 , we see that the identity in (4.1) can be rephrased as

n + 1 ϕ ( x , t ) χ i ( x , t ) d x | t = t 1 t 2 = t 1 t 2 n + 1 ϕ t ( x , t ) χ i ( x , t ) d x d t - t 1 t 2 n + 1 ϕ ( x , t ) h ( x , V t ) d χ E i ( t ) ( x ) d t .

This implies, in particular, that χ i BV loc ( n + 1 × + ) and its (measure) derivative is

χ i ( x , t ) = ( χ E i ( t ) ( x ) , - h ( x , V t ) χ E i ( t ) ( x ) ) d t
(4.2) = ( - ν i ( x , t ) , h ( x , V t ) ν i ( x , t ) ) d n * E i ( t ) d t ,

where = ( , t ) and the identity holds in the sense of ( n + 1 × ) -valued Radon measures on n + 1 × + .

The proof of Theorem 4.1 will be obtained in three main steps, which are the content of Sections 4.14.3, respectively. First, we check that Brakke flow implies L 2 flow and recall the important observation Corollary 4.4 due to Mugnai and Röger [29]. Second, we will prove that χ i BV loc ( n + 1 × + ) for each i { 1 , , N } ; see Proposition 4.5. Then we will characterize the measure-theoretic time-derivative of χ i to show that (4.2) holds.

4.1 Characterization as L 2 flow

In this subsection, we first prove Theorem 2.7 (iv), which we isolate as the following theorem.

Theorem 4.3.

Let { V t } t R + be the Brakke flow defined in Theorem 3.6. Then, for every 0 < T < and for every open and bounded U R n + 1 , the varifolds V t ( U × G ( n + 1 , n ) ) ( t [ 0 , T ) ) are an n-dimensional L 2 flow with generalized velocity vector v ( , t ) = h ( , V t ) .

Proof.

We verify the requirements of Definition 2.3. Conditions (i) and (ii) are satisfied by V t in n + 1 , and thus they are trivially satisfied when the varifolds are restricted to the open set U. Concerning (iii’), we have immediately that

0 T U | h ( , V t ) | 2 d V t 𝑑 t ( max clos U Ω - 1 ) 0 T n + 1 | h ( , V t ) | 2 Ω d V t 𝑑 t < ,

so that h ( , V t ) L 2 ( V t U ; n + 1 ) for a.e. t [ 0 , T ) . Analogously, d V t d t is a Radon measure on U × ( 0 , T ) , and in fact we have

| 0 T U ϕ ( x , t ) d V t ( x ) 𝑑 t | ϕ C 0 ( max clos U Ω - 1 ) n Ω ( Γ 0 ) 0 T exp ( c 1 2 t 2 ) 𝑑 t <

for every ϕ C c ( U × ( 0 , T ) ) . To complete the proof, we are then only left with checking that (2.3) holds with v ( , t ) = h ( , V t ) : with this choice, indeed, the condition in (d’1) is automatically satisfied due to the perpendicularity of mean curvature (see [6, Chapter 5]). Consider first a test function ϕ C c 1 ( n + 1 × ( 0 , T ] ; + ) . Brakke’s inequality (2.1) then implies that

0 T n + 1 ϕ t ( x , t ) + h ( x , V t ) ϕ ( x , t ) d V t ( x ) d t 0 T n + 1 ϕ ( x , t ) | h ( x , V t ) | 2 d V t ( x ) 𝑑 t 0 .

In particular, the assignment

ϕ L ϕ := 0 T n + 1 ϕ t ( x , t ) + h ( x , V t ) ϕ ( x , t ) d V t ( x ) d t

defines a positive linear functional on C c 1 ( n + 1 × ( 0 , T ] ) : hence, L is monotone, that is, L ϕ 1 L ϕ 2 whenever ϕ 1 ϕ 2 everywhere. For every ε ( 0 , 1 ) , let ψ U , ε C c 1 ( n + 1 × ( 0 , T ] ) be the cut-off function

ψ U , ε ( x , t ) = ψ U ( x ) ψ ε ( t )

defined according to the following prescriptions:

  1. 0 ψ U 1 and 0 ψ ε 1 .

  2. ψ U 1 on clos U , and ψ U 0 on dist ( x , clos U ) 1 .

  3. ψ ε 1 on [ ε , T ] , and ψ ε 0 on ( 0 , ε / 2 ] .

  4. ψ U C 0 C and ψ ε C 0 C / ε for a geometric constant C.

Let now ϕ C c 1 ( U × ( 0 , T ) ) be a non-zero function such that spt ( ϕ ) U × [ ε , T ] . For such a function, by the definition of ψ U , ε , it holds

- ψ U , ε ϕ ϕ C 0 ψ U , ε ,

so that the linearity and monotonicity of L yield

| L ϕ | L ψ U , ε ϕ C 0 .

Notice that for such a ϕ the space integration can be restricted to U (as the – continuous – derivatives of ϕ are necessarily zero on U c for every t), and thus the above argument shows that whenever ϕ C c 1 ( U × ( 0 , T ) ) is supported on U × [ ε , T ] , it holds

(4.3) | 0 T U ϕ t ( x , t ) + h ( x , V t ) ϕ ( x , t ) d V t ( x ) d t | L ψ U , ε ϕ C 0 .

We next proceed to estimate L ψ U , ε . We have, setting U 1 := { dist ( x , clos U ) < 1 } ,

L ψ U , ε C ε ε 2 ε V t ( U 1 ) 𝑑 t + 0 T n + 1 | h ( , V t ) | | ψ U | d V t 𝑑 t
C ( max clos U 1 Ω - 1 ) ( n Ω ( Γ 0 ) + ( T n Ω ( Γ 0 ) ) 1 2 ( 0 T n + 1 | h ( , V t ) | 2 Ω d V t 𝑑 t ) 1 2 ) .

We see then that sup ε > 0 L ψ U , ε < : in particular, thanks to (4.3), we conclude that inequality (2.3) holds with v ( , t ) = h ( , V t ) for every ϕ C c 1 ( U × ( 0 , T ) ) , and thus completing the proof. ∎

The following is a simple corollary of Theorem 4.3 and [29, Proposition 3.3].

Corollary 4.4.

Let { V t } t R + be the Brakke flow defined in Theorem 3.6, and let μ be the space-time measure d μ = d V t d t . Then

(4.4) ( h ( x 0 , V t 0 ) 1 ) T ( x 0 , t 0 ) μ

at μ-a.e. ( x 0 , t 0 ) such that the tangent space T ( x 0 , t 0 ) μ exists.

4.2 Existence of measure-theoretic velocities

We have the following proposition.

Proposition 4.5.

Let { E ( t ) } t R + and { V t } t R + be as in Theorem 3.6, and let μ Ω denote the Radon measure μ Ω := Ω d V t d t . Also, for every i { 1 , , N } , set

S ( i ) := { ( x , t ) n + 1 × + : x E i ( t ) } .

Then, for every i, there exists a function u i = u i ( x , t ) L 2 ( μ Ω ) (on R n + 1 × [ 0 , T ] for arbitrary T > 0 ) with spt ( u i ) S ( i ) and such that

(4.5) n + 1 ϕ ( x , t ) χ i ( x , t ) d x | t = t 1 t 2 = t 1 t 2 n + 1 ϕ t ( x , t ) χ i ( x , t ) d x d t + t 1 t 2 n + 1 ϕ ( x , t ) u i ( x , t ) d V t d t

for every 0 t 1 < t 2 < and all test functions ϕ C c 1 ( R n + 1 × R + ) .

Proposition 4.5 readily implies the following corollary.

Corollary 4.6.

For every i { 1 , , N } ,

χ i BV loc ( n + 1 × + ) 𝑎𝑛𝑑 χ i ( x , t ) = ( χ E i ( t ) d t , u i ( x , t ) d V t d t ) .

In particular, since χ i is the indicator function of S ( i ) , the latter is a set of locally finite perimeter in R n + 1 × R + .

Proof of Proposition 4.5.

We shall divide the proof into several steps. At first, we will drop the dependence of the test function on time, which will be introduced again at a later step. The idea of the proof is to obtain an evolution identity for the approximating flow j ( t ) defined in Section 3.2 (where j is the sequence of Theorem 3.6), and then to prove that, as , such identity converges to (4.5). Notice that we can assume without loss of generality that t 1 = 0 , since then the seemingly more general identity can be retrieved by simply taking differences. We will then also rename t 2 = T ( 0 , ) .

Step one: preliminary reductions. By density, it is evidently sufficient to show that identity (4.5) holds when ϕ = ϕ ( x ) is a function in C c 2 ( n + 1 ) . We may assume without loss of generality that | ϕ | Ω . Consider the approximating flow j ( t ) , and fix so large that the flow is defined on the interval [ 0 , T ] . We will deduce the validity of an approximate identity for the approximating flow, with vanishing errors in the limit as . We fix the index i { 1 , , N } , and, for the sake of simplicity in the notation, we drop the corresponding subscripts, so to write ( t ) in place of j ( t ) and E ( t ) in place of E j , i ( t ) . Recalling the construction of ( t ) , we let k T be the integer in { 1 , , j 2 p j } such that T ( ( k T - 1 ) Δ t , k T Δ t ] (with Δ t = Δ t j ), so that E ( T ) = E ( k T Δ t ) = E k T := E j , k T , i . We then have the discretization

(4.6) E ( T ) ϕ ( x ) 𝑑 x - E 0 ϕ ( x ) 𝑑 x = k = 0 k T - 1 { E k + 1 ϕ ( x ) 𝑑 x - E k ϕ ( x ) 𝑑 x } ,

and we can further decompose each summand on the right-hand side of (4.6) as

(4.7) E k + 1 ϕ ( x ) 𝑑 x - E k ϕ ( x ) 𝑑 x = D 1 + D 2 ,

where

D 1 := E k + 1 * ϕ ( x ) 𝑑 x - E k ϕ ( x ) 𝑑 x ,
D 2 := E k + 1 ϕ ( x ) 𝑑 x - E k + 1 * ϕ ( x ) 𝑑 x .

Here, E k + 1 * = int ( f 1 ( E k ) ) for some Lipschitz function f 1 𝐄 vc ( k , j ) satisfying the almost-optimality condition

( f 1 ) k ( Ω ) - k ( Ω ) ( 1 - j - 5 ) Δ vc k ( Ω ) ,

where Δ vc k ( Ω ) := Δ j vc k ( Ω ) , and E k + 1 = f 2 ( E k + 1 * ) , with f 2 ( x ) = x + Δ t h ε ( x , k + 1 * ) ( ε = ε j ). We will now proceed evaluating D 1 and D 2 separately.

Step two: evaluation of D 1 and D 2 . For D 1 , we simply observe that

D 1 = E k + 1 * E k ϕ ( x ) 𝑑 x - E k E k + 1 * ϕ ( x ) 𝑑 x ,

so that, using Definition 3.1 (ii) and (3.1), it holds

(4.8) | D 1 | 2 ϕ C 0 n + 1 ( E k E k + 1 * ) - 2 ϕ C 0 Δ vc k ( Ω ) Δ t Δ t j .

We then proceed with D 2 , and in order to further ease the notation we set f ( x ) := f 2 ( x ) , F ( x ) = f ( x ) - x , as well as g := f - 1 and G ( x ) := g ( x ) - x . We also simply write

E := E k + 1 , E * := E k + 1 * = g ( E ) .

Using that g is a diffeomorphism and changing variable in the second integral in D 2 , we have

D 2 = E ϕ ( x ) 𝑑 x - E ϕ ( g ( x ) ) 𝐉 g ( x ) 𝑑 x
(4.9) = E { ϕ ( x ) - ϕ ( g ( x ) ) } 𝑑 x + E ϕ ( x ) { 1 - 𝐉 g ( x ) } 𝑑 x + E ( ϕ ( x ) - ϕ ( g ( x ) ) ) ( 𝐉 g ( x ) - 1 ) 𝑑 x ,

where 𝐉 g denotes the Jacobian determinant of g. Using that G ( x ) = - F ( g ( x ) ) and that

F ( y ) = Δ t h ε ( y ) , h ε ( ) = h ε ( , k + 1 * ) ,

together with (2.15), we have that if ( spt ϕ ) 1 := { dist ( x , spt ϕ ) 1 } , then

(4.10) | ϕ ( g ( x ) ) - ϕ ( x ) | ϕ C 0 | G ( x ) | χ ( spt ϕ ) 1 ( x ) Δ t ε - 3 χ ( spt ϕ ) 1 ( x ) ,
(4.11) | ϕ ( g ( x ) ) - ϕ ( x ) - ϕ ( x ) G ( x ) | 2 ϕ C 0 | G ( x ) | 2 χ ( spt ϕ ) 1 ( x ) Δ t ε c 2 - 6 χ ( spt ϕ ) 1 ( x ) .

Notice that, in order to apply estimate (2.15), we are using that k + 1 * ( Ω ) k ( Ω ) and we are assuming that ε = ε j is smaller than the ε 1 given by Lemma 2.18 corresponding to the choice

M = sup t [ 0 , T ] j ( t ) ( Ω ) ,

which is bounded by a quantity depending only on T, the initial datum and c 1 due to (3.6). In deducing (4.10) and (4.11), we have also used that, by the definition of g and (2.15), and for all large enough, if x ( spt ϕ ) 1 , then g ( x ) spt ϕ , and we have used that ε can be assumed sufficiently small, so that ε ϕ C 0 and ε 2 ϕ C 0 are both bounded by 1. Next, we also estimate

(4.12) | 𝐉 g ( x ) - 1 | c ( n ) G ( x ) c ( n ) Δ t ε - 4 ,
(4.13) | 𝐉 g ( x ) - 1 - div ( G ( x ) ) | c ( n ) G ( x ) 2 Δ t ε c 2 - 9 ,

where we have used that

G = - [ ( F ) g ] g = - [ ( F ) g ] [ ( f ) - 1 g ] ,

that F = Δ t h ε , and estimate (2.16). We can then conclude from (4.11) and (4.13) that the sum of the first two integrals in (4.9) is

= - E { ϕ ( x ) G ( x ) + ϕ ( x ) div ( G ( x ) ) } 𝑑 x + Err 1 ,

where

(4.14) | Err 1 | Δ t ε c 2 - 10

as soon as ε is small enough to absorb the constants depending on n + 1 ( ( spt ϕ ) 1 ) . On the other hand, using that

ϕ ( x ) G ( x ) + ϕ ( x ) div ( G ( x ) ) = div ( ϕ ( x ) G ( x ) )

and the divergence theorem, we have

- E { ϕ ( x ) G ( x ) + ϕ ( x ) div ( G ( x ) ) } 𝑑 x = - * E ϕ ( x ) G ( x ) ν E ( x ) 𝑑 n ( x )
(4.15) = Δ t * E ϕ ( x ) h ε ( g ( x ) ) ν E ( x ) 𝑑 n ( x ) .

Recall that, in the right-hand side of (4.15), h ε ( ) = h ε ( , k + 1 * ) . We then continue the chain of identities in (4.15) as

Δ t * E ϕ ( x ) h ε ( x , k + 1 ) ν E ( x ) 𝑑 n ( x ) + Err 2 ,

where

Err 2 := Δ t * E ϕ ( x ) ν E ( x ) ( h ε ( g ( x ) , k + 1 * ) - h ε ( x , k + 1 ) ) 𝑑 n ( x ) .

Since k + 1 = f k + 1 * , and f is a diffeomorphism, we have that k + 1 = f k + 1 * , where f denotes the push-forward operator on integral varifolds through f. Setting, for simplicity,

V = k + 1 * and V ^ = k + 1 = f V ,

we can write

(4.16) h ε ( x , V ^ ) - h ε ( g ( x ) , V ) = [ h ε ( x , V ^ ) - h ε ( g ( x ) , V ^ ) ] + [ h ε ( g ( x ) , V ^ ) - h ε ( g ( x ) , V ) ] .

Recall now that the ε-smoothed mean curvature vector h ε ( , V ) is defined by

h ε ( , V ) = - Φ ε Φ ε δ V Φ ε V + ε Ω - 1 .

In particular, using [21, (5.74) and Lemma 4.17], it is easy to estimate the second summand on the right-hand side of (4.16) by

| h ε ( g ( x ) , V ^ ) - h ε ( g ( x ) , V ) | ε c 2 - 2 n - 14 ( V ( Ω ) + V ( Ω ) 2 ) .

For the first summand on the right-hand side of (4.16), instead,

| h ε ( x , V ^ ) - h ε ( g ( x ) , V ^ ) | | g ( x ) - x | h ε C 0 Δ t h ε C 0 h ε C 0 4 ε c 2 - 6

by Lemma 2.18. Thus, we can finally estimate

| Err 2 | Δ t ε c 2 - 2 n - 14 V ^ ( Ω ) ( V ( Ω ) + V ( Ω ) 2 + ε 2 n + 7 )
(4.17) Δ t ε c 2 - 2 n - 14 V ^ ( Ω ) ( V ^ ( Ω ) + V ^ ( Ω ) 2 + ε 2 n + 7 ) ,

where V ^ = k and we used V ( Ω ) = ( f 1 ) k ( Ω ) k ( Ω ) . Concerning the third integral in (4.9), instead, we use (4.10) and (4.12) to estimate

Err 3 := E ( ϕ ( x ) - ϕ ( g ( x ) ) ) ( 𝐉 g ( x ) - 1 ) 𝑑 x

by

(4.18) | Err 3 | Δ t ε c 2 - 8

for all ε sufficiently small. We have thus proved that

(4.19) E k + 1 ϕ ( x ) 𝑑 x - E k ϕ ( x ) 𝑑 x = Δ t ( * E k + 1 ϕ ( x ) ν E k + 1 h ε ( x , k + 1 ) 𝑑 n ( x ) + Err k + 1 ) ,

where

Err k + 1 := ( Δ t ) - 1 ( D 1 + Err 1 + Err 2 + Err 3 )

satisfy estimates (4.8), (4.14), (4.17) and (4.18).

Step three: sum of contributions and limit . Now, we sum the identities (4.19) for k = 0 , , k T - 1 . Introducing back the subscripts j and i, and using that (see (3.5))

j ( t ) = j , k + 1 for all  t ( k Δ t j , ( k + 1 ) Δ t j ] ,

we obtain from (4.6) that

(4.20) E j , i ( T ) ϕ ( x ) 𝑑 x - E 0 , i ϕ ( x ) 𝑑 x = 0 k T Δ t j n + 1 ϕ ( x ) h ε j ( x , j ( t ) ) 𝑑 μ E j , i ( t ) ( x ) 𝑑 t + Err j ,

with

μ E j , i ( t ) = - χ E j , i ( t ) = ν E j , i ( t ) n * E j , i ( t ) .

Moreover, since k T Δ t j T + 1 , we obtain

| Err j | 1 j 0 T + 1 - Δ j vc j ( t ) ( Ω ) Δ t j d t - Δ j vc 0 ( Ω ) j + ε j c 2 - 11 ( T + 1 )
(4.21) + ε j c 2 - 2 n - 14 ( T + 1 ) ( M j 2 + M j 3 + ε j 2 n + 7 ) ,

where M j := sup t [ 0 , T + 1 ] j ( t ) ( Ω ) . Next, Lemma 2.18 implies that

| T k T Δ t j n + 1 ϕ h ε j ( , j ( t ) ) 𝑑 μ E j , i ( t ) | M j ε j c 2 - 3 ,

so that the end-point k T Δ t j in the time integral of (4.20) can be replaced by T. Letting in the identity (4.20), we have that the left-hand side converges to

E i ( T ) ϕ ( x ) 𝑑 x - E i , 0 ϕ ( x ) 𝑑 x

due to (3.10), whereas Err j 0 thanks to (4.21), (3.6) and (3.7). We are then left with studying the limit

(4.22) lim 0 T n + 1 ϕ ( x ) h ε j ( x , j ( t ) ) 𝑑 μ E j , i ( t ) ( x ) 𝑑 t .

To this end, we first observe that, since

h ε ( x , V ) = - Φ ε ( Φ ε δ V Φ ε V + ε Ω - 1 ) ( x ) ,

and using that, by definition, Ω ( x ) Ω ( y ) exp ( c 1 | x - y | ) for all x , y n + 1 together with the properties of the kernel Φ ε , we can then estimate

* E j , i ( t ) Ω | h ε j ( , j ( t ) ) | 2 𝑑 n C n + 1 | Φ ε j δ ( j ( t ) ) | 2 Ω ( y ) Φ ε j j ( t ) + ε j Ω - 1 Φ ε j ( n * E j , i ( t ) ) Φ ε j j ( t ) + ε j Ω - 1 𝑑 y
(4.23) C n + 1 | Φ ε j δ ( j ( t ) ) | 2 Ω ( y ) Φ ε j j ( t ) + ε j Ω - 1 𝑑 y .

Therefore, for every T > 0 and for every ϕ C c ( n + 1 × [ 0 , T ] ) it holds

| 0 T n + 1 ϕ ( x , t ) Ω ( x ) | h ε j ( x , j ( t ) ) d μ E j , i ( t ) ( x ) | 𝑑 t |
( 0 T * E j , i ( t ) Ω | h ε j ( x , j ( t ) ) | 2 𝑑 n 𝑑 t ) 1 2 ( 0 T * E j , i ( t ) Ω | ϕ ( x , t ) | 2 𝑑 n 𝑑 t ) 1 2
(4.24) C ( 0 T n + 1 | Φ ε j δ ( j ( t ) ) | 2 Ω ( y ) Φ ε j j ( t ) + ε j Ω - 1 𝑑 y ) 1 2 ( 0 T * E j , i ( t ) Ω | ϕ ( x , t ) | 2 𝑑 n 𝑑 t ) 1 2 ,

where the last line follows from (4.23). In particular, using (3.6) and (3.7), for any fixed T > 0 ,

Ω h ε j ( , j ( t ) ) d μ E j , i ( t ) d t

are (signed) Radon measures in n + 1 × [ 0 , T ] with uniformly bounded total variation, and we can let σ i denote a subsequential limit as . Since n * E j , i ( t ) j ( t ) , and the latter measures converge, as , to V t for every t + , it is clear from (4.24) that σ i is absolutely continuous with respect to the measure μ Ω = Ω d V t d t . We will let u i denote the Radon–Nikodým derivative of σ i with respect to μ Ω , so that

(4.25) σ i = u i Ω d V t d t

in the sense of measures. Estimate (4.24) also implies that u i L 2 ( μ Ω ) , with

u i L 2 ( μ Ω ) C lim sup ( 0 T n + 1 | Φ ε j δ ( j ( t ) ) | 2 Ω ( y ) Φ ε j j ( t ) + ε j Ω - 1 𝑑 y ) 1 2 .

Hence, we can now calculate the limit in (4.22) (along the aforementioned subsequence, not relabeled), which gives

lim 0 T n + 1 ( ϕ Ω - 1 ) Ω h ε j ( , j ( t ) ) 𝑑 μ E j , i ( t ) 𝑑 t = 0 T n + 1 ϕ u i d V t 𝑑 t .

We have then concluded the existence of u i L 2 ( μ Ω ) such that

E i ( T ) ϕ ( x ) 𝑑 x - E i , 0 ϕ ( x ) 𝑑 x = 0 T n + 1 ϕ ( x ) u i ( x , t ) d V t ( x ) 𝑑 t

for all ϕ C c ( n + 1 ) , that is, we have obtained (4.5) when ϕ does not depend on t.

Step four: the case of time-dependent ϕ. We now consider the general case of a test function

ϕ C c 1 ( n + 1 × [ 0 , T ] ) .

The proof is analogous, with a few minor modifications which take the dependence on t into account. First, formula (4.6) becomes

E ( T ) ϕ ( x , k T Δ t ) 𝑑 x - E 0 ϕ ( x , 0 ) 𝑑 x = k = 0 k T - 1 { E k + 1 ϕ ( x , ( k + 1 ) Δ t ) 𝑑 x - E k ϕ ( x , k Δ t ) 𝑑 x } ,

whereas formula (4.7) becomes

E k + 1 ϕ ( x , ( k + 1 ) Δ t ) 𝑑 x - E k ϕ ( x , k Δ t ) 𝑑 x = D 1 + D 2 + D 3 ,

where

D 1 := E k + 1 * ϕ ( x , k Δ t ) 𝑑 x - E k ϕ ( x , k Δ t ) 𝑑 x ,
D 2 := E k + 1 ϕ ( x , k Δ t ) 𝑑 x - E k + 1 * ϕ ( x , k Δ t ) 𝑑 x ,
D 3 := E k + 1 ϕ ( x , ( k + 1 ) Δ t ) 𝑑 x - E k + 1 ϕ ( x , k Δ t ) 𝑑 x .

By the fundamental theorem of calculus and Fubini’s theorem, we can immediately calculate

D 3 = k Δ t ( k + 1 ) Δ t E k + 1 ϕ t ( x , t ) 𝑑 x 𝑑 t .

The summands D 1 and D 2 are, instead, treated as in the time-independent case, with ϕ replaced by ϕ ( , k Δ t ) . Identity (4.19) then becomes

E k + 1 ϕ ( x , ( k + 1 ) Δ t ) 𝑑 x - E k ϕ ( x , k Δ t ) 𝑑 x
= k Δ t ( k + 1 ) Δ t E k + 1 ϕ t ( x , t ) 𝑑 x 𝑑 t + Δ t ( * E k + 1 ϕ ( x , k Δ t ) ν E k + 1 h ε ( x , k + 1 ) 𝑑 n ( x ) + Err k + 1 ) .

Introducing back the subscripts j and i , we can then write

Δ t * E k + 1 ϕ ( x , k Δ t ) ν E k + 1 h ε ( x , k + 1 ) 𝑑 n ( x )
= k Δ t j ( k + 1 ) Δ t j * E j , i ( t ) ϕ ( x , t ) ν E j , i ( t ) h ε j ( x , j ( t ) ) 𝑑 n ( x ) 𝑑 t + Err ~ k + 1 ,

where

| Err ~ k + 1 | Δ t j t ϕ C 0 ( max spt ϕ Ω - 1 ) k Δ t j ( k + 1 ) Δ t j * E j , i ( t ) Ω | h ε j ( x , j ( t ) ) | 𝑑 n ( x ) 𝑑 t .

Summing over k = 0 , , k T - 1 , we can then replace (4.20) with

E j , i ( T ) ϕ ( x , k T Δ t j ) 𝑑 x - E 0 , i ϕ ( x , 0 ) 𝑑 x
(4.26) = 0 k T Δ t j E j , i ( t ) ϕ t ( x , t ) 𝑑 x 𝑑 t + 0 k T Δ t j n + 1 ϕ ( x , t ) h ε j ( x , j ( t ) ) 𝑑 μ E j , i ( t ) ( x ) 𝑑 t + Err ~ j ,

where Err ~ j contains, with respect to Err j , an additional error term which can be estimated by

| Err ~ j - Err j | C ε j c 2 - 2 t ϕ C 0 ( max spt ϕ Ω - 1 ) M j ( T + 1 ) ,

where we have used (2.15). In particular, when , also Err ~ j 0 . Next, also in this case we can replace k T Δ t j with T in the integrals on both sides of (4.26), paying an additional error term which we estimate by

| T k T Δ t j E j , i ( t ) ϕ t 𝑑 x 𝑑 t | C ε j c 2 t ϕ C 0 n + 1 ( ( spt ϕ ) 1 ) ,

and similarly for others, and they all vanish as . We can then finally conclude that

(4.27) E i ( T ) ϕ ( x , T ) 𝑑 x - E i , 0 ϕ ( x , 0 ) 𝑑 x = 0 T E i ( t ) ϕ t ( x , t ) 𝑑 x 𝑑 t + 0 T n + 1 ϕ ( x , t ) u i ( x , t ) d V t ( x ) 𝑑 t ,

where u i L 2 ( μ Ω ) is defined by (4.25).

Step five: support of u i . To conclude the proof of Proposition 4.5, it only remains to show that spt ( u i ) S ( i ) , where

S ( i ) = { ( x , t ) n + 1 × + : x E i ( t ) } .

Notice that S ( i ) is open (relatively to n + 1 × + ) by [21, Theorem 3.5 (5)]. Suppose now that ( x ^ , t ^ ) S ( i ) . Then either ( x ^ , t ^ ) S ( i ) or ( x ^ , t ^ ) clos ( S ( i ) ) . In the first case, there is a neighborhood ( x ^ , t ^ ) U n + 1 × + such that x E i ( t ) for all ( x , t ) U . In particular, ( x ^ , t ^ ) spt μ ( spt V 0 × { 0 } ) , and since σ i is absolutely continuous with respect to μ, also ( x ^ , t ^ ) spt ( σ i ) = spt ( u i ) .

In the second case, we have that, necessarily, x ^ clos ( E i ( t ^ ) ) . Of course, if x ^ E i ( t ^ ) for some i i , then ( x ^ , t ^ ) S ( i ) , and the above proof implies again that

( x ^ , t ^ ) spt μ ( spt V 0 × { 0 } ) ,

which then gives ( x ^ , t ^ ) spt ( u i ) . Thus, we can assume that x ^ Γ ( t ^ ) clos ( E i ( t ^ ) ) . Since ( x ^ , t ^ ) clos ( S ( i ) ) , there is an open neighborhood ( x ^ , t ^ ) U n + 1 × + such that U clos ( S ( i ) ) = , and we can choose U of the form U = U r ( x ^ ) × [ 0 , r 2 ) if t ^ = 0 , or U = U r ( x ^ ) × ( t ^ - r 2 , t ^ + r 2 ) if t ^ > 0 . In both cases, if ( x , t ) U , then x E i ( t ) . For any function ϕ C c 1 ( U ) , identity (4.27) then gives

0 = U ϕ ( x , t ) u i ( x , t ) d V t ( x ) 𝑑 t = U ϕ ( x , t ) Ω - 1 ( x ) 𝑑 σ i ( x , t ) .

Since ϕ is arbitrary, we deduce then that | σ i | ( U ) = 0 , and thus that ( x ^ , t ^ ) spt ( σ i ) = spt ( u i ) . The proof of the theorem is now complete. ∎

4.3 Boundaries move by their mean curvature

In this subsection, we complete the proof of Theorem 4.1, by achieving the representation for the measure-theoretic time derivative t χ i as specified in Remark 4.2. The crucial step is to prove the following lemma, for which the L 2 flow property of { V t } plays a pivotal role. We shall denote by 𝐩 and 𝐪 the projections of n + 1 × onto its factors, so that 𝐩 ( x , t ) = x and 𝐪 ( x , t ) = t .

Lemma 4.7.

For every i { 1 , , N } , there is a set G i * S ( i ) with H n + 1 ( * S ( i ) G i ) = 0 such that the following assertions hold for every ( x , t ) G i :

  1. The tangent T ( x , t ) μ exists, and T ( x , t ) μ = T ( x , t ) ( * S ( i ) ) .

  2. Relation ( 4.4 ) holds.

  3. x * E i ( t ) and T x V t = T x ( * E i ( t ) ) .

  4. 𝐩 ( ν S ( i ) ( x , t ) ) 0 and ν E i ( t ) ( x ) = | 𝐩 ( ν S ( i ) ( x , t ) ) | - 1 𝐩 ( ν S ( i ) ( x , t ) ) .

  5. T x ( * E i ( t ) ) × { 0 } is a linear subspace of T ( x , t ) μ .

The proof of Lemma 4.7 can be deduced with relatively little effort from arguments already contained in [29], but we include it for the reader’s convenience. Before proceeding, we will need some consequences of the celebrated monotonicity formula by Huisken for MCF. Since said consequences will be needed also later on in the paper, we record them here. First, let us set some notation. For ( y , s ) n + 1 × + , we define the n-dimensional backwards heat kernel ρ ( y , s ) by

(4.28) ρ ( y , s ) ( x , t ) := 1 ( 4 π ( s - t ) ) n 2 exp ( - | x - y | 2 4 ( s - t ) ) for all  t < s  and  x n + 1 ,

as well as the truncated kernel

(4.29) ρ ^ ( y , s ) r ( x , t ) := η ( x - y r ) ρ ( y , s ) ( x , t ) ,

where r > 0 and η C c ( U 2 ; + ) is a radially symmetric function such that η 1 on B 1 , 0 η 1 , | η | 2 and 2 η 4 . The following lemma is [21, Lemma 10.3], and it is a variant of Huisken’s monotonicity formula for MCF.

Lemma 4.8.

There exists c ( n ) > 0 with the following property. For every 0 t 1 < t 2 < s < , y R n + 1 and r > 0 , it holds

V t ( ρ ^ ( y , s ) r ( , t ) ) | t = t 1 t 2 c ( n ) r - 2 ( t 2 - t 1 ) sup t [ t 1 , t 2 ] r - n V t ( B 2 r ( y ) ) .

As a consequence, one has the following lemma [21, Lemma 10.4], which provides a uniform upper bound on mass density ratios.

Lemma 4.9.

For every L > 1 , there exists Λ = Λ ( n , L , Ω , E 0 ( Ω ) ) ( 1 , ) such that

(4.30) sup { r - n V t ( B r ( x ) ) : r ( 0 , 1 ] , x B L ( 0 ) , t [ L - 1 , L ] } Λ .

Proof of Lemma 4.7.

Fix i { 1 , , N } , and observe that, since χ i BV loc ( n + 1 × + ) is the characteristic function of S ( i ) , it holds | χ i | = n + 1 * S ( i ) . We then let g i be the Radon–Nikodým derivative of μ with respect to | χ i | , namely

(4.31) g i ( x , t ) = d μ d | χ i | ( x , t ) := lim r 0 + μ ( B r n + 2 ( x , t ) ) | χ i | ( B r n + 2 ( x , t ) ) ,

where B r n + 2 ( x , t ) is the closed ball with radius r and center ( x , t ) in n + 1 × . By the Lebesgue–Radon–Nikodým theorem, it holds g i ( x , t ) < for n + 1 -a.e. ( x , t ) * S ( i ) , and

(4.32) μ = g i n + 1 * S ( i ) + μ Σ i for a set  Σ i  with  n + 1 ( * S ( i ) Σ i ) = 0 .

On the other hand, it is not difficult to see that

(4.33) μ n + 1 .

Indeed, first notice that, as an immediate corollary of Lemma 4.9, one has

(4.34) Θ * n + 1 ( μ , ( x , t ) ) := lim sup r 0 + μ ( B r n + 2 ( x , t ) ) ω n + 1 r n + 1 < for every  ( x , t ) n + 1 × ( 0 , ) .

Let then A be a set with n + 1 ( A ) = 0 , and set, for Q ,

D Q := { ( x , t ) n + 1 × ( 0 , ) : Θ * n + 1 ( μ , ( x , t ) ) Q } .

Then, by [33, Theorem 3.2], one has

μ ( D Q A ) 2 n + 1 Q n + 1 ( D Q A ) = 0 ,

so that the conclusion μ ( A ) = 0 follows because n + 1 × ( 0 , ) = Q D Q due to (4.34).

Combining (4.32) with (4.33), we immediately have that

(4.35) μ * S ( i ) = g i n + 1 * S ( i ) .

Similarly, by taking the Radon–Nikodým derivative of | χ i | with respect to μ, we see that

d | χ i | d μ ( x , t ) is finite for  μ -a.e.  ( x , t ) .

Since | χ i | μ by Theorem 2.11 (iv) and Corollary 4.6, this is also finite for | χ i | -a.e. ( x , t ) and (4.31) shows that this is equal to 1 / g ( x , t ) , so that, in particular, g ( x , t ) > 0 for n + 1 -a.e. ( x , t ) * S ( i ) . Together with (4.35), the latter information implies that μ * S ( i ) is an ( n + 1 ) -rectifiable measure, so that T ( x , t ) ( μ * S ( i ) ) exists and is equal to T ( x , t ) ( * S ( i ) ) for n + 1 -a.e. ( x , t ) * S ( i ) . On the other hand, by [33, Theorem 3.5], for n + 1 -a.e. ( x , t ) * S ( i ) it holds

(4.36) lim r 0 + μ ( B r n + 2 ( x , t ) * S ( i ) ) r n + 1 = 0 .

If ( x 0 , t 0 ) * S ( i ) is such that (4.36) holds and ϕ C c ( B 1 n + 2 ) is arbitrary, then

lim sup r 0 + | ( n + 1 × + ) * S ( i ) ϕ ( r - 1 ( x - x 0 , t - t 0 ) ) r - ( n + 1 ) 𝑑 μ ( x , t ) |
ϕ C 0 lim r 0 + μ ( B r n + 2 ( x 0 , t 0 ) * S ( i ) ) r n + 1 = 0 .

This shows that

lim r 0 + ϕ ( r - 1 ( x - x 0 , t - t 0 ) ) r - ( n + 1 ) 𝑑 μ ( x , t ) = lim r 0 + * S ( i ) ϕ ( r - 1 ( x - x 0 , t - t 0 ) ) r - ( n + 1 ) 𝑑 μ ( x , t )
= g i ( x 0 , t 0 ) T ( x 0 , t 0 ) ( * S ( i ) ) ϕ ( x , t ) 𝑑 n + 1 ( x , t )

for every ϕ C c ( n + 1 × ) and for n + 1 -a.e. ( x 0 , t 0 ) * S ( i ) . This completes the proof of (i), whereas (ii) is an immediate consequence of (i) and Corollary 4.4, since, on * S ( i ) , μ and n + 1 have the same null sets by (4.35).

Next, to prove (iii) and (iv), we are going to use the coarea formula for sets of finite perimeter, slicing * S ( i ) according to hyperplanes { 𝐪 = t } . Set, for t + ,

( * S ( i ) ) t := * S ( i ) { 𝐪 = t } = { x : ( x , t ) * S ( i ) } .

Then using, e.g., [26, Theorem 18.11] together with the fact that S ( i ) { 𝐪 = t } = E i ( t ) , we have that

(4.37) n ( ( * S ( i ) ) t Δ * E i ( t ) ) = 0

for a.e. t + . Moreover, for such t,

(4.38) 𝐩 ( ν S ( i ) ( x , t ) ) 0 ,
(4.39) ν E i ( t ) ( x ) = 𝐩 ( ν S ( i ) ( x , t ) ) | 𝐩 ( ν S ( i ) ( x , t ) ) |

for n -a.e. x ( * S ( i ) ) t . Let

Z := { t + : ( 4.37 )  fails } ,

and for every t + set

Z t := { x ( * S ( i ) ) t : x * E i ( t )  or (4.38) and (4.39) fail } ,

so that 1 ( Z ) = 0 and, for every t Z , n ( Z t ) = 0 . Consider then the Borel function κ ( x , t ) := χ Z t ( x ) on n + 1 × + . We have

0 = 0 χ + Z ( t ) n ( Z t ) 𝑑 t
= 0 n ( Z t ) 𝑑 t
= 0 * S ( i ) { 𝐪 = t } κ ( x , t ) 𝑑 n ( x ) 𝑑 t
= * S ( i ) κ ( x , t ) | * S ( i ) 𝐪 ( x , t ) | 𝑑 n + 1 ( x , t ) ,

where in the last identity we have used the coarea formula [3, (2.72)] and where * S ( i ) 𝐪 ( x , t ) is the tangential gradient, along * S ( i ) , of 𝐪 at ( x , t ) . Combining now (i) and Corollary 4.4, we see that

( h ( x , V t ) 1 ) T ( x , t ) ( * S ( i ) ) at  n + 1 -a.e.  ( x , t ) * S ( i ) ,

which readily implies that | * S ( i ) 𝐪 ( x , t ) | > 0 for n + 1 -a.e. ( x , t ) * S ( i ) . Hence, it must be κ ( x , t ) = 0 for n + 1 -a.e. ( x , t ) * S ( i ) , thus proving the first parts of (iii) and (iv). At such points, the identity

T x V t = T x ( * E i ( t ) )

is then obtained by repeating the argument in (i) at fixed t.

Finally, (v) is always true at points satisfying (i)–(iv): indeed, if ξ T x ( * E i ( t ) ) , then

( ξ ) ν S ( i ) ( x , t ) = ξ 𝐩 ( ν S ( i ) ( x , t ) ) = | 𝐩 ( ν S ( i ) ( x , t ) ) | ξ ν E i ( t ) ( x ) = 0 .

This completes the proof. ∎

Proof of Theorem 4.1.

By virtue of Proposition 4.5, we only need to show the validity of equation (4.2). Fix i { 1 , , N } . Using that S ( i ) is a set of locally finite perimeter in n + 1 × + , for any 0 t 1 < t 2 < and any ϕ C c 1 ( n + 1 × ( t 1 , t 2 ) ) we have that

n + 1 × + ϕ t ( x , t ) χ i ( x , t ) 𝑑 x 𝑑 t = S ( i ) ϕ t ( x , t ) 𝑑 x 𝑑 t
(4.40) = * S ( i ) ϕ ( x , t ) 𝐪 ( ν S ( i ) ( x , t ) ) 𝑑 n + 1 ( x , t ) .

Let G i be the set of Lemma 4.7. For every ( x , t ) G i , we have

T ( x , t ) μ = ( T x ( * E i ( t ) ) × { 0 } ) span ( h ( x , V t ) 1 ) ,

and thus

(4.41) ν S ( i ) ( x , t ) = 1 1 + | h ( x , V t ) | 2 ( ν E i ( t ) ( x ) - h ( x , V t ) ν E i ( t ) ( x ) ) ,

where we have used that h ( x , V t ) T x V t = T x ( * E i ( t ) ) . In particular, since n + 1 ( * S ( i ) G i ) = 0 , equation (4.40) yields

(4.42) n + 1 × + ϕ t ( x , t ) χ i ( x , t ) 𝑑 x 𝑑 t = - G i ϕ ( x , t ) h ( x , V t ) ν E i ( t ) ( x ) 1 1 + | h ( x , V t ) | 2 𝑑 n + 1 ( x , t ) .

Due to (4.41), at points ( x , t ) G i the coarea factor with respect to the slicing via 𝐪 is precisely

| * S ( i ) 𝐪 ( x , t ) | = 1 1 + | h ( x , V t ) | 2 ,

so that the chain of identities in (4.42) can be continued as

- 0 G i { 𝐪 = t } ϕ ( x , t ) h ( x , V t ) ν E i ( t ) ( x ) 𝑑 n ( x ) 𝑑 t = - 0 * E i ( t ) ϕ ( x , t ) h ( x , V t ) ν E i ( t ) ( x ) 𝑑 n ( x ) 𝑑 t ,

again by the coarea formula. This completes the proof. ∎

5 Ilmanen’s density lower bound

This section contains a proof of the following fact (Theorem 5.3 below), from which we may conclude the proofs of Theorem 2.10 (ii)–(iv): if { V t } t + is as in Proposition 3.6, then there exists a threshold θ 0 > 0 such that, at “most” points (in a sense to be suitably specified) on the support of the evolving varifolds the mass density is not smaller than θ 0 . Such density lower bound is already stated in the work of Ilmanen [18, Section 7.1], and we include here its proof only for the sake of completeness.

Using Lemma 4.8 and Lemma 4.9, we prove next the following lemma, which is a variant of Brakke’s clearing out lemma (cf. [6, Section 6.3], [18, Section 6.1] and [21, Lemma 10.6]). Recalling the notation set in (4.28) and (4.29), we define, for δ > 0 ,

(5.1) ρ ^ y r , δ ( x ) := ρ ^ ( y , t + δ r 2 ) r ( x , t ) = η ( x - y r ) 1 ( 4 π δ r 2 ) n 2 exp ( - | x - y | 2 4 δ r 2 ) .

Lemma 5.1.

For any L > 1 , there exists

δ 0 = δ 0 ( n , L , Ω , 0 ( Ω ) ) ( 0 , 1 )

with the following property: if ( y , t ) B L ( 0 ) × [ L - 1 , L ] and r ( 0 , 1 2 ) are such that

(5.2) V t ( ρ ^ y r , δ 0 ) < 1 2 ,

then ( y , t + δ 0 r 2 ) spt μ .

Proof.

For δ 0 ( 0 , 1 ) to be specified later, suppose towards a contradiction that ( y , t + δ 0 r 2 ) spt μ . Then, by [21, Lemma 10.1 (2)], there is a sequence ( y i , t i ) converging to ( y , t + δ 0 r 2 ) such that V t i 𝐈𝐕 n ( n + 1 ) and y i spt V t i . Since V t i is an integral varifold 𝐯𝐚𝐫 ( M t i , θ t i ) , in any arbitrarily small neighborhood of a point x spt V t i there is a point x ~ M t i where M t i has an approximate tangent plane Tan ( M t i , x ~ ) and the density Θ V t i ( x ~ ) = θ t i ( x ~ ) is an integer. Thus, we can assume without loss of generality that Tan ( M t i , y i ) exists and Θ V t i ( y i ) 1 . Assume also without loss of generality that t i > t for all i, and fix ε > 0 . Then apply Lemma 4.8 with t 1 = t , t 2 = t i , s = t i + ε and y = y i to get

V s ( ρ ^ ( y i , t i + ε ) r ( , s ) ) | s = t t i c ( n ) r - 2 ( t i - t ) sup s [ t , t i ] r - n V s ( B 2 r ( y i ) ) .

Consider first the term

V t i ( ρ ^ ( y i , t i + ε ) r ( , t i ) ) = M t i η ( x - y i r ) ρ ( y i , ε ) ( x , 0 ) θ t i ( x ) 𝑑 n ( x )
M t i B r ( y i ) ρ ( y i , ε ) ( x , 0 ) θ t i ( x ) 𝑑 n ( x ) .

By the change of variable z = ( x - y i ) / ε in the integral on the right-hand side, it is not difficult to see that, in the limit as ε 0 + , the latter converges to

θ t i ( y i ) Tan ( M t i , y i ) ρ ( 0 , 1 ) ( z , 0 ) 𝑑 n ( z ) = θ t i ( y i ) 1 .

We can then conclude

1 V t ( ρ ^ ( y i , t i ) r ( , t ) ) + c ( n ) r - 2 ( t i - t ) sup s [ t , t i ] r - n V s ( B 2 r ( y i ) ) .

We let then i : using that ( y i , t i ) ( y , t + δ 0 r 2 ) and recalling (5.1), from (4.30) we get

1 V t ( ρ ^ y r , δ 0 ) + c ( n ) δ 0 sup s [ t , t + δ 0 r 2 ] r - n V s ( B 2 r ( y ) )
(5.3) V t ( ρ ^ y r , δ 0 ) + c ( n ) δ 0 Λ ( n , L + 1 , Ω , 0 ( Ω ) ) .

Choosing δ 0 such that the second summand in the right-hand side of (5.3) is less than or equal to 1 2 leads to a contradiction with (5.2). ∎

Remark 5.2.

Observe that, since

V t ( ρ ^ y r , δ 0 ) ( 4 π δ 0 ) - n 2 r - n V t ( U 2 r ( y ) ) ,

Lemma 5.1 immediately implies the following: if ( y , t ) B L ( 0 ) × [ L - 1 , L ] is such that, for some r ( 0 , 1 ) ,

(5.4) r - n V t ( U r ( y ) ) < θ 0 := ( 4 π δ 0 ) n 2 2 n + 1 ,

then ( y , t + 4 - 1 δ 0 r 2 ) spt μ .

Theorem 5.3.

For L > 1 , let θ 0 = θ 0 ( n , L + 1 , Ω , E 0 ( Ω ) ) be the number defined in (5.4), and define, for t 0 , the sets

𝒵 0 := { ( y , t ) spt μ ( B L ( 0 ) × [ L - 1 , L ] ) : lim sup r 0 + r - n V t ( U r ( y ) ) < θ 0 } ,
𝒵 t 0 := 𝒵 0 ( n + 1 × { t } ) .

Then there exists G R + with L 1 ( G ) = 0 such that H n - 1 + α ( Z t 0 ) = 0 for every α > 0 and every t R + G .

Proof.

For every θ < θ 0 and every σ ( 0 , 1 ) , define

𝒵 θ , σ 0 := { ( y , t ) spt μ ( B L ( 0 ) × [ L - 1 , L ] ) : r - n V t ( U r ( y ) ) < θ  for every  r ( 0 , σ ) } ,

so that

(5.5) 𝒵 0 = θ < θ 0 , σ ( 0 , 1 ) 𝒵 θ , σ 0 .

Let ( y , t ) 𝒵 θ , σ 0 . For any t ( t , t + 4 - 2 δ 0 σ 2 ] , let r > 0 be such that r 2 = 4 δ 0 - 1 ( t - t ) (so that necessarily r ( 0 , σ / 2 ] ), and let y B γ r ( y ) (with γ ( 0 , 1 ) to be specified) be arbitrary. Then it holds

( y , t ) B L + 1 ( 0 ) × [ ( L + 1 ) - 1 , L + 1 ] .

Furthermore, since U r ( y ) U ( 1 + γ ) r ( y ) , and since ( 1 + γ ) r < 2 r < σ , it holds

r - n V t ( U r ( y ) ) < ( 1 + γ ) n θ .

Hence, for γ small enough (depending on the ratio θ 0 / θ ) it holds r - n V t ( U r ( y ) ) < θ 0 , so that Remark 5.2 implies that

( y , t + 4 - 1 δ 0 r 2 ) = ( y , t ) spt μ .

In particular, ( y , t ) 𝒵 θ , σ 0 . Analogously, if t [ t - 4 - 2 δ 0 σ 2 , t ) and y B γ r ( y ) for r 2 = 4 δ 0 - 1 ( t - t ) and γ small enough, we have that necessarily ( y , t ) 𝒵 θ , σ 0 . Otherwise, it would be

( ( 1 + γ ) r ) - n V t ( U ( 1 + γ ) r ( y ) ) < θ ,

and thus

r - n V t ( U r ( y ) ) < ( 1 + γ ) n θ < θ 0 ,

which would then imply ( y , t + 4 - 1 δ 0 r 2 ) = ( y , t ) spt μ , a contradiction.

We have then concluded the following dichotomy:

either  ( y , t ) 𝒵 θ , σ 0  or  ( y , t ) 𝒵 θ , σ 0  whenever  0 < | t - t | 4 - 2 δ 0 σ 2  and  | y - y | 2 4 γ 2 δ 0 - 1 | t - t | .

Therefore, if ( y , t ) 𝒵 θ , σ 0 , then the truncated double paraboloid

𝒫 ( y , t ) := { | y - y | 2 4 γ 2 δ 0 - 1 | t - t | 4 - 1 γ 2 σ 2 }

intersects 𝒵 θ , σ 0 only in ( y , t ) . Next, setting 2 τ := 4 - 1 γ 2 σ 2 , we consider sets

𝒵 θ , σ , y 0 , t 0 0 := 𝒵 θ , σ 0 ( B 1 ( y 0 ) × [ t 0 - τ , t 0 + τ ] ) , ( y 0 B L ( 0 ) , t 0 [ L - 1 , L ] ) ,

so that a countable union of such sets covers 𝒵 θ , σ 0 . Fix any such set, and call it 𝒵 for the sake of simplicity: it will suffice to show that, setting 𝒵 t := 𝒵 ( n + 1 × { t } ) , it holds n - 1 + α ( 𝒵 t ) = 0 for every α > 0 and a.e. t 0 . Notice that if ( y , t ) 𝒵 , then, by the definition of τ, the set 𝒵 ( { y } × ) is contained in 𝒫 ( y , t ) : in particular, for y B 1 ( y 0 ) the fiber { y } × intersects 𝒵 in at most one point. Let 𝐩 be the coordinate projection 𝐩 ( x , t ) = x , fix δ > 0 , and cover the set 𝐩 ( 𝒵 ) B 1 ( y 0 ) by countably many open balls U r i ( y i ) so that

(5.6) r i δ , i ω n + 1 r i n + 1 2 n + 1 ( B 1 ( y 0 ) ) .

For every center y i of the balls in the covering, let t i be the only point such that ( y i , t i ) 𝒵 , and notice that, as a consequence of the first part of the proof, if ( y , t ) 𝒵 with y U r i ( y i ) , then necessarily | t - t i | < 4 - 1 γ - 2 δ 0 r i 2 . In other words, for δ suitably small, the cylinders

U r i ( y i ) × ( t i - 4 - 1 γ - 2 δ 0 r i 2 , t i + 4 - 1 γ - 2 δ 0 r i 2 )

are a covering of 𝒵 . We can then estimate

t 0 - τ t 0 + τ δ n - 1 + α ( 𝒵 t ) 𝑑 t t 0 - τ t 0 + τ i : | t - t i | < 4 - 1 γ - 2 δ 0 r i 2 ω n - 1 + α r i n - 1 + α d t
i t i - 4 - 1 γ - 2 δ 0 r i 2 t i + 4 - 1 γ - 2 δ 0 r i 2 ω n - 1 + α r i n - 1 + α 𝑑 t
C ( n , α ) γ - 2 δ 0 δ α ,

where we have used (5.6). Letting δ 0 + , we find then

t 0 - τ t 0 + τ n - 1 + α ( 𝒵 t ) 𝑑 t = 0

by monotone convergence. Hence, by taking countable unions,

0 n - 1 + α ( ( 𝒵 θ , σ 0 ) t ) 𝑑 t = 0 ,

with the obvious meaning of the symbols and taking into account that 𝒵 t is empty when t < L - 1 or t > L . The conclusion follows from (5.5). ∎

The following is the immediate corollary of Theorem 5.3, which proves Theorem 2.10 (ii).

Corollary 5.4.

There exists G R + with L 1 ( G ) = 0 such that for every t R + G ,

(5.7) n - 1 + α ( { y n + 1 : ( y , t ) spt μ and  lim r 0 + r - n V t ( U r ( y ) ) = 0 } ) = 0 for all  α > 0 .

In particular, by recalling from Theorem 2.10(i) and Theorem 2.11(iii) that spt V t { x : ( x , t ) spt μ } = Γ ( t ) , the set

Γ ( t ) = n + 1 i = 1 N E i ( t ) = i = 1 N E i ( t )

is H n -equivalent to spt V t for a.e. t 0 , and in fact

dim ( Γ ( t ) spt V t ) n - 1 for a.e.  t 0 .

Since V t 𝐈𝐕 n ( n + 1 ) for a.e. t 0 , there exists a countably n-rectifiable set Γ ~ ( t ) such that V t = 𝐯𝐚𝐫 ( Γ ~ ( t ) , θ t ) and θ t ( x ) = Θ n ( V t , x ) . By definition, we have n ( Γ ~ ( t ) spt V t ) = 0 , and by [33, Theorem 3.5] we have Θ n ( V t , x ) = 0 for n -a.e. x spt V t Γ ~ ( t ) . The latter claim shows that spt V t Γ ~ ( t ) is included in the set appearing in (5.7), and we can conclude that Γ ~ ( t ) is n -equivalent to spt V t . This proves Theorem 2.10 (iii) and (iv).

6 Two-sidedness at unit density point

In this section, we prove Theorem 2.11 (vii) and (viii), which are restated in Propositions 6.3 and 6.4, respectively. To do so, we need to analyze the behavior of approximating flows. The first Lemma 6.1 shows that, for a.e. t, only the reduced boundaries of approximating grains contribute to the limit measure and that the measures coming from “interior boundaries” E j , k ( t ) * E j , k ( t ) vanish in the limit. Roughly speaking, this is due to the measure minimizing property in the length scale of o ( 1 / j 2 ) which has the effect of eliminating the interior boundaries.

Lemma 6.1.

For a.e. t R + and E j ( t ) = { E j , k ( t ) } k = 1 N in Theorem 3.6, we have

(6.1) lim k = 1 N * E j , k ( t ) = 2 μ t .

Proof.

We fix t such that

(6.2) lim Δ j vc j ( t ) ( Ω ) j Δ t j = 0 ,

which holds for a.e. t + due to (3.7), and we drop t for simplicity in the following. It is sufficient to prove that

lim n ( U R spt j k = 1 N * E j , k ) = 0

for arbitrary R 1 since

spt j k = 1 N E j , k ,
n ( k = 1 N E j , k spt j ) = 0 ,
2 j = k = 1 N * E j , k + 2 n spt j k = 1 N * E j , k .

As in [22, Section 4], let r := 1 / ( j ) 5 / 2 . Define

Z := { x U R spt j : j ( B r ( x ) ) c 3 r n } and Z c := U R spt j Z ,

where c 3 is the same constant appearing in [21, Proposition 7.2] and apply it to each ball B r ( x ^ ) with x ^ Z to estimate j ( B r / 2 ( x ^ ) ) : there exists c 4 = c 4 ( n ) as in the claim and an j -admissible function f and r [ r / 2 , r ] such that the following conditions hold:

  1. f ( x ) = x for x n + 1 U r ( x ^ ) .

  2. f ( x ) B r ( x ^ ) for x B r ( x ^ ) .

  3. It holds

    f j ( B r ( x ^ ) ) 1 2 j ( B r ( x ^ ) ) .

  4. For all k,

    n + 1 ( E j , k E ~ j , k ) c 4 ( j ( B r ( x ^ ) ) ) n + 1 n ,

    where { E ~ j , k } k = 1 N = f j .

We may use Lemma 3.5 with C = B r ( x ^ ) and above (i)–(iv) to check that f 𝐄 vc ( j , j ) as long as j is large enough (actually, if exp ( - j / 2 r ) = exp ( - j 3 / 2 / 2 ) < 1 2 ). Then, by the definition of Δ j vc j ( Ω ) , (i) and (iii) as well as max B r ( x ^ ) Ω exp ( 2 c 1 r ) min B r ( x ^ ) Ω , we have

Δ j vc j ( Ω ) f j ( Ω ) - j ( Ω )
( min B r ( x ^ ) Ω ) { exp ( 2 c 1 r ) f j ( B r ( x ^ ) ) - j ( B r ( x ^ ) ) }
( min B r ( x ^ ) Ω ) ( 1 2 exp ( 2 c 1 r ) - 1 ) j ( B r ( x ^ ) )
- min B 2 R Ω 4 j ( B r ( x ^ ) )
(6.3) - min B 2 R Ω 4 j ( B r / 2 ( x ^ ) )

for all sufficiently large j . By the Besicovitch covering theorem, we have a mutually disjoint set of closed balls { B r / 2 ( x ^ i ) } x ^ i Z (whose number is at most c ( n ) r - n - 1 ) such that

(6.4) j ( Z ) 𝐁 ( n ) j ( B r / 2 ( x ^ i ) ) ,

and (6.3) and (6.4) show that

(6.5) j ( Z ) - c ( n ) r n + 1 min B 2 R Ω Δ j vc j ( Ω ) .

Since - Δ j vc j ( Ω ) r n + 1 due to (6.2), the right-hand side of (6.5) converges to 0 as . By (6.5), we need to prove

(6.6) lim n ( Z c k = 1 N * E j , k ) = 0 .

By the Besicovitch covering theorem, we have a set of mutually disjoint balls { B r ( x i ) } with

x i Z c k = 1 N * E j , k

such that

(6.7) n ( Z c k = 1 N * E j , k ) 𝐁 ( n ) i n ( B r ( x i ) Z c k = 1 N * E j , k ) .

Because of the lower bound of the measure j ( B r ( x i ) ) c 3 r n for x i Z c , the number of disjoint balls is bounded by c 3 - 1 r - n j ( B 2 R ) . If we prove that

(6.8) lim sup i r - n n ( B r ( x i ) Z c k = 1 N * E j , k ) = 0 ,

combined with (6.7), we would prove (6.6), ending the proof. Assume for a contradiction that (6.8) were not true and we had a subsequence (denoted by the same index) x Z c such that

(6.9) 0 < α r - n n ( B r ( x ) Z c k = 1 N * E j , k ) .

Consider a rescaling F : n + 1 n + 1 defined by F ( x ) = ( x - x ) / r and consider sequences V := ( F ) ( j ) and F ( E j , k ) ( k = 1 , , N ). This rescaling was discussed in [22, Section 4], and there exist a subsequence (denoted by the same index) and a limit V 0 with the properties stated in [22, Theorem 4.1], which shows that V is a unit density varifold and spt V is a real-analytic minimal hypersurface away from a closed lower-dimensional singularity. Moreover, spt V is two-sided in the following sense. Let E 1 , , E N n + 1 be limits of F ( E j , 1 ) , , F ( E j , N ) (see [22, p. 517]). By the lower-semicontinuity of the BV function, we have * E k V , and { E 1 , , E N } can be defined so that they are mutually disjoint open sets such that n + 1 spt V = k = 1 N E k . The two-sidedness means that, at each regular point x of spt V , there are two distinct indices k 1 , k 2 { 1 , , N } such that x clos E k 1 clos E k 2 (see [22, Lemma 4.8]). In particular, this implies that

(6.10) V = 1 2 k = 1 N * E k .

On the other hand, (6.9) implies that

V ( U 2 ) = V ( U 2 k = 1 N * ( F ( E j , k ) ) ) + V ( U 2 k = 1 N * ( F ( E j , k ) ) )
r - n n ( U 2 r ( x ) Z c k = 1 N * E j , k ) + V ( U 2 k = 1 N * ( F ( E j , k ) ) )
(6.11) α + 1 2 k = 1 N * ( F ( E j , k ) ) ( U 2 ) .

Since

(6.12) * E k ( U 2 ) lim inf * ( F ( E j , k ) ) ( U 2 )

for each k, (6.10)–(6.12) show

V ( U 2 ) lim inf V ( U 2 ) - α V ( B 2 ) - α .

Since V ( U 2 ) = 0 , this is a contradiction. The argument up to this point shows (6.8), which in turn shows the claim (6.1). ∎

The next Lemma 6.2 shows that j ( t ) is locally and subsequentially close to “θ-layered sheets” after appropriate blow-ups, for almost all times and places.

Lemma 6.2.

For L 1 -a.e. t R + and H n -a.e. x spt V t , there exist θ := Θ n ( V t , x ) N ,

T := Tan ( V , x ) 𝐆 ( n + 1 , n ) ,

r 0 + , a subsequence { j } = 1 { j } = 1 and H n -measurable sets W T B r with the following property (after a change of variable, we may assume that x = 0 in the following). Define f ( r ) ( y ) := y / r for y R n + 1 and { E j , 1 , , E j , N } := E j ( t ) . Then we have

(6.13) lim ( f ( r ) ) j ( t ) = 𝐯𝐚𝐫 ( T , θ )

as varifolds,

(6.14) 0 ( B r T - 1 ( a ) i = 1 N E j , i ) = θ

for all a W , and

(6.15) lim sup a W { | T ( x r ) | : x B r T - 1 ( a ) i = 1 N E j , i } = 0 ,
(6.16) lim n ( W ) ω n r n = 1 .

The claims (6.14) and (6.15) also hold with * E j , i in place of E j , i .

Proof.

Without loss of generality, assume that t + satisfies

(6.17) lim inf ( n + 1 | Φ ε j δ ( j ( t ) ) | 2 Ω Φ ε j j ( t ) + ε j Ω - 1 𝑑 x - 1 Δ t j Δ j vc j ( t ) ( Ω ) ) < ,

which holds for a.e. t + by (3.7) and Fatou’s lemma. By the compactness theorem of [21, Theorem 8.6], we may conclude that there exists a converging subsequence in the sense of varifolds { j ( t ) } = 1 (denoted by the same index) and the limit V t 𝐈𝐕 n ( n + 1 ) with μ t = V t . By (3.10), each { E j , i ( t ) } = 1 also converges to E i ( t ) in L loc 1 ( n + 1 ) for i = 1 , , N .

By Corollary 5.4, for a.e. t + , we have

θ t n i = 1 N E i ( t ) = θ t n spt V t = V t ,

where θ t ( x ) := Θ n ( V t , x ) is n -a.e. integer-valued. Note in particular that spt V t as well as i = 1 N E i ( t ) are countably n-rectifiable. We fix such a generic t and subsequently drop the dependence of t.

In the following, we use the same argument in the proof of [21, Theorem 8.6]. Let { j } = 1 be a subsequence (denoted by the same index) such that the quantity in (6.17) is uniformly bounded. Let U n + 1 be a bounded open set. It is sufficient to prove the claim on U. As in [21, pp. 112–113], for each j , q , let A j , q be a set consisting of all x clos U such that

δ ( Φ ε j j ) ( B r ( x ) ) q Φ ε j j ( B r ( x ) )

for all r ( j - 2 , 1 ) , and additionally define

A q := { x clos U : there exist  x A j , q  for infinitely many   with  x x } .

Following the argument in [21, p. 113], one can prove that

V ( U q = 1 A q ) = 0 .

Thus, for n -a.e. x spt V , we have some q such that x A q , and additionally the approximate tangent space exists with multiplicity θ . Without loss of generality, we may assume that x = 0 and write T := Tan ( V , x ) . Since 0 A q , there exists a further subsequence of { j } = 1 (denoted by the same index) such that x j A j , q with lim x j = 0 . Set r = 1 and define f ( r ) ( x ) := x / r and

V j := ( f ( r ) ) j , V ~ j := ( f ( r ) ) ( Φ ε j j ) , V j * := ( f ( r ) ) * j ,

where * j denotes the unit density varifold defined from i = 1 N * E j , i . We may choose a further subsequence with the following properties:

(6.18) lim V j = lim V ~ j = 𝐯𝐚𝐫 ( T , θ ) ,
lim x j r = 0 , lim j - 1 r = 0 ,
(6.19) lim r - n n ( B r j * j ) = 0 .

Note that the choice with (6.19) is possible due to (6.1). We then proceed verbatim as in [21, pp. 114–120] with ν = θ + 1 and d = θ . In particular, we fix λ ( 1 , 2 ) such that λ n + 1 θ < θ + 1 (as in [21, (8.111)]) and we use [21, Lemma 8.5]. In summary, with the same notation as in [21], we obtain a sequence of sets G * * j B ( λ - 1 ) r with the following properties for all large enough :

(6.20) lim r - n j ( B ( λ - 1 ) r G * * ) = 0 ,
(6.21) lim sup { | T ( x r ) | : x G * * } = 0 ,
(6.22) sup a T B ( λ - 1 ) r 0 ( { x G * * : T ( x ) = a } ) θ .

These claims are, respectively, [21, (8.154), (8.156) (stated differently) and (8.159)]. The measurability of G * * is not stated in [21]. However, since each step to estimate j ( B ( λ - 1 ) r G * * ) uses covering arguments, if necessary, we may simply take the complement of these coverings with no change of estimates and obtain a possibly smaller G * * which is a Borel set. We next prove that, writing g ~ ( a ) := 0 ( { x G * * : T ( x ) = a } ) for a T , we have

(6.23) lim n ( { a T B ( λ - 1 ) r : g ~ ( a ) = θ } ) ω n ( λ - 1 ) n r n = 1 .

Note that g ~ ( a ) as above on T is n -measurable (see [11, Lemma 5.8]). To see (6.23), first, by (6.18) and (6.20), we have

(6.24) lim ( f ( r ) ) ( j G * * × 𝐆 ( n + 1 , n ) ) = 𝐯𝐚𝐫 ( T B λ - 1 , θ )

as varifolds. We also have

(6.25) T ( f ( r ) ) ( j G * * × 𝐆 ( n + 1 , n ) ) ( B λ - 1 ) = ( r ) - n T g ~ ( a ) 𝑑 n ( a ) .

Since T commutes with lim and T 𝐯𝐚𝐫 ( T B λ - 1 , 1 ) = 𝐯𝐚𝐫 ( T B λ - 1 , 1 ) , (6.24) and (6.25) show that

(6.26) θ ω n ( λ - 1 ) n = lim ( r ) - n T g ~ ( a ) 𝑑 n ( a ) .

Since g ~ θ due to (6.22) (note also G * * B ( λ - 1 ) r ), equation (6.26) shows

(6.27) 0 = lim ( r ) - n T B ( λ - 1 ) r ( θ - g ~ ) d n lim ( r ) - n n ( T B ( λ - 1 ) r { g ~ θ - 1 } ) .

Then (6.22) and (6.27) show (6.23). Finally, we define

(6.28) W := { a T : g ~ ( a ) = θ } ( T ( j B ( λ - 1 ) r G * * ) T ( j B ( λ - 1 ) r * j ) ) .

Since g ~ is n -measurable, so is W . Due to (6.19), (6.20) and (6.23), we can deduce (6.16). For any a W , T - 1 ( a ) j B ( λ - 1 ) r consists of θ points belonging to G * * * j by (6.28). This proves (6.14) (both with j and * j ). This combined with (6.21) also proves (6.15), and we may conclude the proof after renaming ( λ - 1 ) r as r . ∎

Proposition 6.3.

For V t and { E i ( t ) } i = 1 N in Proposition 3.6, for a.e. t R + , we have

(6.29) n ( { x : θ t ( x ) = 1 } i = 1 N * E i ( t ) ) = 0 .

Proof.

We prove (6.29) for a.e. t such that the conclusion of Lemma 6.2 holds. With such t fixed, we drop t and suppose that (6.29) were not true for a contradiction. We may apply the result of Lemma 6.2 and find a point x spt V i = 1 N * E i with Θ n ( V , x ) = 1 and T := Tan ( V , x ) . Moreover, by the well-known property of the set of finite perimeter, there exists some k { 1 , , N } such that

lim r 0 n + 1 ( E k B r ( x ) ) ω n + 1 r n + 1 = 1 .

Without loss of generality, we may assume that x = 0 , T = { x n + 1 = 0 } and k = 1 . Since χ E j , 1 χ E 1 in L loc 1 as , in choosing the subsequence (denoted by the same index) in Lemma 6.2, we may additionally arrange the choice so that we have

(6.30) lim n + 1 ( E j , 1 B r ) ω n + 1 r n + 1 = 1 .

For simplicity, write the union of the measure-theoretic boundaries as * j := i = 1 N * E j , i . Since j is closed, we remind the reader that

(6.31) * j * j j .

We also define the measure-theoretic interior and exterior of E j , 1 (see [11, Definition 5.13]) by

I := { x n + 1 : lim r 0 n + 1 ( B r ( x ) E j , 1 ) r n + 1 = 0 }

and

O := { x n + 1 : lim r 0 n + 1 ( B r ( x ) E j , 1 ) r n + 1 = 0 } .

We note from the definition that

(6.32) I * j = .

We next use some results in the proof of [11, Theorem 5.23] on the property of measure-theoretic interior and exterior. For each m , k , define

(6.33) { G ( k ) := { x n + 1 : n + 1 ( B r ( x ) O ) ω n r n + 1 3 n + 2  for  0 < r < 3 k } , H ( k ) := { x n + 1 : n + 1 ( B r ( x ) I ) ω n r n + 1 3 n + 2  for  0 < r < 3 k }

and

(6.34) { G ± ( k , m ) := G ( k ) { x : x ± s e n + 1 O  for  0 < s < 3 m } , H ± ( k , m ) := H ( k ) { x : x ± s e n + 1 I  for  0 < s < 3 m } .

Here e n + 1 is the unit vector pointing towards the positive direction of the x n + 1 -axis. These sets have the property (see [11, step 3 of Theorem 5.23]) that

n ( T ( G ± ( k , m ) ) ) = n ( T ( H ± ( k , m ) ) ) = 0

for all k , m . Moreover, for all

a T k , m = 1 T ( G + ( k , m ) G - ( k , m ) H + ( k , m ) H - ( k , m ) )

and with 0 ( T - 1 ( a ) * E j , 1 ) < , if x 1 , x 2 T - 1 ( a ) with T ( x 1 ) < T ( x 2 ) , x 1 I and x 2 O , then there exists x 3 T - 1 ( a ) * E j , 1 such that T ( x 1 ) < T ( x 3 ) < T ( x 2 ) (see [11, step 5 of Theorem 5.23]). Here, x 1 is in the interior and x 2 is in the exterior of E j , 1 over a, and the claim is that there must be a “boundary point” x 3 between these two points. The same claim holds if x 1 O and x 2 I instead.

Let W be the set obtained in Lemma 6.2. By [11, Lemma 5.9], we have n + 1 ( I E j , 1 ) = 0 , so with (6.30) we have

(6.35) lim n + 1 ( I B r ) ω n + 1 r n + 1 = 1 .

Then, by the Fubini theorem and (6.35), we may choose a sequence { b } = 1 + such that b [ r / 3 , r / 2 ] and so that, writing

A + := B r { x n + 1 = b } ,
A - := B r { x n + 1 = - b } ,

we have

(6.36) lim n ( I A + ) n ( A + ) = lim n ( I A - ) n ( A - ) = 1 .

On cylinders T ( A + ) × [ - b , b ] n + 1 , by (6.33), (6.34) and the stated property thereafter, for n -a.e. a W , we have the following property:

(6.37) { if  ( a , s 1 ) I  and  ( a , s 2 ) O  with  - b s 1 < s 2 b , then there exists  s ^ ( s 1 , s 2 )  such that  ( a , s ^ ) * E j , 1 ,

and similarly if ( a , s 1 ) O and ( a , s 2 ) I . On the other hand, we know that T - 1 ( a ) B r * j is a singleton located close to T due to (6.14) and (6.15). Here, we used the fact that (6.14) is satisfied for both j and * j as well as (6.31). We use this fact to get

a W T ( I A + ) T ( I A - ) = : W * .

Note that ( a , r ) and ( a , - r ) are both in I , and ( { a } × [ - r , r ] ) * j is a singleton due to the way W is defined. If ( { a } × [ - r , r ] ) O , then (6.37) implies that there must be at least two points of * E j , 1 in { a } × [ - r , r ] , since both crossings from I to O and the other way around have to happen. Since * j = i = 1 N * E j , i , this is a contradiction. Combined with (6.32), we conclude that for n -a.e. a W * , { a } × [ - r , r ] is a disjoint union of one point of * j and two line segments included in I , with no point of O . Because of (6.16) and (6.36), one also sees

lim n ( W * ) n ( A + ) = 1 .

In particular, W * has a positive n measure in T n × { 0 } and there must be a Lebesgue point a of W * such that ( { a } × [ - r , r ] ) * j is a singleton, say, { ( a , s ) } . Then, by the Fubini theorem and the property of W * ,

r - ( n + 1 ) n + 1 ( O B r ( ( a , s ) ) ) r - ( n + 1 ) n + 1 ( ( B r n ( a ) W * ) × [ s - r , s + r ] )
2 r - n n ( B r n ( a ) W * ) ,

which converges to 0 as r 0 since a is a Lebesgue point of W * in T. Since

n + 1 ( O B r ( x ) ) = n + 1 ( B r ( x ) E j , 1 ) ,

this implies that ( a , s ) I . On the other hand, ( a , s ) * j , a contradiction to (6.32). This concludes the proof. ∎

Proposition 6.4.

Assume N = 2 . For V t and { E i ( t ) } i = 1 2 in Proposition 3.6, for a.e. t R + , we have

θ t ( x ) = { 𝑜𝑑𝑑 for  n -a.e.  x * E 1 ( t ) ( = * E 2 ( t ) ) , 𝑒𝑣𝑒𝑛 for  n -a.e.  x spt V t * E 1 ( t ) .

Proof.

The proof proceeds similarly to the proof of Proposition 6.3, except that we need to localize the argument to each layers. Fix a bounded open set U n + 1 . We may choose a generic t + as before and drop t. Since * E 1 = * E 2 by the definition of the reduced boundary, for n -a.e. x spt V * E 1 , we have

lim r 0 n + 1 ( E i B r ( x ) ) ω n + 1 r n + 1 = 1 for either  i = 1  or  2 .

On the other hand, for n -a.e. x * E 1 , there exists a unit outer normal ν to * E 1 such that, letting

B r + ( - ) ( x ) := { y B r ( x ) : ( y - x ) ν 0 (resp.  0 ) } ,

we have

lim r 0 n + 1 ( E 1 B r + ( x ) ) ω n + 1 r n + 1 = 0 and lim r 0 n + 1 ( E 1 B r - ( x ) ) ω n + 1 r n + 1 = 1 2 .

We use Lemma 6.2, and in the proof of Lemma 6.2 we may additionally assume that the chosen subsequence satisfies

(6.38) lim n + 1 ( E j , i B r ) ω n + 1 r n + 1 = 1 for either  i = 1  or  2

if 0 spt V * E 1 , and

(6.39) lim n + 1 ( E j , 1 B r + ) ω n + 1 r n + 1 = 0 and lim n + 1 ( E j , 1 B r - ) ω n + 1 r n + 1 = 1 2

if 0 * E 1 . Without loss of generality, we may assume i = 1 in (6.38), T = { x n + 1 = 0 } and B r + = B r { x n + 1 0 } in (6.39). By (6.13), we have

(6.40) lim ( r ) - n B 2 r S - T d ( j ) ( x , S ) = 0 .

Set

(6.41) C := { x j B r : Tan ( j , x ) - T 1 10 } .

Given any 0 < δ < r and n -a.e. x j , by the rectifiability of j , there exists 0 < r < δ such that

(6.42) 1 2 1 ω n r n j ( B r ( x ) ) and 1 ω n r n B r ( x ) Tan ( j , x ) - S d ( j ) 1 40 .

Then, for n -a.e. x j B r C , (6.41) and (6.42) show

B r ( x ) T - S d ( j ) Tan ( j , x ) - T j ( B r ( x ) ) - ω n r n 40
ω n r n ( 1 20 - 1 40 )
= ω n r n 40 .

We cover j B r C by such balls and use the Besicovitch covering theorem to show that

(6.43) n ( j B r C ) 40 𝐁 n + 1 B 2 r T - S d ( j ) ,

where 𝐁 n + 1 is the Besicovitch constant. Then (6.40) and (6.43) show

lim ( r ) - n n ( j B r C ) = lim ( r ) - n n ( T ( j B r C ) ) = 0 ,

so that T - 1 ( a ) j B r C for a in a large portion of T B r for large enough . Now recall the property of W in Lemma 6.2. We may redefine W by W T ( j B r C ) and keep the properties (6.14)–(6.16). By this, we additionally have that

Tan ( j , x ) - T 1 10 for all  x T - 1 ( W ) B r * j .

We now proceed similarly to the previous Proposition 6.3, and choose b satisfying (6.36) in the case of (6.38). The case of (6.39) can be handled similarly, so we discuss the former case. For all large , we may choose a Lebesgue point a of W in T such that ( a , - b ) and ( a , b ) are also Lebesgue points of I A - and I A + , respectively. By (6.14) and (6.15), there are - b := u 0 < u 1 < < u θ < u θ + 1 := b such that

k = 1 θ { ( a , u k ) } = T - 1 ( a ) B r * j .

At each point ( a , u k ) , since it is in * E j , 1 , the blow-up of E j , 1 converges to a half-space, with the approximate tangent space having a small slope relative to T due to (6.41). Then, for sufficiently small

0 < δ < min 0 k θ { | u k + 1 - u k | } and k = 1 , , θ ,

we may choose b , k [ δ / 3 , δ / 2 ] so that

n ( T - 1 ( B δ n ( a ) ) { x n + 1 = u k + b , k } O ) ω n δ n 1 - 1 6 θ

and

n ( T - 1 ( B δ n ( a ) ) { x n + 1 = u k - b , k } I ) ω n δ n 1 - 1 6 θ ,

or the inequalities replacing the role of O and I . We may also assume that

n ( B δ n ( a ) W ) ω n δ n 8 9

since a is a Lebesgue point of W , and similarly

n ( T - 1 ( B δ n ( a ) ) I A ± ) ω n δ n 8 9 .

With these properties, we can make sure that, with respect to n , one third of W B δ n ( a ) has the property that, if a ~ is in this set,

( a ~ , ± b ) I , ( a ~ , u k + b , k ) O , ( a ~ , u k - b , k ) I  or vice-versa for  k = 1 , , θ .

By using (6.37), for n -a.e. a ~ as above, there exists some s k ( u k - b , k , u k + b , k ) with ( a ~ , s k ) * E j , 1 for each k = 1 , , θ , and there are no other points of * E j , 1 along the line segment connecting ( a ~ , - b ) and ( a ~ , b ) . Looking at this line segment and the intersection of O and I , since these two endpoints are in I and each ( a ~ , s k ) is sided by O and I , θ has to be necessarily even. This finishes the proof in the case of (6.38). For (6.39), the similar argument results in the situation that ( a ~ , - b ) I and ( a ~ , b ) O , which necessitates that θ is odd. This concludes the proof. ∎

Finally, we comment on the proofs of Theorem 2.12 and 2.13. If V t is a unit density flow in U × ( t 1 , t 2 ) , then θ t ( x ) = 1 for V t -a.e. x U and a.e. t ( t 1 , t 2 ) . By Theorem 2.10 (ii) and (iv), we may assume that V t = 𝐯𝐚𝐫 ( Γ ( t ) , 1 ) , and by Theorem 2.11 (vii), Γ ( t ) may be replaced by i = 1 N * E i ( t ) . Thus, (2.13) follows immediately. To check that (2.7) holds, since (see (2.9) and (2.10))

i j n I i , j ( t ) = 2 n Γ ( t ) = 2 V t ,

the left-hand side of (2.7) is equal to 2 0 T δ V t ( g ) 𝑑 t . Since v i ν i = ( h ν i ) ν i = h for n -a.e. x I i , j ( t ) due to the perpendicularity of the mean curvature vector, the right-hand side of (2.7) is - 2 0 T h g d V t 𝑑 t . Since the generalized mean curvature vector exists for a.e. t > 0 , they are indeed equal. This proves the claim of Theorem 2.12. For Theorem 2.13, under the assumption, one can show that there exists

T 0 = T 0 ( n , n ( Γ 0 ) , r 0 , δ 0 ) > 0

such that (recall (4.28))

Γ 0 ρ ( y , s ) ( x , 0 ) 𝑑 n ( x ) < 2 - δ 0 2

for all y n + 1 and 0 < s T 0 . Then Huisken’s monotonicity formula shows that

n + 1 ρ ( y , s ) ( x , t ) d V t ( x )

is non-increasing on t [ 0 , s ) , and thus < 2 - δ 0 / 2 . For a contradiction, if V t is not unit density on [ 0 , T 0 ] , there would exist some t ( 0 , T 0 ) with V t 𝐈𝐕 n ( n + 1 ) , and y spt V t such that Θ n ( V t , y ) 2 and where T y V t exists. Then one can prove that

lim ϵ 0 + n + 1 ρ ( y , t + ϵ ) ( x , t ) d V t ( x ) = Θ n ( V t , y ) 2 .

Since t + ϵ < T 0 for all small ϵ > 0 , this would be a contradiction. Thus, V t is a unit density Brakke flow on [ 0 , T 0 ] . Once this is proved, by Theorem 2.12, the claim of the BV solution also follows. This is the outline of the proof of Theorem 2.13.

7 Final remarks

7.1 On generalized BV solutions

As explained in Definition 2.4, a BV flow is classically defined as consisting of two objects: N families of sets of finite perimeter E i ( t ) , and velocities v i . From them, naturally, one can define a unit density varifold

V t = 𝐯𝐚𝐫 ( i = 1 N * E i ( t ) , 1 )

and check that (2.7) implies that the generalized mean curvature h ( , V t ) is equal to v i ν i on * E i ( t ) for i = 1 , , N . The generalized BV flow of Theorem 2.11 (vi) involves the accompanying Brakke flow V t in addition to the families of sets of finite perimeter, and one may wonder if the definition makes sense even without the reference to the Brakke flow. In fact, it is interesting to observe that each * E i ( t ) for a.e. t > 0 is C 2 -rectifiable due to Menne’s C 2 -rectifiability theorem [28, Theorem 4.8], and one can define a unique second fundamental form for * E i ( t ) as well as a mean curvature vector by the C 2 -approximability property, independent of V t . The mean curvature vector defined in this sense coincides with h ( , V t ) V t -a.e. on * E i ( t ) . Thus, h ν i on * E i ( t ) is uniquely defined from the C 2 -rectifiability without reference to V t . On the other hand, the summation

h ~ := 1 2 i = 1 N ( h ν i ) ν i

may not correspond, in general, to the generalized mean curvature vector of 𝐯𝐚𝐫 ( i = 1 N * E i ( t ) , 1 ) if there is some non-trivial higher multiplicity portion of V t . For example on 2 , define

E + := { ( x , y ) : y > 0  if  x 0 , y > x 2  if  x > 0 } and E - := { ( x , y ) : y < 0  if  x 0 , y < - x 2  if  x > 0 } .

Then V := 𝐯𝐚𝐫 ( E + , 1 ) + 𝐯𝐚𝐫 ( E - , 1 ) has a bounded generalized mean curvature, while 𝐯𝐚𝐫 ( * E + * E - , 1 ) has a singular first variation at the origin. Note that V has multiplicity = 2 on the negative x-axis. In this sense, formula (2.7) does not hold in general. Note that (2.8) is relevant only when n ( Γ 0 ) is finite, and it follows from Theorem 2.7 (iv) with v i = h ν i . Overall, for the generalized BV solution, it makes sense to consider the pair of sets of finite perimeter and Brakke flow together, unlike the original BV solutions discussed in Definition 2.4.

7.2 MCF with fixed boundary conditions

In [35], given a strictly convex bounded domain U n + 1 with C 2 boundary U , a countably n-rectifiable set Γ 0 U with n ( Γ 0 ) < and an open partition E 0 , 1 , , E 0 , N of U such that Γ 0 = U i = 1 N E 0 , i , existence of a Brakke flow and a family of open partitions with fixed boundary condition is established for the given initial datum. The construction method is along the lines of [21], and we may also carry it out using the volume-controlled Lipschitz maps. If one compares the construction in [35] with that in [21], one sees that differences occur only near the boundary U : more precisely, the approximate smoothed mean curvature vector is damped near the portion of Γ 0 close to U , and there is another step in each epoch – a Lipschitz retraction step (see [35, Section 2.6]). Hence, the proof of the present paper works with no essential change away from U , and (2.12) holds for ϕ C c 1 ( U × + ) ; since the formula does not involve ϕ , by approximation, the same formula holds even for ϕ C 1 ( clos U × [ 0 , T ] ) for arbitrary T > 0 . Since the existence results in [34] are based on [35], the same applies to the solutions discussed in [34].

7.3 A lower bound estimate for extinction time

Proof of Theorem 2.14.

Suppose that n ( Γ 0 ) < , and assume without loss of generality that N is the index of the only grain E 0 , i with infinite volume. Define then the open set E ( t ) := i = 1 N - 1 E i ( t ) . Since ν i h = - ν j h on * E i ( t ) * E j ( t ) for i j n -a.e., after summing over i = 1 , , N - 1 , formula (2.12) gives

| E ( t ) | - | E ( 0 ) | = 0 t * E ( s ) ν h 𝑑 n

for all 0 < t < , where ν is the outer unit normal of * E ( s ) . In particular, if we set v ( t ) := | E ( t ) | , then we have

(7.1) v ( t ) = * E ( t ) ν h 𝑑 n for a.e.  t + .

Next, set a ( t ) := V t ( n + 1 ) . By Brakke’s inequality (2.1), the upper derivative a + ( t ) of a ( t ) satisfies

(7.2) - a + ( t ) | h | 2 d V t Γ ( t ) | h | 2 𝑑 n * E ( t ) | h | 2 𝑑 n for a.e.  t + .

Combining (7.1) and (7.2) then gives the inequality

- v ( t ) ( n ( * E ( t ) ) ) 1 2 ( * E ( t ) | h | 2 𝑑 n ) 1 2 ( - a ( t ) a + ( t ) ) 1 2 = ( - ( a 2 ) + ( t ) 2 ) 1 2

for a.e. t + . Integrating over any interval [ 0 , T ] and using the Cauchy–Schwarz inequality then yields

(7.3) | E ( 0 ) | - | E ( T ) | T 2 ( V 0 ( n + 1 ) 2 - V T ( n + 1 ) 2 ) 1 2 T 2 n ( Γ 0 ) .

Since V t 0 as long as E ( t ) , the extinction time T * is at least equal to the first time when E ( t ) becomes the empty set, so that (7.3) implies

T * 2 ( | E ( 0 ) | n ( Γ 0 ) ) 2 ,

that is, (2.14). ∎


Communicated by Frank Duzaar


Award Identifier / Grant number: 18H03670

Award Identifier / Grant number: 19H00639

Award Identifier / Grant number: 17H01092

Funding statement: The first author is partially supported by the GNAMPA – Gruppo Nazionale per l’Analisi Matematica, la Probabilità e le loro Applicazioni of INdAM. The second author is partially supported by JSPS 18H03670, 19H00639, 17H01092.

A The existence theorem of Kim and Tonegawa revisited

Here we point out places which require a change to 𝐄 vc ( , j ) from 𝐄 ( , j ) in [21, 22]. It turned out that the proofs require no essential change and the only point to be checked is that the same Lipschitz maps used in the proofs satisfy the condition of Definition 3.1 (ii).

A.1 Construction of approximate flows

For the construction of discrete approximate sequence in [21, Section 6], we simply replace 𝐄 ( j , l , j ) by 𝐄 vc ( j , l , j ) and Δ j j , l ( Ω ) by Δ j vc j , l ( Ω ) when f 1 is chosen in [21, (6.9)]. As in [21, (6.10)], if we define

{ E j , l , i } i = 1 N := j , l ,
j , l + 1 * := ( f 1 ) j , l ,
{ E j , l + 1 , i * } i = 1 N = * j , l + 1 ,

by Definition 3.1 (ii) and (3.1), we have for each i = 1 , , N ,

n + 1 ( E j , l + 1 , i * E j , l , i ) { j , l ( Ω ) - j , l + 1 * ( Ω ) } j - ( Δ j vc j , l ( Ω ) ) j .

The change in [21, (6.9)] is also reflected in the estimate [21, (6.4)] and we have

j , l ( Ω ) - j , l - 1 ( Ω ) Δ t j + 1 4 n + 1 | Φ ε j δ ( j , l ) | 2 Ω Φ ε j j , l + ε j Ω - 1 𝑑 x - ( 1 - j - 5 ) Δ t j Δ j vc j , l - 1 ( Ω )
ε j 1 8 + c 1 2 2 j , l - 1 ( Ω ) .

This leads to estimates (3.7) and (3.8).

A.2 Proofs of rectifiability and integrality

For the proofs of rectifiability and integrality of μ t for 1 -a.e. t + , the smallness of Δ j ( t ) is essential in [21, Sections 7 and 8]. The general idea is that, whenever Δ j ( t ) is used in [21], the proofs use contradiction arguments and some appropriate Lipschitz deformations with drastic measure reduction are constructed. Here, “drastic” means that the measure is typically reduced by some factor of the measure itself, so that the reduction is typically much larger than the volume change caused by the deformation. Thus, even if we impose the additional condition, it is satisfied by the same Lipschitz deformations in [21] and we only need to check that it is indeed the case. We point out the following three separate places.

(i) In [21, Proposition 7.2], it is proved that there exists an -admissible function f which reduces the measure ( B r ) by the factor of 1 2 (see (iii)) for some r [ R 2 , R ] when the measure in B R is sufficiently small. Note that (iv) gives the desired estimate on the change of volume of each grain in terms of ( B r ) . Since the radii r of balls used later are typically O ( 1 / j 2 ) , and thus the volume change is O ( r n + 1 ) = O ( r n / j 2 ) , and ( B r ) = O ( r n ) , we can deduce that f belongs to 𝐄 vc ( , j ) . Since the claim of [21, Proposition 7.2] is about the existence of an -admissible function, no change is required in the proof.

(ii) In [21, Theorem 7.3], assumption (iv) should be replaced by the volume-controlled counterpart. In the proof, on [21, pp. 94–95], a Lipschitz map f is defined with the properties stated in the bottom of [21, p. 94] using [21, Proposition 7.2]. Using [21, (7.40)], which gives

( 1 - 2 - 1 2 ) j l Ω ( B r x ( x ) ) j l Ω ( B r x ( x ) ) - ( f x ) j l Ω ( B r x ( x ) ) ,

one can proceed in [21, (7.43)] as

n + 1 ( E i E ~ i ) c 4 k = 1 Λ ( j l ( B ( k ) ) n + 1 n
c 4 c 3 1 n 2 j l 2 k = 1 Λ j l ( B ( k ) )
c 4 c 3 1 n ( min B 3 ( x 0 ) Ω ) - 1 2 ( 1 - 2 - 1 2 ) j l 2 k = 1 Λ ( j l Ω ( B ( k ) ) - f j l Ω ( B ( k ) ) )
= c 4 c 3 1 n ( min B 3 ( x 0 ) Ω ) - 1 2 ( 1 - 2 - 1 2 ) j l 2 ( j l ( Ω ) - f j l ( Ω ) ) .

Thus, for all sufficiently large l, Definition 3.1 (ii) is satisfied. The use of [21, Lemma 4.12] is also justified, and we have f 𝐄 vc ( j l , j l ) .

(iii) Throughout [21, Section 8], the only crucial point that needs to be checked is in [21, Lemma 8.1], which involves the actual construction of a measure-reducing Lipschitz deformation. It proves roughly that when is flat and close to being measure-minimizing within a cylinder of size O ( 1 / j 2 ) , then the measure has to be an integer multiple of discs. The argument proceeds by assuming the contrary. The intuitive picture is that, if does not have a measure close to a multiple of discs, then one can locate a hole which can be expanded horizontally. This would cause a drastic reduction of the measure and lead to a contradiction to the almost measure-minimizing property. More precisely, besides the change of Δ j to Δ j vc throughout, near the end of the proof of [21, Lemma 8.1, p. 106], one needs to check that the “expansion map” f a is in 𝐄 vc ( , E ( r 1 , ρ 1 ) , j ) . We recall that Y T is a set of ν points with

diam Y < j - 2 ([21, Lemma 8.1(3)]),
r 1 < R < j - 2 2 ([21, four lines above (8.7)] and [21, Lemma 8.1(1)]),
ρ 1 = ( 1 + R - 1 r 1 ) ρ < 2 ρ < j - 2 ([21, (8.7)] and [21, Lemma 8.1(1)]),

and, with T 𝐆 ( n + 1 , n ) fixed,

E ( r 1 , ρ 1 ) = { x n + 1 : | T ( x ) | r 1 , dist ( T ( x ) , Y ) ρ 1 } ([21, (8.1)]).

The map f a is defined to be the identity map outside of E ( r 1 , ρ 1 ) , so the change of volume of grains caused by f a is at most

(A.1) n + 1 ( E ( r 1 , ρ 1 ) ) 2 ω n ν r 1 n ρ 1 < 2 ω n ν r 1 n j 2 .

As one can see in [21, (8.67)], the reduction of measure is

(A.2) ( f a ) ( E ( r 1 , ρ 1 ) ) - ( E ( r 1 , ρ 1 ) ) < - 1 2 ( 1 - ζ ) ω n r 1 n ,

and we also have from [21, (8.8)] that

(A.3) ( E ( r 1 , ρ 1 ) ) = ( ν - ζ ) ω n r 1 n .

We need to see the difference with the weight Ω, and since diam E ( r 1 , ρ 1 ) < 4 / j 2 , we have

( f a ) ( Ω ) - ( Ω )
= ( f a ) Ω ( E ( r 1 , ρ 1 ) ) - Ω ( E ( r 1 , ρ 1 ) )
( max E ( r 1 , ρ 1 ) Ω ) ( f a ) ( E ( r 1 , ρ 1 ) ) - ( min E ( r 1 , ρ 1 ) Ω ) ( E ( r 1 , ρ 1 ) )
( min E ( r 1 , ρ 1 ) Ω ) ( e 4 c 1 j 2 ( f a ) ( E ( r 1 , ρ 1 ) ) - ( E ( r 1 , ρ 1 ) ) )
( min E ( r 1 , ρ 1 ) Ω ) { e 4 c 1 / j 2 ( ( f a ) ( E ( r 1 , ρ 1 ) ) - ( E ( r 1 , ρ 1 ) ) ) + ( e 4 c 1 / j 2 - 1 ) ( E ( r 1 , ρ 1 ) }
(A.4) - 1 2 ( 1 - ζ ) ω n ( min E ( r 1 , ρ 1 ) Ω ) e 4 c 1 / j 2 r 1 n + 4 c 1 j 2 e 4 c 1 / j 2 ( ν - ζ ) ω n r 1 n .

In the last line, we used (A.2) and (A.3). Note that the first term of the last line is a negative term of order O ( r 1 n ) , while the change of volume expressed in (A.1) is O ( r 1 n / j 2 ) . Thus, (A.1) and (A.4) give the desired inequality

n + 1 ( E i E ~ i ) n + 1 ( E ( r 1 , ρ 1 ) ) ( ( Ω ) - ( f a ) ( Ω ) ) j

for all sufficiently large j, and we have f a 𝐄 vc ( , E ( r 1 , ρ 1 ) , j ) . The rest of the proof is not affected by the change of volume-controlled deformation.

A.3 Volume change of grains

The motivation of having n + 1 ( E i E ~ i ) < 1 j in [21] is the use in the proof of [21, Lemma 10.10], and in fact it is the only place that this inequality is essentially used to derive any conclusion. In the proof, see the second line from the bottom of [21, p. 134], it is used to make sure that the volume change of grains is small for each discrete time step and the continuity of the labelling of each grain is derived in the end. The similar smallness of volume change is available with the volume-controlled counterparts since j ( t ) ( Ω ) is uniformly bounded for a fixed time interval [ 0 , T ] by (3.6). Thus, the proof can be carried out similarly.

A.4 Changes in the later paper by Kim and Tonegawa

The results from [22] are used in the present paper and the modifications are needed there as well. On the other hand, similarly, one can check that the change to the volume-controlled counterpart does not cause any difficulties. The part which is relevant to the change is [22, Section 4]. In the proof of [22, Lemma 4.2], a Lipschitz retraction map F ^ is used, with the reduction of measure inside of B r R ( z ( ) ) being β ( r R ) n , while the volume of the ball is O ( ( r R ) n + 1 ) (see [22, (4.5) and (4.6)]). Since r = 1 / j 5 / 2 (see [22, just after (4.2)]), one can check that F ^ 𝐄 vc ( j , j ) for all large . The similar argument can be applied to the proofs of [22, Lemmas 4.3 and 4.5]. The proof of [22, Lemma 4.6] uses [21, Lemma 8.1], with the modifications discussed above in Section A.2 (iii). In the proofs of [22, Lemmas 4.7–4.9], the Lipschitz maps reduce the measure in similar manners, and they belong to the volume-controlled counterparts. In particular, all of the results in [22] hold true even with the modifications.

References

[1] W. K. Allard, On the first variation of a varifold, Ann. of Math. (2) 95 (1972), 417–491. 10.2307/1970868Search 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] L. Ambrosio, N. Fusco and D. Pallara, Functions of Bounded Variation and Free Discontinuity Problems, Oxford Math. Monogr., Oxford University, New York, 2000. 10.1093/oso/9780198502456.001.0001Search in Google Scholar

[4] G. Bellettini and S. Y. Kholmatov, Minimizing movements for mean curvature flow of partitions, SIAM J. Math. Anal. 50 (2018), no. 4, 4117–4148. 10.1137/17M1159294Search in Google Scholar

[5] L. Bertini, P. Buttà and A. Pisante, Stochastic Allen-Cahn approximation of the mean curvature flow: Large deviations upper bound, Arch. Ration. Mech. Anal. 224 (2017), no. 2, 659–707. 10.1007/s00205-017-1086-3Search in Google Scholar

[6] K. A. Brakke, The Motion of a Surface by its Mean Curvature, Math. Notes 20, Princeton University, Princeton, 1978. Search in Google Scholar

[7] L. Bronsard and F. Reitich, On three-phase boundary motion and the singular limit of a vector-valued Ginzburg–Landau equation, Arch. Ration. Mech. Anal. 124 (1993), no. 4, 355–379. 10.1007/BF00375607Search in Google Scholar

[8] Y. G. Chen, Y. Giga and S. Goto, Uniqueness and existence of viscosity solutions of generalized mean curvature flow equations, J. Differential Geom. 33 (1991), no. 3, 749–786. 10.4310/jdg/1214446564Search in Google Scholar

[9] D. Depner, H. Garcke and Y. Kohsaka, Mean curvature flow with triple junctions in higher space dimensions, Arch. Ration. Mech. Anal. 211 (2014), no. 1, 301–334. 10.1007/s00205-013-0668-ySearch in Google Scholar

[10] S. Esedoḡlu and F. Otto, Threshold dynamics for networks with arbitrary surface tensions, Comm. Pure Appl. Math. 68 (2015), no. 5, 808–864. 10.1002/cpa.21527Search in Google Scholar

[11] L. C. Evans and R. F. Gariepy, Measure Theory and Fine Properties of Functions, Textb. Math., CRC Press, Boca Raton, 2015. 10.1201/b18333Search in Google Scholar

[12] L. C. Evans and J. Spruck, Motion of level sets by mean curvature. I, J. Differential Geom. 33 (1991), no. 3, 635–681. 10.4310/jdg/1214446559Search in Google Scholar

[13] J. Fischer, S. Hensel, T. Laux and T. M. Simon, The local structure of the energy landscape in multiphase mean curvature flow: Weak-strong uniqueness and stability of evolutions, preprint (2020), https://arxiv.org/abs/2003.05478. Search in Google Scholar

[14] A. Freire, Mean curvature motion of triple junctions of graphs in two dimensions, Comm. Partial Differential Equations 35 (2010), no. 2, 302–327. 10.1080/03605300903419775Search in Google Scholar

[15] Y. Giga and K. Yama-uchi, On a lower bound for the extinction time of surfaces moved by mean curvature, Calc. Var. Partial Differential Equations 1 (1993), no. 4, 417–428. 10.1007/BF01206960Search in Google Scholar

[16] S. Hensel and T. Laux, A new varifold solution concept for mean curvature flow: Convergence of the Allen–Cahn equation and weak-strong uniqueness, preprint (2021), https://arxiv.org/abs/2109.04233. Search in Google Scholar

[17] S. Hensel and T. Laux, Weak-strong uniqueness for the mean curvature flow of double bubbles, preprint (2021), https://arxiv.org/abs/2108.01733. 10.4171/IFB/484Search in Google Scholar

[18] T. Ilmanen, Convergence of the Allen–Cahn equation to Brakke’s motion by mean curvature, J. Differential Geom. 38 (1993), no. 2, 417–461. 10.4310/jdg/1214454300Search in Google Scholar

[19] T. Ilmanen, Elliptic regularization and partial regularity for motion by mean curvature, Mem. Amer. Math. Soc. 108 (1994), no. 520, 1–90. 10.1090/memo/0520Search in Google Scholar

[20] K. Kasai and Y. Tonegawa, A general regularity theory for weak mean curvature flow, Calc. Var. Partial Differential Equations 50 (2014), no. 1–2, 1–68. 10.1007/s00526-013-0626-4Search in Google Scholar

[21] L. Kim and Y. Tonegawa, On the mean curvature flow of grain boundaries, Ann. Inst. Fourier (Grenoble) 67 (2017), no. 1, 43–142. 10.5802/aif.3077Search in Google Scholar

[22] L. Kim and Y. Tonegawa, Existence and regularity theorems of one-dimensional Brakke flows, Interfaces Free Bound. 22 (2020), no. 4, 505–550. 10.4171/IFB/448Search in Google Scholar

[23] T. Laux and F. Otto, Convergence of the thresholding scheme for multi-phase mean-curvature flow, Calc. Var. Partial Differential Equations 55 (2016), no. 5, Article ID 129. 10.1007/s00526-016-1053-0Search in Google Scholar

[24] T. Laux and T. M. Simon, Convergence of the Allen–Cahn equation to multiphase mean curvature flow, Comm. Pure Appl. Math. 71 (2018), no. 8, 1597–1647. 10.1002/cpa.21747Search in Google Scholar

[25] 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

[26] F. Maggi, Sets of Finite Perimeter and Geometric Variational Problems, Cambridge Stud. Adv. Math. 135, Cambridge University, Cambridge, 2012. 10.1017/CBO9781139108133Search in Google Scholar

[27] C. Mantegazza, M. Novaga, A. Pluda and F. Schulze, Evolution of networks with multiple junctions, preprint (2018), https://arxiv.org/abs/1611.08254. Search in Google Scholar

[28] U. Menne, Second order rectifiability of integral varifolds of locally bounded first variation, J. Geom. Anal. 23 (2013), no. 2, 709–763. 10.1007/s12220-011-9261-5Search in Google Scholar

[29] L. Mugnai and M. Röger, The Allen–Cahn action functional in higher dimensions, Interfaces Free Bound. 10 (2008), no. 1, 45–78. 10.4171/IFB/179Search in Google Scholar

[30] W. W. Mullins, Two-dimensional motion of idealized grain boundaries, J. Appl. Phys. 27 (1956), 900–904. 10.1063/1.1722511Search in Google Scholar

[31] M. Röger and R. Schätzle, On a modified conjecture of De Giorgi, Math. Z. 254 (2006), no. 4, 675–714. 10.1007/s00209-006-0002-6Search in Google Scholar

[32] F. Schulze and B. White, A local regularity theorem for mean curvature flow with triple edges, J. Reine Angew. Math. 758 (2020), 281–305. 10.1515/crelle-2017-0044Search in Google Scholar

[33] L. Simon, Lectures on Geometric Measure Theory, Proc. Centre Math. Anal. Aust. National Univ. 3, Australian National University, Canberra, 1983. Search in Google Scholar

[34] S. Stuvard and Y. Tonegawa, Dynamical instability of minimal surfaces at flat singular points, preprint (2020), https://arxiv.org/abs/2008.13728. Search in Google Scholar

[35] S. Stuvard and Y. Tonegawa, An existence theorem for Brakke flow with fixed boundary conditions, Calc. Var. Partial Differential Equations 60 (2021), no. 1, Paper No. 43. 10.1007/s00526-020-01909-zSearch in Google Scholar

[36] K. Takasao and Y. Tonegawa, Existence and regularity of mean curvature flow with transport term in higher dimensions, Math. Ann. 364 (2016), no. 3–4, 857–935. 10.1007/s00208-015-1237-5Search in Google Scholar

[37] Y. Tonegawa, Integrality of varifolds in the singular limit of reaction-diffusion equations, Hiroshima Math. J. 33 (2003), no. 3, 323–341. 10.32917/hmj/1150997978Search in Google Scholar

[38] Y. Tonegawa, A second derivative Hölder estimate for weak mean curvature flow, Adv. Calc. Var. 7 (2014), no. 1, 91–138. 10.1515/acv-2013-0104Search in Google Scholar

[39] Y. Tonegawa, Brakke’s Mean Curvature Flow: An Introduction, Springer Briefs Math., Springer, Singapore, 2019. 10.1007/978-981-13-7075-5Search in Google Scholar

Received: 2021-11-13
Revised: 2022-02-03
Accepted: 2022-03-04
Published Online: 2022-07-23
Published in Print: 2024-01-01

© 2022 Walter de Gruyter GmbH, Berlin/Boston

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

Downloaded on 12.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/acv-2021-0093/html
Scroll to top button