Remanent Magnetic Flux Distribution in Superconducting-Ferromagnetic Layered Heterostructures

Stacks of superconducting tape, made of second-generation coated conductors, present an interesting alternative to bulk superconductors. The layered geometry allows incorporation of other functional, namely ferromagnetic, materials. The ferromagnetic material changes the field configuration within the superconductor and can be used to shape or divert the flux density from its interior. This can be used to optimise the remanent flux density after magnetisation. FEM simulations are employed to test the performance of a variety of superconductor/ferromagnet structures and explain the role of the ferromagnetic material in the composite structure. Cases where trapped field is much higher and lower than remanent magnetisation of the ferromagnetic material are considered, showing that performance can be enhanced or maintained by using ferromagnetic material while reducing the usage of costly superconducting tape.


Introduction
A number of works have demonstrated that the use of ferromagnetic materials can improve the properties of trapped field magnets (TFMs). Most of the work so far has been focused on trapped field magnets made of bulk superconductors (SCs), where ferromagnetic materials (FMs) have shown the ability to tailor the magnetic field profile [1], reduce crossed-field demagnetisation [2] or just increase the trapped field above a bulk superconductor [3]. Other work on second-generation high-temperature superconducting (2G HTS) tapes using numerical modelling showed that magnetic shielding of filaments by ferromagnetic material reduces AC losses in selffield conditions due to decoupling of the filaments and, at the same time, it increases the critical current of the composite [4]. Previous work in the group has shown that a composite TFM made of layers of HTS tape layered with ferromagnetic layers increased the trapped field above the stack as compared to a stack made of only HTS tape [5]. However, the magnetisation was performed using pulsed field magnetisation (PFM) where thermal effects have a significant role; hence, it is not clear if it was the thermal or magnetic properties of the ferromagnetic layers that have led to the improvement in trapped flux. This work aims to use finite element method (FEM) modelling to gauge the possible improvements to trapped field using ferromagnetic materials by simulating "fully saturated" TFM, similar to a condition after field cooling magnetisation (FCM). Several cases are modelled considering bulk FM material outside the stack or FM layered together with SC material with or without considering the anisotropy of the critical current J c (B, θ).

Modelling Framework and Material Properties
The FEM modelling was performed in Comsol Multiphysics 5.3a. An axis-symmetric model using A formulation (magnetic vector potential) was set up to simulate a cylindrical stack.
A stationary (time-independent) study was used and the current density in the superconductor was set according to the E-J power law: where J c (B, θ, T) is the critical current density dependence on temperature, T, the magnetic field magnitude, B, and its direction, θ, determined for a critical current criterion, E 0 . While, E is the average electric field in the stack after some relaxation period, in this case taken as 1 nV/m corresponding to relaxation of~5 min after magnetisation procedure, as determined by time-dependant simulations [6]. For this work, the n value was fixed to 20. The critical current data J c (B, θ, T) was taken from [7], corresponding to SuperOx superconducting tape and the data is shown in Fig. 1a, b. A homogenised model was used, i.e. individual HTS layers were not modelled, instead the superconducting domains carried an engineering critical current density, J e = I c /A, where A is the total cross-section area of the tape that corresponds to the SuperOx tape measured in [7].
Two ferromagnetic materials were considered in this work, namely Ni-W alloy, which is used for substrate in HTS tape made by RABiTS process (e.g. in American Superconductor tape) and a CoFe alloy, trademark Vacoflux 50, which has a very high saturation magnetisation close to 2.3 T. Typical applications for Vacoflux 50 include rotating machinery. Magnetic properties for Ni-W can be found in [8], and data for Vacoflux was found in Comsol materials' library. The magnetisation curves are shown in Fig. 1c.

Modelled Cases
Four cases were considered, as depicted in Fig. 2a. In case 4, layers of both HTS and ferromagnetic material are 100 μm thick, 120 layers in total. In each of the cases, the amount of SC material is the same; hence, stacks in cases 1-3 are 6 mm thick, while in the case of 4-12 mm.
In addition, the effect of ferromagnetic material was determined for several conditions: range of temperatures (30-80 K); anisotropic critical current J c (B, T, θ) vs J c (B, T, θ = 0°); and different width of tape, standard 12 mm and 46 mm wide (currently widest commercially available). Magnetic flux density and total flux 1 mm (3) intrinsic ferromagnetic material and (4) ferromagnetic material layered between HTS layers. The flux density is measured 1 mm above the centre stack surface while the magnetic flux is integrated over the surface above the stack at 1 mm height as shown in b above the stack as shown in Fig. 2b were taken as a merit index for comparison between different cases considered.

