KINETIC.f

本帖最后由 天堂弓箭hit 于 2019-1-22 19:48 编辑

我在做甲烷和氧化镍蒸气重整制氢的气固反应的动力学模拟，其中这个反应的动力学方程如图，动力学方程中k（t）是一个定值，g这个函数也是一个定值，f（x）的定义如图，x是固体转化率，其定义如图，就是（完全氧化的质量-当前固体质量）/（完全氧化质量-完全还原质量）。由于这个反应动力学模型aspen没有（应该是这样的吧），所以要编程，我大概写了一下，感觉可能错误比较多，算x（固体转化率），需要调用固体的质量，和固体的分子量，请问我的程序应该怎么改，用什么语句可以调用固体的质量




我这几天大概改了一下，虽然这个程序还是没有能在我的模拟上成功运算，但是固体数据调用问题可以解决了，对于动力学子程序，如果有其他子物流，想要调用其他子物流摩尔流量时，用SOUT(IDXSUB(n)-1+m)就行了，n代表子物流序号，MIXED为1，CLSOILD为2，NC为3，m为组分序号，就是在物性中你定义的组分的顺序，所以调用固体数据时，让n=2或3就行，具体为几取决于固体是常规固体还是非常规固体

例子：如要调用NiO固体的摩尔流量，可以定义该变量为NNiO（变量名可以变），若NiO为第一个组分，则NNiO=SOUT(IDXSUB(2)-1+1)就表示NiO的摩尔流量
此贴错误欢迎各位大神随时指正

学习学习！谢谢分享！

额，这个方程我还看不大懂，分子量应该是定值，可以边方程的时候直接写个常数的。

谢谢提醒，已修改问题

嗯，这个东西我也没有接触过，期待用户自定义模型用得好的人来答吧。

你好，请问这个N指的是哪个流股里边的摩尔流量呀