Best Practice Advice AC6-15
Vortex ropes in draft tube of a laboratory Kaplan hydro turbine at low load
Application Area 6: Turbomachinery Internal Flow
Application Challenge AC6-15
Best Practice Advice
Key Fluid Physics
The runner rotation leads to fluctuations due to the interaction of runner blades with guide vanes (rotor–stator interaction). An imbalance of the runner or an asymmetry of the stator also leads to pulsations. These pulsations propagate to the spiral case and can interfere with each other. Rotor–stator interaction is determined by a combination of the numbers of guide vanes and runner blades. Different modes are excited depending on their ratio. The intensity of these pulsations is strongly influenced by the radial clearance between the guide vanes and the runner. The frequency of these pulsations is usually high and it is of the order of several tens of the runner rotation frequency. Karman vortices can be formed on the vanes of the stator, guide vanes or runner blades. Fluctuations associated with Karman vortices can resonate with the vane body. If the body is sufficiently flexible, its vibrations may have the opposite effect on the vortex sheet. The profile of the vane trailing edge has a significant impact on the intensity of these pulsations. Thus, this requires a high precision when manufacturing the trailing edge. The fluctuations associated with the Karman vortices have an even higher frequency, exceeding the runner rotation frequency by a factor of 100 or even more. However, the greatest problems are caused by low-frequency pulsations in the draft tube. There are different vibration – generating mechanisms, some of them are associated with the precession of the vortex core (PVC). Precession of the vortex core occurs under partial load (at a discharge rate of 0.5–0.85 of the optimal value), where the flow has a considerable residual swirl. Vortex breakdown in such a flow leads to the formation of a recirculation zone at the flow axis and the rotation of the corkscrew vortex rope around it. The precession of the vortex rope is a serious danger for the hydraulic turbine equipment. Intense flow-induced pulsations lead to strong vibrations of the hydraulic turbine structure that can lead to damage in case of resonance. Pressure pulsations, generated by a precessing vortex rope, may also affect cavitation processes, enhancing cavitation erosion. Cavitation may enhance all the other pulsations.
The complexity of the geometry, curved and bladed regions, tip-clearance and rotor-stator interaction, oscillation of the runner rotational speed which is absent in numerical simulations are some sources of uncertainties which make a high fidelity CFD model difficult to assemble.
Computational Domain and Boundary Conditions
The computational domain included all the elements constituting the turbine system: guide vane, runner and draft tube. At the input, a fixed flow rate was obtained from the experiment. At the output, a fixed pressure was set. This approach was found to produce realistic results.
A specific issue in the computer simulation of hydraulic turbines is the treatment of the rotating runner and the rotor-stator interaction. As discussed above, several approaches can be found in the literature, i.e. dynamic, sliding and moving grid methods and those based on a moving reference frame. In the present work paper, the modelling of the runner rotation was performed in the rotated reference frame for the runner zone. The obtained results are then rotated with the runner rotation speed and as such imposed as the inflow field for the inlet into the draft tube. Earlier test calculations proved that this approach is credible for describing the integral flow characteristics including the dominant flow pulsations. A comparison with the computationally more demanding method using sliding meshes showed that for this type of flows with a focus on draft tube dynamics the results are almost the same.
Discretisation and Grid Resolution
The convective terms should be discretised by second-order schemes, and the CFL number should be kept below 1.0. Grid issues have been discussed above. It is clear that the 2M grid is too coarse for an LES, but interestingly, the LES on a 6M grid produced the main flow features - the pressure pulsations, velocity profiles and the rms of their fluctuations in good agreement with results for 19.3M cells. This relative insensitivity of LES to the grid refinement above 6M is attributed partially to a carefully designed mesh in the shear layers to meet the resolution criterion based on the “shear” scale (see below), but also to the fact that the flow is dominated by large-scale helical vortex structures, unaffected by the small-scale stochastic turbulence that may have remained unresolved on the 6M grid.
The computations have shown that the two linear eddy viscosity models on the 2M grid could not maintain the unsteadiness and failed to reproduce the essential flow features, but on a finer grid of 6 M resulted in unsteady solutions with improved mean-flow patterns, though still failing to capture the precessing twin-helix structures and the pressure pulsations. In contrast, the second-moment closure (Re-stress) model on the same 2M mesh recovers well both the mean velocity field, the rms of its fluctuations, as well as the specific helical double rope structure in accord with experiments and LES. Similar performance is returned by DES.
The origin of the LEVM underperformance is traced to a too high eddy viscosity as a consequence of an excessive (unconditionally positive) production of the turbulent kinetic energy due to the models' inability to account for the turbulent stress anisotropy and the stress-strain phase lag, both naturally accounted for by the RSM. This leads to a much larger modelled kinetic energy and the effective viscosity compared to the RSM. In flows with a weak internal forcing or using a too coarse grid, LEVM can fail to capture flow unsteadiness. If it does, the resolved motion is weak and manifests in a smaller resolved kinetic energy, which diminishes the models' sensitivity and receptivity to internal flow instabilities compared to RSM (and, of course, to DES where the bulk flow is resolved by LES).
It was also shown that the LES on two very different grids, with 6M and 19.3M cells, returned very similar results. This, to some extent surprising finding, is attributed to the fact that the cell size was optimized to be smaller than the Corrsin shear scale, while resolving the fine stochastic scales proved not to be essential for the flow dominated by strong large-scale coherent vortical structures.
Recommendations for Future Work
The focus of future works should be put on the control of the loss-sources in the draft tube, such as decreasing the pressure pulsations, shrinking the size of the vortex rope and the on-axis recirculation region, and even the prevention of the vortex breakdown. The different control techniques should be investigated to decrease the harmful effects of the vortex breakdown and pressure pulsations. The detailed simulations should be extended to a variety of geometries used in hydropower. Flows under different operating conditions and with different geometrical details should be studied.
- Aakti, B., Amstutz, O., Casartelli, E., Romanelli, G., & Mangani, L. (2015). On the performance of a high head Francis turbine at design and off-design conditions. Francis-99 Workshop 1: steady operation of Francis turbines. Journal of Physics: Conference Series, 579(1), 1-12. doi:10.1088/1742-6596/579/1/012010
- Bosioc, A., Tanasa, C., Muntean, S., & Susan-Resiga, R. (2009). 2D LDV measurements of swirling flow in a simplified draft tube. Conference on Modelling Fluid Flow (CMFF'09), the 14th International Conference on Fluid Flow Technologies.
- Chen, C., Nicolet, C., Yonezawa, K., Farhat, M., Avellan, F., Miyagawa, K., & Tsujimoto, Y. (2010). Experimental study and numerical simulation of cavity oscillation in a diffuser with swirling flow. International Journal of Fluid Machinery and Systems, 3(1), 80-90.
- Kuibin P. A., Litvinov I. V., Sonin V. I., Ustimenko A. S., & Shtork S. I. (2016). Modelling inlet flow in draft tube for different regimes of hydro turbine operation. Vestnik Novosibirsk State University. Series: Physics. 11(1), 56-65.
- Launder, B. E., Reece, G. J., & Rodi, W. (1975). Progress in the development of a Reynolds-stress turbulence closure. Journal of Fluid Mechanics, 68(3), 537-566. DOI: https://doi.org/10.1017/S0022112075001814
- Litvinov, I. V., Shtork, S. I., Kuibin, P. A., Alekseenko, S. V., & Hanjalić, K. (2013). Experimental study and analytical reconstruction of precessing vortex in a tangential swirler. International Journal Heat and Fluid Flow, 42, 251-264.
- Minakov, A. V., Platonov, D. V., Dekterev, A. A., Sentyabov, A. V., & Zakharov, A. V. (2015a). The numerical simulation of low frequency pressure pulsations in the high-head Francis turbine. Computers & Fluids, 111, 97-205.
- Minakov, A. V., Platonov, D. V., Dekterev, A. A., Sentyabov, A. V., & Zakharov, A. V. (2015b). The analysis of unsteady flow structure and low frequency pressure pulsations in the high-head Francis turbines. Int. J. of Heat and Fluid Flow, 53, 183-194.
- Minakov, A. V., Sentyabov, A. V., Platonov, D. V., Dekterev, A. A., & Gavrilov, A. A. (2015c). Numerical modeling of flow in the Francis-99 turbine with Reynolds stress model and detached eddy simulation method. Journal of Physics: Conference Series 579(1), doi:10.1088/1742-6596/579/1/012004.
- Minakov A., Platonov D., Litvinov I., Shtork S., Hanjalić K. (2017) Vortex ropes in draft tube of a laboratory Kaplan hydroturbine at low load: An experimental and LES scrutiny of RANS and DES computational models. Journal of Hydraulic Research 55(3), 668-685
- Nicoud, F., & Ducros, F. (1999). Subgrid-scale stress modelling based on the square of the velocity gradient tensor. Flow, Turbulence and Combustion, 62(3), 183-200.
- Picano, S., & Hanjalić, K., (2012). Leray-α regularization of the Smagorinsky-closed filtered equations for turbulent jets at high Reynolds numbers. Flow, Turbulence and Combustion, 89(4), 627-650.
- Smirnov, A., Shi, S., & Celik, I. (2001). Random flow generation technique for large eddy simulations and particle-dynamics modeling. Journal of Fluids Engineering, 123(2), 359-371.
- Susan-Resiga, R. (2008, June). Hydrodynamic design and analysis of a swirling flow generator. Proceedings of the 4th German – Romanian Workshop on Turbomachinery Hydrodynamics (GRoWTH), Stuttgart, Germany.
- Slotnick, J., Khodadoust, A., Alonso, J., Darmofal, D., Gropp, W., Lurie, E., & Mavriplis, D. (2014) CFD Vision 2030 Study: A Path to Revolutionary Computational Aerosciences. NASA/CR–2014-218178.
- Skripkin, S., Tsoy, M., Shtork, S., & Hanjalić, K. (2016). Comparative analysis of twin vortex ropes in laboratory models of two hydro-turbine draft-tubes. Journal of Hydraulic Research, 54(4), 1-11. http://dx.doi.org/10.1080/00221686.2016.1168325.
- Zadravec, M., Basic, S., & Hribersek. (2007). The influence of rotating domain size in a rotating frame of reference approach for simulation of rotating impeller in a mixing vessel. J. of Engin. Science and Technology, 2(2), 126 – 138
Contributed by: A. Minakov [1,2], D. Platonov [1,2], I. Litvinov , S. Shtork , K. Hanjalić  —
 Institute of Thermophysics SB RAS, Novosibirsk, Russia,
 Siberian Federal University, Krasnoyarsk, Russia,
 Delft University of Technology, Chem. Eng. Dept., Holland.
© copyright ERCOFTAC 2019