Central Solenoid Winding Pack Design

Similar documents
High Temperature Superconductors for Fusion at the Swiss Plasma Center

Magnetic Flux Surface Measurements at Wendelstein 7-X

Modelling and Analysis of the JET EP2 Neutral Beam Full Energy Ion Dump Curved End Plate

Temporal Evolution of Temperature and Argon Impurity Density Profiles Observed by X-ray Imaging Spectrometer Measurements at Wendelstein 7-X

JOINTS FOR SUPERCONDUCTING MAGNETS

Electron energy distribution function in the divertor region of the COMPASS tokamak during neutral beam injection heating

EUROFUSION WPJET1-PR(16) CG Albert et al.

The emissivity of W coatings deposited on carbon materials for fusion applications

Qualification Tests and Facilities for the ITER Superconductors

Integrated equilibrium reconstruction and MHD stability analysis of tokamak plasmas in the EU-IM platform

Power Balance Analysis of Wendelstein 7-X Plasmas Using Profile Diagnostics

Thermo-mechanical analyses and ways of optimization of the helium cooled DEMO First Wall under RCC-MRx rules

Development of cable in conduit conductor for ITER CS in Japan

Material, Design, and Cost Modeling for High Performance Coils. L. Bromberg, P. Titus MIT Plasma Science and Fusion Center ARIES meeting

Comparative Transport Analysis of JET and JT-60U Discharges

Turbulent Transport Analysis of JET H-mode and Hybrid Plasmas using QuaLiKiz, TGLF and GLF23

Thermal-hydraulic analysis of LTS cables for the DEMO TF coil using simplifi ed models*

High Temperature Superconductor. Cable Concepts for Fusion Magnets. Christian Barth. \iyit Scientific. ^VI I Publishing

Preliminary system modeling for the EUROfusion Water Cooled Lithium Lead blanket

CFD simulation of the magnetohydrodynamic flow inside the WCLL breeding blanket module

Common approach for thermal-hydraulic calculations

Large Superconducting Conductors and Joints for Fusion Magnets : from Conceptual Design to Test at Full Size Scale

Multiscale modelling of sheath physics in edge transport codes

In situ wavelength calibration of the edge CXS spectrometers on JET

Divertor configuration with two nearby poloidal field nulls: modelling and experiments for EAST and JET tokamaks

Test Results of ITER Conductors in the SULTAN Facility

Which Superconducting Magnets for DEMO and Future Fusion Reactors?

POLITECNICO DI TORINO Repository ISTITUZIONALE

Innovative fabrication method of superconducting magnets using high T c superconductors with joints

Multi-Machine Experiments to Study the Parametric Dependences of Momentum Transport

IEEE TRANSACTIONS ON APPLIED SUPERCONDUCTIVITY, VOL. 14, NO. 2, JUNE A. Nijhuis, Yu. Ilyin, W. Abbas, B. ten Haken, and H. H. J.

Modelling of prompt deposition of tungsten under fusion relevant conditions

Development of miniaturized, spectroscopically assisted Penning gauges for fractional helium and hydrogen neutral pressure measurements

EuCARD-2 Enhanced European Coordination for Accelerator Research & Development. Journal Publication

Massachusetts Institute of Technology 22.68J/2.64J Superconducting Magnets. February 27, Lecture #4 Magnetic Forces and Stresses

Self Field Measurements by Hall Sensors on the SeCRETS Long Sample CICCs in SULTAN

JT-60 SA Toroidal Field coil structural analysis

Studies of Next-Step Spherical Tokamaks Using High-Temperature Superconductors Jonathan Menard (PPPL)

High Temperature Superconductor Cables for DEMO TF-Magnets

Superconducting Magnet Design and R&D with HTS Option for the Helical DEMO Reactor

Flux Motion and Screening Current in High-temperature Superconducting Magnets

Evaluation of the Current Sharing Temperature of the ITER Toroidal Field Model Coil

Critical Current Properties of HTS Twisted Stacked-Tape Cable in Subcooled- and Pressurized-Liquid Nitrogen

Use of High Temperature Superconductors for Future Fusion Magnet Systems

Physics and Engineering Studies of the Advanced Divertor for a Fusion Reactor

Electromagnetic, stress and thermal analysis of the Superconducting Magnet

Design of next step tokamak: Consistent analysis of plasma flux consumption and poloidal field system

