The error message ("FATAL ERROR: Maximum number of iterations exceeded: 100") means that the calculation diverged.
I think you have 2 options to fix the problem.
1.
Changing relaxationFactors for each variable to be smaller
in system/fvSolution file.
2.
Changing gradSchemes and divSchemes to be more stable
in system/fvSchemes file.
I think the following settings may be more stable.
- gradSchemes: "faceLimited Gauss linear 1"
- divSchemes for div grad U: "Gauss linear"
- divSchemes for other (U, k, epsilon...): "bounded Gauss upwind"
Good luck with your work.