首页
发现
课程
培训
文章
案例
问答
需求
服务
行家
赛事
热门搜索
发布
消息
注册
|
登录
首页
/
文章
/
详情
反应速率的UDF实现与软件自带结果对比(一)
余花生
1年前
浏览12483
关注
正文共:
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
最近编辑:1年前
余花生
签名征集中
关注
获赞 217
粉丝 450
文章 322
课程 0
点赞
收藏
0/200
清空
提交
1条评论
绣春刀
签名征集中
1年前
咨询一下大佬,我想做一个氯化铈雾化热解反应,该反应化学方程式自己设置的组分输运模型,但发生不了反应
回复
课程
培训
服务
行家
有问必答第2期-无主题自由问答专场【直播回放】
大学本科课程 现代通信系统概论(上)
金砖大赛工程仿真执行仲裁长赵敏专家:工业软件定义金砖工业未来
ALIAS汽车设计职业培训班-汽车A面从入门到精通(提高班-持续加餐更新)
相关推荐
跟着帮助文档学STAR-CCM 快速入门进阶与精通180讲-全功能模块案例合集
Fluent燃烧与反应流专题9讲:获得燃烧及化学反应流仿真流程及思路
基于Workbench平台Fluent流固热多物理场仿真计算
【Fluent 2023R1系列之八】Fluent组分输运 化学反应 燃烧
最新文章
【机械设计】机械设计的30个机械结构功能,个个眼前一亮
DCDC开关电源中的EMI的原因分析……
Inspire与SimSolid如何“零妥协”重塑增材制造设计
晶体塑性有限元 Abaqus 三维泰森多边形(voronoi模型)插件 V5.0
晶体塑性有限元 Abaqus 三维泰森多边形(voronoi模型)插件 V6.0
热门文章
通用/泛亚/蔚来/弗迪/上海电驱动等百位新能源专家确认发言,易贸三电活动6月苏州见
锂离子电池生产现场异物管控
仿真秀需求交易模块介绍
2024仿真秀618学习狂欢课程特惠大放送
DSP逐波限流(CBC)方案
其他人都在看
2024第四届中国算力之都·宁夏 工程仿真大赛CAE网格赛道【赛事通知】
故障诊断实验台 | BTS100轴承寿命预测故障实验台(16.8w)
固态电池碾压油车,轻松突破1000公里,石油税收大窟窿谁来填?
IAV:采用移动粒子半隐式(MPS)方法模拟电机中的自由流动冷却液的技术
高压条件下的氢气渗透及暴露聚乙烯的破坏-高压氢气装置用高分子材料的性能
VIP会员
学习计划
福利任务
下载APP
联系我们
微信客服
联系客服
人工服务时间为周一至周五的9:30-19:30
非工作时间请在微信客服留言
客服热线:
4000-969-010
邮箱:
service@fangzhenxiu.com
地址:
北京市朝阳区莱锦创意园CN08座
帮助与反馈
返回顶部