phase change udf
If you can't read please download the document
DESCRIPTION
Phase Change UDFTRANSCRIPT
#include "udf.h"#include "sg_mphase.h" #define T_SAT 373 #define LAT_HT 2455e3 DEFINE_SOURCE(vap_src, cell, pri_th, dS, eqn) {Thread *mix_th, *sec_th; real m_dot_v; mix_th = THREAD_SUPER_THREAD(pri_th);sec_th = THREAD_SUB_THREAD(mix_th, 1);/*m_dot_v=0.;*/if(C_T(cell, mix_th)>=T_SAT){m_dot_v = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(C_T(cell, sec_th) - T_SAT)/T_SAT;dS[eqn] = 0.;}if(C_T(cell, mix_th)=T_SAT){m_dot_l = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(C_T(cell, sec_th) - T_SAT)/T_SAT;dS[eqn] = -0.1*C_R(cell, sec_th)*fabs(C_T(cell, sec_th) - T_SAT)/T_SAT; }if(C_T(cell, mix_th)=T_SAT){m_dot = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(C_T(cell, sec_th) - T_SAT)/T_SAT;dS[eqn] = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*LAT_HT/T_SAT;}if(C_T(cell, mix_th)