In this paper, a semi‐discrete rotational velocity‐correction projection method is proposed to solve the Kelvin–Voigt viscoelastic fluid equations. In this method, the rotational velocity‐correction projection method can preserve the divergence free of the velocity
, and the nonlinear equation was linearized. Then, the unconditional stability and optimal convergence order will be provided...