Search

The development of the high-temperature superconductors (HTS) conductors has allowed the development of diverse superconductor devices. Some of these devises, like the power generators and high-field magnets, are classified as large-scale HTS systems, because they are made of hundreds or thousands of turns. Mathematical models are required to address the analysis of these kind of systems. This task cannot be done by means of analytical models, because they are limited to the analysis of simple assemblies. The finite-element models using the H formulation have been extensively used during the last years. Nevertheless, the use of H formulation models to analyze large-scale systems is hindered by the excessive computational load. The recently proposed T-A formulation models have allowed building more efficient models for systems made of HTS tapes. Additionally, the homogenization and multi-scaling methods have been successfully applied in conjunction with the H and T-A formulations, these simplification methods allows reducing the required computational resources. In this article a new simplification method, called densification, is proposed. The strategies emerging from the combined use of the formulations and the simplification methods already mentioned are extensively explored, and the comprehensive validation and comparison of all the resulting strategies is presented.

In recent years, the $H$ formulation of Maxwell's equation has become the de facto standard for simulating the time-dependent electromagnetic behavior of superconducting applications with commercial software. However, there are cases where other formulations are desirable, for example for modeling superconducting turns in electrical machines or situations where the superconductor is better described by the critical state than by a power-law resistivity. In order to accurately and efficiently handle those situations, here we consider two published approaches based on the magnetic vector potential: the $T$-$A$ formulation of Maxwell's equations (with power-law resistivity) and Campbell's implementation of the critical state model. In this contribution, we extend the $T$-$A$ formulation to thick conductors so that large coils with different coupling scenarios between the turns can be considered. We also revise Campbell's model and discuss it in terms of its ability to calculate AC losses: in particular, we investigate the dependence of the calculated AC losses on the frequency of the AC excitation and the possibility of using quick one-step (instead of full cycle) simulations to calculate the AC losses.

The emergence of second generation (2G) high-temperature superconductor (HTS) tapes has favored the development of HTS magnets for their applications in areas such as NMR, MRI and high field magnets. The screening current-induced field and the field drift are two major problems hindering the use of HTS tapes in the mentioned areas. Both problems are caused by the screening current, then it is necessary to have a modeling strategy capable to estimate such phenomena. Thus far, the H formulation has been the most successful and used approach to model medium-size systems (hundreds of tapes). However, its application to large-scale systems is still impaired by excessive computation times and memory requirements. Homogenization and multi-scaling strategies have been successfully implemented to increase the computational efficiency. In this contribution, we show that using the homogenization technique with the recently developed T-A formulation allows reducing the computation time and the amount of memory up to the point that real-time simulations of slow ramping cycles of large-scale systems are possible. The T-A homogeneous model also allows systematically investigating the screening current using numerical simulations.

The emergence of second-generation high temperature superconducting tapes has favored the development of large-scale superconductor systems. The mathematical models capable of estimating electromagnetic quantities in superconductors have evolved from simple analytical models to complex numerical models. The available analytical models are limited to the analysis of single wires or infinite arrays that, in general, do not represent real devices in real applications. The numerical models based on finite element method using the H formulation of the Maxwells equations are useful for the analysis of medium-size systems, but their application in large-scale systems is problematic due to the excessive computational cost in terms of memory and computation time. Then it is necessary to devise new strategies to make the computation more efficient. The homogenization and the multi-scale methods have successfully simplified the description of the systems allowing the study of large-scale systems. Also, efficient calculations have been achieved using the T-A formulation. In the present work, we propose a series of adaptations to the multi-scale and homogenization methods so that they can be efficiently used in conjunction with the T-A formulation to compute the distribution of current density and hysteresis losses in the superconducting layer of superconducting tapes. The computation time and the amount of memory are substantially reduced up to a point that it is possible to achieve real-time simulations of HTS large-scale systems under slow ramping cycles of practical importance on personal computers.

The emergence of second-generation high temperature superconducting tapes has favored the development of large-scale superconductor systems. The mathematical models for superconductors have evolved from simple analytical models to complex numerical models. The available analytical models are just capable of estimating electromagnetic quantities in single wires or small arrays. The numerical models like the H formulation finite element models are useful for the analysis of more medium-size systems, but their application in large-scale systems is problematic due to the excessive cost in terms of memory and the computation time. Then it is necessary to devise new strategies to tackle this challenge. The homogenization and the multi-scale approach are methodologies that have successfully simplified the description of the systems allowing the study of large-scale systems. Recently, efficient calculations have been achieved using the T-A formulation. In this manuscript, we analyze the adequate order of the elements used in the T-A formulation. More importantly, we propose a series of adaptations to the homogenization and multi-scale methodologies in order to allow their implementation to the T-A formulation models. The proposed strategies are validated by means of comparing them with the H reference model of a racetrack coil. The computation time and memory demand are substantially reduced up to a point that makes it possible to perform real-time simulations for slow ramping cycles. At the same time the comparison demonstrates a very good agreement with respect to the reference model, both at a local and global level, is achieved., Comment: arXiv admin note: substantial text overlap with arXiv:1811.09193

