Numerical Considerations About Using Finite-Element Methods to Compute AC Losses in HTS

In this paper, we investigate the influence of various numerical parameters on the precision and the speed of ac loss computations in high-temperature superconductors using the finite-element method. The case considered here is an infinite slab subjected to an external ac magnetic field. This problem can be modeled by a 1-D partial differential equation (diffusion equation). This relatively simple case allowed investigating the influence of the various parameters in a reasonable time. The findings of this work can be used as a starting point for optimizing the numerical settings of more complex models. As the main results, it is shown that choosing first-order elements for approximating the flux density ( B) is the most stable option, and that using high order adaptive time-stepping methods provides good accuracy and fast simulations. A new and simple self-check test for validating the computed ac losses is also proposed. Finally, a detailed analysis about the behavior of the numerical solution near flux/current fronts is provided.