Explicit geometric construction of sparse inverse mass matrices for arbitrary tetrahedral grids

12/02/2020
by   Silvano Pitassi, et al.
0

The geometric reinterpretation of the Finite Element Method (FEM) shows that Raviart Thomas and Nedelec mass matrices map from degrees of freedoms (DoFs) attached to geometric elements of a tetrahedral grid to DoFs attached to the barycentric dual grid. The algebraic inverses of the mass matrices map DoFs attached to the barycentric dual grid back to DoFs attached to the corresponding primal tetrahedral grid, but they are of limited practical use since they are dense. In this paper we present a new geometric construction of sparse inverse mass matrices for arbitrary tetrahedral grids and possibly anisotropic materials, debunking the conventional wisdom that the barycentric dual grid prohibits a sparse representation for inverse mass matrices. In particular, we provide a unified framework for the construction of both edge and face mass matrices and their sparse inverses. Such a unifying principle relies on novel geometric reconstruction formulas, from which, according to a well established design strategy, local mass matrices are constructed as the sum of a consistent and a stabilization term. A major difference with the approaches proposed so far is that the consistent term is defined geometrically and explicitly, that is, without the necessity of computing the inverses of local matrices. This provides a sensible speedup and an easier implementation. We use these new sparse inverse mass matrices to discretize a three dimensional Poisson problem, providing the comparison between the results obtained by various formulations on a benchmark problem with analytical solution.

READ FULL TEXT

page 6

page 30

page 31

page 32

research
09/21/2022

A Yee-like finite element scheme for Maxwell's equations on hybrid grids

A novel finite element method for the approximation of Maxwell's equatio...
research
11/10/2019

MATLAB Implementation of Element-based Solvers

Rahman and Valdman (2013) introduced a vectorized way to assemble finite...
research
06/01/2023

A Yee-like finite-element scheme for Maxwell's equations on unstructured grids

A novel finite element scheme is studied for solving the time-dependent ...
research
11/04/2021

A closed-form multigrid smoothing factor for an additive Vanka-type smoother applied to the Poisson equation

We consider an additive Vanka-type smoother for the Poisson equation dis...
research
12/03/2017

Generalised primal-dual grids for unstructured co-volume schemes

The generation of high-quality staggered unstructured grids for computat...
research
03/28/2018

Investigating the hybrid textures of neutrino mass matrix for near maximal atmospheric neutrino mixing

In the present paper, we have studied that the implication of a large va...
research
01/21/2021

Simple finite elements and multigrid for efficient mass-consistent wind downscaling in a coupled fire-atmosphere model

We present a simple finite element formulation of mass-consistent approx...

Please sign up or login with your details

Forgot password? Click here to reset