水沸腾算例UDF的详细解析

1、【强烈推荐】水沸腾算例UDF勺详细解析强烈推荐】水沸腾算例UDF的详细解析#include udf.h #include sg_mphase.h #define T_SAT 373 #define LAT_HT 1.e3/* 包括常规宏 */* 包括体积分数宏 */* 定义饱和温度 */* 定义水蒸汽潜热 */DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn)/* 定义液相源项 */Thread *mix_th, *sec_th;/* 混合相、第二相 定义计算区域指针 */real m_dot_l;/* 蒸发冷凝速率 定义液相质量转移 */mix_th

2、= THREAD_SUPER_THREAD(pri_th);/* 指向混合区的主相即液相的指针 */sec_th = THREAD_SUB_THREAD(mix_th, 1); /*指向单相控制区的气相的指针 */if(C_T(cell, mix_th)=T_SAT)m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*pri_th) - T_SAT)/T_SAT;*/* 如果液相单fabs(C_T(cell,元的温度高于蒸发温度，液相向气相的质量转移定义源项对质量转移偏导 */dSeqn = -0.1*C_R(cell, pri_th)*f

3、abs(C_T(cell, pri_th) - T_SAT)/T_SAT; /*else m_dot_l =0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(T_SAT-C_T(cell,mix_th)/T_SAT;dSeqn= 0.; /* 于是气相向液相转移，所以液相的质量源项对质量转移的偏导为零*/return m_dot_l;DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) /*定义气相源项 */Thread * mix_th, *pri_th;/* 混合相、第一相 */real m_dot_v;mi

4、x_th = THREAD_SUPER_THREAD(sec_th); /* 指向混合区的第二相即气相的指针 */ pri_th = THREAD_SUB_THREAD(mix_th0, );/* 指向单相控制区的液相的指针，液相为主相 */if(C_T(cell, mix_th)=T_SAT)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; /* 如果混合区单 元的温度高于蒸发温度，液相向气相的质量质量转移 */dSeqn = 0.;/* 由于是液相向气相转移，

5、所以气相的质量源项对来自液相的质量转移的偏导为零 */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;dSeqn = -0.1*C_R(cell, sec_th)*fabs(C_T(cell, sec_th) - T_SAT)/T_SAT; /* 由于是气相向液相转移，所以 气相的质量源项对自身的质量转移的偏导不为零 */return m_dot_v;/* 混合模型能量源项 */* 指向混合区的液相的指针 */* 指向混合区的气相的指针 */DEFINE_S

6、OURCE(enrg_src, cell, mix_th, dS, eqn) 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,pri_th) - T_SAT)/T_SAT;/* 如果混合区dSeqn = -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;dSeqn = -0.1*C_VOF(c