Investigation of the Time Resolution Set Up Used to Compute the Full Load Vortex Rope in a Francis Turbine

Abstract

The flow in a Francis turbine at full load is characterised by the development of an axial vortex rope in the draft tube. The vortex rope often promotes cavitation if the turbine is operated at a sufficiently low Thoma number. Furthermore, the vortex rope can evolve from a stable to an unstable behaviour. For CFD, such a flow is a challenge since it requires solving an unsteady cavitating flow including rotor/stator interfaces. Usually, the numerical investigations focus on the cavitation model or the turbulence model. In the present works, attention is paid to the strategy used for the time integration. The vortex rope considered is an unstable cavitating one that develops downstream the runner. The vortex rope shows a periodic behaviour characterized by the development of the vortex rope followed by a strong collapse leading to the shedding of bubbles from the runner area. Three unsteady RANS simulations are performed using the ANSYS CFX 17.2 software. The turbulence and cavitation models are, respectively, the SST and Zwart models. Regarding the time integration, a second order backward scheme is used excepted for the transport equation for the liquid volume fraction, for which a first order backward scheme is used. The simulations differ by the time step and the number of internal loops per time step. One simulation is carried out with a time step equal to one degree of revolution per time step and five internal loops. A second simulation used the same time step but 15 internal loops. The third simulations used three internal loops and an adaptive time step computed based on a maximum CFL lower than 2. The results show an influence of the time integration strategy on the cavitation volume time history both in the runner and in the draft tube with a risk of divergence of the solution if a standard set up is used.