采用radfrac模拟反应精馏过程，反应选择动力学，但是出现不收敛和质量不守恒的问题..

soaty

RT,我用RADFRAC模拟棕榈酸和异丙醇酯化反应指棕榈酸异丙酯的反应精馏过程，先尝试反应部分选择化学平衡模拟，没有错误，再在相同条件下，采用动力学，由于动力学方程的形式和我查到的表达形式不一致，所以用fortran语言编写了动力学方程导入了aspen中，但是运行的时候出现了不守恒和不收敛的问题，我改了迭代次数和迭代方法，以及蒸馏速率，仍然存在这个问题。求问，是什么原因导致的该问题。会不会是fortran那里的代入有问题？如何解决！！！！！！！ 附件里有fortran的文件，要用simulation engine生成obj文件。fortran采用的应该是6.5的版本，如截图里所示。 采用化学平衡模拟的时候也出现了不收敛，我改变Distillate rate或者是D/F的比值后就可以收敛了，弄成动力学减小D/F并不能收敛，但之前减小改值可以收敛。  
水和异丙醇会共沸，求大神帮助谢谢！


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

tdl522

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

jiaohw910

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

楚天湘水

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

soaty

楚天湘水 发表于 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个字符都不会编译。
今天我的编译器出了问题，没办法调试楼主的程序，楼主修改下那句话，试试吧。

soaty

楚天湘水 发表于 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就是这个样子的

soaty

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

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

楚天湘水

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

tdl522

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

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

soaty

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

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

weihuixiaoniao

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

weihuixiaoniao

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

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

soaty

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

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

weihuixiaoniao

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

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

管家婆

借楼主的程序来学习学习

一直在纠结

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

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

xsahh

{:1106_365:}