Geometrically defined basis functions for polyhedral elements with applications to computational electromagnetics

In the recent years, reformulating the mathematical description of physical laws in an algebraic form using tools from algebraic topology gained popularity in computational physics. Physical variables are defined as fluxes or circulations on oriented geometric elements of a pair of dual interlocked cell complexes, while physical laws are expressed in a metric-free fashion with incidence matrices. The metric and the material information are encoded in the discrete counterpart of the constitutive laws of materials, also referred to as constitutive or material matrices. The stability and consistency of the method is guaranteed by precise properties (symmetry, positive definiteness, consistency) that material matrices have to fulfill. The main advantage of this approach is that material matrices, even for arbitrary star-shaped polyhedral elements, can be geometrically defined, by simple closed-form expressions, in terms of the geometric elements of the primal and dual grids. That is why this original technique may be considered as a “Discrete Geometric Approach” (DGA) to computational physics. This paper first details the set of vector basis functions associated with the edges and faces of a polyhedral primal grid or of a dual grid. Then, it extends the construction of constitutive matrices for bianisotropic media.