pith. sign in

arxiv: 2604.10936 · v1 · submitted 2026-04-13 · 🧮 math.NA · cs.NA

Convergence Analysis of the Hessian Discretisation Method for Fourth Order Semi-linear Elliptic Equations with General Source

Devika Shylaja This is my paper

Pith reviewed 2026-05-10 16:14 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Hessian Discretisation Methodfourth-order elliptic equationsnonconforming finite elementsconvergence analysissemilinear problemsgradient recoveryAdini element
0
0 comments X

The pith

The Hessian Discretisation Method unifies convergence analysis for several numerical schemes on fourth-order semilinear elliptic equations without extra regularity assumptions on the solution.

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

This paper develops a convergence theory for the Hessian Discretisation Method applied to fourth-order semilinear elliptic equations that include a trilinear nonlinearity and a general source term. The method functions as a single abstract framework that encompasses conforming finite elements, nonconforming elements such as the Adini scheme, and gradient recovery techniques. Error estimates with explicit orders of convergence are obtained for the Adini nonconforming finite element method and for gradient recovery methods for the first time. The proofs rest only on four abstract properties of the discretisation together with a companion operator and do not invoke additional regularity on the exact solution.

Core claim

Convergence of the Hessian Discretisation Method for fourth-order semilinear elliptic equations follows from four key properties of the method together with a suitable companion operator, and this yields explicit error bounds that apply uniformly to the Adini nonconforming finite element method and to gradient recovery methods without requiring higher regularity assumptions on the exact solution.

What carries the argument

Hessian Discretisation Method, an abstract framework that encodes discrete approximations to the Hessian through four key properties and a companion operator to obtain unified error analysis.

Load-bearing premise

The chosen numerical schemes satisfy the four key properties of the Hessian Discretisation Method and admit a suitable companion operator.

What would settle it

A numerical test that computes the Adini discrete solution for a fourth-order semilinear problem with a known smooth exact solution and checks whether the error in the discrete norm decreases at the rate claimed by the analysis.

read the original abstract

This paper presents a convergence analysis for the Hessian Discretisation Method (HDM) applied to fourth-order semilinear elliptic equations involving a trilinear nonlinearity and general source, based on two complementary approaches. The HDM serves as a unified framework for the convergence analysis of various numerical schemes, including conforming and nonconforming finite element methods (ncFEMs) and gradient recovery (GR) based methods. Error estimates for the Adini ncFEM and GR methods are derived for the first time, which provide an explicit order of convergence. The analysis relies on four key HDM properties along with a suitable companion operator to establish convergence results. Moreover, a convergence analysis is developed within the HDM framework, which does not require additional regularity assumptions on the exact solution or the assumption that the exact solution is regular. The paper further discusses two significant applications: the Navier--Stokes equations in stream function--vorticity formulation and the von K\'{a}rm\'{a}n equations for plate bending. Numerical experiments are provided to demonstrate the performance of the GR method, Morley, and Adini ncFEMs.

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

3 major / 3 minor

Summary. The paper develops a convergence analysis within the Hessian Discretisation Method (HDM) framework for fourth-order semilinear elliptic equations with trilinear nonlinearity and general source term. It relies on four abstract HDM properties plus a companion operator to obtain error estimates (with explicit orders) for the Adini nonconforming FEM and gradient-recovery methods for the first time, without imposing extra regularity on the exact solution. The framework is applied to the Navier-Stokes equations in stream-function/vorticity form and the von Kármán plate equations; numerical experiments illustrate the performance of the GR, Morley, and Adini schemes.

Significance. If the four HDM properties and companion operator are verified to hold uniformly for the trilinear nonlinearity and general source, the work supplies a unified, regularity-free convergence theory that covers previously unanalyzed discretizations and yields explicit rates. This would be a useful contribution to the numerical analysis of fourth-order nonlinear problems, particularly for the cited applications.

