资源描述:
《fluent-中-模拟蒸发-的UDF.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、#include"udf.h"/H'z7x%F"L8~%{2G&e#include"sg_mphase.h"7s0k+@3K!P1N9H+/_#defineT_SAT373%x;c!E"f)t5O,w#defineLAT_HT1.e36a.I0{7P5o+K3O!C,l;_-e,{('ZDEFINE_SOURCE(liq_src,cell,pri_th,dS,eqn)3l:N5g(l3o{)W-@;e b%b,b Thread*mix_th,*sec_th;%I(p-N(~,S*p(X realm_dot_l;6a*F7L0G"Z&_,s m
2、ix_th=THREAD_SUPER_THREAD(pri_th);+M'D,Y-F9_4i1x4s%M F8F sec_th=THREAD_SUB_THREAD(mix_th,1);&l:i;V7}:h if(C_T(cell,mix_th)>=T_SAT){$l'i1k8t5s/_"K;R0K m_dot_l=-0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*,^'I1r#e$p#_)_ fabs(C_T(cell,pri_th)-T_SAT)/T_SAT;+a(i"@1k;D#J-~ dS
3、[eqn]=-0.1*C_R(cell,pri_th)* M"j H3g+P*e0Z fabs(C_T(cell,pri_th)-T_SAT)/T_SAT;#i,j&~7H6k5Y9e }-o7M2e3e+ X*B1V,
4、$I else{-M3v+t$[#v7Q7E/x m_dot_l=0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*/w7c%B+
5、(Z fabs(T_SAT-C_T(cell,mix_th))/T_SAT; u'R)g0?5['j2
6、Z7'z/{ dS[eqn]=0.;+]"^)g"B,G+~ };u#O/P*U"E5v%F#L ;_0y0P7?+K R'V*y)v;C returnm_dot_l;7m&^#u e1S5o}$Y(A.A+P%{!E)^3[DEFINE_SOURCE(vap_src,cell,sec_th,dS,eqn)0m;`8c0m0i,[{/Y(x:U(V%T2D%j9b Thread*mix_th,*pri_th;.h*`2U1}&z9L4?%S5Z#n realm_dot_v;"]-`/~%_$X0r0r)P8W1i+d*? mix_
7、th=THREAD_SUPER_THREAD(sec_th);0q/b$(S4^"{)
8、!r m pri_th=THREAD_SUB_THREAD(mix_th,0);%S#h7F%{4@1X4K if(C_T(cell,mix_th)>=T_SAT){-@$H;x.?0?"h(x-]+u!u m_dot_v=0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*)
9、"O0V0x4F2O9i4W+c fabs(C_T(cell,mix_th)-T_SAT)/T_SAT;,t'w%r)n8[+J(U(v
10、._ dS[eqn]=0.;0x8q3l#g8D&V/+W*a } e C2Z6X.H8P+A#o ? else{"_,~.F$~/~1A:D m_dot_v=-0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*(j.}0}4Z/[0i'C5D fabs(T_SAT-C_T(cell,mix_th))/T_SAT;+V6J.z$~+f2b/g,E dS[eqn]=-0.1*C_R(cell,sec_th)*%p6} N _8b2Z1h&P fa
11、bs(C_T(cell,sec_th)-T_SAT)/T_SAT;1_%o!?2?8
12、*T })e%`&Z+h!v7w"H"L $L"K;+c)a:I4T/{ returnm_dot_v;.E.?3y p7T"U}#_(w0T2e/u:l8cDEFINE_SOURCE(enrg_src,cell,mix_th,dS,eqn)"?*a!R*V A5k8B{&]7l!q"e/u)v0m(K3p%p)_&Q Thread*pri_th,*sec_th;#~#k'Z+V v)x/N realm_dot;;w1}/D!M;u-D pri_t
13、h=THREAD_SUB_THREAD(mix_th,0);9
14、#H0p2v$S