Validation of Multiphase Flow Modeling in ANSYS CFD

Size: px

Start display at page:

Download "Validation of Multiphase Flow Modeling in ANSYS CFD"

Amice Cobb
5 years ago
Views:

Outline Introduction MPF model validation for

validation MUSIG model Bartolomej testcase (PWR) Lee

transfer (CHT) Summary & Outlook Courtesy by E.

2 Outline Introduction MPF model validation for adiabatic air-water flows Polydisperse MPF model validation MUSIG model Bartolomej testcase (PWR) Lee testcase (BWR) Wall boiling with conjugate heat transfer (CHT) Summary & Outlook Courtesy by E. Krepper (FZD) 2009 ANSYS, Inc. All rights reserved. 2 ANSYS, Inc. Proprietary

3 ANSYS as Part of the German CFD Network in Nuclear Reactor Safety FZ Dresden- ANSYS FZ Karlsruhe Rossendorf Germany GRS Becker Technologies ThAI test facility International Guest Visitors AREVA TU München: TD, Nucl. Energy Univ. Applied Sciences Zittau-Görlitz: IPM Univ. Stuttgart: Nuclear Energy 2009 ANSYS, Inc. All rights reserved. 3 ANSYS, Inc. Proprietary

4 Methodology of CFD Model Development & Validation Experiment phys.-math. Model Complex Geometry 3d CFD Model Complex Flow Conditions Validation Combination with other Models 2009 ANSYS, Inc. All rights reserved. 4 ANSYS, Inc. Proprietary

5 Eulerian MPF Modeling - The Particle Model Mass weighted averaged conservation equations Mass, momentum, energy transport equations for each phase t t r r U k k k k k kl l 1 l k N k r r r P r U U U F I k k k k k k k k k k k Ik F F D F L F WL F TD F V M secondary drag lift wall turbulent mom. transfer lubrication dispersion virtual mass turbulence models for each phase (e.g. k- / k- SST model, 0-eq. disp. phase turb. Model) heat transfer equations for each phase with interfacial transfer closure interfacial i forces need empirical i closure high void fraction effects, bubble induced turbulence, etc ANSYS, Inc. All rights reserved. 5 ANSYS, Inc. Proprietary

6 Lift force, Wall lubrication force & turbulent dispersion Lift force: due to asymmetric wake and deformed asymmetric particle shape Tomiyama C L correlation F C r ( U U ) U C C (Re, Re, Eo) L L G L L G L Wall lubrication force: L L P surface tension prevents bubbles from approaching solid walls Antal, Tomiyama & Frank W.L.F. models F C r U U n n n WL wall G L rel rel W W W Turbulent dispersioni force: 2 C C wall W turbulent dispersion = action of turb. eddies via interphase drag F C r r 3 D tf P F TD F UF UP rp 4 dp rf rp rf (Eo, y/d ) FAD model by Burns et al. (ICMF 04) P 2009 ANSYS, Inc. All rights reserved. 6 ANSYS, Inc. Proprietary

