IEEE TRANSACTIONS ON APPLIED SUPERCONDUCTIVITY, VOL. 24, NO. 1, FEBRUARY (Invited Paper)

Transcription

1 Computation of Losses in HTS Under the Action of Varying Magnetic Fields and Currents Francesco Grilli, Enric Pardo, Member, IEEE, Antti Stenvall, Doan N. Nguyen, Weijia Yuan, and Fedor Gömöry, Member, IEEE (Invited Paper) Abstract Numerical modeling of superconductors is widely recognized as a powerful tool for interpreting experimental results, understanding physical mechanisms, and predicting the performance of high-temperature-superconductor (HTS) tapes, wires, and devices. This is particularly true for ac loss calculation since a sufficiently low ac loss value is imperative to make these materials attractive for commercialization. In recent years, a large variety of numerical models, which are based on different techniques and implementations, has been proposed by researchers around the world, with the purpose of being able to estimate ac losses in HTSs quickly and accurately. This paper presents a literature review of the methods for computing ac losses in HTS tapes, wires, and devices. Technical superconductors have a relatively complex geometry (filaments, which might be twisted or transposed, or layers) and consist of different materials. As a result, different loss contributions exist. In this paper, we describe the ways of computing such loss contributions, which include hysteresis losses, eddy-current losses, coupling losses, and losses in ferromagnetic materials. We also provide an estimation of the losses occurring in a variety of power applications. Index Terms Alternate current (ac) losses, coupling losses, eddy-current losses, hysteresis losses, magnetic materials, numerical modeling. I. INTRODUCTION THE capability of computing ac losses in high-temperature superconductors (HTSs) is very important for designing Manuscript received September 24, 2012; accepted April 16, Date of current version October 23, The work of F. Grilli was supported by the Helmholtz University Young Investigators Program under Grant VH-NG-617. The work of E. Pardo was supported by the Agency for the Structural Funds of the European Union from the Ministry of Education, Science, Research and Sports of the Slovak Republic under Contract The work of A. Stenvall was supported in part by the Academy of Finland under Project , by the Foundation for Technology Promotion in Finland, and by the Emil Aaltonen Foundation. The work of D. Nguyen was supported by the Los Alamos National Laboratory through the Laboratory Directed Research and Development Program under Grant ER. This paper was recommended by Associate Editor P. J. Masson. F. Grilli is with the Institute for Technical Physics, Karlsruhe Institute of Technology, Karlsruhe, Germany ( francesco.grilli@kit.edu). E. Pardo and F. Gömöry are with the Institute of Electrical Engineering, Slovak Academy of Sciences, Bratislava, Slovak Republic. A. Stenvall is with the Tampere University of Technology, Tampere, Finland. D. N. Nguyen is with the Los Alamos National Laboratory, Los Alamos, NM USA. W. Yuan is with the University of Bath, Bath BA2 7AY, U.K. Color versions of one or more of the figures in this paper are available online at Digital Object Identifier /TASC and manufacturing marketable devices such as cables, fault current limiters, transformers, and motors. As a matter of fact, in many cases, the prospected ac loss value is too high to make the application attractive on the market, and numerical calculations can help to find solutions for reducing the losses. In the past decades, several analytical models for computing the losses in superconductor materials have been developed. In general, these models can provide the loss value for a given geometry in given working conditions by means of relatively simple formulas. While these models are undoubtedly fast and very useful for basic understanding of the loss mechanisms and for predicting the loss value for a certain number of geometries and tape arrangements, they have several important limitations, which restrict their usefulness for an accurate estimation of the ac losses in real HTS devices. Numerical models, on the other hand, can overcome these limitations and can simulate geometries and situations of increasing complexity. This increased computing capability comes at the price of more complex software implementation and longer computation times. This paper focuses on methods for computing ac losses by means of analytical and numerical techniques. Technical superconductors consist of different materials in addition to the superconductor itself (metal, buffers and substrates, magnetic materials, etc.), some of which, depending on the operating conditions, can give a significant contribution to the total losses. These loss contributions can be divided into four categories, which are summarized in the following list and schematically illustrated in Fig. 1: 1) Hysteresis losses, which are caused by the penetration and movement of the magnetic flux in the superconducting material; 2) Eddy-current losses, which are caused by the currents induced by a magnetic field and circulating in the normal metal parts of a superconducting tape; 3) Coupling losses, which are caused by the currents coupling two or more superconducting filaments via the normal metal regions separating them; 4) Ferromagnetic losses, which are caused by the hysteresis cycles in magnetic materials. This paper is structured as follows. Sections II V focus each on a different loss contribution, as in the list above. Section VI gives the loss estimation for various power applications. The two appendices contain the details of calculations used in Section II-B IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See for more information.

2 GRILLI et al.: COMPUTATION OF LOSSES IN HTS Fig. 1. Schematic of the different loss contributions in technical HTSs, e.g., ReBCO-coated conductor (tape s cross section displayed): 1) Hysteresis losses in the superconductor; 2) eddy current losses in the normal metal stabilizer; 3) coupling losses between filaments (e.g., through conducting paths due to imperfections of the striation process); and 4) ferromagnetic losses in the substrate. II. HYSTERESIS LOSSES This section is divided into different parts as follows. First, we describe the two most commonly used models for describing a superconductor with the purpose of ac loss computation. Then, we explain how to solve the electromagnetic quantities. Afterward, we illustrate how to compute the ac losses once the electromagnetic quantities are known. Finally, we show the solutions for several particular cases relevant for applications. A. Models for the Superconductor The reason why hard superconductors are able to carry large electrical currents is that the magnetic field (penetrating type-ii superconductors in the form of supercurrent vortices, each carrying the same amount of magnetic flux) is pinned in the volume of the superconducting material. Due to this mechanism, the vortices do not move under the action of the Lorentz force that is pushing them in the direction perpendicular to both the flow of the electrical current and the magnetic field. In addition, a gradient in the density of vortices is reluctant to any rearrangement. Thus, once a dc is established in a hard superconductor, it will also persist after switching the driving voltage off (circulation of an electrical current with less resistance). However, in the ac regime, the vortices must move to follow the change of the magnetic field: The pinning force represents an obstacle, and superseding it is an irreversible process. The accompanied dissipation is called the hysteresis loss in hard superconductors. It is not easy to link the interaction between a current vortex (and the involved pinning centers) and the material s properties that can be used in electromagnetic calculations. However, in the investigation of low-temperature superconductors (LTSs), such as NbTi or Nb 3 Sn, it was found that, for practical purposes, one can obtain very useful predictions utilizing the phenomenological description introduced in [1] and [2], commonly known as the critical-state model (CSM). The model is valid on a macroscopic scale that neglects the details of electrical current distribution in individual vortices and replaces it with an average taken over a large number of vortices. The CSM Fig. 2. Current density distribution in the cross section of a round superconducting wire carrying a current of 0.6I c, taken at different instants of a sinusoidal signal with an amplitude of 0.9I c. states that, in any (macroscopic) part of a hard superconductor, one can find either no electrical current or a current with density J equal to the so-called critical current density J c. In the original formulation, J c is constant, and it fully characterizes the properties of the material in the processes of magnetic field variation. In addition, it was formulated for bodies of high symmetry such as infinitely long cylinders and slabs. Its value is controlled by the magnetic history, according to these two general principles. No current flows in the regions that are not previously penetrated by the flux vortices. In the rest of the superconductor, the flux vortices arrange with a gradient of density that could be of a different direction but always the same magnitude. Since the penetration of vortices is accompanied by a change of local magnetic flux density B that is directly proportional to the density of vortices, the mathematical formulation of the principle of the critical state with constant J c is as follows: { 0, in regions where B =0 J = (1) J c, elsewhere. One should note that the direction of the current density is not defined by this formula. Fortunately, many problems of practical importance can be simplified to a 2-D formulation, e.g., in devices made of straight superconducting tapes or wires, which can be considered infinitely long. The advantage of a 2-D formulation of the critical state is that the current density is always parallel to the longitudinal direction. Therefore, here, we assume that the task is to find a 2-D distribution. We now illustrate the use of (1) for the problem of a round wire carrying a sinusoidal transport current with amplitude of 0.9I c. The wire s critical current I c is simply obtained by multiplying the critical current density J c by the wire s cross section. Fig. 2 shows the current density distributions when the transport current equals 0.6I c, which are taken at three different instants of the sinusoidal signal. The significant difference of the current density distribution is due to the history of the arrangement of the flux vortices in the superconductor during the sinusoidal signal. The solutions shown in Fig. 2 are derived by utilizing (1) and by considering that, due to the circular symmetry and according

3 to Maxwell equations, in the region with B =0also J =0,(no current can flow without generating a magnetic field). In his subsequent work [3], Bean transformed his formulation using the fact that any change in the magnetic field is linked to the appearance of electric field E. Then, he stated that, in a hard superconductor, there exists a limit in the macroscopic current density that it can carry, and any electromotive force, however small, will induce this current to flow. This formulation of critical state can be written as { 0, in regions where E =0in all past history J = (2) J c, elsewhere. Interestingly, the most general J(E) relation of the CSM for infinitely long conductors is [4] { J Jc, if E =0 (3) J = J c, otherwise. For conductors of finite thickness, this equation is equivalent to (2) [39]. The generalization given is essential for thin strips, which present regions where 0 < J <J c. In addition, it is also useful for numerical purposes. Establishing the correspondence between the current density and the electric field is a standard technique in electromagnetic calculations. On the other hand, (2) is multivalued because, for J = J c,anyvalueofe is possible. The value of E is determined by the electromagnetic history of the whole sample. This hinders the direct use of this E(J) relation in calculations. Nevertheless, with this approach, the principal analytical formulas for the estimation of ac losses in LTS devices have been derived [5] [9], and early numerical techniques have been proposed [10]. The CSM still represents the first approximation that allows predicting important features of HTS materials and devices, such as current density and magnetic field profiles. Due to the rapid development of computer technology in both computing power and price affordability, we can now include various refinements of the E(J) relation, which in certain cases are quite important. These include the following: 1) thermal activation leading to substantial flux creep [11]; 2) dependence of J c on the amplitude and orientation of the magnetic field; 3) spatial variation of J c (e.g., due to nonuniformities related to the manufacturing processes). The first item of this list is of fundamental importance because it provides a new formulation of the model for the electromagnetic properties of hard superconductors. The superconductor s behavior described by 1 or 2 is independent of the time derivative of the electromagnetic quantities; this, for example, means that any part of the time axis shown in the upper part of Fig. 2 could be stretched or compressed without influencing the resulting distributions. Only the existence of an electric field E and not its magnitude matters. On the other hand, thermal activation is a process happening on a characteristic time scale. Thus, if thermal activation is taken into account, the current density distribution depends on the rate of change in the case that is illustrated in Fig. 2 on the frequency and on the shape of the waveform of the transport current. Fig. 3. Typical E(J) characteristics for HTS (thin line). Differently from the CSM [thick line, (3)] there is always de/dj > 0 (not visible at low currents because of the used linear scale) caused by the significant flux creep. The CSM is a good approximation in the shaded area. It does not consider the flux-flow regime that would prevail at electric fields much higher than E c. Experimental observations of strongly nonlinear current voltage characteristics of hard superconductors (e.g., [12]) led to a formulation where the link between the current density and the electric field is expressed as J = J c ( E E c ) 1/n (4) where E c is the characteristic electric field (which is usually set equal to 10 4 V/m) that defines the current density J c, and the power exponent n characterizes the steepness of the current voltage curve. The asymptotic behavior for n in practice leads to the same dependence between the current density and the electric field as in (2). Nevertheless, there is one substantial difference: In (4), J depends on the actual value of E at the same given instant. This feature enables the incorporation of hard superconductors in electromagnetic calculations by considering it an electrically conductive and nonmagnetic material with conductivity that depends on the electric field, i.e., σ SC = J E = J c Ec 1/n E 1 n n. (5) This description is valid until the current density reaches the level causing a microscopic driving force on the vortices, which overcomes the pinning force. Then, the movement of vortices enters the so-called flux-flow regime, which can be defined by differential resistivity as [12] ρ FF = ΔE ΔJ ρ B n (6) B c2 where ρ n is the normal state resistivity, and B c2 is the upper critical magnetic field. The different regimes and the interval where the CSM is applicable are schematically illustrated in Fig. 3. In the calculation of critical current or ac losses in HTS materials, it is desirable to take also into account the dependence of the critical current density on the magnetic field and its orientation. In the case of Bi-2223 multifilamentary tapes, the use of a modified Kim s formula [12] was quite successful in reproducing experimental data [13]. In the so-called elliptical