A SUPERCONDUCTING TOKAMAK FUSION TRANSMUTATION OF WASTE REACTOR

Flux-driven multi-channel simulations with the quasilinear gyrokinetic tokamak transport model QuaLiKiz

Evaluation of First Wall Heat Fluxes Due to Magnetic Perturbations for a Range of ITER Scenarios

Thermal optimization of the Helium-Cooled Lithium Lead breeding zone layout design regarding TBR enhancement

Computational Fluid Dynamic analysis of Screw tube relevant for fusion applications

Superconducting Magnet with a Minimal Steel Yoke for the Future Circular Collider Detector

Design Principles of Superconducting Magnets

Simultaneous measurement of critical current, stress, strain and lattice distortions in high temperature superconductors

Inductance and Current Distribution Analysis of a Prototype HTS Cable

Small Spherical Tokamaks and their potential role in development of fusion power

Batavia, Illinois, 60510, USA

New European Accelerator Project EuCARD: Work Package on High Field Magnets

1074 IEEE TRANSACTIONS ON APPLIED SUPERCONDUCTIVITY, VOL. 21, NO. 3, JUNE 2011

High Field HTS SMES Coil

The physics of the heat flux narrow decay length in the TCV scrape-off layer: experiments and simulations

The role of statistical noise in edge plasma transport codes based on kinetic Monte Carlo solvers for neutrals: an analogy with turbulent fluctuations

Challenges on demountable / segmented coil concept for high-temperature superconducting magnet

Magnetic and Cryogenic Design of the MICE Coupling Solenoid Magnet System

Continued Developments in High Magnetic Fields Enabled by Second-Generation High- Temperature Superconductors

Simulation of PHA Soft X-Ray Spectra Expected from W7-X

HiLumi LHC FP7 High Luminosity Large Hadron Collider Design Study. Milestone Report. Cryogenic Scenarios for the Cold Powering System

EDIPO Test Facility User Specification

Application of the ECRH radiation for plasma diagnostics in Wendelstein 7-X

Power loads to misaligned edges in COMPASS

Optimisation of ITER Nb 3 Sn CICCs for coupling loss, transverse electromagnetic load and axial thermal contraction

EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH Laboratory for Particle Physics

Al-Ti-MgB 2 conductor for superconducting space magnets

Nuclear Analysis of the HCLL Blanket Concept for the European DEMO Using the TRIPOLI-4 Monte Carlo Code

Experience in manufacturing a large HTS magnet for a SMES

Experiments on Actuator Management and Integrated Control at ASDEX Upgrade

High Temperature Superconductors for Future Fusion Magnet Systems Status, Prospects and Challenges

Lecture #2 Design Guide to Superconducting Magnet

Measurements of Deuterium Retention and Surface Elemental Composition with Double Pulse Laser Induced Breakdown Spectroscopy

Review of design aspects for high current Nb3Sn conductors

Towards powerful negative ion beams at the test facility ELISE for the ITER and DEMO NBI systems

MHD limits and plasma response in high beta hybrid operations in ASDEX Upgrade

Neutron Emission Spectroscopy Measurements with a Single Crystal Diamond Detector at JET

HTS Roebel and Rutherford Cables for High-Current Applications

Analysis of Coupled Electromagnetic-Thermal Effects in Superconducting Accelerator Magnets

DEVELOPMENT AND PRODUCTION OF SUPERCONDUCTING AND CRYOGENIC EQUIPMENT AND SYSTEMS FOR ACCELERATORS BY IHEP

DEMO Concept Development and Assessment of Relevant Technologies. Physics and Engineering Studies of the Advanced Divertor for a Fusion Reactor

Recent Developments in YBCO for High Field Magnet Applications

RADIATION EFFECTS ON HIGH TEMPERATURE SUPERCONDUCTORS

Status of the Concept Design of CFETR Tokamak Machine

20 T Block Dipole: Features and Challenges

Test Results of a NbTi Wire for the ITER Poloidal Field Magnets: A Validation of the 2-pinning Components Model

Multiphysics Simulations for the design of a Superconducting magnet for proton therapy

Conceptual Design of CFETR Tokamak Machine

Tailoring the charged particle fluxes across the target surface of Magnum-PSI

Upgrade of the Present JET Shape and Vertical Stability Controller