$based on air volume fraction profiles at L/D=59,2 (z=3.$ $flow water velocity J_L [m/s] superficial bubbly flow with near wall void fraction maximum bubbly flow$ $in the transition regime bubbly flow with void fraction maximum at pipe center bubbly bbl flow with$

8 CFX Model Validation MT-Loop & TOPFLOW Test Matrix M01 experimental test series on MT-Loop evaluation based on air volume fraction profiles at L/D=59,2 (z=3.03m) from the sparger system - numerically investigated test case conditions finely disperse bubbly flow water velocity J_L [m/s] superficial bubbly flow with near wall void fraction maximum bubbly flow in the transition regime bubbly flow with void fraction maximum at pipe center bubbly bbl flow with void fraction maximum at pipe center, bimodal slug flow superficial air velocity J_G [m/s] 2009 ANSYS, Inc. All rights reserved. 8 ANSYS, Inc. Proprietary

9 Validation: Bubbly Flows Turbulent Dispersion Force k-eps + RPI TD (0.5) 300 3,00 k-eps + FAD TD. Air Volume Fraction 2,00 SST + RPI TD (0.5) SST + FAD TD Data Model improvement Norm. 1,00 0, Radius, mm 2009 ANSYS, Inc. All rights reserved. 9 ANSYS, Inc. Proprietary

10 Monodispersed Bubbly Flow MT-Loop Test Case FZR-019 FZD-019: J L =1.017 m/s J G =0.004 m/s J G d P =4.8 mm Grace drag Tomiyama lift T./A./F. Wall L. Force FAD Turb. Disp. SST turb. model Sato model t=0.002s 2210 Iterations n [-] d volume fractio normalize Experiment FZR Antal W.L.F., Grid Tomiyama W.L.F., Grid Frank W.L.F., Grid Radius [mm] 2009 ANSYS, Inc. All rights reserved. 10 ANSYS, Inc. Proprietary

11 Monodispersed Bubbly Flow MT-Loop Test Case FZR FZD-052: J L =1.017 m/s J G = m/s 4.0 J G d P =4.4 mm Grace drag Tomiyama lift T./A./F. Wall L. Force FAD Turb. Disp. SST turb. model normalized volume fraction [-] Experiment FZR-052 Antal W.L.F., Grid 2 Tomiyama W.L.F., Grid 2 Frank W.L.F., Grid 2 Sato model t=0.002s Iterations Radius [mm] 2009 ANSYS, Inc. All rights reserved. 11 ANSYS, Inc. Proprietary

TOPFLOW-074 Test Case Conditions from Test

locations: z= 10, 15, 20, 40, 80, 160, 250,

13 TOPFLOW-074 Test Case Conditions from Test Matrix Selection of test case conditions: TOPFLOW-074 test case was subject of validation in the past Superficial velocities: J G = m/s J L =1.017 m/s Wire-mesh sensor measurements at locations: z= 10, 15, 20, 40, 80, 160, 250, 520mm 2009 ANSYS, Inc. All rights reserved. 13 ANSYS, Inc. Proprietary

3d Bubbly Flow Around Obstacle Water Velocity Comparison Comparison CFD Experiment Absolute water velocity distribution in symmetry plane

14 3d Bubbly Flow Around Obstacle Water Velocity Comparison Comparison CFD Experiment Absolute water velocity distribution in symmetry plane Import of exp. data into CFX-Post Pre-interpolation of exp. data to z=0.01m CFD Exp ANSYS, Inc. All rights reserved. 14 ANSYS, Inc. Proprietary

$3d Bubbly Flow Around Obstacle Air Void Fraction Comparison Comparison CFD Experiment Air void fraction$

15 3d Bubbly Flow Around Obstacle Air Void Fraction Comparison Comparison CFD Experiment Air void fraction distribution in symmetry plane CFD Exp ANSYS, Inc. All rights reserved. 15 ANSYS, Inc. Proprietary

16 3d Bubbly Flow Around Obstacle Air Void Fraction Comparison 4) z=40mm 8) z=520mm 8) 7) 6) 3) z=20mm 7) z=250mm 5) 4) 3) 2) z=15mm 6) z=160mm 2) 1) 1) z=10mm 5) z=80mm 2009 ANSYS, Inc. All rights reserved. 16 ANSYS, Inc. Proprietary

$fraction z=520mm$ obstacle y=0mm

18 ANSYS CFX Experimental Data Quantitative Comparison Quantitative data cross sections z=±10, ±15, ±20, ±40, ±80, ±160, ±250, ±520mm: absolute water velocity air volume fraction z=520mm x=-35mm x=+35mm obstacle y=0mm z=-520mm 2009 ANSYS, Inc. All rights reserved. 18 ANSYS, Inc. Proprietary

ANSYS CFX Experimental Data Quantitative Comparison z=-80mm y=0mm Norm. Air Volume Fra action [-] Abso olute Water Ve elocity [m/s] 2.5 2.0 1.5 1.0 0.5 Experiment (z=-80mm) CFX Simulation (z=-80mm) 0.

19 ANSYS CFX Experimental Data Quantitative Comparison z=-80mm y=0mm Norm. Air Volume Fra action [-] Abso olute Water Ve elocity [m/s] Experiment (z=-80mm) CFX Simulation (z=-80mm) x [mm] Experiment (z=-80mm) CFX Simulation (z=-80mm) x [mm] 2009 ANSYS, Inc. All rights reserved. 19 ANSYS, Inc. Proprietary

