闪蒸过程的计算

lng闪蒸量计算
"LNG闪蒸量"是液化天然气（LNG）在一定条件下发生的液化气体的蒸发量，通常是在储存或输送LNG的过程中考虑的参数之一。

计算LNG闪蒸量需要考虑温度、压力等因素，可以使用一些基本的物性和热力学公式。

下面是简化的计算LNG闪蒸量的步骤：
1.计算初始液相体积：确定初始的LNG液相体积（V_initial）。

2.确定储存或输送过程中的温度和压力：确定在储存或输送LNG的过程中的温度（T）和压力（P）。

3.使用物性数据计算饱和蒸汽压：利用LNG的物性数据，如饱和蒸汽压的相关数据，计算在给定温度下LNG的饱和蒸汽压。

4.根据蒸汽压计算蒸发量：使用饱和蒸汽压等数据，可以利用相应的热力学公式计算LNG在给定条件下的闪蒸量。

请注意，具体的计算需要考虑到LNG的具体物性，以及在实际运输和储存中可能存在的各种因素。

理论上，闪蒸量的计算可能涉及到热力学方程、物性数据以及具体系统的工程参数。

在实际应用中，可能需要借助专业软件或者实验数据来进行更准确的计算。

第三节 闪蒸过程的计算2.3 等温闪蒸和部分冷凝过程流程示意图：闪蒸过程的计算方程（MESH ） ⑴物料衡算----M 方程： C 个⑵相平衡--------E 方程： C 个⑶摩尔分率加和式---S 方程： 2个⑷热量平衡式-------H 方程： 1个变量数：3C+8个 (F, F T ,F P ,T,P,V ,L,Q,i i i x y z ,,)方程总数：2C+3个 需规定变量数：C+5个其中进料变量数：C+3个（F, F T ,F P ,i z ）根据其余2个变量的规定方法可将闪蒸计算分为如下五类：11=∑=Ci ix11=∑=Ci iy,...C,i Vy Lx Fz i i i 21 =+=Ci x K y i i i ,...2,1 ==LV F LH VH Q FH +=+表2-4闪蒸计算类型2.3.1 等温闪蒸规定：p 、T计算：Q, V , L,i i x y ,一、汽液平衡常数与组成无关 ()P T f K i ,=已知闪蒸温度和压力，i K 值容易确定，所以联立求解上述（2C+3）个方程比较简单。

