Free Boundary Willmore Flow

• Free boundary Willmore flow is a geometric evolution process where surfaces evolve to minimize bending energy while sliding along a prescribed support.
• It employs the gradient flow of the squared mean curvature, incorporating higher-order boundary conditions to enforce orthogonality and force balance.
• Analytical tools like Hanzawa transformations and maximal regularity estimates ensure short-time well-posedness and convergence to critical configurations.

The free boundary Willmore flow describes the evolution of immersed surfaces or curves under the gradient flow of the Willmore energy, where the surfaces are allowed to meet a prescribed support surface along their boundary, often with additional orthogonality or force-balancing conditions. The Willmore energy quantifies bending by integrating the squared mean curvature, and its geometric gradient flow governs the relaxation toward critical shapes under elastic constraints. When the boundary of the evolving surface—or of a planar curve in one dimension—is not fixed but is instead permitted to slide or evolve along a prescribed subset of the ambient space or a support manifold, the resulting evolution is termed the free boundary Willmore flow.

1. Geometric Framework and Variational Formulation

The Willmore energy for a compact oriented surface $\Sigma$ with (possibly nonempty) boundary is given by

$W(f) = \frac{1}{4} \int_\Sigma H^2 \, d\mu_f,$

where $f:\Sigma \rightarrow \mathbb{R}^3$ is a $W^{4,2}$ immersion, $H$ is the mean curvature, and $d\mu_f$ is the induced area measure (Dall'Acqua et al., 26 Jan 2026). A free boundary constraint involves both positional and geometric conditions imposed on the boundary $\partial\Sigma$:

• The image $f(\partial\Sigma) \subset S$, where $S \subset \mathbb{R}^3$ is a fixed real-analytic support surface.
• The normal vector $\nu_f$ of the immersion at the boundary is orthogonal to the normal of the support: $\langle \nu_f, N^S \circ f \rangle = 0$ along $\partial\Sigma$.

These are typically accompanied by boundary conditions on higher derivatives, derived from variational principles to ensure well-posedness of the evolution and absence of spurious boundary terms in the first variation.

A standard example involves surfaces that are locally close to a half-sphere of small radius (parameterized as $\phi: S^2_+ \rightarrow \mathbb{R}^3$), touching and moving along $\partial\Omega$, the boundary of a smooth domain $\Omega$ (Metsch, 2022). Immersions are required to meet $\partial\Omega$ orthogonally and may be constrained to preserve area or other global quantities.

2. Evolution Equations and Boundary Conditions

The Willmore flow is the $L^2$–gradient flow of the Willmore functional. For a surface without boundary, this takes the form

$$\partial_t f = -\frac{1}{2} \left(\Delta_g H + |A^\circ|^2 H\right) \nu_f,$$

where $\Delta_g$ is the Laplace–Beltrami operator, $|A^\circ|^2$ is the squared norm of the traceless second fundamental form, and $\nu_f$ is the unit normal (Dall'Acqua et al., 26 Jan 2026). For surfaces with boundary meeting a support surface $S$, the evolution is modified to enforce the free boundary constraints.

Boundary conditions for the free boundary Willmore flow include:

• Positional: $f(t,\partial \Sigma) \subset S$ at all $t$.
• Geometric: $\langle \nu_f, N^S \circ f \rangle = 0$ on $\partial\Sigma$.
• Third-order natural boundary: $\partial_\eta H + A^S(\nu_f, \nu_f) H = 0$ along $\partial\Sigma$, where $\eta$ is the outward co-normal and $A^S$ the second fundamental form of $S$ (Metsch, 2022, Dall'Acqua et al., 26 Jan 2026).

When additional constraints such as area preservation are imposed, a Lagrange multiplier appears in the equation, projecting the Willmore operator onto the $L^2$–orthogonal complement of the mean curvature (Metsch, 2022):

$$\langle \partial_t \phi, \nu \rangle = -P_H^\perp[W(\phi)],$$

where $P_H^\perp$ is the $L^2$–orthogonal projection.

For the 1D (curve) case, the elastic (Willmore-plus-length) energy is

$E(\gamma)=\int_\gamma (k^2 + \mu) \, ds,$

with free (Navier-type) boundary conditions ensuring attachment at the prescribed set (e.g., $\gamma_2(y)=0$ for planar curves) and orthogonality of the normal at the endpoints (Diana, 2023).

3. Well-Posedness, Regularity, and Analytical Methods

Short-time existence and uniqueness for free boundary Willmore flows, including those with line tension or area preservation, has been established in various geometric settings (Metsch, 2022, Abels et al., 2014, Diana, 2023). The analytical machinery employs:

• Hanzawa transformations: Representing evolving surfaces as graphs over a fixed reference, which stabilizes the moving boundary and facilitates analysis via fixed domain methods (Abels et al., 2014).
• Maximal regularity estimates in parabolic Hölder or $L^p$ scales, often utilizing Solonnikov's theory for higher-order systems (Diana, 2023).
• Banach manifold and Fredholm theory: Identification of constraint submanifolds of immersions, construction of analytic local slices, and Lyapunov–Schmidt reduction to handle finite-dimensional degeneracies (Dall'Acqua et al., 26 Jan 2026, Metsch, 2022).

For flows in the class of half-spheres sliding on $\partial\Omega$, coupled evolution equations for the shape graph function and barycenter result in a quasilinear parabolic PDE with boundary constraints, with existence shown via implicit function theorems and patching in time (Metsch, 2022).

In the stochastic setting, the Willmore flow of planar curves can be cast as a stochastic Stefan-type free boundary problem, with well-posedness proven for a quasilinear SPDE involving both curvature and time-dependent domain length, and a precise blow-up alternative characterizing singularity formation or domain collapse (Yan, 25 Nov 2025).

4. Stability, Convergence, and Asymptotic Behavior

Stability of free boundary Willmore immersions near critical points is quantified using Łojasiewicz–Simon gradient inequalities. If a critical immersion $f_*$ is a local minimizer, solutions starting sufficiently close to $f_*$ converge to a (possibly reparameterized) free boundary Willmore immersion, with algebraic energy decay rates controlled by the spectrum of the linearized operator (Dall'Acqua et al., 26 Jan 2026).

The convergence of the free boundary flow in the class of small half-spheres is governed by an effective ODE for the barycenter, which, up to small errors, follows a gradient flow on the support surface driven by its mean curvature:

$\frac{d}{d\tau} \xi = \frac{3}{2} \nabla H^S(\xi)$

in the limit of small sphere radius, where $\xi$ parameterizes the barycenter along $S$ (Metsch, 2022). Under nondegeneracy of the critical points of $H^S$, every flow converges to a unique critical configuration determined by these points.

Quantitative rigidity results show that in neighborhoods of free boundary minimal immersions (with $H=0$), any nearby free boundary Willmore immersion must also be minimal, establishing local uniqueness and rigidity (Dall'Acqua et al., 26 Jan 2026).

5. Generalizations, Stochastic and Planar Flows

Extensions of the free boundary Willmore flow include stochastic variants and lower-dimensional analogs:

• Stochastic free boundary Willmore curves are formulated as quasilinear SPDEs for curvature and interval length, allowing random perturbations in force balance and domain evolution. Local strong solutions exist up to a maximal stopping time, with blow-up characterized precisely by curvature singularities or collapse of the interval length (Yan, 25 Nov 2025).
• Planar elastic flows with partial free boundary (Navier conditions) are globally well-posed for nondegenerate initial data, with energy methods ruling out blow-up except for explicit collapse or boundary degeneracy. The analytic setting exploits the structure of boundary terms in the first variation, ensuring natural force balancing at the moving endpoints (Diana, 2023).

6. Applications and Broader Context

The free boundary Willmore flow has significance in modeling the behavior of biological membranes and elastic objects interacting with substrates, where surface flexibility is modulated by bending energy but interactions at the interface are governed by geometric compatibility, force transmission, or energetic wetting and line tension terms (Abels et al., 2014). The area-preserving and line tension modifications accommodate physical constraints such as incompressibility and energetic penalties for contact line motion.

Quantitative stability, convergence, and rigidity results for these flows inform understanding of equilibrium shapes, metastable states, and response to perturbations in both deterministic and stochastic regimes (Dall'Acqua et al., 26 Jan 2026, Metsch, 2022, Yan, 25 Nov 2025). The free boundary framework unifies geometric PDE evolution with variational principles under natural geometric and physical constraints, providing a foundation for rigorous analysis of elastic interface dynamics in various settings.