A Framework for Analysis of DEC Approximations to Hodge-Laplacian Problems using Generalized Whitney Forms
Pith reviewed 2026-05-22 15:12 UTC · model grok-4.3
The pith
Equivalence between DEC cochains and generalized Whitney forms allows FEEC tools to prove convergence rates for Hodge-Laplacian approximations on well-centered meshes.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By establishing that primal and dual cochains coincide with Whitney and generalized Whitney forms, DEC discretizations of the Hodge-Laplacian become amenable to the full machinery of FEEC, which yields rigorous a-priori error bounds with explicit rates that hold uniformly for all k-form degrees on well-centered meshes.
What carries the argument
The equivalence that identifies cochains on primal and dual meshes with Whitney and generalized Whitney forms, thereby transferring FEEC convergence theory to DEC schemes.
If this is right
- Optimal convergence rates hold for the Hodge-Laplacian in full k-form generality on well-centered meshes.
- The observed numerical errors match the theoretically derived rates.
- Superconvergence effects receive a consistent explanation inside the same equivalence framework.
Where Pith is reading between the lines
- The same equivalence argument could be tested on other DEC schemes for different exterior differential operators.
- Extending the analysis beyond well-centered meshes would require either new form definitions or additional correction terms.
Load-bearing premise
The meshes must be well-centered so that the cochains on the primal and dual grids line up exactly with the generalized Whitney forms.
What would settle it
Numerical computation of the Hodge-Laplacian error on a sequence of non-well-centered meshes that shows the predicted convergence rate is not attained.
Figures
read the original abstract
We provide a framework for interpreting Discrete Exterior Calculus (DEC) numerical schemes in terms of Finite Element Exterior Calculus (FEEC). We demonstrate the equivalence of cochains on primal and dual meshes with Whitney and generalized Whitney forms which allows us to analyze DEC approximations using tools from FEEC. We demonstrate the applicability of our framework by rigorously proving convergence with rates for the Hodge-Laplacian problem in full $k$-form generality on well-centered meshes. We also provide numerical results illustrating optimality of our derived convergence rates. Moreover, we demonstrate how superconvergence phenomena can be explained in our framework with corresponding numerical results.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a framework linking Discrete Exterior Calculus (DEC) schemes for Hodge-Laplacian problems to Finite Element Exterior Calculus (FEEC) by establishing an equivalence between cochains on primal/dual meshes and Whitney/generalized Whitney forms. This equivalence is used to transfer FEEC tools and prove rigorous convergence rates for the Hodge-Laplacian in full k-form generality on well-centered meshes. Numerical results are presented to confirm optimal rates and to interpret superconvergence phenomena within the framework.
Significance. If the equivalence preserves the de Rham complex structure and commuting diagram properties without k-dependent degradation, the work provides a valuable bridge between DEC (popular in computational physics) and the mature FEEC theory, enabling rigorous a priori estimates for DEC discretizations that previously lacked them. The full k-generality claim and the explanation of superconvergence via the framework are notable strengths; the well-centered mesh restriction is standard but explicitly acknowledged.
major comments (2)
- [§3.3, Theorem 3.8] §3.3, Theorem 3.8 (equivalence and subcomplex property): The proof that generalized Whitney forms form a subcomplex with the same cohomology as the continuous de Rham complex for arbitrary k relies on the well-centered mesh condition to ensure the interpolation operator commutes with d. However, the argument for k>1 appears to use a recursive construction whose constants may depend on k; this could prevent direct inheritance of the optimal FEEC rates without additional mesh-dependent factors. A explicit bound independent of k is needed to support the central convergence claim.
- [§5, Theorem 5.2] §5, Theorem 5.2 (convergence rates): The a priori error estimate for the Hodge-Laplacian is transferred from FEEC, but the proof sketch does not explicitly verify that the DEC mass matrix and coboundary operator induce the same stability constants as the FEEC inner product for dual meshes when k=2 or 3. If the equivalence only holds up to a mesh-dependent factor, the claimed optimality may not hold uniformly.
minor comments (2)
- [Figures 2 and 4] Figure 2 and 4: The mesh illustrations for well-centered vs. non-well-centered cases would benefit from explicit annotation of the circumcenter locations to clarify the geometric condition used in the equivalence.
- [§3 and §4] Notation: The symbol for the generalized Whitney form (denoted W_k in some places and Ŵ_k in others) should be unified throughout §3 and §4 to avoid confusion.
Simulated Author's Rebuttal
We thank the referee for their thorough review and valuable suggestions. We have carefully considered the major comments and provide point-by-point responses below. Where appropriate, we have revised the manuscript to address the concerns raised.
read point-by-point responses
-
Referee: [§3.3, Theorem 3.8] The proof that generalized Whitney forms form a subcomplex with the same cohomology as the continuous de Rham complex for arbitrary k relies on the well-centered mesh condition to ensure the interpolation operator commutes with d. However, the argument for k>1 appears to use a recursive construction whose constants may depend on k; this could prevent direct inheritance of the optimal FEEC rates without additional mesh-dependent factors. A explicit bound independent of k is needed to support the central convergence claim.
Authors: We appreciate this observation. The recursive construction for higher k in the proof of Theorem 3.8 is based on the definition of generalized Whitney forms, which are constructed inductively. However, because the mesh is well-centered, the circumcenters are inside the simplices, ensuring that the interpolation operator I commutes with the exterior derivative d independently of k. The constants in the estimates arise from the boundedness of the mesh angles and the volume ratios, which are controlled by the well-centeredness condition and do not depend on the form degree k. To make this explicit, we have added a new remark following Theorem 3.8 that provides a k-independent bound on the relevant operator norms, derived from the uniform bound on the mesh quality. This allows direct transfer of the FEEC convergence rates without k-dependent degradation. revision: yes
-
Referee: [§5, Theorem 5.2] The a priori error estimate for the Hodge-Laplacian is transferred from FEEC, but the proof sketch does not explicitly verify that the DEC mass matrix and coboundary operator induce the same stability constants as the FEEC inner product for dual meshes when k=2 or 3. If the equivalence only holds up to a mesh-dependent factor, the claimed optimality may not hold uniformly.
Authors: We agree that explicit verification is important for clarity. The equivalence established in Section 3 shows that the DEC cochains with the primal-dual inner product correspond exactly to the generalized Whitney forms with the L2 inner product, and the coboundary operator matches the exterior derivative. For dual meshes, the mass matrix is defined via the dual volumes, which under the well-centered assumption, preserves the stability constants uniformly. We have expanded the proof of Theorem 5.2 to include a direct comparison of the stability constants for k=2 and k=3, showing they are identical to those in the FEEC setting up to a factor depending only on the mesh regularity (which is fixed). This confirms the optimality of the rates uniformly in k. revision: yes
Circularity Check
No significant circularity; equivalence and FEEC transfer are independently derived.
full rationale
The paper derives the equivalence of cochains to generalized Whitney forms on well-centered meshes within its own analysis, then applies standard FEEC a priori estimates to obtain convergence rates for the Hodge-Laplacian. No load-bearing step reduces by construction to a fitted parameter, self-citation chain, or renamed input; the central claim rests on a self-contained proof of the commuting diagram property and subcomplex structure rather than external self-referential assumptions. This is the typical non-circular outcome for a framework paper that proves its key equivalence explicitly.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Well-centered meshes allow equivalence of cochains on primal and dual meshes with Whitney and generalized Whitney forms.
Forward citations
Cited by 1 Pith paper
-
Convergence of Discrete Exterior Calculus for the Hodge-Dirac Operator
Supplies a convergence proof for DEC discretization of the Hodge-Dirac operator by adapting analysis techniques from a cited paper on Hodge-Laplacian problems.
Reference graph
Works this paper leans on
-
[1]
LocalL2-bounded commuting projections in FEEC.ESAIM Math
Arnold, D., and Guzm´ an, J. LocalL2-bounded commuting projections in FEEC.ESAIM Math. Model. Numer. Anal.55(2021), 2169–2184.https://doi.org/10.1051/m2an/2021054
-
[2]
Arnold, D. N., Falk, R. S., and Winther, R. Finite element exterior calculus, homological techniques, and applications.Acta Numerica15(2006), 1–155
work page 2006
-
[3]
Arnold, D. N., Falk, R. S., and Winther, R. Finite element exterior calculus: from Hodge theory to numerical sta- bility.Bull. Amer. Math. Soc. (N.S.)47(2010), 281–354.https://doi.org/10.1090/S0273-0979-10-01278-4
-
[4]
A note on the Poincar´ e inequality for convex domains.Z
Bebendorf, M. A note on the Poincar´ e inequality for convex domains.Z. Anal. Anwendungen22(2003), 751–756. https://doi.org/10.4171/ZAA/1170. 30 JOHNNY GUZM ´AN, PRATYUSH POTU
-
[5]
Beir˜ ao da Veiga, L. c., Brezzi, F., Marini, L. D., and Russo, A. The virtual element method.Acta Numer.32 (2023), 123–202.https://doi.org/10.1017/S0962492922000095
-
[6]
Bochev, P. B., and Hyman, J. M. Principles of mimetic discretizations of differential operators. InCompatible spatial discretizations, vol.142ofIMA Vol. Math. Appl.Springer, New York, 2006, pp. 89–119.https:// doi-org.revproxy.brown.edu/10.1007/0-387-38034-5_5
-
[7]
PhD thesis, University of Paris-Est, 2014
Bonelle, J.Compatible discrete operator schemes on polyhedral meshes for elliptic and Stokes equations. PhD thesis, University of Paris-Est, 2014
work page 2014
-
[8]
Christiansen, S. H. A construction of spaces of compatible differential forms on cellular complexes.Math. Models Methods Appl. Sci.18(2008), 739–757.https://doi.org/10.1142/S021820250800284X
-
[9]
Christiansen, S. H. Foundations of finite element methods for wave equations of Maxwell type. InApplied wave mathematics. Springer, Berlin, 2009, pp. 335–393.https://doi.org/10.1007/978-3-642-00585-5_17
-
[10]
Topics in structure-preserving discretiza- tion
Christiansen, S. H., Munthe-Kaas, H. Z., and Owren, B. Topics in structure-preserving discretization.Acta Numer.20(2011), 1–119.https://doi.org/10.1017/S096249291100002X
-
[11]
Christiansen, S. H., and Winther, R. Smoothed projections in finite element exterior calculus.Math. Comp.77 (2008), 813–829.https://doi.org/10.1090/S0025-5718-07-02081-9
-
[12]
A 3D discrete duality finite volume method for nonlinear elliptic equations.SIAM J
Coudi` ere, Y., and Hubert, F. A 3D discrete duality finite volume method for nonlinear elliptic equations.SIAM J. Sci. Comput.33(2011), 1739–1764.https://doi.org/10.1137/100786046
-
[13]
Digital geometry processing with discrete exterior calculus
Crane, K., de Goes, F., Desbrun, M., and Schr¨ oder, P. Digital geometry processing with discrete exterior calculus. InACM SIGGRAPH 2013 Courses(New York, NY, USA, 2013), SIGGRAPH ’13, Association for Computing Machinery.https://doi.org/10.1145/2504435.2504442
-
[14]
A glimpse into discrete differential geometry.Notices Amer
Crane, K., and Wardetzky, M. A glimpse into discrete differential geometry.Notices Amer. Math. Soc.64(2017), 1153–1159.https://doi.org/10.1090/noti1578
-
[15]
Discrete exterior calculus for variational problems in computer vision and graphics
Desbrun, M., Hirani, A., and Marsden, J. Discrete exterior calculus for variational problems in computer vision and graphics. In42nd IEEE International Conference on Decision and Control (IEEE Cat. No.03CH37475) (2003), vol.5, pp. 4902–4907 Vol.5
work page 2003
-
[16]
Desbrun, M., Hirani, A. N., Leok, M., and Marsden, J. E. Discrete exterior calculus. arXiv preprint.,https: //arxiv.org/abs/math/0508341, 2005
work page internal anchor Pith review Pith/arXiv arXiv 2005
-
[17]
Droniou, J., Eymard, R., Gallou¨ et, T., Guichard, C., and Herbin, R.The gradient discretisation method, vol.82ofMath´ ematiques & Applications (Berlin) [Mathematics & Applications]. Springer, Cham, 2018.https: //doi.org/10.1007/978-3-319-79042-8
-
[18]
An interpretation of trisk-type schemes from a discrete exterior calculus perspective
Eldred, C., and Bauer, W. An interpretation of trisk-type schemes from a discrete exterior calculus perspective. arXiv preprint.,https://arxiv.org/abs/2210.07476, 2022
-
[19]
Ern, A., and Guermond, J.-L.Theory and practice of finite elements, vol.159ofApplied Mathematical Sciences. Springer-Verlag, New York, 2004.https://doi.org/10.1007/978-1-4757-4355-5
-
[20]
Springer, Cham, [2021]©2021.https://doi.org/10.1007/978-3-030-56341-7
Ern, A., and Guermond, J.-L.Finite elements I—Approximation and interpolation, vol.72ofTexts in Applied Mathematics. Springer, Cham, [2021]©2021.https://doi.org/10.1007/978-3-030-56341-7
-
[21]
Local l2-bounded commuting projections using discrete local problems on alfeld splits
Ern, A., Guzman, J., Potu, P., and Vohralik, M. Local l2-bounded commuting projections using discrete local problems on alfeld splits. arXiv preprint.,https://arxiv.org/abs/2502.03601, 2025
-
[22]
Eymard, R., Gallou¨ et, T., and Herbin, R. Finite volume methods. InHandbook of numerical analysis, Vol. VII, vol.VIIofHandb. Numer. Anal.North-Holland, Amsterdam, 2000, pp. 713–1020.https://doi.org/10.1016/ S1570-8659(00)07005-8
work page 2000
-
[23]
Falk, R. S., and Winther, R. Local bounded cochain projections.Math. Comp.83(2014), 2631–2656.https: //doi.org/10.1090/S0025-5718-2014-02827-5
-
[24]
PhD thesis, University of Texas at Austin, 2011
Gillette, A.Stability of Dual Discretization Methods for Partial Differential Equations. PhD thesis, University of Texas at Austin, 2011
work page 2011
-
[25]
Upwind schemes for scalar advection-dominated problems in the discrete exterior calculus
Griebel, M., Rieger, C., and Schier, A. Upwind schemes for scalar advection-dominated problems in the discrete exterior calculus. InTransport processes at fluidic interfaces, Adv. Math. Fluid Mech. Birkh¨ auser/Springer, Cham, 2017, pp. 145–175
work page 2017
-
[26]
On the convergence of metric and geometric properties of polyhedral surfaces.Geom
Hildebrandt, K., Polthier, K., and Wardetzky, M. On the convergence of metric and geometric properties of polyhedral surfaces.Geom. Dedicata123(2006), 89–112.https://doi.org/10.1007/s10711-006-9109-5
-
[27]
Hirani, A. N. Thesis (Ph.D.)–California Institute of Technology.Discrete exterior calculus. ProQuest LLC, Ann Arbor, MI, 2003.http://gateway.proquest.com/openurl?url_ver=Z39.88-2004&rft_val_fmt=info:ofi/fmt: kev:mtx:dissertation&res_dat=xri:pqdiss&rft_dat=xri:pqdiss:3086864. A FRAMEWORK FOR ANALYSIS OF DEC APPROXIMATIONS 31
work page 2003
- [28]
-
[29]
Mimetic framework on curvilinear quadrilaterals of arbitrary order
Kreeft, J., Palha, A., and Gerritsma, M. Mimetic framework on curvilinear quadrilaterals of arbitrary order. arXiv preprint.,https://arxiv.org/abs/1111.4304, 2011
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[30]
Mohamed, M. S., Hirani, A. N., and Samtaney, R. Discrete exterior calculus discretization of incompressible Navier-Stokes equations over surface simplicial meshes.J. Comput. Phys.312(2016), 175–191.https://doi. org/10.1016/j.jcp.2016.02.028
-
[31]
Mohamed, M. S., Hirani, A. N., and Samtaney, R. Numerical convergence of discrete exterior calculus on arbitrary surface meshes.Int. J. Comput. Methods Eng. Sci. Mech.19(2018), 194–206.https://doi.org/10. 1080/15502287.2018.1446196
-
[32]
Morris, L., Baas, A., Arias, J., Gatlin, M., Patterson, E., and Fairbanks, J. P. Decapodes: A diagrammatic tool for representing, composing, and computing spatialized partial differential equations.Journal of Computational Science81(2024), 102345.https://www.sciencedirect.com/science/article/pii/S1877750324001388
work page 2024
-
[33]
R.Elements of algebraic topology
Munkres, J. R.Elements of algebraic topology. Addison-Wesley Publishing Company, Menlo Park, CA, 1984
work page 1984
-
[34]
Discrete exterior calculus (DEC) for the surface Navier-Stokes equation
Nitschke, I., Reuther, S., and Voigt, A. Discrete exterior calculus (DEC) for the surface Navier-Stokes equation. InTransport processes at fluidic interfaces, Adv. Math. Fluid Mech. Birkh¨ auser/Springer, Cham, 2017, pp. 177– 197
work page 2017
-
[35]
Pavlov, D., Mullen, P., Tong, Y., Kanso, E., Marsden, J. E., and Desbrun, M. Structure-preserving discretization of incompressible fluids.Phys. D240(2011), 443–458.https://doi.org/10.1016/j.physd.2010.10.012
-
[36]
Payne, L. E., and Weinberger, H. F. An optimal Poincar´ e inequality for convex domains.Arch. Rational Mech. Anal.5(1960), 286–292.https://doi.org/10.1007/BF00252910
-
[37]
Schulz, E., and Tsogtgerel, G. Convergence of discrete exterior calculus approximations for Poisson problems. Discrete Comput. Geom.63(2020), 346–376.https://doi.org/10.1007/s00454-019-00159-x
-
[38]
Stern, A., Tong, Y., Desbrun, M., and Marsden, J. E. Geometric computational electrodynamics with variational integrators and discrete differential forms. InGeometry, mechanics, and dynamics, vol.73ofFields Inst. Commun.Springer, New York, 2015, pp. 437–475.https://doi.org/10.1007/978-1-4939-2441-7_19
-
[39]
Vanderzee, E., Hirani, A. N., Guoy, D., and Ramos, E. A. Well-centered triangulation.SIAM J. Sci. Comput. 31(2009/10), 4497–4523.https://doi.org/10.1137/090748214
-
[40]
Wang, M., Jagad, P., Hirani, A. N., and Samtaney, R. Discrete exterior calculus discretization of two-phase incompressible Navier-Stokes equations with a conservative phase field method.J. Comput. Phys.488(2023), Paper No. 112245, 29.https://doi.org/10.1016/j.jcp.2023.112245
-
[41]
Princeton University Press, Princeton, NJ, 1957
Whitney, H.Geometric integration theory. Princeton University Press, Princeton, NJ, 1957
work page 1957
-
[42]
Yavari, On geometric discretization of elasticity, Journal of Mathematical Physics 49 (2008) 022901
Yavari, A. On geometric discretization of elasticity.J. Math. Phys.49(2008), 022901, 36.https://doi.org/ 10.1063/1.2830977
-
[43]
Zhu, C., Christiansen, S. H., Hu, K., and Hirani, A. N. Convergence and stability of discrete exterior calculus for the hodge laplace problem in two dimensions. arXiv preprint.,https://arxiv.org/abs/2505.08966, 2025. AppendixA.Construction of Generalized Whitney forms We follow the work of Christiansen ([8, 9] and [10, Section 5]) to define finite dimensi...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.