On the Design Aspects Affecting Performance of GEM based Detector Development for Plasma Diagnostics

1st Performance Test of the 25 T Cryogen-free Superconducting Magnet

Honghai Song. Yawei Wang, Kent Holland, Ken Schrock, Saravan Chandrasekaran FRIB/MSU & SJTU June 2015, SJTU Xuhui Campus

Transcription:

EUROFUSION WPMAG-CP(16) 15681 R Wesche et al. Central Solenoid Winding Pack Design Preprint of Paper to be submitted for publication in Proceedings of 29th Symposium on Fusion Technology (SOFT 2016) This work has been carried out within the framework of the EUROfusion Consortium and has received funding from the Euratom research and training programme 2014-2018 under grant agreement No 633053. The views and opinions expressed herein do not necessarily reflect those of the European Commission.

This document is intended for publication in the open literature. It is made available on the clear understanding that it may not be further circulated and extracts or references may not be published prior to publication of the original when applicable, or without the consent of the Publications Officer, EUROfusion Programme Management Unit, Culham Science Centre, Abingdon, Oxon, OX14 3DB, UK or e-mail Publications.Officer@euro-fusion.org Enquiries about Copyright and reproduction should be addressed to the Publications Officer, EUROfusion Programme Management Unit, Culham Science Centre, Abingdon, Oxon, OX14 3DB, UK or e-mail Publications.Officer@euro-fusion.org The contents of this preprint and all other EUROfusion Preprints, Reports and Conference Papers are available to view online free at http://www.euro-fusionscipub.org. This site has full search facilities and e-mail alert options. In the JET specific papers the diagrams contained within the PDFs on this site are hyperlinked

Central Solenoid Winding Pack Design R. Wesche, N. Bykovsky, X. Sarasola, K. Sedlak, B. Stepanov, D. Uglietti, P. Bruzzone Ecole Polytechnique Fédérale de Lausanne (EPFL), Swiss Plasma Center (SPC),CH 5232 Villigen PSI, Switzerland The present study aims to reduce the outer radius of the central solenoid (CS) with respect to its nominal size specified by EUROfusion for a maintained CS magnetic flux. A reduced outer CS radius would allow the reduction of the overall size and cost of the DEMO magnet system. The proposed outline design of the winding pack for the CS1 module is based on layer winding. To achieve the same magnetic flux in a CS coil of significantly reduced outer radius the peak magnetic field at the CS conductors needs to be substantially increased. The use of high temperature superconductors is therefore envisaged in the highest field sections of the CS coil. It is planned to use react & wind Nb 3 Sn conductors for intermediate field sections and NbTi at the lowest fields. In order to make a most economic use of the superconductors the proposed winding pack design considers a superconductor grading. Keywords: central solenoid, DEMO, high temperature superconductor, superconductor grading, winding pack E-mail address: rainer.wesche@psi.ch