Stack with No Modifications
As a comparison, a stack with no modifications or ferromagnetic materials was modelled, and the influence of anisotropic critical current density was determined. Figure 3a shows flux density above the stack for 12-and 46-mm-wide tape with and without considering anisotropy of critical current density. The figure shows that the anisotropy does not change the maximum flux density of the stack significantly, especially for 46mm-wide tape; however, Fig. 3b shows that the flux density profile above the stack is changed; hence, the total flux above the stack will also be different. Figure 3b shows that the flux density near the edge of the stack is enhanced, which corresponds to region of the stack where flux density within the stack is closer to θ = 90°, corresponding to a peak in J c (B, T, θ). By considering the anisotropy, the enhancement in flux above the stack for 46-mm tape ranges from 2% at 80 K to 10% at 30 K. This shows that the anisotropic critical current density plays a more important role at low temperatures, which is consistent with J c (B, T, θ) curves in Fig. 1. The results obtained for the stack with no modifications will be used for comparison when ferromagnetic material is added.

Stack with External Ferromagnetic Material
An externally added ferromagnetic material can change field distribution inside and outside the HTS stack. A simple geometry depicted in Fig. 4a was used as a case study to determine optimal shape of the externally added ferromagnetic material, by changing dimensions w 1 , w 2 , h 1 and h 2 . The optimisation was done using a Comsol optimisation module with constraint that h 2 < 6 mm and h 1 , w 1 < 92 mm, so that the top surface of the stack is clear as it would be for most applications and the width and thickness of the FM backplate no larger than twice the stack width-92 mm for 46-mm-wide tape.
It was found that the region of FM material, defined by w 2 , h 2 , does not contribute to increase the flux density above the stack; on the other hand, it was beneficial to continuously increase the w 1 and h 1 . However, as shown in Fig. 5b for 46-mm-diameter stack at 30 K, increasing dimensions of the FM material beyond h 1 , w 1 > 50 mm yields very little improvement. The only exception where the optimal shape was different was the case of 46-mm tape at 80 K, where shape of the J c (B, T, θ) leads to the fact that the dimensions w 1 = 15 mm, h 1 = 92 mm produced slightly higher field compared to h 1 = w 1 = 92 mm (maximum allowed dimensions in the study).Z The improvements of using external FM material with optimised dimensions are detailed in Fig. 5. The results show that stacks made of wide tape benefit much more from the ferromagnetic material, and the improvement is close to 15% at 30 K. While Fig. 4b shows the magnetic field distribution in and around the stack with and without FeCo present. It is worth noting that changing the n-value in the model with fixed h 1 , w 1 = 92 mm in the range n = 10-30 modified the absolute trapped field (higher n value leading to larger field), but the improvement with respect to case with no FM component at 30 K was maintained at 15-18%.

Ferromagnetic Substrate and Layered Stack
Unlike in the case for external ferromagnetic material, it was found that internal ferromagnetic material tends to decrease the flux density outside the material as more flux density is contained within the composite stack itself. Effect is similar for both the layered composite (see Fig. 6b) and a stack where the tape itself is ferromagnetic (homogenised stack with magnetic properties and carrying engineering critical current, akin to the case of ferromagnetic substrate), as shown in Fig. 6a. The flux density outside the stack also decreases with increasing saturation magnetisation of the ferromagnetic material used. Both the current density and magnetic flux density are shown in Fig. 7 for cases 1 and 3 (see Fig. 2a), showing the increased flux density within the stack but not outside it, where it is useful for applications. A lot of detailed theoretical work on magnetisation of thin strips and layered Fig. 6 Flux density above the stacks with and without magnetic properties corresponding to Ni-W and FeCo alloys for 12-and 46-mm stacks (a) and similarly a field above a 60/60 alternating stack of HTS tape and FeCo layers, each 100 μm thick (b) superconductors was done by Brandt [9,10], where a more rigorous time-dependant treatment with comparison of different E-J laws can be found.

Conclusions
The results have shown that a substantial enhancement of field above a trapped field magnet is possible by the use of external ferromagnetic material even at flux densities well above the saturation magnetisation of the ferromagnetic material in question. However, internally added ferromagnetic material by layering or otherwise was detrimental to flux density outside the TFM, and the reduction of flux density was worse with ferromagnetic materials with higher saturation magnetisation.
Lastly, it was determined that it is important to consider the anisotropy of critical current density as it changes the flux density profile and the total amount of magnetic flux above the TFM, although it did not change the maximum flux density above the centre of the stack.
Further work is needed to determine the effect of n-value dependence on temperature and magnetic field and timedependant simulations should shed light on how flux creep is affected by presence of FM. Additional optimisation, especially moving beyond the simple geometry of the FM and tailored for specific J c (B, T, θ) of the superconductor, could lead to higher trapped fields.