A Finite Element Simulation Tool for Predicting Hysteresis Losses in Superconductors Using an H-Oriented Formulation with Cohomology Basis Functions

Currently, modelling hysteresis losses in superconductors is most often based on the H-formulation of the eddy current model (ECM) solved using the finite element method (FEM). In the H-formulation, the problem is expressed using the magnetic field intensity H and discretized using edge elements in the whole domain. Even though this approach is well established, it uses unnecessary degrees of freedom (DOFs) and introduces modelling error such as currents flowing in air regions due to finite air resistivity. In this paper, we present a modelling tool utilizing another H-oriented formulation of the ECM, making use of cohomology of the air regions. We constrain the net currents through the conductors by fixing the DOFs related to the so-called cohomology basis functions. As air regions will be truly non-conducting, DOFs and running times of these nonlinear simulations are reduced significantly as compared to the classical H-formulation. This fact is demonstrated through numerical simulations.

[1]  Weijia Yuan,et al.  Computation of Losses in HTS Under the Action of Varying Magnetic Fields and Currents , 2013, IEEE Transactions on Applied Superconductivity.

[2]  Min Zhang,et al.  3D modeling of high-Tc superconductors by finite element software , 2011 .

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

[4]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[5]  Stephan Foldes.  Fundamental Structures of Algebra and Discrete Mathematics: Foldes/Fundamental , 1994 .

[6]  Antti Stenvall,et al.  Comparison of three eddy current formulations for superconductor hysteresis loss modelling , 2012 .

[7]  Christophe Geuzaine,et al.  Homology and Cohomology Computation in Finite Element Modeling , 2013, SIAM J. Sci. Comput..

[8]  Antti Stenvall,et al.  Development of a three-dimensional finite-element model for high-temperature superconductors based on the H-formulation , 2013 .

[9]  Antti Stenvall,et al.  An H-formulation-based three-dimensional hysteresis loss modelling tool in a simulation including time varying applied field and transport current: the fundamental problem and its solution , 2014 .

[10]  Lauri Kettunen,et al.  Discrete spaces for div and curl-free fields , 1998 .

[11]  S. Dou,et al.  Power-Law Relationship Between Critical Current Density, Microstructure, and the n-Value in MgB2 Superconductor Wires , 2014, Journal of Superconductivity and Novel Magnetism.

[12]  B. Auchmann,et al.  Differential Forms and Boundary Integral Equations for Maxwell-Type Problems , 2014, 1411.2661.

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

[14]  Antti Stenvall,et al.  Programming finite element method based hysteresis loss computation software using non-linear superconductor resistivity and T − φ formulation , 2010 .

[15]  Siddhartha Sen,et al.  Topology and geometry for physicists , 1983 .

[16]  Antti Stenvall,et al.  Scientific Research in the Field of Mesh Method Based Modeling of AC Losses in Superconductors: A Review , 2014 .

[17]  Lauri Kettunen,et al.  A programming interface to the Riemannian manifold in a finite element environment , 2013, J. Comput. Appl. Math..

[18]  Carol S. Woodward,et al.  Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers , 2020, ACM Trans. Math. Softw..

[19]  Linda R. Petzold,et al.  Numerical solution of initial-value problems in differential-algebraic equations , 1996, Classics in applied mathematics.

[20]  Kay Hameyer,et al.  An algorithm to construct the discrete cohomology basis functions required for magnetic scalar potential formulations without cuts , 2002 .

[21]  Charles Leake.  Fundamental Structures of Algebra and Discrete Mathematics , 1994 .

[22]  T. Tarhasaari,et al.  Approaches for Tree-Co-Tree Gauged $T$ –$\varphi$-Formulated Eddy-Current Problem in Superconductor Hysteresis Loss Simulations , 2013, IEEE Transactions on Applied Superconductivity.

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

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

[25]  P. R. Kotiuga.  On making cuts for magnetic scalar potentials in multiply connected regions , 1987 .

[26]  T. A. Coombs,et al.  Numerical solution of critical state in superconductivity by finite element software , 2006 .

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

[28]  Paul W. Gross,et al.  Electromagnetic Theory and Computation: A Topological Approach , 2004 .