1. Introduction The present CS design study is based on the PROCESS data [1] for the 2015 DEMO reference, still valid in 2016, the CAD model of the full tokamak [2] and the reference flat top equilibria [3]. In a previous study [4], the inner and outer radii of the CS coil were varied for a constant magnetic flux generated by the CS coil. In these calculations, the overall current density was kept constant in all layers of each CS module, while different superconductors were considered for the high, medium, and low field sections. The results [4] indicate that the use of high temperature superconductors (HTS) would provide the possibility to increase the peak magnetic field and to generate the same magnetic flux with a CS coil of reduced outer diameter as compared to the PROCESS data [1]. 2. Reference design The CS coil is composed of five modules, namely CSU3, CSU2, CS1, CSL2 and CSL3 [2]. The inner and outer radii of all CS modules are 2.49 and 3.31 m, respectively. The peak magnetic field of the CS coil includes a contribution generated by the six PF coils and the plasma current, and depends on the considered plasma scenario. In the present work, the equilibrium currents I CSU3 = I CSU2 = I CSL2 = 28.07 MA, I CS1 = 57.14 MA, I CSL3 = 20.14 MA, I PF1 =12.38 MA, I PF2 = 4.63 MA, I PF3 = 3.41 MA, I PF4 = 4.34 MA, I PF5 = 3.20 MA, I PF6 = 19.20 MA, I plasma = 0 for pre magnetization [3] are used as a reference leading to a peak magnetic field of 12.36 T in the conductor and a magnetic flux of 307 Vs, generated only by the 5 CS modules, in the central plane of the CS coil. In further considerations, the inner and outer radii of the CS coil have been varied. A maintained magnetic flux of 307 Vs requires an adjustment of the currents in the CS modules. The ratio of the currents in the different CS modules is kept at the values I CS1 :I CSU3 :I CSL3 = 57.14:28.07:20.18 and I CSU3 = I CSU2 = I CSL2. The field contribution of the PF coils has been calculated with the above quoted equilibrium currents. 3. Superconductor properties The superconducting properties of Nb 3 Sn and NbTi are based on the scaling relations and parameters provided in [5]. In case of Nb 3 Sn, the scaling parameters of the JASTEC bronze strand used in the ITER CSJA6 sample and for the CS insert test in 2015 in Naka have been selected. For the react & wind (R&W) Nb 3 Sn conductors, a strain of 0.25% has been retained. The scaling relation for REBa 2 Cu 3 O 7 x (RE 123, RE = rare earth element) HTS is based on the test of the SPC cable. The non copper critical current density obeys the scaling relation [6] with J c A B B p q B B T 0T B B B T 1 T 0 B0 0 1 T 0 0 where B is the magnetic field and the values of the fitting parameters are A = 15.2 MA m 2 T 1, = 3.48, = 1.61, p = 0.54, q = 2.82, B 0 (0) = 170.78 T. The fitting parameters have no physical meaning. Especially B 0 and T 0 should not be considered as the irreversibility field and temperature. The value of A is based on RE 123 tapes with a non copper thickness of 55 m. 4. Dimensioning of graded CS1 module Studies of uniform current density CS coils suggest that a minimum for the outer CS coil radius can be achieved for an inner coil radius of around 1.9 m [4]. The present work considers 10 sub coils, each composed of two layers. In the magnetic field calculations, the overall current densities in the individual sub coils are used omitting the details of conductor trajectories in the winding pack. Furthermore, it is assumed that the same grading as in the CS1 module is also used for the four other CS modules. The peak magnetic field of each sub coil in the central plane of the CS1 module has been used to estimate the current sharing temperature (T cs ) or the ratio of operation to critical current (I op /I c ). For a maintained magnetic flux, a decreased outer CS coil radius leads to higher peak fields and higher overall winding current densities, and hence a larger fraction of stainless steel is required in the winding pack. The area needed for the electrical insulation depends also on the total thickness of the winding pack and the envisaged number of layers. In several iterative cycles, a solution with a self consistent set of parameters needs to be found. The present work considers N = 10 sub coils, each composed of two layers, and an inner radius of 1.9 m for the innermost sub coil. The current in each sub coil is I CS1 k = I CS1 /N, where I CS1 is the total current in the CS1 module and N the number of sub coils. A copper current density J Cu = 120 MA/m 2 has been selected leading to a copper cross section of A Cu k = I CS1 k /J Cu. For HTS sub coils, a non Cu operation current density J sc_k of 80% of J c (T op, B p ) is envisaged, where T op = 4.75 K and B p is the peak magnetic field of the sub coil in question. In case of Nb 3 Sn and NbTi the operation current density is selected in such a way that T cs = 6.25 K leading to a temperature margin of 1.5 K. The resulting non Cu superconductor area in sub coil k is A sc_k = I CS1 k /J sc k. The cross section A cond k is the sum of the superconductor and copper areas. The non steel cross section (cable space) inside the steel conduit is A ns k = A cond k /0.7 taking into account that 30% of the cable space is allocated for voids. The fraction of steel is f steel = A steel k /A CS1 k, where 0 (1) (2)

