水蒸发至湿空气FLUENT UDF教学（内含UDF）

#include "udf.h"

#include "mem.h"

#define MOLAR_MASS_WATER 18.01534

#define MOLAR_MASS_AIR 28.97

DEFINE_PROPERTY(saturation_temp, c, t)

{

real vol = C_VOLUME(c, t);

Thread *pt = THREAD_SUB_THREAD(t, 0);

Thread *st = THREAD_SUB_THREAD(t, 1);

real vf_s = C_VOF(c, st);

real vf_p=C_VOF(c,pt);

real p_mix = C_P(c, t);

Pareal p_op = RP_Get_Real("operating-pressure");

real rho_s = C_R(c, st);

real rho_p=C_R(c,pt);

real p_w;

real t_sat; m_sec = rho_s * vol * vf_s;

real m_pri=rho_p*vol*vf_p;

real m_wv = C_YI(c,st,0) * m_sec;

real m_air = m_sec - m_wv;

real N_wv = m_wv / MOLAR_MASS_WATER;

real N_air = m_air / MOLAR_MASS_AIR;

real N_total = N_wv + N_air;

if(fabs(vf_p-1)<0.001)

{

p_w = C_P(c, t)+ p_op;

}

else if(C_YI(c,st,0)>0.0000001)

{

p_w = (C_P(c, t)+ p_op)* (N_wv / N_total);

}

else

{

p_w=0;

}

C_UDMI(c,t,0)=p_w;

if (p_w>0)

{

t_sat = (1730.63 / (10.196 - log10(p_w))) + 39.724;

}

else

{

t_sat=273.15;

}

C_UDMI(c,t,1)=t_sat-273.15;

C_UDMI(c,t,2)=m_pri/(m_pri+m_sec);

return t_sat;

}

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删