| 查看: 1017 | 回复: 6 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
longma135金虫 (小有名气)
|
[求助]
fluent中添加氧气源项，编译能通过，一计算就出错，附程序，求大牛解答
|
|
|
UDM已经定义，因为程序较多，只贴了氧气的源项，因为其他的源项加入的时候没有错误。 DEFINE_SOURCE(source_o2_g,c,t,dS,eqn) /*氧气的质量源项，定义在气相中*/ { Thread *mixture_t=THREAD_SUPER_THREAD(t); /*混合相的指针*/ Thread *ts=THREAD_SUB_THREAD(mixture_t,1); real Ro13,Ro14,Ro15,Ro16,Ro17,Ro18,Ro19; int n; face_t f; Thread *tf; Thread *t1; cell_t c1; Ro17=0.0; if(C_UDSI(c,ts,4)>1170.0) Ro13=1170.0; else if(C_UDSI(c,ts,4)<730.0) Ro13=730.0; else Ro13=C_UDSI(c,ts,4); Ro14=2500.0*exp(-6420.0/Ro13); Ro15=(Ro14+2.0)/(2.0+2.0*Ro14); /*氧气的化学反应当量系数*/ Ro16=-C_UDMI(c,mixture_t,8)*16.0-C_UDMI(c,mixture_t,9)*16.0-C_UDMI(c,mixture_t,10)*16.0-C_UDMI(c,mixture_t,7)/12.0*Ro15*32.0; c_face_loop(c,t,n) { f=C_FACE(c,t,n); tf=C_FACE_THREAD(c,t,n); if(F_FLUX(f,tf)>0.0) { Ro17+=F_FLUX(f,tf)*C_YI(c,t,5); } else { t1=F_C1_THREAD(f,tf); c1=F_C1(f,tf); Ro17+=-F_FLUX(f,tf)*C_YI(c,t,5); } } Ro18=Ro16+Ro17/C_VOLUME(c,t); return Ro18; dS[eqn]=0; } 添加源项计算后就出现这种问题  [ Last edited by longma135 on 2014-3-13 at 16:12 ] |
回复此楼» 猜你喜欢
0857调剂
已经有3人回复
-
材料学调剂
已经有11人回复
-
材料类求调剂
已经有11人回复
-
高分子化学与物理调剂
已经有11人回复
-
26申博
已经有6人回复
-
328求调剂
已经有5人回复
-
化工专硕348，一志愿985求调剂
已经有8人回复
-
0856求调剂285
已经有8人回复
-
290求调剂
已经有11人回复
-
281求调剂
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
longma135
金虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 1186.6
- 散金: 15
- 帖子: 199
- 在线: 207小时
- 虫号: 2533271
- 注册: 2013-07-05
- 性别: GG
- 专业: 能源化工
6楼2014-11-21 10:56:27
longma135
金虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 1186.6
- 散金: 15
- 帖子: 199
- 在线: 207小时
- 虫号: 2533271
- 注册: 2013-07-05
- 性别: GG
- 专业: 能源化工
2楼2014-03-15 18:14:19
lianghtbuct
木虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 4681.7
- 红花: 1
- 帖子: 211
- 在线: 67.4小时
- 虫号: 2419564
- 注册: 2013-04-15
- 性别: GG
- 专业: 化学反应工程
3楼2014-11-20 21:39:34
longma135
金虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 1186.6
- 散金: 15
- 帖子: 199
- 在线: 207小时
- 虫号: 2533271
- 注册: 2013-07-05
- 性别: GG
- 专业: 能源化工
4楼2014-11-21 10:04:04
10