---
title: Area-Preserving Curve-Shortening Flow
url: https://www.emergentmind.com/topics/area-preserving-curve-shortening-flow-apcsf
type: topic
---

# Area-Preserving Curve-Shortening Flow

Searching arXiv for recent and foundational work on APCSF and closely related geometric-flow results.
The **area-preserving curve-shortening flow** (APCSF) is a nonlocal geometric evolution of planar curves in which the normal velocity is the curvature minus its spatial average, so that the enclosed area is preserved while the length decreases. In its classical closed-curve form, introduced by Gage, the flow is
\[
\partial_t \gamma = (\kappa-\bar\kappa)\,\nu,\qquad 
\bar\kappa(t)=\frac{1}{L(t)}\int_{\gamma_t}\kappa\,ds,
\]
where \(\kappa\) is the signed curvature, \(\nu\) the unit normal, and \(L(t)\) the length of the evolving curve. For simple closed planar curves, \(\int_{\gamma_t}\kappa\,ds=2\pi\), hence \(\bar\kappa=2\pi/L(t)\) [2105.11930]. APCSF is the \(L^2\)-gradient flow of length under an area constraint, and in this sense is the one-dimensional analogue of volume-preserving mean-curvature flow. Its theory spans closed and open curves, free-boundary problems with Neumann conditions, singularity formation and blow-up analysis, self-similar and homothetic solutions, and a growing numerical literature on structure-preserving discretizations [1503.05814], [1701.08948], [2508.00352], [2208.01324], [2605.11721].

## 1. Definition and basic variational structure

For a smooth family of simple closed plane curves \(\gamma:S^1\times[0,T)\to\mathbb{R}^2\), the classical curve-shortening flow is
\[
\partial_t\gamma=\kappa\,\nu,
\]
whereas APCSF inserts the nonlocal correction \(-\bar\kappa\,\nu\) to preserve enclosed area [1701.08948]. In the standard planar orientation convention, \(\bar\kappa\) is the average curvature,
\[
\bar\kappa(t)=\frac{1}{L(t)}\int_{\gamma_t}\kappa\,ds.
\]
For simple closed curves with turning number \(1\), Gauss–Bonnet yields \(\int_{\gamma_t}\kappa\,ds=2\pi\), and therefore
\[
\partial_t X = \Bigl(\kappa-\frac{2\pi}{L(t)}\Bigr)N
\]
in the notation used for Gage’s flow [2105.11930], [2208.01324].

The defining structural identities are area preservation and length dissipation. Under
\[
\partial_t X=\beta N,\qquad \beta=\kappa-\frac{2\pi}{L},
\]
the enclosed area satisfies
\[
\frac{d}{dt}A(t)=0,
\]
while the length satisfies
\[
\frac{d}{dt}L(t)=-\int_\gamma \Bigl(\kappa-\frac{2\pi}{L}\Bigr)^2 ds\le 0
\]
[2105.11930], [2208.01324]. This identifies APCSF as a constrained steepest descent for length at fixed area. A plausible implication is that APCSF should be viewed less as a perturbation of curve-shortening flow than as an isoperimetric relaxation dynamics: it dissipates perimeter while remaining on a fixed-area manifold of immersions.

For closed curves, the curvature evolves by the nonlocal semilinear parabolic equation
\[
\kappa_t=\kappa_{ss}+\kappa^3-\frac{2\pi}{L}\kappa^2
\]
[2105.11930]. The extra quadratic term distinguishes APCSF from classical curve-shortening flow and is the source of both the global area constraint and several analytical complications in nonconvex regimes.

## 2. Closed-curve theory: convexity, star-shaped data, and convergence

The convex closed-curve theory is the classical core of APCSF. For convex initial data, Gage established that the flow preserves enclosed area, decreases length, exists globally, and converges after normalization to a round circle; this standard picture is recalled in later works [2105.11930], [2604.14804]. The circle is the unique stationary solution among simple closed curves with fixed area, since stationarity forces \(\kappa\equiv\bar\kappa\), hence constant curvature.