The modeling of superconducting magnetic bearing (SMB) is of great significance for predicting and optimizing its levitation performance before construction. Although lots of efforts have been made in this area, it still remains some space for improvements. Thus the goal of this work is to report a flexible, fast and trustworthy H-formulation finite element model. First the methodology for modeling and calibrating both bulk-type and stack-type SMB is summarized. Then its effectiveness for simulating SMBs in 2-D, 2-D axisymmetric and 3-D is evaluated by comparison with measurements. In particular, original solutions to overcome several obstacles are given: clarification of the calibration procedure for stack-type and bulk-type SMBs, details on the experimental protocol to obtain reproducible measurements, validation of the 2-D model for a stack-type SMB modeling the tapes real thickness, implementation of a 2-D axisymmetric SMB model, implementation of a 3-D SMB model, extensive validation of the models by comparison with experimental results for field cooling and zero field cooling, for both vertical and lateral movements. The accuracy of the model being proved, it has now a strong potential for speeding up the development of numerous applications including maglev vehicles, magnetic launchers, flywheel energy storage systems, motor bearings and cosmic microwave background polarimeters., Comment: 22 pages, 15 figures

In recent years, commercial HTS superconductors have gained an increasing interest for their use in applications involving large-scale superconductor systems. These systems are typically made from hundreds to thousands of turns of conductors. Due to the large number of turns, the simulations of a whole system can become prohibitive in terms of computing time and load. Therefore, an efficient strategy which does not compromise the accuracy of calculations is needed. Recently, a method, based on a multi-scale approach, showed that the computational load can be lowered by simulating, in detail, only several significant tapes from the system. The main limitation of this approach is the inaccuracy of the estimation of the background magnetic field. To address this issue, we consider the following two complementary strategies. The first strategy consists in the iterative implementation of the multi-scale method. The multi-scale method solves itself a dynamic problem, the iterative implementation proposed here is the iterative application of the multi-scale method, and a dynamic solution is obtained at each iteration. The second strategy is a new interpolation method for current distributions, based on the inverse cumulative density function interpolation technique. With respect to conventional interpolation methods, a more realistic current density distribution is then obtained, which allows for a better estimation of the background magnetic field, and consequently, a better estimation of the hysteresis losses. In contrast with previous works, here we do not focus only on the estimation of the hysteresis losses, but also the estimation of background field and the current density distribution is addressed.

The estimation of the critical current ($I_c$) and AC losses of high-temperature superconductor (HTS) devices through modeling and simulation requires the knowledge of the critical current density ($J_c$) of the superconducting material. This $J_c$ is in general not constant and depends both on the magnitude ($B_{loc}$) and the direction ($\theta$, relative to the tape) of the local magnetic flux density. In principle, $J_c(B_{loc},\theta)$ can be obtained from the experimentally measured critical current $I_c(B_a,\theta)$, where $B_a$ is the magnitude of the applied magnetic field. However, for applications where the superconducting materials experience a local field that is close to the self field of an isolated conductor, obtaining $J_c(B_{loc},\theta)$ form $I_c(B_a,\theta)$ is not a trivial task. It is necessary to solve an inverse problem to correct for the contribution derived from the self field. The methods presented in the literature comprise a series of approaches dealing with different degrees of mathematical regularization, such as the use of brute force or optimization methods to fit the parameters of preconceived non linear formulas. In this contribution, we present a parameter-free method that provides excellent reproduction of experimental data and requires no human interaction or preconception of the $J_c$ dependence with respect to the magnetic field. In particular, it allows going from the experimental data to a ready-to-run $J_c(B_{loc},\theta)$ model in a few minutes.

