MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING DEPARTMENT OF NUCLEAR ENGINEERING 2.64J/22.68J , : & HTS

Transcription

1 MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING DEPARTMENT OF NUCLEAR ENGINEERING 2.64J/22.68J Spring Term 23 May 8, 23 Lecture 1: Protection & HTS Magnets Key magnet issues vs. T op ; Areas of failure in superconducting magnets Key problems: overheating; high voltage; overstressing; over pressure Protection against overheating: self-protecting magnets; normal zone propagation; active protection; detection of quench High voltage driven mode & persistent mode; overstressing Over pressure bath cooled & CICC HTS Magnets & Issues Bi-2223/Ag Data YBCO composite J c ( B 4.2, 1, and 2 K of superconductors 1-D NZP of YBCO: Energy margin & NZP velocity (FSU/NHMFL) Stability experiment & simulation of YBCO at 77 K (MIT) 1-D NZP of YBCO: Energy margin & NZP velocity (ORNL) Y. Iwasa (5/8/3) 1

2 Key Magnet Issues vs Top Cost or Difficulty Protection Conductor Mechanical Stability Cryogenics ~1 T op [K] LTS HTS Y. Iwasa (5/8/3) 2

3 Areas of Failure in Superconducting Magnets In the order of decreasing reported incidents. Insulation Mechanical System performance Conductor External system Coolant Y. Iwasa (5/8/3) 3

4 Key Quench-Induced Problems Overheating High voltage Overstressing Over pressure Y. Iwasa (5/8/3) 4

5 Overheating Energy Conversion B o 2 2 o h cu (T mt )-h cu (T op ) = J/m 3 B o 115 T 1356 K 4.2 K For a uniform energy conversion over the entire winding, a quench-induced T max will remains <1 K, an upper limit to keep the differential thermal expansion stresses negligible. Y. Iwasa (5/8/3) 5

6 Protection Against Overheating 1. Self-Protecting Magnets Normal zone propagation speed, U nzp, fast enough to spread normal zone over most of the winding volume from a localized hot spot initially created. Specifically, a selfprotecting magnet must satisfy: decay ~ a 2 a 1 U nzp 2a 1 2a 2 B o Large magnets are not self-protecting. [U nzp ] HTS << [U nzp ] LTS : HTS magnets not self-protecting. Y. Iwasa (5/8/3) 6

7 An Example B o =8 T B o E mg =1 6 J.6 m Winding volume=.6 m 3 (~2 ft 3 ) Winding % T f [K] Remarks.2 m Absorbing E mg.4 m 1 55 Well below 1 K 5 7 < 1 K 1 13 Barely acceptable 1 58 Unacceptable Y. Iwasa (5/8/3) 7

8 Normal Zone Propagation Velocity Along Wire Axis--No Matrix C n T n t z k n T n z n J 2 (normal state) C s T s t z k s T s z (Superconducting state) T t T z z t U l T z (U l : longitudinal velocity) Normal Superconducting z Y. Iwasa (5/8/3) 8

9 C n U l dt n d z C s U l dt s dz d d z k n d dz k s dt n d z dt s dz Assume k n, k s, and d 2 T n /dz 2 n J 2 near z C n U l dt n dt nj 2 (z ) k s d 2 T s dz 2 C s U l dt s d t (z ) T s (z) = Ae cz T op C C s U l k s Y. Iwasa (5/8/3) 9

10 At z T s () = T c T s (z) = (T c T op ) e ( C s U l / k s )z T op At z k n (dt n /dz) = k s (dt s /dz) k n n J 2 C n U l C s U l (T c T op ) U l J n k n C n C s (T c T op ) Y. Iwasa (5/8/3) 1

11 C s C n C U l J C n k n (T c T op ) U l J U l 1/C [U l ] hts << [U l ] lts In the presence of normal metal matrix: U l C cd J m m k m (T cs T op ) Note that T cs T t Y. Iwasa (5/8/3) 11

12 With C n ( T ), C s ( T ), k n ( T ): U l J C n (T cs ) - n (T cs )k n (T cs ) 1 d k n (T ) T cs C s (T )dt k n (T cs ) dt T op Tcs T cs T op C s (T )d T Note that T cs T t Y. Iwasa (5/8/3) 12

13 1-D Adiabatic NZP Experimental & Simulation Results ~1 cm Slides 13-15: from R.H. Bellis and Y. Iwasa (Cryogenics, Vol ) Y. Iwasa (5/8/3) 13

14 Nb 3 Sn 12 K & 225 A Experimental U l =5 cm/s V 4 V 1 V 2 V 3 Temperature Simulation Simulation 15 ms V 4 25 ms 1 ms T c =16.5 K T cs =13.7 K V 1 V 2 V 3 Y. Iwasa (5/8/3) 14

15 Bi-2223/Ag 14 K & 1 A Experimental [U l ] hts =1 cm/s << [U l ] lts Temperature Simulation Simulation 6 s 1 s T c =93 K T cs =44.4 K Y. Iwasa (5/8/3) 15

16 2. Active Protection Detect-and-Dump Widely used in large magnets. Dissipates most of the stored magnet energy into a dump resistor connected across the magnet terminals. A hot spot heated only over a brief period of current decay. Leads to another criterion for operating current density. L mg ~ V D R D r ( t ) Y. Iwasa (5/8/3) 16

17 L mg di(t ) d t For r(t ) << R D : L mg di(t ) d t [r(t ) + R D ]I(t ) = R D I(t ) = With g k g d g q : A cd C cd (T ) dt dt C cd (T ) A dt m (T ) m A cd I(t ) = I op exp (T) m I 2 (t) C cd (T) dt A m d t [J op ] 2 m exp 2R D t L mg dt R D t L mg m (T ) A m A cd J m (t ) = [J op ] m exp J 2 m (t ) R D t L mg T f C cd (T ) dt Z (T i, T f ) m (T ) T i A m A cd 2 [J op ] m exp 2R d t L mg dt A m A cd 2 L [J op ] mg m 2R D Z (T i, T f ) A m A cd 2 L [J op ] mg m 2R D A m [J op ] 2 E I 2 2 mg op A m m A cd 2V D I op A cd [J op ] m 2 E mg V D I op Y. Iwasa (5/8/3) 17

18 Z (T i, T f ) T f T i C cd (T ) m (T ) dt C cd (T ) C (T ) (T ) cd m m (T ) ~constant Z(T i, T f ) ~T ~1/T T T T i T f T Y. Iwasa (5/8/3) 18

19 Z(T ) Functions: 1: Ag (99.99%) 2: Cu (RRR 2) 3: Cu (RRR 1) 4: Cu (RRR 5) 5: Al (99.99%) Z [1 16 A 2 s/ m 4 ] T [K] Y. Iwasa (5/8/3) 19

20 Protection Criterion [J op ] m A c d A m V D I op Z (T i, T f ) E mg Cryostable Criterion [J op ] cd f p P cd A m q fm m ( A s + A m ) 2 Y. Iwasa (5/8/3) 21

21 Detect-and Dump : Options to improve [J op ] m 1. Increase I op Larger conductor; more expensive for given ka-m Larger current leads--greater heat input Larger ( I B ) force For a given power rating VI, higher I more expensive 2. Increase V D Increased danger of voltage breakdown, particularly > 8 V in helium environment 3. Increase Z( T f ) Better to keep T f below ~1 K-15 K Cu better than Al 4. Reduce E mg Subdivision common practice with HEP & Fusion magnets Y. Iwasa (5/8/3) 22

22 Detection of Non-Recovering Quench di(t ) di(t ) V 1 (t ) L 1 r(t )I(t ) V dt 2 (t ) L 2 dt V out ( t ) = V 2 ( t ) - V 1 ( t ) R 2 L 1 R 1 L 2 V out (t ) = - R 2 R 1 R 2 r(t )I(t ) I( t ) ~ V 1 ( t ) L 1 r( t ) I( t ) o V out ( t ) o R 1 R 2 V 2 ( t ) L 2 Y. Iwasa (5/8/3) 23

23 High Voltage Driven mode magnet, i.e., operated with its power supply generally connected, except during a dump. High voltage, i.e., V D, appears across the magnet terminals. Persistent mode magnet, i.e., operated with its power supply removed, except when the magnet is being charged. Terminal voltage is zero; high internal voltage may be induced. Y. Iwasa (5/8/3) 24

24 Operation Steps Persistent Mode Operation 1. With persistent switch (PS) resistive (heater on), energize the magnet. 2. Turn off the heater; PS becomes superconducting; slowly bring the supply current back to zero. 3. Disconnect the current leads from the magnet terminals. 4. Magnet now operates in persistent mode. Disconnectable joint > > Persistent Switch L mg ~ Heater r ( t ) > Y. Iwasa (5/8/3) 25

25 Persistent Mode Operation (Continued) r ( t ) represents: 1. quench-induced variable resistance; may lead to a high internal voltage OR 2. small constant resistance by imperfect splices; causes a slow field decay (time constant=l mg / r ) in NMR & MRI magnets: B(t ) B e (r / L mg )t B 1 r L mg t (for L mg /r 1hr) Y. Iwasa (5/8/3) 26

26 Internal Voltage In Persistent Mode Magnet No internal voltage in a uniformly (1%) resistive magnet. Internal voltage in a partially resistive magnet. Large voltages can occur when a resistive zone confined to a small section. V L 1% V R ~ r ( t ) 1% No internal voltage Y. Iwasa (5/8/3) 27

27 Internal Voltage in Persistent Mode ~ Extent of r ( t ) 5% 1% 1% 5% 1% Top 1% Middle 1% Top 5% Bottom 5% Bottom 1% Y. Iwasa (5/8/3) 28

28 Internal Voltage In Driven Mode Magnet For R D >> r ( t ), internal voltage is essentially due to the inductive voltage. The maximum inductive occurs at t = and it is equal V D. V L V D (R D >> r ) L mg ~ V D R D r ( t ) Effect of r ( t ) is neglected. Y. Iwasa (5/8/3) 29

29 Overstressing Overstressing generally occurs in a multi-coil system. A quench in one coil increases the current in the rest of the coils this is good for PROTECTION, because this increased current accelerates spreading of normal zones in the rest of the coils by inducing quenches in these coils. This increased current may also cause overstressing in these coils. Y. Iwasa (5/8/3) 3

30 Example: 2-Coil System Overstressing (& Protection) A quench in Coil 1 induces an extra current in Coil Helps spreading out normal zone faster and over a large volume. 2. May overstress Coil 2 I ( t ) Overstress limit I 1 I 2 I o I 1 I o R 1 L 1 I 1 r (t ) ~ M I 1 I 2 R 2 I 2 L 2 I o t I o I 1 I 2 Y. Iwasa (5/8/3) 31

31 Protection Circuit for an NMR Magnet + Y. Iwasa (5/8/3) Within Cryostat 32

32 Over Pressure In every superconducting magnet system, driven, persistent, in which coolant is cryogen (helium, nitrogen), a quench generally causes generation of vapor within the cryostat, raising the cryostat pressure above a safe limit (no greater than ~ 1/3 atm). Cryostat must be equipped with a relief valve or a burst disk to permits a rapid release of the vapor from the cryostat. For CICC there is another problem of having to estimate the internal pressure generated in CICC. Conduit thickness and vent line diameter must be sized according to the maximum quench-induced internal pressure. Y. Iwasa (5/8/3) 33

33 HTS Magnets & Issues Hc2 [ T ] 15 YBCO HTS: YBCO; BSCCO; MgB 2 LTS: Nb-Ti; Nb 3 Sn 1 5 Nb 3 Sn BSCCO 2223 BSCCO 2212 Magnet Grade Conductors: NbTi; Nb 3 Sn; BSCCO 2223 NbTi MgB T [K] Y. Iwasa (5/8/3) 34

34 HTS BSCCO 2223 currently available as magnet grade conductor. Powder in Ag tube BSCCO 2212 being developed but not easily available. YBCO, compared with BSCCO, regarded more promising for electric power devices. MgB2 also regarded promising, because of its low cost. Y. Iwasa (5/8/3) 35

35 Bi-2223 High Current Density Wire: Non-Reinforced Thickness (avg):.21 (+/-.2mm) Width (avg): 4.1 (+/-.2mm) Min. Critical Stress: 75 MPa Min. Critical Strain:.15% Min. Bend Dia: 1 mm Variable Specifications: Min. I c : 115 A 135 A Piece Length: 1 m 2 m Ic/Ic(o) Strain Ic/Ic(o) Strain Stress at 77K (MPa) Slides 35-39: Based on American Superconductor Corp. Data Sheets Y. Iwasa (5/8/3) 36

36 Bi-2223 High Strength Reinforced Tape Stainless Steel Strips Thickness (avg):.31 (+/-.2mm) Width (avg): 4.1 (+/-.2mm) Min. Critical Stress: 265 MPa Min. Critical Strain:.4% Min. Bend Dia: 7 mm Variable Specifications: Min. I c : 115 A 135 A Piece Length: 1 m 3 m Ic/Ic(o) Strain Ic/Ic(o) Stress at 77K (MPa) Strain (%) Y. Iwasa (5/8/3) 37

37 Scaling Ratio to 77K (self-field) in // External Field 6 Scaling Ratio, Ic(T,B)/Ic(77K,) Parallel Magnetic Field (Tesla) Y. Iwasa (5/8/3) 38

38 Scaling Ratio to 77K (self-field) in External Field 6 Scaling Ratio, Ic(T,B)/Ic(77K,) K 35K 5K Perpendicular Magnetic Field (Tesla) Y. Iwasa (5/8/3) 39

39 Scaling Ratio, 77K to 4.2K Bpar, 4.2K Bperp, 4.2K Magnetic Field, Tesla Y. Iwasa (5/8/3) 4

40 1-D NZP in YBCO Composite Energy Margin & NZP Velocity Keithley 21 DMM* Keithley 21 DMM* Keithley 21 DMM* Keithley 2 DMM Keithley 2 DMM Keithley 2 DMM V 3 V 2 V 1 V cernox1 V cernox2 V cernox3 HTS sample NiCr heater HP 6286A DC-P/S Keithley 2 DMM* V tot HP 8112A pulse generator 5 A / 1 mv V(I) HP 6681A DC-P/S Keithley 2 DMM Keithley 21 DMM* GPIB *set on buffer read for increased speed/resolution Slides 45-51: Based on Justin Schwartz ( DOE Meeting July 22) Y. Iwasa (5/8/3) 45

41 Bi2223 tape: Recovery.8 Voltage (V) V 1 V 2. V Time (s) 11 Temperature (K) 1 9 C 2 C Time (s) Y. Iwasa (5/8/3) 46

42 Bi-2223 High Current Density Wire: Non-Reinforced Thickness (avg):.21 (+/-.2mm) Width (avg): 4.1 (+/-.2mm) Min. Critical Stress: 75 MPa Min. Critical Strain:.15% Min. Bend Dia: 1 mm Variable Specifications: Min. I c : 115 A 135 A Piece Length: 1 m 2 m Ic/Ic(o) Strain Ic/Ic(o) Strain Stress at 77K (MPa) Slides 36-4: Based on American Superconductor Corp. Data Sheets Y. Iwasa (5/8/3) 36

43 I op =19 A; I c =22 A 81.1 Temperature (K) C 1 (closest) C 2 C 3 (closest) Time (s) Y. Iwasa (5/8/3) 48

44 19 A & 8.6 K 3.E-2 2.5E-2 Voltage (V) 2.E-2 1.5E-2 1.E-2 5.E-3 Primary Heater pulse V 1 V 2 V 3 V 4 1.E Time (s) Y. Iwasa (5/8/3) 49

45 Energy Margin vs. I/I c (I c = K ) Minimum Quench Energy (mj) Fraction I/I c (%) Y. Iwasa (5/8/3) 5

46 NZP Velocity vs. I/I c Quench Propagation Velocity (cm/s) I/Ic (%) Y. Iwasa (5/8/3) 51

47 Quench/Recovery Experiment & Simulation Bath Forced-flow 2-3 atm 8-81K Y. Iwasa (5/8/3) 52

48 Forced Flow Experimental Setup Pressure Gauge Pressure Regulator Needle Valve Shut/OpenVolume Valve Meter Liquid Nitrogen Water N 2 Gas Tank 77K Hx Sample Holder 77K (Styrofoam Bath) RT Hx High-Pressure (3-atm) Section Y. Iwasa (5/8/3) 53

49 Y. Iwasa (5/8/3) 54

50 New Forced-Flow Sample Holder Y. Iwasa (5/8/3) 55

51 Y. Iwasa (5/8/3) 56

52 Y. Iwasa (5/8/3) 57

53 Y. Iwasa (5/8/3) 58

54 Y. Iwasa (5/8/3) 59

55 Sample 2: V( I ) [ I c =15A ] [ Source-Measured I c > 1A ] 2 15 V [uv] I [A] Y. Iwasa (5/8/3) 6

56 Sample 2: Forced-Flow (~81K; 5cm/s) A I [A] A 45A 65A 71A A V [mv] A 45A 24A Time [s] Y. Iwasa (5/8/3) 61

57 Sample 2: Forced-Flow (~81K; 3cm/s) 2 15 I [A] A 45A 65A 1 8 V [mv] A 2 45A 24A Time [s] Y. Iwasa (5/8/3) 62

58 Sample 2: Forced-Flow A/65A (3cm/s) V [mv] A/65A (5cm/s) Time [s] Y. Iwasa (5/8/3) 63

59 Sample 2: Forced-Flow (5cm/s) & Bath A/71A (5cm/s) V [mv] A/75A (Bath) Time [s] Y. Iwasa (5/8/3) 64

60 Simulation Re 3, turbulent flow Q = h T w =ha s (T w -T f ) h=knu/d Nu =.23Re.8 Pr Re < 3, laminar flow.668 Nu = ( D / L) Re Pr [( D / L) Re Pr] 2/ 3 T(K) Nitrogen Properties at 2 bar ( From Handbook of Cryogenic Engineering, I. G. WeisendΙΙ, 1998) ρ(kg/m 3 ) Cp(j/kg K) µ(pa s) κ(w/m K) Y. Iwasa (5/8/3) 65

61 Sample 2: K & Simulation A 6A T [K] A 2A Time [s] V [mv] A 65A Measurement Simulation 45A 2 45A 24A 24A Y. Iwasa (5/8/3) 66

62 Sample 2: Forced-Flow (3 cm/s) & Simulation 12 V [mv] A 65A Measurement Simulation 45A A 24A 24A T [K] A Entrance Middle Exit 45A 24A Time [s] Y. Iwasa (5/8/3) 67

63 Sample 2: Forced-Flow (5 cm/s) & Simulation A 65A Measurement Simulation 71A V [mv] A 45A 65A 24A 45A A 65A Entrance Middle Exit T [K] A 1 24A Time [s] Y. Iwasa (5/8/3) 68

64 YBCO Samples in a Conduction Cooling Setup The sample was mounted between a Kapton-tapeinsulated Cu-block and a G-1 block. It was affixed to the first- stage cold-head of a Cryomech GB-37 cryocooler for conduction cooling to about 4 K. A copper shield was added around the G-1 block to ensure a better uniform temperature of the tape. A heater was affixed to the copper base,and three PRT Thermometers are placed on the YBCO top surface. Slides 69-78: Based on Winston Lue ( ASC22, August 22 ) Y. Iwasa (5/8/3) 69

65 Over-Current Pulse Induced NZP ] A12 [ 9 I ] A [ 1 I ] K [ 65 T t [s] V 1 V 2 V 3 V 4 V 5 V 6 V 7 V I p = A ] 51 K [ 5 T t [s] t [s] I p = A V 1 V 2 V 3 V 4 V 5 V 6 V 7 V t [s] U [ m V ] U [ m V ] At each operating current, I op an initial over-current pulse was applied. The transient pulse energy was varied by changing the magnitude, I p or the duration of the pulse. The figures show at T= 45 K, and at an I op of 33.7-A and 2-s over-current pulse duration, an I p of A caused a quench while A did not. Distinctive normal zone propagation was observed when there was a quench. The inserts show the temperatures measured by a PRT at the center of the sample. Y. Iwasa (5/8/3) 7

66 Stability Margins in the Range 45 8 K g r12 a m 1 y t 8 i l i 6 b a 4 t S 2 3 ] T cs op ( T T )( I I ) = T + 1 c T [K] op Measured t cop To measure the stability margin, the operating current was set near the lowest zone I c at a particular temperature. At the highest recovery pulse current, the V-I product integrated over the initial pulse divided by the zone volume is used to estimate the stability margin. The measured stability margins ranged from 16 to 12 J/cm 3 as temperature lowers from 8 to 45 K. Comparison with the specific heat integrals indicated better fit to T c than to T cs for the HTS tapes. The higher margins measured at lower temperatures are attributed to the higher conduction cooling. Y. Iwasa (5/8/3) 71

