| 查看: 1197 | 回复: 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 ] |
回复此楼» 猜你喜欢
时间戳变了，能看出什么问题？
已经有20人回复
-
重要消息，中午系统在维护
已经有12人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有20人回复
-
言之凿凿，8月21号(明天)放榜
已经有14人回复
-
只有每年这种时候来逛逛小木虫
已经有21人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有125人回复
-
今天基金会出结果吗？20260819
已经有15人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
longma135
金虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 1186.6
- 散金: 15
- 帖子: 199
- 在线: 207小时
- 虫号: 2533271
- 注册: 2013-07-05
- 性别: GG
- 专业: 能源化工
4楼2014-11-21 10:04:04
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
lianghtbuct
木虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 4681.7
- 红花: 1
- 帖子: 211
- 在线: 67.4小时
- 虫号: 2419564
- 注册: 2013-04-15
- 性别: GG
- 专业: 化学反应工程
5楼2014-11-21 10:18:29
