fortran6.5编写动力学方程有没有步骤？

mhpwhf

醋酸甲酯合成，我编完动力学方程，带入aspen11.1，产物量非常少。我怀疑是方程编写的不对，但我不知道我编写的动力学方程格式对不对。
usrknt.f(8.4KB)

fuwen

这个aspen里面没有吗，还需要用fortran编？

mhpwhf

fuwen 发表于 2013-5-2 13:31

                               
登录/注册后可看大图

这个aspen里面没有吗，还需要用fortran编？

毕业设计要求的

楚天湘水

看了楼主的帖子，里面自己写的语言好像只有r=9863.3*EXP(-48800/8.314/T)，这个应该是反应速率常数K，不是反应速率r,楼主是不是搞错了。你得把你要写的动力学方程贴出来，大家才能帮到你啊。
如果楼主要写fortran的动力学程序，建议楼主参考美国里海大学威廉鲁宾教授和台湾大学余政靖教授合著的书Reactive Distillation Design and Control第九章。
同时http://meng.horse/thread-57962-1-1.html里有关于反应精馏fortran动力学程序的讨论，希望能帮到楼主。

mhpwhf

楚天湘水 发表于 2013-5-2 16:26

                               
登录/注册后可看大图

看了楼主的帖子，里面自己写的语言好像只有r=9863.3*EXP(-48800/8.314/T)，这个应该是反应速率常数K，不是反 ...

r=mcat(k+1xHOAcxMtOH-k-1xMtAcxH2O)

k+1=k0+1exp(-E+1/RT)

k-1=k0-1exp(-E-1/RT)

k+1，k-1 分别为正逆反应的指前因子；E+1，E-1分别为正逆反应的活化能；x 为组分的液相摩尔分数

k0+1=4．944 kmol/(g·s)，k0-1为34．830 kmoL/(g·s)，E+1=48．80 kJ/mol，E-1=58．60 kJ/mol

醋酸进料21.59(mol·h-1),甲醇进料16.85(mol·h-1)         m(Catalyst):m(HOAc)=0.02

楚天湘水

mhpwhf 发表于 2013-5-2 19:19

                               
登录/注册后可看大图

r=mcat(k+1xHOAcxMtOH-k-1xMtAcxH2O)k+1=k0+1exp(-E+1/RT) k-1=k0-1exp(-E-1/RT)k+1，k-1 分别为正逆反应 ...

你得告诉我们你是用的什么反应器或者是采用的反应精馏？这些aspen都有模板，在这个模板的基础上修改。
程序还得你自己写，有了困难后大家再一起讨论。

mhpwhf

楚天湘水 发表于 2013-5-2 20:44

                               
登录/注册后可看大图

你得告诉我们你是用的什么反应器或者是采用的反应精馏？这些aspen都有模板，在这个模板的基础上修改。
程 ...

不好意思，刚学aspen，我模拟的是醋酸甲酯的合成，用RadFrac进行反应精馏。顺便问一下，模板在哪？刚接触aspen，谢谢啦

楚天湘水

mhpwhf 发表于 2013-5-3 08:43

                               
登录/注册后可看大图

不好意思，刚学aspen，我模拟的是醋酸甲酯的合成，用RadFrac进行反应精馏。顺便问一下，模板在哪？刚接触 ...

反应精馏fortran动力学楼主可以参考美国里海大学威廉鲁宾教授和台湾大学余政靖教授合著的书Reactive Distillation Design and Control第九章。

mhpwhf

楚天湘水 发表于 2013-5-3 10:22

                               
登录/注册后可看大图

反应精馏fortran动力学楼主可以参考美国里海大学威廉鲁宾教授和台湾大学余政靖教授合著的书Reactive Dist ...

不好意思，第九章中的fortran程序实在是看不明白

mhpwhf

楚天湘水 发表于 2013-5-3 10:22

                               
登录/注册后可看大图

反应精馏fortran动力学楼主可以参考美国里海大学威廉鲁宾教授和台湾大学余政靖教授合著的书Reactive Dist ...

请问一下，只把动力学编入fortran需要注意什么？