pith. sign in

arxiv: 2507.00193 · v2 · submitted 2025-06-30 · 🧮 math.NA · cs.NA

An energy-stable parametric finite element method for Willmore flow with normal-tangential velocity splitting

Harald Garcke , Robert N\"urnberg , Quan Zhao This is my paper

Pith reviewed 2026-05-19 06:20 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Willmore flowparametric finite element methodenergy stabilitygeometric PDEnormal-tangential splittinghypersurface evolutionmean curvature
0
0 comments X

The pith

A parametric finite element method for Willmore flow achieves unconditional energy stability by splitting normal and tangential velocities.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a fully discrete parametric finite element approximation for Willmore flow of hypersurfaces that allows spontaneous curvature and open boundaries. It relies on a new geometric PDE which evolves the mean curvature to determine the normal velocity while using a separate equation to set the tangential velocity. This structure ensures that the discretization inherits an unconditional energy stability property. Well-posedness of the resulting algebraic system is shown, and the method is tested on curves in the plane and surfaces in space.

Core claim

The central discovery is a weak formulation of a geometric PDE that combines mean curvature evolution with tangential velocity prescription, which upon parametric finite element discretization yields a scheme with proven unconditional stability in the discrete Willmore energy.

What carries the argument

The normal-tangential velocity splitting in the geometric PDE for mean curvature evolution, which decouples the gradient flow direction from tangential motion to enable stable discretization.

Load-bearing premise

The reformulation of Willmore flow via the normal-tangential velocity split must preserve the underlying gradient flow structure for the discrete energy stability to follow unconditionally.

What would settle it

A counterexample or numerical run where the discrete energy fails to decrease or stay non-increasing over time steps would disprove the stability claim.

Figures

Figures reproduced from arXiv: 2507.00193 by Harald Garcke, Quan Zhao, Robert N\"urnberg.

