Simulation of a steady state flash

Transcription

1 Smulaton of a steady state flash Descrpton: Statonary flash smulaton of an Ethanol(1) - Water(2) - mxture Wth followng assumptons: Apart from heater and mass flows, no energy s transferred across the system boundary. Lqud and vapor phase are each deally mxed. Lqud and vapor phase are n equlbrum. The vapor phase behaves as an deal gas and the lqud as an ncompressble flud. VLE descrbed by extended Raoult s law. Vapor pressure calculaton by Antone. Evaporaton enthalpy calculaton by PPDS12. Real flud behavor of the lqud phase s expressed through Wlson s ge model. Excess varables neglected. Wlson and Antone parameter are taken from: Gmehlng, Kolbe; PPDS12 parameter are taken from: Equaton System: 86383: statflash.moseqs IndexSpecfcaton: e[0]86383.nc = 2 Varable Specfcaton: 86386: IV statflash.mosvar Parameter Specfcaton: none. Results Specfcaton: 86388: RE statflash CBzzMath Org.mosvar Herarchcal vew of equatons: Equaton System: 86383: statflash.moseqs Descrpton: Statonary flash wth followng addtonal assumpton to bascflash: 1. Apart from heater and mass flows, no energy s transferred across the system boundary. 2. flash geometry s cylndrc. 1

