利用交错网格simple算法计算顶盖方腔流发散

[我是河滩]

基于有限差分，利用交错网格simple算法计算二维顶盖方腔，层流，计算一直发散。
改变边界条件，相同程序用于poiseuille流和Couette流均收敛，且与解析解能较好符合。因此，怀疑是边界条件的问题，想请教有人知道u网格和v网格的边界条件怎么给吗？
已经改过时间步长和松弛因子但依然发散。

[winsway_zero]

如果能实现顶盖驱动MAC算法，simple也没有问题啊~可以先试试MAC，速度边界条件都是第一类边界条件，其余都为0.只有顶盖有速度。