CERN-ACC-2014-0008 HiLumi LHC FP7 High Luminosity Large Hadron Collider Design Study Deliverable Report THERMAL-ELECTRICAL MODEL OF HILUMI SC LINKS Y. Yang SOTON et al 23 January 2014 The HiLumi LHC Design Study is included in the High Luminosity LHC project and is partly funded by the European Commission within the Framework Programme 7 Capacities Specific Programme, Grant Agreement 284404. This work is part of HiLumi LHC Work Package 6: Cold powering. The electronic version of this HiLumi LHC Publication is available via the HiLumi LHC web site <http://hilumilhc.web.cern.ch> or on the CERN Document Server at the following URL: <http://cds.cern.ch/search?p=CERN-ACC-2014-0008> CERN-ACC-2014-0008 THERMAL-ELECTRICAL MODEL OF HILUMI SC LINKS Doc. Identifier: HILUMILHC-WP6 Del-D6_2_final Date: 23/01/2014 Grant Agreement No: 284404 HILUMI LHC FP7 High Luminosity Large Hadron Collider Design Study S e ve n t h F r a m e w o r k P r o g r a m m e , C a p a c i t i e s S p e c i f i c P r o g r a m m e , R e s e a r c h I n f r a s t r u c t u r e s , C o l l a b o r a t i ve P r o j e c t , D e s i g n S t u d y DELIVERABLE REPORT T HERMAL -E LECTRICAL M ODEL OF H I L UMI SC L INKS D ELIVERABLE : D6.2 Document identifier: HILUMILHC-WP6 Del-D6_2_final Due date of deliverable: End of Month 24 (Oct 2013) Report release date: 23/01/2014 Work package: WP6: Cold Powering Lead beneficiary: University of Southampton Document status: Final Abstract: This document reports the results of the preliminary thermal-electrical model of the superconducting links (SC links) proposed for the cold-powering of superconducting magnetics at P1 and P5 for Hi-Luminosity Upgrades. The thermal-electrical model is concerned with the global characteristics of the superconducting link under both static and transient conditions. The model provides an outline assessment of the electromagnetic environment produced by the SC link and the heating dynamics upon different quench scenarios. Grant Agreement 284404 PUBLIC 1 / 16 THERMAL-ELECTRICAL MODEL OF HILUMI SC LINKS Doc. Identifier: HILUMILHC-WP6 Del-D6_2_final Date: 23/01/2014 Copyright notice: Copyright © HiLumi LHC Consortium, 2014 For more information on HiLumi LHC, its partners and contributors please see www.cern.ch/HiLumiLHC The HiLumi LHC Design Study is included in the High Luminosity LHC project and is partly funded by the European Commission within the Framework Programme 7 Capacities Specific Programme, Grant Agreement 284404. HiLumi LHC began in November 2011 and will run for 4 years. The information herein only reflects the views of its authors and not those of the European Commission and no warranty expressed or implied is made with regard to such information or its use. Delivery Slip Name Partner Date Authored by Y. Yang (SOTON), A. Ballarino (CERN) SOTON 24/11/2013 Edited by Ballarino CERN 26/11/2013 L. Rossi [Project coordinator] CERN 29/11/2013 A. Reviewed by A. Approved by Steering Committee Grant Agreement 284404 23/01/2014 PUBLIC 2 / 16 THERMAL-ELECTRICAL MODEL OF HILUMI SC LINKS Doc. Identifier: HILUMILHC-WP6 Del-D6_2_final Date: 23/01/2014 TABLE OF CONTENTS 1. INTRODUCTION ........................................................................................................................................ 4 2. BASELINE CONFIGURATION OF HI-LUMI SC LINK ....................................................................... 4 3. STATIC THERMAL ELECTRICAL CHARACTERISTICS ................................................................. 5 3.1. 3.2. 4. THERMAL CONSIDERATIONS AT STATIC CONDITIONS .............................................................................5 ELECTROMAGNETIC PROPERTIES AT STATIC CONDITIONS ......................................................................6 TRASIENT THERMAL ELECTRICAL CHARACTERISTICS ........................................................... 7 4.1. HEATING DYNAMICS OF A THERMALLY QUENCHED CABLE ....................................................................8 4.2. HEATING BY INDUCED CURRENT IN THE SC LINK UPON A FAST DISCHARGE OF AN INNER TRIPLETS CIRCUIT ..............................................................................................................................................................10 4.2.1. Field transients due to magnet discharging .................................................................................10 4.2.2. Induced current in the SC link cables by the discharging field transients ...................................11 5. REFERENCES ........................................................................................................................................... 15 ANNEX: GLOSSARY ..........................................................................ERROR! BOOKMARK NOT DEFINED. Grant Agreement 284404 PUBLIC 3 / 16 THERMAL-ELECTRICAL MODEL OF HILUMI SC LINKS Doc. Identifier: HILUMILHC-WP6 Del-D6_2_final Date: 23/01/2014 Executive summary This document reports the results of the preliminary thermal-electrical model of the superconducting links (SC links) proposed for the cold-powering of the superconducting magnets at P1 and P5 for Hi-Luminosity Upgrades. The model provides an outline assessment of the electromagnetic environment produced by the SC link and the heating dynamics upon different transient scenarios. 1. INTRODUCTION Thermal-electrical models for the superconducting links at P1/P5 for LHC High Luminosity Upgrade have been developed to evaluate the global of electromagnetic and thermal behaviours under static and transient conditions. At this stage of the project, the models do not focus on local analysis but on the overall characteristics. The main objectives of the models are to scrutinize the SC link baseline configuration against the required working envelope. The models also aim at identifying the critical aspects in the SC link design for further optimization, detailed local modelling and experimental studies/tests. 2. BASELINE CONFIGURATION OF HI-LUMI SC LINK The baseline configuration of the Hi-Lumi SC links powering the magnets in the insertion regions at P1 and P5 has been proposed [1] to fulfil the requirements of the updated baseline powering layout (Fig. 1) [2]. Fig. 1 Baseline powering layout for the Triplets and D1 at P1 and P5 [2]. Grant Agreement 284404 PUBLIC 4 / 16 THERMAL-ELECTRICAL MODEL OF HILUMI SC LINKS Doc. Identifier: HILUMILHC-WP6 Del-D6_2_final Date: 23/01/2014 The proposed SC link (Fig. 2) has a total current capacity of 150 kA, consisting of six 20 kA cables, seven 3 kA coaxial cables, four 0.4 kA cables and eighteen 0.12 kA cables [2]. The 20 kA cables uses hexagonal configuration of six sub-cables, which are constructed using 18 MgB2 round wires of 0.85 mm diameter twisted around copper stabilizer. The matrix of the MgB2 wire is predominantly Monel. The twist pitch is 100 mm - 300 mm. The same sub-cable is also used to form the inner leg of the 3kA co-axial cable. The outer return leg consists of 18 MgB2 wires and 6 copper wires of the same diameter. Cu MgB2 , = 0.85 mm 18 MgB2 wires = 6.5 mm 20 kA Six cables, = 19.5 mm Concentric 3 kA Seven cables, = 8.4 mm 0.4 kA Four cables 0.12 kA Eighteen cables ext 65 mm Fig. 2 Configuration of the P1/P5 SC link and its cable component [1]. The SC link contains two nominally 17.3 kA circuits for the quadrupoles (Q1+Q3 and Q2a+Q2b respectively) and one 11kA circuit for the dipole D1.In the SC link, these six cables will be identical and designed for transporting a maximum DC current of 20 kA at 25 K. 3. STATIC THERMAL ELECTRICAL CHARACTERISTICS 3.1. THERMAL CONSIDERATIONS AT STATIC CONDITIONS The thermal landscape at static conditions is a simple one. The helium gas cooled SC link operates between 5 K at the cold to the LTS bus-bar and maximum at about 20 K at the warm end near to the current leads. The present cable baseline gives a relative current loading of hence there is zero or negligible heat generation from within the SC link ( ) along the whole length. The forced flow helium gas should maintain a stable temperature gradient of ~0.15-0.2 K/m while removing the heat load from inner cryogenic envelope. Without a sudden interruption of the helium supply or failure of the vacuum envelope, such a static thermal condition is stable and does not interplay with the static electrical operation of the SC link. In theory, there is also a longitudinal conduction heat load along the link. However, the thermal conduction due to the ~ 10 cm2 copper stabilizer with the link only Grant Agreement 284404 PUBLIC 5 / 16 THERMAL-ELECTRICAL MODEL OF HILUMI SC LINKS Doc. Identifier: HILUMILHC-WP6 Del-D6_2_final Date: 23/01/2014 amounts to about 0.2 W at the 4.2 K end, negligible when compared to the helium gas flow for the link. 3.2. ELECTROMAGNETIC PROPERTIES AT STATIC CONDITIONS The main electromagnetic considerations are the magnetic fields produced by the current flow within the link and the inductances of the cables. The distribution of the magnetic field impacts on the critical current of the wires in the sub-cables. Since the low current cables (0.12-0.4 kA) are located inside the six 20kA-class cables, their exposure to the field produced by the high current system is minimal. Furthermore, the coaxial cables do not generate magnetic field to their exterior. Hence the focus is on the 20 kA system. 0T 0.8T Fig. 3 (a) Lines of force and (b) modulus of magnetic field of 6×20 kA cables arranged symmetrically for 3 circuits. The x-y scales are in mm. The six 20 kA cables in the SC link are configured with the cable-pair of one circuit placed next to each other to minimize their self-inductance. The magnetic field produced by the six 20 kA cables with the same ± current is similar to that of a sextupole magnet, as shown by the lines of force and field modulus in Fig. 3 (a) and Fig. 3 (b) respectively. In such a symmetric configuration, the maximum field occurs at the closest proximity between the opposing polarities, both intra-circuit and inter-circuit. For the mechanical stability of the cables, it is not desirable to break the symmetry polarity configuration. If all six cables are fully loaded to 20 kA, the maximum field is at , which is more than 50% above the ( ) maximum field of produced by a single 20 kA cable. With only one circuit powered, the magnetic field of the corresponding cable-pair at ±20 kA reaches a maximum of . Such a high field should be considered if the quadrupoles are not ramped up simultaneously. The calculations assumed that all the wires within a sub-cable carry the same current and no redistribution of current due to the change of magnet field along individual wires within the twist pitch takes place. It is further assumed that the sub-cables of a 20 kA cable also transport the same current. The field distribution within the SC link produced by the six 20kA cables becomes asymmetric (Fig. 4) when the magnets are powered to their nominal strength. In this case, the Grant Agreement 284404 PUBLIC 6 / 16 THERMAL-ELECTRICAL MODEL OF HILUMI SC LINKS Doc. Identifier: HILUMILHC-WP6 Del-D6_2_final Date: 23/01/2014 cables are configured into two pairs of ± 17.3 kA and one pair of ±11 kA. The maximum field, which now appears at the closest proximity between the ± wires of the same circuit, is 0.66T, almost (3.6% excess) proportionally to the reduction of current from 20 kA to 17.3 kA. The field between circuits are further reduced. 0T 0.7T Fig. 4 (a) Lines of force and (b) modulus of magnetic field of SC cables when powered at nominal currents: 2×(±17.3 kA) and 1×(±11kA). The x-y scales are in mm. The inductances of the cables are also obtained as for a single pair of 20 kA cables (i.e. two neighbouring 20 kA cables (+/-) in series) and for the six 20 kA cables arranged symmetrically for the three circuits. For these calculations, the twist pitch of the wires in the sub-cables is not considered as it has minimal impact on the field distribution of cables pairs. The inductance is higher with the current MgB2 wires due to the weakly ferromagnetic monel sheath. The focus on the cable-pair inductances is because the self-inductance of a single straight cable has a logarithmic singularity and is only become defined with a return current leg. The cable/link inductances are clearly very small and should not impose any restriction on the charging characteristics of the circuits. The mutual inductances between two neighbouring cable-pairs are of similar magnitude, again not important for field charging consideration. It should be noted that the wires within sub-cables are twisted, and the current mutually induced by the charging cross-talk between neighbouring circuits near the maximum field regions does not form loops longer than the twist pitch. 4. TRASIENT THERMAL ELECTRICAL CHARACTERISTICS The global dynamic models focus on SC link response to two likely transient events: (a) the quench of a cable and (b) the fast discharge upon quench of a magnet or a superconducting cable in the link. For the former, the main concern is the maximum temperature rise within the quenched cable and its capability to allow the safe discharge of the relevant magnets within the circuit. The latter is concerned with the transient magnetic field generated by the fast discharging circuit and the corresponding induced current heating in the neighbouring circuits. Grant Agreement 284404 PUBLIC 7 / 16 Doc. Identifier: HILUMILHC-WP6 Del-D6_2_final THERMAL-ELECTRICAL MODEL OF HILUMI SC LINKS 4.1. Date: 23/01/2014 HEATING DYNAMICS OF A THERMALLY QUENCHED CABLE The model only considers the adiabatic temperature rise within a cable after a local normal zone has exceeded the minimum propagation zone (MPZ) and an irreversible longitudinal propagation has started. The heat diffusion to the un-quenched section at the propagation front can be safely neglected when the quenched length becomes much longer than MPZ. The model makes no assumption of the causes for triggering a (localised) quench, although they are likely to originate from some local mechanical and/or thermal disturbances. The model further assumes sufficiently small thermal/electrical contact resistances between the MgB2 wires and the copper core to allow full current sharing with the stabilizer. It is noted that the sub-cable is thermally equivalent to the 20 kA cable, although the quench of a single subcable is unlikely not spread to the others. With the assumptions given above, the temperature rise upon a cable quench is given by: ̅ ̅ ( ) ( ) ( ( )) ( )∫ ∫ ∫ ( ) ( ) ( ) ∫ ( 100 Monel 10-1 Cu 10-2 10-8 cm-1 The specific heat of the cable consists of a portion of the superconductor wire and of the copper stabilizer. The superconductor wire is primarily monel and Ni, which has a similar heat capacity characterised by a high electronic heat capacity at low temperatures and a high Debye temperature of ΘD=460K for the phonon heat capacity. For the present model it is sufficient to use simple analytical expressions in J∙mol-1K-1: ( ) ( ) Cv, Jkg-1K-1 Where ( ( )) is known quench load and is determined as function of temperature by the thermal electrical property (density ̅ , heat capacity ̅ and stabilizer resistivity ρ) and the superconducting wire portion λ in the cable. The equation determines the cable temperature ( ) due to resistive (a) 102 heating at a current density of ( ) ( ) for a period t. Here the current density is referred to the total cross101 section area A of the cable. (b) 10-9 RRR 50 100 ) For and monel/Ni: for Cu: . For quench load calculation, cv is converted to J∙kg-1K-1 (Fig.5(a)). With a resistivity between at 10K-300K, monel’s contribution to current sharing is not considered. The copper resistivity [4] as function of temperature depends on the RRR (Fig. 5(b)) with a room temperature value of ~ . 10-10 300 1 10 100 Temperature, K Fig. 5, (a) The relevant specific heat for the SC link cable and (b) resistivity of copper at different RRR as functions of temperature. In case of resistive transition of a cable in the link, a fast discharge of the corresponding magnet chain shall be triggered. While final values for the time constants of the circuits are not yet available, indicative values in the range from 1 s to 3 s were given by the magnet Grant Agreement 284404 PUBLIC 8 / 16 THERMAL-ELECTRICAL MODEL OF HILUMI SC LINKS Doc. Identifier: HILUMILHC-WP6 Del-D6_2_final Date: 23/01/2014 designers. It is noted that the sub-cable and 20 kA main cable have the same quench load due to their identical superconductor fill factors and the nominal current density. The present subcable design has a with a superconducting wires area of 0.102 cm2 with 18 wires of 0.85 mm diameter and a copper cross-section of 12.5 cm2 (about 50% of the available centre core area). The sub-cable has a current density of for the nominal 17.3 kA of the 20 kA cable. The results (Fig. 6) show that at a constant nominal current (2833/sub-cable or ( ) ( )) and room 17.3 kA/cable), the cables will reach 75 K in 1 second ( temperature in 3.5 s. A higher RRR=300 (blue lines) reduces the temperature rise to 50 K and 225 respectively. Although the time for heating is slightly longer for starting at T0=5 K, the assumption of a fully established quench region means the onset temperature is more likely at 20 K or above. 109 (12.5kAcm-2)2(3.5s) (12.5kAcm-2)2(1s) 300K 50K U(T), A2cm-4s 108 75K RRR=100 RRR=300 10 225K 7 T0=20K 10 6 T0=5K 10 100 Temperature, K Fig. 6 Quench load of the baseline SC link cable. Let’s assume that the quench is detected at its onset and the magnet is subsequently discharged at a time constant τ with ( ) . In this case the quench load becomes: ( ( )) ∫ ( ) ( ) The effective constant current heating time for an exponentially discharging current is ( ) which is shorter than the actual time t and terminal temperature at is arrived at . For example, s is equivalent to the full course of discharging at . In this case the cable temperature will not exceed room temperature for time constants of the circuit in the range from 1 to 3 s. The practical indicator for temperature rise is at or , e.g. at for the temperature rises to 220 K with RRR=100 and 150K with RRR=300. For of a reduced time constant of τ = 1s, then cable temperature rise is significantly lower. At the corresponding effective heating time of is only 0. 43s after 1s discharge at τ=1s and Grant Agreement 284404 PUBLIC 9 / 16 Doc. Identifier: HILUMILHC-WP6 Del-D6_2_final THERMAL-ELECTRICAL MODEL OF HILUMI SC LINKS Date: 23/01/2014 the maximum temperature rise will be at 46 K and 30 K for RRR=100 and RRR=300 respectively. Fig. 7 The distribution of field change in the SC link due to the discharging of a quadrupole circuit. 10 105 -1 (a) (b) r/RCable 75 120 60 135 45 150 10-2 dB/dt, Ts-1 90 30 165 10-3 15 180 0 195 10-2 1.5 1.4 1.3 1.2 1.1 1.0 0.9 0.8 0.7 0.6 0.5 345 210 330 225 315 240 10-1 300 255 270 285 Fig. 8 The detailed field change imposed on the wires from the neighbouring cables immediate adjacent to the fast discharging cables. 4.2. HEATING BY INDUCED CURRENT IN THE SC LINK UPON A FAST DISCHARGE OF AN INNER TRIPLETS CIRCUIT 4.2.1. Field transients due to magnet discharging The largest field transient within the SC link occurs upon the fast discharge of quadrupole circuit at 17.3 kA. The present model examines the field change rate ̇ ( ) in the link by using ( ) ( ( )) with ̇ ( ) the field produced by the discharging circuit ( ) is the field produced by the 20 kA cable-pair for the circuit. As shown in where Fig. 7, the largest field change is within the circuit cable-pair itself with a peak value of ( ) 0.77 T. However, a significant field of 0.3 T is imposed on the immediate adjacent ~3 wires in the neighbouring circuits (at 30o and 150o in Fig. 7), as shown Grant Agreement 284404 PUBLIC 10 / 16 Doc. Identifier: HILUMILHC-WP6 Del-D6_2_final THERMAL-ELECTRICAL MODEL OF HILUMI SC LINKS Date: 23/01/2014 ) is the Fig. 8(a). For a discharging time constant of τ=7 s, the initial ̇ ( )= ( o o shown in Fig. 8(b), where 40 mT/s is found near the 30 and 150 regions of the neighbouring cables. The shaded area in Fig. 8(b) represents the region occupied by the cables of the discharging circuit. Note that the field change rate will increase significantly to 0.3 T/s if the quadrupole discharge time constant is decreased to 1s. 4.2.2. Induced current in the SC link cables by the discharging field transients 4.2.2.1. Magnetisation energy of the critical state upon field discharge The superconducting wires in the SC link reacts to the field transients by the induction of “super” current to resist the field change. As the result, work is done by the discharging field to the neighbouring circuits by increasing their Gibbs’ energy , for a discharging field and a critical state magnetisation to retain the flux. The overall energy deposited is proportional volume enclosed by the induced current loop, i.e. . Assuming the super-current penetrates only a portion of the superconducting strands, then the magnetisation is approximately , where is the width of the induced current loop, is the longitudinal length of the induced current loop, and Ic is the critical current of the superconducting wire. The overall magnetization energy of the critical state upon field discharge is When the stored energy is released due to flux diffusion, e.g. flux creep of the critical state or the power-law E-J characteristics of most superconductors, the magnetization energy is converted irreversibly to heat with an effective power of When the flux diffusion time constant becomes comparable to or shorter than the field discharge constant , as often the case for high temperature superconductors (HTS), the heat generation is almost simultaneous: ̇ ( ) 4.2.2.2. Size of the magnetization loops STRAIGHT UNTWISTED FILAMENTS For a superconducting wire with fully coupled filaments, the coupling width is approximately the wire diameter and the induction loop is extended to the whole length of the conductor, hence the energy stored per unit volume (of wire) is If the filaments are completely uncoupled, e.g. by a high matrix resistance, then the coupling width is reduced to the filament diameter . The corresponding energy (per unit volume of wire) by induction for uncoupled filaments is Grant Agreement 284404 PUBLIC 11 / 16 Doc. Identifier: HILUMILHC-WP6 Del-D6_2_final THERMAL-ELECTRICAL MODEL OF HILUMI SC LINKS Date: 23/01/2014 Since most low temperature superconductors (LTS) have very fine filaments, , the induced energy in a filaments-uncoupled wire is much lower, even though the critical current penetration in the fine filaments is higher ( ). Clearly filaments uncoupling is essential to reducing the induction energy upon field change. TWISTED FILAMENTS In practice, however, high matrix resistance for filaments uncoupling is not practical for the cryogenic stabilization, which requires a low normal resistivity for current sharing upon magnet quench. Filaments uncoupling in technical superconductors are achieved by twisting the filaments in a wire at a short pitch so that the coupling loop across the wire diameter is broken as the twisted filaments are displaced from one side of the wire to the other. It should be noted that finer filaments for lower induction energy requires a shorter twist pitch: √ ̇ where is the critical coupling length between filaments separated by a resistive ( ) matrix. The matrix resistance that breaks coupling loop also limits the effective induced ̇ current across the wire which is significantly less than the typical critical current density . However, this coupling current among the twisted filaments leads to a resistive (eddy-current) type of loss known as the coupling current loss ̇ ( ) A successful filament architecture should always have . INDUCTION LOOPS AND CORRESPONDING HEATING IN THE SC LINK CABLES A. Estimated temperature rise due to adiabatic induction heating of different loop sizes We first estimate the upper limit of temperature rise due to adiabatic induction heating by different loop sizes. Although the terminal adiabatic temperature rise is independent of the discharge rate, it is assumed that the flux diffusion is faster than the discharge rate for an instantaneous heat rate per unit length , where the electrical field ̇ and of complete flux penetration in the wire with full critical current ( ) induced. Consider the reference nominal case of a MgB2 round wire of resistivity and a wire critical current ( filaments of diameter with a filament ( ) . There are three distinct coupling widths to consider: 1. Fully uncoupled filaments: 2. Fully coupled filaments: 3. Inter-wire coupling across the Cu stabilizer in the sub-cable: . The thermal mass includes the contributions from the Grant Agreement 284404 PUBLIC ) , matrix (monel) . The wire has 30 12 / 16 Fig. 9 An example of MgB2 round wire used for reference calculation. THERMAL-ELECTRICAL MODEL OF HILUMI SC LINKS Doc. Identifier: HILUMILHC-WP6 Del-D6_2_final Date: 23/01/2014 monel matrix and the copper stabilizer. This large coupling loop requires a negligible resistances in the wire matrix and at the contact between the wire and stabilizer. Similar to the quench load calculation, the temperature rise by induction heating for an exponential magnet discharge of time constant is given by ̅ ( ) ( ) ( ) where is the cross-section heated by the induction current with an instantaneous heat diffusion, and ̅ and ̅ are the weighted average of the specific heat and density respective to the composition of . Thermal diffusion to the copper stabilizer is ignored for case 1-2, i.e. , consistent with a sufficiently high contact resistance to prevent the large coupling loop across the sub-cable. For the same reason, a portion of the stabilizer thermal mass should be included in case 3. 25 20 Temperature, K T0=20K 15 Uncoupled filaments Coupled filaments without Cu Inter-wire coupling with Cu 10 T0=5K 5 10-2 10-1 100 t/ Fig. 10,Adiabatic temperature rise in a MgB2 round wire by inductive current heating due to exposure to the field change by the adjacent cable of a discharging magnet. Three different loops are considered. The differential temperature rise equation is integrated to give an explicit expression of elapsed time to reach temperature : ̅ ̅ ( ) ( ∫ ) ( ) ( ) ( ) Fig. 10 shows the time for temperature rise at for the three cases assuming a wire critical current ( ) linearly reducing with increasing temperature from 590A at 5K to zero at 37K. As expected, fully uncoupled filaments heat up the least, with a maximum temperature rise <2K at the cold end of 5K. Although wire coupling across a sub-cable through the stabilizer has the highest heating power, the temperature rise was tempered by the thermal mass of the copper. An assumed 30% increase of copper limits the cold end temperature rise to 4.5K. In fact, for temperature rise for case 1 and 3 reduces rapidly with increasing temperature, the upper temperature limit of 20K is not exceeded. In contrast, the Grant Agreement 284404 PUBLIC 13 / 16 THERMAL-ELECTRICAL MODEL OF HILUMI SC LINKS Doc. Identifier: HILUMILHC-WP6 Del-D6_2_final Date: 23/01/2014 full coupling of filaments leads to substantial inductive heating within the limited thermal mass of the wire. At cold end (5K), the temperature is increased to above 17.5K, and the warm end to almost 23K. Actually, fully coupled wires initially at (covering about 1/3 of the total Hi-Lumi SC link length at the warm side) will be heated above upper temperature limit of 20K, hence likely to quench or thermally runaway. B. Filament uncoupling by twisting The MgB2 wires in the SC link sub-cables are twisted around the copper stabilizer core at a pitch ~400mm. Even without filament twisting, the effective induction loop length with a wire immediate adjacent to the discharging circuit is , where is the number of wires in close vicinity of the discharging cables and the number of wire in the sub-cable. The coupling length is reduced to the filaments twist pitch if the wire it is shorter. For ̇ at a discharging time constant of 1s, the reference wire of the critical ( ) coupling length at 20K with ( ) . The coupling length at lower temperatures will increase further as increases. Therefore, even if untwisted, the filaments of the reference wire can only be partially-coupled. The coupling current heating through the resistive matrix of partially coupled untwisted filaments is ̇ ( ) ( ) ̇ Which is much less than of fully coupled filaments, hence the conductor will have a negligible temperature rise by the coupling current across the wire matrix. It should be noted that the critical current density for a coupled length of is only ( ) for the reference wire and much lower than the nominal transport current loading of ( ). Hence the contribution by the filament coupling current to E-J heating can be safely ignored within the working temperature range for the SC link cables. Further reduction of inductive heating can be achieved by twisting the filaments within the wire at a pitch less than 90mm, which is the target for the MgB2 round wire. 5. CONCLUSIONS AND FUTURE WORK Static and transient models for the HiLumi SC link cables were established for their outline assessment against the baseline operation conditions. The main conclusions from the modelling analysis are: 1. A maximum field of 0.8-1.0T, corresponds to the powering up of a 20kA single circuit, should be allowed in the wires specification. 2. The self-inductances of the cables are small compared to the magnets in the circuit. Grant Agreement 284404 PUBLIC 14 / 16 THERMAL-ELECTRICAL MODEL OF HILUMI SC LINKS Doc. Identifier: HILUMILHC-WP6 Del-D6_2_final Date: 23/01/2014 3. With a sufficiently good electrical contact to the copper stabilizer, the wires in the subcables undergo moderate temperature rise upon cable thermal quench. For a reference RRR=100 for the stabilizer, the temperature rise is less than 80K if the magnets are discharged at 1s upon quench detection. Slower discharge at 7s leads to a temperature rise to 220K. 4. Upon fast discharge of a quadrupole circuit following a magnet quench (discharging time constant of 7 s), a field change of 0.3T is imposed on the immediately adjacent wires for the neighbouring cables. 5. Fully coupled filaments in these adjacent wires would suffer significant inductive heating leading to potential quench or thermal runaway in cable section at temperatures above 14-15K. 6. With a matrix resistance of and filaments of diameter , the critical filament coupling length at ̇ is much longer than the exposed wire section of 90mm because of the sub-cable twist pitch of 400mm. Therefore the filaments should only be partially coupled. 7. The coupling current loss for a matrix resistance of and filaments -1 of diameter is about 0.001Wm at ̇ , insufficient to cause wire to quench or even to heat up significantly. 8. Filament twist pitch <90mm within the wire will help to further reduce inductive heating and uncouple the filaments. 9. Matrix resistivity between the filaments should not be reduced below at . In-wire stabilizer should be place on the outside of the filaments. Further to the above conclusions, following future studies are recommended: 1. Determine the minimum wire to stabilizer electrical and thermal contact resistances necessary for full current sharing and full transverse thermal diffusion upon thermal quench. 2. Study the minimum quench energy and quench propagation velocity for the sub-cable and cables to establish the criterion for quench detection and the identification of possible quench triggers. Work in this direction has already been done for the SC-Link proposed for LHC P7 [5]. 3. Experimental study and modelling of filaments coupling in MgB2 round wires. 4. Optimize filament architecture and wire configuration to ensure full uncoupling upon imposed field changes from magnet discharge in neighbouring circuits. 5. Study of the e.m. coupling between the different cables to be sure that there is no spurious quench signal due to fast dump of a circuit 6. REFERENCES [1] A. Ballarino, Development of Superconducting Links for the LHC Machine, EEE/CSC & ESAS SUPERCONDUCTIVITY NEWS FORUM (global edition), October 2013 [2] A. Ballarino, Presentation at 4th LHC Parameters and Layout Committee, 26th March 2013, https://indico.cern.ch/conferenceDisplay.py?confId=239311 [3] NIST Cryogenic Materials database: www.cryogenics.nist.gov/MPropsMAY/material%20properties.htm Grant Agreement 284404 PUBLIC 15 / 16 THERMAL-ELECTRICAL MODEL OF HILUMI SC LINKS Doc. Identifier: HILUMILHC-WP6 Del-D6_2_final Date: 23/01/2014 [4] R. L. Powell and F. R. Fickett, Cryogenic Properties of Copper, International Copper Research Association, Dec. 1979. [5] J Spurrell, W Bailey, E A Young, A Ballarino, Y Yang, Cryogenic Stability and Quench Characteristics of Twisted-pair Cables Assemblies of MgB2 tapes, Presented at EUCAS 2013, to be published Grant Agreement 284404 PUBLIC 16 / 16
© Copyright 2024 ExpyDoc