Flat Flow in Geometric Evolutions

• Flat flow is a heterogeneous concept describing variationally constructed weak solutions in curvature-driven evolutions, such as volume preserving mean curvature flow and surface diffusion via discrete minimizing movements.
• It encompasses asymptotic flattening phenomena in noncompact settings, where flows like Yamabe and Ricci converge to flat metrics and preserve asymptotic flatness.
• In discrete differential geometry and continuum mechanics, flat flow signifies piecewise flat approximations on simplicial manifolds and flat profiles in applications like flat diffusers, granular flows, and flat-top solitons.

Flat flow is a heterogeneous term in current mathematical and physical literature. In geometric analysis, it denotes a weak solution constructed by a discrete minimizing movements scheme for curvature-driven evolutions such as the volume preserving mean curvature flow and the surface diffusion equation (Julin et al., 2022, Cicalese et al., 19 Feb 2025). In geometric evolution on noncompact spaces, closely related usage describes flows that preserve asymptotic flatness or converge to flat metrics or flat slices, as in Yamabe flow, Ricci flow, and mean curvature flow on asymptotically flat backgrounds (Chen et al., 2021, Chen, 2019, Kroencke et al., 2019). In discrete differential geometry, piecewise flat Ricci flow evolves edge lengths on simplicial manifolds and uses deficit angles as curvature proxies (Conboye et al., 2016, Conboye, 2018). The adjective also appears in continuum mechanics and nonlinear waves, where it modifies the ambient geometry or profile, as in flat diffusers, flat plates, flat frictional channels, and flat-top solitons (Fedyushkin et al., 2023, Sin et al., 2013, Brodu et al., 2011, Alotaibi et al., 2023). This suggests that “flat flow” is not a single standardized object but a family of usages organized around weak geometric evolutions, asymptotic flattening, piecewise flat discretization, and transport in flat geometries.

1. Flat flow as a variational weak solution

A precise modern meaning of flat flow arises in the volume preserving mean curvature flow (VPMCF). There the evolution law is

$V_t = - (H_{E_t} - \overline{H}_{E_t}),$

where $V_t$ is the normal velocity of $\partial E_t$, $H_{E_t}$ is the mean curvature, and $\overline{H}_{E_t}$ is the average mean curvature enforcing the volume constraint (Julin et al., 2022). Because the equation is nonlocal, does not satisfy a comparison principle, and may develop singularities even in dimension $2$, classical solutions may fail to exist globally in time. The flat flow is therefore defined as a weak solution obtained from a discrete-time minimizing movements scheme.

Given a time step $h>0$, the discrete VPMCF scheme constructs sets $E^h_k$ by minimizing

$$F_h(F, E^h_k) := P(F) + \frac{1}{2h} \int_{F \Delta E^h_k} d_{E^h_k}(x)\,dx$$

among measurable sets $F$ satisfying $V_t$0, where $V_t$1 is the perimeter, $V_t$2 is the signed distance to $V_t$3, and $V_t$4 is the symmetric difference (Julin et al., 2022). Any limit, as $V_t$5, of time-interpolated discrete flows is called a flat flow. In this framework, the term “flat” does not describe zero curvature; it designates a variationally constructed generalized evolution.

An analogous usage appears for the surface diffusion equation

$V_t$6

a fourth-order geometric PDE that preserves the volume of each component while decreasing perimeter (Cicalese et al., 19 Feb 2025). Here flat flow is built through the discrete minimizing movements scheme proposed by Cahn and Taylor. The discrete step minimizes perimeter plus an $V_t$7-type penalization, with an additional tubular-neighborhood constraint used technically to guarantee compactness and existence of minimizers (Cicalese et al., 19 Feb 2025). In both VPMCF and surface diffusion, the terminology identifies a weak evolution generated by an implicit variational scheme rather than by direct classical PDE theory.

2. Consistency, regularity, and uniqueness before singularities

For VPMCF, the central result is a consistency principle. If the initial set $V_t$8 is open, bounded, and satisfies a uniform ball condition, equivalently $V_t$9 regularity, then there exists $\partial E_t$0 such that any flat flow starting from $\partial E_t$1 remains uniformly $\partial E_t$2, satisfying a uniform ball condition with radius $\partial E_t$3 for all $\partial E_t$4 (Julin et al., 2022). The flow becomes instantaneously smooth, with estimates

