Thanks again for your suggestions.
After configuring (and running around 1k iterations) as you suggest, the residuals are still not looking good:
There are still some reversed flow at the inlets as well as high values of turbulent viscosity ratio. The min orthogonal quality after the conversion to poly is 0.41. Do you have any other suggestion that may help improving results?