High temperature superconducting (HTS) bulks or stacks of coated conductors (CCs) can be magnetized to become trapped field magnets (TFMs). The magnetic fields of such TFMs can break the limitation of conventional magnets (<2 T), so they show potential for improving the performance of many electrical applications that use permanent magnets like rotating machines. Towards practical or commercial use of TFMs, effective in situ magnetization is one of the key issues. The pulsed field magnetization (PFM) is among the most promising magnetization methods in virtue of its compactness, mobility and low cost. However, due to the heat generation during the magnetization, the trapped field and flux acquired by PFM usually cannot achieve the full potential of a sample (acquired by the field cooling or zero field cooling method). The multi-pulse technique was found to effectively improve the trapped field by PFM in practice. In this work, a systematic study on the PFM with successive pulses is presented. A 2D electromagnetic-thermal coupled model with comprehensive temperature dependent parameters is used to simulate a stack of CCs magnetized by successive magnetic pulses. An overall picture is built to show how the trapped field and flux evolve with different pulse sequences and the evolution patterns are analyzed. Based on the discussion, an operable magnetization strategy of PFM with successive pulses is suggested to provide more trapped field and flux. Finally, experimental results of a stack of CCs magnetized by typical pulse sequences are presented for demonstration.

The ability of stacks of superconducting tapes to trap large magnetic fields makes them ideal candidates for creating powerful permanent magnets of compact size and mass. Experimentally, several techniques are used to trap the maximum possible field in a given practical application. However, regardless of the magnetization method used, there is a physical limit to the maximum magnetic field that a given superconducting magnet can trap. This limit is given by the geometric design, the particular superconducting material used and the temperature of operation. Knowing the maximum possible trapped field is important for device design as it provides an upper limit for applications such as magnetic bearings or rotating machinery. In this work we present a collection of finite element method (FEM) models in 2D and 3D capable of estimating the maximum trapped field of stacked tape superconducting magnets. The models are computationally fast and can be used to perform parametric studies with ease. For the case of square stacks tape magnets, various sizes are considered and their estimated maximum trapped field is compared with experimental results.

In this work we present a modeling tool designed to estimate the hysteretic losses in the coils of an electric generator with coils made of coated conductor tapes during transient operation. The model is based on a two-stage segregated model approach that allows simulating the electric generator and the current distribution in the superconducting coils using a one-way coupling from the generator to the HTS coils model. The model has two inputs: the rotational speed and the electric load signal. A homogeneous anisotropic bulk model for the coils allows computing the current distribution in the coils. From this distribution, the hysteretic losses are estimated. Beyond the interest on providing an estimate on the global energy dissipation in the machine, in this work we present a more detailed local analysis that allows addressing issues such as coil design, critical current ratting, electric load change rate limits, cryocooler design, identification of quench-prone regions and overall transient performance.

The design of a HTS SMES coil with solenoidal and toroidal geometry is carried out based on a commercially available 2G HTS conductor. A SMES system of practical interest (1 MW / 5 s) is considered. The comparison between ideal toroidal and solenoidal geometry is first discussed and the criteria used for choosing the geometrical parameters of the coils' bore are explained. The design of the real coil is then carried out and the final amount of conductor needed is compared. A preliminary comparison of the two coils in terms of AC loss during one charge discharge cycle is also discussed.

High temperature superconducting (HTS) bulks can be magnetized to become powerful trapped field magnets (TFMs), which are promising for high-performance electrical applications. To magnetize such TFMs, pulsed field magnetization (PFM) is supposed to substitute field cooling (FC) to provide in situ magnetization. However, the heat generation during PFM, which reduces the trapped field has always been an issue, so numerical simulation of the process is important to provide optimal magnetization strategies. In this paper, HTS bulks magnetized by PFM are simulated with an axisymmetric electromagnetic-thermal coupled model based on H-formulation. Influences of important yet difficult-to-characterize parameters of HTS bulks including n values in the E-J Power Law and B0 in the Kim's Law are investigated. Furthermore, controlled magnetic density distribution coils (CMDCs) which generate a non-uniform field are suggested to further improve the trapped field of HTS bulks compared to split coils proposed previously.

Numerical models are powerful tools to predict the electromagnetic behavior of superconductors. In recent years, a variety of models have been successfully developed to simulate high-temperature-superconducting (HTS) coated conductor tapes. While the models work well for the simulation of individual tapes or relatively small assemblies, their direct applicability to devices involving hundreds or thousands of tapes, as for example coils used in electrical machines, is questionable. Indeed the simulation time and memory requirement can quickly become prohibitive. In this article, we develop and compare two different models for simulating realistic HTS devices composed of a large number of tapes: 1) the homogenized model simulates the coil using an equivalent anisotropic homogeneous bulk with specifically developed current constraints to account for the fact that each turn carries the same current; 2) the multi-scale model parallelizes and reduces the computational problem by simulating only several individual tapes at significant positions of the coil's cross-section using appropriate boundary conditions to account for the field generated by the neighboring turns. Both methods are used to simulate a coil made of 2000 tapes, and compared against the widely used H-formulation finite element model that includes all the tapes. Both approaches allow speeding-up simulations of large number of HTS tapes by 1-3 orders of magnitudes, while keeping a good accuracy of the results. Such models could be used to design and optimize large-scale HTS devices., Comment: 17 pages, 11 figures