20 ANSYS CFX Experimental Data Quantitative Comparison z=-20mm y=0mm Norm. Air Volume Fra action [-] Abso olute Water Ve elocity [m/s] Experiment (z=-20mm) CFX Simulation (z=-20mm) Experiment (z=-20mm) x [mm] CFX Simulation (z=-20mm) x [mm] 2009 ANSYS, Inc. All rights reserved. 20 ANSYS, Inc. Proprietary

21 ANSYS CFX Experimental Data Quantitative Comparison z=20mm y=0mm Norm. Air Volume Fra action [-] Abso olute Water Ve elocity [m/s] Experiment (z=20mm) CFX Simulation (z=20mm) x [mm] Experiment (z=20mm) CFX Simulation (z=20mm) x [mm] 2009 ANSYS, Inc. All rights reserved. 21 ANSYS, Inc. Proprietary

ANSYS CFX Experimental Data Quantitative Comparison 3.5 z=80mm 3.0 2.5 y=0mm Norm. Air Volume Fra action [-] Abso olute Water Ve elocity [m/s] 2.0 1.5 1.0 0.

22 ANSYS CFX Experimental Data Quantitative Comparison 3.5 z=80mm y=0mm Norm. Air Volume Fra action [-] Abso olute Water Ve elocity [m/s] Experiment (z=80mm) CFX Simulation (z=80mm) x [mm] Experiment (z=80mm) CFX Simulation (z=80mm) x [mm] 2009 ANSYS, Inc. All rights reserved. 22 ANSYS, Inc. Proprietary

ANSYS CFX Experimental Data Quantitative Comparison z=250mm y=0mm Norm. Air Volume Fra action [-] Abso olute Water Ve elocity [m/s] 2.5 2.0 1.5 1.0 0.5 Experiment (z=250mm) CFX Simulation (z=250mm) 0.

23 ANSYS CFX Experimental Data Quantitative Comparison z=250mm y=0mm Norm. Air Volume Fra action [-] Abso olute Water Ve elocity [m/s] Experiment (z=250mm) CFX Simulation (z=250mm) x [mm] Experiment (z=250mm) CFX Simulation (z=250mm) x [mm] 2009 ANSYS, Inc. All rights reserved. 23 ANSYS, Inc. Proprietary

Polydispersed Bubbly Flow Caused by Breakup &

slug flow: Balance between: coalescence of bubbles

polydisperse bubbly bbl flow counter-current radial

velocity field new population balance model

24 Polydispersed Bubbly Flow Caused by Breakup & Coalescence Transition from disperse bubbly flow to slug flow: Balance between: coalescence of bubbles turbulent bubble breakup bubble size distribution; polydisperse bubbly bbl flow counter-current radial motion of small and large bubbles; more than one velocity field new population balance model (inhomogeneous MUSIG) 2009 ANSYS, Inc. All rights reserved. 24 ANSYS, Inc. Proprietary

25 Inhomogeneous MUSIG Model momentum equations are solved for N gas phases (vel. groups) size fraction equations for M i bubble size classes in each vel. group bubble coalescence and break-up over all M i MUSIG groups N(d P ) breakup/ condensation coalescence/ evaporation d d d d d d d d v 1 v 2 v 3 Velocity group mass transfer d P,krit d P Break up Coalescence 2009 ANSYS, Inc. All rights reserved. 25 ANSYS, Inc. Proprietary

Validation of 3x7 Inhomogeneous MUSIG Model on TOPFLOW-074 good agreement at

$too intensive turbulent dispersion fraction [-] Air volume 25.0 20.00 15.0 10.$ FZR-074, level l I Exp. FZR-074, level L Exp. FZR-074, level O Exp.

0m) CFX 074-A, level C CFX 074-A, level F CFX 074-A, level I CFX 074-A, level L

