{ "schemaVersion": "1.6", "slug": "tetrahedral-moment-stokes-elements", "title": "A moment obstruction and minimal-degree polynomial nonconforming Stokes elements on tetrahedra", "shortTitle": "Why fourth-order tetrahedral Stokes elements need a sextic", "url": "https://evidencepress.org/releases/tetrahedral-moment-stokes-elements/", "oneLine": "Answers the open question of Balazi, Allaire, Jolivet and Omnes: no fifteen quintics complete a fourth-order moment-based Stokes element on tetrahedra, because the 50 moments satisfy a hidden relation on all quintics; one sextic repairs it, and for every order degree n+2 misses exactly (n−1)(n−2)/2 directions while degree n+3 always works.", "abstract": "Nonconforming finite elements for slow viscous (Stokes) flow are often specified by moments: averages of the velocity against polynomials on each face and inside each cell. Balazi, Allaire, Jolivet and Omnes (Numerical Methods for PDEs, 2025) built tetrahedral elements of second and third order this way and left open whether a fourth-order element can be obtained by adding fifteen quintic polynomials to the quartics. This paper shows that it cannot: an alternating face functional vanishes on every quintic but equals 1/34650 on the Vandermonde sextic of the barycentric coordinates, so the 50 moments are dependent on all quintics, and that one sextic repairs the defect, giving a 50-dimensional element independent of vertex numbering. For face moments of degree n and cell moments of degree n−1, the polynomials of degree n+2 invisible to all moments form a seven-dimensional space with an explicit basis of Jacobi and Legendre ridges for every n. Hence degree n+2 suffices only for n=1,2; for n≥3 exactly (n−1)(n−2)/2 moment directions are missing and degree n+3 is necessary and sufficient, for single-cell polynomial spaces containing all polynomials of degree n+1. Numbering-free elements with the fewest top-degree functions are proved to exist for every n, given explicitly for n=4,5 and checked for an explicit family up to n=8. On triangles no extra degree is ever needed; on higher-dimensional simplices sufficiency of degree n+3 is conjectured beyond checked cases. A numerical Stokes test shows rates consistent with fourth-order convergence, with computed discrete inf-sup values 0.40–0.37 on four meshes. AI-assisted, unrefereed candidate; no external peer review, independent reproduction or formal verification is claimed.", "datePublished": "2026-10-22", "dateModified": "2026-10-22", "version": "1.0.0-candidate", "doi": "10.5281/zenodo.23565002", "doiUrl": "https://doi.org/10.5281/zenodo.23565002", "conceptDoi": "10.5281/zenodo.23565001", "pdfUrl": "https://github.com/ipitchford/tetrahedral-moment-stokes-elements/releases/download/v1.0.0-candidate/tetrahedral-moment-stokes-elements-1.0.0-candidate.pdf", "altPdfUrl": null, "altPdfLabel": null, "zenodoUrl": "https://zenodo.org/records/23565002", "repoUrl": "https://github.com/ipitchford/tetrahedral-moment-stokes-elements", "releaseUrl": "https://github.com/ipitchford/tetrahedral-moment-stokes-elements/releases/tag/v1.0.0-candidate", "markdownUrl": "https://evidencepress.org/releases/tetrahedral-moment-stokes-elements/index.md", "bibtexUrl": "https://evidencepress.org/releases/tetrahedral-moment-stokes-elements/cite.bib", "audioUrl": "https://evidencepress.org/assets/audio/tetrahedral-moment-stokes-elements.mp3", "imageUrl": "https://evidencepress.org/assets/og/tetrahedral-moment-stokes-elements.png?v=15178f2110", "coverArtUrl": "https://evidencepress.org/assets/art/tetrahedral-moment-stokes-elements.svg?v=d8c6f047b3", "media": [ { "type": "audio", "url": "https://evidencepress.org/assets/audio/tetrahedral-moment-stokes-elements.mp3", "name": "Plain-English audio briefing", "description": "Communication aid; not additional research evidence.", "transcriptUrl": "https://evidencepress.org/assets/audio/tetrahedral-moment-stokes-elements.txt" } ], "authors": [ "Anonymous" ], "license": "CC0-1.0", "status": "unrefereed-candidate", "verification": { "peerReviewed": false, "independentlyReproduced": false, "formallyVerified": false, "internallyReplayed": true, "detail": "Unrefereed AI-assisted mathematical candidate. Every theorem has a written proof; exact computations, numerical results and conjectures are labelled. Internal AI-assisted review, AI-assisted adversarial review and a verification re-review only. External specialist review, independent reproduction and formal verification have not been established." }, "assurance": [ { "dimension": "availability", "label": "Availability and archiving", "question": "Is the evidence package publicly retrievable from an archive under a persistent identifier?", "state": "passed", "evidenceUrl": "https://zenodo.org/records/23565002", "note": "The versioned research package and manuscript are publicly archived on GitHub and Zenodo. Provider asset parity was checked against frozen local files; availability does not establish correctness or peer review." }, { "dimension": "internalReplay", "label": "Internal replay", "question": "Does the producer’s own pipeline reproduce the stated result from the archived package?", "state": "passed", "note": "verify_all.py regenerates every certificate byte for byte. A clean full replay of the final package in a separate copy reproduced every output; the frozen archive was bound to it by hash, and the checks affected by publication changes were rerun." }, { "dimension": "independentRerun", "label": "Independent rerun", "question": "Has someone else run the supplied implementation and obtained the stated result?", "state": "not-assessed", "note": "No unaffiliated rerun has been established." }, { "dimension": "independentReimplementation", "label": "Independent reimplementation", "question": "Has someone else reached the result from an independent implementation?", "state": "partial", "note": "A separately coded producer-side check (independent_check.py) reconstructs the main ranks, ridges and the value 1/34650; no unaffiliated reimplementation." }, { "dimension": "formalVerification", "label": "Formal verification", "question": "Is a formalised statement machine-checked, and over which trusted base?", "state": "not-assessed", "note": "No proof-assistant formalisation." }, { "dimension": "specialistReview", "label": "Specialist review", "question": "Has a domain specialist assessed the argument?", "state": "not-assessed", "note": "AI-assisted internal and adversarial review is not external specialist review." }, { "dimension": "editorialPeerReview", "label": "Editorial peer review", "question": "Has a journal or venue run peer review to a decision?", "state": "not-assessed", "note": "No external editorial peer review." }, { "dimension": "dataEnvironmentReproducibility", "label": "Data and environment reproducibility", "question": "Are data and computational environment pinned well enough to rebuild?", "state": "partial", "note": "Exact checks need Python 3.12+ with the standard library and SymPy; the numerical test also needs NumPy and SciPy and about 7 GB RAM (versions listed, not locked in a container). Floating-point entries are rounded to four significant digits." } ], "provenance": { "aiGenerated": true, "aiAssisted": true, "generatedBy": [ "AI-assisted mathematical research", "AI-assisted internal review, verification re-review and adversarial review (including OpenAI gpt-5.6-sol)", "Claude Code: publication packaging and release" ], "humanRole": "Publication authority only; no human research contribution inferred.", "reviewStatus": "Internal AI-assisted and adversarial review, not external peer review." }, "attributionBoundary": null, "problem": null, "orderedProductAverage": null, "corrections": [], "publicCorrections": [], "publicCorrectionIndexes": null, "pageStructureVersion": "2.0", "pageStructureVariant": "canonical-theorem", "pageStructureWaivers": [], "recordMaturity": "current", "metadataProvenance": "Prepared from the approved final manuscript and its exact check outputs. Versioned GitHub and Zenodo research archives are checked against the frozen release assets on 22 October 2026.", "grandfatheredAtSchemaVersion": null, "keywords": [ "nonconforming finite elements", "Stokes equations", "unisolvence", "tetrahedral elements", "moment degrees of freedom", "inf-sup stability", "Jacobi polynomials", "Vandermonde polynomial" ], "keyResults": [ "The fourth-order question is answered in the negative: an alternating face functional vanishes on every quintic and equals 1/34650 on the Vandermonde sextic Δ₄ = ∏(λⱼ − λᵢ), so no 15 quintics complete the element (Proposition 3.2, proof; exact).", "For every n the moments are injective on P_{n+1}, their kernel on P_{n+2} is exactly seven-dimensional with an explicit Jacobi/Legendre ridge basis, they miss (n−1)(n−2)/2 directions on P_{n+2} and are surjective on P_{n+3}; the minimum degree is n+2 for n=1,2 and n+3 for n≥3 (Theorem 2.1, proof).", "A 50-dimensional fourth-order element with one sextic, invariant under vertex relabelling and minimal in degree and number of sextic directions (Corollary 5.8, proof); at least (n−1)(n−2)/2 top-degree directions are needed and numbering-free elements attaining this exist for every n (Corollary 5.2, Theorem 5.4, proof), explicit for n=4,5 and checked for Δ₄P_{n−3} up to n=8 (exact computation).", "On triangles no extra degree is ever needed (Corollary 6.3, proof); on d-simplices the kernel dimension and deficit have closed forms (proof), degree n+3 is checked for d=4,5,6 in listed cases (exact computation) and conjectured for d≥4 in general.", "A Stokes test with the fourth-order element shows rates consistent with fourth-order convergence and computed discrete inf-sup values 0.40–0.37 for N=1–4 (numerical); mesh independence comes from the Fortin argument (proof)." ], "reviews": [], "evidencePackage": "Written proofs of every theorem; exact rational and modular rank computations, the obstruction value 1/34650, the Gram matrix, a 14-monomial basis and lean elements for n=3..8, with an independently coded reconstruction of the main ranks; a floating-point Stokes convergence and inf-sup test; every certificate regenerated byte for byte by verify_all.py.", "openProblems": [], "mathObjects": [], "relatedWorks": [ { "citation": "Balazi, Allaire, Jolivet and Omnes (2025), Inf-sup stable non-conforming finite elements on tetrahedra with second- and third-order accuracy. Numerical Methods for Partial Differential Equations 41(6), e70057 (open version hal-04541809v3; source of the question).", "url": "https://doi.org/10.1002/num.70057", "doi": "10.1002/num.70057" }, { "citation": "Crouzeix and Raviart (1973), Conforming and nonconforming finite element methods for solving the stationary Stokes equations I. RAIRO Anal. Numér. 7(R3), 33–75.", "url": "https://doi.org/10.1051/m2an/197307R300331", "doi": "10.1051/m2an/197307R300331" }, { "citation": "Matthies and Tobiska (2005), Inf-sup stable non-conforming finite elements of arbitrary order on triangles. Numerische Mathematik 102, 293–309.", "url": "https://doi.org/10.1007/s00211-005-0648-8", "doi": "10.1007/s00211-005-0648-8" }, { "citation": "Ye and Zhang (2021), A P_{k+2} polynomial lifting operator on polygons and polyhedrons. Applied Mathematics Letters (injectivity of the same face and cell data).", "url": "https://doi.org/10.1016/j.aml.2021.107033", "doi": "10.1016/j.aml.2021.107033" }, { "citation": "Chen, Hu and Zhang (2023), Nonconforming finite element methods of order two and order three for the Stokes flow in three dimensions. Journal of Scientific Computing 97, art. 1.", "url": "https://doi.org/10.1007/s10915-023-02317-z", "doi": "10.1007/s10915-023-02317-z" }, { "citation": "Bohne, Ciarlet and Sauter (2026), Crouzeix–Raviart elements on simplicial meshes in d dimensions. Foundations of Computational Mathematics (arXiv:2407.04361).", "url": "https://doi.org/10.1007/s10208-026-09741-1", "doi": "10.1007/s10208-026-09741-1" }, { "citation": "Ciarlet, Dunkl and Sauter (2018), A family of Crouzeix–Raviart finite elements in 3D. Analysis and Applications 16, 649–691.", "url": "https://doi.org/10.1142/S0219530518500070", "doi": "10.1142/S0219530518500070" }, { "citation": "Sauter and Torres (2023), On the inf-sup stability of Crouzeix–Raviart Stokes elements in 3D. Mathematics of Computation 92, 1033–1059.", "url": "https://doi.org/10.1090/mcom/3793", "doi": "10.1090/mcom/3793" }, { "citation": "Aktaş and Xu (2013), Sobolev orthogonal polynomials on a simplex. International Mathematics Research Notices.", "url": "https://doi.org/10.1093/imrn/rns141", "doi": "10.1093/imrn/rns141" }, { "citation": "Cangiani, Gyrya and Manzini (2016), The nonconforming virtual element method for the Stokes equations. SIAM Journal on Numerical Analysis 54, 3411–3435.", "url": "https://doi.org/10.1137/15M1049531", "doi": "10.1137/15M1049531" } ], "claimAssurance": null, "operatingModel": { "version": "1.0", "workId": "ep-work:tetrahedral-moment-stokes-elements", "attemptIds": [ "ep-attempt:tetrahedral-moment-stokes-elements-publication" ], "aims": [ "science" ], "artifactRoles": [ "research-output", "communication" ], "lineageId": null, "accelerationPrimitives": [ "counterexample-proxy-first", "adversarial-controls", "assurance-vector", "agent-readable-research-object" ], "decisionObject": { "type": "other", "description": "Answer to the open question of Balazi, Allaire, Jolivet and Omnes (2025) on quintic enrichment of the fourth-order tetrahedral moment element, with the exact kernel of face-and-cell moments on P_{n+2} for every n, the sharp minimal degree, minimal enrichment counts with symmetric elements, a d-simplex analysis and an explicit fourth-order element with a numerical Stokes test.", "scope": "Face moments against P_n and cell moments against P_{n-1} on a tetrahedron; single-cell polynomial spaces containing P_{n+1}; no claim about split-cell, nonpolynomial or virtual element spaces, pressure robustness, p-uniform stability or solver speed; d≥4 sufficiency of degree n+3 conjectural beyond checked cases." }, "bottleneckTargeted": [ "assurance", "publication" ], "semanticBridge": { "state": "explicit", "description": "The source question is quoted verbatim from Appendix A of the open version (hal-04541809v3) with Remarks 3.2 and 3.6; the source's notation V_4(K) and Σ_5(K) is mapped to the paper's.", "remainingRisks": [ "The journal version was compared through the open version's numbering.", "Ye and Zhang (2021) and Chen, Hu and Zhang (2023) overlap on injectivity and low-order kernels; the distinctions are stated in the paper.", "The bounded literature check is not exhaustive priority certification." ] }, "humanJudgmentGates": [ "Publication is authorised; scientific endorsement is not inferred.", "External specialist review and independent reproduction remain open." ], "parentLinks": [], "assuranceTarget": { "dimensions": [ "independentRerun", "specialistReview", "formalVerification", "priorityAssessment" ], "nextAction": "Have a finite element specialist check the kernel theorem (Theorem 2.1, Section 4) and the lean-element existence proof (Theorem 5.4).", "claimCeiling": "The exact moment kernel, sharp degree and minimal enrichment statements for face-and-cell moments on tetrahedra within the stated single-cell polynomial scope, with labelled exact, numerical and conjectural extensions. Internal verification is not external assurance." }, "impactClaims": [ { "id": "tetrahedral-moment-stokes-elements-no-impact", "aim": "science", "outcome": "An openly inspectable answer to a published open question on tetrahedral nonconforming Stokes elements.", "setting": "AI-assisted mathematical research", "status": "NO_IMPACT_EVIDENCE", "designClass": "none", "comparator": "None registered.", "estimand": "No research-speed or scientific-impact effect is estimated.", "evidenceRefs": [], "registeredDesignRef": null, "independentAssessment": null } ] }, "researchMetrics": { "schemaVersion": "1.0", "policyUrl": "https://evidencepress.org/api/research-metrics-policy.json", "outcomePolicyUrl": "https://evidencepress.org/api/v2/research-metrics-policy.json", "workLedgerUrl": "https://evidencepress.org/api/work-ledger.json", "claimCeiling": "These records describe scoped workflow inputs and outcomes. They support calibration and reference-class construction, but do not by themselves establish research acceleration, theorem correctness, novelty, independent assurance or impact.", "attempts": [ { "attemptId": "ep-attempt:tetrahedral-moment-stokes-elements-publication", "status": "active", "resultClass": "pending", "metrics": { "schemaVersion": "1.0", "measurementScope": "publication-only", "scopeBoundary": "Only archive packaging, site integration, guarded delivery and canonical readback from registration are prospective. Earlier research, review, revision and extension work are left-censored.", "forecast": { "frozenAt": "2026-10-22T13:22:15Z", "procedureClass": "full-candidate-publication", "targetOutcome": "Accurately scoped GitHub, Zenodo and canonical Evidence Press release with explanatory media and machine readback, if all publication gates pass.", "expectedActiveMinutes": 90, "plausibleLowMinutes": 60, "plausibleHighMinutes": 150, "expectedUnattendedWaitMinutes": 30, "referenceClass": { "label": "Full mathematical candidate publication", "sampleSize": 0, "basis": "Route benchmark planning prior, not empirical calibration." }, "fermiComponents": [ { "component": "Research archive packaging, reuse-bound replay and identities", "count": 1, "lowMinutesPerUnit": 25, "centralMinutesPerUnit": 35, "highMinutesPerUnit": 60, "basis": "One frozen package; the full deterministic suite is reused from a separate clean replay by hash binding, with only changed checks rerun after DOI insertion." }, { "component": "Page, metadata and media integration", "count": 1, "lowMinutesPerUnit": 20, "centralMinutesPerUnit": 30, "highMinutesPerUnit": 45, "basis": "One reader-first release, drafted in parallel with verification." }, { "component": "Guarded deployment and canonical readback", "count": 1, "lowMinutesPerUnit": 15, "centralMinutesPerUnit": 25, "highMinutesPerUnit": 45, "basis": "Mandatory site tests, one merge cycle with ledger carry-forward, zero-cost deployment and public checks." } ], "probabilityPositiveSignal": 0.9, "probabilityTargetClosure": 0.85, "probabilityHorizonMinutes": 300, "assumptions": [ "The approved scientific package remains valid and configured providers remain available." ], "reforecastTriggers": [ "A publication gate reveals a scientific or integrity defect, or a provider outage blocks delivery." ], "stopRule": "Do not weaken assurance, integrity or cost gates to meet the forecast." }, "outcome": null }, "missingnessReason": null } ] } }