Considering the need for cost/performance prediction and optimization of superconducting maglev vehicles, we develop and validate here a 3D finite element model to simulate superconducting linear magnetic bearings. Then we reduce the 3D model to a 2D model in order to decrease the computing time. This allows us to perform in a reasonable time a stochastic optimization considering the superconductor properties and the vehicle operation. We look for the permanent magnet guideway geometry that minimizes the cost and maximizes the lateral force during a displacement sequence, with a constraint on the minimum levitation force. The displacement sequence reproduces a regular maglev vehicle operation with both vertical and lateral movements. For the sake of comparison, our reference is the SupraTrans prototype bearing. The results of the optimization suggest that the bearing cost could be substantially reduced, while keeping the same performances as the initial design. Alternatively, the performances could be significantly improved for the same original cost.

In order to transport sufficiently high current, high-temperature superconductor (HTS) tapes are assembled in cable structures of different forms. In such cables, the tapes are tightly packed and have a strong electromagnetic interaction. In particular, the generated self-field is quite substantial and can give an important contribution in reducing the maximum current the cable can effectively carry. In order to be able to predict the critical current of said cable structures, a static numerical model has been recently proposed. In this contribution, we present in detail the implementation of such models in different programming environments, including finite-element-based and general numerical analysis programs, both commercial an open-source. A comparison of the accuracy and calculation speed of the different implementations of the model is carried out for the case of a Roebel cable. The model is also used to evaluate the importance of choosing a very accurate description of the angular Jc(B) dependence of the superconductor as input for the material's property. The numerical codes, which are open-source, are made freely available to interested users., Comment: The accompanying models descried in this work are available for free download at the HTS modeling workgroup under http://www.htsmodelling.com/?p=917

Nowadays, there is growing interest in using superconducting wires or tapes for the design and manufacture of devices such as cables, coils, rotating machinery, transformers and fault current limiters among others. Their high current capacity has made them the candidates of choice for manufacturing compact and light cables and coils that can be used in the large scale power applications described above. However, the performance of these cables and coils is limited by their critical current, which is determined by several factors, including the conductor's material properties and the geometric layout of the device itself. In this work we present a self-consistent model for estimating the critical current of superconducting devices. This is of large importance when the operating conditions are such that the self-field produced by the current is comparable to the overall background field. The model is based on the asymptotic limit when time approaches infinity of Faraday's equation written in terms of the magnetic vector potential. It uses a continuous E-J relationship and takes the angular dependence of the critical current density on the magnetic flux density into account. The proposed model is used to estimate the critical current of superconducting devices such as cables, coils, and coils made of transposed cables with very high accuracy. The high computing speed of this model makes it an ideal candidate for design optimization.

Different magnet applications require compact high current cables. Among the proposed solutions, the Twisted Stacked Tape Cable (TSTC) is easy to manufacture and has very high tape length usage efficiency. In this kind of cables the tapes are closely packed, so that their electromagnetic interaction is very strong and determines the overall performance of the cable. Numerical models are necessary tools to precisely evaluate this interaction and to predict the cable's behavior, e.g. in terms of effective critical current and magnetization currents. For this purpose, we developed a fully three-dimensional model of a TSTC, which not only takes into account the twisted geometry of these cables, but is also able to account for the contact resistances of the current terminations. The latter can have profound influence on the way the current is partitioned among the tapes, especially on short laboratory prototypes. In this paper, we first use the numerical model to compute the critical current and the magnetization AC loss of a twisted tape, showing the differences with the case of a straight tape. Then, we use it to calculate the current distribution in a TSTC cable, comparing the results with those experimentally obtained on a cable composed of four straight stacked tapes. The results show the ability of the model to simulate twisted conductors and constitutes a first step toward the simulation of TSTC in high-field magnet applications. The presented modeling approach is not restricted to the TSTC geometry, but may be used for any cable configuration with periodical translational symmetry.

