COMSOL Implementation of the H-$\phi$-Formulation With Thin Cuts for Modeling Superconductors With Transport Currents
A. Larry Arsenault, Bruno de Sousa Alves, Frédéric Sirois
Abstract
Despite the acclaimed success of the magnetic field (H) formulation for modeling the electromagnetic behavior of superconductors with the finite-element method, the use of vector-dependent variables in nonconducting domains leads to unnecessarily long computation times. In order to solve this issue, we have recently shown how to use a magnetic scalar potential together with the H-formulation in the COMSOL Multiphysics environment to efficiently and accurately solve for the magnetic field surrounding superconducting domains. However, from the definition of the magnetic scalar potential, the nonconducting domains must be made simply connected in order to obey Ampere's law. In this article, we use thin cuts to apply a discontinuity in Φ and make the nonconducting domains simply connected. This approach is shown to be easily implementable in the COMSOL Multiphysics finite-element program, already widely used by the applied superconductivity community. We simulate three different models in two dimensions and three dimensions using superconducting filaments and tapes, and show that the results are in very good agreement with the H-A and H-formulations. Finally, we compare the computation times between the formulations, showing that the H-Φ-formulation can be up to seven times faster than the standard H-formulation in certain applications of interest.