260 likes | 770 Views
FNST/MASCO/PFC Meeting. Boiling Heat Transfer in ITER First Wall Hypervapotrons. Dennis Youchison, Mike Ulrickson and Jim Bullock Sandia National Laboratories Albuquerque, NM August 6, 2010.
E N D
FNST/MASCO/PFC Meeting Boiling Heat Transfer in ITER First Wall Hypervapotrons Dennis Youchison, Mike Ulrickson and Jim Bullock Sandia National Laboratories Albuquerque, NM August 6, 2010 Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company,for the United States Department of Energy’s National Nuclear Security Administration under contract DE-AC04-94AL85000. 1
Outline • What are hypervapotrons? • Why hypervapotrons? • Geometry optimization • Boiling heat transfer in hypervapotrons • Why CFD? • Benchmarking with HHF test data • CHF prediction 2
Background • Star-CCM+ Version 5.04.006, User Guide, CD-adapco, Inc., New York, NY USA (2010). • S. Lo and A. Splawski, “Star-CD Boiling Model Development”, CD-adapco, (2008). • D.L. Youchison, M.A. Ulrickson, J.H. Bullock, “A Comparison of Two-Phase Computational Fluid Dynamics Codes Applied to the ITER First Wall Hypervapotron,” IEEE Trans. On Plasma. Science, 38 7, 1704-1708 (2010). • Upcoming paper in the 2010 TOFE . 3
Why hypervapotrons? • Advantages: • High CHF with relatively lower pressure drop • Reduction in E&M loads due to thin copper faceplate • Lower Cu/Be interface temperature (no ss liners) • Less bowing of fingers due to thermal loads • Disadvantages: • CuCrZr/SS316LN UHV joint exposed to water 5
What are hypervapotrons? Hypervapotron FW “finger” 6
Two-phase CFD in water-cooled PFCs • Problem: conjugate heat transfer with boiling • • Focus on nucleate boiling regime below critical • heat flux • • Use Eulerian multiphase model in FLUENT & Star-CCM+ • • RPI model (Bergles&Rohsenow) • • Features heat and mass transfer between liquid • and vapor, custom drag law, lift or buoyancy and influence of bubbles on turbulence • CCM+ transitions to a VOF model for the film when vapor fraction is high enough – need to know when to initiate VOF 7
5 MW/m2 400 g/s t=2.05s Velocity distributions Drag on bubbles, lift or buoyancy, changes in viscosity and geometry, all affect the velocity distribution under the heated zone. 2mm-deep teeth and 3-mm spacing optimized to produce a simple reverse eddy in the groove. 8
Star-CCM+ 560 k polyhedra mesh Switches from Eulerian multi-phase mixture to VOF for film boiling. 9
Star-CCM+ Results CCM+ boiling models were benchmarked against US and Russian test data for rectangular channels and hypervapotrons to within 10oC. Case analyzed is a hot “stripe” on a section of the ITER first wall. Surface temperature distribution, t=6.3 s capability to predict CHF from CFD 10
With no boiling, heat transfer is highest under the fins With boiling, the vapor fraction in grooves is 4%-6% on average Star-CCM+ Results Case analyzed is a hot “stripe” on a section of the ITER first wall. The details of the heat transfer change dramatically as boiling ensues. Iso-surface of 2% vapor volume fraction t=6.3 s 11
Star-CCM+ gives same h as Fluentfor nucleate boiling. Heat transfer coefficients increase in grooves where boiling takes place ranging from 12,000 to 13,000 W/m2K. 12
Systematic parameter study performed on rectangular channels – then applied to hypervapotrons. 13
Thermocouple response 3.5 MW/m2 through 6 s Thermocouple response 4.0 MW/m2 through 6 s Russian data ICHF Trip @ 400 C Temperature (C) Not ss yet! Temperature (C) Rectangular channel results 14
Russian HV CHF Mock-up flow 15
Total of 490k poly cells in mesh 3 prism layers Heated area is 100 mm x 48 mm 16
Surface temperature – 6.0 MW/m2, 1 m/s 115 C inlet, 2 MPa 17
CCM+ solid/fluid interface temperatures for 6.0 MW/m2 @6s 18
Thermocouple response through 6 s 4 s for TCs to ss Russian data 20
All flow regimes can exist simultaneously. T: h: • sub-cooled • nucleate to transition boiling • film boiling • sub-cooled 4.0 MW/m2 115 oC, 2 MPa water 1.0 m/s 22
CHF Testing Testing of the HV mock-up T/C (1.5 mm from CuCrZr surface) Water 2 m/s Pabs 10 MW/m2 tpuls 300s Second pulse at 10 MW/m2) ICHF ! 23