Sabine Zaglmayr
Technische Universität Graz
Reduced basis gauging for high-order Maxwell solvers
In order to guartantee unique solvability for the magnetostatic Maxwell
equations additional constraints, so called gauging conditions, have to
be imposed. In particular, we consider Coulomb gauge which enforces
orthogonality to gradient fields and discuss several realizations, e.g.
via Lagrange multipliers.
By a careful construction of high order Nedelec-type finite elements the
Coulomb gauge can be realized by a two step strategy: First, uniqueness
of the discrete magnetostatic problem can be ensured by
eliminating part of the high order basis functions. In a second step the
orthogonality to gradients is restored by postprocessing.
Both subproblems use reduced bases and are better conditioned than the
original problem, so the overall performance can be improved considerably.
The efficiency of our approach is illustrated by numerical examples.