26 Validation of 3x7 Inhomogeneous MUSIG Model on TOPFLOW-074 good agreement at levels A, L through R too fast spreading of the bubble plume from inlet due to too intensive turbulent dispersion fraction [-] Air volume Exp. FZR-074, level A Exp. FZR-074, level C Exp. FZR-074, level F Exp. FZR-074, level l I Exp. FZR-074, level L Exp. FZR-074, level O Exp. FZR-074, level R CFX 074-A, Inlet level (z=0.0m) CFX 074-A, level C CFX 074-A, level F CFX 074-A, level I CFX 074-A, level L CFX 074-A, level O CFX 074-A, level R x [mm] 2009 ANSYS, Inc. All rights reserved. 26 ANSYS, Inc. Proprietary

27 MUSIG Model Extension Basic population balance equations dn(m, r, t) dt Bubble number density New terms n(m,r,t) m( r,t) n(m, r, t) U(m, r, t)n(m, r,t) t r m t B D B D B B C C Size fraction equations Breakup/Coalescence terms t j ( irdf i) ( ir j duif i) SB S B D S B B S C D S C x i m m m m m m i i i i 1 for evaporation i i 1 i 1 i S i = mi m i i i i 1 i for condensation mi mi 1 mi 1 m i i i i 2009 ANSYS, Inc. All rights reserved. 27 ANSYS, Inc. Proprietary

$9 [K] D inj = 1 [mm] Detailed experimental data: Bubble size distribution Radial steam volume fraction distribution Dirk Lucas, Horst-Michael Prasser: Steam bubble condensation in sub-cooled water in$

29 Condensation Test Case P=2 [MPa] J w =1.0 [m/s] J s =0.54 [m/s] T s =214.4 [ C] T w =210.5 [ C] T w =3.9 [K] D inj = 1 [mm] Detailed experimental data: Bubble size distribution Radial steam volume fraction distribution Dirk Lucas, Horst-Michael Prasser: Steam bubble condensation in sub-cooled water in case of co-current vertical pipe flow, Nuclear Engineering and Design, Volume 237, Issue 5, March 2007, Pages ANSYS, Inc. All rights reserved. 29 ANSYS, Inc. Proprietary

30 Physical Model Setup Standard MUSIG & Extended MUSIG 25 bubble size classes 3 velocity ygroups: 0 3 [mm],3 6 [mm], 6 30 [mm] Arranged in accordance with critical Tomiyama bubble diameter for bubble size dependent lift force Break up model: Luo & Svendsen (F B =0.025) Coalescence model: Prince & Blanch (F C =0.05) 2009 ANSYS, Inc. All rights reserved. 30 ANSYS, Inc. Proprietary

31 TOPFLOW Condensation Testcase Inlet BC Inlet Position WLF TD Force Heat Transfer Config 1 D inj = 4mm Config 2 D inj =4 mm Source Wall Source 75 mm F WLF CTD=1.5 Nu=2+0.15Re p 0.8 Pr 0.5 Config 3 D inj = 1 mm Source F mm WLF CTD=1.5 Nu=2+0.15Re 0.8 p Pr ANSYS, Inc. All rights reserved. 31 ANSYS, Inc. Proprietary

32 Results: Vapor Volume Fraction 2009 ANSYS, Inc. All Config rights reserved. 1 Config 32 2 Config 3 ANSYS, Inc. Proprietary

Results: Radial Steam Distribution 2009 ANSYS,

38 CFD Simulation for Fuel Assemblies in Nuclear Reactors Material Properties Wall Boiling & Bulk Condensation Turbulence Conjugate Heat Transfer (CHT) Multiphase Flow Modeling FSI: Stresses & Deformations Validation against Experiments 2009 ANSYS, Inc. All rights reserved. 38 ANSYS, Inc. Proprietary

Condensation Turbulence Conjugate Heat Transfer

Deformations Validation against Experiments 2009

39 CFD Simulation for Fuel Assemblies in Nuclear Reactors Material Properties Wall Boiling & Bulk Condensation Turbulence Conjugate Heat Transfer (CHT) Multiphase Flow Modeling FSI: Stresses & Deformations Validation against Experiments 2009 ANSYS, Inc. All rights reserved. 39 ANSYS, Inc. Proprietary

40 Multiphase Flow Regimes for Boiling Water Flow subcooled flow bubbly flow slug flow annular flow spray flow T ONB OSB wall temperature T sat mean fluid temperature subcooled boiling nucleate boiling (saturated boiling) x 2009 ANSYS, Inc. All rights reserved. 40 ANSYS, Inc. Proprietary

