it can be calculated by using equation based modelling.
therefore, the stream function is:
(psi)xx+(psi)yy=-(vx-uy).
so we can used the left hand term as source term in COMSOL.
I already did that but it was not work correctly something still missing???