A FIT formulation of bianisotropic materials over polyhedral grids

Modeling bianisotropic constitutive equations, i.e., magnetoelectric coupling, in electromagnetics simulation is increasingly important, in particular in metamaterials applications. This paper introduces for the first time such constitutive relationships in the framework of finite integration technique and, furthermore, does so by allowing full generality in the discretization through arbitrary polyhedral grids. The resulting formulation is consistent, stable, and preserves the thermodynamic properties of the bianisotropic constitutive equations because of the energetic approach used to construct the interpolating functions.