A Mixed Finite Element Method for the Dirichlet Vector Laplacian in Three Dimensions

Peiyang Yu; Ralf Hiptmair; Tianwei Yu

arxiv: 2603.14026 · v3 · pith:T322CL2Vnew · submitted 2026-03-14 · 🧮 math.NA · cs.NA

A Mixed Finite Element Method for the Dirichlet Vector Laplacian in Three Dimensions

Ralf Hiptmair , Peiyang Yu , Tianwei Yu This is my paper

Pith reviewed 2026-05-15 11:20 UTC · model grok-4.3

classification 🧮 math.NA cs.NA MSC 65N3065N12

keywords mixed finite elementsvector LaplacianDirichlet boundary conditionsde Rham complexa priori error estimatesCaccioppoli inequalityStokes discretization

0 comments

The pith

Mixed finite element method establishes well-posedness for the three-dimensional Dirichlet vector Laplacian via a non-standard vorticity space.

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

This paper develops and analyzes a mixed finite element method for the vector Laplacian boundary value problem in three dimensions subject to Dirichlet conditions. The Dirichlet boundary data breaks the usual de Rham complex structure, so the vorticity is sought in a non-standard function space that restores well-posedness of the discrete problem. Error estimates are derived by means of a discrete Caccioppoli-type inequality for curl-harmonic functions, yielding (k-1/2)-order convergence in the energy norm on general domains and k-order convergence in L2 on convex domains. The analysis extends an earlier two-dimensional treatment to three-dimensional domains of arbitrary topology and supplies a direct discretization of the Stokes equations in vorticity-velocity-pressure form.

Core claim

The mixed FEEC-type finite element approximation of the three-dimensional vector Laplace problem with Dirichlet boundary conditions is well-posed, and the method admits a priori error bounds of order (k-1/2) in the energy norm on general domains and order k in L2 on convex domains, where k denotes the polynomial degree of the finite element spaces.

What carries the argument

A discrete Caccioppoli-type inequality for discrete curl-harmonic functions in the chosen finite element spaces, which controls the error despite the disruption of the de Rham complex by the Dirichlet condition.

If this is right

The same spaces yield a convergent discretization of the Stokes problem in vorticity-velocity-pressure form.
Optimal L2 convergence is recovered when the domain is convex.
The method extends the two-dimensional mixed finite element analysis of Arnold, Falk and Gopalakrishnan to three dimensions with general topology.
Suboptimal rates are intrinsic on non-convex domains because of the boundary condition.

Where Pith is reading between the lines

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

Similar non-standard spaces may be required for other boundary-value problems that break the de Rham sequence in higher dimensions.
Adaptive refinement near the boundary could potentially restore optimal rates on general domains.
The well-posedness framework might extend directly to time-dependent or nonlinear vector Laplace problems.

Load-bearing premise

The discrete Caccioppoli-type inequality for discrete curl-harmonic functions holds in the chosen finite element spaces on general three-dimensional domains.

What would settle it

Numerical computation on a non-convex domain showing either loss of discrete well-posedness or convergence rates strictly worse than (k-1/2) in the energy norm.

read the original abstract

This work establishes the well-posedness and a priori error analysis for the mixed FEEC-type finite element approximation of the three-dimensional vector Laplace boundary value problem subject to the Dirichlet boundary condition. The Dirichlet condition disrupts the structure of the standard de Rham complex, requiring the vorticity to be sought in a non-standard function space to achieve well-posedness. We derive error estimates that confirm the numerically observed suboptimal convergence rates. In particular, by developing a discrete Caccioppoli-type inequality for discrete curl-harmonic functions, we prove $(k-1/2)$-th order convergence in the energy norm on general domains and $k$-th order convergence in $L^2$ on convex domains, where $k \ge 1$ denotes polynomial degree of the finite element spaces. These results extend the previous two-dimensional analysis developed in [Arnold, D.N., Falk, R.S. and Gopalakrishnan, J., 2012. Mixed finite element approximation of the vector Laplacian with Dirichlet boundary conditions. Mathematical Models and Methods in Applied Sciences, 22(9), p.1250024.]~to three-dimensional domains with general topology. As a direct application, a discretization of the Stokes problem in vorticity-velocity-pressure form is studied.

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

Extends the 2D mixed analysis for the Dirichlet vector Laplacian to 3D with a new discrete Caccioppoli inequality that handles general topology and delivers the observed rates.

read the letter

This paper extends the 2012 Arnold-Falk-Gopalakrishnan 2D mixed finite element analysis for the vector Laplacian with Dirichlet boundary conditions to three dimensions on domains with general topology. The main new piece is a discrete Caccioppoli-type inequality for discrete curl-harmonic functions in the FEEC spaces, which lets them prove well-posedness and the a priori estimates despite the non-standard vorticity space forced by the Dirichlet condition. They recover the (k-1/2) energy-norm rate on general domains and the k rate in L2 on convex domains, and they note that these match the numerics. The Stokes application in vorticity-velocity-pressure form is a direct and useful consequence once the vector Laplacian piece is settled. The setup stays grounded in standard de Rham complex properties plus the new inequality, with no obvious circularity or fitted parameters. The soft spot is the inequality itself: in 3D the space of harmonic fields grows with nontrivial topology, so the proof must control the discrete curl-harmonic component without losing the rate or adding hidden mesh restrictions. The abstract states that it works, but the full derivation needs checking to confirm there are no gaps when the cohomology is nontrivial. This is for numerical analysts who need stable mixed methods for vector problems or Stokes-type systems. A reader who already knows the 2D theory will see a substantive but not revolutionary step forward. It deserves peer review because the claims are concrete, the extension fills a documented gap, and the rates are tied to verifiable numerics.

