| 查看: 2065 | 回复: 4 | |||
| 【悬赏金币】回答本帖问题，作者cfd菜鸡将赠送您 40 个金币 | |||
[求助]
这是我编写湿空气冷凝的udf导入fluent后初始化就无法完成，是不是死循环，求大神指导
|
|||
|
#include "udf.h" #include "sg.h" #include "sg_mphase.h" #include "flow.h" #include "mem.h" #define LAT_HT 2.257e6 #define WALL_NUMBER 3 #define scdomain_index 3 #define prim_index 0 #define index_evap_primary 1 int phase_domain_index; real face_center[ND_ND], cell_center[ND_ND], a[ND_ND], distance; real p_op, p, diffusivity; real cell_vapor_pressure, wall_temper, wall_pressure, wall_mass_fraction; real NV_VEC(A); real area, area_density, mass_transfer_coeff; real water_temper, water_satpressure, water_mass_fraction; real vapor_density; int n; Domain *subdomain; cell_t c; Thread *t; Thread **pt; Thread *tp; Thread *ts; Thread *tf, *tb; face_t fa; face_t fb; real source; DEFINE_SOURCE(vap_src, cell, first, dS, eqn) { Thread *mixer, *sec_th; real m_dot_first; real mass_dot; mixer = THREAD_SUPER_THREAD(first); sec_th = THREAD_SUB_THREAD(mixer, 1); real Db; real water_mass_fiction; real area; real volume; begin_c_loop(cell,mixer) { volume = C_VOLUME(cell, mixer);//网格体积// c_face_loop(cell, mixer, n)//对所有cell面循环// { fa = C_FACE(cell, mixer, n);//获取面指针// tf = C_FACE_THREAD(cell, mixer, n);//获取面所在thread指针// if (THREAD_ID(tf) == WALL_NUMBER) //判断是否为壁面网格// { if (C_VOF(cell, sec_th) == 1.) //判断液体是否充满壁面// mass_dot = 0.0; //质量源为0// else { C_CENTROID(cell_center, cell, mixer);//获取cell中心坐标// F_CENTROID(face_center, fa, tf);//获取面中心坐标// NV_VV(a, =, cell_center, -, face_center);//体心到面心矢量// distance = NV_MAG(a);//体心到面心距离// F_AREA(A, fa, tf);//面积矢量// area = NV_MAG(A);//面积// wall_temper = F_T(fa, tf) - 273.15;//获取壁面温度摄氏度// Db = (4.709 + 0.196 * wall_temper + 0.006122 * wall_temper * wall_temper + 0.0004386 * wall_temper * wall_temper * wall_temper) / 1000;//壁面温度对应饱和含湿量// water_mass_fiction = C_YI(cell, first, 0);//水蒸气质量分数// if (water_mass_fiction > Db) { mass_dot = -C_R(cell, first) * 0.256 * 0.0001 * (water_mass_fiction - Db) / distance * area/volume;//质量源求解公式// dS[eqn] = 0; } } } } } end_c_loop(cell,mixer) return mass_dot; } DEFINE_SOURCE(liq_src, cell, second, dS, eqn) { Thread *mix_th, *fir_th; real mass_dot; real Db; real water_mass_fiction; real area; real volume; mix_th = THREAD_SUPER_THREAD(second); fir_th = THREAD_SUB_THREAD(mix_th, 1); begin_c_loop(cell, mix_th) { volume = C_VOLUME(cell, mix_th);//网格体积// c_face_loop(cell, mix_th, n) { fa = C_FACE(cell, mix_th, n); tf = C_FACE_THREAD(cell, mix_th, n); if (THREAD_ID(tf) == WALL_NUMBER) { if (C_VOF(cell, fir_th) == 1.) mass_dot = 0.0; else { C_CENTROID(cell_center, cell, mix_th); F_CENTROID(face_center, fa, tf); NV_VV(a, =, cell_center, -, face_center); distance = NV_MAG(a); F_AREA(A, fa, tf);//面积矢量// area = NV_MAG(A);//面积// wall_temper = F_T(fa, tf) - 273.15; Db = (4.709 + 0.196 * wall_temper + 0.006122 * wall_temper * wall_temper + 0.0004386 * wall_temper * wall_temper * wall_temper) / 1000; water_mass_fiction = C_YI(cell, second, 0); if (water_mass_fiction > Db) { mass_dot = C_R(cell, second) * 0.256 * 0.0001 * (water_mass_fiction - Db) / distance * area / volume; dS[eqn] = 0; } } } } } end_c_loop(cell, mixer) return mass_dot; } DEFINE_SOURCE(enrg_src, cell, mixer, dS, eqn) { Thread* first, * sec_th; real m_dot; real mass_dot; first = THREAD_SUB_THREAD(mixer, 0); sec_th = THREAD_SUB_THREAD(mixer, 1); real Db; real water_mass_fiction; real area; real volume; begin_c_loop(cell, mixer) { c_face_loop(cell, mixer, n) { fa = C_FACE(cell, mixer, n); tf = C_FACE_THREAD(cell, mixer, n); if (THREAD_ID(tf) == WALL_NUMBER) { if (C_VOF(cell, sec_th) == 1.) mass_dot = 0.0; else { C_CENTROID(cell_center, cell, mixer); F_CENTROID(face_center, fa, tf); NV_VV(a, =, cell_center, -, face_center); distance = NV_MAG(a); F_AREA(A, fa, tf);//面积矢量// area = NV_MAG(A);//面积// wall_temper = F_T(fa, tf) - 273.15; Db = (4.709 + 0.196 * wall_temper + 0.006122 * wall_temper * wall_temper + 0.0004386 * wall_temper * wall_temper * wall_temper) / 1000; water_mass_fiction = C_YI(cell, first, 0); if (water_mass_fiction > Db) { mass_dot = -C_R(cell, first) * 0.256 * 0.0001 * (water_mass_fiction - Db) / distance * area / volume; dS[eqn] = 0; } } } } } end_c_loop(cell, mixer) return mass_dot * LAT_HT; |
回复此楼» 猜你喜欢
085701环境工程专硕求调剂
已经有6人回复
-
0817化学工程319求调剂
已经有10人回复
-
伙伴们，祝我生日快乐吧
已经有5人回复
-
274求调剂0856材料化工
已经有10人回复
-
337一志愿华南理工0805材料求调剂
已经有6人回复
-
（081700）化学工程与技术-298分求调剂
已经有7人回复
-
材料与化工304求B区调剂
已经有5人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有9人回复
-
一志愿北京化工大学材料与化工296分求调剂
已经有17人回复
-
0703化学一志愿211 总分320求调剂
已经有3人回复
mopsite
木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 8854.7
- 散金: 424
- 红花: 2
- 帖子: 2579
- 在线: 615小时
- 虫号: 172296
- 注册: 2006-01-19
- 专业: 传热传质学
2楼2021-05-12 14:28:33
3楼2021-09-28 21:22:09
4楼2021-09-28 21:22:25
5楼2021-10-01 09:45:33
