arxiv: 2511.16070 · v2 · submitted 2025-11-20 · 🧮 math.NA · cs.NA

Optimal error analysis of an interior penalty virtual element method for fourth-order singular perturbation problems

Fang Feng , Yuanyi Sun , Yue Yu This is my paper

Pith reviewed 2026-05-17 21:27 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords interior penalty virtual element methodfourth-order singular perturbationoptimal error estimatesuniform convergenceboundary layersvirtual element methodssingularly perturbed problems

0 comments

The pith

The interior penalty virtual element method achieves optimal uniform error estimates for fourth-order singular perturbation problems even with boundary layers.

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

This paper shows that an interior penalty virtual element method for fourth-order singular perturbation problems attains optimal-order convergence that remains uniform with respect to the small perturbation parameter. Prior analyses had established only a suboptimal half-order uniform rate; the current work upgrades this to the full optimal rate by a more precise error analysis. The results continue to hold when thin boundary layers appear in the solution. A reader would care because such problems arise in many applications where the perturbation parameter can be arbitrarily small, and uniform optimality removes the need for special layer-adapted meshes.

Core claim

We demonstrate that the proposed IPVEM in fact achieves optimal and uniform error estimates, even in the presence of boundary layers. The theoretical results are substantiated through extensive numerical experiments, which confirm the validity of the error estimates and highlight the method's effectiveness for singularly perturbed problems.

What carries the argument

The interior penalty virtual element method equipped with penalty terms that weakly enforce inter-element continuity and control the fourth-order operator, enabling an error analysis that separates the contribution from smooth regions and from boundary layers while remaining independent of the perturbation parameter.

If this is right

Error bounds of full optimal order hold uniformly in the perturbation parameter.
The method converges optimally on quasi-uniform meshes without layer adaptation.
Numerical experiments confirm the predicted rates for a wide range of perturbation values.
The analysis applies directly to solutions containing boundary layers.

Where Pith is reading between the lines

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

The same technique may yield uniform optimality for virtual element discretizations of other high-order singularly perturbed equations.
Related nonconforming or discontinuous Galerkin methods could be re-examined for hidden optimal uniform rates.
The result suggests that computational cost for stiff fourth-order problems can be reduced by avoiding parameter-dependent mesh refinement.

Load-bearing premise

The analysis assumes the IPVEM formulation and penalty terms are exactly as defined in the cited prior works, with solution regularity sufficient away from layers.

What would settle it

A set of numerical tests on successively refined meshes showing that the observed convergence order falls below the expected optimal rate once the perturbation parameter becomes sufficiently small would falsify the uniform optimality claim.

read the original abstract

In recent studies \cite{ZZ24, FY24}, the Interior Penalty Virtual Element Method (IPVEM) has been developed for solving a fourth-order singular perturbation problem, with uniform convergence established in the lowest-order case concerning the perturbation parameter. However, the resulting uniform convergence rate is only of half-order, which is suboptimal. In this work, we demonstrate that the proposed IPVEM in fact achieves optimal and uniform error estimates, even in the presence of boundary layers. The theoretical results are substantiated through extensive numerical experiments, which confirm the validity of the error estimates and highlight the method's effectiveness for singularly perturbed problems.

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.

Desk Editor's Note private letter to a colleague

This paper improves the uniform convergence analysis for IPVEM on fourth-order singular perturbation problems from half-order to optimal order.

read the letter

The core advance is a sharper error analysis that delivers optimal uniform rates in the energy norm for the interior penalty virtual element method on fourth-order singularly perturbed problems, where ZZ24 and FY24 only reached half-order. The authors keep the same method but tighten the consistency and approximation estimates so the constants stay independent of the perturbation parameter ε. Numerical experiments on a range of ε values and mesh sizes back this up, showing the expected rates without visible degradation near layers. That is a useful incremental result for anyone already using VEM on biharmonic-type problems with boundary layers. The main soft spot is the reliance on solution regularity away from the layers; if the proof does not explicitly split the regular part from the layer part in the interpolation and consistency terms, an ε-dependent factor could still sneak in and cap the rate. The abstract does not spell out the precise mesh assumptions or the exact form of the penalty terms beyond referencing the prior papers, so a referee would need to check those details carefully. Overall the work is for specialists in virtual element methods and singularly perturbed PDEs. It is solid enough on its own terms to merit a serious referee rather than a desk reject.

