有大侠做过vof的气液两相流中的open channel flow问题吗？

linjun0528

我想用steady-state implicit 公式模拟一个二维的自由表面流动，下面是水，水深20m,长200米，速度1m/s，上面是air, 空气高度5m。我做了一个长方形的流动域，x方向是水流方向，在fluent中选用2D,segrated,steady,vof(implicit)，在solve-solution定义了body force weighted,,PISO,2-order scheme.并这样定义了边界条件，在长方形流域的左边界分成上下两边，上面的边用作air inlet(velocity inlet:1m/s),下面的边用作water inlet(velocity-inlet:1m/s),同样在长方形流域右边定义air out(pressure outlet)和water outlet(velocity inlet:1m/s),出口高度与进口相同，长方形流域上边定义 pressure-inlet,长方形流域下边定义为wall.然后我定义了水相和气相，最后pacth水相，但是不收敛。也试过右边定义velocity inlet:-1m/s，也不行。
最近拿到fluent6.2，发现vof中新增了计算open channel flow的功能，而且能做steady问题，但是我怎么试都很快发散；unsteady问题相对好一点，但是时间步长取到0.0005还是会在气液界面处中段（或进口位置）出现较大流速，我想可能是两个问题：1、上界面的边界取法；2、k\e值，我取k=0.0025,e=0.00002934因为经常会出现提示：turbulent viscosity limited to viscosity ratio of 1.0000000e+005 in……,这句话不知是什么意思，确实我对这个问题中怎么取水里半径不清楚？
各位大侠请多指教，有经验的能否上传一个成功的case?

linjun0528

水相计算有点像了，但就是不收敛，入口处的水气界面出现了较大流速,各位给点建议，谢谢！

whuxyl

在王福军-<计算流体动力学分析-CFD软件原理与应用>书中有类似的例子:含自由表面的河流跌坎分析.在本论坛中我曾贴过这个例子的CAS文件,你可以搜索一下