热力学状态方程的几种数值求解方法

给出了4种热力学状态方程求摩尔体积的求解算法,讨论了这些算法的数学原理、收敛特性和具体迭代格式,并以Peng-Robinson 方程为例进行了具体计算和比较.主要结论:简单迭代法收敛稳定但收敛速度慢;Aitken-Steffensen 方法简单、稳定且和Newton 方法具有相同的二阶收敛;Newton 方法速度较快但算法比较复杂,需要计算导数;加速Newton 方法速度最快但最复杂.推荐优先使用Aitken-Steffensen 方法,其次是加速Newton 方法.