Am I wrong to think that the pressure coefficient from the output of a transient incompressible solver (pimpleFoam) should be calculated as: $$C_p = \frac{p-p0}{0.5 U_{\infty}^2}$$ as opposed to $$C_p = \frac{p-p0}{0.5 \rho U_{\infty}^2}$$ with the consideration that the pressure values are kinematic pressure (p/rho)? The mean Cp value I get at the stagnation point of a suspended square cylinder is getting higher than 1. Any help is appreciated.