Figure 2.1
Figure 2.1. Figure 2.1: Sketch of Γ(t) with boundary ∂Γ(t), as well as the three unit vectors ⃗τ, ⃗µ and ⃗ν which form a basis of R 3 . 2.1. The Willmore flow. For simplicity we define ⃗x(t) = ⃗x(·, t) and κ(·, t) = κ(t). In this work, we consider the energy contribution (2.5) Eκ(κ(t), ⃗x(t)) = 1 2 Z Γ(t) (κ − κ) 2 dHd−1 , κ ∈ R, where κ is the so-called spontaneous curvature and dHd−1 is the (d − 1)-dimensional Hausdorff measu… view at source ↗
Figure 5.1
Figure 5.1. Figure 5.1: Numerical errors for an expanding/shrinking circle with [PITH_FULL_IMAGE:figures/full_fig_p014_5_1.png] view at source ↗
Figure 5.2
Figure 5.2. Figure 5.2: [κ = −2] Numerical errors in the evolution of an initial circle segment un￾der Navier boundary conditions (left panel) and clamped boundary conditions (right panel), where ∆t = ( 2 5 h 5 ) 2 . We fix κ = −2 and consider both the cases of Navier boundary and clamped boundary conditions. As expected, a second-order convergence rate for the numerical solutions can be observed as well in [PITH_FULL_IMAGE:fi… view at source ↗
Figure 5.3
Figure 5.3. Figure 5.3: [κ = −0.5] Evolution of an initial circle segment towards the steady state (red line) under different boundary conditions. We plot Γm at times t = 0, 2, 4, · · · , 20, 50. On the bottom are plots of the discrete energy and the mesh ratio Rm. Example 3: In this example, we consider the evolution of an open curve under different types of boundary conditions. We also monitor the mesh ratio R m = maxσ∈T m Hd… view at source ↗
Figure 5.4
Figure 5.4. Figure 5.4: [κ = −2] Evolution of an initial circle segment towards the steady state (red line) under different boundary conditions. We plot Γm at times t = 0, 0.1, · · · , 0.6, 2. On the bottom are plots of the discrete energy and the mesh ratio Rm. -1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 Navier -1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 clamped -1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 mixed 0 0.2 0.4 0.6 0.8 1 t 0 5 10 15 20 Navier clamped m… view at source ↗
Figure 5.5
Figure 5.5. Figure 5.5: [κ = 2] Evolution of an initial circle segment towards the steady state (red line) under different boundary conditions. We plot Γm at times t = 0, 0.05, · · · , 0.3, 1. On the bottom are plots of the discrete energy and the mesh ratio Rm. Here the initial setting is the same as that in Example 2 for a circle segment. We fix J = 128 and ∆t = 10−3 and conduct an experiment with κ = −0.5, which will [PITH_… view at source ↗
Figure 5.6
Figure 5.6. Figure 5.6: [κ = 0] Evolution of an initial cigar shape of dimension 4×1×1 towards the steady state, where we plot Γm at t = 0, 1. On the bottom is a plot of the discrete energy. lead to the curve expanding. The numerical results are reported in [PITH_FULL_IMAGE:figures/full_fig_p017_5_6.png] view at source ↗
Figure 5.7
Figure 5.7. Figure 5.7: [κ = 0] Evolution of an initial thin torus towards a Clifford torus. We plot Γ m at t = 0, 1. On the bottom is a plot of the discrete energy. 4 × 1 × 1. We fix ∆t = 10−3 and choose κ = 0. Here J = 9216, K = 4610. The numerical results are shown in [PITH_FULL_IMAGE:figures/full_fig_p018_5_7.png] view at source ↗
Figure 5.8
Figure 5.8. Figure 5.8: [κ = −2] Evolution for a torus, where we plot Γm at t = 0, 1. On the bottom is a plot of the discrete energy. 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 t 50 100 150 200 energy [PITH_FULL_IMAGE:figures/full_fig_p019_5_8.png] view at source ↗
Figure 5.9
Figure 5.9. Figure 5.9: [κ = 1] Evolution for a torus, where we plot Γm at t = 0, 0.5. On the bottom is a plot of the discrete energy [PITH_FULL_IMAGE:figures/full_fig_p019_5_9.png] view at source ↗
Figure 5.10
Figure 5.10. Figure 5.10: Evolution of an initial sphere cap under Navier boundary conditions. Left [PITH_FULL_IMAGE:figures/full_fig_p020_5_10.png] view at source ↗
Figure 5.11
Figure 5.11. Figure 5.11: Evolution of an initial sphere cap under clamped boundary conditions. Left [PITH_FULL_IMAGE:figures/full_fig_p020_5_11.png] view at source ↗
Figure 5.12
Figure 5.12. Figure 5.12: Evolution of an initial sphere cap from a standard quadruple bubble with [PITH_FULL_IMAGE:figures/full_fig_p021_5_12.png] view at source ↗
Figure 5.13
Figure 5.13. Figure 5.13: Evolution of an initial spherical cap from a standard quadruple bubble with [PITH_FULL_IMAGE:figures/full_fig_p021_5_13.png] view at source ↗
Figure 5.14
Figure 5.14. Figure 5.14: Evolution of an initial cigar shape of dimension of 8 [PITH_FULL_IMAGE:figures/full_fig_p022_5_14.png] view at source ↗
read the original abstract

We propose and analyze an energy-stable fully discrete parametric approximation for Willmore flow of hypersurfaces in two and three space dimensions. We allow for the presence of spontaneous curvature effects and for open surfaces with boundary. The presented scheme is based on a new geometric partial differential equation (PDE) that combines an evolution equation for the mean curvature with a separate equation that prescribes the tangential velocity. The mean curvature is used to determine the normal velocity within the gradient flow structure, thus guaranteeing an unconditional energy stability for the discrete solution upon suitable discretization. We introduce a novel weak formulation for this geometric PDE, in which different types of boundary conditions can be naturally enforced. We further discretize the weak formulation to obtain a fully discrete parametric finite element method, for which well-posedness can be rigorously shown. Moreover, the constructed scheme admits an unconditional stability estimate in terms of the discrete energy. Extensive numerical experiments are reported to showcase the accuracy and robustness of the proposed method for computing Willmore flow of both curves in $\mathbb{R}^2$ and surfaces in $\mathbb{R}^3$.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 4 minor

Summary. The manuscript develops a fully discrete parametric finite element scheme for Willmore flow of hypersurfaces in R^2 and R^3, including spontaneous curvature and open surfaces with boundary. The central construction reformulates the gradient flow by introducing an auxiliary evolution equation for the mean curvature H together with a separate transport equation for the tangential velocity component. A weak formulation is derived whose discretization yields an unconditional discrete energy dissipation identity; well-posedness of the resulting nonlinear algebraic system is proved via a fixed-point argument that uses the energy bound to obtain a uniform a-priori estimate without mesh or time-step restrictions. Extensive numerical experiments are presented for curves and surfaces.

Significance. If the stability and well-posedness results hold, the work supplies a useful advance for long-time simulation of curvature-driven geometric flows. The unconditional energy stability, obtained directly from the discrete dissipation identity without artificial parameters, and the rigorous fixed-point solvability argument are genuine strengths. The method's ability to accommodate different boundary conditions through the weak form further increases its applicability. These features distinguish the contribution from many existing schemes that require CFL-type restrictions.

minor comments (4)
  1. [Introduction] The introduction would benefit from a short paragraph contrasting the new normal-tangential splitting with prior parametric approaches for Willmore flow (e.g., those based on direct discretization of the fourth-order equation), to make the novelty of the auxiliary H-equation clearer.
  2. [Section 3] In the weak formulation, the precise function space for the tangential velocity test functions should be stated explicitly so that readers can verify that the integration-by-parts step leading to the energy identity remains valid for all admissible test functions.
  3. [Section 5] Numerical experiments: convergence tables reporting observed rates in the L^2 and H^1 norms for successively refined meshes would strengthen the accuracy claims; currently only qualitative plots are shown.
  4. A few minor notation inconsistencies appear (e.g., the discrete energy is sometimes denoted E_h and sometimes E^n); a single consistent symbol throughout would improve readability.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive summary of the manuscript, for highlighting the unconditional energy stability and the rigorous fixed-point solvability argument as genuine strengths, and for recommending minor revision. We appreciate the recognition that the weak formulation naturally accommodates different boundary conditions and that the method is applicable to both curves and surfaces.

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper reformulates Willmore flow via an auxiliary evolution equation for mean curvature H combined with a separate transport equation for tangential velocity. The weak form is obtained by direct multiplication and integration by parts on this PDE, after which the discrete energy dissipation identity follows immediately upon substituting the discrete velocity as test function. Well-posedness of the resulting nonlinear algebraic system is established by a fixed-point argument that uses only the a-priori energy bound; no parameter fitting, self-referential definition, or load-bearing self-citation is required. The construction preserves equivalence to the original gradient flow without reducing any claimed prediction or stability result to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the validity of the newly proposed geometric PDE as a reformulation of Willmore flow and on standard results from finite element theory for geometric evolution equations; no additional free parameters or invented physical entities are introduced.

axioms (2)
  • domain assumption Existence and uniqueness results for weak solutions of the continuous geometric PDE
    Invoked to justify the subsequent discretization and stability analysis.
  • standard math Standard approximation properties of parametric finite element spaces for hypersurfaces
    Used to establish well-posedness of the fully discrete scheme.

pith-pipeline@v0.9.0 · 5725 in / 1364 out tokens · 47553 ms · 2026-05-19T06:20:56.251581+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

37 extracted references · 37 canonical work pages

  1. [1]

    Balzani and M

    N. Balzani and M. Rumpf , A nested variational time discretization for parametric Willmore flow, Interfaces Free Bound., 14 (2012), pp. 431–454

    work page 2012

  2. [2]

    Bao and Y

    W. Bao and Y. Li, An energy-stable parametric finite element method for the planar Willmore flow, SIAM J. Numer. Anal., 63 (2025), pp. 103–121

    work page 2025

  3. [3]

    J. W. Barrett, H. Garcke, and R. N ¨urnberg, A parametric finite element method for fourth order geometric evolution equations , J. Comput. Phys., 222 (2007), pp. 441–467

    work page 2007

  4. [4]

    , On the parametric finite element approximation of evolving hypersurfaces in R3, J. Comput. Phys., 227 (2008), pp. 4281–4307

    work page 2008

  5. [5]

    , Parametric approximation of Willmore flow and related geometric evolution equations , SIAM J. Sci. Comput., 31 (2008), pp. 225–253

    work page 2008

  6. [6]

    , Parametric approximation of surface clusters driven by isotropic and anisotropic sur- face energies, Interfaces Free Bound., 12 (2010), pp. 187–234

    work page 2010

  7. [7]

    , Computational parametric Willmore flow with spontaneous curvature and area differ- ence elasticity effects, SIAM J. Numer. Anal., 54 (2016), pp. 1732–1762

    work page 2016

  8. [8]

    , Stable variational approximations of boundary value problems for Willmore flow with Gaussian curvature, IMA J. Numer. Anal., 37 (2017), pp. 1657–1709

    work page 2017

  9. [9]

    , Parametric finite element approximations of curvature driven interface evolutions , Handb. Numer. Anal. (Andrea Bonito and Ricardo H. Nochetto, eds.), 21 (2020), pp. 275– 423. An energy-stable parametric FEM for Willmore flow 25

    work page 2020

  10. [10]

    Blatt, A singular example for the Willmore flow , Analysis (Munich), 29 (2009), p

    S. Blatt, A singular example for the Willmore flow , Analysis (Munich), 29 (2009), p. 407–430

    work page 2009

  11. [11]

    Bretin, S

    E. Bretin, S. Masnou, and E. Oudet , Phase-field approximations of the willmore functional and flow, Numer. Math., 131 (2015), pp. 115–171

    work page 2015

  12. [12]

    Clarenz, U

    U. Clarenz, U. Diewald, G. Dziuk, M. Rumpf, and R. Rusu , A finite element method for surface restoration with smooth boundary conditions , Comput. Aided Geom. Design, 21 (2004), pp. 427–445

    work page 2004

  13. [13]

    T. A. Davis, Algorithm 832: UMFPACK V4.3—an unsymmetric-pattern multifrontal method , ACM Trans. Math. Software, 30 (2004), pp. 196–199

    work page 2004

  14. [14]

    Droske and M

    M. Droske and M. Rumpf, A level set formulation for willmore flow , Interfaces Free Bound., 6 (2004), pp. 361–378

    work page 2004

  15. [15]

    Q. Du, C. Liu, and X. Wang , A phase field approach in the numerical study of the elastic bending energy for vesicle membranes , J. Comput. Phys., 198 (2004), pp. 450–468

    work page 2004

  16. [16]

    Duan and B

    B. Duan and B. Li , New artificial tangential motions for parametric finite element approxi- mation of surface evolution , SIAM J. Sci. Comput., 46 (2024), pp. A587–A608

    work page 2024

  17. [17]

    B. Duan, B. Li, and Z. Zhang , High-order fully discrete energy diminishing evolving surface finite element methods for a class of geometric curvature flows , Ann. Appl. Math, 37 (2021), pp. 405–436

    work page 2021

  18. [18]

    Dziuk, Computational parametric Willmore flow , Numer

    G. Dziuk, Computational parametric Willmore flow , Numer. Math., 111 (2008), pp. 55–80

    work page 2008

  19. [19]

    Dziuk, E

    G. Dziuk, E. Kuwert, and R. Sch ¨atzle, Evolution of elastic curves in Rn: Existence and computation, SIAM J. Math. Anal., 33 (2002), pp. 1228–1245

    work page 2002

  20. [20]

    Franken, M

    M. Franken, M. Rumpf, and B. Wirth , A phase field based PDE constrained optimization approach to time discrete Willmore flow , Int. J. Numer. Anal. Model., 10 (2013), pp. 116– 138

    work page 2013

  21. [21]

    Garcke, R

    H. Garcke, R. N ¨urnberg, and Q. Zhao , Arbitrary Lagrangian-Eulerian finite element ap- proximations for axisymmetric two-phase flow, Comput. Math. Appl., 155 (2024), pp. 209– 223

    work page 2024

  22. [22]

    arXiv:2503.23152, 2025

    , Stable fully discrete finite element methods with BGN tangential motion for Willmore flow of planar curves . arXiv:2503.23152, 2025

    work page arXiv 2025

  23. [23]

    Hu and B

    J. Hu and B. Li, Evolving finite element methods with an artificial tangential velocity for mean curvature flow and Willmore flow , Numer. Math., 152 (2022), pp. 127–181

    work page 2022

  24. [24]

    Kov´acs, B

    B. Kov´acs, B. Li, and C. Lubich, A convergent evolving finite element algorithm for Willmore flow of closed surfaces , Numer. Math., 149 (2021), pp. 595–643

    work page 2021

  25. [25]

    Kuwert and R

    E. Kuwert and R. Sch ¨atzle, The Willmore flow with small initial energy , J. Differ. Geom., 57 (2001), pp. 409–441

    work page 2001

  26. [26]

    , Gradient flow for the Willmore functional , Comm. Anal. Geom., 10 (2002), p. 307–339

    work page 2002

  27. [27]

    C. M. Elliott and H. Fritz, On approximations of the curve shortening flow and of the mean curvature flow based on the DeTurck trick , IMA J. Numer. Anal., 37 (2017), pp. 543–603

    work page 2017

  28. [28]

    F. C. Marques and A. Neves, Min-max theory and the Willmore conjecture, Ann. Math., 179 (2014), pp. 683–782

    work page 2014

  29. [29]

    U. F. Mayer and G. Simonett, A numerical scheme for axisymmetric solutions of curvature- driven free boundary problems, with applications to the Willmore flow , Interfaces Free Bound., 4 (2002), pp. 89–109

    work page 2002

  30. [30]

    Olischl ¨ager and M

    N. Olischl ¨ager and M. Rumpf , Two step time discretization of Willmore flow , in Mathe- matics of Surfaces XIII, E. R. Hancock, R. R. Martin, and M. A. Sabin, eds., Berlin, 2009, Springer, pp. 278–292

    work page 2009

  31. [31]

    Rumpf, J

    M. Rumpf, J. Sassen, and C. Smoch , A hybrid minimizing movement and neural network approach to Willmore flow . arXiv:2502.14656, 2025

    work page arXiv 2025

  32. [32]

    R. E. Rusu , An algorithm for the elastic flow of surfaces , Interfaces Free Bound., 7 (2005), pp. 229–239

    work page 2005

  33. [33]

    Schlierf, Spontaneous curvature effects of the Helfrich flow: singularities and convergence, Comm

    M. Schlierf, Spontaneous curvature effects of the Helfrich flow: singularities and convergence, Comm. Partial Diff. Equ., 50 (2025), p. 441–476

    work page 2025

  34. [34]

    Schmidt and K

    A. Schmidt and K. G. Siebert , Design of Adaptive Finite Element Software: The Finite Element Toolbox ALBERTA, vol. 42 of Lecture Notes in Computational Science and En- gineering, Springer-Verlag, Berlin, 2005

    work page 2005

  35. [35]

    Simonett, The Willmore flow near spheres , Differ

    G. Simonett, The Willmore flow near spheres , Differ. Integral Equ., 14 (2005), p. 1005–1014

    work page 2005

  36. [36]

    T. J. Willmore, Riemannian geometry, Oxford University Press, 1993

    work page 1993

  37. [37]

    Q. Zhao, W. Jiang, and W. Bao , An energy-stable parametric finite element method for simulating solid-state dewetting , IMA J. Numer. Anal., 41 (2021), pp. 2026–2055

    work page 2021