A major development concerns nonconvex initial curves. The paper “Star-shaped Centrosymmetric Curves under Gage’s Area-preserving Flow” proves that smooth embedded star-shaped curves that are also centrosymmetric evolve smoothly for all time, become convex in finite time, and converge to a circle as \(t\to\infty\) [2105.11930]. The argument combines a polar-graph formulation
\[
X(\theta,t)=r(\theta,t)(\cos\theta,\sin\theta)
\]
with a quasilinear parabolic PDE for the radial function,
\[
r_t = \frac{1}{g^2} r_{\theta\theta} -\frac{2}{r g^2}\,r_\theta^2 -\frac{r}{g^2} +\frac{2\pi g}{r L},
\qquad g=\sqrt{r^2+r_\theta^2},
\]
and uses preservation of star-shapedness and centrosymmetry to obtain curvature bounds and global existence [2105.11930].

That centrosymmetry assumption was later removed. The paper “Star-shaped Curves under Gage's Area-preserving Flow and the CSF” proves the folklore conjecture that smooth embedded star-shaped initial curves evolve globally under Gage’s area-preserving flow and converge to a circle [2412.18102]. Its key mechanism is Dittberner’s singularity analysis, mediated through the Chow–Liou–Tsai turning-angle invariant \(\Theta(X_0)\). The paper shows that star-shaped embedded curves satisfy \(\Theta(X_0)>-\pi\), which is sufficient to exclude finite-time singularities under APCSF [2412.18102]. This result substantially enlarges the known nonconvex global-existence class.

At the same time, star-shapedness itself is not preserved in full generality. The same work constructs a star-shaped “flying wing” example showing that Gage’s area-preserving flow may lose star-shapedness during the evolution [2412.18102]. This also yields a negative answer to Mantegazza’s open problem on whether curve-shortening flow always preserves star-shapedness [2412.18102]. Thus, for APCSF, global regularity and preservation of a given geometric class can diverge sharply.

## 3. Free-boundary APCSF with Neumann conditions

A substantial branch of the theory concerns **open curves** whose endpoints lie on a support curve \(\Sigma\) and satisfy a Neumann free-boundary condition. In this setting, the evolution is
\[
\partial_t\gamma(p,t)=(\kappa(p,t)-\bar\kappa(t))\,\nu(p,t),
\qquad p\in(a,b),
\]
with endpoints constrained to \(\Sigma\) and tangent orthogonality at the boundary [1503.05814], [1701.08948]. For a convex closed support curve \(\Sigma\), the boundary condition is
\[
\tau(a,t)=-\nu_\Sigma(\gamma(a,t)),\qquad \tau(b,t)=\nu_\Sigma(\gamma(b,t)),
\]
and for the line support \(\Sigma=\{y=0\}\), it becomes
\[
\tau(a,t)=e_2,\qquad \tau(b,t)=-e_2
\]
[1701.08948].

Since the evolving curve is not closed, “area preservation” is defined relative to a boundary arc of \(\Sigma\). One attaches the short piece \(\sigma_t\subset\Sigma\) joining the endpoints and defines the assembled closed curve \(\gamma_t+\sigma_t\). Its oriented area is
\[
A(\gamma_t+\sigma_t)
=\frac12\int_{\gamma_t}(p^1\,dp^2-p^2\,dp^1)
+\frac12\int_{\sigma_t}(p^1\,dp^2-p^2\,dp^1),
\]
and APCSF preserves this quantity while still decreasing the length of \(\gamma_t\) [1503.05814], [1701.08948].

The 2015 paper on APCSF with Neumann free boundary conditions establishes a global-existence and subconvergence theory under geometric smallness assumptions on the initial curve. For strictly convex, embedded initial curves lying outside a convex support domain and satisfying a suitable length restriction, the flow exists for all time and subconverges smoothly to an arc of a circle meeting \(\Sigma\) orthogonally and enclosing the same area with \(\Sigma\) as the initial curve [1503.05814]. The limiting stationary states are exactly circular arcs with orthogonal contact to the support.

This free-boundary theory already displays a dichotomy that later becomes central: either the flow exists globally and approaches a circular equilibrium, or it develops a finite-time singularity. The 2017 singularity paper sharpened this dichotomy considerably [1701.08948].

## 4. Singularity formation, type II blow-up, and grim-reaper limits