Numerical models for computing the effective critical current of devices made of HTS tapes require the knowledge of the Jc(B,theta) dependence, i.e. of the way the critical current density Jc depends on the magnetic flux density B and its orientation theta with respect to the tape. In this paper we present a numerical model based on the critical state with angular field dependence of Jc to extract the Jc(B,theta) relation from experimental data. The model takes into account the self-field created by the tape, which gives an important contribution when the field applied in the experiments is low. The same model can also be used to compute the effective critical current of devices composed of electromagnetically interacting tapes. Three examples are considered here: two differently current rated Roebel cables composed of REBCO coated conductors and a power cable prototype composed of Bi-2223 tapes. The critical currents computed with the numerical model show good agreement with the measured ones. The simulations reveal also that several parameter sets in the Jc(B,theta) give an equally good representation of the experimental characterization of the tapes and that the measured Ic values of cables are subjected to the influence of experimental conditions, such as Ic degradation due to the manufacturing and assembling process and non-uniformity of the tape properties. These two aspects make the determination of a very precise Jc(B,theta) expression probably unnecessary, as long as that expression is able to reproduce the main features of the angular dependence. The easiness of use of this model, which can be straightforwardly implemented in finite-element programs able to solve static electromagnetic problems, is very attractive both for researchers and devices manufactures who want to characterize superconducting tapes and calculate the effective critical current of superconducting devices.

Use of 2G HTS coated conductors in several power applications has become popular in recent years. Their large current density under high magnetic fields makes them suitable candidates for high power capacity applications such as stacks, coils, magnets, cables and current leads. For this reason, modeling and simulation of their electromagnetic properties is very desirable in the design and optimization processes. For many applications, when symmetries allow it, simple models consisting of 1D or 2D representations are well suited for providing a satisfying description of the problem at hand. However, certain designs such as racetrack coils and finite-length or non-straight stacks, do pose a 3D problem that cannot be easily reduced to a 2D configuration. Full 3-D models have been developed, but their use for simulating superconducting devices is a very challenging task involving a large-scale computational problem. In this work, we present a new method to simulate the electromagnetic transient behavior of 2G HTS stacks and coils. The method, originally used to model stacks of straight superconducting tapes or circular coils in 2D, is now extended to 3D. The main idea is to construct an anisotropic bulklike equivalent for the stack or coil, such that the geometrical layout of the internal alternating structures of insulating, metallic, superconducting and substrate layers is reduced while keeping the overall electromagnetic behavior of the original device. Besides the aforementioned interest in modeling and simulating 2G HTS coated conductors, this work gives a further step towards efficient 3D modeling and simulation of superconducting devices for large scale applications.

In this contribution we develop a refined numerical model of pancake coils assembled from a coated conductor Roebel cable, which includes the angular dependence of the critical current density $J_c$ on the magnetic field and the actual (three-dimensional) shape of the current lead used to inject the current. Previous works of ours indicate that this latter has an important influence on the measured value of the AC losses. For the simulation of the superconductor, we used two alternative models based on different descriptions of the superconductor's properties and implemented in different mathematical schemes. For the simulation of the current lead we use a full three-dimensional finite-element model. The results of the simulation are compared with measurements and the main issues related to the modeling and the measurement of Roebel coils are discussed in detail.

A homogenization method to model a stack of second generation (2G) High Temperature Superconducting (HTS) tapes under AC applied transport current or magnetic field has been obtained. The idea is to find an anisotropic bulk equivalent for the stack, such that the geometrical layout of the internal alternating structures of insulating, metallic, superconducting and substrate layers is "washed" out while keeping the overall electromagnetic behavior of the original stack. We disregard assumptions upon the shape of the critical region and use a power law E-J relationship allowing for overcritical current densities to be considered. The method presented here allows for a computational speedup factor of up to 2 orders of magnitude when compared to full 2-D simulations taking into account the actual dimensions of the stacks without compromising accuracy.

First experimental observations of the ferromagnetic shielding effect in high-Tc superconducting coated conductors were carried out. Experimental results were compared to simulations calling upon finite-element calculations based on the H-formulation of Maxwell equations to model superconducting strips with ferromagnetic shields. Samples of copper-stabilized coated conductors were electroplated with nickel shields and afterwards characterized. Both externally applied oscillating transverse magnetic fields as well as transport currents were studied. Having observed promising gains with respect to the reduction of ac losses in both cases, we further investigated the potential of ferromagnetic shielding. The numerical model was able to reproduce and also predict experimental results very well and will serve as an indispensable tool to determine the potential of soft ferromagnetic materials to significantly reduce hysteretic losses., Comment: 4 pages, 3 figures, 1 table