4 GRILLI et al.: COMPUTATION OF LOSSES IN HTS Fig. 4. Simplified model describing the nonuniformity of the current density across the tape s width (in the direction of the coordinate x) by the piecewise linear profile of J c(x) defined by (8). The tape width is w s,andtheprofileis symmetric with respect to the tape center x =0. approximation, J c depends on the components of the magnetic field parallel and perpendicular to the flat face of the tape by means of the anisotropy parameter γ (with γ>1), i.e., J c (B,B )= ( 1+ J c0 B 2 +B2 /γ2 ) β. (7) B 0 The cutoff value at low fields and the rapidity of the J c reduction are given by B 0 and β, respectively. An alternative way to express the asymmetric scaling of J c with respect to the perpendicular and parallel magnetic field components is by introducing the angular dependence of the characteristic field B 0 [14]. With the introduction of artificial pinning centers in the superconductor, the simple description given by the four-parameter expression (7) is no longer adequate; more elaborated expressions are necessary to describe the experimentally observed behavior of the samples [15]. The issue of possible variations of the superconductor s critical current density with respect to spatial coordinates is quite important in the case of HTSs. The longitudinal uniformity is assessed by the determination of critical current variations along the conductor length. Since the loss at full penetration is proportional to the critical current, the change of ac losses due to longitudinal nonuniformity is only marginal for commercial tapes, which typically present a fluctuation of the critical current below 10%. On the contrary, the nonuniformity across the tape s width, which is investigated in detail for Bi-2223 multifilamentary tapes in [16] and [17], has been proposed as explanation for anomalous loss behavior of coated conductor tapes [18] [20]. This effect is of particular importance for lowloss configurations, such as bifilar coils used in resistive fault current limiters. Due to the manufacturing process, the critical current density is lower at the tape s edges than in the center, where it is usually quite uniform. As the first approximation is to take this nonuniformity into account, a symmetric piecewise linear profile, as shown in Fig. 4, can be used. Its mathematical form is J c (x) =J c0 { 1, if x hx ws 2 1 (1 h y ) 2 x ws h x 1 h x, if x>h x ws 2. (8) B. How to Solve the Electromagnetic Quantities The main step to calculate the ac losses in superconductors is to solve the electromagnetic state variable. Depending on the simulation method, the state variable can be current density J(t, r), magnetic field H(t, r), the pair composed of the current vector potential and the magnetic scalar potential T(t, r) Ω(t, r) [21], [22], or vector potential A(t, r), where t and r are the time and the position vector, respectively. Additionally, scalar potential φ(t, r) may be an extra variable. For formulations using the vector potential, the gauge is usually the Coulomb gauge (see Appendix B). All the methods reviewed in this paper assume that the displacement current is negligible, i.e., J D/ t, where D is the displacement vector. Once the state variable is solved for the time-varying excitation, the power loss can be calculated by using the methods in Section II-C. The time-varying excitation can be cyclic or not, such as a ramp increase in dc magnets. It can be transport current I(t), applied field H a (t, r) (which can be nonuniform and with an orientation varying in time; see Section II-D for particular solutions), or a combination of both. 1) Cross-Sectional Methods: In the following, we outline the methods that solve the state variable for an arbitrary combination of the current and the applied field. Other methods optimized for particular situations and their results are summarized in Section II-D. These methods solve the cross section of infinitely long conductors and bodies with cylindrical symmetry; hence, they are mathematically 2-D problems. Simulations for bidimensional surfaces and 3-D bodies are treated in Sections II-B.3 and II-B.4, respectively. A superconductor with a smooth E(J) relation can be solved by calculating different electromagnetic quantities. Brandt s method [23], [24], which is generalized for simultaneous currents and fields by Rhyner and Yazawa et al. [25], [26], solves J and φ in the superconductor volume. Alternatively, finiteelement method (FEM) models solve H [20], [27] [29], T, and Ω [30] [33], or A and φ [34] [36] in a finite volume containing the superconductor and the surrounding air, forming the simulation volume. All FEM models require setting the boundary conditions of the state variable on the boundary of the simulation volume. For most of the methods, the boundary conditions are for asymptotic values; hence, the simulation volume is much larger than the superconducting one. However, for the implementation of Kajikawa et al., the boundary conditions only require to contain the superconductor volume [28], thus reducing the simulation volume and computation time. In this sense, Brandt s method is also advantageous because it only simulates the superconductor volume. On the other hand, FEM simulations are suitable for commercial software, which simplifies the implementation and analysis. The computing time for all the methods aforementioned dramatically increases for a power-law E(J) relation with a high exponent. Variational methods may also be applied to solve problems with a smooth E(J) relation [4], [37] [39], although they are mostly used for the CSM. Critical-state calculation models are ideal to simulate superconductors with a large exponent in the power-law E(J) relation. Moreover, they are also usually faster than the simulations with a smooth E(J) relation for relatively low n. This improvement in speed may justify the sacrifice in accuracy caused by using the critical-state approximation. However, the CSM cannot describe relaxation effects or overcurrent situations.

5 Most of the existing methods for the general current and applied field imply variational methods. They were first proposed by Bossavit [4], although their most important contribution is from Prigozhin, who developed the J formulation [38], [40]. Later, Badia and Lopez proposed the H formulation [41]. An alternative numerical implementation to minimize the functional in the J formulation and to set the current constraints is the minimum magnetic energy variation (MMEV) (this method sets the current constraints directly, not through the electrostatic potential). The general formulation is described in [42] and [43], although it was first introduced in [44] for magnetization cases. The early works on MMEV implicitly assume that the current fronts penetrate monotonically in a half-cycle. This occurs in some practical cases but not in general cases [42]. In addition, FEM models with the A formulation can solve the critical state situation, as shown by Gömöry et al. [45], [46], by means of the A variation method. This technique is inspired by Campbell s A formulation for superconductors with both reversible and irreversible contributions to pinning [47], [48]. Another method based on the A formulation was developed by Barnes et al. [49], who merged the FEM technique with the critical-state constraint on J, J J c. Methods assuming a smooth E(J) relation can also approximate the CSM by considering an E(J) relation with piecewise linear segments or a power law with a high exponent [28], [37], [50]. All the given methods, both for a smooth E(J) relation and the CSM, can take the magnetic field dependence of the critical current density into account and, in principle, position dependence. 2) On the A φ and T Ω Formulations: For completeness, we next discuss the meaning of the scalar and vector potentials for the A φ and T Ω formulations. Further details can be found in a dedicated paper in this issue [51]. First, we consider the A φ formulation and, later, the T Ω formulation. In the Coulomb gauge, scalar potential φ is the electrostatic potential created by the electrical charges (see Appendix B). In practice, the scalar potential needs to be taken into account in the A φ formulation in the following cases. The first case is when there is a net transport current. In the Coulomb gauge, the current distribution creates a nonzero A/ t at the current-free core (see Fig. 2 for a typical current distribution), whereas the electric field E = A/ t φ vanishes. Therefore, it is necessary to have a certain φ in order to compensate for A/ t. One may include φ in A as a gauge, i.e., A = A + φ, where A is in the Coulomb gauge. As a result, A =0at the current-free core. However, φ has still to be taken into account in the boundary conditions far away from the superconductor. Explicitly, A φ and A μ 0 I ln r/(2π)+ φ for 3-D and infinitely long 2-D geometries, respectively, whereas in the Coulomb gauge, the boundary conditions are A =0and A μ 0 I ln r/(2π), respectively (r is the radial distance from the superconductor electrical center). In this case, φ is the electrostatic potential generated by the current source in order to keep a constant current I. The physical source of φ is the surface electrical charges at both ends of the superconductor wire and on the lateral surfaces. Second, φ is necessary when a multifilamentary body with isolated filaments is submitted to a changing applied magnetic field. In order to ensure a zero net current in all filaments, scalar potential φ between the filaments is needed. The reason is as follows. The electric field in the current-free zone is zero (see Fig. 10 for a typical current distribution in this case). In order to achieve zero electric field there, it is necessary to have φ because A/ t is nonzero at the current-free zone. The cause is that there is a net flux between the filaments, which is equal to A dr. A similar situation appears when the filaments are connected by a normal conducting material. Finally, the scalar potential is necessary for the general 3-D case and for 2-D thin surfaces. A general body submitted to a magnetic field created by a round winding produces A with cylindrical symmetry (in the Coulomb gauge). Then, if the superconducting body does not match this cylindrical symmetry, e.g., because it has corners, a certain φ that corrects the direction of E is needed, in order to allow currents parallel to the superconductor surface [48]. This is clear when E and J are parallel, but a certain φ should also be present for an arbitrary tensor relation between E and J. The meaning of the T and Ω potentials is as follows. The T potential is such that J = T. Although it follows that J = H, T is not necessarily H. Since (T H) = 0, the quantity T H can be written as a gradient of a scalar function, i.e., T H = Ω and H = T Ω. Note that T is subjected to gauge invariance. That is, T = T ϕ generates the same current density J = T = T, where ϕ is any scalar function. As a consequence, the meaning of Ω depends on the gauge of T. When the gauge function ϕ is zero, H = 2 Ω; hence, Ω is the magnetic scalar potential due to the magnetic pole density. For this gauge, T is the magnetic field created by the currents. In addition, other gauges are also used in computations [21], [22]. 3) Methods for 2-D Surfaces: An obviously mathematically 2-D geometry is thin films with finite length. The advantage of flat 2-D surfaces is that the current density can be written as a function of a scalar field, i.e., J(x, y) = ẑg(x, y), where the z-direction is perpendicular to the surface [52]. g is the z component of the current potential T since J = T. It could be also regarded as an effective density of magnetic dipoles M. For a smooth E(J) relation, this quantity can be computed by numerically inverting the integral equation for the vector potential in (57) [52] or the Biot Savart law [53]. Naturally, FEM models based on the T formulation are useful for these surfaces [54], which is the same with those based on the H formulation [55], [56]. Variational principles can solve the situation of the CSM [57], [58]. A more complex situation is when the surfaces are bent in three dimensions, such as twisted wires, spiral cables, and Roebel cables (this latter with important simplifications). The models developed so far use the FEM with the T [54], [59], [60] and H formulations [56]. 4) Methods for Fully 3-D Problems: We define a fully 3-D situation as the case when the body of study is a 3-D domain, and the problem cannot be reduced to a mathematically 2-D domain. Thus, we exclude problems with cylindrical symmetry, as well as thin surfaces with 3-D bending (see Section II-B3).

6 GRILLI et al.: COMPUTATION OF LOSSES IN HTS For a fully 3-D situation (and flat surfaces in 3-D bending), the magnetic field is not perpendicular to the current density in certain parts of the sample. Then, the modeling should take flux cutting [61] into account. This can be done by using the dependence of J c on the angle between the electric field, i.e., E, and the current density, i.e., J [41], [62], [63]. Experimental evidence reveals elliptical dependence [64]. The implementation of 3-D models can be done as follows. In principle, variational methods can describe the CSM in 3-D [4], although this has not yet been brought into practice. FEM models have successfully calculated several 3-D situations, although in all cases, the model assumes that J is parallel to E and a smooth E(J) relation. These calculations are for the T Ω [32], A φ [32], [65], [48], or H [55], [56] formulations. An alternative model based on the H formulation that does not use the Garlekin method for the FEM can be found in [66]. References for examples of 3-D calculations are in Section II-D12. C. How to Calculate AC Losses Once the Electromagnetic Quantities are Known This section deduces several formulas to calculate ac losses in superconductors once the electromagnetic quantities are known. These can be obtained by any method described in Section II-B. First, Section II-C1 justifies that the density of local dissipation is p = J E. Section II-C2 outlines how to calculate the ac losses for the CSM. Section II-C3 presents simplified formulas for the magnetization ac losses. Finally, Section II-C4 discusses how to calculate the ac losses from the energy delivered by the power sources. 1) Local Instantaneous Dissipation in Type-II Superconductors: The local instantaneous power dissipation p in a type-ii superconductor is p = J E. (9) In the following, we summarize the mechanisms that lead to this relation [12], [67], [68]. The driving force per unit length F d on a vortex is F d = J Φ 0, where Φ 0 is the vortex flux quantum, with value π /e, and the direction of Φ 0 follows the direction of the vortex. Then, the local rate of work per unit volume p on the vortices is p =(J Φ 0 n) v =(Φ 0 n v) J (10) where n and v are the vortex density and velocity, respectively. From electromagnetic analysis, it can be shown that the moving vortices with speed v create an electric field (see Appendix A), i.e., E = Φ 0 n v = B v. (11) The given fields E and B are the average in a volume containing several vortices. Finally, by inserting (11) into (10), we obtain (9). Equation (9) also applies to normal conductors [6, pp. 258]. Then, this equation is also valid to calculate the loss due to the eddy and coupling currents. However, the loss mechanism in normal conductors and superconductors is very different. In normal conductors, E creates work on the charge carriers and induces movement in them, and this movement manifests in J. In contrast, in superconductors, J is the quantity that exerts work and induces movement on the vortices. This vortex movement manifests in E. For a given local J, E does not change in time, which means that v is constant. This is because several dissipation mechanisms cause a viscous flow of the vortices [12], [68]. For power applications, the relation between E and J does not depend on the speed of the change of the current density. Kötzler et al. experimentally showed that the response to ac fields for a wide range of frequencies, i.e., 3 mhz to 50 MHz, is consistent with a unique E(J) relation of the studied sample [68], [70]. 2) Application to the CSM: The application of (9) for the local loss dissipation requires to know J and E. Once the state variable is known, the current density can be found by means of one of the following relations: H = J, T = J, and ( A) =J. Then, by assuming a smooth E(J) relation for the superconductor, E is simply obtained from that relation. For the CSM, the relation between E and J is multivalued; hence, E cannot be found in this way. For the CSM, E is calculated from E = A/ t φ. (12) The vector potential in the Coulomb gauge is obtained from J by using (57). For this gauge, φ becomes the electrostatic potential (see Appendix B). The electrostatic potential can be solved as a variable during the computation of J [38], [39] or can be calculated from A from the following relation [42], [46]: φ(t) = A(t, r 0) t (13) where r 0 is a position in the superconductor where J < J c ; hence, E =0 [see (3)]. In long conductors and bodies with cylindrical symmetry, the single component of J, i.e., J, vanishes at the current-free core and at the current fronts (boundaries between J =+J c and J = J c ). In the thin-film approximation, J <J c in the under-critical region. It should be noted that the time derivative in (13) is the partial time derivative; thus, the time dependence of r 0 should not be taken into account in evaluating the time derivative. From (9) and (13), the instantaneous power loss density is [ A(t, r0 ) p(t, r) =J t A(t, r) t ]. (14) The given equation can be simplified for monotonic penetration of the current fronts, constant J c, and infinitely long conductors or cylindrical symmetry. The deduction in the following follows the concepts from [5], [71], and [42]. Next, we present a deduction assuming any combination of a transport current and a uniform applied field proportional to each other. With monotonic penetration of the current fronts, the current density in the reverse and returning curves in an ac cycle, i.e., J rev and

