The flow field \mathbf{u} computed in the first time step seems to be calculated incorrectly and contain (mostly) NaN values, see picture below.
\mathbf{u}
NaN