Stochastic Surface Diffusion Flow

• Stochastic surface diffusion flow is a fourth-order geometric evolution of planar curves driven by curvature diffusion and random (Stratonovich) forcing.
• The analysis leverages arc-length parameterization and a fixed-domain SPDE reformulation to establish local well-posedness and investigate blow-up scenarios.
• The framework equates the evolution to a stochastic one-phase Stefan problem, coupling curvature dynamics with boundary motion to explain singularity formation.

Stochastic surface diffusion flow refers to the evolution of planar curves under a fourth-order geometric evolution law perturbed by stochastic forcing, typically modeled by a Stratonovich noise term. The canonical deterministic surface diffusion flow governs the evolution of a curve’s normal velocity as minus the Laplacian of the curvature; its stochastic generalization incorporates both geometric nonlinearities and randomness, and is analytically formulated as a stochastic one-phase Stefan problem in the curvature and the curve’s arc length. The principal object of interest is the coupled evolution of curvature $k(s, t)$ and length $L(t)$ for a family of immersed curves $\gamma_t$ on $\mathbb{R}^2$, parameterized by arc-length $s \in [0, L(t)]$.

1. Geometric Formulation and Stochastic Extension

The stochastic surface diffusion flow yields a dynamical system for the curvature-length pair $(k, L)$ with Stratonovich noise added to the classical surface diffusion drift. For an immersed planar curve $\gamma_t$ with curvature $k(s, t)$ and instantaneous length $L(t)$, the evolution reads: $$\begin{aligned} d k(s, t) &= \Big[-\partial_{ss} (\partial_{ss}k + \tfrac{1}{2}k^3) - k^2(\partial_{ss}k + \tfrac{1}{2}k^3)\Big] dt + (\partial_{ss} + k^2) \circ dW(s, t), \ dL(t) &= -\int_0^{L(t)} k(s, t)(\partial_{ss}k + \tfrac{1}{2}k^3)\, ds\, dt - \int_0^{L(t)} k(s, t) \circ dW(s, t) \end{aligned}$$ for $s \in [0, L(t)]$, with initial conditions $k(s, 0) = k_0(s)$ and $L(0) = L_0 > 0$. Here, $W(s, t)$ represents either a one-dimensional Brownian motion or an $\ell^2$-cylindrical Wiener process along the curve. For closed curves with $\int_0^{L(t)}k\, ds = 2\pi$, the Itô form simplifies to: $$\begin{aligned} dk &= \left[-\partial_{ss}(\partial_{ss}k + \tfrac{1}{2}k^3) - k^2(\partial_{ss}k + \tfrac{1}{2}k^3) + k^3\right] dt + k^2 dW_t, \ dL &= \int_0^{L(t)} k (\partial_{ss}k + \tfrac{1}{2}k^3)\, ds\, dt - 2\pi\,dW_t. \end{aligned}$$ This reflects the interaction between high-order geometric flows and stochastic forcing.

2. Arc-Length Parameterization and Frenet–Serret Framework

Any parameterization $F(u, t)$ of the curve yields a local speed $v(u, t) = |\partial_u F|$, and arc-length $s(u, t) = \int_0^u v(r, t)\, dr$, so that differentiation by arc-length satisfies $\partial_s = (1/v)\partial_u$. The geometric properties of the curve, specifically curvature $k$ and normal vector $N$, adhere to the standard Frenet–Serret relations.

Under normal-velocity flows $\partial_t F = V N$, this structure implies evolution equations: $$\partial_t k = \partial_{ss}V + k^2 V, \qquad \partial_t L = -\int_0^{L} k V ds$$ which dictate the interactions between curvature, arc-length, and their respective flows. Stochastic surface diffusion flow thus extends these deterministic relations by incorporating random perturbations along both the bulk and the boundary.

3. Stochastic One-Phase Stefan Problem Equivalence

The evolution system for $(k, L)$ is equivalent to a stochastic one-phase Stefan (free-boundary) problem:

• The “bulk” evolution for curvature $k$ on the domain $[0, L(t)]$:

$\partial_t k = \partial_{ss}V + k^2 V$

• The boundary motion law (Stefan condition):

$\partial_t L = -\int_0^{L} k V ds$

with $V$ embodying both deterministic drift ($-\partial_{ss}k - \tfrac{1}{2}k^3$) and stochastic forcing. Notably, both the PDE for curvature and the SDE for boundary carry stochastic terms. This construction contrasts with classical two-phase Stefan problems, where boundary motion is typically deterministic.

4. Quasilinear Parabolic Reformulation via Fixed-Domain SPDE–SDE

