本算法迭代速度快，一般只需2～6次即可求解，而且精度较高。
图2已知温度T和熵s求解其他参数的程序框图
Fig 2 The flow diagram for calculating the other
parameters under given temperature T and entropy s
3 算例
已知某天然气干线天然气组分如下表所示，输气站中某一压气机的入口状态工况：（1）P=1000K Pa，T=144.869 K；（2）P=1000 K Pa，T=266.048 K；（3）P=10000 K Pa，T=251.947 K；（4）P=1000 K Pa，T=323.100 K；求压气机出口处天然气的熵值。
表1 天然气组分表（%）
| CH4 |
C2H6 |
C3H8 |
i-C4 |
n-C4 |
iC5 |
nC5 |
nC6 |
nC7 |
nC8 |
N2 |
CO2 |
H2O |
| 86.43 |
1.83 |
0.49 |
0.12 |
0.13 |
0.06 |
0.05 |
0.05 |
0.09 |
0.02 |
0.62 |
10.1 |
0.01 |
The Gascomposition table（%）
table1 经过计算得到在四种工况下的新算法和传统算法计算结果比较如下表2：
表2 计算结果的比较
| 已知条件 |
传统算法结果 s(kJ.k-1.mol-1.K-1) |
新算法结果 s* (kJ.k-1.mol-1.K-1) |
迭代次数 |
| P/ (K Pa) |
T (K) |
| 1000 |
144.869 |
225.000 |
225.000 |
4 |
| 1000 |
266.048 |
250.000 |
250.000 |
3 |
| 10000 |
251.947 |
220.000 |
220.000 |
5 |
| 10000 |
323.100 |
235.000 |
235.000 |
3 |
Comparisonof calculation results
table24 结论
以上介绍了两种种计算压气机实际工况的方法。传统算法是很精确的一种算法，但是在计算过程中运用了大量的计算公式，需要求解很多复杂的方程组，往往求解一个参数就需要花费一个小时甚至更长时间。并且在每一次求解时都会产生一定的误差。而在新算法中，运用了C语言编程技术，经过几次迭代就可以求得结果，并且这个程序对任何工况下物性参数的求解可通用，方便、简单、省时，并且通过与传统算法的计算结果比较误差也很小，因此这种结合C语言编程技术的算法计算速度快、精度高，在实际应用中更有优势，更广泛。
参考文献
1 王树立：输气管道设计与管理。化学工业出版社（北京），2006。第2页。
2 姚光镇：输气管道设计与管理。石油大学出版社（东营），1989。第107页。
3 吴玉国，陈保东：BWRS方程在天然气物性计算中的应用，油气储运，2003年，第22卷，第10期。第22～24页。
4 谭浩强：C语言程序设计，清华大学出版社（北京），1998。76～78页。
2/2 首页 上一页 1 2 |