$\partial E_t$5

and any flat flow agrees with the classical solution as long as the latter exists (Julin et al., 2022). A direct consequence is uniqueness and smoothness up to the first singular time.

The proof strategy is geometric and discrete. A two-point function method propagates the uniform ball condition at the discrete level, while higher-order regularity is derived from discrete analogues of classical evolution identities (Julin et al., 2022). The same minimizing-movements approach also applies to standard mean curvature flow and yields an alternative proof of consistency that does not rely on the comparison principle (Julin et al., 2022).

For surface diffusion in dimension three, the Cahn–Taylor flat flow is likewise consistent with smooth evolution. If the initial set is sufficiently regular, the discrete scheme converges to the unique smooth solution of the equation, and the convergence holds up to the maximal smooth existence time (Cicalese et al., 19 Feb 2025). The approximating discrete flows are shown to remain uniformly regular, first in $\partial E_t$6 and then in $\partial E_t$7, while the associated height functions converge in $\partial E_t$8 to the classical solution (Cicalese et al., 19 Feb 2025). The significance of these results lies in the fact that both VPMCF and surface diffusion lack the direct compactness and comparison mechanisms available for lower-order scalar parabolic equations.

3. Flatness as an asymptotic state under geometric evolution

A second major usage concerns flows on asymptotically flat manifolds or flows converging to flat limiting geometries. For the Yamabe flow

$\partial E_t$9

starting from an asymptotically flat manifold, long-time existence holds and asymptotic flatness is preserved (Chen et al., 2021). The decisive criterion for convergence is the Yamabe invariant

$H_{E_t}$0

If $H_{E_t}$1, the flow converges in weighted Hölder norms to the unique asymptotically flat, scalar-flat metric $H_{E_t}$2 in the conformal class of $H_{E_t}$3; if $H_{E_t}$4, the flow does not converge (Chen et al., 2021). When the scalar curvature is nonnegative and integrable, the ADM mass at time infinity drops by the limit of the total scalar curvature along the flow (Chen et al., 2021). Related results establish that the ADM mass is well-defined and monotone non-increasing under Yamabe flow on asymptotically flat manifolds, and that in dimensions $H_{E_t}$5 or $H_{E_t}$6 it is invariant under the flow (Cheng et al., 2011). A separate global existence theorem gives unique global Yamabe flow on asymptotically flat manifolds of order $H_{E_t}$7, with preservation of the asymptotically flat structure and, for $H_{E_t}$8, preservation of ADM mass under suitable hypotheses (Ma, 2021).

Ricci flow provides an analogous flattening mechanism. On asymptotically flat manifolds of dimension $H_{E_t}$9, if the scale-invariant $\overline{H}_{E_t}$0-norm of curvature is sufficiently pinched relative to the inverse Sobolev constant, then the Ricci flow exists for all time and converges in weighted Hölder spaces to the flat Euclidean metric on $\overline{H}_{E_t}$1 (Chen, 2019). In particular, the initial manifold must have been diffeomorphic to $\overline{H}_{E_t}$2 (Chen, 2019). On compact manifolds, Ricci flow also yields a strengthened version of the Gromov–Ruh theorem: the pointwise condition $\overline{H}_{E_t}$3 can be replaced by the weaker integral condition

$\overline{H}_{E_t}$4

and manifolds satisfying it are diffeomorphic to infranil manifolds (Chen et al., 2022). In this sense, Ricci flow turns integral near-flatness into classical almost-flat pinching.

In Lorentzian product geometry, mean curvature flow of uniformly spacelike graphs in $\overline{H}_{E_t}$5, with $\overline{H}_{E_t}$6 asymptotically flat and $\overline{H}_{E_t}$7, exists for all times and converges uniformly in all derivatives to the flat slice $\overline{H}_{E_t}$8 when the initial graph is asymptotic to that slice at infinity (Kroencke et al., 2019). The flattening mechanism is enforced through barriers at infinity, gradient estimates, and higher-derivative control (Kroencke et al., 2019). A further stability theorem states that an ALE Ricci-flat manifold that is linearly stable and integrable is dynamically stable under Ricci flow, so any nearby Ricci flow exists for all time and converges modulo diffeomorphism to a nearby ALE Ricci-flat metric; ALE Calabi–Yau manifolds satisfy these hypotheses (Deruelle et al., 2017).

