文档介绍:Forpersonaluseonlyinstudyandresearch;mercialuse薃对飞昂公司水-水蒸汽蒸发过程UDF程序段的解释薄#include""‘包括常规宏衿#include""‘包括体积分数宏CVOF(C,T)莆#_SAT373‘定义蒸发温度100℃薆#‘定义蒸发潜热J/Kg蚃芀肈莅DEFINE_SOURCE(liq_src,cell,pri_th,dS,eqn)‘液相质量源项UDF螃{蚁蒆Thread*mix_th,*sec_th;‘定义计算区线指针肄realm_dot_l;‘定义液相质量转移kg/()袃螈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=-*C_VOF(cell,pri_th)*C_R(cell,pri_th)*蚆fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;袆羃dS[eqn]=-*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;蚀‘定义源项对质量转移偏导,莇蚄}肃else{肀m_dot_l=*C_VOF(cell,sec_th)*C_R(cell,sec_th)*袅fabs(T_SAT-C_T(cell,mix_th))/T_SAT;蒃‘如果指向混合区液相的单元温度小于蒸发温度,气相向液相的质量转移,液相得膃蒁dS[eqn]=0.;‘由于是气相向液相转移,所以液相的质量源项对质量转移的偏导为零薇蒆}芃薈returnm_dot_l;艿芅}莂罿DEFINE_SOURCE(vap_src,cell,sec_th,dS,eqn)‘气相质量源项UDF螇羄{蒂莀Thread*mix_th,*pri_th;葿realm_dot_v;螃蒂mix_th=THREAD_SUPER_THREAD(sec_th);‘指向混合区的第二相即气相的指针螁袇pri_th=THREAD_SUB_THREAD(mix_th,0);指向单相控制区的液相的指针,液相为主相螆薂if(C_T(cell,mix_th)>=T_SAT)‘如果混合区单元的温度高于蒸发温度,液相向气相的质量质量转移袈蕿{薅m_dot_v=*C_VOF(cell,pri_th)*C_R(cell,pri_th)*蚂fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;艿dS[eqn]=0.;‘由于是液相向气相转移,所以气相的质量源项对来自液相的质量转移的偏导为零肆}莃else{螂m_dot_v=-*C_VOF(cell,sec_th)*C_R(cell,sec_th)*虿fabs(T_SAT-C_T(cell,mix_th))/T_SAT;螈‘如果指向混合区的单元温度小于蒸发温度,气相向液相的质量转移,气相失莆袂dS[eqn]=-*C_R(cell,sec_th)*fabs(C_T(