具体步骤如下： 1. 输出变量求解将E---方程：代入M —方程： 消去i y ，得到： 将L=F-V 代入上式：汽化率代入（2-66）式，得到：Ci VK V F Fz x iii ,...2,1 =+-=(2-66))1(1-+=i ii K z x ψ（2-67） Ci x K y i i i ,...2,1 ==,...C ,i Vy Lx Fz i i i 21 =+=C i x VK Lx Fz i i i i ,...2,1 =+=FV /=ψ将（2-67）和（2-68）式代入S---方程，得到：两式相减，得：0)1(1)1()(=-+-=∑i ii K z K f ψψ--------------------------闪蒸方程0)1(1)1()(=-+-=∑i ii K z K f ψψ （2-71）)1(1-+=i ii K z x ψ i i i x K y = F=V+L L V F LH VH Q FH +=+通过闪蒸方程（2-71）求出汽化率ψ后，由（2-67）和（2-68）式可分别求出i i y x 和，进而由总物料衡算式（2-64）可求出V 和L,由热量衡算式（2-65）可求出Q汽化率ψ的迭代： 设ψ初值，计算:)(ψf可采用Newton-Raphson 法迭代ψ：（2-68）)1(1-+==i ii i i i K z K x K y ψ1)1(11=-+∑=Ci i iK z ψ（2-69）（2-70）1)1(11=-+∑=Ci i ii K z K ψ0)1(1)1()(1=-+-=∑=Ci i ii K z K f ψψ（2-71）2. Q 的计算L V F LH VH Q FH +=+Q-----吸热为正，移热为负H-----混合物的摩尔焓对于理想混合：3. 判断闪蒸过程是否可行的方法 方法一：已知T 、P对i Z 进行泡点计算：∑==-=Ci i i B Z K T f 101)( 试差泡点B T对i Z 进行露点计算：∑=⎪⎪⎭⎫⎝⎛=-=Ci i i D K Z T f 101)( 试差露点D T 判断：若D B T T T 闪蒸问题成立方法二：对i Z 在T 、P 下进行露点计算：对i Z 在T 、P 下进行泡点计算：—i Ci P T Li iL Ci P T Vi iV H Hx H Hy H ∑∑====1),(1),(纯组分摩尔焓判断：若 同时成立，闪蒸问题有解闪蒸过程计算框图：开始打印 BD BT T T T --=ψ输入T,P,F,Z ()()∑-+-=)1(11i i i k k Z f ψψ计算计算泡点B T []打印，结束−→−<YF εψ)(计算露点D T []22)1(1)1()('-+-∑-=i i i k k Z f ψψ)(')(1ψψψψf f k k -=+汽液平衡常数与组成有关的闪蒸计算 对i i y x ,,ψ分层迭代：开始给定F,Z,P,T估计初值x,y ψ计算()i i i i y x P T k K ,,,=),(p T F k i =打印过冷液体−→−>YB T T 过热蒸汽−→−<YD T T 由（2-67），（2-68）计算x,y 归一化i i y x ,比较：估计和归一化值 比较：k k ψψ和)1(+如果不直接迭代，重新估计x,y 值 由Rachford-Rice 方程迭代()1+k ψ思考题1、相平衡关系可用几种方法来表达。

闪蒸过程计算范文闪蒸过程是一种常见的汽化方式，其原理是通过降低液体的压力，使其在常温下迅速汽化。

闪蒸过程通常用于蒸发器、炉窑等设备中，可以通过瞬间释放液体中的挥发性成分，以提高产品的质量和产量。

下面将详细介绍闪蒸过程的计算方法。

在闪蒸过程中，液体从高压状态变成低压状态，即蒸发而不烧开。

首先，我们需要确定液体的性质，包括饱和蒸汽压、饱和温度等参数。

这些参数可以通过查阅相关的物性数据手册或者进行实验测定获得。

闪蒸过程的计算通常可以分为以下几个步骤：1.确定闪蒸前的液体状态：计算液体的焓和焓差。

液体的焓可以通过查表或者使用热力学软件计算得到。

焓差可以通过闪蒸前和闪蒸后的饱和蒸汽压和温度的差值计算得到。

2.确定闪蒸后的蒸汽质量：根据闪蒸前的液体状态和闪蒸后的压力条件，使用蒸汽表或者蒸发焓公式计算得到闪蒸后的蒸汽的温度和质量。

蒸汽质量也可以通过查表或者使用计算软件进行计算。

3.确定闪蒸过程中的液体和蒸汽的流量：根据设备的流量要求和液体的属性，可以通过质量守恒方程和能量守恒方程计算得到闪蒸过程中的液体流量和蒸汽流量。

液体流量可以通过闪蒸后的液体密度和闪蒸后的液体速度计算得到。

蒸汽流量可以通过闪蒸后的蒸汽密度和蒸汽速度计算得到。

4.确定闪蒸后的冷却效果：闪蒸后生成的蒸汽需要通过冷凝器进行冷却处理。

冷却效果取决于冷凝器的设计参数和冷却介质的属性。

可以通过热传导原理进行计算，确定冷却传热的速率和效果。

5.对闪蒸过程进行总体性能评估：可以通过计算液体和蒸汽的能量平衡，计算设备的热效率和质量效率，并分析各种因素对闪蒸过程的影响。

需要注意的是，闪蒸过程中液体的压力骤降会引起液体的蒸发，液体的温度也会随之降低。

在实际操作中，需要考虑设备的材料和结构强度，避免由于压力骤降引起的瞬间膨胀和爆炸事故的发生。

此外，在进行闪蒸过程计算时，还需要考虑实际操作中的一些细节，如设备的压力损失、管道的阻力、流体的黏性等因素。

这些因素会对闪蒸过程的计算结果产生一定的影响，因此在实际操作中需要进行精确的计算和合理的安全考虑。

专业班级化工三班学号09040403 指导教师刘瑾姓名实验日期化工分离工程上机题目一：等温闪蒸过程计算：1.闪蒸过程以及流程闪蒸是连续单级蒸馏过程。

该过程使进料混合物部分汽化或冷凝得到含易挥发组分较多的蒸汽和含难挥发组分较多的液体。

2.过程的数学描述（2C+3个方程）3.计算的分类4. 计算框图：1. 试问乙酸甲酯（1）0.33、丙酮（2）0.34和甲醇（3）0.33在50℃和85.46kPa 条件下是否存在闪蒸现象。

如果存在闪蒸，计算汽液组成和汽化分率。

已知参数：s s3),,,(y x P T f K i∑∑∑===ΛΛ-Λ-=ck cj kjj kjk c j ij j i x x x 111)(ln 1ln γ例如：⎥⎦⎤⎢⎣⎡+Λ+Λ+Λ++Λ+Λ+Λ+-Λ+Λ+-=33223113233221121331221113312211)ln(1ln x x x x x x x x x x x x x x x γ1.闪蒸问题核实：用安托尼公式计算P i s用wilson 公式计算γi用 计算K i用 1>∑iiK z 判断闪蒸是否成立。

2.ψ的迭代公式：∑=-+--=ci ii i K z K d df 122)]1(1[)1()(ψψψ3. 计算x i 、y i :4.归一化x,y ： ∑==ci iii xx x 1∑==ci iii yy y 1核算温度Option Base 1Dim a(), b(), c(), v(3,3), l(3, 3), xy()Dim x(), y(), z(), k(), fq(), dq(), pis(), f, dPP K si i i γ=1>∑i i z K )()()()()1(/)()(k k k k k ΨΨf Ψf ΨΨd d -=+)1(1-+=i ii K Ψz x )1(1-+=i ii i K Ψz K yPrivate Sub Command1_Click()a(1) = 16.1295: a(2) = 16.6513: a(3) = 18.5875b(1) = 2601.92: b(2) = 2940.46: b(3) = 3626.55c(1) = -56.15: c(2) = -35.93: c(3) = -34.29z(1) = 0.33: z(2) = 0.34: z(3) = 0.33t = 323: p = 641v(1, 1) = 1: v(1, 2) = 1.1816: v(1, 3) = 0.52297: v(2, 1) = 0.71891: v(2, 2) = 1v(2, 3) = 0.50878: v(3, 1) = 0.57939: v(3, 2) = 0.97531: v(3, 3) = 1For i = 1 To 3pis(i) = Exp(a(i) - b(i) / (t + c(i)))l(i, 1) = (z(1) * v(1, i) / (z(1) + z(2) * v(1, 2) + z(3) * v(1, 3)))l(i, 2) = z(2) * v(2, i) / (z(1) * v(2, 1) + z(2) + z(3) * v(2, 3))l(i, 3) = z(3) * v(3, i) / (z(1) * v(3, 1) + z(2) * v(3, 2) + z(3))y(i) = Exp(1 - Log(z(1) * v(i, 1) + z(2) * v(i, 2) + z(3) * v(i, 3)) - (l(i, 1) + l(i, 2) + l(i, 3))) k(i) = y(i) * pis(i) / pNextIf z(1) * k(1) + z(2) * k(2) + z(3) * k(3) >1 and z(1) / k(1) + z(2) / k(2) + z(3) / k(3) >1 then Print “存在闪蒸现象，可继续计算”q = 0.5: o = 1While (q - o) / o >= 0.0001 Or (q - o) / o <= -0.0001For j = 1 To 3fq(j) = (k(j) - 1) * z(j) / (1 + q * (k(j) - 1))dq(j) = -((k(j) - 1) ^ 2 * z(j) / (1 + q * (k(j) - 1)) ^ 2)Nextf = fq(1) + fq(2) + fq(3)d = dq(1) + dq(2) + dq(3)o = qq = o - f / dWendFor h = 1 To 3x(h) = z(h) / (1 + q * (k(h) - 1))xy(h) = k(h) * z(h) / (1 + q * (k(h) - 1))NextFor m = 1 To 3x(m) = x(m) / (x(1) + x(2) + x(3))xy(m) = xy(m) / (xy(1) + xy(2) + xy(3))Nextn = (z(1) - x(1)) / (y(1) - x(1))Print”气化率分率为：”;nPrint "液相组成：乙酸甲酯："; x(1) ，"丙酮"; x(2) ，"甲醇：";x(3)Print "气相组成：乙酸甲酯:"; xy(1)，"丙酮:"; xy(2)，"甲醇:"; xy(3)ElsePrint “不存在闪蒸’EndEnd ifEnd Sub第二题：核实闪蒸温度Option Base 1Dim a(), b(), c(), xy()Dim x(), y(), z(), k(), fq(), dq(), pis(), f, dPrivate Sub Command1_Click()a(1) = 20.7936: a(2) = 20.7665b(1) = 2788.15: b(2) = 2911.32c(1) = -52.36: c(2) = -56.51z(1) = 0.5: z(2) =0.5t =333.15: p = 42500y(1) = Exp(0.315559 / (1 + 0.315559 * z(1) / 0.402295 / z(2)) ^ 2)y(2) = Exp(0.402295 / (1 + 0.402295 * z(2) / 0.315559 / z(1)) ^ 2)For i = 1 To 2k(i) = y(i) * pis(i) / pNext计算组成及汽化分率If z(1) * k(1) + z(2) * k(2) + z(3) * k(3) >1 and z(1) / k(1) + z(2) / k(2) + z(3) / k(3) >1 then Print “存在闪蒸，可继续计算”q = 0.5: o = 1While (q - o) / o >= 0.0001 Or (q - o) / o <= -0.0001For j = 1 To 2fq(j) = (k(j) - 1) * z(j) / (1 + q * (k(j) - 1))dq(j) = -((k(j) - 1) ^ 2 * z(j) / (1 + q * (k(j) - 1)) ^ 2)Nextf = fq(1) + fq(2)d = dq(1) + dq(2)o = qq = o - f / dWendFor h = 1 To 2x(h) = z(h) / (1 + q * (k(h) - 1))xy(h) = k(h) * z(h) / (1 + q * (k(h) - 1)) NextFor m = 1 To 2x(m) = x(m) / (x(1) + x(2))xy(m) = xy(m) / (xy(1) + xy(2))Nextn = (z(1) - x(1)) / (y(1) - x(1))Print "x1="; x(1): Print "x2="; x(2) ElsePrint “不存在闪蒸现象，不可再计算’EndEnd ifEnd Sub例2-5Dim x(1 To 3) As SingleDim p As SingleDim t As SinglePublic Function ln(q As Single) As Singleln = Log(q) / 0.434End FunctionPrivate Sub Command1_Click()For i = 1 To 3x(i) = InputBox("请输入" & i & "组分的液相摩尔含量", "数据输入")Next iEnd SubPrivate Sub Command2_Click()p = InputBox("混合气体总压", "数据输入")End SubPrivate Sub Command4_Click()t = InputBox("输入温度初始值（开尔文温度）", "数据输入")End SubPrivate Sub Command3_Click()Dim vml(1 To 3) As Singlevml(1) = 73.52: vml(2) = 89.57: vml(3) = 97.79Dim y1 As Single, y2 As SingleDoDim yd(1 To 3) As SingleFor i2 = 1 To 3yd(i2) = 1Next i2Dim y(1 To 3, 1 To 3) As SingleDim ym(1 To 3, 1 To 3) As Singley(1, 2) = 5741.401: y(2, 1) = -2722.056: y(1, 3) = -1226.628: y(3, 1) = 2698.313: y(2, 3) = -1696.533: y(3, 2) = 11322.895For i = 1 To 3For j = 1 To 3If i <> j Thenym(i, j) = (vml(j) / vml(i)) * Exp(-y(i, j) / (8.314 * t))End IfIf i = j Thenym(i, j) = 0End IfNext jNext iDim hd(1 To 3) As Singlehd(1) = Exp(1 - ln(ym(1, 2) * x(2) + ym(1, 3) * x(3)) - (ym(2, 1) * x(2) / (ym(2, 1) * x(1) + ym(2, 3) * x(3)) + ym(3, 1) * x(3) / (ym(3, 1) * x(1) + ym(3, 2) * x(2))))hd(2) = Exp(1 - ln(ym(2, 1) * x(1) + ym(2, 3) * x(3)) - (ym(1, 2) * x(1) / (ym(1, 2) * x(2) + ym(1, 3) * x(3)) + ym(3, 2) * x(3) / (ym(3, 1) * x(1) + ym(3, 2) * x(2))))hd(3) = Exp(1 - ln(ym(3, 1) * x(1) + ym(3, 2) * x(2)) - (ym(1, 3) * x(1) / (ym(1, 2) *x(2) + ym(1, 3) * x(3)) + ym(2, 3) * x(2) / (ym(2, 1) * x(1) + ym(2, 3) * x(3))))Dim pscs(1 To 3, 1 To 3) As SingleDim ps(1 To 3) As Singlepscs(1, 1) = 14.6363: pscs(1, 2) = 2940.46: pscs(1, 3) = -35.93pscs(2, 1) = 14.5836: pscs(2, 2) = 3150.42: pscs(2, 3) = -36.55pscs(3, 1) = 14.1366: pscs(3, 2) = 2790.5: pscs(3, 3) = -57.15For m = 1 To 3ps(m) = Exp(pscs(m, 1) - pscs(m, 2) / (pscs(m, 3) + t))Next mDim tc(1 To 3) As Single, tr(1 To 3) As SingleDim pc(1 To 3) As Single, pr(1 To 3) As SingleDim w(1 To 3) As Singletc(1) = 508.1: tc(2) = 535.6: tc(3) = 523.25pc(1) = 4701.5: pc(2) = 4154.33: pc(3) = 3830.09w(1) = 0.309: w(2) = 0.329: w(3) = 0.363Dim b(1 To 3) As Single, b0(1 To 3) As Single, b1(1 To 3) As SingleDim zd(1 To 3) As SingleDim yds(1 To 3) As SingleFor n = 1 To 3tr(n) = t / tc(n)b0(n) = 0.083 - 0.422 / ((tr(n)) ^ 1.6)b1(n) = 0.139 - 0.172 / ((tr(n)) ^ 4.2)b(n) = 8.314 * tc(n) * (b0(n) + w(n) * b1(n)) / pc(n)zd(n) = 1 + p * b(n) / (8.314 * t)yds(n) = Exp(2 * b(n) * ps(n) / (8.314 * t) - ln(zd(n)))Next nDim k(1 To 3) As SingleDim my(1 To 3) As Singley1 = 0For i1 = 1 To 3k(i1) = hd(i1) * ps(i1) * yds(i1) * Exp(vml(i1) * (p - ps(i1)) / (8.314 * t)) / (p * yd(i1))my(i1) = k(i1) * x(i1)y1 = y1 + my(i1)Next i1For i5 = 1 To 3my(i5) = my(i5) / y1yd(i5) = my(i5) * yds(i5)Next i5Do例2-9Option Base 1Dim a(3, 3), x(3), pis(3), vil(3), xi(3)Dim r, t, p1, p2, f1Private Sub Command1_Click()a(1, 1) = 1: a(2, 1) = 0.71891: a(3, 1) = 0.57939a(1, 2) = 1.1816: a(2, 2) = 1: a(3, 2) = 0.97513a(1, 3) = 0.52297: a(2, 3) = 0.50878: a(3, 3) = 1pis(1) = 78.049: pis(2) = 81.818: pis(3) = 55.581vil(1) = 83.77: vil(2) = 76.81: vil(3) = 42.05r = 8.314: t = 323.16x(1) = 0.33: x(2) = 0.34: x(3) = 0.33For i = 1 To 3xi(i) = x(i)Nextf1 = 0.5Doffirst = f1xi(1) = xxx(1)xi(2) = xxx(2)xi(3) = xxx(3)fend = f2()If Abs(ffirst - fend) < 0.001 ThenExit DoElsef1 = fendLoopx(1) = xi(1)x(2) = xi(2)x(3) = xi(3)Print x(1), x(2), x(3), fendEnd SubPublic Function ln(x)ln = Log(x) / Log(2.71828)End FunctionPublic Function gama(i)gama = Exp(1 - ln(ai(i)) - (xi(1) * a(1, i) / ai(1) + xi(2) * a(2, i) / ai(2) + xi(3) * a(3, i) / ai(3)))Print xi(1)End FunctionPublic Function ai(i)ai = xi(1) * a(i, 1) + xi(2) * a(i, 2) + xi(3) * a(i, 3)End FunctionPublic Function ki(i)ki = gama(i) * pis(i) * Exp(vil(i) * (p1 - pis(i)) / 1000 / r / t) / p()End FunctionPublic Function p()p = pis(1) * x(1) + pis(2) * x(2) + pis(3) * x(3)End FunctionPublic Function fai()fai = (ki(1) - 1) * xi(1) / (1 + f1 * (ki(1) - 1)) + (ki(2) - 1) * xi(2) / (1 + f1 * (ki(2) - 1)) + _(ki(3) - 1) * xi(3) / (1 + f1 * (ki(3) - 1))End FunctionPublic Function df1df()df1df = -((ki(1) - 1) ^ 2 * xi(1) / (1 + f1 * (ki(1) - 1)) + ki(2) - 1) ^ 2 * xi(2) / (1 + f1 * (ki(2) - 1))(ki(3) - 1) ^ 2 * xi(3) / (1 + f1 * (ki(3) - 1))End FunctionPublic Function f2()f2 = f1 - fai() / df1df()End FunctionPublic Function xxx(i)xxx = xi(i) / (1 + f1 * (ki(i) - 1))End Function。

第三节 闪蒸过程的计算2.3 等温闪蒸和部分冷凝过程流程示意图：闪蒸过程的计算方程（MESH ） ⑴物料衡算----M 方程： C 个⑵相平衡--------E 方程： C 个⑶摩尔分率加和式---S 方程： 2个⑷热量平衡式-------H 方程： 1个变量数：3C+8个 (F, F T ,F P ,T,P,V ,L,Q,i i i x y z ,,)方程总数：2C+3个 需规定变量数：C+5个其中进料变量数：C+3个（F, F T ,F P ,i z ）根据其余2个变量的规定方法可将闪蒸计算分为如下五类：11=∑=Ci ix11=∑=Ci iy,...C,i Vy Lx Fz i i i 21 =+=Ci x K y i i i ,...2,1 ==LV F LH VH Q FH +=+表2-4闪蒸计算类型2.3.1 等温闪蒸规定：p 、T计算：Q, V , L,i i x y ,一、汽液平衡常数与组成无关 ()P T f K i ,=已知闪蒸温度和压力，i K 值容易确定，所以联立求解上述（2C+3）个方程比较简单。

具体步骤如下： 1. 输出变量求解将E---方程：代入M —方程： 消去i y ，得到： 将L=F-V 代入上式：汽化率代入（2-66）式，得到：Ci VK V F Fz x iii ,...2,1 =+-=(2-66))1(1-+=i ii K z x ψ（2-67） Ci x K y i i i ,...2,1 ==,...C ,i Vy Lx Fz i i i 21 =+=C i x VK Lx Fz i i i i ,...2,1 =+=FV /=ψ将（2-67）和（2-68）式代入S---方程，得到：两式相减，得：0)1(1)1()(=-+-=∑i ii K z K f ψψ--------------------------闪蒸方程0)1(1)1()(=-+-=∑i ii K z K f ψψ （2-71）)1(1-+=i ii K z x ψ i i i x K y = F=V+L L V F LH VH Q FH +=+通过闪蒸方程（2-71）求出汽化率ψ后，由（2-67）和（2-68）式可分别求出i i y x 和，进而由总物料衡算式（2-64）可求出V 和L,由热量衡算式（2-65）可求出Q汽化率ψ的迭代： 设ψ初值，计算:)(ψf可采用Newton-Raphson 法迭代ψ：（2-68）)1(1-+==i ii i i i K z K x K y ψ1)1(11=-+∑=Ci i iK z ψ（2-69）（2-70）1)1(11=-+∑=Ci i ii K z K ψ0)1(1)1()(1=-+-=∑=Ci i ii K z K f ψψ（2-71）2. Q 的计算L V F LH VH Q FH +=+Q-----吸热为正，移热为负H-----混合物的摩尔焓对于理想混合：3. 判断闪蒸过程是否可行的方法 方法一：已知T 、P对i Z 进行泡点计算：∑==-=Ci i i B Z K T f 101)( 试差泡点B T对i Z 进行露点计算：∑=⎪⎪⎭⎫⎝⎛=-=Ci i i D K Z T f 101)( 试差露点D T 判断：若D B T T T 闪蒸问题成立方法二：对i Z 在T 、P 下进行露点计算：对i Z 在T 、P 下进行泡点计算：—i Ci P T Li iL Ci P T Vi iV H Hx H Hy H ∑∑====1),(1),(纯组分摩尔焓判断：若 同时成立，闪蒸问题有解闪蒸过程计算框图：开始打印 BD BT T T T --=ψ输入T,P,F,Z ()()∑-+-=)1(11i i i k k Z f ψψ计算计算泡点B T []打印，结束−→−<YF εψ)(计算露点D T []22)1(1)1()('-+-∑-=i i i k k Z f ψψ)(')(1ψψψψf f k k -=+汽液平衡常数与组成有关的闪蒸计算 对i i y x ,,ψ分层迭代：开始给定F,Z,P,T估计初值x,y ψ计算()i i i i y x P T k K ,,,=),(p T F k i =打印过冷液体−→−>YB T T 过热蒸汽−→−<YD T T 由（2-67），（2-68）计算x,y 归一化i i y x ,比较：估计和归一化值 比较：k k ψψ和)1(+如果不直接迭代，重新估计x,y 值 由Rachford-Rice 方程迭代()1+k ψ思考题1、相平衡关系可用几种方法来表达。