pith. sign in

arxiv: 2604.18406 · v1 · submitted 2026-04-20 · 🧮 math.NA · cs.NA

Virtual element methods for a quad-curl problem on general planar domains

Susanne C. Brenner , Li-yeng Sung , Jai Tushar This is my paper

Pith reviewed 2026-05-10 03:41 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords virtual element methodquad-curl problemHodge decompositionpolygonal domainsdivergence-free fieldsfinite element analysisnumerical PDEsplanar domains
0
0 comments X

The pith

Virtual element methods for the quad-curl problem on general polygonal domains are constructed and analyzed via the Hodge decomposition of divergence-free vector fields.

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

The paper develops virtual element methods to solve a quad-curl problem, a fourth-order vector PDE, on arbitrary polygonal domains rather than requiring structured meshes. The approach rests on discretizing the Hodge decomposition to enforce the divergence-free constraint inside the virtual element spaces. If the discretization preserves stability and approximation quality, the resulting schemes deliver reliable numerical solutions for problems on complex geometries where triangles or quadrilaterals are difficult to fit. Readers care because the method extends standard finite-element tools to more flexible polygonal meshes while retaining convergence guarantees. Numerical experiments in the paper confirm the expected error behavior.

Core claim

Virtual element methods based on the Hodge decomposition of divergence-free vector fields can be designed and analyzed for the quad-curl problem on general planar polygonal domains, with the discrete spaces chosen so that the decomposition is preserved at the algebraic level and the resulting schemes satisfy optimal error estimates.

What carries the argument

Virtual element spaces that discretize the Hodge decomposition of divergence-free fields to enforce the constraint while allowing arbitrary polygonal meshes.

If this is right

  • The methods achieve optimal-order convergence in the natural energy norm for the quad-curl problem.
  • The schemes remain stable on meshes composed of arbitrary convex and non-convex polygons.
  • The discrete Hodge decomposition guarantees that the computed solutions are exactly divergence-free at the algebraic level.
  • The same framework supplies a practical way to treat other higher-order curl problems on non-standard meshes.

Where Pith is reading between the lines

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

  • The same virtual-element construction could be adapted to three-dimensional polyhedral domains for analogous higher-order problems.
  • The approach suggests a route to handle quad-curl equations in electromagnetism or fluid models on domains with re-entrant corners without special mesh refinement.
  • It opens the possibility of combining virtual element spaces with other mimetic discretizations for coupled multiphysics problems.

Load-bearing premise

The Hodge decomposition of divergence-free vector fields admits a stable and accurate discretization inside virtual element spaces on general polygonal domains.

What would settle it

Numerical experiments on a sequence of successively refined polygonal meshes that fail to recover the predicted convergence rates in the appropriate norms or that produce unstable solutions would falsify the analysis.

Figures

Figures reproduced from arXiv: 2604.18406 by Jai Tushar, Li-yeng Sung, Susanne C. Brenner.

Figure 2.1
Figure 2.1. Figure 2.1: Domains with Betti numbers 0, 1 and 2 [PITH_FULL_IMAGE:figures/full_fig_p002_2_1.png] view at source ↗
Figure 6.1
Figure 6.1. Figure 6.1: Structured Voronoi mesh (left) and unstructured Voronoi mesh (right) For the other experiments where the exact solutions are not available, we use nested polytopal meshes to facilitate the computation of the differences between the solutions on consecutive levels. These nested meshes are obtained by the following procedure. (i) We create a random initial Voronoi mesh on the domain Ω. (ii) We connect the … view at source ↗
Figure 6.2
Figure 6.2. Figure 6.2: Refinement procedure: refining a polygon into quadrilaterals (left) and uniform refinement of a quadrilateral (right). Experiment 1 We solve (1.1) on the unit square (0, 1)2 with β = γ = 0, and we use both structured and nonstrutured Voronoi meshes. The manufactured solution is given by u = curl ϕ, where ϕ(x) = sin3 (πx1) sin3 (πx2). The numerical results for the virtual element methods are presented in … view at source ↗
Figure 6.3
Figure 6.3. Figure 6.3: Nested meshes on Ω = (0, 1)2 . Since the exact solution is not available, we displayed in [PITH_FULL_IMAGE:figures/full_fig_p022_6_3.png] view at source ↗
Figure 6.4
Figure 6.4. Figure 6.4: Nested meshes on Ω = (−1, 1)2 \ [PITH_FULL_IMAGE:figures/full_fig_p023_6_4.png] view at source ↗
Figure 6.5
Figure 6.5. Figure 6.5: Nested meshes on Ω = (0, 1)2 \ [1/4, 3/4]2 [PITH_FULL_IMAGE:figures/full_fig_p024_6_5.png] view at source ↗
Figure 6.6
Figure 6.6. Figure 6.6: Nested meshes on Ω = (0, 1)2 \ [PITH_FULL_IMAGE:figures/full_fig_p025_6_6.png] view at source ↗
read the original abstract