4. Piecewise flat Ricci flow and discrete curvature

A third meaning of flat flow is explicitly discrete. Piecewise flat Ricci flow approximates Ricci flow on triangulated $\overline{H}_{E_t}$9-manifolds by evolving the lengths of edges in a simplicial complex (Conboye et al., 2016). Curvature is concentrated on hinges, which in dimension three are edges, through the deficit angle

$2$0

where the sum is over incident dihedral angles (Conboye et al., 2016). The discrete scalar curvature at a vertex $2$1, sectional curvature orthogonal to an edge $2$2, and Ricci curvature along $2$3 are assembled from these deficit angles and dual volumes, yielding the edge evolution law

$2$4

or its normalized form (Conboye et al., 2016, Conboye, 2018).

This framework was tested on diverse manifolds, including the $2$5-sphere, $2$6-cylinder, Nil geometry, Gowdy manifolds, a three-torus embedded in Euclidean four-space, and a perturbation of a flat three-torus (Conboye et al., 2016, Conboye, 2018). A flat metric is characterized by the vanishing of all deficit angles, and in simulations the average magnitudes of Ricci curvature and deficit angles decay to zero when the limit is flat (Conboye, 2018). For Nil and Gowdy manifolds, the piecewise flat flow converges to known smooth Ricci flow solutions; for the embedded and perturbed tori, the flow is asymptotically flat, with minimum and maximum coordinate-curve lengths approaching one another under the evolution (Conboye, 2018).

Stability is a separate issue. The original piecewise flat Ricci flow exhibits exponential numerical instability on conventional cubic and skew triangulations, even for flat manifolds that should remain stationary under the flow (Conboye, 2023). Linearization gives

$2$7

and the coefficient matrix has a positive eigenvalue for these triangulations; in the cubic case the instability rate is $2$8 for blocks of volume $2$9 (Conboye, 2023). The remedy is to enforce “flat blocks,” meaning that body-diagonal lengths are constrained so that the associated deficit angles vanish at each step (Conboye, 2023). With this adaptation, the largest eigenvalue becomes zero and the spurious exponential mode disappears, which is essential for reliable convergence to smooth Ricci flow behavior (Conboye, 2023).

5. Flat sides and interface regularity in curvature flows

Another established usage concerns curvature flows with a flat side. For the scalar curvature flow

$h>0$0

the flat side is separated from the strictly convex region by an interface $h>0$1 (Jang et al., 2018). Under suitable initial non-degeneracy conditions,

$h>0$2

with $h>0$3, the interface propagates with finite and non-degenerate speed until the flat side vanishes (Jang et al., 2018). The estimates

$h>0$4

hold uniformly near the interface, and the analysis yields optimal decay estimates of curvatures, an Aronson–Bénilan-type lower bound

$h>0$5

and Hölder regularity of curvature ratios up to the free boundary (Jang et al., 2018). The solution is smooth up to the interface for all time until the flat side disappears (Jang et al., 2018).

For the $h>0$6-Gauss curvature flow with flat side,

$h>0$7

the flat side persists for some time when $h>0$8, whereas for $h>0$9 immediate strict convexity occurs (Huang et al., 2024). The main regularity theorem addresses the hypersurface near the interface rather than only the interface itself. Writing the evolving hypersurface as a graph $E^h_k$0 and defining

$E^h_k$1

the authors prove that if $E^h_k$2, then $E^h_k$3 is $E^h_k$4 up to the interface for $E^h_k$5 (Huang et al., 2024). If $E^h_k$6, then $E^h_k$7 is smooth up to the interface; otherwise the $E^h_k$8 regularity is optimal (Huang et al., 2024). The degeneracy is controlled through a hodograph transformation and weighted Hölder/Schauder theory adapted to the singular behavior near the flat part (Huang et al., 2024).

6. Flat geometries and flat profiles in fluid mechanics, granular media, and nonlinear waves

