Numerical integration of an arbitrary initial boundary problem for Navier-Stokes equations in is considered. The method uses conditionally stable semi-Lagrangian advection scheme for invicid fluxes with high order interpolation on unstructured mesh. The fractural step method on physical processes is used to project pressure field. In order to increase time step integration the BFECC method with limiter TVD correction is used. Flow over a Cylinder for 3D case is compared to the experimental and simulation data.