fluent 中 模拟蒸发 的UDF

mix_th) - T_SAT)/T_SAT; dS[eqn] = 0.; //‘由于是液相向气相转移，所以气相的质量源项对来自液相的质量转移的偏导为零 } else { m_dot_v =

-0.1*C_VOF(cell, sec_th)*C_R(cell,

sec_th)* fabs(T_SAT-C_T(cell,mix_th))/T_SAT;//‘如果指向混合区的单元温度小于蒸发温度，气相向液相的质量转移，气相失 dS[eqn] = -0.1*C_R(cell, sec_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;//由于是气相向液相转移，所以气相的质量源项对自身的质量转移的偏导不为零 } return

m_dot_v;}DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn) //‘混合模型能量源项UDF{ Thread *pri_th, *sec_th; real m_dot; pri_th =

THREAD_SUB_THREAD(mix_th, 0);//‘指向混合区的液相的指针 sec_th = THREAD_SUB_THREAD(mix_th, 1); //‘指向混合区的气相的指针 if(C_T(cell, mix_th)>=T_SAT) //‘如果混合区的单元温度高于蒸发温度。质量转移由液相向气相转移，吸热，质量转移量前有负号{ m_dot = -0.1*C_VOF(cell, pri_th)*C_R(cell,

pri_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; dS[eqn] = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT; } else { // 相反，气相向液相转移则放热 m_dot = 0.1*C_VOF(cell, sec_th)*C_R(cell,

sec_th)* fabs(T_SAT-C_T(cell,mix_th))/T_SAT; dS[eqn] = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;} return LAT_HT*m_dot; //气化潜热与质量转移率相乘得能量源项W/m3}

Injection

The droplet type is available when heat transfer is being modeled and at least two chemical species are active or the nonpremixed or partially premixed combustion model is active. You should use the ideal gas law to define the gas-phase density when you select the droplet type

#include \

#include \#define T_SAT 373

#define LAT_HT 2.253e+6 (水潜热）

DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) {

Thread *mix_th, *sec_th; real m_dot_l;

mix_th = THREAD_SUPER_THREAD(pri_th); sec_th = THREAD_SUB_THREAD(mix_th, 1);

if(C_T(cell, mix_th)>=T_SAT) 如果液相的温度大于饱和温度 {

m_dot_l = -4220*C_VOF(cell, pri_th)*C_R(cell, pri_th)*

fabs(C_T(cell, mix_th) - T_SAT)/ LAT_HT; （4220为水的比热，LAT_HT为水潜热）

dS[eqn] = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;

（dS这相我不明白什么意思，但估计前面的0.1也应该换成4220）

}

else { 应该是气态转换成液体，需要计算气态所含多余的热量} m_dot_l = 2028*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(T_SAT-C_T(cell,sec_th))/ LAT_HT; （2028蒸汽比热）

dS[eqn] = 0.; （不明白什么意思） }

return m_dot_l; }

DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) {

Thread * mix_th, *pri_th; real m_dot_v;

mix_th = THREAD_SUPER_THREAD(sec_th); mix-th指针代表气相 pri_th = THREAD_SUB_THREAD(mix_th, 0); pri-th指针代表液相

if(C_T(cell, mix_th)>=T_SAT) 如果气相温度大于373（估计这里需要改成液相的） {

m_dot_v = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*

fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; 计算转换成气态的水，公式应该不对，后面一些也是这样的问题，

当液相和饱和蒸气相平衡共存时，在单位时间内从单位液体表面蒸发出去的分子的数目（蒸发的分子通量 ）应该和在单位时间内凝结到单位液体表面的分子的数目（凝结的分子通量 ）相等。如果忽略饱和蒸气分子碰到液体表面时的反射，认为碰到液体的分子都要凝结，就可以用饱和蒸气的分子通量 J来表示在单位时间内凝结到单位液体表面的分子的数目（凝结的分子通量 ）。u可以称为液体的蒸发速率，其实就是在单位时间内从单位液体表面蒸发出去的液体的体积，它等于蒸发的分子通量除以液体的分子数密度，也等于凝结的分子通量除以液体的分子数密度。设饱和蒸气可以视为理想气体，饱和蒸气的压强为 p、热力学温度为T 、分子的平均速率为v、分子数密度为ng、密度为ρg，分子的质量为m，摩尔质量为μ，液体的分子数密度为nl、密度为ρl，普适气体常量为R，圆周率为π，则有 u=J/nl=ngv/(4nl) =mngv(4mnl) =ρgv/(4ρl) =[μp/(RT)]

=[8RT/(πμ)]^(1/2)/(4ρl) =(p/ρl)[μ/(2πRT)]^(1/2)

对于20℃的水，已知T=293K,p=2.34*10^3Pa,ρl =1.00*10^3kg*m-3，μ=18.0*10^-3 kg*mol-1，R=8.31J*mol^-1*K-1. 将这些数据代入上式得 u =(p/ρl)[μ/(2πRT)]^(1/2)

?=(2.34*10^3/1.00*10^3)[18.0*10^-3/(2π*8.31*293)]^1/2 m*s^-1

?≈2.54*10^-3m*s^-1 ≈15.2 cm*min^-1.