Numerical simulation of the effect of rod bowing on critical heat flux