H method for solving 3D eddy current problems

A novel finite element model, which uses magnetic field intensity H as a state variable in the whole solution region of 3D eddy current problems, is presented. A surface integral on the interface, which leads to asymmetric equations in discretization, is considered in this formulation. The formulation uses consistently linear edge elements and flux continuity is ensured by imposing ∇·B=0 with the penalty method. Numerical results of the problem No. 3 of the Team Workshop are provided, showing the validity of the algorithm