We design and analyze virtual element methods for a quad-curl problem on general polygonal domains that are based on the Hodge decomposition of divergence-free vector fields. Numerical results that corroborate the theoretical analysis are also presented.

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

Summary. The manuscript claims to design and analyze virtual element methods (VEM) for a quad-curl problem on general polygonal domains based on the Hodge decomposition of divergence-free vector fields. It supplies explicit local space definitions, projection operators, a discrete de Rham sequence commuting with the continuous one, derives a priori error estimates via adapted VEM arguments, and presents numerical tests on convex and non-convex meshes that confirm optimal convergence.

Significance. This is a solid contribution to numerical PDE methods, particularly for extending VEM to fourth-order problems on general domains. The Hodge decomposition enables effective discretization of the div-free constraint while preserving stability and approximation properties. Credit is due for the explicit constructions and the commuting diagram, which facilitate the error analysis, as well as the numerical corroboration of the theory. Such work is useful for applications involving higher-order curl operators.

minor comments (3)
  1. [Abstract] The abstract is concise but could benefit from including the specific error estimates or convergence orders to give readers a quicker sense of the results.
  2. [Section 3] The definition of the local virtual element spaces would be clearer with an explicit example for a sample polygon, such as a quadrilateral or pentagon.
  3. [Numerical experiments] It would be helpful to include a table listing the observed L2 and H1 error norms and rates for different degrees of approximation.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our work and the recommendation of minor revision. The summary accurately captures the main contributions of the manuscript regarding the design, analysis, and numerical validation of virtual element methods for the quad-curl problem based on Hodge decomposition.

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper constructs VEM spaces for the quad-curl problem by discretizing the standard Hodge decomposition of divergence-free fields on polygonal domains. It supplies explicit local space definitions, commuting discrete de Rham sequences, projection operators, and stability estimates that are independent of the target error bounds. The a priori analysis then follows from consistency arguments adapted to the fourth-order operator, without any reduction of the central claims to self-definitions, fitted parameters renamed as predictions, or load-bearing self-citations that presuppose the result. The approach remains self-contained against external benchmarks of VEM theory and Hodge decomposition.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

No free parameters, axioms, or invented entities are described in the abstract; the work appears to build on existing Hodge decomposition and virtual element frameworks without introducing new postulates.

pith-pipeline@v0.9.0 · 5318 in / 1038 out tokens · 33100 ms · 2026-05-10T03:41:31.496828+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

