pith. sign in

arxiv: 2606.25934 · v1 · pith:YAAHVJESnew · submitted 2026-06-24 · 🧮 math.NA · cs.NA

Error estimates for A-stable backward difference full discretizations of Willmore flow of closed surfaces

Nils Bullerjahn , Bal\'azs Kov\'acs This is my paper

Pith reviewed 2026-06-25 19:34 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Willmore flowbackward difference methodsurface finite elementserror estimatesA-stabilityG-stabilitygeometric PDE
0
0 comments X

The pith

Optimal-order H¹-norm error estimates hold for A-stable BDF1 and BDF2 full discretizations of Willmore flow on closed surfaces.

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

The paper proves optimal-order error estimates in the H¹ norm for full discretizations of Willmore flow on closed two-dimensional surfaces. The methods combine surface finite elements of polynomial degree at least two in space with A-stable backward difference time stepping of order one or two. The proof proceeds from a stability analysis that uses energy estimates on the anti-symmetric structure of the second-order system, together with G-stability and multiplier techniques that include a new upper bound. A reader would care because the result justifies the use of these standard time integrators for a nonlinear geometric evolution equation while preserving the expected convergence rate.

Core claim

The authors prove that A-stable backward difference full discretizations of order 1 and 2 for Willmore flow of closed surfaces achieve optimal-order H¹-norm error estimates. The numerical method evolves surface finite elements of polynomial degree at least two in space together with backward difference time stepping. The convergence analysis rests on a stability analysis that uses energy estimates from the anti-symmetric structure of the second-order system, together with G-stability and multiplier techniques including a new upper bound.

What carries the argument

Stability analysis based on energy estimates exploiting the anti-symmetric structure of the second-order system, combined with Dahlquist's G-stability and the multiplier techniques of Nevanlinna and Odeh with a new upper bound.

If this is right

  • The full discretization converges at the optimal rate in the H¹ norm for both time orders considered.
  • The estimates apply directly to the coupled second-order system after surface finite element discretization.
  • Numerical experiments are expected to confirm the predicted rates on closed surfaces.
  • The analysis covers polynomial degrees at least two without order reduction in time.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The multiplier bound may extend the reach of G-stability arguments to other nonlinear geometric flows that admit an anti-symmetric formulation.
  • If the same structural properties hold for related higher-order systems, the approach could support error analysis of BDF methods beyond order two.
  • The result suggests that standard A-stable integrators remain reliable for long-time simulation of Willmore flow without custom stabilization.

Load-bearing premise

The stability analysis exploiting the anti-symmetric structure combined with G-stability and multiplier techniques with a new bound remains valid for the coupled Willmore system and its surface finite element discretization.

What would settle it

A numerical test on a sphere with a known exact solution where the observed H¹ error for the BDF2 scheme falls short of the predicted rate would falsify the claim.

read the original abstract

A proof of optimal-order $H^1$-norm error estimates is given for $A$-stable backward difference full discretizations (of order 1 and 2) of Willmore flow for closed two-dimensional surfaces. The numerical method discretizes a coupled system of evolution equations by evolving surface finite elements of polynomial degree at least two in space and backward difference method of order 1 or 2 in time. The convergence analysis is based on a stability analysis, based on energy estimates exploiting the anti-symmetric structure of the second-order system, in combination with Dahlquist's $G$-stability and the multiplier techniques of Nevanlinna and Odeh, with a new upper bound in the spirit of Dahlquist. Numerical experiments illustrate and complement the theoretical results.

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 / 2 minor

Summary. The paper claims to prove optimal-order H¹-norm error estimates for A-stable backward difference full discretizations (BDF1 and BDF2) of Willmore flow on closed two-dimensional surfaces. The method uses surface finite elements of polynomial degree at least two in space to discretize a coupled second-order system, with the convergence analysis relying on energy estimates that exploit the anti-symmetric structure, combined with Dahlquist G-stability and Nevanlinna-Odeh multiplier techniques together with a new upper bound.

Significance. If the analysis holds, the result supplies the first rigorous optimal-order error bounds for BDF time discretizations of Willmore flow, extending G-stability and multiplier techniques to a nonlinear geometric PDE setting. This is valuable for the numerical analysis of fourth-order geometric flows, where such stability tools are not routine.

