-We performed numerical simulation complemented by experiments of flow in a 60:1 scaled-down laboratory model of a Kaplan hydroturbine at a part load of about 40% nominal flowrate, using two eddyviscosity-(EVM) and a Reynolds stress (RSM) RANS models (realizable k-, k- SST, LRR) and detachededdy-simulations (DES), all on 2M (million) and 6M grid cells, as well as large-eddy simulations (LES) on 6M and 19.3M grids. Unlike the linear EVMs, the RSM, DES (on 2M grid), and LES (on 6M grid) reproduced well the mean velocity components, the rsm of their fluctuations and pressure pulsations in the diffusor draft tube. Despite relatively coarse meshes and insufficient resolution of the near-wall region, LES, DES and RSM also reproduced well the intrinsic flow unsteadiness and the dominant flow structures, all capturing a twin-rope pattern and the associated pressure pulsations in the draft tube. The EVM URANS on 2M grid could not maintain unsteady solution, but did so on a finer grid of 6 M, resulting in improved mean-flow patterns, but still failing to capture the twin-helix structures and the associated pressure pulsations.