A second order turbulence model based on a Reynolds stress approach for two-phase boiling flow. Part 1: Application to the ASU-annular channel case

摘要

High-thermal performance PWR (pressurized water reactor) spacer grids require both low pressure loss and high critical heat flux (CHF) properties. Numerical investigations on the effect of angles and position of mixing vanes and to understand in more details the main physical phenomena (wall boiling, entrainment of bubbles in the wakes, recondensation) are required.rnIn the field of fuel assembly analysis or design by means of CFD codes, the overwhelming majority of the studies are carried out using two-equation eddy viscosity models (EVM), especially the standard K-ε model, while the use of Reynolds Stress Transport Models (RSTM) remains exceptional.rnBut extensive testing and application over the past three decades have revealed a number of shortcomings and deficiencies in eddy viscosity models. In fact, the K-ε model is totally blind to rotation effects and the swirling flows can be regarded as a special case of fluid rotation. This aspect is crucial for the simulation of a hot channel in a fuel assembly. In fact, the mixing vanes of the spacer grids generate a swirl in the coolant water, to enhance the heat transfer from the rods to the coolant in the hot channels and to limit boiling.rnFirst, we started to evaluate computational fluid dynamics results against the AGATE-mixing experiment: single-phase liquid water tests, with Laser-Doppler liquid velocity measurements upstream and downstream of mixing blades. The comparison of computed and experimental azimuthal (circular component in a horizontal plane) liquid velocity downstream of a mixing vane for the AGATE-mixing test shows that the rotating flow is qualitatively well reproduced by CFD calculations but azimuthal liquid velocity is underestimated with the K-ε model.rnBefore comparing performance of EVM and RSTM models on fuel assembly geometry, we performed calculations with a simpler geometry, the ASU-annular channel case. A wall function model dedicated to boiling flows is also proposed.