minor comments (2)
  1. [Abstract] The abstract states that the new upper bound is 'in the spirit of Dahlquist,' but the precise statement and derivation of this bound should be isolated in a dedicated lemma or subsection to allow readers to verify its application to the coupled surface system.
  2. The description of the surface finite element discretization (degree at least two) and the handling of the evolving metric in the error analysis would benefit from an explicit statement of the interpolation estimates used, even if standard.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive summary, recognition of the significance of the first rigorous optimal-order error bounds for BDF discretizations of Willmore flow, and the recommendation of minor revision. No major comments are provided in the report, so we have no specific points requiring detailed rebuttal or revision at this stage.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The derivation rests on a stability analysis that combines the anti-symmetric structure of the second-order system with Dahlquist G-stability and Nevanlinna-Odeh multiplier techniques, plus a new upper bound derived in the paper itself. Error estimates then follow from standard consistency arguments for the surface finite element / BDF discretization. No step reduces a claimed result to a fitted parameter, self-defined quantity, or load-bearing self-citation; the cited stability tools are external and the new bound is presented as an original contribution within the work.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review; no explicit free parameters, invented entities, or non-standard axioms are stated. The proof implicitly relies on standard assumptions from finite element theory for parabolic systems and geometric flows.

pith-pipeline@v0.9.1-grok · 5665 in / 1082 out tokens · 24302 ms · 2026-06-25T19:34:56.756006+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