A CS1 k includes the area of the electrical insulation except the ground insulation. The outer radius of sub coil k is: ACS 1 k ro k ri k hcs 1 2tg Here r i k is the inner radius of sub coil k, h CS1 the height of the CS1 module and t g = 8 mm [5] is the thickness of the ground insulation. It is supposed that each of the sub coils contributes B p /N to the peak magnetic field and the field at the CS outer radius is zero leading to B(r o k ) = B p (1 k/n) at the outer radius of sub coil k. The calculation of the radial and the hoop stress in a superconducting solenoid with different sub coil overall current densities follows the procedure described in [7]. (3) 7 2460.43 / 2498.27 2502.27 / 2540.11 8 2544.11 / 2588.41 2592.41 / 2636.71 9 2640.71 / 2678.93 2682.93 / 2721.15 10 2725.15 / 2762.27 2766.27 / 2803.38 Fig. 2 shows the overall current densities in the 10 sub coils. For CS1 1 and CS1 2 the use of HTS conductors is envisaged, while for CS1 3 to CS1 7 the use of R&W Nb 3 Sn ( 0.25% strain) is proposed. In the lowest field sub coils CS1 8 to CS1 10, the overall winding current densities are based on the use of NbTi. 5. Outline design of CS1 winding pack A sketch of the lower part of the CS1 winding pack (WP) is presented in Fig. 1. The WP design is based on 10 sub coils, each composed of two layers. The total number of conductors equals 20 78 = 1560. Each of the 10 sub coils carries a current of 8.36 MA leading to a conductor current of 53.6 ka. The proposed WP dimensions lead to B p as high as 18.15 T and a flux of 307 Vs exclusively generated by the 5 CS modules. Fig. 2. Overall sub coil winding current densities. Fig. 1. Sketch of the lower part of the CS1 WP. The WP including ground insulation is 5710 mm high. The inner and outer radii are 1892 and 2813.38 mm, respectively. Both, the conductor and the layer insulation are 1 mm thick. The radial positions of the conductors without insulation are listed in Table 1. Table 1. Inner and outer radii of conductors without insulation. Subcoil Inner layer r i / r o (mm) Outer layer r i / r o (mm) 1 1902.00 / 1947.88 1951.88 / 1997.76 2 2001.76 / 2046.82 2050.82 / 2095.87 3 2099.87 / 2144.49 2148.49 / 2193.10 4 2197.10 / 2238.52 2242.52 / 2283.94 5 2287.94 / 2327.60 2331.60 / 2371.26 6 2375.26 / 2413.85 2417.85 / 2456.43 Fig. 3. Variation of the total magnetic field (solid symbols) and the field generated only by the CS coil (open symbols) in the central plane of the CS1 module. Fig. 3 shows the numerically calculated values of the total magnetic field and the field generated only by the CS coil in the central plane of the CS1 module. The variation of the radial and the hoop stress in the central plane of the CS1 module is presented in Fig. 4. The peak value of the hoop stress is reached at the inner radius of the CS1 winding pack. From inner to outer radius, it decreases from 428.46 to 252.03 MPa. To limit the peak value of the hoop stress to 660 MPa in the steel the

fraction of steel needs to be as large as 64.92%. For the time being, the fraction of stainless steel is kept at the same value in all sub coils. The radial stress is always negative (compressive) and reaches a maximum absolute value of 16.6 MPa at a radius of about 2.44 m. 6. Proposed conductor design In the proposed conductor design, J Cu k and the area of voids have been adjusted according to geometrical needs. In Fig. 5, a sketch of the typical conductor layout is presented (L cond = 71 mm, r 4 = 3 mm). The resulting inner radius of the cable space is Finally, the length of the cable space is L cable = L cond 2t steel. For the two innermost sub coils (first 4 layers), the use of RE 123 (REBa 2 Cu 3 O 7 x with RE a rare earth element) HTS is envisaged. In Fig. 6, a sketch of the HTS cable for the innermost sub coil at the highest field is presented. Twelve strands of 7.5 mm diameter are wound around a copper core of 6.2 mm diameter. In each strand, 42 RE 123 tapes of 4 mm width and 95 m thickness are embedded. The HTS conductor properties are listed in Table 2. Table 2. Main characteristics of the 53.6 ka RE 123 cables for the innermost sub coils 1 & 2. Subcoil N str D str (mm) N tapes L cable d cable (mm 2 ) d core (mm) 1 12 7.5 504 46.32 21.20 6.20 2 13 7.0 455 46.52 20.58 6.58 Fig. 4. Variation of hoop and radial stress in the central plane of the CS1 module. Fig. 6. Sketch of the RE 123 cable for sub coil CS1 1. Fig. 5. Sketch of the typical conductor layout. The length L cond of the conductors without insulation is 71 mm in all 20 layers, whereas the conductor thickness d cond varies. r 3 Lcond d cond 2 Lcond dcond Lcond dcond 4 r 2 4 Asteel Here A steel is the cross section of the stainless steel conduit. The thickness of the cable space is d cable = 2r 3. The wall thickness of the steel conduit is t steel = d cond /2 r 3. (4) Fig. 7. Schematic layout of the Nb 3 Sn cable for sub coil 3. In the intermediate field sub coils 3 to 7, the use of R&W Nb 3 Sn conductors is planned. For the strain in the Nb 3 Sn filaments, a design value of 0.25% has been selected. In all considered Nb 3 Sn strands, the Cu : non Cu ratio equals 1. In Fig. 7, a sketch of the Nb 3 Sn flat cable for sub coil 3 is presented. Thirteen cable

