Programming finite element method based hysteresis loss computation software using non-linear superconductor resistivity and T − φ formulation

Due to the rapid development of personal computers from the beginning of the 1990s, it has become a reality to simulate current penetration, and thus hysteresis losses, in superconductors with other than very simple one-dimensional (1D) Bean model computations or Norris formulae. Even though these older approaches are still usable, they do not consider, for example, multifilamentary conductors, local critical current dependency on magnetic field or varying n-values. Currently, many numerical methods employing different formulations are available. The problem of hysteresis losses can be scrutinized via an eddy current formulation of the classical theory of electromagnetism. The difficulty of the problem lies in the non-linear resistivity of the superconducting region. The steep transition between the superconducting and the normal states often causes convergence problems for the most common finite element method based programs. The integration methods suffer from full system matrices and, thus, restrict the number of elements to a few thousands at most. The so-called T ? formulation and the use of edge elements, or more precisely Whitney 1-forms, within the finite element method have proved to be a very suitable method for hysteresis loss simulations of different geometries. In this paper we consider making such finite element method software from first steps, employing differential geometry and forms.

[1]  Luciano Martini,et al.  Integral equations for the current density in thin conductors and their solution by the finite-element method , 2008 .

[2]  Naoyuki Amemiya,et al.  AC loss reduction of YBCO coated conductors by multifilamentary structure , 2004 .

[3]  Pascal Tixador,et al.  Thermal-electromagnetic modeling of superconductors , 2007 .

[4]  A. Campbell,et al.  A new method of determining the critical state in superconductors , 2007 .

[5]  Hybrid and integral formulations for 3D eddy current problems , 1997 .

[6]  Adel Razek,et al.  Computation of 3-D electromagnetic field using differential forms based elements and dual formulations , 1996 .

[7]  M. Majoros,et al.  Numerical Estimation of AC Loss in MgB2 Wires in Self-field Condition , 2008 .

[8]  P. Bruzzone.  The index n of the voltage–current curve, in the characterization and specification of technical superconductors , 2004 .

[9]  O. Bíró.  Edge element formulations of eddy current problems , 1999 .

[10]  G. Deschamps.  Electromagnetics and differential forms , 1981, Proceedings of the IEEE.

[11]  Bertrand Dutoit,et al.  Finite-element modelling of YBCO fault current limiter with temperature dependent parameters , 2007 .

[12]  Lawrence F. Shampine,et al.  Analysis and implementation of TR-BDF2 , 1996 .

[13]  W. T. Norris,et al.  Calculation of hysteresis losses in hard superconductors carrying ac: isolated conductors and edges of thin sheets , 1970 .

[14]  Luciano Martini,et al.  Development of an edge-element model for AC loss computation of high-temperature superconductors , 2006 .

[15]  C. F. Hempstead,et al.  RESISTIVE STATES OF HARD SUPERCONDUCTORS , 1964 .

[16]  Lauri Kettunen,et al.  Computation of AC losses in high-temperature superconductors , 1998 .

[17]  H. Kardestuncer,et al.  Finite element handbook , 1987 .

[18]  Brandt.  Thin superconductors in a perpendicular magnetic ac field: General formulation and strip geometry. , 1994, Physical review. B, Condensed matter.

[19]  D. Arnold,et al.  Finite element exterior calculus, homological techniques, and applications , 2006, Acta Numerica.

[20]  D. Larbalestier,et al.  Resistive transitions and the origin of the n value in superconductors with a Gaussian critical‐current distribution , 1993 .

[21]  H. Hashizume,et al.  Numerical analysis of electromagnetic phenomena in superconductors , 1992 .

[22]  Leonid Prigozhin.  Regular ArticleThe Bean Model in Superconductivity: Variational Formulation and Numerical Solution , 1996 .

[23]  Enric Pardo,et al.  Magnetic properties of arrays of superconducting strips in a perpendicular field , 2003 .

[24]  R. Mikkonen,et al.  Quench current in conduction-cooled HTS magnets , 2003 .

[25]  Luciano Martini,et al.  AC losses in thin superconductors: the integral equation method applied to stacks and windings , 2009 .

[26]  C. P. Bean.  Magnetization of hard superconductors , 1962 .

[27]  H. London,et al.  Alternating current losses in superconductors of the second kind , 1963 .

[28]  N. Amemiya,et al.  Electromagnetic field analysis of rectangular superconductor with large aspect ratio in arbitrary orientated magnetic fields , 2005, IEEE Transactions on Applied Superconductivity.

[29]  Leonid Prigozhin,et al.  The Bean Model in Superconductivity , 1996 .

[30]  Lauri Kettunen,et al.  Formulation of the eddy current problem in multiply connected regions in terms of h , 1998 .