| 查看: 1246 | 回复: 3 | |||
[求助]
UDF近壁面cell源项编写 已有1人参与
|
|
各位大神，目前我想向所有贴近壁面第一层cell中添加能量源项，代码如下，但是加载后运行结果与不加载是一样的，大神能不能看看问题出现在了哪里？ 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)*1000000; dS[eqn]=0; //Enrg_Source=2; Message("The value of Enrg_Source is: %f\n", Enrg_Source); } else { Enrg_Source=0; dS[eqn]=0; } return Enrg_Source; } end_f_loop(f, t) } |
回复此楼» 收录本帖的淘帖专辑推荐
 博士科研新路历程 |
» 猜你喜欢
211本科材料化工求调剂
已经有4人回复
-
272分材料子求调剂
已经有43人回复
-
270求调剂
已经有7人回复
-
271求调剂
已经有13人回复
-
283求调剂 086004考英二数二
已经有16人回复
-
085410 273分调剂
已经有4人回复
-
一志愿厦大0856，306求调剂
已经有3人回复
-
304求调剂（085602，过四级，一志愿985）
已经有28人回复
-
调剂
已经有4人回复
-
070300化学279求调剂
已经有19人回复
2楼2017-01-17 22:02:42
xingfuww
专家顾问 (正式写手)
-
专家经验: +78 - 仿真EPI: 2
- 应助: 46 (小学生)
- 金币: 6631.4
- 散金: 52
- 红花: 29
- 帖子: 542
- 在线: 284小时
- 虫号: 1776059
- 注册: 2012-04-24
- 性别: GG
- 专业: 化学工程及工业化学
- 管辖: 仿真模拟
3楼2017-01-18 13:28:56
4楼2017-01-18 15:55:36