59 extracted references · 59 canonical work pages

  1. [1]

    Adams and J.J.F

    R.A. Adams and J.J.F. Fournier. Sobolev Spaces (Second Edition). Academic Press, Amsterdam, 2003

    work page 2003

  2. [2]

    Ahmad, A

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

    work page 2013

  3. [3]

    Assous and M

    F. Assous and M. Michaeli. Hodge decomposition to solve singular static Maxwell’s equations in a non-convex polygon. Appl. Numer. Math. , 60:432–441, 2010

    work page 2010

  4. [4]

    Beir˜ ao da Veiga, F

    L. Beir˜ ao da Veiga, F. Brezzi, L.D. Marini, and A. Russo. H(div) and H(curl)-conforming virtual element methods. Numer. Math., 133:303–332, 2016

    work page 2016

  5. [5]

    Beir˜ ao da Veiga, F

    L. Beir˜ ao da Veiga, F. Brezzi, L.D. Marini, and A. Russo. The virtual element method. Acta Numer., 32:123–202, 2023

    work page 2023

  6. [6]

    Beir˜ ao da Veiga, C

    L. Beir˜ ao da Veiga, C. Lovadina, and A. Russo. Stability analysis for the virtual element method.Math. Models Methods Appl. Sci. , 27:2557–2594, 2017

    work page 2017

  7. [7]

    Beir˜ ao da Veiga, F

    L. Beir˜ ao da Veiga, F. Brezzi, A. Cangiani, G. Manzini, L.D. Marini, and A. Russo. Basic principles of virtual element methods. Math. Models Methods Appl. Sci. , 23:199–214, 2013

    work page 2013

  8. [8]

    Bernardi, V

    C. Bernardi, V. Girault, F. Hecht, P.-A. Raviart, and B. Rivi` ere. Mathematics and Finite Element Discretizations of Incompressible Navier-Stokes Flows. Society for Industrial and Applied Mathematics, Philadelphia, 2025

    work page 2025

  9. [9]

    D. Biskamp. Magnetic Reconnection in Plasmas. Cambridge University Press, Cambridge, 2000

    work page 2000

  10. [10]

    Bramble and S.R

    J.H. Bramble and S.R. Hilbert. Estimation of linear functionals on Sobolev spaces with applications to Fourier transforms and spline interpolation. SIAM J. Numer. Anal. , 7:113–124, 1970

    work page 1970

  11. [11]

    Brenner, C

    S.C. Brenner, C. Cavanaugh, and L.-Y. Sung. A Hodge decomposition finite element method for the quad-curl problem on polyhedral domains. J. Sci. Comp. , 100:80, 2024

    work page 2024

  12. [12]

    Brenner, C

    S.C. Brenner, C. Cavanaugh, and L.-Y. Sung. A Hodge decomposition finite element method for an elliptic Maxwell boundary value problem on general polyhedral domains.Comput. Methods Appl. Math., 25:823–848, 2025

    work page 2025

  13. [13]

    Brenner, J

    S.C. Brenner, J. Cui, Z. Nan, and L.-Y. Sung. Hodge decomposition for divergence-free vector fields and two-dimensional Maxwell’s equations. Math. Comp., 81:643–659, 2012

    work page 2012

  14. [14]

    Brenner, J

    S.C. Brenner, J. Cui, and L.-Y. Sung. Multigrid methods based on Hodge decomposition for a quad-curl problem. Comput. Methods Appl. Math. , 19:215–232, 2019

    work page 2019

  15. [15]

    Brenner, J

    S.C. Brenner, J. Gedicke, and L.-Y. Sung. Hodge decomposition for two-dimensional time-harmonic Maxwell’s equation: impedance boundary condition. Math. Methods Appl. Sci. , 40:370–390, 2017

    work page 2017

  16. [16]

    Brenner, Q

    S.C. Brenner, Q. Guan, and L.-Y. Sung. Some estimates for virtual element methods. Comput. Methods Appl. Math., 17:553–574, 2017

    work page 2017

  17. [17]

    Brenner and L.R

    S.C. Brenner and L.R. Scott. The Mathematical Theory of Finite Element Methods (Third Edition ). Springer-Verlag, New York, 2008

    work page 2008

  18. [18]

    Brenner, J

    S.C. Brenner, J. Sun, and L.-Y. Sung. Hodge decomposition methods for a quad-curl problem on planar domains. J. Sci. Comput. , 73:495–513, 2017

    work page 2017

  19. [19]

    Brenner and L.-Y

    S.C. Brenner and L.-Y. Sung. Virtual element methods on meshes with small edges or faces. Math. Models Methods Appl. Sci. , 28(7):1291–1336, 2018

    work page 2018

  20. [20]

    Cakoni, D

    F. Cakoni, D. Colton, P. Monk, and J. Sun. The inverse electromagnetic scattering problem for anisotropic media. Inverse Problems, 26:074004, 2010

    work page 2010

  21. [21]

    Cakoni and H

    F. Cakoni and H. Haddar. A variational approach for the solution of the electromagnetic interior trans- mission problem for anisotropic media. Inverse Probl. Imaging , 1:443–456, 2007

    work page 2007

  22. [22]

    Cangiani, Z

    A. Cangiani, Z. Dong, E.H. Georgoulis, and P. Houston. hp-Version Discontinuous Galerkin Methods on Polygonal and Polyhedral Meshes . Springer, Cham, 2017

    work page 2017

  23. [23]

    S. Cao, L. Chen, and X. Huang. Error analysis of a decoupled finite element method for quad-curl problems. J. Sci. Comput. , 90:Paper No. 29, 25, 2022

    work page 2022

  24. [24]

    Chac´ on, A.N

    L. Chac´ on, A.N. Simakov, and A. Zocco. Steady-state properties of driven magnetic reconnection in 2D electron magnetohydrodynamics. Phys. Rev. Lett., 99:235001, 2007. 30 SUSANNE C. BRENNER, LI-YENG SUNG, AND JAI TUSHAR

    work page 2007

  25. [25]

    G. Chen, W. Qiu, and L. Xu. Analysis of an interior penalty DG method for the quad-curl problem. IMA J. Numer. Anal. , 41:2990–3023, 2021

    work page 2021

  26. [26]

    L. Chen, X. Huang, and C. Zhang. Distributional finite element curl div complexes and application to quad curl problems. SIAM J. Numer. Anal. , 63:1078–1104, 2025

    work page 2025

  27. [27]

    P.G. Ciarlet. The Finite Element Method for Elliptic Problems . North-Holland, Amsterdam, 1978

    work page 1978

  28. [28]

    Dahmen, P

    W. Dahmen, P. Oswald, and X.Q. Shi. C 1-hierarchical bases. J. Comput. Appl. Math. , 51:37–56, 1994

    work page 1994

  29. [29]

    M. Dauge. Elliptic Boundary Value Problems on Corner Domains, Lecture Notes in Mathematics 1341. Springer-Verlag, Berlin-Heidelberg, 1988

    work page 1988

  30. [30]

    Di Pietro and J

    D.A. Di Pietro and J. Droniou. The hybrid high-order method for polytopal meshes . Springer, Cham, 2020

    work page 2020

  31. [31]

    Dupont and R

    T. Dupont and R. Scott. Polynomial approximation of functions in Sobolev spaces. Math. Comp. , 34:441–463, 1980

    work page 1980

  32. [32]

    F¨ uhrer, P

    T. F¨ uhrer, P. Herrera, and N. Heuer. A DPG method for the quad-curl problem.Comput. Math. Appl. , 149:221–238, 2023

    work page 2023

  33. [33]

    Girault and P.-A

    V. Girault and P.-A. Raviart. Finite Element Methods for Navier-Stokes Equations. Theory and Algo- rithms. Springer-Verlag, Berlin, 1986

    work page 1986

  34. [34]

    Grisvard

    P. Grisvard. Elliptic Problems in Non Smooth Domains . Pitman, Boston, 1985

    work page 1985

  35. [35]

    H. Guo, M. Zhang, Q. Zhang, and Z. Zhang. Unfitted finite element method for the quad-curl interface problem. Adv. Comput. Math. , 51:Paper No. 3, 30, 2025

    work page 2025

  36. [36]

    Q. Hong, J. Hu, S. Shu, and J. Xu. A discontinuous Galerkin method for the fourth-order curl problem. J. Comput. Math. , 30:565–578, 2012

    work page 2012

  37. [37]

    K. Hu, Q. Zhang, and Z. Zhang. Simple curl-curl-conforming finite elements in two dimensions. SIAM J. Sci. Comput. , 42:A3859–A3877, 2020

    work page 2020

  38. [38]

    X. Huang. Nonconforming finite element Stokes complexes in three dimensions. Sci. China Math. , 66:1879–1902, 2023

    work page 1902

  39. [39]

    Huang and C

    X. Huang and C. Zhang. Robust mixed finite element methods for a quad-curl singular perturbation problem. J. Comput. Appl. Math. , 451:Paper No. 116117, 19, 2024

    work page 2024

  40. [40]

    J. Li, S. Mao, C. Wang, and Z. Zhang. An efficient finite element method for the quad-curl problem. SIAM J. Numer. Anal. , 63:2272–2295, 2025

    work page 2025

  41. [41]

    P. Monk. Finite Element Methods for Maxwell’s Equations . Oxford University Press, New York, 2003

    work page 2003

  42. [42]

    Monk and J

    P. Monk and J. Sun. Finite element methods for Maxwell’s transmission eigenvalues. SIAM J. Sci. Comput., 34:B247–B264, 2012

    work page 2012

  43. [43]

    L. Mu, J. Wang, and X. Ye. Weak Galerkin finite element methods on polytopal meshes. Int. J. Numer. Anal. Model., 12:31–53, 2015

    work page 2015

  44. [44]

    Neˇ cas.Direct Methods in the Theory of Elliptic Equations

    J. Neˇ cas.Direct Methods in the Theory of Elliptic Equations . Springer, Heidelberg, 2012

    work page 2012

  45. [45]

    J. Sun. A mixed FEM for the quad-curl eigenvalue problem. Numer. Math., 132:185–200, 2016

    work page 2016

  46. [46]

    J. Sun, Q. Zhang, and Z. Zhang. A curl-conforming weak Galerkin method for the quad-curl problem. BIT, 59:1093–1114, 2019

    work page 2019

  47. [47]

    Z. Sun, J. Cui, F. Gao, and C. Wang. Multigrid methods for a quad-curl problem based on C 0 interior penalty method. Comput. Math. Appl. , 76:2192–2211, 2018

    work page 2018

  48. [48]

    Talischi, G.H

    C. Talischi, G.H. Paulino, A. Pereira, and I.F.M. Menezes. PolyMesher: a general-purpose mesh gen- erator for polygonal elements written in Matlab. Struct. Multidiscip. Optim. , 45:309–328, 2012

    work page 2012

  49. [49]

    C. Wang, J. Cui, and Z. Sun. A linear edge finite element method for quad-curl problem. Comput. Math. Appl., 176:99–108, 2024

    work page 2024

  50. [50]

    C. Wang, J. Wang, and S. Zhang. Weak Galerkin finite element methods for quad-curl problems. J. Comput. Appl. Math. , 428:Paper No. 115186, 13, 2023

    work page 2023

  51. [51]

    Wang and S

    C. Wang and S. Zhang. Stabilizer-free weak Galerkin methods for quad-curl problems on polyhedral meshes without convexity assumptions. J. Comput. Appl. Math. , 481:Paper No. 117279, 2026

    work page 2026

  52. [52]

    L. Wang, H. Li, Q. Zhang, and Z. Zhang. H(curl2)-conforming triangular spectral element method for quad-curl problems. J. Comput. Appl. Math. , 459:Paper No. 116362, 21, 2025. VIRTUAL ELEMENT METHODS FOR A PLANAR QUAD-CURL PROBLEM 31

    work page 2025

  53. [53]

    L. Wang, W. Shan, H. Li, and Z Zhang. H(curl 2)-conforming quadrilateral spectral element method for quad-curl problems. Math. Models Methods Appl. Sci. , 31:1951–1986, 2021

    work page 1951

  54. [54]

    L. Wang, M. Zhang, and Q. Zhang. Fully H (gradcurl)-nonconforming finite element method for the singularly perturbed quad-curl problem on cubical meshes. IMA J. Numer. Anal. , 46:271–298, 2026

    work page 2026

  55. [55]

    Zhang and Z

    B. Zhang and Z. Zhang. A new family of nonconforming elements with H(curl)-continuity for the 3D quad-curl problem. Commun. Comput. Phys. , 33:1069–1089, 2023

    work page 2023

  56. [56]

    Zhang, L

    Q. Zhang, L. Wang, and Z. Zhang. H(curl2)-conforming finite elements in 2 dimensions and applications to the quad-curl problem. SIAM J. Sci. Comput. , 41:A1527–A1547, 2019

    work page 2019

  57. [57]

    J. Zhao, K. Deng, C. Wang, and B. Zhang. Error analysis on the mixed finite element method for a quad-curl problem with low-order terms in three dimensions. Appl. Numer. Math. , 222:17–31, 2026

    work page 2026

  58. [58]

    Zhao and B

    J. Zhao and B. Zhang. The curl-curl conforming virtual element method for the quad-curl problem. Math. Models Methods Appl. Sci. , 31:1659–1690, 2021

    work page 2021

  59. [59]

    Zheng, Q

    B. Zheng, Q. Hu, and J. Xu. A nonconforming finite element method for fourth order curl equations in R3. Math. Comp., 80:1871–1886, 2011. Susanne C. Brenner, Department of Mathematics and Center for Computation and Tech- nology, Louisiana State University, Baton Rouge, LA 70803, USA Email address : brenner@math.lsu.edu Li-yeng Sung, Department of Mathemat...

    work page 2011