7 J ret, are related to the one in the initial curve, i.e., J in,[2],[3], [72] as J rev (u) =J in,m 2J in [(1 u)/2] J ret (u) = J in,m +2J in [(1 + u)/2] (15) where J in,m is the current distribution at the end of the initial curve, and u is defined as follows. For a transport current, u(t) is such that I(t) =u(t)i m, where I m is the amplitude of the transport current. If there is no transport current but the applied field is created by a winding with a linear current field characteristic, we define u(t) as H a (t, r) =u(t)h m (r), where H m (r) is the maximum applied field at position r. This relation is also valid for the simultaneous action of a transport current and an in-phase applied field. For monotonic current penetration, there is always at least one point, which is called a kernel, where J =0 in the whole cycle. Then, r 0 in (14) becomes constant. Integrating the loss power density of (14) in the volume and integrating by parts for the time variable, we find that Q = dt V J t [A(t, r 0) A(t, r)] d 3 r (16) where we took into account that J in each position changes once in each half-cycle and that this change is instantaneous from J c to J c (or vice versa). The time derivative of this discontinuous dependence is J(t, r)/ t = 2J c δ[t t c (r)], where δ is the Dirac delta, and t c (r) is the time of the shift of J at position r. The time integration results in Q =4 J in,m (r) {A [t c (r), r 0 ] A [t c (r), r]} d 3 r. (17) V The relation of (15) is also valid for the vector potential A in the Coulomb gauge. Then A(r,t)=±A in,m (r) 2A in [r, (1 u(t)) /2] (18) where the top sign is for the reverse curve, and the bottom sign is for the returning curve. The current density shifts its sign at the current front (see Fig. 2). For the initial curve, the current front encloses the flux-free region, where B = J =0; hence, A is uniform. Then, A at the current front is the same as in the kernel. By taking this into account and using (18), the loss per cycle from (17) becomes Q =4 J in,m (r)[a in,m (r 0 ) A in,m (r)] d 3 r. (19) V For a pure transport current, the given equation turns into Q =4J c [A in,m (r 0 ) A in,m (r)] d 3 r (20) V because, at the zone where J in,m vanishes, A in,m (r) = A in,m (r 0 ). Equation (20) corresponds to the original formula from Norris [5]. The general equation (19) reduces to the one given by Rhyner for a uniform applied field [71]. 3) Magnetization Loss: When the only time-dependent excitation is the applied magnetic field H a (t, r), the expression for the ac losses per cycle can be simplified. Next, we take into account the ac losses in the superconductor only. In the following, we deduce several formulas for the magnetization loss. From (9), the density of the instantaneous power dissipation in the whole superconductor is p = J s E (21) where J s is the current distribution in the superconductor. Then, the instantaneous power dissipation is P = V J s E d 3 r, where V is any volume containing the superconductor, and d 3 r is the volume differential. The loss per cycle Q is the time integration in one cycle. Using H s = J s, E = B/ t, the differential vector relation ( H s ) E = (H s E)+( E) H s, and the divergence theorem, we find that P = V H s B t d3 r + V H s E ds (22) where H s is the magnetic field created by the superconductor currents, V is the surface of the volume V, and ds is the surface differential. When the volume V contains not only the superconductor but also the sources of the applied magnetic field (the current or magnetic pole densities) and V approaches infinity, the second term of (22) vanishes. This can be seen as follows. Let us take V as a sphere of radius r. With large r, the largest contribution to the fields is their dipolar component. Then, H s 1/r 3, and since E = A/ t φwith A 1/r 2 and φ 1/r 3, it follows that H s E 1/r 5. Since the surface only increases as r 2, the integration vanishes with r. Then, the instantaneous power loss is P = H s V B t d3 r (23) where V is the whole space. Note that (23) for the total power loss is general, although its integrand is no longer the local rate of power dissipation. The ac losses per cycle Q as a response to a periodic ac field can be further simplified. If the whole system is in the void, B = μ 0 H = μ 0 (H a + H s ), and H s B/ t = μ 0 [H s H a / t +(1/2) H s 2 / t]. When integrating in one cycle, the second term trivially vanishes, resulting in Q = μ 0 dt H s V H a t d 3 r. (24) Note that the applied field H a (t, r) in the given equation is not necessarily uniform, and it may rotate during the cycle. In addition, the volume integral in (24) is not the instantaneous power loss. For uniform H a, we can use the fact that V H s d 3 r = (1/2) V r J s d 3 r = m; hence, (24) becomes Q = μ 0 m H a dt (25) t

8 GRILLI et al.: COMPUTATION OF LOSSES IN HTS where H a may still rotate in one cycle. For an applied field in a constant direction, the ac losses per cycle are Q = μ 0 m z dh a (26) where we choose the direction of the z-axis as the direction of the applied field. Next, we rewrite (24) in a more useful way for calculations using the J or A formulations. Using μ 0 H s = A s,differential vector equality ( A s ) H a = (A s H a )+ ( H s ) A a, and the divergence theorem, we find that Q = dt J s V A a t d 3 r + V A s H a t ds. (27) Using the same arguments as for (22), the term with the surface integral vanishes. Note that, even if the applied field may be uniform in the superconductor, far away from its sources, it has to decay with the distance. Then, (27) becomes Q = dt J s V A a t d 3 r. (28) For magnetic materials originated by a microscopic density of magnetic dipoles M, it is practical to use formulas independent of the loss mechanism associated to the change of the local density of dipoles because they can be applied to any magnetic material. The variation of the free energy of a magnetic material is δf = V M δh a [6, pp. 116]. The loss per cycle is the integration of the free energy in a cycle, i.e., Q = μ 0 dt V M H a t d 3 r. (29) From this formula, it is trivial to reproduce (25) and (26) for uniform applied fields. This formula can be also written in terms of the total magnetic field H, i.e., Q = μ 0 dt t V M H/ t. For soft magnetic materials, M is practically parallel to H. By taking this into account and neglecting the loss due to a rotation in the magnetic field, the loss per cycle is Q μ 0 V MdH = μ 0 V q(h m,h bias ) d 3 r. (30) The given equation assumes that the magnitude of the local magnetic field oscillates with amplitude H m around H bias. The quantity q(h m,h bias ) is the density of hysteresis loss corresponding to this amplitude, which can be experimentally obtained. If the applied magnetic field oscillates with no dc component, H bias =0for the whole volume. The density of hysteresis loss can be also written as a function of the magnetic flux density, i.e., B m and B bias. For infinitely long wires and tapes, the power loss and the loss per cycle per unit length, i.e., P l and Q l, follow (23), (24), (28) (30), after the replacement of the volume integration with the cross-sectional surface integration. The same applies to (25) and (26) but with the replacement of the magnetic moment m with the magnetic moment per unit length m l. 4) AC Losses From the Point of View of the Power Source: By continuity of the energy flow, the loss per cycle dissipated in the superconductor equals the energy per cycle supplied by the current sources. In general, a superconductor may be under the effect of a transport current or a magnetic field created by an external winding, which implies two current sources [74], [75]. If the model solves a single conductor in pure transport current or a complete superconducting coil, there is only one power source. The voltage delivered by the power source corresponds to the difference of the electrostatic potential Δφ. Then, the ac losses per cycle in the coil are Q = IΔφ dt. (31) For the CSM, the drop of the electrostatic potential in each turn can be calculated from the scalar potential from (13) and by multiplying by the length of the conductor [76], [43]. The total voltage drop is the sum of the drop in all the turns. Although (31) for the loss per cycle is valid, the instantaneous power dissipation in the coil not necessarily equals the instantaneous power delivered by the source, i.e., I(t)Δφ(t). The cause is the (nonlinear) inductance of the coil. This can be seen as follows. From (12) and (21), the instantaneous power ac losses are ( P (t) = J s φ + J s A ) d 3 r. (32) t V Since there are no external sources of magnetic field, A in the Coulomb gauge only depends on J s through (57). As a consequence, V J s A/ t d 3 r = V J s/ t A d 3 r. Then, V J s A/ t d 3 r is the time derivative of the magnetic energy U =(1/2) V J s A d 3 r. Thus, the time integral in one cycle of the second term in (32) vanishes. The remaining term is (31) because, in straight conductors, φ is uniform [77], as well as in each turn of circular coils [76]. If there is also an applied magnetic field, the total loss per cycle from (12) and (21) becomes Q = IΔφ dt dt J s V A a t d 3 r (33) where we follow the same steps as the reasoning after (32). The second term of (33) is the magnetization loss from (28). That formula is equivalent to (24), which reduces to (25) for uniform applied fields in the superconductor. Under this condition, the total loss in the superconductor is Q = IΔφ dt μ 0 m H a t dt. (34) The first term is usually called transport loss, whereas the second term is called magnetization loss. If the applied magnetic field and the transport current are in phase, the transport and magnetization losses correspond to the loss covered by the transport and magnetization sources, respectively. This is because the sources are only coupled inductively; hence, the energy transfer from one source to the other in one cycle vanishes. Equation (34) is universal; thus, it is also valid when

9 the sources are not in phase to each other. However, in that case, the two terms in (34) cannot be always attributed to the loss covered by each source separately [75]. D. Solutions for Particular Cases The purpose of this section is to gather a reference list of the hysteresis loss classified by topics. The number of published papers in the field is vast. Although we made the effort of making justice to the bulk of published works, a certain degree of omission is inevitable, and we apologize in advance to the concerned authors. Readers interested in a particular topic may go directly to the section of interest according to its title. In the following, we do not regard the case of [78]. A review on that problem can be found in another paper of this issue [78]. 1) AC Applied Magnetic Field: The response to an applied magnetic field strongly depends on the geometry of the superconducting body, as discussed in [8] and [9]. The first works were for geometries that are infinitely long in the field direction because there are no demagnetizing effects. Bean introduced the CSM for slabs [1] and cylinders [3], and calculated the magnetization and the ac losses in certain limits. Later, Goldfarb et al. and Clem obtained formulas for the ac susceptibility for any field amplitude [79], [80]. The formula for the ac losses per unit volume is Q v = 2μ 0H 3 m 3H p if H m <H p = μ 0H p (6H m 4H p ) 3 if H m H p 2μ 0 H p H m if H m H p. (35) In the given equation, H m is the amplitude of the applied magnetic field, and H p = J c w/2, where w is the thickness of the slab. To obtain this equation, we used the fact that the ac losses per unit volume Q v are related to the imaginary part of the ac susceptibility χ as Q v = μ 0 πhmχ 2. The imaginary part of the ac susceptibility (and the loss factor Q/Hm)presentsa 2 peak at an amplitude, i.e., H peak = 2 3 J cw. (36) This formula is useful to obtain J c from magnetization ac loss measurements. There are extensive analytical results for the critical state with magnetic field dependence of the critical current density J c (B) for infinite geometries. Most papers only calculate the current distribution and magnetization in slabs [81], [82]. The most remarkable analytical work for the ac losses is from Chen and Sanchez [83], who studied rectangular bars. Comparisons between the CSM and a power-law E(J) relation are in [84] for a slab and in [85] for a cylinder. The latter paper also finds that the frequency for which the ac susceptibility from the CSM and a power-law E(J) relation is the most similar, at least for field amplitudes above that of the peak in χ. This frequency is f c = E c /(2πμ 0 a 2 J c ) (37) where a is the radius of the cylinder. This frequency is based on an extension of the scaling law by Brandt [86]. This scaling law applied to χ reads ( ) χ PL(H m,f)=χ PL H m C 1 n 1,Cf (38) where χ PL is χ for a power-law E(J) relation, C is any arbitrary constant, and n is the exponent of the power law. Combining the scaling law in (38) with (37), one can find a relation between χ for the CSM, i.e., χ CSM, and that for the power law, i.e., χ PL(H m,f)=χ PL [ χ CSM ( fc H m f [ ( fc H m f ) 1 n 1,fc ] ) ] 1 n 1. (39) The first geometries with demagnetizing effects to be solved were thin strips and cylinders. For these geometries, conformal mapping techniques allowed to find analytical solutions for the CSM. Halse and Brandt et al. obtained the case of a thin strip in a perpendicular applied field [6], [72], although an earlier work from Norris developed the conformal mapping technique for a transport current [5]. The ac losses per unit volume are Q v = μ 0 wj c H m ( 2 x ln cosh x tanh x ) μ 0 wj c H m if H m H c (40) where x is defined as x = H m /H c with H c = J c d/π, and d and w are the thickness and width of the strip, respectively. For the same dimension w and low H m, the ac losses per unit volume in a slab are smaller than in a thin film. However, the limit of high applied field amplitudes is the same for the strip and the slab, i.e., Q v μ 0 wj c H m [from (35) and (40)]. The peak of χ (and Q/H 2 m)isatanamplitude[87] H peak =0.7845J c d. (41) The ac losses in a thin disk with constant J c were obtained by Clem and Sanchez [88]. The solutions of the current and field distributions (as well as the magnetization) are analytical, although the expression for ac losses is in integral form. In general, to solve the geometries with intermediate thickness requires numerical techniques. In the following, we first summarize the works on geometries for wires in a transverse applied field (rectangular, circular, elliptical, and tubular); then, the works on several of these wires parallel to each other (multiple conductors); afterward, the works on the case of wires under an arbitrary angle with the surface of the wires; and finally, the works on the solutions for several cases with rotational symmetry subjected to a magnetic field in the axial direction (cylinders, rings, spheres, spheroids, etc.). Rectangular wires were studied in detail by Brandt, who calculated the current distribution and other electromagnetic quantities for a power-law E(J) relation [23]. For the CSM, Prigozhin computed the current distribution [38], and