major comments (3)
  1. [§3] §3 (HDM properties for the semilinear problem): The four key properties (consistency, stability, approximation, and the companion-operator relation) are invoked to pass from the linear biharmonic case to the trilinear nonlinearity, yet the manuscript does not supply a uniform-in-nonlinearity verification that the consistency error remains controlled by the same mesh-size powers when the right-hand side is merely in L² or H^{-2}. This verification is load-bearing for the claim of “no additional regularity assumptions.”
  2. [Theorem 4.1] Theorem 4.1 (error estimate for Adini ncFEM): The proof reduces the nonlinear residual to the linear HDM estimate via the companion operator, but the Lipschitz constant of the trilinear form is allowed to depend on the discrete solution; no a-priori bound independent of h is derived that would justify passing to the limit without extra regularity on the exact solution.
  3. [§5.2] §5.2 (gradient-recovery method): The construction of the companion operator for the GR scheme is stated to satisfy the required commutation property, but the proof that this operator remains bounded independently of the trilinear term (when the source is general) is omitted; the explicit convergence order claimed in the abstract therefore rests on an unverified step.
minor comments (3)
  1. [§2] The notation for the discrete trilinear form and the companion operator should be introduced once in §2 and used consistently thereafter; several ad-hoc symbols appear only in the proofs.
  2. [§6] Figure 1 (numerical convergence plots) lacks error bars or tabulated rates; adding a table of observed orders for the three methods would strengthen the experimental section.
  3. [§6] The statement that the analysis applies to “general source” should be accompanied by an explicit example (e.g., f ∈ H^{-2} but not in L²) in the numerical section to illustrate the regularity-free claim.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the detailed and insightful comments on our manuscript. We appreciate the recognition of the potential contribution of our unified HDM framework for fourth-order semilinear problems. Below, we provide point-by-point responses to the major comments and outline the revisions we plan to make to strengthen the paper.

read point-by-point responses

  1. Referee: [§3] §3 (HDM properties for the semilinear problem): The four key properties (consistency, stability, approximation, and the companion-operator relation) are invoked to pass from the linear biharmonic case to the trilinear nonlinearity, yet the manuscript does not supply a uniform-in-nonlinearity verification that the consistency error remains controlled by the same mesh-size powers when the right-hand side is merely in L² or H^{-2}. This verification is load-bearing for the claim of “no additional regularity assumptions.”

    Authors: We agree that explicit uniform verification is important for the claim of no additional regularity. In Section 3, the consistency error for the semilinear term is estimated using the general source term and the properties of the companion operator, which ensures control by the same powers of h as in the linear case. However, to address this concern directly, we will add a new lemma in the revised manuscript that explicitly verifies the uniformity with respect to the nonlinearity, showing that the consistency error bound holds independently of the discrete solution under the given assumptions on the source term. revision: yes

  2. Referee: [Theorem 4.1] Theorem 4.1 (error estimate for Adini ncFEM): The proof reduces the nonlinear residual to the linear HDM estimate via the companion operator, but the Lipschitz constant of the trilinear form is allowed to depend on the discrete solution; no a-priori bound independent of h is derived that would justify passing to the limit without extra regularity on the exact solution.

    Authors: In the proof of Theorem 4.1, the dependence of the Lipschitz constant on the discrete solution is handled by first establishing an a priori bound on the discrete solution using the stability property of the HDM (Property 2) and the coercivity of the bilinear form, combined with the boundedness of the trilinear form for solutions in the discrete space. This bound is independent of the mesh size h because it relies on the continuous embedding and the general source term in L² or H^{-2}. We will revise the proof to explicitly derive and state this a priori estimate before applying the Lipschitz continuity, thereby justifying the passage to the limit without requiring extra regularity on the exact solution. revision: partial

  3. Referee: [§5.2] §5.2 (gradient-recovery method): The construction of the companion operator for the GR scheme is stated to satisfy the required commutation property, but the proof that this operator remains bounded independently of the trilinear term (when the source is general) is omitted; the explicit convergence order claimed in the abstract therefore rests on an unverified step.

    Authors: We acknowledge that the boundedness of the companion operator for the gradient-recovery scheme, independent of the trilinear term, was not fully detailed in Section 5.2. The construction ensures the commutation property, and boundedness follows from the recovery operator's properties and the general source. In the revised version, we will provide the complete proof of this boundedness, demonstrating that it holds uniformly and supports the explicit convergence order without additional regularity assumptions. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation rests on independent HDM properties and companion operator

full rationale

