pith. sign in

arxiv: 2607.02179 · v1 · pith:HAFDJMW7new · submitted 2026-07-02 · 🧮 math.NA · cs.NA

Unconditional Optimal Error Estimates and Energy Stability for a Linearly Implicit Mass-Lumped Projection Finite Element Method for the Harmonic Map Flow

Yongyong Cai , Xingwei Yang This is my paper

Pith reviewed 2026-07-03 07:50 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords harmonic map flowfinite element methodmass lumpingnodal projectionenergy stabilityoptimal error estimatesunconditional convergence
0
0 comments X

The pith

A linearly implicit mass-lumped projection method for harmonic map flow yields optimal error estimates without any time-step restriction relative to mesh size.

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

The paper develops a scheme that first solves a linear mass-lumped predictor derived from a cross-product form of the harmonic map equation and then applies a nodal projection to restore the exact unit-length constraint at every node. On Cartesian rectangular and cuboidal meshes this projection is nonexpansive in the discrete Dirichlet energy, producing an unconditional discrete energy decay law. The same contraction property together with quadrature estimates and edge-based cancellation identities yields optimal rates O(Δt + h²) in the discrete ℓ^∞(L²) norm and O(Δt + h) in the discrete ℓ²(H¹) seminorm for smooth solutions, with no coupling condition imposed between time step and mesh size.

Core claim

The linearly implicit mass-lumped predictor followed by nodal projection satisfies an exact nodal orthogonality relation that places the auxiliary solution on or outside the unit sphere, rendering the projection well-defined; on tensor-product meshes the projection is nonexpansive in the discrete energy, and the projected error contracts in the discrete L² inner product, enabling unconditional energy stability and the stated optimal convergence rates without time-mesh coupling.

What carries the argument

The nodal projection applied after the mass-lumped linear predictor, which contracts the projected error in the discrete L² norm and is nonexpansive in the discrete Dirichlet energy on Cartesian tensor-product meshes.

If this is right

  • The discrete solution satisfies the unit-length constraint exactly at every finite-element node for all time steps.
  • The scheme produces a monotonically nonincreasing discrete energy without any restriction on the ratio of time step to mesh size.
  • The predictor step is equivalent to both a corrected mass-lumped discretization and a tangent-plane scheme.
  • Error bounds hold in the stated norms for any sufficiently smooth exact solution of the continuous harmonic map flow.

Where Pith is reading between the lines

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

  • The edge-based cancellation identities introduced for the error analysis may transfer to other projection-based discretizations of geometric evolution equations.
  • If a comparable nonexpansiveness property can be established on unstructured meshes, the unconditional stability result would extend beyond tensor-product grids.
  • The absence of a time-step restriction suggests the method could remain stable and accurate for long-time simulations on fine spatial meshes where explicit or coupled schemes become impractical.

Load-bearing premise

The nodal projection is nonexpansive in the discrete Dirichlet energy on Cartesian rectangular and cuboidal tensor-product meshes.

What would settle it

A computation on a rectangular mesh with fixed large time step and successively refined spatial meshes that shows either discrete energy increase or observed convergence rates falling below O(Δt + h²) in L².

read the original abstract