10 GRILLI et al.: COMPUTATION OF LOSSES IN HTS Pardo et al. discussed the ac susceptibility and presented tables in [87]. Circular and elliptical wires are geometries often met in practice. The main papers for the CSM are the following. The earliest work is from Ashkin, who numerically calculated the current fronts and the ac losses in round wires [10]. Later, several authors calculated the elliptical wire by analytical approximations, either by assuming elliptical flux fronts [8], [89], [90] or very low thickness [91]. Numerical solutions beyond these assumptions are in [40], [90], and [92], with the latter containing tables for ac susceptibility. Earlier work already solved this situation for a power-law E(J) relation [93], [94]. The most complete work on tubular wires is from Mawatari [95]. This paper presents analytical solutions for the current and flux distribution, as well as the ac losses, in a thin tubular wire. The current distribution in a tube with finite thickness had previously been published in [40]. For the hysteresis loss in multiple conductors, it is necessary to distinguish between the uncoupled and fully coupled cases. These cases are the following. If several superconducting conductors (wires or tapes) are connected to each other by a normal conductor, either at the ends or along their whole length, coupling currents appear, and as a consequence, the ac losses per cycle depend on the frequency (see Section IV). The ac losses become hysteretic not only for the limit of high interconductor resistance (or low frequency) but also for low resistance (or high frequency), as shown by measurements [96] and simulations [54]. These limits are the uncoupled and fully coupled cases that we have introduced earlier. It is possible to simulate these two cases by setting different current constraints: For the uncoupled case, the net current is zero in all the conductors, whereas the fully coupled one allows the current to freely distribute among all the tapes [97]. For thin strips, there exist analytical solutions. Mawatari found the ac losses in vertical and horizontal arrays of an infinite number of strips [98] for the uncoupled case, and Ainbinder and Maksimova solved two parallel strips for both the uncoupled and fully coupled cases in [99]. The study of an intermediate number of tapes requires numerical calculations. Simulations for vertical arrays can be found in [100] and [101] for a powerlaw E(J) relation and in [97] for the CSM. Complicated cross sections for the fully coupled case can be found in [102]. Matrix arrays for both the coupled and uncoupled cases are in [97] and [103]. A wire with uncoupled filaments is in [104]. A systematic study of the initial susceptibility in multifilamentary tapes with uncoupled filaments is in [105]. All the given results for single and multiple conductors are for applied magnetic fields parallel to either the thin or wide dimension of the conductor cross section. Solutions for intermediate angles can be found in the following papers. Mikitik et al. analytically studied a thin strip, but they did not calculate the ac losses [106], [107]. The ac losses in single conductors are in [25] and [108] [110] for a power-law E(J) relation and in [46] and [111] for the CSM. Solutions for arrays of strips for both a power-law E(J) relation and the CSM are in [112]. The latter reference also reports that the ac losses in a thin strip with magnetic field dependence of J c may depend on both the parallel and perpendicular components of the applied field; therefore, the ac loss analysis cannot be always reduced to only the variation on the perpendicular component of the applied field. Next, we list the works for bodies with rotational symmetry. Apart from the analytical solutions in [88] and [113] for thin disks, the main works for cylinders subjected to a uniform applied field are the following. Brandt studied in detail several electromagnetic properties of cylinders with finite thickness in [24] and the ac susceptibility in [114] for a power-law E(J) relation. Later, Sachez and Navau studied the same geometry for the CSM [44], and Chen et al. discussed the ac susceptibility and publish tables [92]. For cylindrical rings, the ac susceptibility and ac losses were calculated in [86] for thin rings and, in [115], for finite thickness. Pioneering works on the current distribution and ac losses in spheres and spheroids in the CSM can be found in [116] [118]. Finally, the current distribution for several bodies of revolution can be found in [38] and [40]. All the numerical models discussed in Section II-A can solve the general situation of a conductor with an arbitrary cross section, either infinitely long or with rotational symmetry. 2) AC Transport Current: The effect of an alternating transport current in a superconducting wire (or tape) also depends on the wire geometry. The simplest situation is a round wire because of symmetry. London calculated this situation in the CSM [2]. London ideated the CSM in parallel to Bean [1]. The ac losses per unit length Q l in a round wire is 2πQ l μ 0 I 2 c =(2 i m )i m +2(1 i m )ln(1 i m ) i3 m 3 if i m 1 (42) where the expression on the left is adimensional, i m = I m /I c, and I m is the current amplitude. Norris found that this formula is also valid for wires with an elliptical cross section (elliptical wires), which is based on complex mathematical arguments [5]. He also developed the technique of conformal mapping to calculate the current distribution in thin strips [5]. The resulting ac losses in a thin strip are 2πQ l μ 0 I 2 c =2 [ (1 + i m )ln(1+i m )+(1 i m )ln(1 i m ) i 2 ] m i4 m 3 if i m 1. (43) Later, Clem et al. calculated the magnetic flux outside the strip [119]. They found the key result that it is necessary to use C- shaped loops closing at least the width of the strip in order to measure properly the ac loss. Contrary to the case of an applied magnetic field, thin films cannot be simplified as a slab with critical current density penetrating in a uniform depth across the film width. The formula for a slab from Hancox [120] can only be applied to estimate the ac losses due to the current penetration from the wide surfaces of the film, also known as top and bottom losses. Curiously, the formula for a slab exactly corresponds to the limit of low current amplitudes in (42). Numerical methods are usually necessary to calculate the ac losses for other geometries or an E(J) relation.

11 For an elliptical wire, Amemiya et al. compared the ac loss for the CSM and a power-law E(J) relation [93]. Later, Chen and Gu further discussed this kind of comparison for round wires in [121], where tables for the ac losses are presented. They also found a scaling law for the ac losses, as for the magnetization case. The scaling law for the loss per unit length Q l,pl is Q l,pl Im 2 (I m,f)= Q ( ) l,pl Im 2 I m C 1 n 1,Cf (44) where I m is the amplitude of the current, and C is an arbitrary constant. A study of the effect of the field dependence in an elliptical wire can be found in [94]. The effect of nonhomogeneity in the cross section of a round wire and a thin strip are discussed in [18] and [122], respectively. In summary, degraded superconductor at the edges of the conductor increases the ac loss and vice versa. A similar effect appears if the strips are thinner close to the edges [123]. The main results for rectangular wires with finite thickness are the following. Norris was the first to numerically calculate the current fronts and ac losses for the CSM [124]. Later, several authors published more complete works [87], [125], [126], where [87] also presents tables and a fitting formula. Results for a power-law E(J) relation can be found in [127]. The following papers calculate the current distribution and ac losses in multiple wires and tapes connected in parallel; thus, the current can freely distribute among the conductors (the case of wires connected in series corresponds to coils, which is discussed in Section II-D6). The only analytical solution is for two coplanar thin strips in the critical state [99]. Numerical results that are always in the critical state for Bi-2223 multifilamentary tapes with several geometries are in [128] and [129]. The effect of arranging rectangular tapes in horizontal, vertical, and rectangular arrays is studied in [130]. A comparison between the CSM and a power-law E(J) relation in a matrix array of coated conductors is in [103]. Valuable earlier work using a power law for several kinds of multifilamentary wires and tapes is in [102] and [131]. We review cables transporting alternating current in Sections II-D4 and II-D5. 3) Simultaneous Alternating Transport Current and Applied Field: This section outlines the results for superconducting tapes under the action of both an alternating transport current and a uniform applied field. First, we summarize the case of inphase current and field and, afterward, with an arbitrary phase shift. For a current and a field in phase with each other, analytical solutions only exist for slabs and strips in the CSM. Carr solved the situation of a slab and presented a formula for the ac losses [132]. (The reader can find the same formula in SI units in [42].) Later, Brandt and Indenbom and Zeldov et al. simultaneously obtained the magnetic field and current distributions for a thin strip [133], [134], and Schonborg deduced the ac losses from them [135]. However, as discussed in the original two papers [133], [134], the formulas are only valid when the critical region monotonically penetrates in the initial curve. This is ensured at least for the high-current low-field regime, which appears for I I c tanh(h a /H c ) with H c = J c d/π. The range Fig. 5. Schematic view of a two-layer HTS power cable and of 2-D and 3-D regions to be considered for simulation. Figure taken from [59] and reproduced with permission from Institute of Physics (IOP). of applicability of the slab and strip formulas compared with numerical calculations is discussed in [42]. In that reference, the current penetration process in a rectangular wire with finite thickness is studied in detail. However, the earliest numerical calculations are for a power-law E(J) relation [26]. Elliptical wires with a power law are computed and discussed in detail by several authors [93], [27], [35], [136]. The latter reference also compares different shapes: round, elliptical, and rectangular wires. The main features for the thin geometry of a coated conductor is discussed in [31]. In electrical machines, the current and the magnetic field are not always in phase with each other. Analytical solutions for an arbitrary phase shift only exist for slabs [137] and strips [138] in the CSM and for low applied magnetic fields. These calculations predict that the ac loss maximum is when the alternating current and the field are in phase with each other. However, simulations for slabs [139] and elliptical wires [75], [140] show that, for large applied fields, the ac loss maximum is at intermediate phase shifts for both a power-law E(J) relation and the CSM. Experiments confirm this fact [75], [140]. 4) Power-Transmission Cables: The cross section of stateof-the-art power-transmission cables is of superconducting tapes lying on a cylindrical former, in order to minimize the ac loss. These tapes are spirally wound on the former, allowing the cable to bend (see Fig. 5). To estimate the loss, the cross section is often approximated as a cylindrical shell (or monoblock) because its solution in the CSM is analytical [120]. However, other analytical approximations may also be applicable [141]. Most of the models for the real cross section assume that the tapes are straight, neglecting the effect of the spiral shape; this corresponds to considering just the cross section in Fig. 5(a) and to assuming that it extends straight along the cable s axis. The error of this approximation is negligible for single-layer cables. Moreover, due to symmetry, the cross section in the computations can be reduced to a single circular sector [142] [144]. Simulations for tapes with a ferromagnetic substrate show that its presence always increases the ac losses [145], [146]. In multilayer cables, the twist pitch of the tapes strongly influences the current sharing between layers [147]. One reasonable approximation is to assume that the twist pitch is optimum, so that all the tapes carry the same current independently of the layer that they belong to [148] [150]. All the given numerical calculations are for a smooth E(J) relation, although it is also possible to assume the CSM [151]. These computations show that nonuniformity in the tapes increases ac losses [151].

12 GRILLI et al.: COMPUTATION OF LOSSES IN HTS Fig. 6. Typical structure of a Roebel cable composed of meander-shaped strands obtained from ReBCO-coated conductors. Only half of the stands are shown for clarity. Analytical solutions are more recent than numerical ones, due to the complexity of the deductions. These have been found for cables made of thin strips in the CSM. They show that, if the tapes are bent and, hence, the cable cross section is a tube with slits, the ac losses are lower than if the tapes are straight, forming a polygon [152]. The real spiral geometry is more complicated to compute. The first step has been to calculate the critical current [153]. AC loss calculations for cables made of Bi-2223 tapes require 3-D models [see Fig. 5(b)], resulting in high computing complexity, so that only coarse meshes are considered [32], [154]. More recent models for coated conductors improve the accuracy by reducing the problem to a mathematically 1-D [155] or 2-D [55], [59] one. 5) Roebel Cables: Roebel cables are compact transposed cables for windings that require high currents (see Fig. 6). Due to transposition, coupling loss is minimized. The complicated structure of Roebel cables can be simplified with a 2-D approximation. If the transposition length is much larger than the cable width, the loss in the crossed strands is much smaller than in the straight parts. As a consequence, the cable can be modeled as a matrix of tapes (see [103] and [156] for more details). For an applied magnetic field, the computations distinguish between the uncoupled and fully coupled cases, as for multiple tapes (see Section II-D1), not only if the applied field is perpendicular to the strand surface [103] [156] but also for an arbitrary angle [112]. The equivalent can be also done for the transport loss, with the interpretation that the uncoupled and fully coupled cases are for transposed and untransposed strands [103], [156], [157]. The uncoupled case requires that the net current in all the strands is the same, whereas the fully coupled one allows the current to distribute freely among strands. Until now, only two papers have taken the real 3-D bending into account [60], [158]. The work by Nii et al. is based on the T formulation and assumes the strands to be infinitely thin [60]. The authors calculate the current density and ac loss distribution in the uncoupled case for one example of transport current, applied magnetic field, and a combination of both. That work shows that there appear points of high loss density close to the crossing strands. This increases the ac loss comparing to cross-sectional models. The increase is important for low ac applied fields (and no transport current) because the current penetration in the crossed strands is much larger than in the straight parts. Otherwise, cross-sectional models do not introduce an error larger than 10% concerning the total loss. Zermeno et al. used the H formulation to develop a full 3-D electromagnetic model of a Roebel cable [158]. The use of periodic boundary conditions allowed simulating a periodic cell, keeping the number of degrees of freedom at a manageable level in this way. This model also reveals high-loss-density regions near the crossing strands in the case of applied perpendicular field. In the transport current case, the losses are mostly localized along the long edges of the strands, which supports the possibility of using a 2-D model for this case. Results are also compared with experimental data, showing a fairly good agreement; the observed mismatch for the magnetization losses is probably due to different reasons: for example, the fact that the experimental setup does not grant the perfect uncoupling of the strands (assumed in the simulations) and that the model does not take into account the dependence of J c on the magnetic field. 6) Coils: The ac losses in the tape (or wire) that composes a superconducting coil is subjected to both a transport current and a magnetic field created by the rest of the turns. Since real coils are made of hundreds or thousands of turns, computations are complex. Early calculations computed only the critical current, which already provided valuable information [159]. One possible simplification for the ac losses is to approximate that the effect of the whole coil in a certain turn is the same as an applied magnetic field. This applied field is computed by assuming that the current density is uniform in the rest of the turns. Afterward, the ac losses in the turn of study are estimated by either measurements in a single tape [160], [161] or by numerical calculations [162]. The problem of this approximation is that the neighboring turns shield the magnetic field from the coil [97], [101]. This situation is the worst if the winding consists of stacks of pancake coils because the whole pancake shields the magnetic field, as it is the case for stacks of tapes [97], [101]. There are several ways to take into account the nonuniform current distribution in all the turns simultaneously. The first way is to assume infinite stacks of tapes, describing only the ac losses at the center of a single pancake coil, either by analytical calculations for the CSM [163] or by numerical ones for a smooth E(J) relation [164]. The second is to simply solve all the turns simultaneously for single pancake coils (pancakes) [39], [43], [76], [165], [166], double pancakes [167], and stacks of pancakes [112]. For this, it is possible to calculate the circular geometry [43], [76], [112], [165], [166] or assume stacks of infinitely long tapes [39], [167]. Finally, for closely packed turns, it is possible to approximate the pancake cross section as a continuous object [39], [168]. Semi-analytical approaches assume that the region with critical current is the same for all the turns [169], [170]. This condition can be relaxed by assuming a parabolic shape of the current front [168]. However, advanced numerical simulations do not require any assumption about the shape of these current fronts [39]. The presence of ferromagnetic materials strongly influences the ac losses in the superconductor. Magnetic diverters reduce

