封闭容器，储液晃荡，为什么不收敛‘大家帮帮忙

我行

在模型中设定舱壁静止、液体运动的初始状态.液舱长25 m ,高16 m ,液面高度12 m. 用粘性不可压的水模拟LNG,用可压缩的空气模拟自由液面上的舱内环境. 一个大气压取为101. 325 kPa ,液体的初始速度为10m/s，时间步长取为2 ms. 采用VOF 法进行模拟计算. 计算过程中采用标准k2ε湍流模型,假定区域内的流场是非定常湍流
编写了液体初始速度的udf，如下；
#include "udf.h"
DEFINE_INIT(my_init_function,mixture_domain)
{
Thread **pt;
Thread *cell_threads;
cell_t cell;
mp_thread_loop_c(cell_threads,mixture_domain,pt)
{
  begin_c_loop_all(cell,cell_threads)
  {
            C_U(cell,cell_threads)=10.0;/*速度，可以修改，注意要小数点*/
  }
  end_c_loop_all(cell,cell_threads)
}
}
但是残差图就是不收敛
下面附有msh、cas文件





[ 本帖最后由 我行 于 2009-3-29 02:19 编辑 ]

我行

大家帮帮忙

我行

大家帮帮忙呀