We propose and analyze a linearly implicit mass-lumped finite element method for the heat flow of harmonic maps into the unit sphere. The method consists of a linear predictor followed by a nodal projection and therefore preserves the unit-length constraint exactly at all finite element nodes. The predictor is derived from a cross-product reformulation of the equation and is shown to be equivalent to a mass-lumped discretization of the original formulation with a correction term enforcing nodal orthogonality, as well as to a tangent plane scheme. A key ingredient is the consistent use of the discrete inner product in both the mass and stiffness terms. This yields a nodal orthogonality relation implying that the auxiliary solution lies on or outside the unit sphere at every node. Consequently, the projection is well defined and the projected error satisfies a contraction property in the discrete \(L^2\)-norm. On Cartesian rectangular and cuboidal tensor-product meshes, the nodal projection is also nonexpansive in a discrete Dirichlet energy, which gives an unconditional discrete energy dissipation law. For sufficiently smooth solutions, we prove optimal error estimates without any coupling condition between the time step and the mesh size: the method converges with order \(O(\Delta t+h^2)\) in \(\ell^\infty(0,T;L^2)\) and order \(O(\Delta t+h)\) in \(\ell^2(0,T;H^1)\). The proof combines the projected-error contraction, quadrature consistency estimates, edge-based cancellation identities, and a bootstrap argument for controlling nonlinear terms. Numerical experiments confirm the predicted convergence rates and the discrete energy decay.

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 proposes a linearly implicit mass-lumped projection finite element method for the heat flow of harmonic maps into the unit sphere. The scheme consists of a linear predictor derived from a cross-product reformulation of the PDE, followed by a nodal projection that enforces the unit-length constraint exactly at all nodes. The predictor is shown to be equivalent to a mass-lumped discretization of the original formulation with a correction term and to a tangent-plane scheme. Consistent use of the discrete inner product yields a nodal orthogonality relation, implying the auxiliary solution lies on or outside the unit sphere, so that the projection is well-defined and the projected error contracts in the discrete L² norm. On Cartesian rectangular and cuboidal tensor-product meshes the nodal projection is nonexpansive in the discrete Dirichlet energy, producing an unconditional discrete energy dissipation law. For sufficiently smooth solutions the method is proved to converge at the optimal rates O(Δt + h²) in ℓ^∞(0,T; L²) and O(Δt + h) in ℓ²(0,T; H¹) with no coupling restriction between time step and mesh size; the proof combines projected-error contraction, quadrature consistency, edge-based cancellation identities, and a bootstrap argument to control the nonlinear terms. Numerical experiments confirm the predicted rates and the discrete energy decay.

Significance. If the stated results hold, the contribution is significant: it supplies a linearly implicit scheme that exactly preserves the pointwise constraint at nodes, achieves unconditional energy stability on the indicated meshes, and delivers optimal convergence rates without any Δt–h coupling. Such unconditional stability and optimal rates are uncommon for constrained geometric evolution equations and would be of direct interest to the numerical analysis of harmonic map flows and related geometric PDEs. The combination of mass lumping, discrete-inner-product consistency, projected-error contraction, and bootstrap control constitutes a technically coherent proof strategy that advances the literature on structure-preserving discretizations.

minor comments (2)
  1. [Abstract] Abstract, paragraph on energy stability: the nonexpansiveness property and resulting unconditional energy law are stated only for Cartesian rectangular and cuboidal tensor-product meshes; a brief sentence clarifying that this mesh restriction does not affect the error estimates would improve readability.
  2. The abstract refers to 'edge-based cancellation identities' without a forward reference to the relevant lemma or equation; adding such a pointer would help readers trace the argument.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive evaluation and recommendation to accept the manuscript. The summary accurately captures the key contributions regarding the linearly implicit scheme, nodal constraint preservation, unconditional energy stability on tensor-product meshes, and optimal convergence rates without time-step restrictions.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The derivation chain proceeds from the cross-product reformulation of the harmonic map flow, the linearly implicit predictor, nodal projection step, and consistent discrete inner-product identities. These yield the nodal orthogonality relation, projected-error contraction in discrete L2, quadrature consistency, edge-based cancellation, and bootstrap control of nonlinear terms. The nonexpansiveness of the nodal projection (used for unconditional energy stability) is established as a mesh-specific property on Cartesian rectangular/cuboidal tensor-product meshes within the analysis, not imported via self-citation or ansatz. No parameter is fitted to data and then relabeled as a prediction; no uniqueness theorem or ansatz is smuggled through prior self-work; the O(Δt + h²) / O(Δt + h) rates follow from the stated hypotheses without reduction to inputs by construction. The argument is internally self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The paper introduces no free parameters or invented entities. It relies on standard finite-element approximation theory and smoothness assumptions for the error analysis.