13 the ac losses in pancake coils [171], [172], whereas a ferromagnetic substrate increases it [167], [173]. 7) Noninductive Windings: Noninductive windings are the typical configurations for resistive superconducting fault current limiters (SFCLs). The reason is the low inductance and low ac losses. Modeling of this situation starts with two parallel tapes on top of each other transporting opposite current (antiparallel tapes). For thin strips, it is possible to obtain approximated analytical formulas because the ac losses are dominated by the penetration from the edges [5]. However, if the thickness of the rectangular wire is finite, numerical computations show that, for low current amplitudes, the ac losses is dominated by the flux penetration from the wide surfaces [124]. This contribution of the ac losses is called top and bottom loss [174]. For practical superconductors, the aspect ratio is high (at least 15 and 1000 for Bi-2223 and ReBCO tapes, respectively). For this cases, the ac losses for antiparallel tapes is much lower than if the tapes are placed far away from each other. Most of this ac loss reduction is lost if the tapes are not well aligned [175]. Degradation at the edges of the coated conductors also worsens the low-loss performance of the antiparallel configuration [176]. Bifilar windings (windings made of antiparallel tapes) reduce even further the ac loss if they are either made of elliptical wires [177] or coated conductors [174], [176]. The most favorable situation concerning ac losses is for bifilar pancake coils [176]. For strips in this configuration, there exist valuable analytical formulas [174]. There also exist analytical formulas for the equivalent of a bifilar coil but with all the strips on the same plane [178]. The ac losses for this configuration are larger than if the tapes are far away from each other, which is opposite to bifilar pancake coils. 8) Rotating Fields: Modeling the effect of rotating fields is very different if the magnetic field is either perpendicular to the current density or not. For infinitely long wires or tapes in a transverse applied magnetic field, the local current density is always perpendicular to the magnetic field. All the numerical methods from Section II- B1 can solve this situation. The main feature of the current distribution is that it becomes periodic only after the first halfcycle that follows from the initial curve, as seen for cylinders [40], square bars [40], and rectangular bars [37]. This feature has an effect on the magnetization [179]. If the sample has been previously magnetized, a rotating field partially demagnetizes the sample after a few cycles [179], [180]. If the magnetic field is not perpendicular to the current density in part of the ac cycle, there appear flux cutting mechanisms in addition to vortex depinning; thus, the situation is more complex. This situation can be described by a double CSM with two characteristic critical current densities: one for depinning and another one for flux cutting [181]. This model allows calculating analytical solutions of the whole electromagnetic process, including the ac losses, for a slab submitted to a rotating magnetic field parallel to its surface [182]. An extension of this model takes into account a continuous variation of the critical current density as a function of the orientation of the electric field with respect to the current density [41], [183]. This model requires numerical methods to solve the current distribution [41], [183], [184]. Experimental evidence suggests elliptical dependence between the electric field and the critical current density [64]. 9) Pulsed Applied Field: This section first summarizes the main modeling issues for bulk applications (slabs and cylinders) and, afterward, for coils in pulse mode. Calculations for slabs show that, essentially, the response to a pulsed applied field is very similar to a periodic ac field [185], [186]. If the pulse starts from a zero applied field, then increases to a certain maximum H m, and finally decreases to zero again, the response for the periodic case in the increasing and decreasing parts of the pulses is the same as the initial and reverse curves for a periodic ac field, respectively. For the CSM, the response is exactly the same, although for a smooth E(J) relation, the instantaneous current and flux distribution depends on the shape of the pulse, as seen for slabs [185]. In case that there is a bias dc field superimposed to the pulse, the movement of the current fronts is exactly the same as for an ac field, except for the first increase in current. In addition, as a consequence of the dc field, there is a region with critical current density inside the sample. This has been analytically studied in [186]. The geometry of finite cylinders is more realistic to model bulk samples, but it requires numerical computations. The process of the current and flux penetration has been studied for the CSM [187] and for a power-law E(J) relation [188]. The latter work also takes into account the self-heating effects caused by the changing magnetic field. For the design of dc magnets, it is essential to consider the dissipation in pulse mode [189], [190]. There are several ways to model the loss in pulse mode. The first one is to use the measured ac losses of the cable in a pulsed applied field as a function of the peak field. The total energy loss is calculated from computations of the magnetic field of the coil by assuming uniform current density in the superconductor and by taking the corresponding measured loss in the cable into account [189]. The second way is to assume the slab model for the wire [190], which is a good approximation for large magnetic fields. The last way is to calculate the detailed current distribution in the whole magnet, as explained in Section II-D6. 10) AC Ripple on DC Background: The problems that a researcher can encounter in the modeling of the ac losses in a dc background are very different if the dc background is a transport current or an applied magnetic field. The case of a dc applied magnetic field perpendicular to the ac excitation is outlined in Section II-D11. The situation of an ac applied magnetic field superimposed to a dc component in the same direction is basically the same as pulsed applied fields (see Section II-D9). After the increase in the applied field to the peak, the process of flux penetration is the same as a pure ac field. The only difference is that the critical current density is lower because of the presence of the dc component, as seen for slabs [186] and stacks of tapes in perpendicular applied field [191]. This concept is exploited to estimate the magnetization ac losses in complex coils due to current fluctuations around the working current. One possibility is to use analytical approximations for the ac losses in one wire [192].

14 GRILLI et al.: COMPUTATION OF LOSSES IN HTS An interesting case is when the dc applied field is in the same direction as a certain alternating transport current. For this situation, the ac transport current overwhelms the dc magnetization, approaching to the reversible value after a few cycles [193]. These calculations require advanced models because the magnetic field has a parallel component to the current density (see Section II-D8). For a wire with a dc transport current and a transverse applied ac field, the most interesting feature is the appearance of a dynamic resistance [194]. This dynamic resistance appears for amplitudes of the applied magnetic field above a certain threshold. Analytical formulas exist for slabs with either a constant J c [194] or a field-dependent one [195]. This case has also been numerically calculated for stacks of tapes, in the modeling of a pancake coil [191]. However, the applied fields are small compared with the penetration field of the stack; hence, the dynamic resistance may not be present. A dc transport current may increase the loss due to an ac oscillation, as shown by analytical approximations for slabs and strips in the CSM [196]. These formulas also take an alternating applied field into account. Extended equations including also the effect of a dc magnetic field can be found in [197], where they are applied to optimize the design of a dc powertransmission cable. 11) AC Applied Magnetic Field Perpendicular to DC Background Field: Given a superconductor magnetized by a dc magnetic field, a relatively small ac magnetic field perpendicular to the dc one strongly reduces the magnetization (collapse of magnetization), as experimentally shown in [198]. The explanation is different if the sample is slab-like, with both applied fields parallel to the surface or if either the dc field is perpendicular to the sample surface or the sample has a finite thickness. For a slab with both applied fields parallel to the surface, the collapse of the magnetization can only be explained by flux cutting mechanisms [198]. By assuming the double CSM, solutions can be analytically found for a limited amplitude of the ac field [198]. General solutions for continuous angular dependence between the critical current density and the electric field require numerical methods [62], [199]. For infinitely long wires in a transverse applied field, the current density is always perpendicular to the magnetic field independently of the angle with the surface. Then, flux cutting cannot occur. The collapse of magnetization is explained by the flux penetration process. This can be seen by physical arguments and analytical solutions for thin strips [200] and by numerical calculations for rectangular wires [201]. Collapse of magnetization is also present when the ac applied field is parallel to the infinite direction of a strip (and the dc one is perpendicular to the surface) [202]. 12) 3-D Calculations: In this paper, we consider 3-D calculations as those where the superconductor is a 3-D object that cannot be mathematically described by a 2-D domain. A summary of the requirements and existing models for the 3-D situation is in Section II-B4. The calculated examples are power-transmission cables made of tapes with finite thickness [32], [56], [154], multifilamentary twisted wires [55], two rectangular prisms connected by a normal metal [32], a permanent magnet moving on a thick superconductor [66], a finite cylinder in a transverse applied field [48], a cylinder with several holes [65], and an array of cylinders or rectangular prisms [56]. III. EDDY CURRENTS Eddy-current computations go beyond problems involving superconductivity. They are important for modeling electric machines, induction heaters, eddy-current brakes, electromagnetic launching, and biomedical apparatus, to name just a few applications [203]. Thus, we begin this section by shortly reviewing the key research done by the computational electromagnetism community. Then, we discuss the specific numerical methods to solve eddy-current problems and derive losses. Finally, we review some of the numerical loss computations related to nonsuperconducting and nonmagnetic but conducting parts present in technical superconductors. A. Eddy Currents and Conventional Electromagnetism Numerical simulations of eddy-current problems in conventional computational electromagnetism have been under investigation since the early 1970s [204], [205]. With the exception of some work done at the very beginning [206], [207], most of the early research concentrated on 3-D computations [205], [208], [209], and theoretical background for numerical eddycurrent simulations is mostly established in 3-D space [210] [213]. The need for 3-D modeling comes from practical electrotechnical devices, such as electric motors. It is necessary to model them in three dimensions to guarantee accurate simulations, which also include end effects [209], [214], [215]. In order to compute the eddy currents, the majority of researchers solved partial differential equations with the FEM [216]. However, since the computing capacity in the 1980s was relatively low, hybrid methods, combining the good sides of integral equation and FEMs, have also been widely utilized from various perspectives [208], [217] [220]. The trend to analyze only 3-D geometries in conventional eddy-current problems characteristically differs from the approaches to model hysteresis losses in superconductors by treating them as conventional electric conductors with nonlinear power-law resistivity in two dimensions [51]. Typically, a 2-D analysis of superconductors is adequate to solve a wide class of problems, but in certain cases, a new and more accurate analysis also requires modeling in three dimensions in the case of superconductors [59], [221]. The analysis of a 2-D geometry results in simpler formalism and program implementation than the analysis in a 3-D one since, for a solution of eddy-current problem in two dimensions, it is not mandatory to solve all the components of a vector field [214]. However, edge elements for vector-valued unknowns have been used for more than twenty years [222] [225], and their usability in 3-D eddy-current problems based on the vector potential is well known and demonstrated [226], [227], as well as in 2-D problems based on the H-formulation [228]. Nowadays, there are several publicly available or commercial software packages for eddy-current simulations

15 [229] [232], and the visualization in postprocessing, skills to use the software, and the computation time are the major differences between 2-D and 3-D eddy-current simulations. However, in 3-D eddy-current computations of conventional conductor materials, such as copper, not much attention is typically paid to the conductivity, which is expected to be constant. Naturally, this is a modeling decision since, for example, copper suffers from anisotropic magnetoresistance [233], [234]. Thus, the problem is closer to the problem of utilizing power law in 3-D geometries for superconductor resistivity than what the literature related to computations suggests. Anyhow, the good correspondence between the simulations and the experiments still suggests that the approximation of constant conductivity does not lead to severe errors in loss computations as it would in the case of superconductor simulations [235], [236]. B. Computation Methods for Eddy Currents The methods to compute eddy-current losses from solved current density and/or magnetic field distributions do not differ from the methods used to solve superconductor hysteresis loss from a solution obtained with an eddy-current solver utilizing nonlinear resistivity for the superconductor. Thus, the main methods are the utilization of Joule heating and domain integrals and the integration of Poynting vector through surfaces surrounding the domain under consideration [228], [237], [238]. However, when one needs to separate the magnetization loss in the superconductor from the eddy-current loss, it is difficult to use the Poynting vector approach since these two domains are typically in direct contact with each other. The main numerical methods for eddy-current simulations are the FEM, the integral equation method (IEM), and their combinations [203]. IEMs can be combined with FEMs in two ways: The boundary terms are computed either by the volume integration method [239] or by the boundary integration method [240]. Integration methods typically suffer from dense system matrices, but fast multipole acceleration considerably increases their speed [241]. However, to the best of our knowledge, these sophisticated methods are not available in commercial software packages. All these methods can work with several formulations. Within these formulations, there is often the possibility to make different modeling decisions: For example, in the so-called T Ω formulation, the field T can be solved by Biot Savart integration [227] or by using the so-called cotree gauge [242]. Table I summarizes the main formulations used for eddycurrent simulations with the most typical methods. Several other formulations with slight variations have been reviewed in [222]. C. Computation of Eddy-Current Losses in Superconductors Although the primary concern for the computation of losses in technical superconductors is the hysteresis loss of superconducting materials, a numerical estimation of eddy-current losses is also important. Furthermore, methods to reduce eddycurrent losses, at the expense of stability, are however well known, and low ac loss conductors have thus increased matrix TABLE I MOST COMMON FORMULATIONS AND SOLUTION METHODS FOR NUMERICAL EDDY-CURRENT COMPUTATIONS resistivity [243], [244]. The contribution of eddy-current loss at power frequencies (below 200 Hz) to the total loss of a superconductor [245], [246] or a superconducting coil [247] is generally reported to be low, and many studies neglect these losses [248], [249]. This applies to both the silver-sheathed Bibased conductors [245], [247] as well as to copper-stabilized layered YBCO-coated conductors [250], [163]. For a 1-cm-wide YBCO tape having YBCO layer thickness of 0.9 μm, copper stabilization of 50 μm, and self-field critical current of 260 A at 77 K, the hysteresis loss at 60-Hz sinusoidal current with peak current of critical current is 380 mw/m. The corresponding eddy-current loss in the copper stabilizer can be estimated to be in the order of 5 mw/m [250]. However, if the tape is brought to 4.2 K, the copper resistivity still decreases, and the eddy-current loss increases. The eddy-current power per unit length P eddy in a copper stabilizer layer of a YBCO tape can be estimated as [163] P eddy = 4μ2 0 π twf 2 ρ I2 c h(i) (45) where t, w, and ρ are the copper stabilizer thickness, width, and resistivity, respectively. h is a complex function of normalized operation current i. Ati =1, it receives a value of 0.03 and, at i =0.5, avalueof10 5 (see details in [163]). Thus, from this formula, we can conclude that the eddy-current loss in the stabilizer or the matrix is proportional to the second power of the frequency and inversely proportional to the resistivity of the normal metal; in addition, it also largely depends on the current amplitude. If one wants to go a bit deeper in the analysis, one should note that a realistic operation current is considerably below the critical current; thus, the eddy-current losses can be on the level of hysteresis losses. In Bi-based Ag-sheathed tapes, the contribution of eddy-current losses is from 30% down to 1% when the operation current increases from 0.1I c to 0.8I c at a frequency of 60 Hz and a temperature of 77 K [251]. In addition, it should be noted that, at low currents, the filament arrangement can change the contribution of eddy-current losses at least from 5% to 30% of the overall losses. Thus, one should keep in mind that eddy-current losses in the matrix cannot be totally forgotten, although in many cases, they seem to give little contribution to the overall losses at power frequencies [252]. Eddy-current losses can be estimated from the overall loss measurements since their frequency dependence is different from hysteretic and resistive losses. The measured total average

