| 查看: 3552 | 回复: 17 | ||
[求助]
好心人帮帮忙吧，为什么这个UDF导入fluent以后，在边界条件处没有显示~~~ 已有1人参与
|
||
|
这是一个关于膜渗透的过程，利用UDF进行质量源项的添加。具体的程序贴出来了，问题是导入fluent以后，DEFINE_ADJUST定义的宏表现出来了，可是DEFINE_SOURCE 定义的函数在边界条件那里没有，我在fluent中设定了memory是2个，其他的就没有再改，是不是漏掉什么了？大家帮我看看吧~~ #include"udf.h" #include"mem.h" #define MW_N2 28.0 #define MW_C2H4 28 #define Q_N2 0 #define Q_C2H4 0.0000002435 #define P_outer 2160000 #define P_inner 120000 #define feed_domain 17 #define perment_domain 18 real molefrac( real frac, real molarone, real molartwo) { real a = frac / molarone ; real b = (1-frac)/molartwo ; return a/(a+b) ; } DEFINE_ADJUST(filter_adjust, domain) { cell_t c1, c0; Thread* wall_thread; Thread* wall_shadow_thread ; Thread* t0; Thread*t1; face_t wall = 27; face_t wall_shadow = 2; real A[ND_ND] ; real x_C2H4_outer, x_N2_outer, x_C2H4_inner, x_N2_inner, At, flux_C2H4, flux_N2, tot_flux ; int memb_id = 17 ; wall_thread= Lookup_Thread(domain,wall); tot_flux = 0 ; begin_f_loop (wall, wall_thread) { c0 = F_C0(wall,wall_thread); wall_shadow_thread = THREAD_SHADOW(wall_thread); wall_shadow = F_SHADOW(wall,wall_thread); c1 = F_C0(wall_shadow,wall_shadow_thread); t0 = F_C0_THREAD(wall,wall_thread); t1 =THREAD_T0(wall_shadow_thread); x_C2H4_outer = molefrac(C_YI(c0, t0, 0), MW_C2H4, MW_N2); x_N2_outer = 1.0 - x_C2H4_outer; x_C2H4_inner = molefrac(C_YI(c1, t1, 0), MW_C2H4, MW_N2); x_N2_inner = 1.0 - x_C2H4_inner; F_AREA(A,wall, wall_thread); At = NV_MAG(A); flux_C2H4 = Q_C2H4 * ( x_C2H4_outer * P_outer - x_C2H4_inner * P_inner); flux_N2 = Q_N2 * ( x_N2_outer * P_outer - x_N2_inner * P_inner); tot_flux = tot_flux + MW_C2H4*flux_C2H4*At + MW_N2*flux_N2*At; C_UDMI(c0, t0, 0) = -(MW_C2H4*flux_C2H4)*At/C_VOLUME(c0,t0); C_UDMI(c1, t1, 1) = (MW_C2H4*flux_C2H4)*At/C_VOLUME(c1,t1); C_UDMI(c0, t0, 2) = -(MW_N2*flux_N2)*At/C_VOLUME(c0,t0); C_UDMI(c1, t1, 3) = (MW_N2*flux_N2)*At/C_VOLUME(c1,t1); } end_f_loop (wall,wall_thread) } DEFINE_SOURCE(gas1, c, t, dS, eqn) { Thread* tf; int n ; real source=0.0; Domain *domain = Get_Domain(feed_domain); face_t wall = 27; Thread* wall_thread=Lookup_Thread(domain,wall); Thread* wall_shadow_thread = THREAD_SHADOW(wall_thread); c_face_loop( c, t, n) {tf = C_FACE_THREAD(c,t,n); if(tf==wall_thread) return source=C_UDMI(c, t, 0); } return source; } DEFINE_SOURCE(gas2, c, t, dS, eqn) { Thread* tf; int n ; real source=0.0; Domain *domain = Get_Domain(perment_domain); face_t wall = 27; Thread* wall_thread=Lookup_Thread(domain,wall); Thread* wall_shadow_thread = THREAD_SHADOW(wall_thread); c_face_loop( c, t, n) {tf = C_FACE_THREAD(c,t,n); if(tf==wall_shadow_thread) return source=C_UDMI(c,t,1); } return source; } 有点长，可是很简单。求求大家帮我看看吧 |
回复此楼» 猜你喜欢
333求调剂
已经有10人回复
-
307求调剂
已经有3人回复
-
考研化学308分求调剂
已经有4人回复
-
298求调剂
已经有7人回复
-
263求调剂
已经有5人回复
-
276求调剂。有半年电池和半年高分子实习经历
已经有6人回复
-
求调剂材料学硕080500，总分289分
已经有21人回复
-
北科281学硕材料求调剂
已经有8人回复
-
08工学调剂
已经有9人回复
-
291 求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 新手求救啊！Fluent非稳态计算 怎么样才能在每一个时间步上调用一次UDF啊？
已经有15人回复
- fluent中用udf编写三维速度边界条件出现错误，大家来指导一下啊
已经有9人回复
- 用fluent模拟消声器的入口边界条件为什么是速度入口条件？
已经有10人回复
- Fluent 压力出口边界条件设置
已经有6人回复
- Fluent中压力入口边界条件怎么设置
已经有21人回复
- fluent中定义两个出口边界条件
已经有10人回复
- FLUENT源项UDF非稳态迭代的一个疑问请教
已经有8人回复
- 关于FLUENT中共轭传热中的couple边界条件问题
已经有12人回复
- fluent中壁面热边界条件的问题
已经有16人回复
- 有关fluent的udf中的hook
已经有6人回复
- 【讨论】Fluent边界条件及虚拟质量力
已经有6人回复
- 【求助】FLUENT中怎么同时定义压力和流量边界条件啊
已经有18人回复
- 【求助】fluent计算中，udf编译通过，初始化出错
已经有11人回复
4楼2012-05-18 07:28:00
030301029
捐助贵宾 (正式写手)
- 应助: 76 (初中生)
- 金币: 2203.4
- 散金: 306
- 红花: 8
- 帖子: 890
- 在线: 122.3小时
- 虫号: 563419
- 注册: 2008-05-25
- 专业: 化学反应工程
9楼2012-05-22 09:24:26
2楼2012-05-18 01:37:51
3楼2012-05-18 07:10:00
5楼2012-05-18 09:29:55
6楼2012-05-18 10:28:37
7楼2012-05-21 19:28:41
030301029
捐助贵宾 (正式写手)
- 应助: 76 (初中生)
- 金币: 2203.4
- 散金: 306
- 红花: 8
- 帖子: 890
- 在线: 122.3小时
- 虫号: 563419
- 注册: 2008-05-25
- 专业: 化学反应工程
8楼2012-05-22 09:22:44
030301029
捐助贵宾 (正式写手)
- 应助: 76 (初中生)
- 金币: 2203.4
- 散金: 306
- 红花: 8
- 帖子: 890
- 在线: 122.3小时
- 虫号: 563419
- 注册: 2008-05-25
- 专业: 化学反应工程
10楼2012-05-22 10:02:42
