HIGH HEAT FLUX AND CRITICAL HEAT FLUX (CHF) TESTS IN SUPPORT TO ITER HHFC

Transcription

1 HIGH HEAT FLUX AND CRITICAL HEAT FLUX (CHF) TESTS IN SUPPORT TO ITER HHFC by F.Escourbiac With the contributions of : J.Schlosser, N.Vignal, V.Cantone, M.Missirlian, M.Richou, R.Guigon * Association EURATOM-CEA, France J.L. Jouvelot, S.Constans, I.Bobin * AREVA-NP industry, Le Creusot, France J.Boscary * Association FZJ, Germany M.Merola * ITER Organization, Cadarache, France B. Riccardi * Fusion For Energy, Barcelona, Spain Experimental work was supported by the European Communities Slide 1

2 INTRODUCTION Vplasma= 25 m3 In powerfull fusion machines, the HHFC have to sustain HHF of MW/m2 R=2.715m R=2.22m When cooled by water in the subcooled boiling regime, typically -40 bar, C, the critical heat f lux is a concern In addition with thermal design rules on armour or joint temperatures (not reported in this talk), it is necessary to estimate the margin to critical heat flux for any design of HHFC ITER VT Objective of this talk : to come back on the notion of CHF to review the EU studies and to assess the CHF margin of the ITER Divertor CEA-EURATOM Slide 2

3 OUTLINE THE CHF TESTS A LONG CHF TESTING STORY CHF RESULTS HEAT TRANSFER EVALUATION CHF PREDICTION CHF margin : APPLICATION TO THE ITER DIVERTOR CONCLUSIONS Slide 3

4 OUTLINE THE CHF TESTS A LONG CHF TESTING STORY CHF RESULTS HEAT TRANSFER EVALUATION CHF PREDICTION CHF MARGIN : APPLICATION TO THE ITER DIVERTOR CONCLUSIONS Slide 4

5 THE CHF TESTS Principle Sweeping frequency: up to 10 khz (b) 90% of the power in a < 1.5 mm disc. diameter (a) 4 mm 4 mm (c) Main facility for CHF tests: FE200 (electron beam 200 kw) in Le Creusot, France (flexible pressurized water loop : bar, up to 6 kg/s, C) Slide 5

6 THE CHF TESTS Increase power step by step * measure Absorbed Heat Flux (AHF) PW Coolant circuit High Heat Flux 23 2 * Infrared monitoring 790 C Up to sudden temperature excursion T Pyr2 ( C) Shot t (s) Low-temperature pyrometers vs time during one shot (absorbed flux in MW/m² is given in red) (destruction of component in a few 100 ms if gun not stopped) 520 Infra-red picture at steady state 970 C 525 Infra-red picture, transient OUTPUT : Absorbed CHF (ACHF) Slide 6

7 THE CHF TESTS CHF is often used for Burn-out, ACHF, ICHF, WCHF More precisely, when Burn-out is detected during HHF testing : * ACHF : Absorbed Critical Heat Flux is measured by calorimetry ACHF = QC p (T out -T in )/(w x L) Heat Flux Heat Flux * ICHF : Incident Critical Heat Flux is derived from ACHF taking account radiation (width) Tin L (length) Tout ICHF = ACHF + Radiated Heat Flux * WCHF : Wall Critical Heat Flux is calculated by Finite Element modeling WCHF = ICHF x Pf ICHF s WCHF 0 WHF=F(s) s(mm) πr Slide 7

8 OUTLINE THE CHF TESTS A LONG CHF TESTING STORY CHF RESULTS HEAT TRANSFER EVALUATION CHF PREDICTION CHF MARGIN : APPLICATION TO THE ITER DIVERTOR CONCLUSIONS Slide 8

9 A LONG CHF TESTING STORY Investigation of various smooth/swirl tubes diameters 1996 Investigations of cooling concept (swirl rod insert, double smooth and swirl tubes, HV tubes) CEA 3a-5b CEA 3a-3,4,5 CEA 3a-1,2 AF1 ID/OD 11/16 AF3 ID/OD 18.6/ CEA 3e-1,2,3 1.5 CEA 3b 2.5 CEA 3f-1,2 2.5 DSM2 10 ID DST2/DST4 Tr = 2 and 4 HV1/HV (copper alloy tubes : Glidcop, CuCrZr) Slide 9

10 A LONG CHF TESTING STORY 1997 : Study of twist ratio and tape thickness ST2 ST : monoblock PRODIV1 4.5 CFC material 23 Tape thick.: 2 mm Twist ratio: 2, 3 or 4 23 Tape thick.: 0.8 mm Twist ratio: 2, 3 or CFC monoblock, Glidcop tube Tape thickness: 2 mm Twist ratio: : Study on hypervapotron as an alternative solution to the monoblock concept : investigation of CFC monoblocks mock-ups width and tube diameter /12 ID/OD 10/12 ID/OD 10/12 ID/OD 12/15 ID/OD > 200 CHF values have been obtained Slide 10