16 GRILLI et al.: COMPUTATION OF LOSSES IN HTS loss power P tot of a superconductor under periodic excitation can be expressed as P tot = P }{{} res + P hyst + P eddy (46) }{{}}{{} RIrms 2 fq hyst f 2 where P res and P hyst are the average contributions of resistive and hysteretic losses in the sample, respectively; R is the sample s resistance along resistive paths where part of the current flows; I rms is the RMS value of the excitation current; and Q hyst is the superconductor hysteresis loss over one cycle. Thus, the loss in one cycle Q cycle can be expressed as [253] Q cycle (f) =Q hyst + P res f + C eddy f (47) where C eddy is a frequency-independent eddy-current loss constant (in joules per hertz). The constants C hyst, P res, and P eddy can be extracted from measurements performed with the same amplitudes but varying frequencies. In YBCO- and Bi-based tapes, materials with high conductivity, such as silver, need to be used to allow high-quality superconducting regions. Then, the eddy-current loss in these materials can pose a real challenge. In [253], the authors measured Bi-2223/Ag tapes with filaments from 1 to 55 and extracted the C eddy factor of (47). They studied frequency range of Hz. Additionally, FEM modeling was performed to estimate the losses in the silver matrix. The modeling resulted in a prediction for eddy-current loss at least four times smaller than that extracted from simulations using (47). However, the measured eddy-current loss had negligible contribution to Q cycle when the frequency was below 200 Hz, and the operation current was below 0.3 of the critical one. IV. COUPLING LOSSES There is no wide consensus about what loss components are to be put under the term coupling loss [254]. This is because coupling losses and eddy-current losses have similar nature since they both arise from currents in nonsuperconducting parts of the conductor [255]. Additionally, filament coupling influences the current penetration profiles inside a superconductor; thus, it changes the hysteresis loss [103]. In this paper, with the term coupling loss, we indicate the loss that arises in normal conducting regions when a current loop going from a filament to another via normal metal is formed due to a time-varying applied magnetic field. Here, we discuss first about performing numerical simulations from which it is possible to extract the coupling loss. Then, we review how filamentary coupling affects hysteresis losses and how this can be simulated. A. Simulating Filament Coupling Effects Coupling loss simulations that take into account the real conductor geometry need in general to be performed with models having three spatial dimensions [256]. Important characteristics for the simulations that represent real-world cases are, at a minimum, the modeling of superconducting domain, the modeling of the normal metal matrix, and the modeling of the contact resistance between these two [255], [257]. For estimating coupling loss, there is no easy way to use analytical formulas, which work in general cases. The reason for this is that the actual filament distribution largely determines the coupling loss. Additionally, separation of coupling loss from other loss components, such as eddy-current loss, is difficult from measurement results since these losses have similar nature. However, in general, one can say that a short twist pitch in conductors results in low coupling loss, and low transverse resistivity results in increased coupling loss. For a given conductor, the coupling loss increases with the second power of frequency. Analytical approaches to coupling losses are reviewed in [258]. So far, FEM-like tools have been less frequently used to model coupling losses than hysteresis losses in superconductors. There are three main reasons for this. First, when the magnetic field is not perpendicular to the current density, it is uncertain what kind of a model or a constitutive law should be adopted to model currents in the superconductor [259]. This case is usual in 3-D geometries. The CSM does not work in this case; thus, it is doubtful to use also its power-law approximation. Although power law works in some simulations of simple geometries in three dimensions [32], [260], [262] [264], experimental verification of these studies are missing. However, currently, there is ongoing interest in solving this problem [64], [183], [264]. Second, in three dimensions, it is not possible to obtain the same discretization level and, thus, accuracy, as in two dimensions. Since superconductor modeling is highly nonlinear, even solving small problems takes a considerable amount of time [32]. In three dimensions, realworld problems with low accuracy still typically require several tens of thousands of degrees of freedoms. Finally, methods to reduce coupling losses are very well known for filamentary superconductors [265]. These methods include twisting the conductor, increasing the matrix transverse resistivity, and using highly resistive barriers around filaments. Thus, particularly in LTSs, the coupling losses are small, and from the perspective of wire development and applications, the establishment of methods for coupling loss simulations has not been crucial. However, utilization of these manufacturing methods with a filamentary HTS, such as Bi-based or MgB 2 wires, is still an unknown issue. Additionally, these methods are not applicable for striated YBCO conductors as such, and there are many recent developments aiming at reducing YBCO ac losses by filamentarizing the end product [257], [266] [269]. This effort must also include considerations on coupling losses; thus, their modeling is important since, even in the 10 Hz regime, the coupling loss in striated YBCO can be important [270]. An alternative way of modeling coupling losses in a program that solves directly Maxwell s field equations is the utilization of equivalent circuit models [269]. This kind of models have been mainly developed for LTS cables, such as cable-in-conduit conductors (CICC) [271]. These models are fairly complicated (as shown in Fig. 7) and require measurement of contact resistance between many points. However, they are still applicable to coupling loss simulations of Roebel cables [272], which also have a loose cable structure such as many CICCs [273].

17 Fig. 7. Circuit analysis model for coupling loss computation of CICC. Reprinted with permission from Elsevier [271]. Fig. 8. Schematic view of the thin strip approximation of a striated YBCO tape. One magnetization current loop is shown for the case of uncoupled filaments. In addition, when FEM-type simulations target high accuracy, the laborious characterization of interface contact resistance becomes mandatory anyway. When coupling losses are modeled in thin conductors by directly solving the fields, it is possible to take 2-D approaches. It can be estimated that there is no current flow from the bottom of the tape to the top or vice versa [274]. This modeling decision is generally called thin strip approximation [257] (see Fig. 8). Striated YBCO conductors are examples of geometries that can be simulated with such approximation. Thus, although the tape can have complicated 3-D geometry, at every cross section, in local coordinates, one current density component always vanishes [59]. Since superconductor hysteresis losses are computed simultaneously to coupling losses, we will review this method more deeply in the following. B. Influence of Filamentary Coupling on Hysteresis Losses In 2-D hysteresis loss simulations of a multifilamentary conductor or a cable consisting of several conductors, it is necessary to make a modeling decision concerning filament coupling effects when a magnetic field is applied. However, typical programs consider only two extreme cases: Either the filaments are perfectly coupled or not coupled at all [103]. Full coupling means that the filaments interact together similarly to a single larger filament. Uncoupled situation means that filaments operate completely individually, although the fields that they generate can alter current distributions in the other filaments [274], [275]. Three-dimensional simulations are typically mandatory to study the intermediate cases, e.g., the effect of matrix resistivity and filament-to-matrix contact resistance on the current distribution in a superconductor [55]. Naturally, these 3-D simulations also consider actual sample sizes and do not require infinite long samples. However, for high-aspect-ratio tapes, such as YBCO, it is possible to assume that the current density is always parallel to the tape s wide surface. Then, the problem reduces to a 2-D one, which allows different solutions with intermediate coupling where the resistivity between the filaments is varied [274]. Naturally, in this model, coupling currents can only flow from the edges of the superconductor to the normal metal, which limits the use of this model to the study of in-plane coupling currents (see Fig. 8). The field problem in this thin strip approximation is formulated in terms of current density vector potential such that J = (nt ), where n is the unit normal to the computation surface, and T is a scalar field. The governing equation includes also integral terms since the field generated by the surface currents is computed by integral equations. Power-law resistivity has been used to model the electric behavior of the superconductor [274]. In [274], the influence of the resistance between striations of a multifilamentary YBCO tape on hysteresis loss and coupling loss was studied with the thin strip approximation. The striated YBCO tape had 20 filaments. When the frequency was below 5000 Hz, coupling loss increased with decreasing resistance between the filaments. Roughly, one order of magnitude increase in contact resistance between filaments decreased the coupling loss by one order of magnitude. In addition, the variation of hysteresis loss as a function of filament-tofilament contact resistance was studied. The difference of losses in the fully coupled case and the uncoupled case was a bit more than one order of magnitude. Losses computed with different values of contact resistance fell between these two limits. At 50 Hz, the coupling loss was still, in every case, much lower than the hysteresis loss, and the most important factor of the filamentary coupling was the change in hysteresis loss and not the amount of coupling loss. Naturally, these results depend largely on the number of filaments as well as on the presence of a stabilization layer. Additionally, in practical applications, the end of the YBCO tape needs to be soldered to current lead. Then, the filaments are well coupled, and introduction of filament transposition becomes mandatory. In Figs. 9 and 10, we present typical current distributions for the two extreme case of fully coupled and uncoupled filaments, when an external magnetic field is applied. In these 2-D simulations, the current density is a scalar quantity, and it flows perpendicular to the modeling plane. We had two filaments, both having a radius of 0.5 mm and separated by 1 mm. The matrix eddy-current losses were neglected in this study. We applied vertical magnetic flux density 50 sin(2πft) mt. We had field independent J c of 127 A/mm 2. We utilized A V formulation for the simulation [36]. As shown, in the case of uncoupled filaments, the current density distributions look similar to one filament alone at the applied field. For the coupled case, the interaction of filaments is clearly visible. This leads to higher magnetization of the superconductor and, thus, to higher losses. Similar simulations for the cross section of YBCO Roebel cable were presented in [103].

18 GRILLI et al.: COMPUTATION OF LOSSES IN HTS Fig. 9. Current penetration simulation of the coupled case. Red and blue represent currents toward and from the plane, respectively. Black arrows represent magnetic flux density vectors. Gray background presents the mesh for numerical computation. View is displayed from the time when applied field reaches its maximum. Fig. 10. Current penetration simulation of the uncoupled case. Details correspond to those in Fig. 9. V. I NFLUENCE OF MAGNETIC PARTS The presence of magnetic materials in superconducting systems is not uncommon because magnetic materials can be a part of superconducting wires or a component of superconducting devices. YBCO-coated conductors manufactured by the Rolling-assisted Biaxially Textured Substrates (RABiTS) technique typically have a ferromagnetic substrate. First-generation HTS tapes (Bi-2223) or MgB 2 wires can be manufactured with magnetic stabilizer/reinforcing layers. Superconductors and magnetic cores are the primary components of many power applications, such as generators, transformers, and motors. The magnetic components can influence the ac losses in superconducting systems in two ways: 1) They change the magnetic profile inside the superconductor (thereby affecting the superconductor losses); and 2) they contribute additional hysteretic losses to the system. Therefore, the prediction of the effect of magnetic components on the total ac losses in a superconducting system is essential and has attracted significant effort from the research community. In general, two types of approaches can be distinguished: analytical and numerical models. Analytical models provide relatively simple formulas to compute the losses in the superconductor; their main limitations concern the type of analyzed Fig. 11. Schematic view of the cross section of a YBCO-coated conductor with a magnetic substrate. Not drawn to scale. geometries of the superconductor and/or the magnetic parts, the physical properties of the material (typically, CSM for the superconductor and infinite magnetic permeability for the magnetic material), and the fact that they neglect the losses in the magnetic material. On the other hand, numerical methods can overcome these limitations, although they generally require a larger computational effort. In the following, we shortly review the main methods found in the literature. All of these models are usually solved for the 2-D cross section of a conductor or a device, which is perpendicular to the direction of the transport current. For example, these models can be used to calculate the magnetic and current density distribution in the xy plane for the cross section of a superconductor with ferromagnetic substrate, as shown in Fig. 11. A. Analytical Models Genenko et al. pioneered analytical calculations of ferromagnetic/superconductor heterostructures by investigating the effect of an infinite magnetic barrier placed next to a current-carrying superconducting strip [276]. It was found that the effects on the current density distribution inside the superconductor mostly depend on the barrier s shape and distance from the superconductor and only marginally on its magnetic permeability. The method was later used to propose configurations aiming at increasing the transport capability [277] and reducing the ac losses [278] of superconducting strips. The aforementioned models, however, consider magnetic domains that are not embedded in the structure of the tape or wire, and as such, they are not suitable to describe technical conductors used in applications, e.g., coated conductors with magnetic substrate or MgB 2 wires with magnetic material. An analytical model to calculate the magnetic distribution and ac loss of an isolated YBCO tape with magnetic substrate of infinite permeability μ was proposed in [279]. This analytical model derives formulas to calculate the magnetic field, current distribution, and ac loss in a superconducting strip with a magnetic substrate (SC/FM strip conductors) for some special simple cases. In the model, a 2-D magnetic field, H =(H x,h y ) was considered as a complex field H(ζ) =H ay + ih ax, where ζ is a complex variable ζ = x + iy. The primary idea of the model is to solve the generalized Biot Savart law for a SC/FM strip [279]: H(ζ) =H ay + ih ax )+ 1 2π a a dx K z(x )+iσ m (x ) ζ x (48)