41 Flows with Subcooled Boiling (DNB) RPI-Wall Boiling Model Mechanistic wall heat partioning model: q Wall q F q Q q E convective heat flux q F A1 hf ( TW TL) quenching heat flux q Q A2 hq ( TW TL) evaporation heat flux q m (h h ) E G L y A 2 A1 A 2 Convective heat flux u Quenching heat flux m * * * m m t 2009 ANSYS, Inc. All rights reserved. 41 ANSYS, Inc. Proprietary

42 RPI-Wall Boiling Model Submodels for Model Closure Submodels for closure of RPI wall boiling model: Nucleation site density: Lemmert & Chawla, User Defined Bubble departure diameter: Tolubinski & Kostanchuk, Unal, Fritz, User Defined Bubble detachment frequency: Terminal rise velocity over Departure Diameter, User Defined Bubble waiting time: Proportional to Detachment Period, User Defined Quenching heat transfer: Del Valle & Kenning, User Defined Turbulent Wall Function for liquid convective heat transfer coefficient Correlation for bulk flow mean bubble diameter required: e.g. Kurul & Podowski correlation via CCL Supported combination of wall boiling & CHT in the solid GGI & 1:1 solid-fluid interfaces 2009 ANSYS, Inc. All rights reserved. 42 ANSYS, Inc. Proprietary

45 The Bartolomej Test Case R = 7.7 mm Variable Value Z= 2 m q=0.5 57MW/m 2 P 4.5MPa R 7.7 mm G in 900 kg/(s m2) q q 0.57MW/m2 Subcooling 58.2 K G in =900 kg/(s m 2 ) 2009 ANSYS, Inc. All rights reserved. 45 ANSYS, Inc. Proprietary

46 Multiphase Flow Model Steam-Water 2-phase flow: Water: continuous phase Water Steam: disperse bubbles (particle model) Material properties (EOS): IAPWS-IF97 water - water steam property tables Modified law for interfacial area Kurul & Podowski type bulk bubble diameter: d B =f(t sub ) Accounting for higher volume fraction of the steam phase Turbulence Model SST turbulence model for continuous phase 0-eq. disperse phase turb. model + Sato bubble induced turbulence 2009 ANSYS, Inc. All rights reserved. 46 ANSYS, Inc. Proprietary

47 Inter-Phase Mass, Momentum and Energy Transfer Mass transfer model Thermal Phase Change Model (bulk boiling/condensation model) RPI wall boiling model Momentum transfer models Grace drag FAD turbulent dispersion force Tomiyama lift force Wall lubrication force (none, Antal, Tomiyama) Heat transfer models Water: Thermal Energy Water Steam: Saturation temperature Two resistance model Ranz Marshall correlation for bubble heat transfer 2009 ANSYS, Inc. All rights reserved. 47 ANSYS, Inc. Proprietary

48 Numerical Grids Validation on mesh hierarchy with regular refinement factor of 4 (2d meshes) Grid Grid1 Grid2 Grid3 # Nodes (uniform) 20x150 40x300 80x600 Max y Δt [s] x ANSYS, Inc. All rights reserved. 48 ANSYS, Inc. Proprietary

53 Comparison to Experimental Data - Parameter & Model Variation Influence of wall heat flux: Influence of wall lubrication force model: 2009 ANSYS, Inc. All rights reserved. 53 ANSYS, Inc. Proprietary

55 Lee et al. (2008) Testcase Axially symmetric circular annulus Radial dimensions Inner radius of outer tube: R = mm Outer radius of inner tube: R 0 = 9.5 mm r Outlet R R0 R C Measuring Plane (for experimental and numerical Results) Core radius: R C = 3/4 R 0 Annulus width: 9.25 mm Axial dimensions Total heating section height: ht L T = 1670 mm Distance between inlet and measuring plane: L M = 1610 mm L M L T Heated Wall Adiabatic Wall Radial Position: R P Dimensionless, radial distance from inner tube (R P = 0) to outer tube (R P = 1) across the annulus: R P r R R R 0 0 Inner Tube (Heating Rod) Annulus Outer Tube 2009 ANSYS, Inc. All rights reserved. 55 ANSYS, Inc. Proprietary Axis Inlet z

