I’ve already partially succeeded, with very low radial velocity of the inner cylinder,
the Couette flow can be recovered quite well, but the system looses density perpetually.
I have also understood that interpolated bounce-back does not guarantee mass conservation. Have you benchmarked your code by simulating some standard examplex?