Numerical modelling of high temperature superconducting tapes and cables

This Ph.D. thesis is focused on the numerical modelling of high-Tc superconductors (HTS) at the operating temperature of 77 K (liquid nitrogen). The purpose of numerical modelling is to precisely calculate the current and field distributions inside HTS devices (tapes, cables) and the corresponding AC losses, which are one of the most important limiting factors for a large-scale application of such materials. From the electrical point of view, superconductors are characterized by a strongly non-linear voltage-current relation, which defines the transition from superconducting to normal state. In the case of HTS, the steepness of this transition is smoother than for low-Tc superconductors (LTS), so that the commonly used critical state model (CSM) gives a too simplified representation of their electromagnetic behaviour and can be used for a qualitative description only. In this thesis the finite element method (FEM) has been used for precisely computing the current and field distributions as well as for evaluating the AC losses in HTS devices. The superconducting transition is modelled with a power-law relation, E(J) = Ec(J/Jc)n, which is derived from the fit of transport measurements. Firstly, the results obtained with the software package FLUX3D on multi-filamentary tapes have been validated by means of a comparison with the ones obtained by another software package (FLUX2D). The results from FLUX2D, having already been successfully compared with experimental measurements within the framework of two previous Ph.D. theses at LANOS, have been used as reference. Secondly, two power-law models, which take into account the spatial variation of the critical current density inside HTS tapes and its strongly anisotropic dependence on the magnetic field, have been implemented in FLUX3D. In most cases, this latter dependence is extremely important, since the transport capacity of the superconductor is considerably reduced (and its power loss sensibly increased) by the presence of a magnetic field. Afterwards, FEM modelling of HTS tapes has been extended to cables. HTS cables are in general composed by different layers of several tapes and have a quite complex three-dimensional structure: in fact, the layers are wound around a central cylindrical support with different pitch lengths and relative winding orientations, in order to obtain a uniform repartition of the transport current among the layers, which minimizes the total AC losses. For overcoming the difficulties of a direct 3D FEM simulation, a simple electrical model, which allows to find the optimal pitch lengths and whose results are the input data for a 2D FEM evaluation of the AC losses, has been developed. FEM computations have also been used to investigate the influence of the non-uniformity of the tape properties (contact resistance, critical current, power index) on the global loss performance of a single-layer HTS cable. As an alternative to FEM computations, an equivalent circuit model of HTS cables has been utilized. It describes the cable from the macroscopic point of view and allows to compute the current repartition among the layers and the corresponding AC losses, without the necessity of having detailed information about individual tapes. In the framework of the European project BIG-POWA, I have collaborated to the assembling process of a HTS power-link at the Pirelli Labs, in the Milan region, Italy. Finally, 3D simulations have been performed to extensively study the coupling effect between two superconducting filaments via the resistive matrix, which is a typical case where analytical solutions exist for very peculiar geometries and physical conditions only. FEM simulations have been utilized to study the dependence of the filament coupling on the physical and geometrical parameters of the conductors.

Quench behavior of high-temperature superconductor tapes for power applications: a strategy toward resilience

Nicolò Riva

High-Temperature Superconductors (HTS) can be superconducting in liquid nitrogen 77 K, holding immense promises for our future. They can enable disruptive technologies such as nuclear fusion, lossless power transmission, cancer treatment devices, and technologies for future transportation. In the past years, the numerical models to describe the electrical resistivity of REBCO commercial tapes for devices working near and above the critical current, have been shown to be not accurate or very empirical. The resistivity in this regime, in fact, is not very well known. The lack of this knowledge is a significant issue in developing quality simulation tools. The major challenge in retrieving such properties lies in the fact that when I>Ic, heating effects, and thermal instabilities can quickly destroy the conductor if nothing is done to protect it. Moreover, due to the current sharing between the layers, it is difficult to know the amount of current carried by the superconducting layer and its resistivity. The present work aims to understand better the overcritical current regime combining ultra-fast pulsed current measurements performed on HTS REBCO based coated conductors with Finite Element Modeling. The experimental activities were carried out mostly at EPFL and in part at PM and KIT. The modeling activities were carried out between EPFL and KIT. The major result is a resistivity relationship describing the overcritical current regime to be used in numerical simulations of REBCO tapes. The first part of the thesis illustrates a post-processing method based on the so-called Uniform Current (UC) model to estimate the REBCO material's resistivity in the overcritical from experimental measurements. Pulsed current measurements as short as 15 us and with current magnitude up to 5 Ic were performed in liquid nitrogen bath 77 K on samples from various manufacturers, without damaging the tapes. The second part of the thesis discusses a post-processing method based on regularization of data to treat the experimental measurements extracted in the overcritical current regime. The output of this technique is a look-up table that can be shared with interested partners and used in numerical modeling afterward. The third part of the thesis presents the overcritical current model (rho-\eta\beta): a mathematical relationship of the overcritical current regime based on measurements performed between 77 K and 90 K and in self-field conditions. We compare such models with the power-law model, and we provide a short discussion of the fitting parameters and their typical values. The last part of the thesis discusses the overcritical current model, based on experimental measurements obtained as outlined above. The model was validated experimentally and used to show that for the case of a superconducting fault current limiter when the power-law model is used to model its electro-thermal response, the device quenches faster than with the overcritical model. In conclusion, this work can help optimize the use of superconductors and, consequently, the stabilizer. More interestingly, it opens the study of the overcritical current regime, a new exciting aspect of REBCO commercial tapes. This project has received funding from the Swiss Federal Office of Energy SFOE under grant agreement SI/500193-02.

EPFL2021

Numerical Modelling of a HTS Cable

Bertrand Dutoit, Francesco Grilli

The finite element method has been employed in order to simulate the behaviour of a HTS cable for transport current applications. The E-J model includes an anisotropic dependence of the critical current density Jc and power index n on the local magnetic field, whose magnitude is non-negligible and determines the effective critical current of the cable. The cable consists of four electrically insulated layers, each of them composed by 20 superconducting tapes, providing a total critical current of about 4 kA. In order to obtain a uniform repartition of the current among the layers, a sufficiently high contact resistance has been inserted in the electric circuit. The ac losses and the field distribution have been computed. A comparison is also made with a simple electric model of a HTS cable.

2003

Finite element method analysis of the coupling effect between superconducting filaments of different aspect ratio

Bertrand Dutoit, Francesco Grilli, Yong Yang

Using 3D finite element method (FEM) modelling, the present work investigates the coupling effect between two finite superconductors through a resistive matrix. This effect, related to the finite length of the conductors, is typically three-dimensional and cannot be analyzed by the widely utilized 2D models. Superconductors are modelled with the non linear power law E=Ec(J/Jc)^n, which has been implemented in the FEM software. The main focus is to demonstrate the feasibility of such calculations and to establish the correlation between the coupling effect and the aspect ratio of the conductor cross-section, in order to extend the existing theory, which is only precise for superconductors of infinite slabs, fully penetrated by the magnetic field. The effect of other parameters such as the conductor length and the gap between the superconducting filaments is also considered. The latter is a parameter which does not feature in the approximate theories.

2003

Quench behavior of high-temperature superconductor tapes for power applications: a strategy toward resilience

Nicolò Riva

EPFL2021