56 Geometry and Mesh Geometry & Mesh generation: Figure 4: Generated one layer mesh Figure 1: 25 segment of the geometry Annulus Core Shell Figure 2: 1 simulated segment of the geometry symmetry symmetry Figure 3: 1 segment with one-layer mesh example symmetry symmetry 2009 ANSYS, Inc. All rights reserved. 56 ANSYS, Inc. Proprietary

57 Mesh Hierarchy Mesh Name Grid 01 (coarse) Grid 02 (medium) Grid 03 (fine) Domains ( 1 = HFO, 2 = CHT ) * No. of Nodes 1: : : : : : No. of Elements 1: 20x150 1: 40x300 1: 80x600 (hexahedra) 2: 40x150 2: 80x300 2: 160x600 y + max (at 1 st node near wall) Set16 ~84 ~41 ~24 Set25 ~88 ~45 ~25 Tstep Δt [s] Set Set ANSYS, Inc. All rights reserved. 57 ANSYS, Inc. Proprietary

58 Selection of Extreme/Limiting Testcase Conditions Concentrating on 2 (out of 12) datasets: Set 25 (least of all steam) Set 16 (most of all steam) Parameter comparison Set No.* q [kw m^-2] G [kg m^-2s] T in [ C] P in [kpa] ANSYS, Inc. All rights reserved. 58 ANSYS, Inc. Proprietary

59 Required Parameter Modifications in Comparison to PWR Conditions Found that submodels need modifications for BWR conditions (see also Tu&Yeoh, Anglart et al., Krepper, Koncar): 1. Bulk bubble diameter (BBD) Kurul & Podowski d B,max wall modified d B law d B,max 2. Bubble departure diameter (BDD) Tolubinski & Kostanchuk d W ~0.5mm max. const. bubble dept. diam. d W =1mm - 3mm 3. A 2 - Wall area fraction influenced by steam bubbles default increased up to ANSYS, Inc. All rights reserved. 59 ANSYS, Inc. Proprietary

60 BBD & BDD Modifications Test Matrix Overview Trying to systematically increase Bubble Departure Diameter to investigate its influence on Heat Flux to Vapor (Q V ) profile Test series with increasing BDD starting from d W.max 0.5 mm 1 mm; 2 mm; 3 mm T&K * 4.0 BDD Tolubinsky & Kostanchuk BDD User defined d W [mm] K&P yes - bbdmod01-1 = const. bbdmod02-2 = const. bbdmod03-3 = const ANSYS, Inc. All rights reserved. 60 ANSYS, Inc. Proprietary

61 BBD Modification / Set 25: Gas Volume z = 1610 [mm] Set25 : Bulk Diameter Modification Comparison: Gas Volume Fraction (r G ) r G [ ] Radial Position (R p ) [ ] G3_K&P G3_dbmod01 G3_dbmod02 G3_dbmod03 Experimental Data 2009 ANSYS, Inc. All rights reserved. 61 ANSYS, Inc. Proprietary

62 BDD Modification / Set 25: Gas Volume z = 1610 [mm] Set25 Bubble Departure Diameter Modification Comparison: Gas Volume Fraction (r G ) r G [ ] Radial Position (R p ) [ ] G3_K&P G3_dbmod02_bddmod01 bddmod01 G3_dbmod02_bddmod02 bddmod02 G3_dbmod02_bddmod03 Experimental Data 2009 ANSYS, Inc. All rights reserved. 62 ANSYS, Inc. Proprietary

63 A 2F Limiter Modification: Results Set 25, Gas Volume z=1610[mm] rg [ ] 0,45 0, ,35 0,30 0,25 0, ,15 0,10 0,05 0,00 Set25 A 2F Mod Comparison: Gas Volume Fraction (r G ) ,1 02 0,2 03 0,3 04 0,4 05 0,5 06 0,6 07 0,7 08 0,8 09 0,9 1 G3_K&P G3_dbmod02_bddmod02_A2F02 Experimental Data Radial Position (R p ) [ ] G3_dbmod02_bddmod01_A2F02 G3_dbmod02_bddmod03_A2F ANSYS, Inc. All rights reserved. 63 ANSYS, Inc. Proprietary

