[sara2115]

Hi.
I have written a code for flow field around a square cylinder.
The base of my code is written in vorticity stream function approach.
For calculation of stream function on the cylinder I have used "single pressure magnitude" for multi-connected domains. now it works properly. but after I decided to insert two slot on the top and down edge of square it dose not work any more. the down slot sucks the fluid and the upper slot blow the fluid by constant velocity. I was wondering is there any one who know how i can determine the right magnitude of stream function on the cylinder?
thanks