The paper “Singularities of the area preserving curve shortening flow with a Neumann free boundary condition” gives explicit geometric criteria guaranteeing finite-time singularity formation for APCSF outside a convex support curve or at a line [1701.08948]. For initial curves satisfying \(L_0<d_\Sigma\), with \(d_\Sigma\) the minimum width of the support curve, and with total curvature in a specified interval
\[
(2l-2)\pi < \int_{\gamma_0}\kappa\,ds < 2l\pi,
\]
the flow must develop a finite-time singularity if either
\[
A(\gamma_0+\sigma_0)<0,
\]
or
\[
A(\gamma_0+\sigma_0)>0,\qquad
\frac{L_0^2}{A(\gamma_0+\sigma_0)} \le \pi\,\frac{(2l-1)^2}{l}
\]
[1701.08948].

A central result is the classification of any such finite-time singularity as **type II**, not type I. In the paper’s formulation, if \(T_{\max}<\infty\), then
\[
\max_{p\in[a,b]}|\kappa|(p,t)\to\infty\qquad (t\to T_{\max}),
\]
and moreover
\[
\max_{p\in[a,b]} \bigl(|\kappa|^2(p,t)(T_{\max}-t)\bigr)
\]
is unbounded as \(t\to T_{\max}\) [1701.08948]. This is a notable contrast with classical curve-shortening flow, where finite-time singularities of embedded convex curves are type I shrinkers. Here the area-preserving constraint shifts the blow-up mechanism toward translators.

The blow-up analysis uses Hamilton’s rescaling. Choosing spacetime points \((p_j,t_j)\) where
\[
Q_j=|\kappa(p_j,t_j)|\to\infty,
\]
one rescales by
\[
\tilde\gamma_j(p,\tau)=Q_j\Bigl(\gamma\Bigl(p,\frac{\tau}{Q_j^2}+t_j\Bigr)-\gamma(p_j,t_j)\Bigr).
\]
Because the average-curvature term scales as \(Q_j^{-1}\bar\kappa\), it disappears in the limit, and the blow-up limit solves the pure curve-shortening flow
\[
\partial_\tau\tilde\gamma_\infty=\tilde\kappa_\infty\,\tilde\nu_\infty
\]
[1701.08948]. For convex initial data satisfying the singularity criterion, the limit is either a full grim reaper or, at a boundary singularity, a half grim reaper meeting a line orthogonally [1701.08948].

The line-supported case is handled by reflection. Reflecting the free-boundary curve across the line \(y=0\) produces a closed curve \(\delta_t\), and Escher–Ito’s theory for closed APCSF can then be invoked. In this reflected setting, the turning index is always an odd integer, and finite-time singularity follows if either the reflected enclosed area is negative or an isoperimetric condition involving the index holds [1701.08948]. This makes the free-boundary APCSF theory structurally parallel to the closed immersed-curve theory.

## 5. Self-similar solutions and shrinkers of APCSF

Although APCSF is area-preserving for embedded simple closed curves of turning number \(1\), homothetic solutions can still exist in the broader immersed setting because the enclosed signed area may vanish. The paper “Shrinkers of the area-preserving curve-shortening flow: Existence and saddle-point property” studies such homothetic evolutions and identifies their governing equation [2508.00352].

If
\[
X(\cdot,t)=p_*+\psi(t)\bigl(X(\cdot,0)-p_*\bigr)
\]
is a homothetic APCSF evolution, then after normalization the initial curve satisfies
\[
\kappa=\langle x,\mathsf{n}\rangle+\bar\kappa[x]
\]
[2508.00352]. This differs from the classical self-shrinker equation for curve-shortening flow by the nonlocal additive term \(\bar\kappa[x]\). A direct consequence is that APCSF shrinkers must have zero signed enclosed area:
\[
\vol[x]=0
\]
[2508.00352].

