分离工程 第二章 单级平衡过程
分离工程第二章单级平衡过程
第二章 单级平衡过程
Chapter2 Single stage balance process
单级平衡分离是指两相经一次紧密接触达到平衡后随即引离的过程，由于平衡两相的组成不同，因而可起到一个平衡级的分离作用。其相平衡用于阐述混合物分离原理、传质推动力和设计计算。
§2-1汽液相平衡（Vapour-liquid phase equilibrium ）
所谓相平衡是指两个或两个以上的相处于平衡状态。“平衡”的意思是指：在宏观上系统的性质随时间而改变的趋势已达到零。而所谓“相”，是指任何数量的物质在其所占据的空间内宏观性质是均匀一致的，没有不连续的地方。一定数量的物质，即使被分割成若干部分，但只要它们的性质和组成完全一样，则可把它们称为一个“均相”。若有两个或两个以上的均相，虽然它们互相紧密接触，但它们各自的性质并不随时间而改变，通常，就用“相平衡“这一词来表达这一状态。
对于双组分系统温度—压力—组成的平衡关系，常常利用实验来测得，而多组分系统的相平衡关系用实验方法来测定就比较复杂。随着相平衡理论研究的深入，对双组分和多组分系统的汽液相平衡已建立了一些定量的关系式，利用这些关系式，它只需用少量的双组分的实验数据，这就大大地减轻了实验工作量。 一、 汽液相平衡关系（Vapour-liquid phase equilibrium ） 1. 基本关系式
相平衡条件：组分i 在汽液两相中的化学位相等，L i V i μμ= 也可表示为组分i 在汽液两相中的逸度相等，
L i
V i f f ??= i V i V i i V i V i y f y P f 0??γφ== i
L i L i i L i L i x f x P f 0??γφ== i L i i V i x P y P φφ??= 或 i
i L i i V i x f y P γφ0?= ()
RT
P P V P f S
i i S i
S
i L
i
-=ex p 0φ
2. 相平衡常数i k
① 定义
工程中常用相平衡常数来表示相平衡关系
i
i
i x y k =
i k 表示了i 组分在平衡的汽液两相中的分配情况，俗称分配系数。
②i k 的计算方法 ⑴状态方程法
V i
L i
i
i
i x y k φφ??==
只要给出组分i 的汽液两相的分逸度系数，即可求出i k ，而已知i k 则可由()i i y x 求与之相平衡的()i i x y 。
L i φ?和V i
φ?均可用状态方程来计算，但该状态方程必须同时适用于汽液两相，常见的有SRK 、PR 和BWR 方程，此法适用于中压下，液相非理想性不是很强的烃类系统。
⑵活度系数法
V i
L
i i i i i P f x y k φγ?0== 用于只能计算V i
φ?的状态方程，如维里方程，RK 方程，而i γ则由活度系数模型来计算。该法用于压力不高，液相非理想性强的系统。
a. 可凝性组分基准态逸度 当1→i x 时，1→i γ，()
RT
P P V P f S
i i S i
S
i L
i
-=ex p 0φ
即纯液体i 在T 、P 下的逸度等于饱和蒸汽压乘以两个校正系数，S i φ为校正处于饱和蒸汽压对理想气体的偏差，指数校正（Poynting ）因子是校正压力偏离饱和蒸汽压的影响。
b. 不凝性组分基准态逸度 当0→i x 时，
1*
→i γ，i
L i
x L
i
x f H f i ?lim
0→==
或亨利定律：i
L i Hx f =?
③i k 与ij α的关系
相对挥发度ij α的定义是i ，j 两组分的相平衡常数之比。
固有分离因子：j
i j i j j i
i j i ij x x y y x y x y k k =
==α 若1=ij α，表示汽液两相中i 、j 两组分的浓度之比相等，因此不能用一般的精馏来分离，ij α值越大，两相平衡后的比值
∑=
i
ij
ij
i x k α
α ∑=
i
ij
i
ij i x x y α
α
∑
=ij
i
ij i y k αα ij
i ij
i i y y x αα∑=
二、 汽液平衡的分类与计算（Classification and calculation of vapour-liquid phase
equilibrium constant ） 1. 汽液平衡的分类
基本关系式：i
i S i S i i i L i i V i x P x f y P γφγφ==0? 忽略Poynting 因子，S
S L p f φ=0
2. 汽液平衡常数的计算 ① 逸度系数的计算
计算逸度系数可以通过气体状态方程，PVT 实验数据，普遍化关系，可根据情况任选一种，其基本关系为： ⑴纯组分的逸度系数
?-=P
i i P
dP Z 0
)
1(ln φ 当2≥r V ,)(Pr ln 10
B B Tr
i ωφ+=
6.100422083.0r T B -
= 42
1
172.0139.0r
T B -= 当2 )(ln ,ln ,10r r i i P T f =φφ可查图。 ⑵二元混合物的逸度系数 ?-=P dP Z i i )1(?ln φ a. 维里方程 RT P B B y M j ij j i ??? ?????-???? ??=∑2?ln φ ∑∑=i j ij j i M B y y B 10B B RT BP C C ω+= ij B 用混合规则求，即先求ij Cij Cij p T ω,,。见 P34和化工热力学。 b. RK 方程 ??? ??-??? ?????+-???? ??+++--+-=∑RT PV b V b V b V b RT b a V b V b RT a y b V b b V V m m m m m m m m i m m m m m k ik k m m i m m m i ln ln )ln(2ln ?ln 2 2323φ ∑=i i i m b y b ∑∑=j ij j i i m a y y a 书中给出了Vdw Equ.的计算，见P30。 ② 活度系数的计算 若1>i γ，称为对拉乌尔定律有正偏差，1 大多数非烃类混合物，都表现为正偏差，少数由于有缔合现象或属于电解质溶液，能发生负偏差。一般可以液相中分子间的作用力来估计偏差的正负。 目前，还不能完满地定量计算活度系数，对于双组分系统，主要通过实验测定，也可采用某些公式计算，对于三组分或更多组分的非理想溶液，实验数据很少，活度系数主要用公式估算。 ()∑=i i i E x RT G γln () j n P T i E i n RT nG ,,ln ????????=γ 只要知道E G 的数学模型，就可通过组分i 的摩尔数n i 求偏导数得到i γ的表达式。 ⑴二元溶液的活度系数 a. Van-Laar Equ. 2 22111212 11lg ??? ? ??+= x A x A A γ 2 11222121 21lg ??? ? ??+= x A x A A γ 讨论 1）A 的物理意义 ∞→==110 12lg lg lim 1γγx A ∞→==2 20 21lg lg lim 2γγx A 2）由2112,lg A A i →γ 2 112 2112lg lg 1lg ???? ? ?+=γγγx x A 2 221 1221lg lg 1lg ??? ? ? ?+=γγγx x A 3）当A A A ==2112时，此二元系统称对称系统，方程可变为单参数的对称方程： 2 222121lg Ax x A ==γ；2121212lg Ax x A ==γ 4）当02112==A A 时，1=i γ为理想体系； 当0,02112<>A A 时，1>i γ为正偏差非理想体系 A 可用来判别实际溶液与理想溶液的偏离度。 b. Margules Equ. ()[]12211122 2 12lg A A x A x -+=γ ()[]21122212122lg A A x A x -+=γ 10 12lg lim 1γ→=x A 20 21lg lim 2γ→=x A 1 2 12 21212lg 2lg x x x x A γγ+-= Van Laar 方程和Margules 方程有悠久的历史，仍有实用价值，特别是定性分析方面。 优点：数学表达式简单；容易从活度系数数据估计参数； 非理想性强的二元混合物包括部分互溶物系，也经常能得到满意的结果。 缺点：不能用二元数据正确推断三元系的活度系数。不能用于多元系相平衡计算。 c. Wilson Equ. ∑∑ ∑===--=n j n k n j kj j ki k ij j i x x x 1 1 1 ln 1ln λ λλγ ?? ? ???+++-+-=2211212122111221)ln(1ln x x x x x x x x i λλλλγ )exp(RT g g V V ii ij L j L i ij -- = λ 注意：1）当01→x 时，21121ln 1ln λλγ--=∞；当02→x 时，12212 ln 1ln λλγ--=∞ 2）当1,12112>>λλ时，为理想体系； 当1,12112<>λλ时，为负偏差非理想体系；当1,12112<>λλ时，为正偏差非理想体系； λ可用来判别实际溶液与理想溶液的偏离度。 3）ji ij g g =，但ji ij λλ≠ 该方程考虑了T 、P 的影响，汽液平衡计算有较高的精度；适用范围很广，不能用于液液平衡的计算，不能用于部分互溶系统。 d. NRTL 方程和UNIQUAC 方程 根据局部摩尔分率的概念建立的。NRTL 方程能进行汽液平衡和液液平衡的计算；但方程中每对二元系多了第三参数α12；Uniquac 方程有NRTL 的优点，但数学表达式最复杂；适用于分子大小相差悬殊的混合物。 e. UNIFAC 基团贡献法0=E S 正规溶液 ⑵三元溶液的活度系数 a. Margules Equ. ()[][]+-++-+=)(22lg 13311132312211122 2 1A A x A x A A x A x γ [])21()(2)(21233231331132132132x C A A x A A x A A A x x ---+-+-++ )]21()(2123323x C A A x ---+ 2 13 3132231221A A A A A A C -+-+-= ij A 为有关的双组分溶液之端值常数，可查阅手册，顺序轮回替换下标，用2代1，用3代2，1代3，便可求2γ及3γ。 b. Wilson Equ. )exp(RT g g V V ii ij L j L i ij -- = λ ? ? ????++++++++-++-=3322311313233221121 2133122111331221)ln(1ln x x x x x x x x x x x x x x x i λλλλλλλλλλγ 只需查出有关三个二元溶液的威尔逊参数，计算结果较好，虽较繁复，随着计算机的使用，目前应用较多。 c.NRTL 及UNIFAC 法 ③i k 的简化计算：P —T —K 图 相平衡常数是温度、压力和汽液组成的函授，无论用状态方程还是用活度系数模型，其计算工作量都很大，必须借助于计算机辅助计算。烃类物系在化工中十分重要，其行为接近理想情况，可仅考虑P 、T 对i k 的影响，迪普里斯特（Depriester ）以BWR 方程为基础，经广泛的实验和理论推算，作出了轻烃类的PTK 列线图，见P38图2-1，这些图虽然没有假设理想溶液这个条件，但在图上所示的有限的压力范围内，组成对K 值的影响很小，仍然把K 看成是T 、P 的函数，平均误差为8-15%。适用于0.8-1Mpa （绝对压）以下的较低压区域。例2-1，2。 §2-2 多组分物系的泡点和露点计算（Calculation about bubble-point and dew-point of multicomponent system ） 泡点温度（压力）是在恒压（温）下加热液体混合物，当液体混合物开始汽化出现第一个气泡时的温度（压力），简称泡点。 露点温度（压力）是在恒压（温）下冷却气体混合物，当气体混合物开始冷凝出现第一个液滴时的温度（压力），简称露点。 根据泡、露点的概念，精馏塔塔顶温度即为对应塔顶汽相组成的露点，塔釜温度即为对应塔釜液相组成的泡点。由于气体中出现的露点液和液体中出现的泡点气泡，其物质量仅为原混合物的极为微小部分，因而物系形成了第一个液滴和第一个气泡后，原有物系组成并未改变。在露点和泡点下均会出现汽液两相，可以通过泡点、露点的计算，了解在该温度下的汽液平衡组成。 一、 泡点温度和压力的计算（Calculation of temperature and pressure at bubble-point ） 1. 泡点计算与有关方程 已知)(,T P x i ，求，泡点压力）泡点温度，（P T y i ,, C C f =+-=2π 即只要给定C 个变量，整个系统就规定了，可以利用相平衡关系计算。 ① 相平衡关系：i i i x k y = ② 浓度总和式：1=∑=∑i i i x k y ③ 相平衡常数关联式：),,,(i i i y x P T f k = 给定P （或T ）和C-1个x i ，则上述方程有唯一解。 2. 计算方法 ① 平衡常数与组成无关的泡点计算 :),(P T f k i = a.PTK 列线图； b.)ex p(1 T C B A P k i i i i -- = )(,,T f C B A i i i 物性，= ⑴手算 思路：设i i i B Y i i i i i T P x k y T T x k x k k P T ==?→?≤-∑→∑→??? →?，求已知ξ1)()( ） N ξ为试差的允许偏差，手算中一般取0.01~0.001。 按初设温度T 所求得的i i x k ∑值若大于1，表明所设温度偏高，B T T >，反之若小于1，则表明所设温度偏低。这是因为P 一定，T 升高，k 增加。 如何调整T （P ）值，为避免盲目性，加速试差过程的收敛可用下法： i iG G i G i G i i x k x k k k x k α∑=∑ =∑)( G 表示对i i x k ∑值影响最大的组分，iG α表示i 组分对G 组分的相对挥发度，在一定的温度范围内const iG ≈α。 const x k k i i G =∑1 亦即对各次试差：1)1 ()1( -∑=∑m i i G m i i G x k k x k k m 为试差序号，为使第m 次试差时1)(=∑m i i x k 1)( -∑=m i i G Gm x k k k 由该Gm k 便可以从K 图读出第m 次试差时应假设的温度m T 值，按上述方法通常经过2-3次试算便可求得解。 对压力的迭代还可采用下面的方式： i i r r x k P P ∑=+B 1B ，此式可自动调整压力P 。 P47例2-3。 ⑵计算机计算 a ．)exp(1T C B A P P p k i i i S i i --== Hayden 提出适合于烃类的经验方程。 3423213 1T b T b T b b T k i i i i i +++= 常见物质的mi b 的数据可查表，见“多元汽液平衡与精馏”。 b ．迭代方法，Newton-Raphson 迭代求解。 目标函数：0)(=x F ，求根， 迭代公式：) () (1 m m m m x F x F x x '-=+ 泡点计算目标函数：01)(=-∑=i i x k T F ，泡点方程 dT dk x T F i i ∑=')( ) () (1m m m m T F T F T T '- =+ 自学P48例2-4。 Richmond 迭代法： ) () (2)()(22 1m m m m m m T f T f T f T f T T '''-'- =+ ②平衡常数与组成有关的泡点计算 ),,,(i i i y x P T f k =，用于系统非理想性较强时,需计算混合物中i 组分逸度系数或 活度系数，而V i φ?是i y 的函数,i γ是i x 的函数，在i y 或i x 未求得之前无法求得V i φ?或i γ值，于是计算时还需对V i φ?，i γ进行试差,对泡点计算,由于已知i x ，除需送代泡点温度或压力外，还需对V i φ?进行试差。 i i S i S i i P P k φ γφ?= 设i i i B Y S S T P x k y T T x k x k y k P P T ==?→?≤-∑→→=→??→????→?=，设已知ξφγφφ1?,,)(1?)( N 由于i k 主要受温度影响，且i k ln 与T 1近似线性关系，故判别收敛判据为： 0ln )1(=∑=i i x k T G 计算框图见P50图2-5，P53图2-6 参考文献：“极性物系多组分分离泡点计算”，齐鲁石油化工，91（3），208 P54图2-7给出了状态方程法计算泡点压力的框图。 天然气化工，1995年第1期P52介绍了“一种新的泡点计算方法”该法使目标函数与自变函数的关系变为单调函数的线性关系,即使目标函数线性化，从而消除了极值点和拐点。由于形式简单，技术也简单，收敛非常迅速，一般经过2-4次迭代即可达到收敛精度，且初值范围大。 自学例2-5，2-6。露点温度和压力的计算（Calculation of temperature and pressure at dew –point ） 1.露点计算与有关方程 已知)(,T P y i ，求，露点压力）露点温度，（P T x i ,, ①相平衡关系：i i i k y x = ②浓度总和式：1=∑ =∑i i i k y x ③相平衡常数关联式：),,,(i i i y x P T f k = 2.计算方法 ① 平衡常数与组成无关的泡点计算 ⑴手算 思路：设i i i D Y i i i i i T P k y x T T k y k y k P T ==?→?≤-∑→∑→??? →?，求已知ξ1)()( ） N ξ为试差的允许偏差，手算中一般取0.01-0.001。 按初设温度T 所求得的01<-∑i x ，表明所设温度偏高，D T T >，反之若 01>-∑i x ，则表明所设温度偏低。达到允许误差时，所设温度即为所求的露点。 调整T 或P ：1)(-∑?=m i i G Gm k y k k ⑵计算机计算 露点方程：01)(=-∑ =i i k y T F ， )( )(2 dT dk k y x T F i i i i ∑-=' )() (1m m m m T F T F T T '-=+ 或0ln )(=∑=i i k y T φ i i k y T f T ∑ '=') ()(φ ) (ln 1T f x x T T i i m m '∑∑- =+ 收敛 速度快。 ②平衡常数与组成有关的露点计算 求算露点，由于已知i y ，则需迭代露点温度或压力外，还需对i γ进行试差。 设T T x k y x k P P T D Y i i i i i i S i i S i T P i =?→?≤-∑→→=→??→????→?=ξγφφγ1,?,)(1)(设已知 比较不变 N 调整T （P ） 自学例题2-7。 §2-3闪蒸过程的计算（Flash calculation ） 平衡蒸馏又称闪蒸是连续、稳态的单级蒸馏过程，该过程使进料混合物经加热、冷却或降压部分汽化或冷凝得到含易挥发组分较多的蒸汽和含难挥发组分较多的液体，即一个平衡级过程。在分离流程中常遇到的部分汽化和冷凝，绝热闪蒸以及部分互溶系统精馏塔顶蒸汽冷凝后的分层过程均属这类单级分离过程。 一、 等温闪蒸（Isothermal Flash ） 1. 混合物的相态 在进行闪蒸计算时，需判断混合物在指定温度和压力下是否处于两相区，为此需对进料作如下检验。 = 1 B T T = 进料处于泡点，0=ψ i i Z k ∑ > 1 T >B T 可能为汽液两相区，ψ>0 < 1 T = 1 D T T = 进料处于露点，1=ψ i i k Z ∑ > 1 T < 1 T >D T 进料为过热蒸汽 计算时，所指定的温度应在泡露点温度之间，这样才会出现汽液两相，否则只会是单相不必进行等温闪蒸的计算。若)(i i k Z ∑< 1，所指定的温度高于露点温度。若i i Z k ∑< 1，所指定的温度低于泡点温度。则所指定的温度下不可能实现闪蒸。 只有)(i i k Z ∑> 1和i i Z k ∑< 1时，混合物始处于汽液两相区（0<ψ<1），需要进行等温闪蒸的计算。 此外可分别用泡露点方程计算闪蒸压力下的泡露点温度，若D T >T >B T ，则闪蒸成立。 2. 等温闪蒸过程 气体混合物的部分冷暖和液体混合物的部分汽化如图所示，两类计算： ① 已知：P T Z F i ,,,，（指定工艺条件） 求：Q x L y V i i ,,,, 液化率：F L e == 总加入量液化的量 汽化率：F V ==ψ总加入量 液化的量 1=ψ+e ，该情况俗称等温闪蒸。 ②已知：)(,,,e P Z F i ψ，指定汽化率求操作条件，即求：Q L V x y T i i ),(,,, 3. 基本计算公式 由于进、出口量是相等的及出口处的汽液两相达到平衡，因而可通过物料平衡和相平衡来计算。 ① 组分的物料衡算式 i i i Lx Vy FZ L V F +=+=, ②相平衡关系式 i i i x k y = ③相平衡常数关联式 ),,,(y x P T f k i = ④浓度总和式 1,1=∑=∑i i y x 关系方程 )1()1(ψ-+ψ=+-= i i i i i k Z e k e Z x ) 1(1)1(-ψ+=+-= i i i i i i i k Z k e k e Z k y 4. 计算方法 ① 汽液平衡常数与组成无关 ⑴第一类计算 思路：结束 设已知 ?→? ????=∑=∑→y i i i T P y x x 11, 要注意e 、Ψ的初值的选取不能使i i i Z y x =或。 当指定i i Z k P T 和被指定)(,求解e 或Ψ时，由于需求解高度非线性方程，为避免试差计算的盲目性，可采用计算机计算中常用的迭代求根法，如Newton-Raphson 法，按该法的迭代公式对Ψ可写出 ) () (1k k k k F F ψ'ψ- ψ=ψ+ 目标函数 1) 1(11)(--ψ+∑ =-∑=ψi i i k Z x F [] 2 )1(1) 1()(-ψ+-∑ -=ψ'i i i k k Z F 迭代公式 []2 1 )1(1)1(1 ) 1(1-ψ+-∑ --ψ+∑ + ψ=ψ+i i i i i k k k k Z k Z 按指定的T 、P 求出i k 值并假设Ψ的初值1ψ，便可以计算)(),(ψ'ψF F ，并求出下次迭代用的2ψ值依次重复进行直至)(ψF 小于指定的允许偏差ψε为止。 同样需注意Ψ的初值的选取，以免发散，可取 B D B T T T T --= ψ为初值，此外，应选择适宜的迭代方法，而同一迭代方法若采用不同 的目标函数对收敛速度和稳定性有很大影响。目标函数可任选取，如0=∑-∑i i x y 1)1(1=-ψ+∑ =∑i i i k Z x ， 1) 1(1=-ψ+∑=∑i i i i k Z k y 0)1(1)1()()(2 =-ψ+-∑=-∑=ψi i i i i k k Z x y F [] 2 2 )1(1)1()(-ψ+-∑ -=ψ'i i i k k Z F ) () (1ψ'ψ- ψ=ψ+F F k k 该式接近于一线性单调收敛函数，其收敛性能甚佳，对初值1ψ的选择并无特殊要求，无论在Ψ根的左方或右方均能迅速求得解。 例题2-8。 计算框图： 开始 输入Z F P T ,,, 计算),(P T F k i = 计算泡点B T 打印过冷液体?→? >Y B T T 计算露点D T ?→? B T T -= ψ (ψF []?→?<ψY F ε)( 打印，结束 ) () (1ψ'ψ- ψ=ψ+F F k k [] 2 2 1(1)1()(-ψ+-∑ -=ψ'i i i k k Z F ⑵第二类计算 思路：结束设已知 已知?→? ????=∑→=∑→ψy i i i P y x x 11 为有助于收敛，选用Rachford-Rice 方程核实假定值是否正确为佳。此外，可用下式估计T ： ) (1) ()(1k k G k G T f d T k T k ?+= + G k 为基准组分的平衡常数，d 为阻尼因子，0.1≤d 在上述两类计算中加热或冷凝的热量，均可由热量衡算得到。 L V F LH VH FH Q --= 设为理想溶液，则i i F Z H H ?∑= i i V H y H ∑=， i i L h x H ∑= “三组分平衡闪蒸中v k Z --的关联式”，石油化工，1993，（3）：170 “多元平衡闪蒸计算的新方法”，石油化工，1988，17（4）：217 ② 汽液平衡常数与组成有关的闪蒸计算 当i k 不仅是温度和压力的函数而且还是组成函数时，计算更繁复，由于 ),,,(y x P T f k i =，因此需先按),(P T f k i =初估i i y x ,，迭代求Ψ至收敛，再估算新的 一组x 和y ，并计算k ，重新迭代Ψ，直至x 和y 没有变化为止。 设→??????=-=-→=?→???????=-=-→=→ψ0101),,,(101),(i i i y i i i y x y x P T f k y x P T f k 结束 调Ψ 调Ψ P63图2-9 ，自学例2-9 二、绝热闪蒸（Adiabatic flash ） 1. 等焓节流过程及其计算内容 流量为F ，组成为i Z 的液相进料于压力1P ，温度1T 下经节流阀在绝热情况下减压至2P ，于节流过程中部分液体依靠进料本身携带的热量汽化，而系统温度降至2T ，由于上 述节流过程是在绝热情况下进行的（0=Q ），因而节流前后混合物的焓21H H =，节流生成的汽相和剩余液相分离时假设达到平衡。 计算内容：已知：211,,,,P T P Z F i ；求：i i y x L V T ,,,,2 等焓节流和等温闪蒸都能产生汽液平衡相，属单级平衡分离，但所用的分离剂不同。 2. 计算公式 组分的物料衡算式 i i i Lx Vy FZ L V F +=+=, 相平衡关系式 i i i x k y = 相平衡常数关联式 ),,,(y x P T f k i = 浓度总和式 1,1=∑=∑i i y x 焓平衡方程 21H H = i i H Z H ∑=1， i i V H y H ∑=， i i L h x H ∑= L V L V H H H LH VH FH )1(,22ψ-+ψ=∴+= L V L L V L H H H H H H H H --=--= ψ12 关系方程 )1()1(ψ-+ψ=+-= i i i i i k Z e k e Z x ) 1(1)1(-ψ+=+-=i i i i i i i k Z k e k e Z k y 3. 计算方法 ① 作图法 思路：先假设一个2t （2t 应在B D T T 与之间），这样可按闪蒸来计算产生的汽液两相组成和量，即得出ψ-T 关系，闪蒸曲线，然后再由进出料热焓相等的原则来校核2t ，即等焓平衡线，交点为2t ，Ψ。 ⑴闪蒸曲线（ψ-T 曲线） a 、 试差求2P 下的B D T T ,； b 、 在B D T T ,之间设22,t t '''…，试差求21,ψψ…。 ⑵等焓平衡线 a 、 在上述B D T T ,之间所设22,t t '''…下，由求出的 i i y x ,并根据22,T P 查出iL iV H H ,，求出L V H H ,，L V L H H H H --= ψ1*1 b 、作*2ψ-T 等焓平衡线，交点即为2t ，Ψ，由此求出i i y x L V ,,, ② 手算或计算机计算 思路：?→?=-→??→?→ψ???→???→?→-y H H i i R N T H H H y x k H T P L V 0,1 22,211,12迭代设结束 本计算为两层迭代循环，内层为Ψ循环，用Newton-Raphson 迭代求解汽化分率Ψ，收敛精度005.0)(<ψF 。外层则为2T 循环，用正割法调整节流后的温度，需二点温度： ) ()() ()(2222222T f T f T f T T f T T '-'''''-'''= '''；0)(122=-=H H T f 收敛精度 005.0) (1 12<-H H H 211,,,,P P T Z F i i i i H Z H P T f H ∑==111),,( 1 设2T ' =ψ T n + 1 END 该法适合于宽沸程混合物的闪蒸计算，所谓宽沸程混合物，是指构成混合物的各组分的挥发度相差悬殊，在很宽的温度范围内，易挥发组分主要在蒸汽相中，而难挥发组分主要留在液相中。进料热焓的增加将使平衡温度升高，但对汽液流率V和L几乎无影响。因此宽沸程的热量衡算更主要地取决于温度，而不是Ψ。根据序贯算法迭代变量的排列原则，最好是使内层循环中迭代变量的收敛值将是下次内层循环运算的最佳初值，因此的宽沸程闪蒸，由于Ψ对T的取值不敏感，所以Ψ作为内层迭代变量较为合理。见P65图2-10。 此外，估计新的闪蒸温度，除用正割法外，仍可用前述的牛顿法。 对于窄沸程闪蒸，由于各组分的沸点相近，因而热量衡算主要受汽化潜热的影响，反映在受汽相分率的影响，改变进料热焓会使汽液相流率发生变化，而平衡温度没有太明显的变化。显然，应该通过热量衡算计算Ψ（即V和L），解闪蒸方程式确定闪蒸温度。并且，由于收敛T值对Ψ值不敏感，故应在内层循环迭代T，外层循环迭代Ψ。其窄沸程绝热闪蒸的收敛方案见P66图2-11。 此外，书中还介绍了同时收敛法，请自学。 “多组分多级绝热闪蒸闪蒸过程的研究”，化学工程，1994，（1）：19