Referee Report

2 major / 2 minor

Summary. The paper develops a mixed finite-element exterior calculus (FEEC) discretization for the three-dimensional vector Laplacian subject to Dirichlet boundary conditions. It establishes well-posedness by placing the vorticity in a non-standard space that restores the required de Rham structure, derives a discrete Caccioppoli-type inequality for curl-harmonic fields, and obtains a priori error estimates of order k-1/2 in the energy norm on general domains and order k in L² on convex domains (k ≥ 1). The analysis extends the 2012 Arnold-Falk-Gopalakrishnan 2D result to three dimensions with arbitrary topology and includes an application to the Stokes problem in vorticity-velocity-pressure form.

Significance. If the central technical step is rigorous, the work supplies the first complete well-posedness and error theory for mixed FEEC approximations of the Dirichlet vector Laplacian in 3D. The explicit handling of non-trivial cohomology and the confirmation of the observed suboptimal rates are valuable for applications in incompressible flow and electromagnetism; the extension from the 2D reference is a clear technical advance.

major comments (2)

[§3] §3 (discrete Caccioppoli inequality): the proof that the inequality holds for discrete curl-harmonic functions in the chosen spaces on domains with non-trivial topology is load-bearing for the (k-1/2) energy-norm rate. The argument must explicitly control the harmonic component arising from the Dirichlet condition and the enlarged cohomology; any gap here would invalidate the claimed rates on general domains.
[Theorem 4.2] Theorem 4.2 (energy-norm estimate): the constant in the (k-1/2) bound appears to depend on the discrete harmonic projection; please verify that this dependence remains uniform with respect to mesh size and does not degrade the rate when the first Betti number is positive.

minor comments (2)

[Abstract] The abstract states that the rates 'confirm the numerically observed suboptimal convergence rates' but does not cite the section containing the numerical experiments; add a parenthetical reference.
[§2] Notation for the non-standard vorticity space (introduced after Eq. (2.7)) should be used consistently in all subsequent statements of the discrete problem and error estimates.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below with clarifications and proposed revisions to strengthen the presentation.

read point-by-point responses

Referee: [§3] §3 (discrete Caccioppoli inequality): the proof that the inequality holds for discrete curl-harmonic functions in the chosen spaces on domains with non-trivial topology is load-bearing for the (k-1/2) energy-norm rate. The argument must explicitly control the harmonic component arising from the Dirichlet condition and the enlarged cohomology; any gap here would invalidate the claimed rates on general domains.

Authors: We appreciate the referee highlighting this critical step. In the proof of the discrete Caccioppoli inequality in Section 3, the harmonic component is controlled explicitly via the discrete Hodge decomposition in the non-standard vorticity space, which incorporates the Dirichlet boundary conditions and the enlarged cohomology to handle non-trivial topology. The bound on the harmonic fields follows from the finite-dimensionality of the cohomology groups and the stability of the discrete de Rham complex, yielding a mesh-independent constant. To make this control fully transparent, we will revise the manuscript by adding a dedicated remark and expanded steps in the proof of the inequality. revision: partial
Referee: [Theorem 4.2] Theorem 4.2 (energy-norm estimate): the constant in the (k-1/2) bound appears to depend on the discrete harmonic projection; please verify that this dependence remains uniform with respect to mesh size and does not degrade the rate when the first Betti number is positive.

Authors: In Theorem 4.2 the constant depends on the discrete harmonic projection, but Lemmas 3.1–3.3 establish that this projection is stable with a bound independent of the mesh size h, including when the first Betti number is positive. The uniformity follows from the finite-dimensional cohomology and the properties of the FEEC spaces under the discrete Hodge decomposition. Consequently the (k-1/2) rate is preserved. We will insert an explicit remark in the statement and proof of Theorem 4.2 confirming this mesh-independent bound. revision: partial

Circularity Check

0 steps flagged

No significant circularity; discrete Caccioppoli inequality derived independently

full rationale

The paper's core contribution is a new derivation of the discrete Caccioppoli-type inequality for curl-harmonic functions in 3D mixed FEEC spaces under Dirichlet conditions, used to obtain the stated convergence rates. This step is presented as an original analysis extending the 2D Arnold-Falk-Gopalakrishnan result (different authors, 2012) rather than reducing to it by construction, self-definition, or load-bearing self-citation. No parameters are fitted and then renamed as predictions; the well-posedness and error estimates follow from standard de Rham complex properties plus the newly proved inequality on general domains. The derivation chain remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The analysis relies on the standard de Rham complex in 3D and on the existence of finite element spaces that preserve discrete exact sequences; the key new ingredient is the discrete Caccioppoli inequality developed inside the paper.

axioms (2)

standard math The de Rham complex structure holds in three dimensions and is disrupted by Dirichlet boundary conditions
Invoked to motivate the non-standard space for vorticity.
domain assumption Suitable finite element spaces exist that satisfy discrete de Rham properties
Required for the mixed FEEC formulation to be well-defined.

pith-pipeline@v0.9.0 · 5526 in / 1314 out tokens · 61545 ms · 2026-05-15T11:20:10.376437+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking echoes

?

echoes
ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

extension to three-dimensional domains with general topology... discrete Caccioppoli-type inequality for discrete curl-harmonic functions... H1, H2 harmonic spaces characterized by Betti numbers b1, b2

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.