重力场计算

Hoijuon

计算重力场出现图中的情况 请问造成边界处和刚体处静压的可能原因？

静压

相关设置：
2D（3D也会出现这样的情况）瞬态 重力已加载
标准k-epsilon
单相 水
操作压力点在右上方 具体值根据水下压强计算
边界条件：中间的是固定wall  四周的是压力出口 UDF指定如下（上表面定义为压力零点）

DEFINE_PROFILE(pressure_profile,t,i)
{
    real x[ND_ND];    /* this will hold the position vector */
    real y;
    face_t f;
    begin_f_loop(f,t)
      {
       F_CENTROID(x,f,t);
       y = x[1];
           if(y<0.020)
           {
         F_PROFILE(f,t,i) = 0.0+998.2*9.81*(0.020-y);
           }
           else
           {
             F_PROFILE(f,t,i) = 0.0;
           }
      }
    end_f_loop(f,t)
}

SIMPLE算法
其他为默认

THANKS