2 Connected Equatons: Eq: 86378: StatEnergyBalance.mosequ (usng Nota: 8639: not - Desc.: Energy balance (statonary) 0 = F F,n h F,n h B,n F B,n h D,n F D,n + Q (1) Eq: 86379: StatComponentBalance.mosequ (usng Nota: 8639: not Desc.: Component balance (statonary) 0 = x F F F,n x B F B,n y D F D,n (2) Connected EQ-Systems: 86375: bascflash.moseqs 8638: tankgeometry.moseqs Connecton Level (1) EQ-Systems connected to 86383: statflash.moseqs: Equaton System: 86375: bascflash.moseqs Descrpton: basc model of a flash wth followng assumptons: 1. lqud and vapor phase are each deally mxed. Lqud and vapor phase are n equlbrum. 2. The vapor phase behaves as an deal gas and the lqud as an ncompressble flud. 3. VLE descrpton by extended Raoult s law. Vapor pressure calculaton by Antone 5. Actvty calculatons n the lqud phase by Wlson s g E -Model Connected Equatons: 2

3 Eq: 86350: sumrelatonvapor.mosequ (usng Nota: 8639: not - Desc.: Molar summaton relaton (vapor) NC 1 = =1 y D (3) Eq: 86353: molarenthalpyvapor.mosequ (usng Nota: 8639: not - Desc.: Molar enthalpy (vapor) NC h D,n = =1 y D h D,n + h E,D,n () Eq: 86356: molarenthalpyfeed.mosequ (usng Nota: 8639: not - Desc.: Molar enthalpy (feed) NC h F,n = =1 x F h F,n + h E,F,n (5) Eq: 8636: molarevapenthalpylqud.mosequ (usng Nota: 8639: not Desc.: Molar evaporaton enthalpy (lqud, pure component), source: h B,LV,n = R T c (A P P DS12, (1 T T c ) BP P DS12, (1 T T c + C P P DS12, (1 T T c ) + D P P DS12, (1 T T c ) 2 3 ) 2 + E P P DS12, (1 T T c ) 6 ) Eq: 86368: molarenthalpycomponentvapor.mosequ (usng Nota: 8639: not Desc.: Molar enthalpy (vapor, pure component), assumptons: deal gas (6) h D,n = h o,n + C cp, 3 + A cp, (T T o ) + B cp, 2 ((T ) 3 (T o ) 3 ) + D cp, ((T ) 2 (T o ) 2 ) ((T ) (T o ) ) (7) 3

4 Eq: 86367: molarenthalpycomponentfeed.mosequ (usng Nota: 8639: not Desc.: Molar enthalpy (feed, pure component), assumptons: ncompressble flud deal gas behavor of correspondng vapor phase (pressure ndependency of h n V-phase) h F,n = h o,n + A cp, (T F T o ) + B cp, ((T F ) 2 (T o ) 2 ) 2 + C cp, ((T F ) 3 (T o ) 3 ) + D cp, 3 ((T F ) (T o ) ) + v L,n (p F p F,LV ) h F,LV,n (8) Eq: 8637: actvtycoeffwlsonparameterlqud.mosequ (usng Nota: 8639: not Desc.: Parameter calculaton (lqud) based on Wlson s g E -model (lqud), source: Gmehlng, Kolbe S.239 α B = NC =1 vl,n v L,n v L,n exp( λ T ) (9) Eq: 86372: volumeflowlqud.mosequ (usng Nota: 8639: not - Desc.: Volume flow (lqud) F B,n = F B,v ρ B,n (10) Eq: 86363: pressuredrop.mosequ (usng Nota: 8639: not - Desc.: Pressure drop (Feed > flash) p = p F p (11) Eq: 86361: volumelqud.mosequ (usng Nota: 8639: not - Desc.: Volume n flash (lqud) NC HU L,v = ( v L,n x B + v L,n,E ) HU L,n 1000 =1 (12)

5 Eq: 86351: sumrelatonlqud.mosequ (usng Nota: 8639: not - Desc.: Molar summaton relaton (lqud) NC 1 = =1 x B (13) Eq: 86352: molarenthalpylqud.mosequ (usng Nota: 8639: not - Desc.: Molar enthalpy (lqud) NC h B,n = =1 x B h B,n + h B,n,E (1) Eq: 86359: molarholdupcomponent.mosequ (usng Nota: 8639: not Desc.: Molar component holdup wthn flash HU n = x B HU L,n + y D HU V,n (15) Eq: 86357: levellqud.mosequ (usng Nota: 8639: not Desc.: Lqud level calculaton L L = HU L,v π (d) 2 (16) Eq: 86369: vaporpressureantonelqud.mosequ (usng Nota: 8639: not Desc.: Vapor pressure by Antone equaton (lqud, pure component), source: Kolbe, Mehlng p.59 p LV,B = (10) 3 (10) A Antone, B Antone, C Antone, +(T ) (17) Eq: 86371: molardenstylqud.mosequ (usng Nota: 8639: not - Desc.: Molar densty wth no excess volume (lqud) NC 1 = ρ B,n =1 x B v L,n (18) 5

6 Eq: 86370: vaporpressureantonefeed.mosequ (usng Nota: 8639: not Desc.: Vapor pressure by Antone equaton (feed, pure component), source: Kolbe, Mehlng p.59 p LV,F = (10) 3 (10) A Antone, B Antone, C Antone, +(T F ) (19) Eq: 86362: volume.mosequ (usng Nota: 8639: not Desc.: Volume of flash HU v = HU L,v + HU V,v (20) Eq: 86365: molarevapenthalpyfeed.mosequ (usng Nota: 8639: not Desc.: Molar evaporaton enthalpy (feed, pure component), source: h F,LV,n = R T c (A P P DS12, (1 T F ) BP P DS12, (1 T F T c T c ) CP P DS12, (1 T F ) + D P P DS12, (1 T F T c T c ) 2 + E P P DS12, (1 T F ) 6 ) T c (21) U Eq: 86358: nternalenergy.mosequ (usng Nota: 8639: not - Desc.: Internal energy wthn flash (22) = HU L,n (h B,n p ( NC =1 xb v L,n + v L,n,E )) + HU V,n (h D,n R T z) (10) 6 Eq: 86360: volumevapor.mosequ (usng Nota: 8639: not - Desc.: Volume n flash (vapor) HU V,v = HU V,n R T z p 1000 (23) 6

7 Eq: 8635: sumrelatonfeed.mosequ (usng Nota: 8639: not - Desc.: Molar summaton relaton (feed) NC 1 = =1 x F (2) Eq: 86373: actvtycoeffwlsonlqud.mosequ (usng Nota: 8639: not Desc.: Actvty coeffcent calculatons by Wlson s g E -model (lqud) γ B = 1 x B + α B (1 x B ) α B exp((1 x B ) ( x B + α B (1 x B ) NC =1 αb α B ( NC =1 αb α B) xb + (1 x B ))) (25) Eq: 86355: VLEextendedRaoult.mosequ (usng Nota: 8639: not - Desc.: Extended Raoult s law for VLE n flash y D = x B pb,lv p γ B (26) Eq: 86366: molarenthalpycomponentlqud.mosequ (usng Nota: 8639: not Desc.: Molar enthalpy (lqud, pure component), assumptons: ncompressble flud deal gas behavor of correspondng vapor phase (pressure ndependency of h n V-phase) h B,n = h o,n + A cp, (T T o ) + B cp, ((T ) 2 (T o ) 2 ) 2 + C cp, ((T ) 3 (T o ) 3 ) + D cp, 3 ((T ) (T o ) ) + v L,n (p p B,LV ) h B,LV,n (27) 7

8 Equaton System: 8638: tankgeometry.moseqs Descrpton: Tank geometry Connected Equatons: Eq: 86381: volumetank.mosequ (usng Nota: 8639: not - Desc.: Flash s volume HU v = A L (28) Eq: 86382: ratodameterlengthtank.mosequ (usng Nota: 8639: not Desc.: dameter to heght rato of tank r D,L = d L (29) Eq: 86380: crosssectonarea.mosequ (usng Nota: 8639: not - Desc.: Flash s cross secton area A = π (d)2 (30) Equaton nstances: Eq: 86378: StatEnergyBalance.mosequ (usng Nota: 8639: not. Descrpton: Energy balance (statonary). 0 = e0.f F,n e0.h F,n e0.h B,n e0.f B,n e0.h D,n e0.f D,n + e0.q (31) Eq: 86379: StatComponentBalance.mosequ (usng Nota: 8639: not. Descrpton: Component balance (statonary). 0 = e0.x F =1 e0.f F,n e0.x B =1 e0.f B,n e0.y D =1 e0.f D,n (32) 0 = e0.x F =2 e0.f F,n e0.x B =2 e0.f B,n e0.y D =2 e0.f D,n (33) Eq: 86350: sumrelatonvapor.mosequ (usng Nota: 8639: not. Descrpton: Molar summaton relaton (vapor). 1 = (e0.y D =1 + e0.y D =2) (3) 8

9 Eq: 86351: sumrelatonlqud.mosequ (usng Nota: 8639: not. Descrpton: Molar summaton relaton (lqud). 1 = (e0.x B =1 + e0.x B =2) (35) Eq: 8635: sumrelatonfeed.mosequ (usng Nota: 8639: not. Descrpton: Molar summaton relaton (feed). 1 = (e0.x F =1 + e0.x F =2) (36) Eq: 86355: VLEextendedRaoult.mosequ (usng Nota: 8639: not. Descrpton: Extended Raoult s law for VLE n flash. e0.y D =1 = e0.x B =1 e0.pb,lv =1 e0.γ B =1 e0.p (37) e0.y D =2 = e0.x B =2 e0.pb,lv =2 e0.γ B =2 e0.p (38) Eq: 86352: molarenthalpylqud.mosequ (usng Nota: 8639: not. Descrpton: Molar enthalpy (lqud). e0.h B,n = (e0.x B =1 e0.h B,n =1 + e0.xb =2 e0.h B,n =2 ) + e0.hb,e,n (39) Eq: 86353: molarenthalpyvapor.mosequ (usng Nota: 8639: not. Descrpton: Molar enthalpy (vapor). e0.h D,n = (e0.y D =1 e0.h D,n =1 + e0.yd =2 e0.h D,n =2 ) + e0.hd,e,n (0) Eq: 86356: molarenthalpyfeed.mosequ (usng Nota: 8639: not. Descrpton: Molar enthalpy (feed). e0.h F,n = (e0.x F =1 e0.h F,n =1 + e0.xf =2 e0.h F,n =2 ) + e0.he,f,n (1) Eq: 86357: levellqud.mosequ (usng Nota: 8639: not. Descrpton: Lqud level calculaton. e0.l L = e0.hu L,v e0.π (e0.d) (2) (2) Eq: 86358: nternalenergy.mosequ (usng Nota: 8639: not. Descrpton: Internal energy wthn flash. e0.u = e0.hu L,n (e0.h B,n e0.p ((e0.x B =1 e0.vl,n =1 + e0.xb =2 e0.vl,n =2 ) + e0.ve,l,n )) + e0.hu V,n (e0.h D,n (10) (6) Eq: 86359: molarholdupcomponent.mosequ (usng Nota: 8639: not. Descrpton: Molar component holdup wthn flash. (3) e0.hu n =1 = e0.x B =1 e0.hu L,n + e0.y D =1 e0.hu V,n () 9

10 e0.hu n =2 = e0.x B =2 e0.hu L,n + e0.y D =2 e0.hu V,n (5) Eq: 86360: volumevapor.mosequ (usng Nota: 8639: not. Descrpton: Volume n flash (vapor). e0.hu V,v = e0.hu V,n e0.r e0.t e0.z e0.p 1000 (6) Eq: 86361: volumelqud.mosequ (usng Nota: 8639: not. Descrpton: Volume n flash (lqud). e0.hu L,v = ((e0.v L,n =1 e0.xb =1 + e0.v L,n =2 e0.xb =2) + e0.v E,L,n ) e0.hu L,n 1000 (7) Eq: 86362: volume.mosequ (usng Nota: 8639: not. Descrpton: Volume of flash. e0.hu v = e0.hu L,v + e0.hu V,v (8) Eq: 86363: pressuredrop.mosequ (usng Nota: 8639: not. Descrpton: Pressure drop (Feed > flash). e0. p = e0.p F e0.p (9) Eq: 8636: molarevapenthalpylqud.mosequ (usng Nota: 8639: not. Descrpton: Molar evaporaton enthalpy (lqud, pure component), source: e0.h B,LV,n =1 = e0.r e0.t=1 c (e0.a P P DS12,=1 (1 e0.t e0.t=1 c ) ( 3) 1 + e0.b P P DS12,=1 (1 e0.t e0.t=1 c ) ( 2 3) + e0.cp P DS12,=1 (1 e0.t e0.t=1 c ) + e0.d P P DS12,=1 (1 e0.t e0.t=1 c ) (2) + e0.e P P DS12,=1 (1 e0.t e0.t=1 c ) (6) ) (50) e0.h B,LV,n =2 = e0.r e0.t=2 c (e0.a P P DS12,=2 (1 e0.t e0.t=2 c ) ( 3) 1 + e0.b P P DS12,=2 (1 e0.t e0.t=2 c ) ( 2 3) + e0.cp P DS12,=2 (1 e0.t e0.t=2 c ) + e0.d P P DS12,=2 (1 e0.t e0.t=2 c ) (2) + e0.e P P DS12,=2 (1 e0.t e0.t=2 c ) (6) ) (51) 10

11 Eq: 86365: molarevapenthalpyfeed.mosequ (usng Nota: 8639: not. Descrpton: Molar evaporaton enthalpy (feed, pure component), source: e0.h F,LV,n =1 = e0.r e0.t c =1 (e0.a P P DS12,=1 (1 e0.t F + e0.b P P DS12,=1 (1 e0.t F e0.t c =1 e0.t c =1 ) ( 1 3) ) ( 2 3) + e0.cp P DS12,=1 (1 e0.t F e0.t=1 c ) + e0.d P P DS12,=1 (1 e0.t F e0.t=1 c ) (2) + e0.e P P DS12,=1 (1 e0.t F ) (6) ) e0.t c =1 (52) e0.h F,LV,n =2 = e0.r e0.t c =2 (e0.a P P DS12,=2 (1 e0.t F + e0.b P P DS12,=2 (1 e0.t F e0.t c =2 e0.t c =2 ) ( 1 3) ) ( 2 3) + e0.cp P DS12,=2 (1 e0.t F e0.t=2 c ) + e0.d P P DS12,=2 (1 e0.t F e0.t=2 c ) (2) + e0.e P P DS12,=2 (1 e0.t F ) (6) ) e0.t c =2 (53) Eq: 86366: molarenthalpycomponentlqud.mosequ (usng Nota: 8639: not -. Descrpton: Molar enthalpy (lqud, pure component), assumptons:. ncompressble flud deal gas behavor of correspondng vapor phase (pressure ndependency of h n V-phase) e0.h B,n =1 = e0.hn,o =1 + e0.a cp,=1 (e0.t e0.t o =1) + e0.b cp,=1 2 ((e0.t ) (2) (e0.t o =1) (2) ) + e0.c cp,=1 3 + e0.d cp,=1 ((e0.t ) () (e0.t o =1) () ) + e0.v L,n =1 (e0.p e0.pb,lv =1 ) e0.h B,LV,n =1 ((e0.t ) (3) (e0.t o =1) (3) ) (5) 11

12 e0.h B,n =2 = e0.hn,o =2 + e0.a cp,=2 (e0.t e0.t o =2) + e0.b cp,=2 2 ((e0.t ) (2) (e0.t o =2) (2) ) + e0.c cp,=2 3 + e0.d cp,=2 ((e0.t ) () (e0.t o =2) () ) + e0.v L,n =2 (e0.p e0.pb,lv =2 ) e0.h B,LV,n =2 ((e0.t ) (3) (e0.t o =2) (3) ) (55) Eq: 86368: molarenthalpycomponentvapor.mosequ (usng Nota: 8639: not -. Descrpton: Molar enthalpy (vapor, pure component), assumptons:. deal gas e0.h D,n =1 = e0.hn,o =1 + e0.a cp,=1 (e0.t e0.t o =1) + e0.b cp,=1 2 ((e0.t ) (2) (e0.t o =1) (2) ) + e0.c cp,=1 3 + e0.d cp,=1 ((e0.t ) () (e0.t o =1) () ) e0.h D,n =2 = e0.hn,o =2 + e0.a cp,=2 (e0.t e0.t o =2) + e0.b cp,=2 2 ((e0.t ) (2) (e0.t o =2) (2) ) + e0.c cp,=2 3 + e0.d cp,=2 ((e0.t ) () (e0.t o =2) () ) ((e0.t ) (3) (e0.t o =1) (3) ) (56) ((e0.t ) (3) (e0.t o =2) (3) ) (57) Eq: 86367: molarenthalpycomponentfeed.mosequ (usng Nota: 8639: not -. Descrpton: Molar enthalpy (feed, pure component), assumptons:. ncompressble flud deal gas behavor of correspondng vapor phase (pressure ndependency of h n V-phase) e0.h F,n =1 = e0.hn,o =1 + e0.a cp,=1 (e0.t F e0.t o =1) + e0.b cp,=1 2 ((e0.t F ) (2) (e0.t o =1) (2) ) + e0.c cp,=1 3 + e0.d cp,=1 ((e0.t F ) () (e0.t o =1) () ) + e0.v L,n =1 (e0.pf e0.p F,LV =1 ) e0.hf,lv,n =1 ((e0.t F ) (3) (e0.t o =1) (3) ) (58) 12

13 e0.h F,n =2 = e0.hn,o =2 + e0.a cp,=2 (e0.t F e0.t o =2) + e0.b cp,=2 2 ((e0.t F ) (2) (e0.t o =2) (2) ) + e0.c cp,=2 3 + e0.d cp,=2 ((e0.t F ) () (e0.t o =2) () ) + e0.v L,n =2 (e0.pf e0.p F,LV =2 ) e0.hf,lv,n =2 ((e0.t F ) (3) (e0.t o =2) (3) ) (59) Eq: 86369: vaporpressureantonelqud.mosequ (usng Nota: 8639: not -. Descrpton: Vapor pressure by Antone equaton (lqud, pure component), source: Kolbe, Mehlng p.59. ( ) e0.p B,LV =1 = (10) (3) (10) e0.p B,LV =2 = (10) (3) (10) ( e0.b e0.a Antone,=1 Antone,=1 e0.c Antone,=1 +(e0.t ) e0.b e0.a Antone,=2 Antone,=2 e0.c Antone,=2 +(e0.t ) Eq: 86370: vaporpressureantonefeed.mosequ (usng Nota: 8639: not. Descrpton: Vapor pressure by Antone equaton (feed, pure component), source: Kolbe, Mehlng p.59. ( ) e0.p F,LV =1 = (10)(3) (10) e0.p F,LV =2 = (10)(3) (10) ( e0.b e0.a Antone,=1 Antone,=1 e0.c Antone,=1 +(e0.t F ) e0.b e0.a Antone,=2 Antone,=2 e0.c Antone,=2 +(e0.t F ) ) ) (60) (61) (62) (63) Eq: 86371: molardenstylqud.mosequ (usng Nota: 8639: not. Descrpton: Molar densty wth no excess volume (lqud). 1 = e0.ρ B,n (e0.x B =1 e0.v L,n =1 + e0.xb =2 e0.v L,n =2 ) (6) Eq: 86372: volumeflowlqud.mosequ (usng Nota: 8639: not. Descrpton: Volume flow (lqud). e0.f B,n = e0.f B,v e0.ρ B,n (65) Eq: 86373: actvtycoeffwlsonlqud.mosequ (usng Nota: 8639: not. Descrpton: Actvty coeffcent calculatons by Wlson s g E -model (lqud). e0.γ B 1 =1 = e0.x B =1 + e0.αb =1 (1 e0.xb =1 ( ) exp (1 e0.x B e0.α=1 B =1) ( e0.x B =1 + e0.αb =1 (1 e0.xb =1 ) (e0.α=1 B + e0.αb =2 ) ) e0.αb =1 ((e0.α=1 B + e0.αb =2 ) e0.αb =1 ) e0.xb =1 + (1 e0.xb =1 )) (66) 13

14 e0.γ B 1 =2 = e0.x B =2 + e0.αb =2 (1 e0.xb =2 ( ) exp (1 e0.x B e0.α=2 B =2) ( e0.x B =2 + e0.αb =2 (1 e0.xb =2 ) (e0.α=1 B + e0.αb =2 ) ) e0.αb =2 ((e0.α=1 B + e0.αb =2 ) e0.αb =2 ) e0.xb =2 + (1 e0.xb =2 )) (67) Eq: 8637: actvtycoeffwlsonparameterlqud.mosequ (usng Nota: 8639: not. Descrpton: Parameter calculaton (lqud) based on Wlson s g E -model (lqud), source: Gmehlng, Kolbe S.239. e0.α B =1 = (e0.vl,n =1 + e0.vl,n =2 ) ( ) e0.vl,n =1 e0.λ=1 e0.v L,n exp e0.t =1 e0.α B =2 = (e0.vl,n =1 + e0.vl,n =2 ) ( ) e0.vl,n =2 e0.λ=2 e0.v L,n exp e0.t =2 (68) (69) Eq: 86380: crosssectonarea.mosequ (usng Nota: 8639: not. Descrpton: Flash s cross secton area. e0.a = e0.π (e0.d) (2) (70) Eq: 86381: volumetank.mosequ (usng Nota: 8639: not. Descrpton: Flash s volume. e0.hu v = e0.a e0.l (71) Eq: 86382: ratodameterlengthtank.mosequ (usng Nota: 8639: not. Descrpton: dameter to heght rato of tank. e0.r D,L = e0.d e0.l (72) Varable Specs 86386: IV statflash.mosvar Desgn varables e0.a Antone,=1 = e0.a Antone,=2 = e0.a P P DS12,=1 = e0.a P P DS12,=2 = e0.a cp,=1 = e0.a cp,=2 = e0.b Antone,=1 = e0.b Antone,=2 =

15 e0.b P P DS12,=1 = e0.b P P DS12,=2 = e0.b cp,=1 = e0.b cp,=2 = e0.c Antone,=1 = e0.c Antone,=2 = e0.c P P DS12,=1 = e0.c P P DS12,=2 = e0.c cp,=1 = 8.386E 5 e0.c cp,=2 = 1.058E 5 e0.d P P DS12,=1 = e0.d P P DS12,=2 = e0.d cp,=1 = E 9 e0.d cp,=2 = 3.59E 9 e0.e P P DS12,=1 = e0.e P P DS12,=2 = e0.f F,n = 1.75 e0.l = 0.5 e0.l L = 0.25 e0.r = 8.31 e0.t = e0.t F = e0.t=1 c = e0.t=1 o = e0.t=2 c = 67.3 e0.t=2 o = e0. p = e0.λ =1 = e0.λ =2 = e0.π = e0.d = 0.16 e0.h B,E,n = 0.0 e0.h D,E,n = 0.0 e0.h E,F,n = 0.0 e0.h n,o =1 = e0.h n,o =2 = e0.p = e0.v E,L,n = 0.0 e0.v L,n =1 = 5.869E 5 e0.v L,n =2 = 1.807E 5 e0.x F =1 = 0.15 e0.z = 1.0 Iteraton varables e0.a =

16 e0.f B,n = 1.0 e0.f B,v = 1.0 e0.f D,n = 1.0 e0.hu L,n = 1.0 e0.hu L,v = 1.0 e0.hu V,n = 1.0 e0.hu V,v = 1.0 e0.hu v = 1.0 e0.hu=1 n = 1.0 e0.hu=2 n = 1.0 e0.q = 1.0 e0.u = 1.0 e0.α=1 B = 1.0 e0.α=2 B = 1.0 e0.γ=1 B = 1.0 e0.γ=2 B = 1.0 e0.ρ B,n = 1.0 e0.h B,n = 1.0 e0.h D,n = 1.0 e0.h F,n = 1.0 e0.h B,LV,n =1 = 1.0 e0.h B,n =1 = 1.0 e0.h D,n =1 = 1.0 e0.h F,LV,n =1 = 1.0 e0.h F,n =1 = 1.0 e0.h B,LV,n =2 = 1.0 e0.h B,n =2 = 1.0 e0.h D,n =2 = 1.0 e0.h F,LV,n =2 = 1.0 e0.h F,n =2 = 1.0 e0.p F = 1.0 e0.p B,LV =1 = 1.0 e0.p F,LV =1 = 1.0 e0.p B,LV =2 = 1.0 e0.p F,LV =2 = 1.0 e0.r D,L = 1.0 e0.x B =1 = 1.0 e0.x B =2 = 1.0 e0.x F =2 = 1.0 e0.y=1 D = 1.0 e0.y=2 D =

17 Results 86388: RE statflash CBzzMath Org.mosvar Desgn varables e0.a Antone,=1 = e0.a Antone,=2 = e0.a P P DS12,=1 = e0.a P P DS12,=2 = e0.a cp,=1 = e0.a cp,=2 = e0.b Antone,=1 = e0.b Antone,=2 = e0.b P P DS12,=1 = e0.b P P DS12,=2 = e0.b cp,=1 = e0.b cp,=2 = e0.c Antone,=1 = e0.c Antone,=2 = e0.c P P DS12,=1 = e0.c P P DS12,=2 = e0.c cp,=1 = 8.386E 5 e0.c cp,=2 = 1.058E 5 e0.d P P DS12,=1 = e0.d P P DS12,=2 = e0.d cp,=1 = E 9 e0.d cp,=2 = 3.59E 9 e0.e P P DS12,=1 = e0.e P P DS12,=2 = e0.f F,n = 1.75 e0.l = 0.5 e0.l L = 0.25 e0.r = 8.31 e0.t = e0.t F = e0.t=1 c = e0.t=1 o = e0.t=2 c = 67.3 e0.t=2 o = e0. p = e0.λ =1 = e0.λ =2 = e0.π = e0.d = 0.16 e0.h B,E,n = 0.0 e0.h D,E,n = 0.0 e0.h E,F,n =

18 e0.h n,o =1 = e0.h n,o =2 = e0.p = e0.v E,L,n = 0.0 e0.v L,n =1 = 5.869E 5 e0.v L,n =2 = 1.807E 5 e0.x F =1 = 0.15 e0.z = 1.0 Iteraton varables e0.a = e0.f B,n = e0.f B,v = E 5 e0.f D,n = e0.hu L,n = e0.hu L,v = e0.hu V,n = E e0.hu V,v = e0.hu v = e0.hu=1 n = e0.hu=2 n = e0.q = e0.u = e0.α=1 B = e0.α=2 B = e0.γ=1 B = e0.γ=2 B = e0.ρ B,n = e0.h B,n = e0.h D,n = e0.h F,n = e0.h B,LV,n =1 = e0.h B,n =1 = e0.h D,n =1 = e0.h F,LV,n =1 = e0.h F,n =1 = e0.h B,LV,n =2 = e0.h B,n =2 = e0.h D,n =2 = e0.h F,LV,n =2 = e0.h F,n =2 = e0.p F = e0.p B,LV =1 =

19 e0.p F,LV e0.p B,LV e0.p F,LV =1 = =2 = =2 = e0.r D,L = 0.32 e0.x B =1 = e0.x B =2 = e0.x F =2 = 0.85 e0.y=1 D = e0.y=2 D = Notaton 8639: not flash.mosnot Base lne symbols A Area [m]; Parameter B Parameter C Parameter D Parameter E Parameter F Flow[mol/s] HU Holdup I Integral (by tme) K Equlbrum constant [1]; Control constant L Heght [m] M Molar mass [g/mol] Q Heat flow [W] R Ideal gas constant [J/mol/K] T Temperature [K]; Tme constant [s] U Internal Energy [J] F Flow dfference L Level dfference [m] Q Heat flow dfference [J/s] T Temperature dfference [K] p Pressure dfference [Pa] α Wlson Parameter [1] γ Actvty coeffcent [1] λ Interacton parameter (Wlson Modell) π Natural constant π [1] ρ Densty [mol/m] d Dameter [m] h Specfc enthalpy[j/mol]; p Pressure [Pa] r rato [1] t Tme [s] 19

20 v Specfc molar volume [mol/m] x Mol fracton (lqud) [1] y Mol fracton (vapor) [1] z Compressblty factor [1] Superscrpts B D E F L LV V c max n o v Bottom Dstllate Excess Feed Lqud; Heght Lqud/Vapor Vapor Crtcal Maxmum Molar Reference volumetrc Subscrpts Antone P P DS12 cp Antone equaton PPDS12 equaton ( molar heat capacty Indces 1..N C Komponentenndex 20