比较最近两次计算的∑y是否相等,如果不相等,则归一化所有计算得到的y,并 以此值、设定温度T和已知压力P共同计算φ",返回V;如果相等,进入(ⅤI V考察∑y是否等于1?若∑y≠1,则重新给定温度r值,返回(I;若∑y,=1, 则计算结束,最新给定的温度就是相平衡温度T,最后一次归一化的y即为汽相组成 另外,如果∑y>1,重新设定温度T时,应该有所降低,相反,∑y<1时,应提高温 开始 停止 输入P及各x 打印结果T,y 假定T令各”=1 否 计算P,,1 调整T ∑y 否 计算全部y 所有y归一化 ∑y是 计算全部ψ 否改变 是否 计算∑y 否 第一次迭代 图5-9泡点温度组成计算框图 3)露点压力和组成计算(DEWP) 已知:平衡温度T,液相组成y,y… 求:平衡压力P,汽相组成x,x2…xN 将式(5-14)做恒等变形,得: pyi g (i=1.2…N) pr (520) 由式(5-16)~(5-19)的分析可知,已知温度T就直接可得到p、g的值,若假定 了P,联合已知的温度T和液相组成,可计算¢的值。那么,等式两边就是关于x的迭
11 VI. 比较最近两次计算的 ∑ i y 是否相等,如果不相等,则归一化所有计算得到的 yi,并 以此值、设定温度 T 和已知压力 p 共同计算 v φi ˆ ,返回 V;如果相等,进入(VII); VII. 考察∑ i y 是否等于 1?若∑ ≠ 1 i y ,则重新给定温度 T 值,返回(II);若∑ = 1 i y , 则计算结束,最新给定的温度就是相平衡温度 T,最后一次归一化的 yi 即为汽相组成。 另外,如果∑ > 1 i y ,重新设定温度 T 时,应该有所降低,相反,∑ < 1 i y 时,应提高温 度。 图 5-9 泡点温度组成计算框图 3) 露点压力和组成计算(DEWP) 已知:平衡温度 T,液相组成 y1,y2 ···yN 求 :平衡压力 p,汽相组成 x1,x2 ···xN 将式(5-14)做恒等变形,得: i s i s i v i i i p py x φ γ φ ˆ = (i = 1, 2 … N) (5-20) 由式(5-16)~(5-19)的分析可知,已知温度 T 就直接可得到 s pi 、 s φi 的值,若假定 了 p,联合已知的温度 T 和液相组成 yi,可计算 v φi ˆ 的值。那么,等式两边就是关于 xi 的迭
代。计算框图如图5-10所示。具体的计算步骤为: 1.假设P;因为y与x有关,第一次试算时,先令所有组分的y1=1 I1用 Antonie方程计算p';选择适当的EOS计算g"; I.选择适当的EOS计算"; Ⅳ.用式(520)计算所有的x,进而计算∑x 若是第一次迭代,则归一化所有计算得到的x,并以此值和已知的温度T共同计算y; 若不是第一次迭代,跳到(VID; Ⅵ将最新的x和、p'、、已知的汽相组成y、设定的压力P,代入式(5-20),计 算所有的x,进而计算∑x; VI.返回(V) VI比较最近两次计算的∑x是否相等,如果不相等,则归一化最近一次计算得到的x, 并以此值和已知的温度T共同计算y,返回(Ⅵ);如果相等,进入(IX); IX考察∑x是否等于1?若∑x≠1,则重新给定平衡压力P值,返回(Ⅲ):若 ∑x,=1,则计算结束,最新给定的压力就是相平衡压力P,最后一次归一化的x即 为液相组成 另外,如果∑x>1,压力P的设定值应该有所降低,相反,∑x1<1时,应提高压力设 定值
12 代。计算框图如图 5-10 所示。具体的计算步骤为: I. 假设 p;因为 i γ 与 xi 有关,第一次试算时,先令所有组分的γ i = 1; II. 用 Antonie 方程计算 s pi ;选择适当的 EOS 计算 s φi ; III. 选择适当的 EOS 计算 v φi ˆ ; IV. 用式(5-20)计算所有的 xi,进而计算∑ i x ; V. 若是第一次迭代,则归一化所有计算得到的 xi,并以此值和已知的温度 T 共同计算 i γ ; 若不是第一次迭代,跳到(VIII); VI. 将最新的 i γ 和 v φi ˆ 、 s pi 、 s φi 、已知的汽相组成 yi、设定的压力 p,代入式(5-20),计 算所有的 xi,进而计算∑ i x ; VII. 返回(V); VIII. 比较最近两次计算的∑ i x 是否相等,如果不相等,则归一化最近一次计算得到的 xi, 并以此值和已知的温度 T 共同计算 i γ ,返回(VI);如果相等,进入(IX); IX. 考察 ∑ i x 是否等于 1?若 ∑ ≠ 1 i x ,则重新给定平衡压力 p 值,返回(III);若 ∑ = 1 i x ,则计算结束,最新给定的压力就是相平衡压力 p,最后一次归一化的 xi 即 为液相组成。 另外,如果∑ > 1 i x ,压力 p 的设定值应该有所降低,相反, ∑ < 1 i x 时,应提高压力设 定值
开始 停止 输入T及各y 设P,令所有y=1 打印结果P 计算p,剪 是 计算全部 调整P ∑x 否 所有x归化 计算全部x 计算全部y 否变化 是否 计算Σx 次迭 图5-10露点压力组成计算框图 4)露点温度和组成计算(DEWT) 已知:平衡压力P,汽相组成y,y2…W 求:平衡温度T,液相组成x,x2…x 这类计算与前面的计算大同小异,这里只给出具体的计算框图,如图5-11所示,读者 可自行给出计算步骤,并与框图对照比较
13 图 5-10 露点压力组成计算框图 4) 露点温度和组成计算(DEWT) 已知:平衡压力 p,汽相组成 y1,y2 ···yN 求 :平衡温度 T,液相组成 x1,x2 ···xN 这类计算与前面的计算大同小异,这里只给出具体的计算框图,如图 5-11 所示,读者 可自行给出计算步骤,并与框图对照比较