stages, composed of 18 superconducting strands and one copper wire, form a flat cable. The residual resistivity ratio (RRR) equals 100 for copper in the strands and RRR = 400 for the copper in the solid mixed 95% Cu and 5% Cu Ni matrix. In order to reduce the ac loss the two layers of sub cables are separated by a thin stainless steel strip. The assembled cable is shown in Fig. 8. The cooling channels on the left and right hand sides are beneficial for welding of the two halves of the stainless steel conduit. The peak magnetic fields reach 14.48, 12.64, 10.81, 8.97 and 7.14 T in sub coils 3, 4, 5, 6 and 7, respectively. The design of the lowest field sub coils 8, 9 and 10 is based on the use of NbTi. The layout of the NbTi conductors is analogous to that of the Nb 3 Sn conductors. A design value of T cs 6.25 K has been selected for all NbTi and Nb 3 Sn sub coils. In the NbTi sub coils 8, 9 and 10, the peak magnetic fields reach 5.3, 3.47 and 1.63 T, respectively. The main properties of the Nb 3 Sn and NbTi conductors are listed in Table 3. 10 12 0.71 6 sc + 1Cu 12.5 3.75 7. Conclusions A winding pack design for a CS1 module composed of 10 sub coils with superconductor grading is presented. The use of HTS in the highest field sub coils allows us to maintain the reference magnetic flux in a CS coil of reduced outer diameter by means of an increased peak magnetic field. Next, a variation of the steel fraction in different sub coils is envisaged requiring a more detailed mechanical analysis including the vertical loads. Acknowledgments The technical support of the Paul Scherrer Institute (PSI) is gratefully acknowledged. This work has been carried out within the framework of the EUROfusion Consortium and has received funding from Euratom research and training programme 2014 2018 under grant agreement No 633053. The views and opinions expressed herein do not necessarily reflect those of the European Commission. References Fig. 8. Assembled Nb 3 Sn conductor for sub coil 3. Table 3. Main properties of 53.6 ka Nb 3 Sn and NbTi cables. Subcoil N sub D str (mm) Cabling pattern L sc d sc (mm 2 ) 3 13 1.21 18 sc + 1 Cu 37.8 10.38 4 16 0.87 18 sc + 1 Cu 33.44 7.50 5 20 1.10 6 sc + 1Cu 32.4 5.74 6 18 0.97 6 sc + 1Cu 25.6 5.11 7 18 0.82 6 sc + 1Cu 21.7 4.31 8 14 1.06 18 sc + 1 Cu 35.6 9.14 9 12 1.00 6 sc + 1Cu 17.6 5.21 [1] R. Wenninger, Reference Design 2015 April (EU DEMO1 2015) PROCESS Full Output. https://idm.euro fusion.org/?uid=2mdkfh [2] B. Meszaros, EU DEMO1 2015 DEMO Tokamak Complex, (CAD model) https://idm.euro fusion.org/?uid=2d3fbf [3] R. Ambrosino and R. Albanese, Reference flat top equilibria for DEMO with aspect ratio 3.1 https://idm.euro fusion.org/?uid=2aq5gp [4] R. Wesche et al., Winding pack proposal for the TF and CS coils of European DEMO, IEEE Trans. Appl. Supercond. 26 (3) (2016) 4200405. [5] K. Sedlak et al., Common operating values for DEMO magnets design for 2016. https://idm.euro fusion.org/?uid=2mmdtg [6] D. Uglietti et al., Test of 60 ka coated conductor cable prototypes for fusion magnets, Supercond. Sci. Technol. 28 (2015) 124005. [7] H. Burkhard, Magnetic and stress fields of cylindrical non uniform current density superconducting coils, Appl. Phys. 6 (1975) 357 362.