Steady state and non-linear harmonic (NLH) flow simulations were performed within the framework of the Francis-99 project in order to assess the capacity of the NLH method to capture the main pressure fluctuations associated with the rotor-stator interactions between the distributor and the runner of the turbine. This paper focusses on the methodology developed to obtain harmonic solutions and presents preliminary results from the simulations using the flow solver NUMECA FineTURBO on intermediate grid level meshes. Comparisons of the first simulations to experimental data reveal good agreement concerning the predicted pressure amplitudes notably at high load operating condition.