IAEA CRP benchmark of ROCOM PTS test case for the use of CFD in reactor design using the CFD-Codes ANSYS CFX and TrioCFD

Abstract Over the last 15 years, considerable effort has been expended in assembling the available information on the use of CFD in the nuclear reactor safety field. Typical application areas here are heterogeneous mixing and heat transfer in complex geometries, buoyancy-induced natural and mixed convection, etc., with specific reference to Nuclear Reactor Safety (NRS) accident scenarios such as Pressurized Thermal Shock (PTS), boron dilution, hydrogen build-up in containments, thermal fatigue and thermal striping issues, etc. The development, verification and validation of CFD codes in respect to Nuclear Power Plant (NPP) design necessitates further work on the complex physical modelling processes involved, and on the development of efficient numerical schemes needed to solve the basic equations. Therefore, a set of ROCOM CFD-grade test data were made available to set up an International Atomic Energy Agency (IAEA) benchmark, relating to PTS scenarios. The benchmark deals with the injection of the relatively cold Emergency Core Cooling (ECC) water, which can induce buoyancy-driven stratification. Data obtained from the PTS experiment were compared in the study presented here with predictions obtained from the CFD software packages ANSYS CFX and TrioCFD. In addition a test case without buoyancy forces was selected to show the influence of density differences. The results of the two test cases and the numerical calculations show that mixing efficiency is strongly influenced by buoyancy effects. At higher mass flow rates without density differences the injected slug propagates in the circumferential direction around the core barrel. Buoyancy forces reduce this azimuthal propagation. The ECC water falls in an almost vertical path and reaches the lower downcomer sensor below the inlet nozzle. Therefore, density effects play an important role during natural convection with ECC injection in PWRs. Both CFD codes were able to predict the observed flow patterns and mixing phenomena.

[1]  Frank-Peter Weiss,et al.  Coolant mixing in a PWR - de-boration transients, steam line breaks and emergency core cooling injection - experiments and analyses , 2002 .

[2]  F. Nicoud,et al.  Subgrid-Scale Stress Modelling Based on the Square of the Velocity Gradient Tensor , 1999 .

[3]  H. Reichardt,et al.  Vollständige Darstellung der turbulenten Geschwindigkeitsverteilung in glatten Leitungen , 1951 .

[4]  P. Quemere,et al.  Verification and validation considerations regarding the qualification of numerical schemes for LES for dilution problems , 2010 .

[5]  B. Launder,et al.  Progress in the development of a Reynolds-stress turbulence closure , 1975, Journal of Fluid Mechanics.

[6]  P.-E. Angeli,et al.  FVCA8 Benchmark for the Stokes and Navier–Stokes Equations with the TrioCFD Code—Benchmark Session , 2017 .

[7]  Thomas Höhne,et al.  Experimental and Numerical Modeling of Transition Matrix From Momentum to Buoyancy-Driven Flow in a Pressurized Water Reactor , 2009 .

[8]  Frank-Peter Weiss,et al.  Coolant Mixing in a Pressurized Water Reactor: Deboration Transients, Steam-Line Breaks, and Emergency Core Cooling Injection , 2003 .

[9]  Soeren Kliem,et al.  Experiments at the mixing test facility ROCOM for benchmarking of CFD codes , 2008 .

[10]  H. Prasser,et al.  A new electrode-mesh tomograph for gas–liquid flows , 1998 .

[11]  Soeren Kliem,et al.  Fluid mixing and flow distribution in the reactor circuit, measurement data base , 2005 .

[12]  D. Kuzmin,et al.  High-resolution FEM-TVD schemes based on a fully multidimensional flux limiter , 2004 .

[13]  P. Raviart,et al.  Conforming and nonconforming finite element methods for solving the stationary Stokes equations I , 1973 .

[14]  Soeren Kliem,et al.  Hydrodynamic phenomena in the downcomer during flow rate transients in the primary circuit of a PWR , 2007 .

[15]  Thomas Höhne,et al.  Modeling of a buoyancy-driven flow experiment at the ROCOM test facility using the CFD codes CFX-5 and Trio_U , 2006 .