首页
发现
课程
培训
文章
案例
问答
需求
服务
行家
赛事
热门搜索
发布
消息
注册
|
登录
首页
/
文章
/
详情
反应速率的UDF实现与软件自带结果对比(一)
余花生
10月前
浏览10967
关注
正文共:
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
最近编辑:10月前
余花生
签名征集中
关注
获赞 195
粉丝 372
文章 306
课程 0
点赞
收藏
0/200
清空
提交
1条评论
绣春刀
签名征集中
10月前
咨询一下大佬,我想做一个氯化铈雾化热解反应,该反应化学方程式自己设置的组分输运模型,但发生不了反应
回复
课程
培训
服务
行家
2023仿真知识周(二):学犀牛和Grasshopper年薪30到60万是如何练成的
LS-DYNA偏心不耦合的临自由面岩石双孔微差爆破(JH-2本构)
有问必答第2期-无主题自由问答专场【直播回放】
LS-DYNA的PBM-SPH-FEM耦合模拟边坡爆破岩块抛掷
相关推荐
跟着帮助文档学STAR-CCM 快速入门进阶与精通180讲-全功能模块案例合集
ansys fluent表达式功能实现二维表格插值(以锂电池充电map为例)
【Fluent 2023R1系列之八】Fluent组分输运 化学反应 燃烧
Fluent燃烧与反应流专题9讲:获得燃烧及化学反应流仿真流程及思路
结构及流体仿真杂谈:UDF、燃烧、耦合项目及案例讲解
最新文章
一款关于人体建模的仿真软件:AnyBody
综述 | 设备状态监测中处理工业数据分布不平衡的重采样技术(下)
应力、应变、位移、变形
储层裂缝研究方法
射频类公司汇总(苏州篇)
热门文章
盘点·近十年来国外各公司推出的碳纤维产品
几种常见的热仿真软件
ABAQUS中Cohesive粘聚力模型的2种定义方式(附案例操作步骤)
仿真工作者必须知道的15款开源软件!
“卡脖子”的国产工业软件,现状如何?路又在何方?
其他人都在看
Abaqus分析常见问题及解决方法(2):零主元和过约束
STAR CCM 案例|电池包散热
Abaqus分析常见问题及解决方法(3):负特征值(Negative Eigenvalue)
电磁场仿真 | ChatGPT请回答,我想和你聊聊
仿真笔记——ANSYS APDL命令汇总(收藏备用)
VIP会员
学习
福利任务
兑换礼品
下载APP
联系我们
微信客服
联系客服
人工服务时间为周一至周五的9:30-19:30
非工作时间请在微信客服留言
客服热线:
4000-969-010
邮箱:
service@fangzhenxiu.com
地址:
北京市朝阳区莱锦创意园CN08座
帮助与反馈
返回顶部