首页
发现
课程
培训
文章
案例
问答
需求
服务
行家
赛事
热门搜索
发布
消息
注册
|
登录
首页
/
文章
/
详情
反应速率的UDF实现与软件自带结果对比(一)
余花生
7月前
浏览9938
关注
正文共:
1606
字
11
图
预计阅读时间:
5
分钟
1 前言
前期推送了多个关于燃烧和化学反应的案例,从应用范围来讲,通用有限速率模型适应更多的应用场景。比如在氢能行业,关于燃料电池电化学的模拟、关于重整制氢的模拟等等,就需要考虑详细的化学反应。因此,我们觉得有必要对详细化学反应的模拟做一些更深入的探讨,接下来我们以一个具体的例子,先来看看FLUENT自身模型求解详细化学反应的关键内容,然后用UDF来模拟化学反应并与软件自身计算结果进行对比,以便获得更深刻认识,今天我们先进行第一个环节。
2 反应案例
我们考虑以下的化学反应,这是甲烷蒸汽重整制氢的反应之一,反应焓等于
摩尔加权
的反应物和生成物的生成焓之差,本反应的反应焓为正直,说明是一个吸热反应。
FLUENT材料物性的标准状态焓就是材料的生成焓,单位为J/kmol,因此本反应的反应焓转化为FLUENT的反应焓就是206.3×10
6
J/kmol。这个反应的字面理解就是每消耗1kmol的CH4或者每消耗1kmol的H2O或者每生成1kmol的CO或者每生成3kmol的H2,需要吸收206.3×10
6
J的能量。当然,也可以理解为每消耗1/3kmol的CH4或者每消耗1/3kmol的H2O或者每生成1/3kmol的CO或者每生成1kmol的H2,需要吸收1/3×206.3×10
6
J的能量,也即是前文所述的反应焓是摩尔加权计算结果,注意数值上一定要对应。通用有限速率模型中,FLUENT通过在能量方程中增加由于化学反应产生的源项来计算反应热,式5.10中各项的单位为
S
h,rxn
:W/m3,
h
j
:J/kmol、M
j
:kg/kmol、
R
j
:kg/m
3
s。
该式也等效于
ΔH
(J/kmol)×
R
(kmol/m
3
s)。
我们假定反应速率如下,反应速率式的各个变量、各个指数的单位与FLUENT的单位吻合。另外,根据反应式,我们可以得出各个组分的摩尔速率关系为:
r
CH4
=
r
H2O
=
r
CO
=
r
I
,
r
H2
=
3r
I
。也就是各组分
摩尔速率
之比就等于化学计量数之比,这是化学反应基本常识。
在化学反应实验中,通常是得出某个物质的反应速率,也就得到了整个反应以及其他物质的反应速率。FLUENT通过反应摩尔速率式乘以各反应组分的化学计量数,等到该组分的反应摩尔速率,注意反应物的速率为负(表示消耗)、生成物的速率为正(表示生成)。
根据上面的反应速率式,在FLUENT中输入对应的参数。
我们设定反应入口的组分为H2O和CH4,并且摩尔比重为0.6和0.4。其他的相关设置我们就不再细说了,读者朋友可以参考以前的案例。这里需要指出的是,为了能够顺利自启动反应,我们给反应入口设定一个高的温度773.15K。
接下来我们对计算结果进行平衡检查。
首先,检查一下质量守恒情况,进出口满足质量守恒,这是最基础的,没有这个平衡,整个模拟计算可执行一票否决。
然后,我们看一下进出口的各组分的质量分数,进出口的质量分数差值乘以上述总的质量流量就得到各组分在本反应中的质量消耗量(反应物)和质量生成量(生成物),单位为kg/s。再除以各组分的摩尔质量(kg/kmol),就得到了摩尔消耗量和摩尔生成量,单位为kmol/s。各物质的摩尔消耗量和摩尔生成量分别为(CH4):1.1884E-08kmol/s、(H2O):1.1893E-08kmol/s、(CO):1.1885E-08kmol/s和(H2):3.5725E-08kmol/s,满足化学反应计量数1:1:1:3的关系。
最后,检查一下能量守恒情况。根据以上的反应速率,乘以反应焓就可以得到该反应吸收的总热流,比如用CH4的反应速率计算能量如下,结果为2.452W。我们读取FLUENT的能量报表,可以看出反应热为-2.450W,与手动计算值一样(负号表示吸热)。
反应吸热通过壁面来补充,因此进出口的温度基本不变,但是这不代表局部温度不变,化学反应剧烈的地方由于吸热作用,温度一定会降低(见以下温度云图)。
来源:仿真与工程
登录后免费查看全文
立即登录
Fluent
燃烧
化学
燃料电池
UDF
通用
材料
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-01-16
最近编辑:7月前
余花生
签名征集中
关注
获赞 176
粉丝 312
文章 295
课程 0
点赞
收藏
0/200
清空
提交
1条评论
绣春刀
签名征集中
7月前
咨询一下大佬,我想做一个氯化铈雾化热解反应,该反应化学方程式自己设置的组分输运模型,但发生不了反应
回复
课程
培训
服务
行家
总布置工程师入门第四阶段:新手小白需要掌握的office办公软件PPT、Word、Excel
表面化学反应模拟
总布置工程师入门第二阶段:整车开发流程及我们在其中扮演的角色
反应堆内粒子流动
相关推荐
结构及流体仿真杂谈:UDF、燃烧、耦合项目及案例讲解
汽车外气动计算和动画制作
基于FLUENT/UDF编写圆柱电芯各向异性导热系数
Fluent及CFD-Post后处理专题:获得仿真后处理常用技巧和学术报告科研数据展示能力
Fluent燃烧与反应流专题9讲:获得燃烧及化学反应流仿真流程及思路
最新文章
趋势ℱ制造产业升级 碳纤维行业迎来高速发展机遇
技术ℱ当碳纤维邂逅3D打印技术
复材ℱ三维编织复合材料研究现状及应用
应用ℱ碳纤维加固是怎么施工的?这里有个实例
『专题』日本碳纤维的研发现状
热门文章
ABAQUS中Cohesive粘聚力模型的2种定义方式(附案例操作步骤)
几种常见的热仿真软件
STAR CCM 案例|电池包散热
盘点·近十年来国外各公司推出的碳纤维产品
工业互联网中数字孪生系统的机理+数据融合建模方法
其他人都在看
某车型中高速工况热制动车身抖动分析与控制
Abaqus分析常见问题及解决方法(2):零主元和过约束
Abaqus分析常见问题及解决方法(3):负特征值(Negative Eigenvalue)
分不清楚中国动车组型号?看完本文秒懂(还有技术参数哦)
电磁场仿真 | ChatGPT请回答,我想和你聊聊
VIP会员
学习
福利任务
兑换礼品
下载APP
联系我们
微信客服
联系客服
人工服务时间为周一至周五的9:30-19:30
非工作时间请在微信客服留言
客服热线:
4000-969-010
邮箱:
service@fangzhenxiu.com
地址:
北京市朝阳区莱锦创意园CN08座
帮助与反馈
返回顶部