Referee Report

2 major / 2 minor

Summary. The paper develops an interior penalty virtual element method (IPVEM) for fourth-order singular perturbation problems. It claims that, contrary to prior half-order uniform rates in ZZ24 and FY24, the method achieves optimal uniform error estimates independent of the perturbation parameter ε, even with boundary layers present. The theoretical results are supported by extensive numerical experiments confirming the rates.

Significance. If the optimal uniform estimates are rigorously established, the work would meaningfully advance virtual element methods for singularly perturbed higher-order PDEs by removing the suboptimality seen in earlier analyses. The combination of improved theory and numerical validation strengthens the case for IPVEM in layer-dominated regimes.

major comments (2)

[§4, Theorem 4.1] §4 (Error Analysis), particularly the consistency and approximation estimates leading to Theorem 4.1: the proof invokes generic H^s regularity away from layers but does not explicitly decompose the solution into regular and layer parts when bounding the interpolation and penalty terms. Without this split, layer derivatives of order O(ε^{-k/2}) can reintroduce ε-dependent factors, undermining the claimed uniformity and optimality of the rates.
[§3.2] §3.2 (IPVEM formulation): the penalty parameter choice and stabilization terms are stated to follow ZZ24/FY24 exactly, yet the improved rates rest on a sharper analysis of these terms under layer effects; the manuscript does not provide a separate lemma isolating the layer contribution to the consistency error.

minor comments (2)

[Tables 1-2] Table 1 and Table 2: the reported error tables would benefit from an additional column showing the observed rate as ε varies to make the uniformity visually immediate.
[§2] Notation: the symbol for the virtual element space and the broken Sobolev norm are introduced without a dedicated notation table; a short summary table would improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We address the major comments point by point below, indicating where revisions will be made to strengthen the presentation while preserving the core results.

read point-by-point responses

Referee: [§4, Theorem 4.1] §4 (Error Analysis), particularly the consistency and approximation estimates leading to Theorem 4.1: the proof invokes generic H^s regularity away from layers but does not explicitly decompose the solution into regular and layer parts when bounding the interpolation and penalty terms. Without this split, layer derivatives of order O(ε^{-k/2}) can reintroduce ε-dependent factors, undermining the claimed uniformity and optimality of the rates.

Authors: We appreciate the referee's observation on the structure of the proof. The current analysis employs approximation properties of the virtual element space and consistency estimates that are already formulated in norms robust to the singular perturbation; the penalty terms are scaled to absorb the large layer derivatives without introducing ε-factors into the final bound. Nevertheless, to make the uniformity fully explicit and address the concern directly, we will revise Section 4 to include an explicit decomposition of the solution into a smooth component and a boundary-layer component. Separate estimates for each part will then be derived, confirming that layer contributions remain controlled uniformly. This addition clarifies the argument without changing the stated theorem. revision: yes
Referee: [§3.2] §3.2 (IPVEM formulation): the penalty parameter choice and stabilization terms are stated to follow ZZ24/FY24 exactly, yet the improved rates rest on a sharper analysis of these terms under layer effects; the manuscript does not provide a separate lemma isolating the layer contribution to the consistency error.

Authors: The formulation in §3.2 indeed adopts the same parameter choices as ZZ24 and FY24. The improvement in rates arises from the refined bounding technique applied to the consistency error within the proof of Theorem 4.1, where layer effects are handled through the specific scaling of the interior-penalty terms. While we have not isolated this contribution in a standalone lemma, the estimates are carried out explicitly in the main argument. To improve readability, we will insert a short remark (or brief paragraph) immediately after the formulation in §3.2 that outlines how the layer contributions to consistency are controlled uniformly; this does not require a new numbered lemma but makes the sharper analysis more transparent. revision: partial

Circularity Check

0 steps flagged

Minor self-citation to prior IPVEM formulation; optimal uniform error rates derived via independent mathematical analysis

full rationale

