Multiphase CFD Applied to Steam Condensation Phenomena in the Pressure Suppression Pool

Transcription

1 Multiphase CFD Applied to Steam Condensation Phenomena in the Pressure Suppression Pool Marco Pellegrini NUPEC STAR Japanese Conference 2016 Yokohama, Japan June 9 th 2016

2 NUCLEAR PLANTS AFFECTED BY THE 3.11 EARTHQUAKE 2 Operating reactor March 11 th 2011 Under inspection Onagawa ~ 130 km JMA seismic intensity Fukushima Daiichi Fukushima Daini

3 STATION BLACK OUT 3 R/B High Pressure Alternate Cooling system Courtesy of A. Obonai, Tohoku Electric Power CO Reactor Core Isolation Cooling system T-quencher Experiment at SIET, Italy (IAE) RCIC quencher Experiment at SIET, Italy (IAE)

4 DIRECT CONTACT CONDENSATION IN S/C 4 MAKE-UP WATER SYSTEMS DIRECT CONTACT CONDENSATION R/B RPV pressure [MPa(abs)] 8.E+06 7.E+06 7.E+06 7.E+06 7.E+06 7.E+06 6.E+06 6.E Time(s) Injection point Computation by A. Buccio (IAE), 2016 ~ 30 m

5 EULERIAN TWO-PHASE FLOW 5 Instantaneous representation Average representation Heat flux Heat flux Energy equation Source terms,

6 HEXAHEDRAL MESH APPLIED TO A SPHERE 6 D D/16 D/32 D/64 D/128 Volume Fraction D Area Density Magnitude of Volume Fraction Gradient 6/10/2016 NURETH-16, Hyatt Regency, Chicago

7 HEXAHEDRAL MESH APPLIED TO A SPHERE 7 D D/16 D/32 D/64 D/128 Error between the computed and theoretical area Error [%] 20.0% 18.0% 16.0% 14.0% 12.0% 10.0% 8.0% 6.0% 4.0% 2.0% 0.0% d/8 d/16 d/32 d/64 d/128 ~ 9% error with large refinement 6/10/2016 NURETH-16, Hyatt Regency, Chicago

8 POLYHEDRAL MESH APPLIED TO A SPHERE 8 D/8 D/16 D/32 Error between the computed and theoretical area Error [%] 3.50% 3.00% 2.50% 2.00% 1.50% 1.00% 0.50% 0.00% d/16 d/32 d/64 d/128 ~ 2.5% error with large refinement 6/10/2016 NURETH-16, Hyatt Regency, Chicago

9 DOMAIN AND MESH STRATEGIES 9 Small nozzle diameter D = 2 mm Large nozzle diameter D = 210 mm Mesh elements: 305,067

10 DOMAIN AND MESH STRATEGIES 10 Small nozzle diameter D = 2 mm Large nozzle diameter D = 210 mm Mesh elements: 305,067 D/16 Mesh elements: 405,067

11 MESH SENSITIVITY MESH x1.25 MESH x1.0 MESH x0.750 MESH x IFA [Mesh X 0.5] IFA [Mesh X 0.75] IFA [Mesh x1.0] IFA [Mesh x1.25] Interfacial area Interfacial area [cm 2 ] Time [ms]

12 DIRECT CONTACT CONDENSATION: CHUGGING 12 Single hole pipe Pressure sensor In recent experiment we employed transparent pipes to visualize the bubble phenomenology during direct contact condensation

13 EXPERIMENTAL EVIDENCE T POOL = C water level 0.2 kg/s 13 Pool temperature [ C] 2.8 m 1.24 m Steam reaching point m 6/10/2016 Severe Accident Mitigation and Research Collaboration

14 DIRECT CONTACT CONDENSATION: CHUGGING-2 14 Multi hole pipe In recent experiment we employed transparent pipes to visualize the bubble phenomenology during direct contact condensation pressure sensor

15 RAYLEIGH-TAYLOR INSTABILITY 15 Accelerating flow field steam steam Psteam water P steam Pwate r P water A Pstea m water Pwater P steam < P water Final terms for area growth a i n t t e t t s k 1 2 6/10/2016 NURETH-16, Hyatt Regency, Chicago

16 IMPLEMENTATION INTO STAR-CCM+ 16 Compressible steam flow Compressible steam flow Record amplitude length at previous time step

17 LARGE NOZZLE DIAMETER: POOLEX 17 pressure outlet velocity inlet T = 106 C v = m/s adiabatic walls D/16 WATER Incompressible Constant properties k-ε standard Temperature = 62 ºC STEAM Compressible Time step = Courant number limited Stopping criteria at interfacial mass transfer (1% of inlet mass flow rate) Mesh elements: 405,067

18 EFFECT OF RTI MODELIZATION 18 Pressure monitor

19 VOLUME FRACTION 19 Minimum area model Steam flow Rayleigh-Taylor Instability Model Steam flow T pool = 62 ºC T pool = 62 ºC

20 Tanskanen, Ph.D. Thesis EXP RTI model No RTI model

21 EFFECT OF MISPREDICTION OF CHUGGING 21 Prediction of oscillating bubble creates thermal stratification in the pool Chugging is responsible for very large mixing in the pool 6/10/2016 NURETH-16, Hyatt Regency, Chicago

22 SMALL NOZZLE DIAMETER: CLERX ET AL. 22 D/16 WATER Incompressible Constant properties k-ε standard Temperature = 25 ºC STEAM Compressible Time step = Courant number limited Stopping criteria at interfacial mass transfer (1% of inlet mass flow rate) Mesh elements: 405,067

23 VOLUME FRACTION FIELD ms 0.6 ms 0.9 Clerx et al., ms 1.5 ms 1.8 ms Bubble implosion is less than 2 ms in the experiment at it appears immediately Minimum area model Rayleigh-Taylor Instability Model

24 CLERX ET AL. EXPERIMENT 24 Clerx et al., 2009 Minimum area model RTI Model Penetration Length [mm] Clerx Experiment BL + no RTI RTI Time [ms]

25 PREDICTION OF TEMPERATURE DISTRIBUTION 25 Minimum area model Measured temperature field RTI Model Clerx et al., 2009

26 THE CHALLENGE OF ACCIDENT COMPUTATION 26 R/B R/B R/B accident time scale [days] Unit 1 vent pipes Unit 2 RCIC Unit 3 RCIC

27 UNIT 1 Fukushima Daiichi power plant what are the conditions at this moment? UNIT 2 UNIT 3 Courtesy of S. Mizokami, TEPCO 27