Hi everyone,

I have been trying to create on open exit boundary in my matlab code. I have looked at the paper on Neumann/Dirchlet boundary conditions and tried all kinds but have had no success so far.

My mesh velocity directions are…

% 8 4 7

% \ | /

% \ | /

% \ | /

% 1 <-------o-------> 3

% / |

% / |

% / |

% 5 2 6

f3(:,N)=f3(:,N-1);

f6(:,N)=f6(:,N-1);

f7(:,N)=f7(:,N-1);

and so i have been trying boundary conditions such as those just above based on what i have read on the Neumann boundary condition about creating a ghost and having the same velocities at these nodes (N-1) and the boundary nodes (N). If anbody has any help or ideas then i would really appreciate it.

Thanks. James.