应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 用UDF编写化学反应速率模拟结果不合实际
111/2返回列表12下一页
查看: 2368  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

BBKtime

新虫 (小有名气)

[求助] 用UDF编写化学反应速率模拟结果不合实际 已有3人参与

目前由于刚开始使用UDF，参照fluent帮助文档内的reaction rates中的实例把其中的程序复制粘贴后通过function hooks命令load到程序当中，除管道尺寸外，其它参数均参照文档中进行设置（由于刚开始学习，主要是自己瞎倒腾，一些设置自己也不清楚）。
 \"用UDF编写化学反应速率模拟结果不合实际\"
模拟的问题见上图，管道底部和右端为多孔介质区，化学反应也发生在该区域，管道其余部分不发生反应。其中UDF的代码为
#include "udf.h"
#define K1 2.0e-2
#define K2 5.
DEFINE_VR_RATE(vol_reac_rate,c,t,r,mole_weight,species_mf,rate,rr_t)
{
    real s1=species_mf[0];
    real mw1=mole_weight[0];
    if(FLUID_THREAD_P(t) && THREAD_VAR(t).fluid.porous)
      *rate=K1*s1/pow((1.+K2*s1),2.0)/mw1;
    else
      *rate=0.;
    *rr_t=*rate;
}
fluent帮助文件中给出的反应物的质量分数图为
 \"用UDF编写化学反应速率模拟结果不合实际-1\"

而我模拟出的结果图为
 \"用UDF编写化学反应速率模拟结果不合实际-2\"
最终得到的模拟结果与感觉与实际情况相差太远，有没有好心人给看看都哪方面出了问题，在此提前拜谢了！一个人默默地钻研真是心累，希望能得到大家的帮助吧
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : porous-reaction.cas
    • 2016-10-21 21:07:34, 467.86 K

» 猜你喜欢
1楼 2016-10-21 21:10:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

BBKtime

新虫 (小有名气)

用UDF编写化学反应速率模拟结果不合实际-3
1.jpg
赞一下 回复此楼
2楼2016-10-21 21:18:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

BBKtime

新虫 (小有名气)
其中第一张图是帮助文档中给出的模拟结果，第二张是自己模拟出的不切实际的结果
用UDF编写化学反应速率模拟结果不合实际-4
3.jpg


用UDF编写化学反应速率模拟结果不合实际-5
2.jpg
赞一下 回复此楼
3楼2016-10-21 21:20:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Xaiver

新虫 (初入文坛)

【答案】应助回帖

感谢参与，应助指数 +1
速度场就不对吧，流体在催化剂层流动没有阻力，整个管子就是简单的层流。
而且，定义的化学反应，为什么计量系数和速率指数都是0。
赞一下 回复此楼
4楼2016-10-22 21:57:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

BBKtime

新虫 (小有名气)
引用回帖:
4楼: Originally posted by Xaiver at 2016-10-22 21:57:14
速度场就不对吧，流体在催化剂层流动没有阻力，整个管子就是简单的层流。
而且，定义的化学反应，为什么计量系数和速率指数都是0。

非常感激您的回复，我不太清楚用udf给定了化学反应速率之后计量系数和速率指数还需要设置么？按照fluent帮助文档中的那个例子设置了多孔介质区的惯性阻力系数，但并没有得到相应的结果。不知道是哪方面设置的问题，还望您能给看下，拜托了

发自小木虫IOS客户端
赞一下 回复此楼
5楼2016-10-23 10:50:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dai26520

专家顾问 (著名写手)
当时看过这个算例，不过没有详细做过，按着设置做，想来不会差别这么大

发自小木虫Android客户端
赞一下 回复此楼
个性只是一个代号，做好自己就行
6楼2016-10-23 18:15:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

BBKtime

新虫 (小有名气)
引用回帖:
6楼: Originally posted by dai26520 at 2016-10-23 18:15:58
当时看过这个算例，不过没有详细做过，按着设置做，想来不会差别这么大

这个帮助文件中给出的设置信息较少，算例中给出的信息都在fluent中进行了设置，其它的设置按照自己的理解，不知道是哪方面出了问题，还望您能够指点迷津
赞一下 回复此楼
7楼2016-10-23 22:44:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

星星在我家门

木虫 (著名写手)

热爱互联网
我也看到了这个，正打算做做，楼主是做哪方面反应的呢？
赞一下 回复此楼
于无声听惊雷
8楼2017-05-12 22:03:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

banallove

新虫 (初入文坛)
楼主这个问题解决了没有呢，您的case文件首先没有对多孔介质单独分区吧，这样就没法设置阻力系数。我也在做这个，就是反应的问题不太明白
赞一下 回复此楼
9楼2017-08-01 15:05:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huwanpeng2

金虫 (职业作家)

【答案】应助回帖

看到COMSOL软件写一个表达式即可

发自小木虫Android客户端
赞一下 回复此楼
skloche.ecust
10楼2017-08-01 23:57:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 BBKtime 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