[sara2115]

HI everyone.
I've wrote a code to simulate velocity field over a square cylinder.
I've used vorticity-stream function solver at 0 angle of attack.(not immersed)
but now I don't know how to specify outflow boundary condition for different angles of attack.
I used a rectangular as my domain.
please help me.
thank you for reading my thread.

[mettler]

I did this for a finite difference method (except I did flow over a step). My exit boundary condition was @x=L, Streamfunction(@x=L) = Streamfunction(@x=L-1). I did this based on assuming that the flow after the step was fully developed and not changing. My vorticity BC was the same. The results were good.

[sara2115]

Originally Posted by mettler

I did this for a finite difference method (except I did flow over a step). My exit boundary condition was @x=L, Streamfunction(@x=L) = Streamfunction(@x=L-1). I did this based on assuming that the flow after the step was fully developed and not changing. My vorticity BC was the same. The results were good.

Thank you for your reply Mettler.
but in my case I have vortex shedding which makes it more complicated.
By the way, did free stream velocity have an angle of attack in your case? or it was simply horizontal?

[mettler]

it had no angle of attack, but why can't you extend the rectangular grid long enough to where the vortex shedding has dissipated and the flow is fully developed?

[sara2115]

Originally Posted by mettler

it had no angle of attack, but why can't you extend the rectangular grid long enough to where the vortex shedding has dissipated and the flow is fully developed?

Dear Mettler.
Well, I did so.
but still there are lot's of problems.
for example, even if the flow behind the cylinder is fully developed, it still has its angle of attack, so I can't impose the boundary condition you used for outflow.
there are other difficulties in dealing with far field boundaries.
thank you for your reply.