11 OUTLINE THE CHF TESTS A LONG CHF TESTING STORY CHF RESULTS HEAT TRANSFER EVALUATION CHF PREDICTION CHF MARGIN : APPLICATION TO THE ITER DIVERTOR CONCLUSIONS Slide 11

12 CHF RESULTS Influence of thermalhydraulics parameters ICHF (MW/m²) Série9 Série8 ST2 5 m/s ST2 10m/s ST2 15m/s ST4 5m/s ST4 10m/s ST4 15m/s 15 m/s 10 m/s 5 m/s ST4 CHF with T sat -T bulk (K) (subcooling) * Means : CHF with Tsat (water pressure) * Means : CHF if Tbulk (water temperature) ICHF (MW/m²) Série8 Série7 Série9 HV1 5m/s HV1 10m/s HV1 15m/s HV3 5m/s HV3 10m/s HV3 15m/s 15 m/s 10 m/s 5 m/s HV T sat -T bulk (K) CHF with water flow rate (The lines show the CEA-Tong75 correlation described later) Slide 12

13 CHF RESULTS Influence of the swirl tape design : twist ratio CHF results Mock-up Geometry Twist Ratio ICHF (MW/m 2 ) flat peaked Thick tapes 2 ST ST ST Thin tapes 2 ST ST ST (3.5 MPa, Tsub,out=100 C, V=12m/s, ID = 10 mm) CHF if twist ratio Twist ratio 2 Twist ratio 3 Twist ratio 4 Slide 13

14 CHF RESULTS Influence of the swirl tape design : thickness CHF results Mock-up Geometry Twist Ratio ICHF (MW/m 2 ) flat peaked Thick tapes 2 ST ST ST Thin tapes 2 ST ST ST (3.5 MPa, Tsub,out=100 C, V=12m/s, ID = 10 mm) CHF if twist ratio CHF with thickness tape mm Twist ratio 2 Twist ratio 3 Twist ratio mm Slide 14

15 CHF RESULTS Influence of the incident heat flux axial profile CHF results Mock-up Geometry Twist Ratio ICHF (MW/m 2 ) flat peaked Thick tapes 2 ST ST ST Thin tapes 2 ST ST ST (3.5 MPa, Tsub,out=100 C, V=12m/s, ID = 10 mm) CHF if twist ratio CHF with thickness tape CHF with peaked axial heat flux Incident Heat Flux (MW/m²) mm Twist ratio 2 Twist ratio 3 Twist ratio mm 20 ITER 15 Transient Reference Case 10 FE200 Peaked Profile 5 X(cm) Slide 15

16 CHF RESULTS Influence of geometry ICHF ( MW/m²) Subcooling ~120 C Axial Velocity ~4 m/s peak. unif Width (mm) ICHF (MW/m 2 ) Width (mm) Subcooling ~100 C Axial velocity ~10 m/s ID/OD 10/12 uniform ID/OD 10/12 peaked ID/OD 12/15 uniform CHF if width /12 ID/OD 10/12 ID/OD 10/12 ID/OD 12/15 ID/OD Slide 16

17 OUTLINE THE CHF TESTS A LONG CHF TESTING STORY CHF RESULTS HEAT TRANSFER EVALUATION CHF PREDICTION CHF MARGIN : APPLICATION TO THE ITER DIVERTOR CONCLUSIONS Slide 17

18 HEAT TRANSFER EVALUATION Development of correlations Convective regime : Sieder Tate Nu = f ST x Re 0.8 H Pr 1/3 (µ b /µ w ) 0.14 h con = k Nu / D h, f ST =1.15 if twisted tape; 1 if not Subcooled boiling regime : Thom-CEA T w - T sat = (10-6 Φ w ) e (-P.10-5 /87) 1/2.8 Comparison between the different correlations using B. & R. method Slide 18

19 HEAT TRANSFER EVALUATION Validation of subcooled boiling regime correlation with swirl tubes 100 ST22 ST22R ST24 2 TH7 thermocouples TH8 Jens & Lottes 1/4 Thom-CEA 1/2.8 wall heat flux ( MW/m 2 ) 10 ST082 Thom 1/2 1 T wall - T sat (K) Slide 19

20 HEAT TRANSFER EVALUATION Validation of full set of correlations with swirl tubes 2 TH7 thermocouples TH Comparison between thermocouples and FE calculations temperature ( C ) TH7 2 mm from heated surface TH8 7 mm from heated surface 150 Step wise increase of EB gun power up to CHF Int. HHFC Workshop, Dec , 2008, incident UCSD, heat flux ( CA, MW / m 2 ) Slide 20

21 HEAT TRANSFER EVALUATION Validation of heat transfer evaluation : hypervapotron 6.6 Tsurf ( C) TSEFEY-M 1500 FE MW/m² T( C) bar, 120 C, 9 m/s AHF (MW/m²) Tth ( C) Slide 21

22 OUTLINE THE CHF TESTS A LONG CHF TESTING STORY CHF RESULTS HEAT TRANSFER EVALUATION CHF PREDICTION CHF MARGIN : APPLICATION TO THE ITER DIVERTOR CONCLUSIONS Slide 22

