We present an adaptive wavelet method to integrate the velocity–vorticity formulation of the two-dimensional Navier–Stokes equations coupled with a penalisation technique to handle efficiently solid boundaries of arbitrary shape. We demonstrate the validity of this method, called coherent vortex simulation (CVS), to compute the flow around an impulsively started cylinder at high Reynolds number and compare the results with a classical vortex method.