19 Fig. 12. Normalized imaginary part of ac susceptibility as a function of normalized applied magnetic field. where H a = H axˆx + H ay ŷ is a uniform applied magnetic field, K z (x) is the sheet current in the superconducting strip, and σ(x) is the effective sheet magnetic charge in the ferromagnetic substrate. Both K z (x) and σ(x) can be determined from Faraday s laws [6, Eqs. (4) and (5)] and must satisfy the following constraints: The total net current in the superconducting strip must be equal to the transport current, and the total net magnetic charge in the substrate is zero. In order to be able to derive the explicit results for the current and magnetic field distribution, the model must assume that the ferromagnetic substrate is an ideal soft magnet with infinite permeability, i.e., μ. With this approximation, the outside of the magnetic substrate has only the perpendicular component on its surface. This determines the boundary condition at the surface of a ferromagnetic substrate. First, the model was solved for a SC/FM strip conductor in the ideal Meissner state and exposed to either a perpendicular or parallel magnetic field, or a transport current. In the ideal Meissner state, the perpendicular magnetic field component on the surface of the superconducting strip must be zero. Based on the CSM, formulas for the case of the SC/FM strip conductor exposed in a perpendicular dc or ac field were also derived. In ac perpendicular applied field H a (t) =H 0 cos(ωt), the time dependence of magnetic moment per unit length m y (t) is derived explicitly and can be expressed as the Fourier series to determine the complex ac susceptibility χ n + iχ n at the nth harmonics. The superconducting ac loss Q generated per unit length of a SC/FM strip in an ac perpendicular field can be therefore calculated from the imaginary part of the first harmonic χ 1 as follows [88]: Q(H 0 )=πμ 0 H 2 0 χ (H 0 ) (49) Fig. 12 plots the normalized imaginary part χ 1/πa 2 (a is the half width of the conductor) for the fundamental frequency (n =1) of the ac susceptibility as a function of normalized applied field H 0 /j c d s (j c and d s are the critical current density and thickness of the superconducting strip, respectively) for strip superconductors with ferromagnetic substrate (SC/FM) and without magnetic substrate (SC/NM). Experimental results reported by Suenaga et al. in [280] indicated that, in perpendicular applied field, the ferromagnetic hysteresis loss in the Ni W substrate of a RABiTS tape would be significantly lower than that in the superconducting layer and can be therefore negligible in the total magnetization losses. The model given predicted that magnetization losses in the SC/FM conductor are higher than those in the SC/NM counterpart at low applied fields, i.e., μ 0 H 0 < 0.14j c d s.for a conductor with j c = A/m 2 and d s =2.5 μm, the crossing point μ 0 H 0 =0.14j c d s =4.9 mt is consistent with the experimental results reported in [280] for a YBCO tape with a Ni W substrate. That quantitative agreement suggests that the model with infinite permeability can be used to quantitatively predict ac losses for a superconducting strip with a magnetic substrate of large but finite permeability exposed to a perpendicular applied magnetic field. However, this model currently does not apply for other important cases, which represent better for the real operating conditions of HTSs in devices, such as a SC/FM strip carrying an ac transport current or/and exposed to a perpendicular/parallel applied magnetic field. In addition, as already pointed out, this model does not provide a way to predict the ferromagnetic losses in the magnetic substrate, which may be significant in some situations [281], [282]. Thus, the accurate prediction of the total ac losses in an SC/FM system is difficult to be obtained with this analytical model, and numerical methods are therefore needed. B. Numerical Models Two numerical models have been widely used to compute losses in superconducting devices with magnetic parts, and they both employ the FEM. While those FEM models are supposed to be suitable with soft magnetic materials whose μ(h) is considered to be approximately the same for both the upper and the lower branches of a B H loop, they may be also applied for a system with hard magnetic materials by using μ(h) approximately determined from the virgin magnetization curve. An accurate model for hard magnetic materials that is capable of considering the hysteresis in B H loop would be very complicated and has not been developed yet. 1) H-Formulation FEM Model: The H-formulation FEM model to compute ac losses for superconducting systems with magnetic materials was first proposed in [281] and [282]. This model is capable of taking the field dependence of permeability and ferromagnetic loss of the magnetic materials into account. With respect to the standard version of the H-formulation [29], [228], the writing of Faraday s equation need to be modified to take into account that μ r depends on H and, hence, on time. As a consequence, μ r cannot be taken out of the time derivative in the governing equation (μ r μ 0 H) + ρ J =0 (50) t and the equation needs to be expanded. The details of the new formulations can be found in [281] and [282]. The hysteresis loss in the magnetic material Q fe is the area of the B H loop and is usually a function of B m, which is the maximum magnetic field of that loop. The hysteretic loss

20 GRILLI et al.: COMPUTATION OF LOSSES IN HTS Fig. 13. FEM and experimental results for the ac loss components in a 1-cmwide RABiTS YBCO tape (I c = 330 A) carrying an ac transport current. at each location inside the magnetic material can be therefore calculated from the maximum value of the induction field B m at that location during an ac cycle. Therefore, based on the Q fe (B m ) relation fitted from the experimental data, the total magnetic loss q fe can be estimated by taking the integral of Q fe (x, y) over the cross section A SUB of the magnetic material [281], [282], i.e., q fe = Q fe (B m (x, y))dxdy. (51) A SUB This model is quite versatile and can be used for HTS devices with complicated geometries, as long as the field dependence of μ r and Q fe for the magnetic materials in the device are known. In fact, the model has been used extensively to calculate ac losses for superconducting systems using RABiTS-coated conductors with an Ni W substrate. For this type of conductors, the field dependence of μ r and Q fe of the substrate material was determined experimentally and published in [283]. The functions fitting the experimentally observed field dependence of μ r and Q fe were entered in the model and used for the calculation of ac losses [281], [282]. As an example, Fig. 13 shows the comparison between calculated and measured total transport losses in a RABiTS tape as a function of transport current. In the simulation, the superconducting loss in the YBCO layer and ferromagnetic loss in the substrate were calculated separately, and the sum of these loss components provides the total selffield loss in this sample. In general, the modeled total self-field loss agrees well with that obtained from the measurement. The good agreement between simulations and experiments was also observed for other cases of RABiTS conductors with Ni W substrates, such as the total ac loss in parallel applied field [281], the total ac losses for stacks of two parallel tapes [282], bifilar coils [284], racetrack coils [167], pancake coils [285], [286], or cables [287]. 2) A-Formulation FEM Method: This method is based on the following general equation: ( 1 μ A ) = j = σ ( ) A t + φ (52) where A is the magnetic vector potential, and φ is the scalar potential. For infinitely long conductors and circular rings, φ is constant in the cross section. Miyagi proposed a 2-D edgebased hexahedral finite-element modeling technique to solve that equation with constant and variable μ. At the beginning, (52) was solved for conductors with a substrate of constant μ, but a considerable discrepancy between calculated results and experimental data was observed [288]. The model was therefore improved to take the field dependence of permeability into account for better accuracy [289]. The authors used the Newton Raphson method to handle the nonlinearity of permeability of the substrate. The detailed algorithm and techniques used to solve (52) with nonlinear σ(j) and μ(b) can be found in [289]. This model, however, has not been validated against experiment. Although authors did not use this method to predict the ferromagnetic loss in the substrate material, this can be done by a similar way as presented in the H-formulation technique. Another numerical technique to solve (52) was introduced in [45]. This method is based on the CSM and two assumptions: 1) There exists a neutral zone inside the conductor where the current density j (or the electric field E) vanishes during the entire cycle; and 2) the magnetic flux monotonically penetrates from the surface when the applied current or magnetic field increases monotonically. Because the electric field is zero, the vector potential A c in the neutral zone must be uniform and satisfies φ = A c / t. Then, the authors introduced a new variable, the so-called calibrated vector potential A = A A c, which is zero in the neutral zone. The electric field can then be calculated by the following simple equation: E = A / t. (53) Since A in the neutral zone is zero, A (x, y, t) is the magnetic flux per unit length between the neutral zone and the point (x, y). Thus, based on the CSM, the current density from one time step can be analytically calculated from the calibrated vector potential of the previous time step. The ferromagnetic hysteretic loss in the magnetic material was estimated by a similar technique, as described in the H-formulation method. The details of the calculation procedures can be found in [45]. The authors employed the model to calculate ac losses in a Bi-2223 conductor with an elliptical cross section for two cases: 1) a conductor carrying an ac transport current; and 2) a conductor exposed to a perpendicular applied ac field. The edges of the conductor are covered with ferromagnetic material (Ni) for ac loss reduction. As an example, Fig. 14 depicts the calculated evolution of magnetic and current density profiles in that conductor when the applied field changes from 5 to 50 mt. The theoretical results of ac losses in this conductor were then compared with the experimental data [45]. The quite good agreement between the theoretical prediction and the experimental data has validated this numerical method, and it can provide another option to calculate ac loss and electromagnetic behavior in a SC/FM system [290]. 3) Other FEM-Based Models: Other numerical models, based on the FEM, have recently been developed to study the interaction of superconducting and magnetic materials.

21 TABLE II IDEAL AND PRACTICAL CARNOT SPECIFIC POWER FOR OPERATING TEMPERATURE FROM 4.2 TO 273 K TABLE III SUMMARY OF THE ESTIMATED HEAT LOAD OF HTS DEVICES AT THE OPERATING TEMPERATURE Fig. 14. Evolution of magnetic and current density profiles calculated for the initial increase in the perpendicular magnetic field. From top to bottom, the applied fields are 5, 12, 22 and 50 mt, respectively. Density of electrical current is characterized by the darkness of the color, the magnetic field is represented by the lines of constant vector potential. Figure taken from [45] and reprinted with permission from IOP. Genenko et al. [291] studied the influence of the shape of magnetic shields to reduce the ac losses in a superconducting strip. It was found that not only the shape of the shield but also the magnetic permeability of the magnetic material plays a key role for that purpose if the size of the shield is comparable to or less than that of the superconductor. The same group also developed a model to simulate superconducting strips with magnetic substrate in applied ac field based on an electrostatic magnetostatic equivalence [292]. Farinon et al. extended their adaptive resistivity algorithm [293] to take into account magnetic materials [294]. VI. AC LOSSES IN HTS POWER APPLICATIONS Calculation of ac losses is essential for HTS power applications since the design of the cooling systems needs to take into account the thermal loads. The cooling system maintains the required cryogenic temperatures during the operation of HTS devices. AC loss is a contribution source of the heat load; thus, the importance of ac loss lies in two aspects: The calculation is essential for the engineering design of cooling systems, and the calculation acts as a guide to explore methods for minimizing ac losses. This section discusses the impact of the ac losses in the most important power applications and the main contribution of the ac loss for each case. A. Cryogenic Considerations One important reason why HTSs are very attractive is that they can operate without liquid helium, which is a very ex- pensive type of cryogen. Both BSCCO- and YBCO-based superconductors can operate at 77 K; the price per liter of liquid nitrogen is times lower than that of liquid helium depending on the purchase location [295]. The latent heat of liquid nitrogen is 161 kj/l at one atmosphere [296]. In order to achieve a larger critical current and higher magnetic field, HTS devices are operated at temperatures lower than 77 K with the aid of cryocoolers. One important parameter in determining the efficiency of cryocoolers is Carnot efficiency. The inverse of the Carnot efficiency is called the Carnot specific power and is the number of watts required at ambient temperature to provide 1 watt of refrigeration at the lower operating temperature. Column 2 in Table II presents the ideal Carnot specific power (watt input per watt lifted) [297]. However, commercial refrigerators run only at a fraction of Carnot efficiency; thus, column 3 gives a more realistic value for the specific power (watt input at 303 K per watt lifted at T op ). Thus, in the cryogenic system design, the cooling power required for removing the heat load due to the ac losses can be calculated using Table II. Table III summarizes the estimated heat load of HTS devices at the operating temperature [297]. B. Techniques for AC Loss Reduction Ac loss calculation can be used to find solutions to reduce ac losses of HTS power devices. As it has been described in this paper, ac losses have different components: hysteretic, ferromagnetic, eddy-current, and coupling current losses. Hysteresis losses can be reduced by decreasing the filament size [243]. Eddy-current losses can be

