| 查看: 2019 | 回复: 4 | ||||
[交流]
UDF源项编写中如何调用壁面温度
|
|
楼主新人，目前想写一个液滴的能量源项，但是在源项中不知道如何调用壁面温度，代码如下，可以编译，但无法运行。 DEFINE_SOURCE(liq_enrg_src, c, sec_th, dS, eqn) { Thread *mix_th, *pri_th, *t; Domain *domain=Get_Domain(1); int zone_ID=18;/*wall ID*/ real m_dot_d=0; real kd=0; real Cnt=0; real a; real b; real d=0; real vg=0; real Reg=0; real tr=0; real Enrg_Source=0; face_t f; mix_th = THREAD_SUPER_THREAD(sec_th); pri_th = THREAD_SUB_THREAD(mix_th,0); t=Lookup_Thread(domain,zone_ID); begin_f_loop(f, t) { Cnt=(1.0-C_VOF(c,pri_th))*C_R(c,sec_th); a=Cnt/C_R(c,pri_th); if(a<0.3) { kd=0.18*sqrt(0.026/(C_R(c,pri_th)*0.008)); } else { kd=pow(a,-0.85)*sqrt(0.026/(C_R(c,pri_th)*0.008)); } m_dot_d=kd*Cnt; dS[eqn]=0; vg=sqrt(pow(C_U(c,pri_th),2)+pow(C_V(c,pri_th),2)+pow(C_W(c,pri_th),2)); Reg=C_R(c,pri_th)*vg*0.008/C_MU_EFF(c,pri_th); d=0.00796*0.026/C_R(c,pri_th)/pow(C_VOF(c,pri_th)*vg,2)*pow(Reg,2/3)*pow(C_R(c,pri_th)/C_R(c,sec_th),-1/3)*pow(C_MU_EFF(c,pri_th)/C_MU_EFF(c,sec_th),2/3); tr=3.14*sqrt(C_R(c,sec_th)*d*d*d/16/0.026); b=pow(18*pow(C_K_L(c,pri_th)*tr,3)*C_R(c,pri_th)*LAT_HT*pow(m_dot_d,5)/pow(C_R(c,sec_th),4)/C_MU_EFF(c,pri_th)/(1.0-C_VOF(c,pri_th))/(F_T(f,t)-T_SAT),0.25); Enrg_Source=(F_T(f,t)-T_SAT)*b*4/0.008; if(C_T(c,mix_th)<T_SAT) { return 0; } else { return Enrg_Source; } } |
回复此楼» 收录本帖的淘帖专辑推荐
 博士科研新路历程 |
» 猜你喜欢
302分求调剂 一志愿安徽大学085601
已经有5人回复
-
材料与化工371求调剂
已经有4人回复
-
313求调剂
已经有3人回复
-
求调剂
已经有6人回复
-
材料工程302分求调剂
已经有5人回复
-
生物学308分求调剂（一志愿华东师大）
已经有4人回复
-
304求调剂（085602，过四级，一志愿985）
已经有8人回复
-
复试调剂
已经有6人回复
-
求
已经有4人回复
-
求调剂
已经有8人回复
2楼2016-12-27 16:16:32
3楼2016-12-27 23:44:14
4楼2016-12-28 10:46:57
|
这个是我重新写的，不知是否已经实现我所说的功能，另我这个源项加进去是要模拟液滴与壁面之间的对流换热，所以我的想法是给壁面液体设置一个正能量源项，给气体设置一个负的能量源项。 DEFINE_SOURCE(liq_enrg_src, c, sec_th, dS, eqn) { Thread *mix_th, *pri_th, *t; Domain *domain=Get_Domain(1); int zone_ID=18;/*wall ID*/ real m_dot_d=0; real kd; real Cnt; real a; real Area[ND_ND]; real b; real d; real vg; real Reg; real tr; real Enrg_Source=0; face_t f; mix_th = THREAD_SUPER_THREAD(sec_th); pri_th = THREAD_SUB_THREAD(mix_th,0); t=Lookup_Thread(domain,zone_ID); begin_f_loop(f, t) { if(c == F_C0(f,t)) { Cnt=(1.0-C_VOF(c,pri_th))*C_R(c,sec_th); a=Cnt/C_R(c,pri_th); if(a<0.3) { kd=0.18*sqrt(0.026/(C_R(c,pri_th)*0.008)); } else { kd=pow(a,-0.85)*sqrt(0.026/(C_R(c,pri_th)*0.008)); } m_dot_d=kd*Cnt; vg=sqrt(pow(C_U(c,pri_th),2)+pow(C_V(c,pri_th),2)+pow(C_W(c,pri_th),2)); Reg=C_R(c,pri_th)*vg*0.008/C_MU_L(c,pri_th); d=0.00796*0.026/C_R(c,pri_th)/pow(C_VOF(c,pri_th)*vg,2)*pow(Reg,2/3)*pow(C_R(c,pri_th)/C_R(c,sec_th),-1/3)*pow(C_MU_L(c,pri_th)/C_MU_L(c,sec_th),2/3); tr=3.14*sqrt(C_R(c,sec_th)*d*d*d/16/0.026); b=pow(18*pow(C_K_L(c,pri_th)*tr,3)*C_R(c,pri_th)*LAT_HT*pow(m_dot_d,5)/pow(C_R(c,sec_th),4)/C_MU_EFF(c,pri_th)/(1.0-C_VOF(c,pri_th))/(F_T(f,t)-T_SAT),0.25); F_AREA(Area,f,t); Enrg_Source=(F_T(f,t)-T_SAT)*b/C_VOLUME(c,mix_th)/(1-C_VOF(c, pri_th))*NV_MAG(Area); dS[eqn]=0; //Enrg_Source=2; return Enrg_Source; } } end_f_loop(f, t) } |
5楼2016-12-28 13:28:51