The paper states that its error estimates and regularity-free convergence follow from four key HDM properties plus a companion operator applied to the semilinear problem. These are treated as given inputs that the chosen schemes (Adini ncFEM, GR) are asserted to satisfy; the subsequent analysis derives bounds from them without redefining the properties in terms of the derived error orders, without fitting parameters to the target quantities, and without load-bearing self-citations that close the loop. No equation reduces the output to the input by construction, and the framework is presented as self-contained once the abstract properties are granted.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Convergence depends on four unspecified key HDM properties and existence of a companion operator; no free parameters or new entities introduced in the abstract.

axioms (2)
  • domain assumption Four key HDM properties hold for the discretisations considered
    Invoked to establish convergence results for the schemes.
  • domain assumption A suitable companion operator exists
    Used together with the four properties to prove convergence.

pith-pipeline@v0.9.0 · 5496 in / 1236 out tokens · 29475 ms · 2026-05-10T16:14:15.533934+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

31 extracted references · 31 canonical work pages

  1. [1]

    D. Adak, D. Mora, S. Natarajan, and A. Silgado. A virtual element discretization for the time dependent navier–stokes equations in stream-function formulation.ESAIM: Mathematical Modelling and Numerical Analysis, 55(5):2535–2566, 2021

    work page 2021

  2. [2]

    Blum and R

    H. Blum and R. Rannacher. On the boundary value problem of the biharmonic operator on domains with angular corners.Math. Methods Appl. Sci., 2(4):556–581, 1980

    work page 1980

  3. [3]

    S. C. Brenner. A two-level additive Schwarz preconditioner for nonconforming plate elements. Numer. Math., 72(4):419–447, 1996

    work page 1996

  4. [4]

    S. C. Brenner, M. Neilan, A. Reiser, and L.-Y. Sung. A𝐶0 interior penalty method for a von K´arm´an plate.Numer. Math., 135(3):803–832, 2017

    work page 2017

  5. [5]

    S. C. Brenner, L.-Y. Sung, H. Zhang, and Y. Zhang. A Morley finite element method for the displacement obstacle problem of clamped Kirchhoff plates.J. Comput. Appl. Math., 254:31–42, 2013

    work page 2013

  6. [6]

    F. Brezzi. Finite element approximations of the von K ´arm´an equations.RAIRO Anal. Num ´er., 12(4):303–312, 1978

    work page 1978

  7. [7]

    Brezzi, J

    F. Brezzi, J. Rappaz, and P.-A. Raviart. Finite-dimensional approximation of nonlinear problems. I. Branches of nonsingular solutions.Numer. Math., 36(1):1–25, 1980/81

    work page 1980

  8. [8]

    Carstensen, G

    C. Carstensen, G. Mallik, and N. Nataraj. A priori and a posteriori error control of discontinuous Galerkin finite element methods for the von K´arm´an equations.IMA J. Numer. Anal., 39(1):167–200, 2019

    work page 2019

  9. [9]

    Carstensen, N

    C. Carstensen, N. Nataraj, G. C. Remesan, and D. Shylaja. Unified a priori analysis of four second- order fem for fourth-order quadratic semilinear problems.Numerische Mathematik, 154(2):323–368, 2023

    work page 2023

  10. [10]

    M. E. Cayco and R. A. Nicolaides. Finite element technique for optimal pressure recovery from stream function formulation of viscous flows.Mathematics of Computation, 46(174):371–377, 1986

    work page 1986

  11. [11]

    P. G. Ciarlet.The finite element method for elliptic problems. Access Online via Elsevier, 1978

    work page 1978

  12. [12]

    P. G. Ciarlet.Mathematical elasticity. Vol. II, volume 27 ofStudies in Mathematics and its Applica- tions. North-Holland Publishing Co., Amsterdam, 1997. Theory of plates

    work page 1997

  13. [13]

    Droniou, R

    J. Droniou, R. Eymard, T. Gallou¨et, C. Guichard, and R. Herbin.The gradient discretisation method, volume 82 ofMath ´ematiques & Applications (Berlin) [Mathematics & Applications]. Springer, Cham, 2018

    work page 2018

  14. [14]

    Droniou, B

    J. Droniou, B. Lamichhane, and D. Shylaja. The Hessian discretisation method for fourth order linear elliptic equations.J. Sci. Comput., 78(3):1405–1437, 2019

    work page 2019

  15. [15]

    Droniou, N

    J. Droniou, N. Nataraj, and D. Shylaja. Hessian discretisation method for fourth-order semi-linear elliptic equations: applications to the von k ´arm´an and navier–stokes models.Advances in Compu- tational Mathematics, 47(2):20, 2021. Published online: 16 February 2021. REFERENCES 27

    work page 2021

  16. [16]

    L. C. Evans.Partial differential equations, volume 19 ofGraduate Studies in Mathematics. American Mathematical Society, Providence, RI, 1998

    work page 1998

  17. [17]

    E. L. Foster, T. Iliescu, and Z. Wang. A finite element discretization of the stream function formulation of the stationary quasigeostrophic equations of the ocean.Computer Methods in Applied Mechanics and Engineering, 261–262:105–117, 2013

    work page 2013

  18. [18]

    Girault and P.-A

    V. Girault and P.-A. Raviart.Finite element approximation of the Navier-Stokes equations, volume 749 ofLecture Notes in Mathematics. Springer-Verlag, Berlin-New York, 1979

    work page 1979

  19. [19]

    Grisvard.Singularities in boundary value problems, volume 22 ofRecherches en Math ´ematiques Appliqu´ees [Research in Applied Mathematics]

    P. Grisvard.Singularities in boundary value problems, volume 22 ofRecherches en Math ´ematiques Appliqu´ees [Research in Applied Mathematics]. Masson, Paris; Springer-Verlag, Berlin, 1992

    work page 1992

  20. [20]

    B. P. Lamichhane. A mixed finite element method for the biharmonic problem using biorthogonal or quasi-biorthogonal systems.J. Sci. Comput., 46(3):379–396, 2011

    work page 2011

  21. [21]

    B. P. Lamichhane. A finite element method for a biharmonic equation based on gradient recovery operators.BIT, 54(2):469–484, 2014

    work page 2014

  22. [22]

    Lions.Quelques m ´ethodes de r ´esolution des probl `emes aux limites non lin ´eaires

    J.-L. Lions.Quelques m ´ethodes de r ´esolution des probl `emes aux limites non lin ´eaires. Dunod; Gauthier-Villars, Paris, 1969

    work page 1969

  23. [23]

    Lovadina, D

    C. Lovadina, D. Mora, and I. Vel ´asquez. A virtual element method for the von k ´arm´an equations. ESAIM: Mathematical Modelling and Numerical Analysis, 55(2):533–560, 2021

    work page 2021

  24. [24]

    Mallik and N

    G. Mallik and N. Nataraj. Conforming finite element methods for the von K ´arm´an equations.Adv. Comput. Math., 42(5):1031–1054, 2016

    work page 2016

  25. [25]

    Mallik and N

    G. Mallik and N. Nataraj. A nonconforming finite element approximation for the von K ´arm´an equations.ESAIM Math. Model. Numer. Anal., 50(2):433–454, 2016

    work page 2016

  26. [26]

    T. Miyoshi. A mixed finite element method for the solution of the von K ´arm´an equations.Numer. Math., 26(3):255–269, 1976

    work page 1976

  27. [27]

    D. Shylaja. Improved𝐿 2 and𝐻 1 error estimates for the Hessian discretization method.Numer. Methods Partial Differential Equations, 2020

    work page 2020

  28. [28]

    D. Shylaja. Numerical analysis of optimal control problems governed by fourth-order linear elliptic equations using the hessian discretization method.Optimal Control Applications and Methods, 45:2027 – 2050, 2022

    work page 2027

  29. [29]

    Shylaja and S

    D. Shylaja and S. Kumar. Morley type virtual element method for von k´arm´an equations.Advances in Computational Mathematics, 50(94), 2024

    work page 2024

  30. [30]

    Zhang, Y

    B. Zhang, Y. Yang, and M. Feng. A𝐶0-weak galerkin finite element method for the two-dimensional navier-stokes equations in stream-function formulation.Journal of Computational Mathematics, 38(2):310–336, 2020

    work page 2020

  31. [31]

    Zhang and M

    X. Zhang and M. Feng. A mixed virtual element method for the two-dimensional navier-stokes equations in stream-function formulation.Computers & Mathematics with Applications, 175:119– 137, 2024

    work page 2024