64 Grid Independency: Results Set 25, Gas Volume z = 1610 [mm] Set25 New Grid Comparison: Gas Volume Fraction (r G ) rg [ ] 0,40 0, ,30 0,25 0,20 0,15 0,10 0,05 0,00 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1 Radial Position (R p ) [ ] G1_dbmod02_bddmod02_A2F02 G3_dbmod02_bddmod02_A2F02 Experimental Data G2_dbmod02_bddmod02_A2F02 G4_dbmod02_bddmod02_A2F ANSYS, Inc. All rights reserved. 64 ANSYS, Inc. Proprietary

66 Heat Source in Solid Material & Conjugate Heat Transfer Prediction HFO (Heat Flux Only): Fluid Domain (Annulus) area specific heat flux boundary condition # Outlet Fluid Domain Inlet CHT (Conjugated Heat Transfer): Fluid Domain (Annulus) + Solid Domain (Non-Heated Rod Shell) + + Solid Domain (Heated Rod Core) volume specific heat source # Outlet Solid Domain # Solid Domain Fluid Domain Inlet 2009 ANSYS, Inc. All rights reserved. 66 ANSYS, Inc. Proprietary

67 The RPI Wall Boiling Model: Lee et al. Testcase with CHT Specific energy source in solid material, Set25 (equiv. to q Wall ): E Core = [W/m 3 ] Temperature and Steam VF distribution in vertical plane 2009 ANSYS, Inc. All rights reserved. 67 ANSYS, Inc. Proprietary

68 The RPI Wall Boiling Model: Lee et al. Testcase with CHT Set25 & CHT: Water temperature monitors x W =1.5mm, z=83.5mm, T L z=835mm T L Outlet Iterations T L Intlet Iterations 2009 ANSYS, Inc. All rights reserved. 68 ANSYS, Inc. Proprietary

69 The RPI Wall Boiling Model: Lee et al. Testcase with CHT Set25 & CHT: Grid independence for temperature z=1610[mm] Heated Core Fluid Domain Unheated Cladding 2009 ANSYS, Inc. All rights reserved. 69 ANSYS, Inc. Proprietary

71 New R&D Consortium R&D Initiative: Modeling, Simulation & Experiments for Boiling Processes in Fuel Assemblies of PWR Karlsruhe Inst. of Technology (KIT) ANSYS Germany TUD, Dept. Fluid Mechanics FZ Dresden/ Rossendorf TUM, Dept. Thermodynamics TUD, Dept. Nucl. Eng. Univ. Bochum, Dept. Energy Systems Univ. Appl. Sciences Zittau/ Görlitz TUD Medical Faculty 2009 ANSYS, Inc. All rights reserved. 71 ANSYS, Inc. Proprietary

72 Modeling, Simulation & Experiments for Boiling Processes in Fuel Assemblies of PWR Ultrafast electron beam X-ray CT of fuel rod bundle in titanium pipe on FZD: Images by courtesy of U. Hampel, FZD 2009 ANSYS, Inc. All rights reserved. 72 ANSYS, Inc. Proprietary

73 Modeling, Simulation & Experiments for Boiling Processes in Fuel Assemblies of PWR Wall boiling simulation in a 3x3 rod bundle with spacer grid: Wall superheat T W -T Sat 2009 ANSYS, Inc. All rights reserved. 73 ANSYS, Inc. Proprietary

74 Summary & Outlook Overview on ANSYS CFD multiphase flow model development and validation Continuous effort in model improvement, R&D Emphasis in validation on BPG, comparison to data, geometry & grid independent modeling High interoperability of physical models Outlook: Ongoing & customer driven CFD model development Research cooperation with Industry & Academia More & more complex MPF phenomena Coupling of wall boiling model to inhomogeneous MUSIG Extension of the wall heat partitioning in wall boiling model 2009 ANSYS, Inc. All rights reserved. 74 ANSYS, Inc. Proprietary

Development and Validation of the Wall Boiling Model in ANSYS CFD

Development and Validation of the Wall Boiling Model in ANSYS CFD Th. Frank, C. Lifante, A.D. Burns PBU, Funded CFD Development ANSYS Germany Thomas.Frank@ansys.com E. Krepper, R. Rzehak FZ Dresden-Rossendorf