求UDF高手帮忙修改动量源项

lixinjun911

本人正在做风雨流场的多相流，因为要考虑在建筑壁面处雨相可以穿透，才能求出雨在壁面的抓取率（catch ratio），多方求助得到一份udf,希望能够将非稳态改成稳态计算，烦请各位懂UDF的帮助一下，谢谢大家。
  F_AREA(A,f,tf);
    m_A=sqrt(A[0]*A[0]+A[1]*A[1]+A[2]*A[2]);

    if ((A[2]/m_A*(-1.0)==1.0)&&(C_W(c,t)==0.0))
                                            
{  
V_z=-9.8*CURRENT_TIMESTEP;
V_x=C_U(c,t);
V_y=C_V(c,t);
}
else
{V_y=C_V(c,t);
V_x=C_U(c,t);
V_z=C_W(c,t);
  }
if (fabs(V_x)<1.0e-50)
{V_x=0.0;}
                       
if (fabs(V_y)<1.0e-50)
{V_y=0.0;}       
if (fabs(V_z)<1.0e-50)
{V_z=0.0;}       
   m_v=sqrt(V_x*V_x+V_y*V_y+V_z*V_z);
if (m_v==0.0)
{m_v=1.0;}
alfa=A[0]/m_A*V_x/m_v+A[1]/m_A*V_y/m_v+A[2]/m_A*V_z/m_v;
beta_x=A[0]/m_A;
beta_y=A[1]/m_A;
beta_z=A[2]/m_A;
   if((A[0]/m_A*V_x/m_v+A[1]/m_A*V_y/m_v+A[2]/m_A*V_z/m_v)>=0.0)
{
source =source-C_R(c,t)*C_VOF(c,t)/CURRENT_TIMESTEP*alfa;
C_UDMI(c,tm,0) =C_UDMI(c,tm,0)+ source*C_U(c,t);
C_UDMI(c,tm,1) =C_UDMI(c,tm,1)+ source*C_V(c,t);
C_UDMI(c,tm,2) =C_UDMI(c,tm,2)+ source*C_W(c,t);
C_UDMI(c,tm,3) =C_UDMI(c,tm,3)+ source*C_U(c,tp);
C_UDMI(c,tm,4) =C_UDMI(c,tm,4)+ source*C_V(c,tp);
  C_UDMI(c,tm,5) =C_UDMI(c,tm,5)+ source*C_W(c,tp);
  }
  else
  {
  source = source+0.0;
}
期待有人能够解答，如何将公式中的时间项修改，使得计算变成稳态，非常感谢。

mayo_0

把时间相关的宏替换成1或者其他值即可。

lixinjun911

好的 试试  谢谢啊