分离工程2.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 、 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * ① 松弛法的基本方程 F为进料量，若各板的存液量相等，各板的松弛系数也可取等值，但对不同组分有时可用不同的松弛系数值 称为松弛系数 松弛法基本方程 ② 计算步骤 ⑴假定一组初值 选 直至稳定为止 判断 ⑵由松弛法基本方程计算新的 （设为恒摩尔流） 或 设一组各板的 ），泡点试差得 （设 设一组各板的 ② 计算步骤 ⑶由S方程求 不成立则重复⑵～⑸式 ⑹判断 ⑸由H方程计算 计算 ， ⑷由 ③ 优缺点 优点:适用于各种复杂的精馏过程，对初值选定没有严格要求，且某一组分的量发生变化时，对计算影响不大，所以收敛很稳定 但必须使用计算机，且收敛速度也太慢，当其他方法无法解时，才用松弛法 改进:石川、平田等人进行了修正，但收敛速度慢 由丁惠华提出的新松弛法则在Rose松弛法的基础上，经过适当处理允许采用较大的松弛系数，克服了收敛速度慢的缺点，可用于非理想物系的精馏计算 2.3 精馏塔的模拟计算方法的改进 各种平衡级精馏模拟方法的区别仅在于联解非线性的MESH方程组F X 0时采用的方法和步骤互异 现行的严格算法优缺点： 1 方程解耦法 该法简单，易编程，占用内存少，但若组份间交互作用较显著和流程较复杂，此方法误差较大或不收敛 2 同时校正法。 即Newton-Raphson N-R 法充分考虑了组份间的交互作用，适用于各种类型的复杂精馏过程。但它只具局域收敛性，若无好的初值就难以收敛 2.3 精馏塔的模拟计算方法的改进 3 松弛法。该法的稳定性高，不受初值影响，但只在最初几次迭代内收敛较快，越接近真实解收敛越缓慢。 各种方法都存在一些缺陷，需开发一种可靠而有效的通用型算法 解决方法: 将现有的算法进行组合，相互取长补短，使其适合于求解复杂的问题 2.3.1松弛法与N-R法的联合算法 2.3.1.1研究方法 先利用稳定性最好的松弛法进行计算，待精馏过程变量落入N-R法的收敛域之后，就作为高质量的初值由N-R方法快速地进行收敛计算 2.3.1.2需解决的问题 1 需要有切实可行的方法加快松弛法的收敛速度 2.3.1松弛法与N-R法的联合算法 2 设计可靠的决策系统使计算过程能及时自动地由松弛法切换到N-R法 3 需要寻找更通用而有效的数学求解方法 2.3.1.3解决策略 1 将松驰因子的选取方式进行了改进，即采取每步迭代自动更换最优μ值的方法，使松驰法很快就会得出能进入N-R法收敛域的迭代变量值，大大加快了松驰法的收敛速度 2.3.1松弛法与N-R法的联合算法 2 利用一个有效的开关不等式及时地将计算过程切换到N-R方法，判断改进松驰法所产生的初值是否进入N-R法的收敛域 3 联合算法在转入N-R法进行计算时，引入了求解带外有非零元素的近似三对角 块状 矩阵的修正的H-S算法，可成功地应用于求解复杂内联塔系的H-R线性方程组 吴松涛，江青茵，曹志凯.精馏模拟计算的改进新松弛-泡点法，计算机与应用化学，2005， 22（8）：655-658 孙德慧，新松弛法与泡点法结合用于低温含氢体系的模拟计算， 吉林化工学院学报，2000年02期 参考文献 2.3.2 自由度N C十2 同时校正法 2.3.2.1同时校正法的优缺点 该法充分考虑了组分间的交互作用，适用于各种类型的复杂精馏过程，且当初值较接近于解时收敛迅速 但在计算雅可比矩阵各元素时，由于描述精馏过程方程组的高度非线性，通常需用差分法求定一阶偏导数，运算工作量和存储量非常庞大，计算准确性差，限制了该算法联立求解大规模精馏问题的能力 2.3.2 自由度N C十2 同时校正法 2.3.2.2解决策略 1 对原算法的自由度N（2C＋1）数学模型进行化简，代之以自由度N（C＋2）数学模型为核心的压缩型块三对角线方程组联立算法，从根本上降低模型的方程数和变量数 2 对相平衡常数K和焓值H进行多项式回归，用解析法代替差分法求定一阶偏导数，可减少内存单元数和一阶偏导数的运算工作量，改善联立求解的效率和稳定性 2.3.3 二对角矩阵法 2.3.3.1三对角矩阵法的优缺点 占用计算机内存少、无需计算偏导数矩阵 对强非理想体系和沸点范围宽的体系，有时计算不收敛 原因: 由于收敛前由物料衡算三对角矩阵解出的组成 圆整后的组成又不满足物料衡算的要求 不能满足 2.3.3 二对角矩阵法 2.3.