Francesco Bonaldi,Daniele A. Di Pietro,Jerome Droniou,Kaibo Hu
Abstract:We develop in this work the first polytopal complexes of differential forms. These complexes, inspired by the Discrete De Rham and the Virtual Element approaches, are discrete versions of the de Rham complex of differential forms built on meshes made of general polytopal elements. Both constructions benefit from the high-level approach of polytopal methods, which leads, on certain meshes, to leaner constructions than the finite element method. We establish commutation properties between the interpolators and the discrete and continuous exterior derivatives, prove key polynomial consistency results for the complexes, and show that their cohomologies are isomorphic to the cohomology of the continuous de Rham complex.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a new discrete framework for dealing with differential forms in the polyhedral method of arbitrary dimensions and arbitrary approximation orders. Specifically, the author aims to combine the Discrete de Rham complex with the Virtual Element Method (VEM) to construct a discrete differential form complex that can work on general polyhedral meshes.
### Core problems of the paper
1. **Combining two successful research directions**: Finite Element Differential Forms and arbitrary - order polyhedral methods. Both of these are very successful fields in numerical analysis, but they usually rely on conforming meshes, while polyhedral methods can work on more general meshes.
2. **Maintaining the homology structure**: Ensure that the homology space of the discrete complex is isomorphic to that of the continuous de Rham complex. This is crucial for designing stable numerical schemes because these homology structures are closely related to the topological features of the problem.
3. **Developing a compatible framework**: Although polyhedral techniques introduce some non - conformity, they can still be used to develop a compatible framework. This enables the discrete complex to maintain important mathematical properties such as polynomial consistency and compatibility.
### Specific objectives
- **Constructing the discrete de Rham complex**: The author proposes two discrete de Rham complexes of arbitrary dimensions and arbitrary approximation orders, based on the Discrete de Rham method (DDR) and the Virtual Element Method (VEM) respectively. These constructions do not require an explicit global conforming differential form space but use a fully discrete space that contains polynomial components on mesh cells.
- **Defining the discrete exterior derivative and the discrete potential**: By combining polynomial components to mimic the Stokes formula, the discrete exterior derivative and the discrete potential are constructed. The design of these discrete operators revolves around two key properties: polynomial consistency and compatibility.
- **Applying to practical problems**: These discrete complexes can be applied to various models, such as magnetostatics, Stokes equations, and Yang - Mills equations, etc., and can remain robust when physical parameters change.
### Conclusion
The main contribution of this paper is to provide a new framework for constructing discrete differential form complexes on general polyhedral meshes. This method not only supports a wider range of mesh types but also can maintain important mathematical and physical properties, thus providing a powerful tool for solving complex partial differential equations.
### Key formulas
- Definition of the discrete exterior derivative:
\[
d_h: \mathcal{P}^{-}_{k}\Lambda^{\ell}(\mathcal{T}_h) \to \mathcal{P}^{-}_{k}\Lambda^{\ell + 1}(\mathcal{T}_h)
\]
- Definition of the discrete potential:
\[
\Pi_h: \mathcal{P}^{-}_{k}\Lambda^{\ell}(\mathcal{T}_h) \to \mathcal{P}_{k}\Lambda^{\ell}(\mathcal{T}_h)
\]
These formulas ensure the polynomial consistency and compatibility of the discrete complex, thereby guaranteeing the stability and accuracy of the numerical scheme.