Error from solver k_epsilon_eq.f, solver crashes

Hello,
I am simulating a single flow case. The case works well without k-e turbulent model. But it reports error (solver crash!error from solve_k_epsilon.f:375, error divide by zero). I also tried several k-e model parameters, but all failed. There is a post about this problem but the discussions did not provide solution. Anyone could help me with this?
bonuli2.zip (24.5 MB)

I assume the model runs without the k-e model on?

The k-e model is selected.