谁能帮我解释一下FLUENT-HELP湿蒸汽udf

mumu8888

[这个贴子最后由mumu8888在 2008/05/20 11:03am 第 1 次编辑]

[br][br][以下内容由 mumu8888 在 2008年05月20日 09:50am 时添加] [br]
real
wetst_satP(real T)
{
  real psat  ;
  real SUM=0.0;
  real pratio;
  real F  ;
  real a1  = -7.41924200      ;
  real a2  =  2.97210000E-01  ;
  real a3  = -1.15528600E-01  ;
  real a4  =  8.68563500E-03  ;
  real a5  =  1.09409899E-03  ;
  real a6  = -4.39993000E-03  ;
  real a7  =  2.52065800E-03  ;
  real a8  = -5.21868400E-04  ;
  if (T > ws_Tc) T = ws_Tc ;
  F   = ws_aaa*(T - ws_TPP)   ;
  SUM = a1 + F*(a2+ F*(a3+ F*(a4+ F*(a5+ F*(a6+ F*(a7+ F*a8)))))) ;
  pratio = (ws_Tc/T - 1.0)*SUM;
  psat   = ws_Pc *exp(pratio) ;
  return psat; /*Pa */
}


real
wetst_satT(real P, real T)
{
  real  tsat  ;
  real  dT, dTA,dTM,dP,p1,p2,dPdT;
  int  i ;
    for (i=0; i<25; ++i)
      {
if (T > ws_Tc) T = ws_Tc-0.5;
p1= wetst_satP(T)        ;
p2= wetst_satP(T+0.1)    ;
dPdT = (p2-p1)/0.1 ;
dP = P - p1 ;
dT = dP/dPdT   ;
dTA = fabs(dT);
dTM = 0.1*T  ;
if (dTA > dTM)  dT=dT*dTM/dTA ;
T = T + dT;
if (fabs(dT) < TEMP_eps*T) break;
      }
  tsat = T;
  return tsat; /*K */
}
real
wetst_eosP(real rho, real T)
{
  real P  ;
  P = rho* Rgas_v * T ;
  return P; /*Pa */
}