soaty 发表于 2013-5-22 23:49:59

采用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并不能收敛,但之前减小改值可以收敛。
水和异丙醇会共沸,求大神帮助谢谢!


已经解决了!感谢楚天湘水的帮助!

tdl522 发表于 2013-5-22 23:49:59

@楚天湘水
兄弟啊,来看看这个帖子可以解决不。

jiaohw910 发表于 2013-5-22 23:49:59

{:1106_365:}反应精馏以前也模拟过一次,貌似也是出了一些问题

楚天湘水 发表于 2013-5-22 23:49:59

@tdl522 ,呵呵,惭愧啊,其实我也懂的不多。
楼主能不能简单介绍下这个反应,有没有共沸生成?还有能不能把动力学方程贴出来,这样方便大家讨论

soaty 发表于 2013-5-22 23:49:59

楚天湘水 发表于 2013-5-23 09:11
@tdl522 ,呵呵,惭愧啊,其实我也懂的不多。
楼主能不能简单介绍下这个反应,有没有共沸生成?还有能不能 ...

我把动力学传到附件里了,还有异丙醇和水是共沸体系。麻烦你看一下是什么原因。另外fortran语言不是我编写的,是我的导师编的,所以这一块的问题并不是很清楚。

楚天湘水 发表于 2013-5-22 23:49:59

看了楼主的动力学子程序,发现有这样一句话RATNET=RATE*REALB(1)*1.D-03,但我打开楼主的BKP文件里面并没有看到realb这个值,楼主的动力学子程序是从RAMTBE的那个子程序改过来,可能没有读懂那个程序,这个地方没有修改。
如果楼主真的需要用到动力学子程序的话,建议楼主把手册认真读个几遍吧,还有楼主写fortran程序的时候最好用个fortran语言的编辑器,这样可以避免很多语法错误。感觉楼主的PARAMETER(K_IPA=1)这句话跑到前面去了,aspen中用的fortran还是老版本,前面6个字符都不会编译。
今天我的编译器出了问题,没办法调试楼主的程序,楼主修改下那句话,试试吧。

soaty 发表于 2013-5-22 23:49:59

楚天湘水 发表于 2013-5-23 09:58
看了楼主的动力学子程序,发现有这样一句话RATNET=RATE*REALB(1)*1.D-03,但我打开楼主的BKP文件里面并没有 ...

谢谢,我也感觉可能是FORTRAN的问题,不过我完全不懂这一块,是导师编的。这个fortran的版本也是导师电脑上的。那我先安装新版本的fortran试试可否。

楚天湘水 发表于 2013-5-22 23:49:59

soaty 发表于 2013-5-23 10:03
谢谢,我也感觉可能是FORTRAN的问题,不过我完全不懂这一块,是导师编的。这个fortran的版本也是导师电脑 ...

不是安装新版本fortran的问题,aspen中的fortran就是这个样子的

soaty 发表于 2013-5-22 23:49:59

楚天湘水 发表于 2013-5-23 10:16
不是安装新版本fortran的问题,aspen中的fortran就是这个样子的

哦~明白您的意思了。那这个问题怎么解决啊。不懂。我以为是我安装的fortran的问题呢。

楚天湘水 发表于 2013-5-22 23:49:59

和楼主讨论了下,子程序已经修改好了,不收敛是因为共沸的问题,大家都来帮帮忙吧。

tdl522 发表于 2013-5-22 23:49:59

楚天湘水 发表于 2013-5-23 10:47
和楼主讨论了下,子程序已经修改好了,不收敛是因为共沸的问题,大家都来帮帮忙吧。

感谢兄弟修改子程序,现已调收敛。

soaty 发表于 2013-5-22 23:49:59

楚天湘水 发表于 2013-5-23 10:47
和楼主讨论了下,子程序已经修改好了,不收敛是因为共沸的问题,大家都来帮帮忙吧。

感谢帮忙,下午的时候已经可以运行错误了。只是很奇怪,生成物水和酯均是从塔底出来,而塔顶出来的水很少,理论上水应该从塔顶出来啊。 我猜想是否与进料温度有关。略微试了下,结果还是不行。恳请回复~~

weihuixiaoniao 发表于 2013-5-22 23:49:59

请问你用那个版本的fortran编译器可以吗?是不是版本不能太高啊?我用的inter visual fortran2011,总是出现如下图的错误:

weihuixiaoniao 发表于 2013-5-22 23:49:59

weihuixiaoniao 发表于 2014-2-27 17:37
请问你用那个版本的fortran编译器可以吗?是不是版本不能太高啊?我用的inter visual fortran2011,总是出 ...

我用的aspen软件是2010年出的aspen plus v7.2

soaty 发表于 2013-5-22 23:49:59

weihuixiaoniao 发表于 2014-2-27 17:43
我用的aspen软件是2010年出的aspen plus v7.2

恩。我用的fortran是以前版本的,好像是6.5的。aspen我用的v7.1

weihuixiaoniao 发表于 2013-5-22 23:49:59

weihuixiaoniao 发表于 2014-2-27 17:43
我用的aspen软件是2010年出的aspen plus v7.2

刚才那个问题解决了,谢谢了,可以配套

管家婆 发表于 2013-5-22 23:49:59

借楼主的程序来学习学习

一直在纠结 发表于 2013-5-22 23:49:59

weihuixiaoniao 发表于 2014-3-6 23:29
刚才那个问题解决了,谢谢了,可以配套

我在做编译的时候也出现了您这个问题,请问您是如何解决的?我用的是V8.4

xsahh 发表于 2013-5-22 23:49:59

{:1106_365:}
页: [1]
查看完整版本: 采用radfrac模拟反应精馏过程,反应选择动力学,但是出现不收敛和质量不守恒的问题..