67 Minimum Propagation Currents 2 15 ] A [ 1 I ] 12 A [ 9 I 75 7 ] K 65 [ 6 T t [s] V 1 V 2 V 3 V 4 V 5 V 6 V 7 V I op = 27.4 A ] K [ 8 T t [s] t [s] I op = 28.4 A V 1 V 2 V 3 V 4 V 5 V 6 V 7 V ] V m [ U ] V m [ U At 45.4 K, the sample was subjected to an initial overcurrent pulse of 118 A for 2 s while increasing the operating propagation in 1 A steps. The figures show recovery at an operating current of 27.4 A and propagation at 28.4 A - a minimum propagation current. The initial heat input to zone 4 was about 4.6 J or 35.4 J/cm 3 - about 3 times the stability margin t [s] Y. Iwasa (5/8/3) 72

68 Minimum Propagation Current vs. Temperature 25 ] A [ 2 I m p 3 15 measured extrapolated fitted average T [K] Minimum propagation currents were also determined by extrapolating the velocity versus current plots to zero velocity. Minimum propagation currents of 12 to 27 A were observed at 8.8 to 45.4 K. The existence of minimum propagation current is thought to be the result of conduction cooling. Y. Iwasa (5/8/3) 73

69 i t 1 a g a 8 p o r 6 p h c n e u Q I i op [A] t 3 a g I =22.3-A T=76.4-K 25 c2 a I =2.4-A p c3 o 2 I =18.7-A c7 r p v h c 1 n e u 5 Q I c2 =68.8 A I c6 =64.9 A I c7 =61.5 A v 2 v 6 v 7 v 3-2 v 7-8 T=45.8 K I [A] op NZP Velocity vs. I op NPZ always starts from the middle (zone 4) of this sample. NZP velocities were found to be increasing linearly with the current. No more than 2 mm/s of NZP velocity was measured when I op was set below I c. The outer zones which received better cooing resulted in lower NZP velocity and higher minimum propagation current. Y. Iwasa (5/8/3) 74

70 i t a g a p o r p h c n e u Q v 2 v 7 NZP Velocity vs. Temperature I op =3.3 A T [K] 2 1 I 7 c / p I o NZP velocities were measured at five different temperatures. NZP velocities at an operating current of 3.3 A were found at the different temperatures and plotted in the figure. NZP velocity vs. temperature shows a similar dependence as the ratio of operating to critical current. Y. Iwasa (5/8/3) 75

71 Modified Adiabatic Theories of NZP Velocity 1. Constant material properties (Dresner) v A I I mp A C p T c k T op 2. Variable material properties (Zhao and Iwasa) v B I I mp A T c T op k n 1 dk C s dt C n n k n dt n T c T op C s dt 3. Constant material properties with a different transition temperature (Bellis and Iwasa) v C I I m p A k n C n C s T t T op T t T cs 1 2 T c T cs Y. Iwasa (5/8/3) 76

72 i t 2 a g a p15 o r p 1 h c n e 5 u Q NZP velocities: Experiment & Modified Adiabatic Theories v A v B v C v m (Fig. 9) I op =3.3 A T [K] Inclusion of a minimum propagation current is necessary to achieve better agreement with the data. The agreement of the two Iwasa s theories indicated the validity of a new transition point for HTS, Fair agreement is found between the data the modified theory of Iwasa. Y. Iwasa (5/8/3) 77

73 ORNL Conclusion Stability margins ranged from 16 to 12 J/cm 3 as temperature lowers from 8 to 45 K. Comparison with the specific heat integrals indicated better fit to T c than to T cs for the HTS tapes. Distinctive NZP in YBCO was observed at each of the measured temperatures. Minimum propagation current as a function of temperature was measured. Its existence is thought to be the result of conduction cooling. NZP velocities of no more than 2 mm/s were measured. Velocity vs. temperature data show a similar dependence as the ratio of operating to critical current and agrees with a modified adiabatic theory. Y. Iwasa (5/8/3) 78