An intrinsic finite element framework for scalar- and vector-valued partial differential equations on general manifolds
Abstract
We present an intrinsic finite element framework for the numerical approximation of scalar- and vector-valued partial differential equations on general manifolds described by atlases of charts. Weak formulations and their discretisation are expressed exclusively in the flat parametric space of the manifold, where the exact geometry enters only through the metric tensor that is evaluated at quadrature points. Conforming finite element spaces of arbitrary order for the whole de Rham complex are obtained on multi-chart atlases through purely topological inter-chart gluing. We develop a nodal change of basis construction for vector-valued spaces on multi-chart atlases, where nodal degrees of freedom couple charts with incompatible coordinate systems through transmission maps. Inter-chart continuity holds at the nodes with jumps of the order of the approximation error, and the resulting discrete field is exactly tangent to the manifold. Optimal convergence rates for this method are observed numerically for the surface Stokes problem; no a priori error analysis is provided. Numerical experiments on the cubed sphere manifold in two and three dimensions for the Hodge Laplacian problems demonstrate the intrinsic framework is free of the geometric consistency errors for a sufficient quadrature degree, and that intrinsic assembly is substantially cheaper than its extrinsic counterpart. Further considering the rotating shallow water equations that include orography as a perturbation of the geometry demonstrates that mass is conserved exactly, and energy is conserved up to the time discretisation error.
Related
- GridapGeosciences.jl: A Julia finite element package for partial differential equations on general manifolds
- Differential Harnack Estimates and Frequency Monotonicity for Filtration Equations on Riemannian Manifolds
- Bounds-Constrained Finite Element Approximation of Time-Dependent Partial Differential Equations
- Rough differential equations on manifolds via natural bundles
- The Priority-Axis Z-Algebra: A Unified Framework of Axioms, Structure Theory, and Intrinsic Differential Geometry
BibTeX
@misc{tambyah2026intrinsic,
title = {{An intrinsic finite element framework for scalar- and vector-valued partial differential equations on general manifolds}},
author = {Tamara A. Tambyah and Alberto F. Martín and David Lee and Santiago Badia},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.09592},
eprint = {2609.09592},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.09592},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - An intrinsic finite element framework for scalar- and vector-valued partial differential equations on general manifolds
AU - Tamara A. Tambyah
AU - Alberto F. Martín
AU - David Lee
AU - Santiago Badia
PY - 2026
DA - 2026/09/09
PB - arXiv
DO - 10.48550/arxiv.2609.09592
UR - https://arxiv.org/abs/2609.09592
AB - We present an intrinsic finite element framework for the numerical approximation of scalar- and vector-valued partial differential equations on general manifolds described by atlases of charts. Weak formulations and their discretisation are expressed exclusively in the flat parametric space of the manifold, where the exact geometry enters only through the metric tensor that is evaluated at quadrature points. Conforming finite element spaces of arbitrary order for the whole de Rham complex are obtained on multi-chart atlases through purely topological inter-chart gluing. We develop a nodal change of basis construction for vector-valued spaces on multi-chart atlases, where nodal degrees of freedom couple charts with incompatible coordinate systems through transmission maps. Inter-chart continuity holds at the nodes with jumps of the order of the approximation error, and the resulting discrete field is exactly tangent to the manifold. Optimal convergence rates for this method are observed numerically for the surface Stokes problem; no a priori error analysis is provided. Numerical experiments on the cubed sphere manifold in two and three dimensions for the Hodge Laplacian problems demonstrate the intrinsic framework is free of the geometric consistency errors for a sufficient quadrature degree, and that intrinsic assembly is substantially cheaper than its extrinsic counterpart. Further considering the rotating shallow water equations that include orography as a perturbation of the geometry demonstrates that mass is conserved exactly, and energy is conserved up to the time discretisation error.
ER -