axioms (2)
  • standard math Finite element spaces satisfy standard approximation and inverse inequalities used in error analysis.
    Invoked implicitly for the O(h²) and O(h) rates in the abstract.
  • domain assumption The exact solution is sufficiently smooth for the bootstrap argument and nonlinear term control to close.
    Explicitly stated as the condition under which the optimal rates hold.

pith-pipeline@v0.9.1-grok · 5826 in / 1448 out tokens · 31158 ms · 2026-07-03T07:50:43.387221+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

132 extracted references · 132 canonical work pages · 1 internal anchor

  1. [1]

    Gui, Xinping and Li, Buyang and Wang, Jilu , TITLE =. SIAM J. Numer. Anal. , FJOURNAL =. 2022 , NUMBER =

    work page 2022

  2. [2]

    Bramble, J. H. and Hilbert, S. R. , TITLE =. SIAM J. Numer. Anal. , FJOURNAL =. 1970 , PAGES =

    work page 1970

  3. [3]

    Li, Buyang and Sun, Weiwei , TITLE =. SIAM J. Numer. Anal. , FJOURNAL =. 2013 , NUMBER =

    work page 2013

  4. [4]

    and Bartels, S\"oren and Feng, Xiaobing and Prohl, Andreas , TITLE =

    Barrett, John W. and Bartels, S\"oren and Feng, Xiaobing and Prohl, Andreas , TITLE =. SIAM J. Numer. Anal. , FJOURNAL =. 2007 , NUMBER =

    work page 2007

  5. [5]

    Nie, Yi Yong and Thom\'ee, Vidar , TITLE =. IMA J. Numer. Anal. , FJOURNAL =. 1985 , NUMBER =

    work page 1985

  6. [6]

    Du, Qiang and Liu, Shuang and Yang, Jiang , TITLE =. J. Sci. Comput. , FJOURNAL =. 2025 , NUMBER =

    work page 2025

  7. [7]

    ESAIM Math

    Akrivis, Georgios and Bartels, S\"oren and Ruggeri, Michele and Wang, Jilu , TITLE =. ESAIM Math. Model. Numer. Anal. , FJOURNAL =. 2026 , NUMBER =

    work page 2026

  8. [8]

    Acta Mech

    Li, Xiwei and Zhang, Hanjie and Wang, Dongdong , TITLE =. Acta Mech. Sin. , FJOURNAL =. 2022 , NUMBER =

    work page 2022

  9. [9]

    ESAIM Math

    Gui, Xinping and Jiang, Yiwei and Wang, Jilu and Yin, Lihong , TITLE =. ESAIM Math. Model. Numer. Anal. , FJOURNAL =. 2025 , NUMBER =

    work page 2025

  10. [10]

    Bartels, S\"oren and Kov\'acs, Bal\'azs and Wang, Zhangxian , TITLE =. IMA J. Numer. Anal. , FJOURNAL =. 2024 , NUMBER =

    work page 2024

  11. [11]

    Bartels, S\"oren , TITLE =. Math. Comp. , FJOURNAL =. 2016 , NUMBER =

    work page 2016

  12. [12]

    A convergent implicit finite element discretization of the

    Ba. A convergent implicit finite element discretization of the. SIAM J. Numer. Anal. , FJOURNAL =. 2008 , NUMBER =

    work page 2008

  13. [13]

    Chen, Yun Mei and Struwe, Michael , TITLE =. Math. Z. , FJOURNAL =. 1989 , NUMBER =

    work page 1989

  14. [14]

    Chen, Yunmei and Lin, Fang Hua , TITLE =. J. Geom. Anal. , FJOURNAL =. 1998 , NUMBER =

    work page 1998

  15. [15]

    Di Fratta, Giovanni and Pfeiler, Carl-Martin and Praetorius, Dirk and Ruggeri, Michele , TITLE =. Comput. Methods Appl. Math. , FJOURNAL =. 2023 , NUMBER =

    work page 2023

  16. [16]

    Hu, Qiya and Tai, Xue-Cheng and Winther, Ragnar , TITLE =. SIAM J. Numer. Anal. , FJOURNAL =. 2009 , NUMBER =

    work page 2009

  17. [17]

    An, Rong and Gao, Huadong and Sun, Weiwei , TITLE =. SIAM J. Numer. Anal. , FJOURNAL =. 2021 , NUMBER =

    work page 2021

  18. [18]

    Chang, Kung-Ching and Ding, Wei Yue and Ye, Rugang , TITLE =. J. Differential Geom. , FJOURNAL =. 1992 , NUMBER =

    work page 1992

  19. [19]

    2026 , eprint=

    Stability and error analysis of fully discrete original energy-dissipative and length-preserving scheme for the Landau-Lifshitz-Gilbert equation , author=. 2026 , eprint=

    work page 2026

  20. [20]

    [2021] 2021 , PAGES =

    Ern, Alexandre and Guermond, Jean-Luc , TITLE =. [2021] 2021 , PAGES =

    work page 2021

  21. [21]

    Nonlinear partial differential equations in differential geometry (

    Struwe, Michael , TITLE =. Nonlinear partial differential equations in differential geometry (

    work page

  22. [22]

    and Dean, Joseph P

    Baratta, Igor A. and Dean, Joseph P. and Dokken, J. 2023 , howpublished =

    work page 2023

  23. [23]

    Bartels, S\"oren and Palus, Christian and Wang, Zhangxian , TITLE =. SIAM J. Numer. Anal. , FJOURNAL =. 2023 , NUMBER =

    work page 2023

  24. [24]

    Bartels, S\"oren , TITLE =. SIAM J. Numer. Anal. , FJOURNAL =. 2005 , NUMBER =

    work page 2005

  25. [25]

    Monk, P. B. and Vacus, O. , TITLE =. SIAM J. Numer. Anal. , FJOURNAL =. 1999 , NUMBER =

    work page 1999

  26. [26]

    A convergent linear finite element scheme for the

    Ba. A convergent linear finite element scheme for the. Electron. Trans. Numer. Anal. , FJOURNAL =. 2015 , PAGES =

    work page 2015

  27. [27]

    An, Rong and Sun, Weiwei , TITLE =. IMA J. Numer. Anal. , FJOURNAL =. 2022 , NUMBER =

    work page 2022

  28. [28]

    and Rodrigo, Jos\'e

    Robinson, James C. and Rodrigo, Jos\'e. The three-dimensional. 2016 , PAGES =

    work page 2016

  29. [29]

    2001 , PAGES =

    Prohl, Andreas , TITLE =. 2001 , PAGES =

    work page 2001

  30. [30]

    E, Weinan and Wang, Xiao-Ping , TITLE =. SIAM J. Numer. Anal. , FJOURNAL =. 2000 , NUMBER =

    work page 2000

  31. [31]

    Discrete Contin

    Alouges, Fran cois , TITLE =. Discrete Contin. Dyn. Syst. Ser. S , FJOURNAL =. 2008 , NUMBER =

    work page 2008

  32. [32]

    Cai, Yongyong and Chen, Jingrun and Wang, Cheng and Xie, Changjian , TITLE =. Math. Methods Appl. Sci. , FJOURNAL =. 2023 , NUMBER =

    work page 2023

  33. [33]

    Akrivis, Georgios and Feischl, Michael and Kov\'acs, Bal\'azs and Lubich, Christian , TITLE =. Math. Comp. , FJOURNAL =. 2021 , NUMBER =

    work page 2021

  34. [34]

    Cai, Yongyong and Chen, Jingrun and Wang, Cheng and Xie, Changjian , TITLE =. J. Comput. Phys. , FJOURNAL =. 2022 , PAGES =

    work page 2022

  35. [35]

    Chen, Jingrun and Wang, Cheng and Xie, Changjian , TITLE =. Appl. Numer. Math. , FJOURNAL =. 2021 , PAGES =

    work page 2021

  36. [36]

    Yang, Lei and Chen, Jingrun and Hu, Guanghui , TITLE =. J. Comput. Phys. , FJOURNAL =. 2021 , PAGES =

    work page 2021

  37. [37]

    Yang, Xingwei and Huang, Pengzhan and He, Yinnian , TITLE =. Math. Comput. Simulation , FJOURNAL =. 2026 , PAGES =

    work page 2026

  38. [38]

    Cheng, Qing and Shen, Jie , TITLE =. SIAM J. Sci. Comput. , FJOURNAL =. 2023 , NUMBER =

    work page 2023

  39. [39]

    Guti\'errez-Santacreu, Juan Vicente and Restelli, Marco , TITLE =. SIAM J. Numer. Anal. , FJOURNAL =. 2017 , NUMBER =

    work page 2017

  40. [40]

    Bartels, S\"oren and Lubich, Christian and Prohl, Andreas , TITLE =. Math. Comp. , FJOURNAL =. 2009 , NUMBER =

    work page 2009

  41. [41]

    Li, Xiaoli and Shen, Jie and Zheng, Nan , TITLE =. J. Comput. Phys. , FJOURNAL =. 2026 , PAGES =

    work page 2026

  42. [42]

    An, Rong and Li, Yonglin and Sun, Weiwei , TITLE =. IMA J. Numer. Anal. , FJOURNAL =. 2025 , NUMBER =

    work page 2025

  43. [43]

    d'Aquino, Massimiliano and Serpico, Claudio and Miano, Giovanni , TITLE =. J. Comput. Phys. , FJOURNAL =. 2005 , NUMBER =

    work page 2005

  44. [44]

    Bartels, S\"oren and Ko, Joy and Prohl, Andreas , TITLE =. Math. Comp. , FJOURNAL =. 2008 , NUMBER =

    work page 2008

  45. [45]

    Bartels, S\"oren and Prohl, Andreas , TITLE =. SIAM J. Numer. Anal. , FJOURNAL =. 2006 , NUMBER =

    work page 2006

  46. [46]

    Alouges, Fran cois and Jaisson, Pascal , TITLE =. Math. Models Methods Appl. Sci. , FJOURNAL =. 2006 , NUMBER =

    work page 2006

  47. [47]

    Cimr\'ak, Ivan , TITLE =. Arch. Comput. Methods Eng. , FJOURNAL =. 2008 , NUMBER =

    work page 2008

  48. [48]

    Recent developments in the modeling, analysis, and numerics of ferromagnetism , JOURNAL =

    Kru. Recent developments in the modeling, analysis, and numerics of ferromagnetism , JOURNAL =. 2006 , NUMBER =

    work page 2006

  49. [49]

    and Joly, P

    Cohen, G. and Joly, P. and Roberts, J. E. and Tordjman, N. , TITLE =. SIAM J. Numer. Anal. , FJOURNAL =. 2001 , NUMBER =

    work page 2001

  50. [50]

    2008 , PAGES =

    Guo, Boling and Ding, Shijin , TITLE =. 2008 , PAGES =

    work page 2008

  51. [51]

    Differential Integral Equations , FJOURNAL =

    Carbou, Gilles and Fabrie, Pierre , TITLE =. Differential Integral Equations , FJOURNAL =. 2001 , NUMBER =

    work page 2001

  52. [52]

    Guo, Bo Ling and Hong, Min Chun , TITLE =. Calc. Var. Partial Differential Equations , FJOURNAL =. 1993 , NUMBER =

    work page 1993

  53. [53]

    Nonlinear Anal

    Alouges, Fran cois and Soyeur, Alain , TITLE =. Nonlinear Anal. , FJOURNAL =. 1992 , NUMBER =

    work page 1992

  54. [54]

    A Lagrangian formulation of the gyromagnetic equation of the magnetization field , author=. Phys. Rev. , volume=

    work page

  55. [55]

    Alouges, Fran cois , TITLE =. SIAM J. Numer. Anal. , FJOURNAL =. 1997 , NUMBER =

    work page 1997

  56. [56]

    Kov\'acs, Bal\'azs and Li, Buyang and Lubich, Christian , TITLE =. Numer. Math. , FJOURNAL =. 2019 , NUMBER =

    work page 2019

  57. [57]

    and Osher, Stanley J

    Vese, Luminita A. and Osher, Stanley J. , TITLE =. SIAM J. Numer. Anal. , FJOURNAL =. 2002 , NUMBER =

    work page 2002

  58. [58]

    Physikalische Zeitschrift der Sowjetunion , volume=

    On the theory of the dispersion of magnetic permeability in ferromagnetic bodies , author=. Physikalische Zeitschrift der Sowjetunion , volume=

    work page

  59. [59]

    2023 , month =

    Firedrake User Manual , author =. 2023 , month =

    work page 2023

  60. [60]

    Le, Kim-Ngan and Page, Marcus and Praetorius, Dirk and Tran, Thanh , TITLE =. Appl. Anal. , FJOURNAL =. 2015 , NUMBER =

    work page 2015

  61. [61]

    Zhang, Guo-Dong and He, Xiaoming and Yang, Xiaofeng , TITLE =. J. Comput. Phys. , FJOURNAL =

    work page

  62. [62]

    Gao, Huadong , TITLE =. SIAM J. Numer. Anal. , FJOURNAL =. 2014 , NUMBER =

    work page 2014

  63. [63]

    and Mulder, W

    Geevers, S. and Mulder, W. A. and van der Vegt, J. J. W. , TITLE =. SIAM J. Sci. Comput. , FJOURNAL =. 2018 , NUMBER =

    work page 2018

  64. [64]

    Feischl, Michael and Tran, Thanh , TITLE =. SIAM J. Numer. Anal. , FJOURNAL =. 2017 , NUMBER =

    work page 2017

  65. [65]

    Ding, Qianqian and Mao, Shipeng , TITLE =. J. Sci. Comput. , FJOURNAL =. 2020 , NUMBER =

    work page 2020

  66. [66]

    Appl. Numer. Math. , FJOURNAL =. 2021 , issn =

    work page 2021

  67. [67]

    Math. Comput. Simul. , FJOURNAL =. 2026 , issn =

    work page 2026

  68. [68]

    2003 , publisher=

    Sobolev Spaces , author=. 2003 , publisher=

    work page 2003

  69. [69]

    Coey, J. M. D. History of Magnetism and Basic Concepts. Handbook of Magnetism and Magnetic Materials. 2020

    work page 2020

  70. [70]

    Characteristics of strong ferromagnetic Josephson junctions with epitaxial barriers , author =. Phys. Rev. B , volume =. 2005 , month =

    work page 2005

  71. [71]

    2003 , month =

    Monk, Peter , title =. 2003 , month =

    work page 2003

  72. [72]

    Bohn, Jan and Feischl, Michael and Kov\'acs, Bal\'azs , TITLE =. Comput. Methods Appl. Math. , FJOURNAL =. 2023 , NUMBER =

    work page 2023

  73. [73]

    Le, Kim-Ngan and Tran, Thanh , TITLE =. Comput. Math. Appl. , FJOURNAL =. 2013 , NUMBER =

    work page 2013

  74. [74]

    Hu, Shuaifei and Zhao, Guomei and An, Rong , TITLE =. Comput. Math. Appl. , FJOURNAL =. 2022 , PAGES =

    work page 2022

  75. [75]

    Cimr\'ak, Ivan , TITLE =. Math. Methods Appl. Sci. , FJOURNAL =. 2007 , NUMBER =

    work page 2007

  76. [76]

    Guo, Boling and Su, Fengqiu , TITLE =. J. Partial Differential Equations , FJOURNAL =. 2001 , NUMBER =

    work page 2001

  77. [77]

    Boling, Guo and Fengqiu, Su , TITLE =. J. Math. Anal. Appl. , FJOURNAL =. 1997 , NUMBER =

    work page 1997

  78. [78]

    Cimr\'ak, Ivan , TITLE =. J. Math. Anal. Appl. , FJOURNAL =. 2007 , NUMBER =

    work page 2007

  79. [79]

    Numer. Math. , FJOURNAL =. 1980 , issn =

    work page 1980

  80. [80]

    Japan J. Appl. Math. , FJOURNAL =. 1985 , author =

    work page 1985

Showing first 80 references.