Another thing, it's quite clear your epsilon gets bounded but that doesn't help, the solver explodes at solving epsilon equation anyway. If the mesh doesn't ruin your day, I'm quite sure y+ somewhere in your domain is well below 30. Again, use yPlus function object and check if that's true