【求教】DES计算时，连续方程残差居高不下，不知何故？（多图）

luoke

当我采用2d-DES计算绕柱流时，连续方程残差居高不下，不知什么原因，
模型是：计算域0.1×0.8m，圆柱d＝0.01m，来流7m/s
我试着将网格加密从几万到几十万，残差降低很小，
调整压力的under－relaxation number也试过，有一点作用，但是如果通过大幅度调整松弛因子，我觉得还是不可取。
不知众位大侠有何高招？请赐教，谢谢！

aguang

可能是流体方程或边界条件设的不对所造成的吧

luoke

左侧圆弧与上下两侧均为速度入口，右端为压力出口，圆柱为壁面
请问这样设置边界条件合理吗？