23 CHF PREDICTION F TONG75H calculated from CEA-TONG75 correlation Bo Hcrit =1.84 (D H /D o ) Re H -0.6 [ (P/Pc) 1.8 Re H -0.5.Ja] F TONG75H = Bo Hcrit.ρ.V.i fg (Re H = VD H / ν taking account swirl presence) WCHF : estimated with application of a corrective factor extrapolated from tests WCHF = C f. F TONG75H C f : depends on the geometry, the promoter and the heat flux profile (range 1,2 2,2, precision +/- 20%) CHF margin = WCHF /max.whf max.whf being calculated with design values Slide 23

24 OUTLINE THE CHF TESTS A LONG CHF TESTING STORY CHF RESULTS HEAT TRANSFER EVALUATION CHF PREDICTION CHF MARGIN : APPLICATION TO THE ITER DIVERTOR CONCLUSIONS Slide 24

25 CHF MARGIN : APPLICATION TO THE ITER DIVERTOR CHF is a concern near the strike point on the CFC monoblocks (lower part of the VT) Heat loads : * Nominal : 10 MW/m², 400 sec. * Transient : 20 MW/M², 10 sec. Design Erosion during lifetime Potential Imperfection of bonding (defect) CFC or W θ θ Cu CuCrZr ID/OD 10/12 mm θ : Location, θ : Extension CFC, W or Cu : Interface Slide 25

26 APPLICATION TO ITER DIVERTOR CHF margin estimation max WHF (MW/m²) (radiation included) CFC IHF Defect CHF thickness Margin 18 mm 7 mm max Wall Heat Flux (loading 20 MW/m²; 10 sec.) CFC 7 mm ; loc. 45 CFC 7 mm ; loc. 0 CFC 18 mm ; loc. 45 CFC 18 mm ; loc Extension defect θ 10 MW/m², SS No 3,5 yes 3,8 20 MW/m², 10 No 2,2 sec yes 2,0 10 MW/m², SS No 3,5 yes 3,3 20 MW/m², 10 No 1,6 sec yes 1,4 Defect centred at 0 can protect from CHF! (K) OUTPUT : Lowest margin to CHF = 1.4 Slide 26

27 OUTLINE THE CHF TESTS A LONG CHF TESTING STORY CHF RESULTS HEAT TRANSFER EVALUATION CHF PREDICTION CHF MARGIN : APPLICATION TO THE ITER DIVERTOR CONCLUSIONS Slide 27

28 CONCLUSIONS The actively cooled plasma facing components have to sustain HHF of MW/m2 Cooled by water in the subcooled boiling regime, the critical heat flux is a concern Several CHF studies in Europe during since more 15 years investigated various cooling concepts at various cooling conditions. More than 200 values of CHF have been obtained. Critical Heat Flux prediction and heat transfer correlations were developed during these studies This methodology already used for the design of the HHFC of water cooled controlled fusion machines such as Tore Supra and W7-X Stellerator was applied to ITER Vertical target divertor, this gives a CHF margin of 1.4, the current state of the art Slide 28

29 Additional information Slide 29

30 ZOOM ON NUKIYAMA S BOILING CURVE Different design, different limits NUKIYAMA S BOILING CURVE Swirl tube ACHF ~ local CHF ( no transition boiling, except at low velocity [Marshall]) Hypervapotron ACHF > local CHF ( there is transition boiling in normal conditions) Forced Convection Wall Heat Flux Subcooled Nucleate Boiling Partial Boiling ONB Fully Developed Boiling local CHF CHF Wall Temperature Transition Boiling Film Boiling Swirl Hypervapotron tube Slide

31 FE200 FACILITY ELECTRON BEAM : 180 kw at steady state (200 kw, t<10s) advanced sweeping system (10 khz, +/-10 ) 0.1 sec. to 6000 sec. VACUUM CHAMBER 8 m3, mockups up to 2 m WATER LOOP : heat exchanger 0 kw Pressurized water loop (3.9 MPa) 50 to 2 C, 6 kg/s, 0.7 MPa DIAGNOSTICS : IR, CCD devices 2 pyrometers thermocouples data acquisition Co-operated by CEA Framatome Since 1992 in Le Creusot, France more than 100 components tested Slide 31

32 DEFINITIONS i fg : latent heat of vaporization [J.kg -1 ] P c : critical pressure of water (221 bars) Ja = (ρ f /ρ g ).[(i sat -i)/i fg ]: Jacob number µ: dynamic viscosity [kg.m -1.s -1 ] D 0 : reference diameter (12.7 mm = 0.5 inch) Slide 32

33 COMPARISON OF SWIRL TUBE AND HYPERVAPOTRON CHF (MW/m²) (P = 35 bar, Dtsubout = 120 C) DP (bar/m) Higher pressure drop for swirl tube but it is about only 21% of the total pressure drop Both swirl and HV tubes meet the ITER requirements kg/s/m 2 Slide 33