Transformation to a fixed domain employs the rescaling $s = r L(t)$, $r \in \mathbb{T} = \mathbb{R}/\mathbb{Z}$, and function $f(r, t) = k(rL(t), t)$. By Itô–Wentzell calculus, the problem becomes a coupled system in $E = X \times \mathbb{R}$ with $X = L^q(\mathbb{T})$: $$dX_t + A(X_t) X_t dt = F(t, X_t)\,dt + G(X_t)\, dW_t$$ where $X_t = (f(\cdot, t), L(t))^T$. The various operators are:

• $A(f, L)$: fourth-order elliptic operator

$$A(f, L) = \begin{pmatrix} L^{-4} \partial_{rrrr} + (\tfrac{5}{2}L^{-2}f^2 - 2\pi^2 r^2/L^2)\partial_{rr} & 0 \ 0 & 0 \end{pmatrix}$$

• $F(t, f, L)$: collects lower-order drift terms and nonlocal interactions, such as cubic $f$ nonlinearities and integral coupling.
• $G(f, L)$: multiplicative noise operator, e.g., in the one-dimensional case

$$G(f, L)(r) = \left(f(r)^2 - \frac{2\pi r}{L}\partial_r f(r), -2\pi\right)^T$$

The analysis proceeds in a UMD-Banach space scale: $X = L^q(\mathbb{T}) \times \mathbb{R}$, $X_1 = W^{4, q}(\mathbb{T}) \times \mathbb{R}$, with critical interpolation spaces $$X_p = (X, X_1)_{1 - 1/p, p} = B^{4-4/p}_{q, p}(\mathbb{T}) \times \mathbb{R}$$.

5. Local Well-posedness, Maximal Regularity, and Blow-up Criterion

Application of stochastic maximal regularity theory, particularly the Agresti–Veraar framework (2022), establishes:

• For all initial data $$(k_0, L_0) \in B^{4 - 4/p}_{q, p} \times (0, \infty)$$ with $1 > 4/p + 1/q$, there exists a unique local $L^p$-maximal strong solution $(f, L, \sigma)$, with $\sigma > 0$ almost surely, satisfying

$$(f, L) \in L^p(\Omega; H^{\theta, p}([0, \sigma_n]; H^{4(1 - \theta), q} \times \mathbb{R})) \cap L^p(\Omega; C([0, \sigma_n]; B^{4 - 4/p}_{q, p} \times \mathbb{R}))$$

for each $n \in \mathbb{N}$ and $\theta \in [0, 1/2)$, with stopping times

$$\sigma_n = \inf \left\{ t < \sigma : \|f(t)\|_{B^{4 - 4/p}_{q, p}} + |L(t)| > n \text{ or } L(t) \not\in (1/n, n) \right\}$$

• Blow-up alternative: if $\sigma < \infty$ almost surely then

$$\limsup_{t \to \sigma^-} \left( \|f(t)\|_{B^{4 - 4/p}_{q, p}} + L(t) + L(t)^{-1} \right) = \infty$$

This characterizes the maximal existence interval for solutions in terms of their geometric and analytical growth.

6. Geometric Blow-up Scenarios

Blow-up of the stochastic surface diffusion flow manifests through two distinct geometric mechanisms:

• Divergence of curvature: $\|k(\cdot, t)\|_{W^{4, q}} \to \infty$ at some point, corresponding to the development of a singularity such as a corner or cusp in the curve.
• Shrinking of the curve: $L(t) \to 0$ in finite time, corresponding to collapse of the curve to a point. These scenarios delimit the boundary of well-posed evolution and reflect critical transitions in the underlying geometry of the flow.

7. Energy Behavior and Regularity Considerations

Deterministically, Willmore (free-elastic) flow represents the $L^2$ gradient descent of the bending energy $E = \frac{1}{2}\int k^2 ds$. In the stochastic setting, the noise injects fluctuations in the normal direction; however, a net drift toward lower bending energy persists in the average sense. The absence of a closed stochastic energy-dissipation law is compensated by a priori maximal regularity estimates produced via the Agresti–Veraar machinery, which bound $\|f\|_{H^{4, q}}$ and $L(t)$ in $L^p$ spaces.

The stochastic Stefan formulation is novel in that both bulk motion (the PDE for curvature) and boundary motion (the SDE for arc length) incorporate Stratonovich forcing, with the boundary SDE coupled intrinsically to the evolution of curvature. This stands in contrast to classical two-phase stochastic Stefan problems, where the boundary law is deterministic. The general framework is robust and extends to fourth-order surface-diffusion (curve-diffusion) flows with minimal modification (Yan, 25 Nov 2025).