The paper cites ZZ24 and FY24 solely for the IPVEM formulation and penalty terms, then performs a fresh error analysis to obtain optimal (rather than half-order) uniform rates. No step reduces the claimed rates to a fit, a self-citation chain, or a redefinition of inputs; the derivation rests on standard consistency, approximation, and stability arguments plus solution regularity away from layers. Numerical experiments supply external corroboration, keeping the central claim self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard assumptions from virtual element and finite element theory for fourth-order problems; no free parameters or new entities are introduced in the abstract.

axioms (1)

domain assumption Standard assumptions on virtual element spaces, penalty parameters, and solution regularity for IPVEM applied to singularly perturbed fourth-order equations
Invoked to obtain the uniform error bounds in the presence of boundary layers.

pith-pipeline@v0.9.0 · 5399 in / 1105 out tokens · 38871 ms · 2026-05-17T21:27:39.870344+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

14 extracted references · 14 canonical work pages

[1]

Ahmad, A

B. Ahmad, A. Alsaedi, F. Brezzi, L.D. Marini, and A. Russo . Equivalent projectors for virtual element methods. Comput. Math. Appl. , 66(3):376–391, 2013

work page 2013
[2]

Beirão Da Veiga, F

L. Beirão Da Veiga, F. Brezzi, A. Cangiani, G. Manzini, L. D. Marini, and A. Russo. Basic principles of virtual element methods. Math. Models Meth. Appl. Sci. , 23(1):199–214, 2013

work page 2013
[3]

Beirão Da Veiga, F

L. Beirão Da Veiga, F. Brezzi, L. D. Marini, and A. Russo. T he Hitchhiker’s guide to the virtual element method. Math. Models Meth. Appl. Sci. , 24(8):1541–1573, 2014

work page 2014
[4]

S. C. Brenner. Poincaré-Friedrichs inequalities for pi ecewise H 1 functions. SIAM J. Numer. Anal. , 41(1):306–324, 2003

work page 2003
[5]

S. C. Brenner and L. R. Scott. The Mathematical Theory of Finite Element Methods . Springer, New York, 2008

work page 2008
[6]

Brezzi, A

F. Brezzi, A. Buﬀa, and K. Lipnikov. Mimetic ﬁnite diﬀere nces for elliptic problems. M2AN Math. Model. Numer. Anal. , 43(2):277–295, 2009

work page 2009
[7]

Chen and J

L. Chen and J. Huang. Some error analysis on virtual eleme nt methods. Calcolo, 55(1):5, 2018

work page 2018
[8]

M. Chen, J. Huang, and X. Huang. A nonconforming ﬁnite ele ment method with the use of Nitsche’s technique for a strain gradient elastic model. J. Sci. Comput. , 104(1):Paper No. 2, 29, 2025

work page 2025
[9]

Fang and Y

F. Fang and Y. Yu. A modiﬁed interior penalty virtual elem ent method for fourth-order singular perturbation problems. J. Sci. Comput. , 101:Paper No. 21, 32, 2024

work page 2024
[10]

T. K. Nilssen, X. Tai, and R. Winther. A robust nonconfor ming H 2-element. Math. Comput. , 70(234):489–505, 2001. 10

work page 2001
[11]

J. Qiu, F. Wang, M. Ling, and J. Zhao. The interior penalt y virtual element method for the fourth-order elliptic hemivariational inequality. Commun. Nonlinear Sci. Numer. Simul. , 127(4644807):Paper No. 107547, 17, 2023

work page 2023
[12]

Talischi, G

C. Talischi, G. H. Paulino, A. Pereira, and Ivan F. M. Men ezes. Polymesher: a general-purpose mesh generator for polygonal elements written in Matlab. Struct. Multidiscip. Optim. , 45(3):309–328, 2012

work page 2012
[13]

Zhang and J

B. Zhang and J. Zhao. The virtual element method with int erior penalty for the fourth-order singular perturbation problem. Commun. Nonlinear Sci. Numer. Simul. , 133:Paper No. 107964, 2024

work page 2024
[14]

J. Zhao, S. Mao, B. Zhang, and F. Wang. The interior penal ty virtual element method for the biharmonic problem. Math. Comp. , 92(342):1543–1574, 2023. 11

work page 2023