22 GRILLI et al.: COMPUTATION OF LOSSES IN HTS reduced by increasing the stabilizer resistivity. Coupling current losses can be reduced by increasing the matrix resistivity. Ferromagnetic losses can be eliminated by avoiding ferromagnetic materials. There are a few general methods for reducing ac losses of HTS applications. Controlling the amplitude and magnetic field orientation are effective ways for reducing ac losses [298], [299]. Using flux diverters or shielding layers are common means to reduce the magnetic field, which in turn reduces ac losses. Since HTS-coated conductors produce significantly higher ac losses in a perpendicular magnetic field, reducing the perpendicular component is an effective way to lower the ac losses. The large cross-sectional aspect ratio of HTS-coated conductors leads to a large hysteresis loss in external fields. Filamentization of the coated conductor proved to be effective in reducing hysteresis losses [300], [301]. Sumption et al. pointed out that the losses for a coated conductor subdivided into ten filaments is shown to be reduced significantly [301]. However, filamentarization introduces coupling losses, caused by the relatively good conductivity of nonsuperconducting barriers and remnant superconducting bridges between neighboring strips [302]. Another contribution of losses in a striated conductor comes from magnetic flux coupling between filaments because magnetic field distribution inside the superconducting filaments is influenced by the field concentration in the nonsuperconducting region [302]. C. Electrical Power Applications In this section, ac losses in various applications of HTS power applications are briefly reviewed. Typical techniques for reducing losses in different applications are described. 1) Rotating Machines: In an ac synchronous machine, superconducting coils can either be used as armature coils or rotor coils. If used as stator coils, superconductors will carry an ac and will be exposed to a rotating field, thus producing large ac losses [303], [304]. If used as the rotor coils, the magnetic field of the superconducting coils is phase-locked with the stator field. Under balanced load conditions, the rotor experiences a dc field from the armature. However, nonsynchronous disturbances, which are usually a result of unbalanced loads, transient events caused by system faults, and noise from the grid, can produce nonsynchronous ac effects that affect the superconducting coils on the rotor. The induced time-varying fields in the armature at frequencies other than the synchronous frequency induce compensating ac in the field windings, which in turn causes ac losses in the superconducting coils on the rotor. In a synchronous superconducting wind generator design, the superconducting coil windings on the rotor will also experience external ripple field from the power electronics excitation ripple and wind turbulence. This effect is investigated in [305]. In the Siemens 4-MW demonstration motor experiment [306], superconducting coils are used as rotor coils, whereas conventional copper coils are used in the armature. The experiment shows that the associated rotor losses, including ac losses, are smaller than 120 W, which would cost 10-kW compressor power; thus, the cooling system would decrease the system efficiency by 10 kw/4 MW =0.25%, which is quite small for a large-scale application [306]. This power requirement is encouraging because it shows that the cooling system to remove the heat load is not a significant barrier for a large superconducting machine. 2) Power-Transmission Cables: In HTS powertransmission cables, superconducting tapes are usually helically wound around a cylindrical former. A very simple but useful model for calculating ac loss is the monoblock model in which the cable is considered as a tube of superconductor [141]. Although the monoblock model does not take into account the interaction between phases and the geometry of the tapes, it constitutes a frequently used engineering estimation for large-scale applications [307]. Another useful model is the electrical circuit model [308], [309], which models cables with resistance and inductance connected in series and in parallel. The resistance values are calculated from Norris equation, whereas the inductance values are calculated from self-inductance and mutual inductance formulas. This model is able to take into account the interaction between the conductor layers, the insulation layer, and the shielding layer. Clem et al. proposed that the ac loss of HTS powertransmission cables could be broken down into six different mechanisms [310]: Bean model losses as in a Bean slab, losses in helically wound phase conductors arising from finite gaps between adjacent tapes, losses arising from the polygonal configuration of the conductor cross section if the tapes do not conform to the cylindrical former, flux cutting loss, losses when the current exceeds the critical current, and losses in the possible magnetic substrates. The key parameters affecting ac losses in power-transmission cables are the number of conductors, the conductor s width, the gap size between adjacent tapes, and the radius of the former and the current distribution. One common design of superconducting cables is the monolayer design. AC losses in a monolayer cable do not significantly increase over one single tape since most of the flux lines lie outside the superconducting layers in the cable. It has been found that using narrower tapes and/or decreasing gaps between adjacent tapes are effective solutions for ac loss reduction [257], [311]. Due to the mechanical requirements and architecture of the conductor, a monolayer cable consisting of mm-wide tape would have minimized ac losses [312], [313]. In a two-layer counter-wound cable using coated conductors, gap, and polygonal losses, flux transfer losses in imbalanced two-layer cables and ferromagnetic losses for conductors with magnetic substrates are the major contribution to ac losses [310]. Current imbalance and ferromagnetic losses can be minimized by orienting the substrates of the inner winding inward and the outer winding outward [310]. In multilayer cables, the ac losses behave differently. It has been found that ac losses in cables with many layers are hardly affected by gaps between adjacent coated conductors, as well as a lateral critical current distribution [313]. 3) Compact High-Current Cables: Different from powertransmission cables, superconductors are often required to be

23 parallel or stack placed to create a high-current-density cable for applications in high-field magnets and busbars. Several kiloampere to tens of kiloampere cables are expected in these applications. Twisting and transposition are typical methods for reducing ac losses in parallel and stack-placed superconducting cables [103], [243], [301]. This is because magnetic flux coupling can be reduced by transposing the conductors. However, it is not easy for HTSs be transposed as LTS wires since, often, HTSs are made into flat tapes. Recently, significant effort has been devoted to develop HTS Roebel cables, Rutherford cables, and other types of a twisted stack of coated conductors for high-current applications [314] [316]. Filamentarization is another technique to reduce ac losses in highcurrent cables [252]. In addition, decreasing the aspect ratio (width/thickness) of the cable by assembling stacks of strands instead of individual strands in cables can also reduce ac losses [156]. 4) Fault Current Limiters: For resistive-type SFCLs, it is preferable to contain a long length of superconductors in a small volume. One type of SFCLs uses an array or a matrix of superconducting tapes [317] [319]. Fo. For ac loss calculation, a stack of superconductors can be modeled for this type of SFCLs [319]. Another design is to use a noninductive superconducting coil [174], [284]. To calculate ac losses of this type of SFCL, HTS coils need to be studied [174], [284]. In the dc-biased iron-core-type SFCL, superconducting coils carry dc; thus, there are no ac losses [320]. To calculate losses in the shielded iron-core-type SFCL, one needs to analyze the ac losses of a superconducting bulk cylinder or a tape. An effective way to reduce ac losses in a resistive SFCL is the bifilar design [174]. This design can significantly reduce the inductance and also ac losses by decreasing the magnetic field volume [321]. Another technique for reducing ac losses is to use the two-in-hand structured wire [284]. The loss of a coil using two-in-hand structure wire is only half of that generated in an isolated tape per unit length. The ferromagnetic losses can be minimized by orienting the substrates of the inner winding layer inward and the outer winding outward (face-toface configuration [310]). In a commercial 12-kV resistive SFCL, ac losses account for 60% of the total thermal load [321]. 5) Transformers: From the point of view of engineering design, ac losses in transformers can be regarded as ac losses in HTS coils. AC losses come from the transport losses in the primary and secondary windings and the magnetization losses due to the magnetic fields, in particular, the leakage magnetic field in the end region of superconducting windings. One technique for reducing losses is to use a compensation coil to cancel the unwanted magnetic field, particularly the perpendicular component. Heydari et al. discussed how to use auxiliary windings to reduce the leakage magnetic field, hence, the ac losses [322]. Another way of reducing leakage magnetic field is by using flux diverters to shape the field. Al-Mosawi et al. examined the validity of this method and proved that ac losses are reduced by 40% in a 10-kVA superconducting transformer using Bi-2223 conductors [323]. Low-ac-loss conductor designs, e.g., Roebel cables, are also proposed to be used in transformers for carrying a large transport current [315]. In an example study, it is estimated that the losses only reduces the efficiency by 0.2% in a 800-kVA transformer using BSCCO tapes, with a penalty factor of 20 taken into account [324], [325]. This shows that, such as in a large superconducting machine, the cooling power is only a fraction of the total output power of superconducting transformers. 6) Superconducting Energy Storage: The energy storage devices in superconducting energy storage systems (SMES) are superconducting magnets. The ac loss associated with changing currents and magnetic fields are nearly zero if the magnets operate in storage mode. The two significant continuous energy losses for HTS magnets are the joint resistance of HTSs [326] and the normal resistance in the power electronics transistors. However, ac losses may occur in the system as a result of variation in the operating current and magnetic field due to the charging or discharging of energy [327]. The ac loss per unit time might be significantly higher than in normal conditions during a discharge process since a SMES is expected to react quickly in response to a grid incident, i.e., a quick discharge to maintain grid stability [328] [330]. In certain discharge circumstances, the ac loss per unit time will be much larger than a cryocooler capacity. An example study shows that, for 1 s of the discharging time in a 600-kJ SMES system, the total ac loss is W, which is a large amount of heat, increasing the system temperature by 3 K [331]. Thus, the temperature rise during a discharge process should be analyzed in the design work. The ac loss calculation is also essential for a SMES cyclic efficiency estimation. VII. CONCLUSION AND OUTLOOK We have presented a review of the literature on ac loss computation in HTS tapes, wires, and devices. Two-dimensional calculations by different methods are now fully developed to the point that ac losses in superconductors with geometries of arbitrary complexity can be quickly computed with good accuracy. Magnetic materials can be also taken into account, and at least with certain models, their nonlinear magnetic characteristics can be modeled. In addition, situations different from pure ac excitations (such as ramps, pulses, and combination of ac and dc field and currents) can be handled. We are confident to say that all the situations relevant for realistic applications can be modeled in 2-D. The most challenging step for the scientific community will be the extension to 3-D calculations. The challenge resides not only in the numerical difficulties related to 3-D modeling (choice of the mesh, increase in memory requirements, and computation time, etc.) but also in the choice of an appropriate constitutive relation for modeling the different types of superconductors in 3-D and in its implementation in the numerical codes. AC losses in different electrical power applications are investigated and summarized. It can be concluded that ac loss calculation in such applications is important for cooling system design and examining innovative methods for heat load reduction.

24 GRILLI et al.: COMPUTATION OF LOSSES IN HTS Fig. 15. Vortex moving at velocity v creates an electric field E v due to its magnetic field B v. The cylinder represents the magnetic core of the vortex. The wire figure is an arbitrary integration volume (see Appendix A). APPENDIX A ELECTRIC FIELD CREATED BY MOVING VORTICES This appendix outlines the main steps to obtain the macroscopic electric field created by moving vortices, i.e., (11). The electric field created by a moving vortex, i.e., E v,at a certain speed v is as follows. This electric field is E v = A v / t φ v, where A v is the vector potential in the Coulomb gauge, and φ v is the scalar potential. The time dependence is due to the movement of the vortex. Then, A v (t, r) = A v [r r v (t)], where r v is the position of the center of the vortex. Thus, E v =(v )A v φ v. By means of differential vector relations, this equation becomes E v = B v v + (v A v φ v ). For a long straight vortex, the average in a volume V of height l and base surface S (see Fig. 15) containing a portion of the vortex is E av,v =[1/(lS)] V E v d 3 r, where the volume integral is E v d 3 r = B v vd 3 r + (v A φ v )ds. (54) V V In the given equation, V is the surface of the volume V.The second term in the equation vanishes. This is because of two reasons. First, both A v and φ v vanish far away from the vortex center. Second, the integral contributions for the portion close to the vortex center cancel by symmetry. This reasoning also applies if the volume V is the cell of the vortex lattice. Both the vector and scalar potentials created by the vortices in the lattice vanish on the border of the cell [68], [332]. The integral of v A v and φ v on the top and bottom surfaces also vanishes by symmetry. Then, the electric field created by one vortex is E av,v =(Φ 0 v)/s; hence, for n vortices per unit surface, the average electric field is V E =Φ 0 n v = B v. (55) Equation (55) for the average electric field assumes that the dimensions of the base of the volume (see Fig. 15) are larger than the vortex separation. In practice, the vortex separation is often much smaller than the sample dimensions. Thus, E is approximately the local electric field in a differential volume of the sample. For thin films in a perpendicular field, the vortices become Pearl vortices [333], but the above analysis is still valid. APPENDIX B PROPERTIES OF THE VECTOR AND SCALAR POTENTIALS AND COULOMB GAUGE Several calculation methods assume the Coulomb gauge for the vector and scalar potentials A and φ. In this paper, the discussions in Sections II-C2 and II-C4 assume the Coulomb gauge. For the reader s convenience, this appendix lists the main properties of the Coulomb gauge and summarizes their deduction. Since B =0, the flux density can be written as a rotational of function A as B = A. This defines the vector potential A. From the differential Ampère s law with negligible displacement current, i.e., B = μ 0 J, the vector potential is A(t, r) = μ 0 4π V J(t, r ) r r d3 r + Ψ(t, r) (56) where the function Ψ(t, r) is an arbitrary scalar function. The equation above can be deduced by differential vector analysis, as in [69, pp ]. The choice of this arbitrary function defines the gauge. The Coulomb gauge is defined as Ψ=0. Then, the vector potential in Coulomb gauge, i.e., A c,is A c (t, r) = μ 0 4π V J(t, r ) r r d3 r. (57) As a consequence, A c =0because J =0. In addition, A c approaches zero as r approaches infinity, except if there is current at infinity. Another property of A in the Coulomb gauge is that it only depends on the current density. The scalar potential is defined as follows. Faraday s law and B = A imply that (E + A/ t) =0. Then, E + A/ t can be written as the gradient of a scalar function φ; hence, E = A/ t φ. Inserting this equation into the Gauss law, i.e., E = ρ/ɛ 0 and using equation (56) results in 2 (φ + Ψ/ t) = ρ/ɛ 0, where ρ is the charge density and Ψ is the function that defines the gauge of the vector potential. The unique solution of this Laplace equation, under the assumption that φ + Ψ/ t vanishes at the infinity, is φ(t, r) = 1 4πɛ 0 V ρ(t, r ) r r d3 r Ψ t. (58) The Coulomb gauge for A is defined as Ψ=0. With this gauge, the scalar potential becomes the electrostatic potential as follows: φ c (t, r) = 1 4πɛ 0 V ρ(t, r ) r r d3 r. (59) For infinitely long conducting wires, the vector potential in the Coulomb gauge is the integral along the length of the wire of (57), resulting in A c (t, r) = μ 0 J(t, r )ln r r d 3 r. (60) 2π V