36 extracted references · 27 canonical work pages

  1. [1]

    Kov´ acs, B., Li, B., Lubich, C.: A convergent evolving finite element algorithm for Willmore flow of closed surfaces. Numer. Math. 149(3), 595–643 (2021) https://doi.org/10.1007/s00211-021-01238-z

    work page doi:10.1007/s00211-021-01238-z 2021

  2. [2]

    Al. I. Cu za

    Willmore, T.: Note on embedded surfaces. An. S ¸ti. Univ. “Al. I. Cu za” Ia¸ si Sect ¸. I a Mat. (N.S.) 11B, 493–496 (1965)

    1965

  3. [3]

    Oxford Science Publications, p

    Willmore, T.: Riemannian Geometry. Oxford Science Publications, p. 318. The Clarendon Press, Oxford University Press, New York (1993)

    1993

  4. [4]

    theory and possible experiments

    Helfrich, W.: Elastic properties of lipid bilayers: Theory and possible experiments. Zeitschrift f¨ ur Naturforschung C 28(11-12), 693–703 (1973) https://doi.org/10.1515/znc-1973-11-1209

    work page doi:10.1515/znc-1973-11-1209 1973

  5. [5]

    Elliott, C., Stinner, B.: Modeling and computation of two phase geom et- ric biomembranes using surface finite elements. J. Comput. Phys. 229(18), 6585–6612 (2010) https://doi.org/10.1016/j.jcp.2010.05.014

    work page doi:10.1016/j.jcp.2010.05.014 2010

  6. [6]

    Barrett, J., Garcke, H., N¨ urnberg, R.: Numerical computation s of the dynam- ics of fluidic membranes and vesicles. Phys. Rev. E 92, 052704 (2015) https://doi.org/10.1103/PhysRevE.92.052704

    work page doi:10.1103/physreve.92.052704 2015

  7. [7]

    Chen, Y., Lowengrub, J., Shen, J., Wang, C., Wise, S.: Efficient ener gy stable schemes for isotropic and strongly anisotropic Cahn-Hilliard s ystems with the Willmore regularization. J. Comput. Phys. 365, 56–73 (2018) https://doi.org/10.1016/j.jcp.2018.03.024

    work page doi:10.1016/j.jcp.2018.03.024 2018

  8. [8]

    In: Geometric Analysis

    Marques, F., Neves, A.: Min-max theory and a proof of the Willmore conjecture. In: Geometric Analysis. IAS/Park City Math. Ser., vol. 22, pp. 277– 300. Amer. Math. Soc., Providence, RI (2016). https://doi.org/10.1090/pcms/022/06

    work page doi:10.1090/pcms/022/06 2016

  9. [9]

    Interfaces Free Bound

    Mayer, U.F., Simonett, G.: A numerical scheme for axisymmetric so lutions of curvature-driven free boundary problems, with applications to th e Willmore flow. Interfaces Free Bound. 4(1), 89–109 (2002) https://doi.org/10.4171/IFB/54 32

    work page doi:10.4171/ifb/54 2002

  10. [10]

    In : Evolution Equations: Applications to Physics, Industry, Life Sciences and Ec onomics (Levico Terme, 2000)

    Mayer, U.F., Simonett, G.: Self-intersections for Willmore flow. In : Evolution Equations: Applications to Physics, Industry, Life Sciences and Ec onomics (Levico Terme, 2000). Progr. Nonlinear Differential Equations Appl., vol. 55 , pp. 341–348. Birkh¨ auser, Basel (2003)

    2000

  11. [11]

    , title =

    Rusu, R.: An algorithm for the elastic flow of surfaces. Interfa ces Free Bound. 7(3), 229–239 (2005) https://doi.org/10.4171/IFB/122

    work page doi:10.4171/ifb/122 2005

  12. [12]

    Dziuk, G.: Computational parametric Willmore flow. Numer. Math. 111(1), 55– 80 (2008) https://doi.org/10.1007/s00211-008-0179-1

    work page doi:10.1007/s00211-008-0179-1 2008

  13. [13]

    Barrett, J., Garcke, H., N¨ urnberg, R.: A parametric finite elem ent method for fourth order geometric evolution equations. J. Comput. Phys. 222(1), 441–462 (2007) https://doi.org/10.1016/j.jcp.2006.07.026

    work page doi:10.1016/j.jcp.2006.07.026 2007

  14. [14]

    Barrett, J., Garcke, H., N¨ urnberg, R.: Parametric approximation of Willmore flow and related geometric evolution equations. SIAM J. Sci. Comput. 31(1), 225–253 (2008) https://doi.org/10.1137/070700231

    work page doi:10.1137/070700231 2008

  15. [15]

    Kov´ acs, B., Li, B., Lubich, C.: A convergent evolving finite elemen t algorithm for mean curvature flow of closed surfaces. Numer. Math. 143(4), 797–853 (2019) https://doi.org/10.1007/s00211-019-01074-2

    work page doi:10.1007/s00211-019-01074-2 2019

  16. [16]

    Bullerjahn, N., Kov´ acs, B.: Error estimates for full discretiza tion of Cahn-Hilliard equation with dynamic boundary conditions. IMA J. Numer. Anal. 46(2), 713– 757 (2026) https://doi.org/10.1093/imanum/draf009

    work page doi:10.1093/imanum/draf009 2026

  17. [17]

    BIT 18(4), 384–401 (1978) https://doi.org/10.1007/BF01932018

    Dahlquist, G.: G-stability is equivalent to A-stability. BIT 18(4), 384–401 (1978) https://doi.org/10.1007/BF01932018

    work page doi:10.1007/bf01932018 1978

  18. [18]

    Nevanlinna, O., Odeh, F.: Multiplier techniques for linear multi- step methods. Numer. Funct. Anal. Optim. 3(4), 377–423 (1981) https://doi.org/10.1080/01630568108816097

    work page doi:10.1080/01630568108816097 1981

  19. [19]

    Lubich, C., Mansour, D., Venkataraman, C.: Backward differenc e time discretiza- tion of parabolic differential equations on evolving surfaces. IMA J. Numer. Anal. 33(4), 1365–1385 (2013) https://doi.org/10.1093/imanum/drs044

    work page doi:10.1093/imanum/drs044 2013

  20. [20]

    Akrivis, G., Lubich, C.: Fully implicit, linearly implicit and implicit-explicit backward difference formulae for quasi-linear parabolic equations. Numer. Math. 131(4), 713–735 (2015) https://doi.org/10.1007/s00211-015-0702-0

    work page doi:10.1007/s00211-015-0702-0 2015

  21. [21]

    Kov´ acs, B., Li, B., Lubich, C., Power Guerra, C.: Convergence o f finite elements on an evolving surface driven by diffusion on the surface. Numer. Ma th. 137(3), 643–689 (2017) https://doi.org/10.1007/s00211-017-0888-4 33

    work page doi:10.1007/s00211-017-0888-4 2017

  22. [22]

    Acta Numer

    Deckelnick, K., Dziuk, G., Elliott, C.: Computation of geometric par tial dif- ferential equations and mean curvature flow. Acta Numer. 14, 139–232 (2005) https://doi.org/10.1017/S0962492904000224

    work page doi:10.1017/s0962492904000224 2005

  23. [23]

    Progre ss in Nonlinear Differential Equations and Their Applications

    Ecker, K.: Regularity Theory for Mean Curvature Flow. Progre ss in Nonlinear Differential Equations and Their Applications. Birkh¨ auser, Boston (2012)

    2012

  24. [24]

    Advances in Design and Control

    Walker, S.: The Shapes of Things: A Practical Guide to Differentia l Geometry and the Shape Derivative. Advances in Design and Control. SIAM, So ciety for Industrial and Applied Mathematics, Philadelphia (2015)

    2015

  25. [25]

    Finite elements for the

    Dziuk, G.: Finite elements for the Beltrami operator on arbitrar y surfaces. In: Par- tial Differential Equations and Calculus of Variations. Lecture Note s in Math., vol. 1357, pp. 142–155. Springer, Berlin (1988). https://doi.org/10.1007/BFb0082865

    work page doi:10.1007/bfb0082865 1988

  26. [26]

    Demlow, A.: Higher-order finite element methods and pointwise er ror estimates for elliptic problems on surfaces. SIAM J. Numer. Anal. 47(2), 805–827 (2009) https://doi.org/10.1137/070708135

    work page doi:10.1137/070708135 2009

  27. [27]

    A cta Numer

    Dziuk, G., Elliott, C.M.: Finite element methods for surface PDEs. A cta Numer. 22, 289–396 (2013) https://doi.org/10.1017/S0962492913000056

    work page doi:10.1017/s0962492913000056 2013

  28. [28]

    Solving Ordinary Differential Equations

    Hairer, E., Wanner, G.: Solving Ordinary Differential Equations. I I, 2nd edn. Springer Series in Computational Mathematics, vol. 14, p. 614. Spr inger, Berlin (1996). https://doi.org/10.1007/978-3-642-05221-7 . Stiff and differential- algebraic problems

    work page doi:10.1007/978-3-642-05221-7 1996

  29. [29]

    BIT 60(1), 93–99 (2020) https://doi.org/10.1007/s10543-019-00768-1

    Akrivis, G., Katsoprinakis, E.: Maximum angles of A( ϑ)-stability of backward difference formulae. BIT 60(1), 93–99 (2020) https://doi.org/10.1007/s10543-019-00768-1

    work page doi:10.1007/s10543-019-00768-1 2020

  30. [30]

    Akrivis, G., Feischl, M., Kov´ acs, B., Lubich, C.: Higher-order linearly implicit full discretization of the Landau–Lifshitz–Gilbert equation. Math. Com p. 90(329), 995–1038 (2021)

    2021

  31. [31]

    Contri, A., Kov´ acs, B., Massing, A.: Error analysis of BDF 1–6 tim e-stepping methods for the transient Stokes problem: velocity and pressure estimates. SIAM J. Numer. Anal. 63(4), 1586–1616 (2025) https://doi.org/10.1137/23M1606800

    work page doi:10.1137/23m1606800 2025

  32. [32]

    arXiv (2025) 2502.03847 [math.NA]

    Bullerjahn, N.: Error estimates for full discretization by an almo st mass conser- vation technique for Cahn–Hilliard systems with dynamic boundary co nditions. arXiv (2025) 2502.03847 [math.NA]

    Pith/arXiv arXiv 2025

  33. [33]

    Akrivis, G., Katsoprinakis, E.: Backward difference formulae: Ne w multipliers and stability properties for parabolic equations. Math. Comp. 85(301), 2195–2216 (2015) 34

    2015

  34. [34]

    Akrivis, G., Chen, M., Yu, F., Zhou, Z.: The energy technique for t he six-step BDF method. SIAM J. Numer. Anal. 59(5), 2449–2472 (2021) https://doi.org/10.1137/21M1392656

    work page doi:10.1137/21m1392656 2021

  35. [35]

    ar Xiv:2605.14035, https://go.upb.de/ellFEM

    Kov´ acs, B., Lantelme, M.F.R.: ℓFEM: A fast and loop-free MATLAB implemen- tation of isoparametric bulk and surface finite elements (2026). ar Xiv:2605.14035, https://go.upb.de/ellFEM

    Pith/arXiv arXiv 2026

  36. [36]

    SIAM Rev

    Persson, P., Strang, G.: A simple mesh generator in Matlab. SIAM Rev. 46(2), 329–345 (2004) https://doi.org/10.1137/S0036144503429121 35

    work page doi:10.1137/s0036144503429121 2004