Hi Dears
I have confused a n important question: for walls of channel boundary condition Can we implement Zou-He with(u=v=0) instead of bounce back : (f(i)=f(opposite(i))??
6 2 5
\ | /
3----0----1
/ |
7 4 8
for example for upper wall Zou-He BC yields:
f(2)=f(4)
f(5)=f(7)-0.5*(f(1)-f(3))
f(6)=f(8)+0.5*(f(1)-f(3))
–buonce back yields:
f(2)=f(4)
f(5)=f(7)
f(6)=f(8)
I tried two above methods for a uniq problem. the results are different (with the first BC, flow reachs fully developed in shorter length)and it seems that the first method is more accurate. what is your opponion??
What about half bounce back? I don’t know how to programmed half bounce back. Does anyone have a part of code showing how to imlement half bounce back??