I understand that you are not using the forcing scheme from Lamorgese but a velocity sort of feedback scheme every time step. I was wondering where you got the idea for that implementation and whether it was I guess physically correct