Dra. Camaño-Valenzuela, Jessika

A method classically used in the lower polynomial degree for the construction of a finite element basis of the space of divergence-free functions is here extended to any polynomial degree for a bounded domain without topological restrictions. The method uses graphs associated with two differential operators: the gradient and the divergence, and selects the basis using a spanning tree of the first graph. It can be applied for the two main families of degrees of freedom, weights and moments, used to express finite element differential forms.

We propose and analyze an efficient algorithm for the computation of a basis of the space of divergence-free Raviart–Thomas finite elements. The algorithm is based on graph techniques. The key point is to realize that, with very natural degrees of freedom for fields in the space of Raviart–Thomas finite elements of degree r+1r+1 and for elements of the space of discontinuous piecewise polynomial functions of degree r≥0r≥0, the matrix associated with the divergence operator is the incidence matrix of a particular graph. By choosing a spanning tree of this graph, it is possible to identify an invertible square submatrix of the divergence matrix and to compute easily the moments of a field in the space of Raviart–Thomas finite elements with assigned divergence. This approach extends to finite elements of high degree the method introduced by Alotto and Perugia (Calcolo 36:233–248, 1999) for finite elements of degree one. The analyzed approach is used to construct a basis of the space of divergence-free Raviart–Thomas finite elements. The numerical tests show that the performance of the algorithm depends neither on the topology of the domain nor or the polynomial degree r.