采用radfrac模拟反应精馏过程,反应选择动力学,但是出现不收敛和质量不守恒的问题..
本帖最后由 soaty 于 2013-5-26 00:41 编辑RT,我用RADFRAC模拟棕榈酸和异丙醇酯化反应指棕榈酸异丙酯的反应精馏过程,先尝试反应部分选择化学平衡模拟,没有错误,再在相同条件下,采用动力学,由于动力学方程的形式和我查到的表达形式不一致,所以用fortran语言编写了动力学方程导入了aspen中,但是运行的时候出现了不守恒和不收敛的问题,我改了迭代次数和迭代方法,以及蒸馏速率,仍然存在这个问题。求问,是什么原因导致的该问题。会不会是fortran那里的代入有问题?如何解决!!!!!!! 附件里有fortran的文件,要用simulation engine生成obj文件。fortran采用的应该是6.5的版本,如截图里所示。 采用化学平衡模拟的时候也出现了不收敛,我改变Distillate rate或者是D/F的比值后就可以收敛了,弄成动力学减小D/F并不能收敛,但之前减小改值可以收敛。
水和异丙醇会共沸,求大神帮助谢谢!
已经解决了!感谢楚天湘水的帮助!
@楚天湘水
兄弟啊,来看看这个帖子可以解决不。 {:1106_365:}反应精馏以前也模拟过一次,貌似也是出了一些问题 @tdl522 ,呵呵,惭愧啊,其实我也懂的不多。
楼主能不能简单介绍下这个反应,有没有共沸生成?还有能不能把动力学方程贴出来,这样方便大家讨论 楚天湘水 发表于 2013-5-23 09:11
@tdl522 ,呵呵,惭愧啊,其实我也懂的不多。
楼主能不能简单介绍下这个反应,有没有共沸生成?还有能不能 ...
我把动力学传到附件里了,还有异丙醇和水是共沸体系。麻烦你看一下是什么原因。另外fortran语言不是我编写的,是我的导师编的,所以这一块的问题并不是很清楚。
看了楼主的动力学子程序,发现有这样一句话RATNET=RATE*REALB(1)*1.D-03,但我打开楼主的BKP文件里面并没有看到realb这个值,楼主的动力学子程序是从RAMTBE的那个子程序改过来,可能没有读懂那个程序,这个地方没有修改。
如果楼主真的需要用到动力学子程序的话,建议楼主把手册认真读个几遍吧,还有楼主写fortran程序的时候最好用个fortran语言的编辑器,这样可以避免很多语法错误。感觉楼主的PARAMETER(K_IPA=1)这句话跑到前面去了,aspen中用的fortran还是老版本,前面6个字符都不会编译。
今天我的编译器出了问题,没办法调试楼主的程序,楼主修改下那句话,试试吧。 楚天湘水 发表于 2013-5-23 09:58
看了楼主的动力学子程序,发现有这样一句话RATNET=RATE*REALB(1)*1.D-03,但我打开楼主的BKP文件里面并没有 ...
谢谢,我也感觉可能是FORTRAN的问题,不过我完全不懂这一块,是导师编的。这个fortran的版本也是导师电脑上的。那我先安装新版本的fortran试试可否。 soaty 发表于 2013-5-23 10:03
谢谢,我也感觉可能是FORTRAN的问题,不过我完全不懂这一块,是导师编的。这个fortran的版本也是导师电脑 ...
不是安装新版本fortran的问题,aspen中的fortran就是这个样子的 楚天湘水 发表于 2013-5-23 10:16
不是安装新版本fortran的问题,aspen中的fortran就是这个样子的
哦~明白您的意思了。那这个问题怎么解决啊。不懂。我以为是我安装的fortran的问题呢。
和楼主讨论了下,子程序已经修改好了,不收敛是因为共沸的问题,大家都来帮帮忙吧。 楚天湘水 发表于 2013-5-23 10:47
和楼主讨论了下,子程序已经修改好了,不收敛是因为共沸的问题,大家都来帮帮忙吧。
感谢兄弟修改子程序,现已调收敛。
楚天湘水 发表于 2013-5-23 10:47
和楼主讨论了下,子程序已经修改好了,不收敛是因为共沸的问题,大家都来帮帮忙吧。
感谢帮忙,下午的时候已经可以运行错误了。只是很奇怪,生成物水和酯均是从塔底出来,而塔顶出来的水很少,理论上水应该从塔顶出来啊。 我猜想是否与进料温度有关。略微试了下,结果还是不行。恳请回复~~ 请问你用那个版本的fortran编译器可以吗?是不是版本不能太高啊?我用的inter visual fortran2011,总是出现如下图的错误: weihuixiaoniao 发表于 2014-2-27 17:37
请问你用那个版本的fortran编译器可以吗?是不是版本不能太高啊?我用的inter visual fortran2011,总是出 ...
我用的aspen软件是2010年出的aspen plus v7.2 weihuixiaoniao 发表于 2014-2-27 17:43
我用的aspen软件是2010年出的aspen plus v7.2
恩。我用的fortran是以前版本的,好像是6.5的。aspen我用的v7.1 weihuixiaoniao 发表于 2014-2-27 17:43
我用的aspen软件是2010年出的aspen plus v7.2
刚才那个问题解决了,谢谢了,可以配套 借楼主的程序来学习学习 weihuixiaoniao 发表于 2014-3-6 23:29
刚才那个问题解决了,谢谢了,可以配套
我在做编译的时候也出现了您这个问题,请问您是如何解决的?我用的是V8.4 {:1106_365:}
页:
[1]