The paper connects APCSF shrinkers to the theory of **\(\lambda\)-curves**, namely solutions of
\[
\kappa=\langle x,\mathsf{n}\rangle+\lambda.
\]
In tangential polar coordinates, the curvature \(\kappa(\theta)\) of a \(\lambda\)-curve satisfies the ODE
\[
\kappa''=\frac{1}{\kappa}-\kappa+\lambda,
\]
with first integral
\[
E=\frac12(\kappa')^2+\frac{\kappa^2}{2}-\lambda\kappa-\log\kappa
\]
[2508.00352]. Closedness is controlled by a semi-period map \(\Theta_\lambda(E)\); rational values of \(\Theta_\lambda(E)/\pi\) correspond to closed \(n\)-symmetric curves with tangent turning index \(m\) [2508.00352].

Using this framework, the paper proves existence of non-circular APCSF shrinkers for every coprime pair \((m,n)\) satisfying
\[
\frac12<\frac{m}{n}<1.
\]
For each such rational \(m/n\), there exists a non-circular \(n\)-symmetric APCSF shrinker with tangent turning index \(m\) [2508.00352]. This is an APCSF analogue of the Abresch–Langer theory for classical curve-shortening flow, although the classification is presently partial rather than complete.

The same paper also establishes a **saddle-point property**. Small normal perturbations of a non-circular APCSF shrinker bifurcate into two distinct dynamical regimes: one side evolves toward an \(m\)-fold cover of a circle, while the other develops a singular curve with \(n\) cusps in finite time [2508.00352]. This suggests that APCSF shrinkers organize nearby dynamics in much the same way as Abresch–Langer curves do for curve-shortening flow.

## 6. Variants, nonstandard area-preserving flows, and affine analogues

The term “area-preserving curve-shortening flow” is most often reserved for Gage’s nonlocal flow with Lagrange multiplier \(\bar\kappa\), but nearby constructions also appear in the literature and illuminate what is specific to APCSF.

A 2025 paper studies a different nonlocal area-preserving curvature flow with normal velocity
\[
\partial_t X=(\kappa-\lambda(t))\,\nu,\qquad 
\lambda(t)=\frac{L^2(t)-2\pi A(t)}{A(t)}.
\]
This flow also preserves area, decreases length, preserves convexity, and converges smoothly to a circle for strictly convex initial data [2502.16409]. It is not the classical APCSF because the Lagrange multiplier depends on \(L^2-2\pi A\) rather than average curvature. The comparison is useful because it isolates which analytical features are robust under nonlocal perimeter-constrained flows and which depend specifically on the average-curvature structure.

There are also non-Euclidean analogues. The paper “A fourth-order area-preserving curve flow in centro-equiaffine geometry” constructs a local, fourth-order, area-preserving flow in which the velocity is
\[
\partial_t\gamma=(1-\sigma\mathcal{K})\,\mathcal{N}=\sigma_{ss}\,\mathcal{N},
\]
using the affine Minkowski identity \(\sigma_{ss}+\sigma\mathcal{K}=1\) [2604.14804]. This flow preserves Euclidean enclosed area, exists globally for smooth strictly convex origin-symmetric initial curves, and converges smoothly to a round circle up to \(\mathrm{SL}(2)\) [2604.14804]. It is not Euclidean APCSF: it is fourth-order, local rather than nonlocal, and affine rather than Euclidean in structure. Nonetheless, it clarifies a broader theme: area preservation can be enforced either through a nonlocal Lagrange multiplier, as in APCSF, or through a local geometric identity built into the velocity law.

## 7. Numerical discretization and structure preservation

Numerical work on APCSF focuses on preserving the flow’s defining geometry at the discrete level: area conservation, perimeter decrease, convexity preservation, and mesh quality.

The 2022 paper “A convexity-preserving and perimeter-decreasing parametric finite element method for the area-preserving curve shortening flow” develops a semi-discrete PFEM based on Dziuk’s mass-lumped formulation [2208.01324]. In parametric form, the continuous APCSF is written as
\[
X_t=\frac{1}{|X_\xi|}\partial_\xi\!\left(\frac{X_\xi}{|X_\xi|}\right)
-\frac{2\pi}{L(t)}\Bigl(\frac{X_\xi}{|X_\xi|}\Bigr)^\perp,
\]
and the semi-discrete scheme is proved to preserve two fundamental geometric properties for initially convex curves: discrete convexity and perimeter decrease [2208.01324]. The paper also establishes an \(O(h)\) \(H^1\)-error estimate and reports approximately first-order convergence in the \(H^1\)-seminorm and close to second-order behavior in \(L^2\) and velocity errors in numerical tests [2208.01324].

More recent work emphasizes higher-order time integration and constraint handling. The 2026 dual-SAV PFEM framework treats APCSF as a constrained \(L^2\)-gradient flow of length with an area constraint enforced by a Lagrange multiplier, while using separate auxiliary variables for geometric energy and mesh regularization [2605.11721]. For one global area constraint, the nonlinear algebraic system per time step is reduced to dimension \(K+1=2\), independently of the number of mesh vertices [2605.11721]. Numerical experiments for APCSF show near machine-precision area conservation together with significant mesh-quality improvement under tangential regularization [2605.11721].

A complementary 2026 paper on arbitrary-order structure-preserving discretizations does not treat APCSF explicitly, but it develops a continuous Petrov–Galerkin framework that preserves area or volume laws for curvature-driven flows at arbitrary order in space and time [2605.20371]. It shows that conservation laws can be replicated discretely by introducing suitable auxiliary variables and test functions, and it explicitly notes that the same machinery can be adapted to APCSF by introducing a scalar Lagrange multiplier and an area constraint tested with \(1\) [2605.20371]. This suggests a route to genuinely high-order, structure-preserving APCSF solvers.

## 8. Relation to curve-shortening flow and broader geometric-flow theory

APCSF is best understood in comparison with classical curve-shortening flow. Under curve-shortening flow, convex embedded curves shrink to round points in finite time and type I self-similar shrinkers govern singularity models. Under APCSF, the area constraint prevents extinction and instead drives the evolution toward stationary constant-curvature states when global existence holds [2105.11930], [2412.18102]. For convex or suitably controlled nonconvex closed curves, this means convergence to a circle; for free-boundary problems, to a circular arc [1503.05814].

Yet the area constraint also creates new singular behavior. In free-boundary APCSF, finite-time singularities are type II and have grim-reaper rather than shrinking-soliton blow-up models [1701.08948]. In immersed closed settings, APCSF admits non-circular homothetic shrinkers with zero signed area [2508.00352]. These two facts together show that APCSF blends elliptic equilibrium selection with singular translator dynamics in a way that has no direct analogue in the simplest embedded-curve theory for classical curve-shortening flow.

A recurring misconception is that “area-preserving” implies benign long-time behavior for all smooth initial data. The literature does not support that view. For convex closed data, the behavior is global and circular; for star-shaped embedded data, global existence still holds [2412.18102]; but for general immersed or free-boundary data, finite-time singularities can occur, and explicit criteria are known [1701.08948], [2508.00352]. Another misconception is that APCSF is merely curve-shortening flow with a harmless lower-order correction. The classification of type II blow-ups and the existence of non-circular APCSF shrinkers indicate that the nonlocal term qualitatively changes both singularity theory and the global phase portrait.

A plausible implication of the current literature is that APCSF occupies an intermediate position between unconstrained second-order flows and higher-order area-preserving flows such as curve diffusion. It retains a strong maximum-principle component, especially for curvature and convexity arguments, but already exhibits the rich dynamical structure typical of constrained geometric evolution: nonlocality, multiple equilibria in immersed classes, and delicate interactions between topology, turning number, and singularity formation.

## 9. Outlook

Several directions remain active. The first is classification. For closed embedded star-shaped curves, global existence is now well understood [2412.18102], but for general immersed curves the complete APCSF singularity landscape remains only partially charted. The recent theory of APCSF shrinkers suggests an Abresch–Langer-type classification indexed by rational turning-number data, but uniqueness up to similarity is still conjectural [2508.00352].

A second direction concerns free boundaries and support geometry. The convex-support Neumann theory is comparatively mature [1503.05814], [1701.08948], but more general support manifolds, anisotropies, and higher-dimensional analogues remain open. The reflection argument at a straight line suggests that some closed-curve theory can be imported to boundary problems, yet the boundary singularity structure is genuinely richer.

A third direction is discretization. Recent PFEM and dual-SAV frameworks show that discrete area preservation, perimeter decrease, convexity preservation, and mesh regularization can be achieved simultaneously in practical schemes [2208.01324], [2605.11721]. The arbitrary-order structure-preserving methodology developed for related curvature flows indicates that high-order APCSF solvers with exact discrete conservation laws are now technically plausible [2605.20371].

In aggregate, APCSF has evolved from a constrained variant of curve-shortening flow into a distinct subject within geometric analysis: one in which nonlocality, isoperimetric optimization, singularity formation, and discrete structure preservation all play essential roles. Its central equations are simple, but the theory they generate is not; it includes global convergence theorems, translator blow-ups, free-boundary grim reapers, non-circular homothetic shrinkers, and numerically nontrivial conservation laws, all tied together by the single principle of preserving area while shortening curve length.

Source: https://www.emergentmind.com/topics/area-preserving-curve-shortening-flow-apcsf