Outside geometric analysis, the adjective “flat” frequently specifies the ambient geometry or the waveform rather than the solution concept. In a flat diffuser, numerical simulations of the incompressible Navier–Stokes equations show that weak periodic vibration of the inlet velocity can symmetrize an asymmetric laminar viscous flow (Fedyushkin et al., 2023). The domain is a plane diffuser bounded by two arcs with opening angle $E^h_k$9 and length $$F_h(F, E^h_k) := P(F) + \frac{1}{2h} \int_{F \Delta E^h_k} d_{E^h_k}(x)\,dx$$0, with inlet condition

$$F_h(F, E^h_k) := P(F) + \frac{1}{2h} \int_{F \Delta E^h_k} d_{E^h_k}(x)\,dx$$1

constant outlet pressure, and no-slip walls (Fedyushkin et al., 2023). Even amplitudes less than $$F_h(F, E^h_k) := P(F) + \frac{1}{2h} \int_{F \Delta E^h_k} d_{E^h_k}(x)\,dx$$2 of the inlet velocity restore symmetry, independently of vibration frequency within the investigated range, and the simulations reproduce Richardson’s “annular effect,” namely near-wall peaks in the time-averaged longitudinal velocity under harmonic forcing (Fedyushkin et al., 2023).

For unsteady reversed stagnation-point flow over a flat plate, similarity reduction of the incompressible Navier–Stokes equations produces a third-order nonlinear equation

$$F_h(F, E^h_k) := P(F) + \frac{1}{2h} \int_{F \Delta E^h_k} d_{E^h_k}(x)\,dx$$3

with boundary conditions $$F_h(F, E^h_k) := P(F) + \frac{1}{2h} \int_{F \Delta E^h_k} d_{E^h_k}(x)\,dx$$4, $$F_h(F, E^h_k) := P(F) + \frac{1}{2h} \int_{F \Delta E^h_k} d_{E^h_k}(x)\,dx$$5, and $$F_h(F, E^h_k) := P(F) + \frac{1}{2h} \int_{F \Delta E^h_k} d_{E^h_k}(x)\,dx$$6 (Sin et al., 2013). The study shows that no steady solution exists with these conditions, emphasizing the fundamentally unsteady character of the reversed stagnation-point problem (Sin et al., 2013). Viscous terms are negligible only in the outer region; they remain essential near the plate for satisfying the no-slip condition and for describing separation (Sin et al., 2013).

In shallow granular flows down flat frictional channels, the flat base generates a basal rolling layer of slipping grains that acts as an effective bumpy base for the overlying bulk (Brodu et al., 2011). Above that rolling layer, the flow obeys Bagnold-type scaling, while at larger inclination angles the system transitions from unidirectional layered flow to a convective regime with steady counter-rotating longitudinal vortices and an inverted density profile (Brodu et al., 2011). The rheological transition appears as a discontinuity in the $$F_h(F, E^h_k) := P(F) + \frac{1}{2h} \int_{F \Delta E^h_k} d_{E^h_k}(x)\,dx$$7 constitutive relation, and the convective regime is interpreted as a Rayleigh–Bénard type instability driven by the hotter granular basal layer (Brodu et al., 2011).

In nonlinear wave dynamics, flat-top solitons governed by a nonlinear Schrödinger equation with cubic–quintic nonlinearity display unidirectional flow through two reflectionless potential wells of slightly different depths (Alotaibi et al., 2023). The phenomenon is restricted to a finite velocity window, which is wider for shallower wells and narrower for wider flat-top solitons (Alotaibi et al., 2023). Transport is quantified by reflection, transmission, and trapping coefficients satisfying $$F_h(F, E^h_k) := P(F) + \frac{1}{2h} \int_{F \Delta E^h_k} d_{E^h_k}(x)\,dx$$8, and the results show that the soliton’s width controls how strongly it resolves the asymmetry of the double-well potential (Alotaibi et al., 2023).

Across these settings, “flat” may refer to a weak variational evolution, an asymptotically flat ambient space, a piecewise flat discretization, a flat side in a free-boundary problem, or a flat channel or flat-top profile in continuum mechanics and nonlinear optics. The common feature is not a single formal definition but the repeated appearance of flattening, flat asymptotics, or flat geometric structure as an organizing principle of the flow.