[715] | 1 | #*------------------------------------------------------------------- |
---|
| 2 | * EMSO Model Library (EML) Copyright (C) 2004 - 2007 ALSOC. |
---|
| 3 | * |
---|
| 4 | * This LIBRARY is free software; you can distribute it and/or modify |
---|
| 5 | * it under the therms of the ALSOC FREE LICENSE as available at |
---|
| 6 | * http://www.enq.ufrgs.br/alsoc. |
---|
| 7 | * |
---|
| 8 | * EMSO Copyright (C) 2004 - 2007 ALSOC, original code |
---|
| 9 | * from http://www.rps.eng.br Copyright (C) 2002-2004. |
---|
| 10 | * All rights reserved. |
---|
| 11 | * |
---|
| 12 | * EMSO is distributed under the therms of the ALSOC LICENSE as |
---|
| 13 | * available at http://www.enq.ufrgs.br/alsoc. |
---|
| 14 | * |
---|
| 15 | *---------------------------------------------------------------------- |
---|
| 16 | * File containg models of columns: distillation, stripping, absorbers |
---|
| 17 | * rectifier, .... |
---|
| 18 | * |
---|
| 19 | * The default nomenclature is: |
---|
| 20 | * Type_Column_reboilertype_condensertyper |
---|
| 21 | * |
---|
| 22 | * where: |
---|
| 23 | * Type = refluxed or reboiled or section |
---|
| 24 | * Column = Stripping, Absorption, Rectifier, Distillation |
---|
| 25 | * Reboiler type (if exists) = kettle or thermosyphon |
---|
| 26 | * Condenser type (if exists) = with subccoling or without subcooling |
---|
| 27 | * |
---|
| 28 | *----------------------------------------------------------------------- |
---|
| 29 | * Author: Paula B. Staudt |
---|
| 30 | * $Id: column.mso 511 2008-05-12 17:25:33Z paula $ |
---|
| 31 | *---------------------------------------------------------------------*# |
---|
| 32 | |
---|
| 33 | using "tray"; |
---|
| 34 | using "reboiler"; |
---|
| 35 | using "condenser"; |
---|
| 36 | using "mixers_splitters/splitter"; |
---|
| 37 | using "tank"; |
---|
| 38 | using "pressure_changers/pump"; |
---|
| 39 | |
---|
[771] | 40 | Model ControlSection |
---|
| 41 | |
---|
| 42 | ATTRIBUTES |
---|
| 43 | Pallete = false; |
---|
| 44 | Brief = "Control variables for column section."; |
---|
| 45 | |
---|
| 46 | PARAMETERS |
---|
| 47 | outer NComp as Integer (Brief="Number of components"); |
---|
| 48 | |
---|
| 49 | Tindicator_TrayNumber as Integer (Brief="Temperature Indicator in the tray"); |
---|
| 50 | Pindicator_TrayNumber as Integer (Brief="Pressure Indicator in the tray"); |
---|
| 51 | |
---|
| 52 | end |
---|
| 53 | |
---|
[715] | 54 | Model InitializeSection |
---|
| 55 | |
---|
| 56 | ATTRIBUTES |
---|
| 57 | Pallete = false; |
---|
[719] | 58 | Brief = "Initial conditions for column section."; |
---|
[715] | 59 | |
---|
| 60 | PARAMETERS |
---|
| 61 | outer NComp as Integer (Brief="Number of components"); |
---|
| 62 | |
---|
| 63 | TopTemperature as temperature (Brief = "Tray Temperature at Column Top", Default = 300); |
---|
| 64 | BottomTemperature as temperature (Brief = "Tray Temperature at Column Bottom", Default = 330); |
---|
[721] | 65 | LevelFraction as fraction (Brief = "Tray Level Fraction", Default = 0.5); |
---|
[715] | 66 | |
---|
| 67 | TopComposition(NComp) as fraction (Brief = "Component Molar Fraction at Column Top", Default = 0.30); |
---|
| 68 | BottomComposition(NComp) as fraction (Brief = "Component Molar Fraction at Column Bottom", Default = 0.30); |
---|
| 69 | |
---|
| 70 | end |
---|
| 71 | |
---|
[721] | 72 | Model InitializeStage |
---|
| 73 | |
---|
| 74 | ATTRIBUTES |
---|
| 75 | Pallete = false; |
---|
| 76 | Brief = "Initial conditions for Packed column section."; |
---|
| 77 | |
---|
| 78 | PARAMETERS |
---|
| 79 | outer NComp as Integer (Brief="Number of components"); |
---|
| 80 | |
---|
| 81 | TopStageTemperature as temperature (Brief = "Tray Temperature at Column Top", Default = 300); |
---|
| 82 | BottomStageTemperature as temperature (Brief = "Tray Temperature at Column Bottom", Default = 300); |
---|
| 83 | |
---|
| 84 | TopStageComposition(NComp) as fraction (Brief = "Component Molar Fraction at Top"); |
---|
| 85 | BottomStageComposition(NComp) as fraction (Brief = "Component Molar Fraction at Bottom"); |
---|
| 86 | LiquidMolarHoldup as mol (Brief="Molar liquid holdup", Default=0.01); |
---|
| 87 | |
---|
| 88 | end |
---|
| 89 | |
---|
[715] | 90 | |
---|
[787] | 91 | Model SectionColumnBasic |
---|
[775] | 92 | |
---|
[715] | 93 | ATTRIBUTES |
---|
[819] | 94 | Pallete = false; |
---|
[787] | 95 | Brief = "Model of a basic column section."; |
---|
[819] | 96 | Info = |
---|
[787] | 97 | "Model of a basic column section containing a vetor of TRAYS numbered from the top-down."; |
---|
[715] | 98 | |
---|
[775] | 99 | PARAMETERS |
---|
[787] | 100 | outer PP as Plugin (Brief="External Physical Properties", Type="PP"); |
---|
| 101 | outer NComp as Integer (Brief="Number of components"); |
---|
[819] | 102 | |
---|
[787] | 103 | NumberOfFeeds as Integer (Brief="Number of Feed Trays",Default=3,Protected=true); |
---|
| 104 | FeedTrayLocation(NumberOfFeeds) as Integer (Brief="Feed tray Location",Default=2); |
---|
[775] | 105 | NumberOfTrays as Integer (Brief="Number of trays", Default=8); |
---|
| 106 | LiqSideTrayIndex(NumberOfTrays) as Integer (Brief="Liquid Side Tray Index", Default=0,Hidden=true); |
---|
| 107 | VapSideTrayIndex(NumberOfTrays) as Integer (Brief="Vapour Side Tray Index", Default=0,Hidden=true); |
---|
[787] | 108 | |
---|
[775] | 109 | LiquidSideStreamLocation as Integer (Brief="Liquid Side Stream Location", Default=2); |
---|
| 110 | VapourSideStreamLocation as Integer (Brief="Vapour Side Stream Location", Default=2); |
---|
[816] | 111 | Pi as constant (Brief="Pi Number",Default=3.14159265, Symbol = "\pi",Hidden=true); |
---|
[831] | 112 | Gconst as acceleration (Brief="Gravity Acceleration",Default=9.81,Hidden=true); |
---|
[775] | 113 | Mw(NComp) as molweight (Brief="Component Mol Weight",Hidden=true); |
---|
[831] | 114 | zero_flow as flow_mol (Brief = "Stream Flow closed",Default = 0, Hidden=true); |
---|
| 115 | low_flow as flow_mol (Brief = "Low stream Flow",Default = 1E-6, Hidden=true); |
---|
[819] | 116 | |
---|
[775] | 117 | VapourFlowModel as Switcher (Valid = ["Reepmeyer", "Feehery_Fv", "Roffel_Fv", "Klingberg", "Wang_Fv", "Elgue"], Default = "Reepmeyer"); |
---|
[831] | 118 | LiquidFlowModel as Switcher (Valid = ["Francis", "Wang_Fl", "Olsen", "Feehery_Fl", "Roffel_Fl"], Default = "Francis"); |
---|
[819] | 119 | |
---|
[816] | 120 | TrayDiameter as length (Brief="Tray Diameter",Default=1.600); |
---|
| 121 | TraySpacing as length (Brief="Tray Spacing",Default=0.600); |
---|
| 122 | Fraction_HoleArea as fraction (Brief="Fraction of the active area that is occupied by the holes with respect to the total tray area",Default=0.10); |
---|
| 123 | Fraction_DowncomerArea as fraction (Brief="Fraction of the downcomer area with respect to the total tray area",Default=0.20); |
---|
[831] | 124 | WeirLength as length (Brief="Weir length", Default = 1); |
---|
| 125 | WeirHeight as length (Brief="Weir height", Default= 0.05); |
---|
| 126 | TrayLiquidPasses as positive (Brief="Number of liquid passes in the tray", Lower = 1,Default=1); |
---|
| 127 | HeatSupply as heat_rate (Brief="Rate of heat supply",Default = 0); |
---|
[819] | 128 | AerationFraction as Real (Brief="Aeration fraction", Default = 1); |
---|
| 129 | DryPdropCoeff as Real (Brief="Dry pressure drop coefficient", Default= 0.60); |
---|
[775] | 130 | |
---|
[819] | 131 | PlateArea as area (Brief="Plate area = Atray - Adowncomer",Protected=true); |
---|
| 132 | TrayVolume as volume (Brief="Total Volume of the tray",Protected=true); |
---|
| 133 | HolesArea as area (Brief="Total holes area",Protected=true); |
---|
[787] | 134 | |
---|
[819] | 135 | FeeheryCoeff as Real (Brief="Feeherys correlation coefficient", Unit='1/m^4', Default=1,Hidden=true); |
---|
| 136 | ElgueCoeff as Real (Brief="Elgues correlation coefficient", Unit='kg/m/mol^2', Default=1,Hidden=true); |
---|
| 137 | OlsenCoeff as Real (Brief="Olsens correlation coefficient", Default=1,Hidden=true); |
---|
[816] | 138 | |
---|
[831] | 139 | VapourFlow as Switcher (Brief="Flag for Vapour Flow condition",Valid = ["on", "off"], Default = "off",Hidden=true); |
---|
| 140 | LiquidFlow as Switcher (Brief="Flag for Liquid Flow condition",Valid = ["on", "off"], Default = "off",Hidden=true); |
---|
[775] | 141 | |
---|
| 142 | SET |
---|
[787] | 143 | VapSideTrayIndex(VapourSideStreamLocation) =1; |
---|
| 144 | LiqSideTrayIndex(LiquidSideStreamLocation) =1; |
---|
[775] | 145 | Mw = PP.MolecularWeight(); |
---|
[831] | 146 | zero_flow = 0 * 'kmol/h'; |
---|
| 147 | low_flow = 1E-6 * 'kmol/h'; |
---|
[775] | 148 | |
---|
[818] | 149 | PlateArea = 0.25*Pi*(TrayDiameter^2)*(1-Fraction_DowncomerArea); |
---|
[816] | 150 | TrayVolume = 0.25*Pi*(TrayDiameter^2)*TraySpacing; |
---|
| 151 | HolesArea = 0.25*Pi*(TrayDiameter^2)*Fraction_HoleArea; |
---|
| 152 | |
---|
[715] | 153 | VARIABLES |
---|
[819] | 154 | |
---|
| 155 | INITIALIZATION as InitializeSection (Brief="Column Model Initialization"); |
---|
| 156 | CONTROL as ControlSection (Brief="Column Model Control"); |
---|
| 157 | |
---|
[787] | 158 | out TCI as control_signal (Brief="Temperature Indicator", Protected = true, PosX=1, PosY=0.73); |
---|
| 159 | out PCI as control_signal (Brief="Pressure Indicator", Protected = true, PosX=0, PosY=0.24); |
---|
[819] | 160 | |
---|
| 161 | in LiquidInlet as stream (Brief="Liquid Inlet in the section", PosX=0.70, PosY=0); |
---|
| 162 | out VapourOutlet as vapour_stream (Brief="Vapour Outlet in the section", PosX=0.30, PosY=0); |
---|
[806] | 163 | |
---|
[819] | 164 | out VapourDrawOff as vapour_stream (Brief="Vapour Draw Off Port in the section", PosX=1, PosY=0.35,Protected = true); |
---|
| 165 | out LiquidDrawOff as liquid_stream (Brief="Liquid Draw Off Port in the section", PosX=1, PosY=0.65,Protected = true); |
---|
[787] | 166 | |
---|
[819] | 167 | TRAYS(NumberOfTrays) as tray (Brief="Number of trays in the Column Section",Protected=true); |
---|
[806] | 168 | VapourDrawOffFlow as flow_mol (Brief="Vapour Draw Off Stream Molar Flow Rate"); |
---|
| 169 | LiquidDrawOffFlow as flow_mol (Brief="Vapour Draw Off Stream Molar Flow Rate"); |
---|
| 170 | MurphreeEff as Real (Brief="Murphree efficiency for All Trays",Lower=0.01,Upper=1); |
---|
[787] | 171 | |
---|
[806] | 172 | in VapourInlet as stream (Brief="Vapour Inlet in the section", PosX=0.30, PosY=1); |
---|
| 173 | out LiquidOutlet as liquid_stream (Brief="Liquid Outlet in the section", PosX=0.70, PosY=1); |
---|
[787] | 174 | |
---|
[819] | 175 | LiquidConnector as stream (Brief="Liquid connection in the middle TRAYS", PosX=0.75, PosY=1,Hidden=true); |
---|
| 176 | VapourConnector as stream (Brief="Vapour connection in the middle TRAYS", PosX=0.55, PosY=0,Hidden=true); |
---|
| 177 | |
---|
[715] | 178 | CONNECTIONS |
---|
[775] | 179 | TRAYS([2:NumberOfTrays]).OutletVapour to TRAYS([1:NumberOfTrays-1]).InletVapour; |
---|
[819] | 180 | TRAYS([1:NumberOfTrays-1]).OutletLiquid to TRAYS([2:NumberOfTrays]).InletLiquid; |
---|
[787] | 181 | |
---|
| 182 | LiquidConnector to TRAYS(1).InletLiquid; |
---|
| 183 | VapourConnector to TRAYS(NumberOfTrays).InletVapour; |
---|
| 184 | |
---|
[775] | 185 | INITIAL |
---|
[819] | 186 | |
---|
[775] | 187 | for i in 1:NumberOfTrays do |
---|
| 188 | |
---|
| 189 | "The initial temperature of the TRAYS" |
---|
| 190 | TRAYS(i).OutletLiquid.T = INITIALIZATION.TopTemperature+(INITIALIZATION.BottomTemperature-INITIALIZATION.TopTemperature)*((i-1)/(NumberOfTrays-1)); |
---|
| 191 | |
---|
| 192 | "The initial Level of the TRAYS" |
---|
[816] | 193 | TRAYS(i).Level = INITIALIZATION.LevelFraction*WeirHeight; |
---|
[775] | 194 | end |
---|
| 195 | |
---|
| 196 | for i in 1:NComp-1 do |
---|
| 197 | |
---|
| 198 | for j in 1:NumberOfTrays do |
---|
| 199 | |
---|
| 200 | "The initial composition of the TRAYS - Normalized" |
---|
| 201 | TRAYS(j).OutletLiquid.z(i) = INITIALIZATION.TopComposition(i)/sum(INITIALIZATION.TopComposition) +(INITIALIZATION.BottomComposition(i)/sum(INITIALIZATION.BottomComposition)-INITIALIZATION.TopComposition(i)/sum(INITIALIZATION.TopComposition) )*((j-1)/(NumberOfTrays-1)); |
---|
| 202 | end |
---|
| 203 | |
---|
| 204 | end |
---|
| 205 | |
---|
[715] | 206 | EQUATIONS |
---|
[775] | 207 | |
---|
| 208 | "Tray Temperature Indicator" |
---|
| 209 | #TCI*'K' = TRAYS(min([NumberOfTrays, CONTROL.Tindicator_TrayNumber])).OutletVapour.T; |
---|
| 210 | TCI*'K' = TRAYS(CONTROL.Tindicator_TrayNumber).OutletVapour.T; |
---|
| 211 | |
---|
| 212 | "Tray Pressure Indicator" |
---|
| 213 | #PCI*'atm' = TRAYS(min([NumberOfTrays, CONTROL.Pindicator_TrayNumber])).OutletVapour.P; |
---|
| 214 | PCI*'atm' = TRAYS(CONTROL.Pindicator_TrayNumber).OutletVapour.P; |
---|
| 215 | |
---|
| 216 | for i in [1:NumberOfTrays] do |
---|
| 217 | |
---|
| 218 | "Murphree Efficiency" |
---|
| 219 | TRAYS(i).OutletVapour.z = MurphreeEff * (TRAYS(i).yideal - TRAYS(i).InletVapour.z) + TRAYS(i).InletVapour.z; |
---|
| 220 | |
---|
| 221 | "Level of clear liquid over the weir" |
---|
[816] | 222 | TRAYS(i).Level = TRAYS(i).ML*TRAYS(i).vL/PlateArea; |
---|
[775] | 223 | |
---|
| 224 | "Geometry Constraint" |
---|
[816] | 225 | TrayVolume = TRAYS(i).ML* TRAYS(i).vL + TRAYS(i).MV*TRAYS(i).vV; |
---|
[775] | 226 | |
---|
| 227 | "Energy Holdup" |
---|
[816] | 228 | TRAYS(i).E = TRAYS(i).ML*TRAYS(i).OutletLiquid.h + TRAYS(i).MV*TRAYS(i).OutletVapour.h - TRAYS(i).OutletLiquid.P*TrayVolume; |
---|
[775] | 229 | |
---|
| 230 | "Energy Balance" |
---|
| 231 | diff(TRAYS(i).E) = ( TRAYS(i).Inlet.F*TRAYS(i).Inlet.h + TRAYS(i).InletLiquid.F*TRAYS(i).InletLiquid.h + TRAYS(i).InletVapour.F*TRAYS(i).InletVapour.h- TRAYS(i).OutletLiquid.F*TRAYS(i).OutletLiquid.h - TRAYS(i).OutletVapour.F*TRAYS(i).OutletVapour.h |
---|
[816] | 232 | -TRAYS(i).VapourSideStream.F*TRAYS(i).VapourSideStream.h - TRAYS(i).LiquidSideStream.F*TRAYS(i).LiquidSideStream.h + HeatSupply ); |
---|
[775] | 233 | |
---|
| 234 | switch LiquidFlow |
---|
| 235 | case "on": |
---|
| 236 | switch LiquidFlowModel |
---|
[831] | 237 | case "Francis": |
---|
[775] | 238 | "Francis Equation" |
---|
[819] | 239 | TRAYS(i).OutletLiquid.F*TRAYS(i).vL = 1.84*'1/s'*WeirLength*((TRAYS(i).Level-(AerationFraction*WeirHeight))/(AerationFraction))^2; |
---|
[775] | 240 | |
---|
| 241 | case "Wang_Fl": |
---|
[819] | 242 | TRAYS(i).OutletLiquid.F*TRAYS(i).vL = 1.84*'m^0.5/s'*WeirLength*((TRAYS(i).Level-(AerationFraction*WeirHeight))/(AerationFraction))^1.5; |
---|
[775] | 243 | |
---|
| 244 | case "Olsen": |
---|
[816] | 245 | TRAYS(i).OutletLiquid.F / 'mol/s'= WeirLength*TrayLiquidPasses*TRAYS(i).rhoL/sum(Mw*TRAYS(i).OutletVapour.z)/(0.665*OlsenCoeff)^1.5 * ((TRAYS(i).ML*sum(Mw*TRAYS(i).OutletLiquid.z)/TRAYS(i).rhoL/PlateArea)-WeirHeight)^1.5 * 'm^0.5/mol'; |
---|
[775] | 246 | |
---|
| 247 | case "Feehery_Fl": |
---|
[816] | 248 | TRAYS(i).OutletLiquid.F = WeirLength*TRAYS(i).rhoL/sum(Mw*TRAYS(i).OutletLiquid.z) * ((TRAYS(i).Level-WeirHeight)/750/'mm')^1.5 * 'm^2/s'; |
---|
[775] | 249 | |
---|
| 250 | case "Roffel_Fl": |
---|
[831] | 251 | TRAYS(i).OutletLiquid.F = 2/3*TRAYS(i).rhoL/sum(Mw*TRAYS(i).OutletLiquid.z)*WeirLength*(TRAYS(i).ML*sum(Mw*TRAYS(i).OutletLiquid.z)/(PlateArea*1.3)/TRAYS(i).rhoL)^1.5*sqrt(2*Gconst/ |
---|
[816] | 252 | (2*(1 - 0.3593/'Pa^0.0888545'*abs(TRAYS(i).OutletVapour.F*sum(Mw*TRAYS(i).OutletVapour.z)/(PlateArea*1.3)/sqrt(TRAYS(i).rhoV))^0.177709)-1)); #/'(kg/m)^0.0888545/s^0.177709'; |
---|
[775] | 253 | end |
---|
[819] | 254 | when TRAYS(i).Level < (AerationFraction *WeirHeight) switchto "off"; |
---|
[775] | 255 | |
---|
| 256 | case "off": |
---|
[831] | 257 | |
---|
[775] | 258 | "Low level" |
---|
[831] | 259 | TRAYS(i).OutletLiquid.F = zero_flow; |
---|
| 260 | |
---|
[836] | 261 | when TRAYS(i).Level > (AerationFraction * WeirHeight) switchto "on"; |
---|
[831] | 262 | |
---|
[775] | 263 | end |
---|
| 264 | |
---|
| 265 | switch VapourFlow |
---|
| 266 | case "on": |
---|
| 267 | switch VapourFlowModel |
---|
| 268 | case "Reepmeyer": |
---|
[819] | 269 | TRAYS(i).InletVapour.F*TRAYS(i).vV = sqrt((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)/(TRAYS(i).rhoV*DryPdropCoeff))*HolesArea; |
---|
[775] | 270 | |
---|
| 271 | case "Feehery_Fv": |
---|
[831] | 272 | TRAYS(i).InletVapour.F = TRAYS(i).rhoV/PlateArea/FeeheryCoeff/sum(Mw*TRAYS(i).OutletVapour.z) * sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-(TRAYS(i).rhoV*Gconst*TRAYS(i).ML*TRAYS(i).vL/PlateArea))/TRAYS(i).rhoV); |
---|
[775] | 273 | |
---|
| 274 | case "Roffel_Fv": |
---|
[831] | 275 | TRAYS(i).InletVapour.F^1.08 * 0.0013 * 'kg/m/mol^1.08/s^0.92*1e5' = (TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)*1e5 - (AerationFraction*sum(TRAYS(i).M*Mw)/(PlateArea*1.3)*Gconst*1e5) * (TRAYS(i).rhoV*HolesArea/sum(Mw*TRAYS(i).OutletVapour.z))^1.08 * 'm^1.08/mol^1.08'; |
---|
[775] | 276 | |
---|
| 277 | case "Klingberg": |
---|
[831] | 278 | TRAYS(i).InletVapour.F * TRAYS(i).vV = PlateArea * sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-TRAYS(i).rhoL*Gconst*TRAYS(i).Level)/TRAYS(i).rhoV); |
---|
[775] | 279 | |
---|
| 280 | case "Wang_Fv": |
---|
[831] | 281 | TRAYS(i).InletVapour.F * TRAYS(i).vV = PlateArea * sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-TRAYS(i).rhoL*Gconst*TRAYS(i).Level)/TRAYS(i).rhoV*DryPdropCoeff); |
---|
[775] | 282 | |
---|
| 283 | case "Elgue": |
---|
[816] | 284 | TRAYS(i).InletVapour.F = sqrt((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)/ElgueCoeff); |
---|
[775] | 285 | end |
---|
[831] | 286 | when TRAYS(i).InletVapour.F < low_flow switchto "off"; |
---|
[775] | 287 | |
---|
| 288 | case "off": |
---|
[831] | 289 | TRAYS(i).InletVapour.F = zero_flow; |
---|
| 290 | |
---|
[836] | 291 | when TRAYS(i).InletVapour.P > TRAYS(i).OutletVapour.P switchto "on"; |
---|
[831] | 292 | |
---|
[775] | 293 | end |
---|
| 294 | |
---|
| 295 | end |
---|
| 296 | |
---|
[721] | 297 | "Liquid Inlet Flow" |
---|
| 298 | LiquidConnector.F = LiquidInlet.F; |
---|
[715] | 299 | |
---|
[721] | 300 | "Liquid Inlet Temperature" |
---|
[715] | 301 | LiquidConnector.T = LiquidInlet.T; |
---|
[721] | 302 | |
---|
| 303 | "Liquid Inlet Pressure" |
---|
[715] | 304 | LiquidConnector.P = LiquidInlet.P; |
---|
[721] | 305 | |
---|
| 306 | "Liquid Inlet Composition" |
---|
[715] | 307 | LiquidConnector.z = LiquidInlet.z; |
---|
[721] | 308 | |
---|
| 309 | "Liquid Inlet Vapour Fraction" |
---|
[715] | 310 | LiquidConnector.v = LiquidInlet.v; |
---|
[721] | 311 | |
---|
| 312 | "Liquid Inlet Enthalpy" |
---|
[715] | 313 | LiquidConnector.h = LiquidInlet.h; |
---|
[819] | 314 | |
---|
[721] | 315 | "Vapour Inlet Flow" |
---|
| 316 | VapourConnector.F = VapourInlet.F; |
---|
| 317 | |
---|
| 318 | "Vapour Inlet Temperature" |
---|
[715] | 319 | VapourConnector.T = VapourInlet.T; |
---|
[721] | 320 | |
---|
| 321 | "Vapour Inlet Pressure" |
---|
[715] | 322 | VapourConnector.P = VapourInlet.P; |
---|
[721] | 323 | |
---|
| 324 | "Vapour Inlet Composition" |
---|
[715] | 325 | VapourConnector.z = VapourInlet.z; |
---|
[721] | 326 | |
---|
| 327 | "Vapour Inlet Vapour Fraction" |
---|
[715] | 328 | VapourConnector.v = VapourInlet.v; |
---|
[721] | 329 | |
---|
| 330 | "Vapour Inlet Enthalpy" |
---|
[715] | 331 | VapourConnector.h = VapourInlet.h; |
---|
[819] | 332 | |
---|
[721] | 333 | "Liquid Outlet Flow" |
---|
[735] | 334 | LiquidOutlet.F = TRAYS(NumberOfTrays).OutletLiquid.F; |
---|
[721] | 335 | |
---|
| 336 | "Liquid Outlet Temperature" |
---|
[735] | 337 | LiquidOutlet.T = TRAYS(NumberOfTrays).OutletLiquid.T; |
---|
[721] | 338 | |
---|
| 339 | "Liquid Outlet Pressure" |
---|
[735] | 340 | LiquidOutlet.P = TRAYS(NumberOfTrays).OutletLiquid.P; |
---|
[721] | 341 | |
---|
| 342 | "Liquid Outlet Composition" |
---|
[735] | 343 | LiquidOutlet.z = TRAYS(NumberOfTrays).OutletLiquid.z; |
---|
[819] | 344 | |
---|
[721] | 345 | "Vapour Outlet Flow" |
---|
[735] | 346 | VapourOutlet.F = TRAYS(1).OutletVapour.F; |
---|
[721] | 347 | |
---|
| 348 | "Vapour Outlet Temperature" |
---|
[735] | 349 | VapourOutlet.T = TRAYS(1).OutletVapour.T; |
---|
[721] | 350 | |
---|
| 351 | "Vapour Outlet Pressure" |
---|
[735] | 352 | VapourOutlet.P = TRAYS(1).OutletVapour.P; |
---|
[721] | 353 | |
---|
| 354 | "Vapour Outlet Composition" |
---|
[735] | 355 | VapourOutlet.z = TRAYS(1).OutletVapour.z; |
---|
[819] | 356 | |
---|
[721] | 357 | "Vapour Draw Off Stream - Flow" |
---|
| 358 | VapourDrawOff.F*VapSideTrayIndex= TRAYS.VapourSideStream.F; |
---|
[715] | 359 | |
---|
[721] | 360 | "Vapour Draw Off Stream - Temperature" |
---|
[715] | 361 | VapourDrawOff.T = TRAYS(VapourSideStreamLocation).VapourSideStream.T; |
---|
[721] | 362 | |
---|
| 363 | "Vapour Draw Off Stream - Pressure" |
---|
[715] | 364 | VapourDrawOff.P = TRAYS(VapourSideStreamLocation).VapourSideStream.P; |
---|
[721] | 365 | |
---|
| 366 | "Vapour Draw Off Stream - Composition" |
---|
[715] | 367 | VapourDrawOff.z = TRAYS(VapourSideStreamLocation).VapourSideStream.z; |
---|
| 368 | |
---|
[721] | 369 | "Vapour Draw Off Stream" |
---|
| 370 | VapourDrawOffFlow = VapourDrawOff.F; |
---|
[819] | 371 | |
---|
[721] | 372 | "Liquid Draw Off Stream - Flow" |
---|
[715] | 373 | LiquidDrawOff.F*LiqSideTrayIndex= TRAYS.LiquidSideStream.F; |
---|
[721] | 374 | |
---|
| 375 | "Liquid Draw Off Stream - Temperature" |
---|
[715] | 376 | LiquidDrawOff.T = TRAYS(LiquidSideStreamLocation).LiquidSideStream.T; |
---|
[721] | 377 | |
---|
| 378 | "Liquid Draw Off Stream - Pressure" |
---|
[715] | 379 | LiquidDrawOff.P = TRAYS(LiquidSideStreamLocation).LiquidSideStream.P; |
---|
[721] | 380 | |
---|
| 381 | "Liquid Draw Off Stream - Composition" |
---|
[715] | 382 | LiquidDrawOff.z = TRAYS(LiquidSideStreamLocation).LiquidSideStream.z; |
---|
| 383 | |
---|
[721] | 384 | "Liquid Draw Off Stream" |
---|
[715] | 385 | LiquidDrawOffFlow = LiquidDrawOff.F; |
---|
| 386 | |
---|
| 387 | end |
---|
[719] | 388 | |
---|
[787] | 389 | Model Section_Column as SectionColumnBasic |
---|
[719] | 390 | |
---|
| 391 | ATTRIBUTES |
---|
[715] | 392 | Pallete = true; |
---|
[787] | 393 | Icon = "icon/SectionColumn"; |
---|
| 394 | Brief = "Model of a column section."; |
---|
| 395 | Info = |
---|
| 396 | "== Model of a column section containing == |
---|
| 397 | * NumberOfTrays TRAYS. |
---|
| 398 | * One Feed Inlet. |
---|
| 399 | |
---|
| 400 | == Specify == |
---|
| 401 | * the feed stream of each tray (Inlet); |
---|
| 402 | * the Murphree eficiency for each tray Emv; |
---|
| 403 | * the InletLiquid stream of the top tray; |
---|
| 404 | * the InletVapour stream of the bottom tray. |
---|
| 405 | |
---|
| 406 | == Initial Conditions == |
---|
| 407 | * the TRAYS temperature (OutletLiquid.T); |
---|
| 408 | * the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F); |
---|
| 409 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray. |
---|
| 410 | "; |
---|
| 411 | |
---|
| 412 | SET |
---|
| 413 | NumberOfFeeds = 1; |
---|
| 414 | |
---|
| 415 | VARIABLES |
---|
| 416 | |
---|
| 417 | in FeedTray as stream (Brief="Feed stream", PosX=0, PosY=0.55); |
---|
| 418 | |
---|
| 419 | EQUATIONS |
---|
| 420 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 421 | # Equating Feed Tray Variables to Trays Variables |
---|
| 422 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
[791] | 423 | for i in 1:NumberOfTrays do |
---|
[787] | 424 | |
---|
[791] | 425 | if i equal FeedTrayLocation(1) then |
---|
[787] | 426 | |
---|
[791] | 427 | "FeedTrayTop Inlet Flow" |
---|
| 428 | FeedTray.F= TRAYS(i).Inlet.F; |
---|
[787] | 429 | |
---|
[791] | 430 | "FeedTrayTop Inlet Temperature" |
---|
| 431 | FeedTray.T = TRAYS(i).Inlet.T; |
---|
[787] | 432 | |
---|
[791] | 433 | "FeedTrayTop Inlet Pressure" |
---|
| 434 | FeedTray.P = TRAYS(i).Inlet.P; |
---|
[787] | 435 | |
---|
[791] | 436 | "FeedTrayTop Inlet Composition" |
---|
| 437 | FeedTray.z = TRAYS(i).Inlet.z; |
---|
[787] | 438 | |
---|
[791] | 439 | "FeedTrayTop Inlet Vapour Fraction" |
---|
| 440 | FeedTray.v = TRAYS(i).Inlet.v; |
---|
| 441 | |
---|
| 442 | "FeedTrayTop Inlet Enthalpy" |
---|
| 443 | FeedTray.h = TRAYS(i).Inlet.h; |
---|
| 444 | |
---|
| 445 | else |
---|
| 446 | |
---|
[832] | 447 | "Inlet Tray - Flow Sealed" |
---|
| 448 | 0*'mol/h'= TRAYS(i).Inlet.F; |
---|
| 449 | |
---|
[791] | 450 | "Inlet Tray - Temperature" |
---|
[832] | 451 | 300*'K' = TRAYS(i).Inlet.T; |
---|
| 452 | |
---|
[791] | 453 | "Inlet Tray - Pressure" |
---|
[832] | 454 | 1*'atm' = TRAYS(i).Inlet.P; |
---|
| 455 | |
---|
[791] | 456 | "Inlet Tray - Composition" |
---|
[832] | 457 | 0.1 = TRAYS(i).Inlet.z; |
---|
| 458 | |
---|
[791] | 459 | "Inlet Tray - Vapour Fraction" |
---|
[832] | 460 | 0 = TRAYS(i).Inlet.v; |
---|
| 461 | |
---|
[791] | 462 | "Inlet Tray - Enthalpy" |
---|
[832] | 463 | 0*'J/mol' = TRAYS(i).Inlet.h; |
---|
[791] | 464 | |
---|
[787] | 465 | end |
---|
| 466 | |
---|
[791] | 467 | end |
---|
| 468 | |
---|
| 469 | end |
---|
| 470 | |
---|
[787] | 471 | Model Section_Column2 as SectionColumnBasic |
---|
| 472 | |
---|
| 473 | ATTRIBUTES |
---|
| 474 | Pallete = true; |
---|
| 475 | Icon = "icon/SectionColumn"; |
---|
| 476 | Brief = "Model of a column section."; |
---|
| 477 | Info = |
---|
| 478 | "== Model of a column section containing == |
---|
| 479 | * NumberOfTrays TRAYS. |
---|
| 480 | * Two Feed Inlets. |
---|
| 481 | |
---|
| 482 | == Specify == |
---|
| 483 | * the feed stream of each tray (Inlet); |
---|
| 484 | * the Murphree eficiency for each tray Emv; |
---|
| 485 | * the InletLiquid stream of the top tray; |
---|
| 486 | * the InletVapour stream of the bottom tray. |
---|
| 487 | |
---|
| 488 | == Initial Conditions == |
---|
| 489 | * the TRAYS temperature (OutletLiquid.T); |
---|
| 490 | * the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F); |
---|
| 491 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray. |
---|
| 492 | "; |
---|
| 493 | |
---|
| 494 | SET |
---|
| 495 | NumberOfFeeds = 2; |
---|
| 496 | |
---|
| 497 | VARIABLES |
---|
| 498 | |
---|
| 499 | in FeedTrayOne as stream (Brief="Feed stream", PosX=0, PosY=0.45); |
---|
| 500 | in FeedTrayTwo as stream (Brief="Feed stream", PosX=0, PosY=0.65); |
---|
| 501 | |
---|
| 502 | EQUATIONS |
---|
| 503 | |
---|
| 504 | for i in 1:NumberOfTrays do |
---|
| 505 | |
---|
| 506 | if i equal FeedTrayLocation(1) then |
---|
| 507 | |
---|
| 508 | "FeedTrayTop Inlet Flow" |
---|
| 509 | FeedTrayOne.F= TRAYS(i).Inlet.F; |
---|
| 510 | |
---|
| 511 | "FeedTrayTop Inlet Temperature" |
---|
| 512 | FeedTrayOne.T = TRAYS(i).Inlet.T; |
---|
| 513 | |
---|
| 514 | "FeedTrayTop Inlet Pressure" |
---|
| 515 | FeedTrayOne.P = TRAYS(i).Inlet.P; |
---|
| 516 | |
---|
| 517 | "FeedTrayTop Inlet Composition" |
---|
| 518 | FeedTrayOne.z = TRAYS(i).Inlet.z; |
---|
| 519 | |
---|
| 520 | "FeedTrayTop Inlet Vapour Fraction" |
---|
| 521 | FeedTrayOne.v = TRAYS(i).Inlet.v; |
---|
| 522 | |
---|
| 523 | "FeedTrayTop Inlet Enthalpy" |
---|
| 524 | FeedTrayOne.h = TRAYS(i).Inlet.h; |
---|
| 525 | |
---|
| 526 | else if i equal FeedTrayLocation(2) then |
---|
| 527 | |
---|
| 528 | "FeedTrayBottom Inlet Flow" |
---|
| 529 | FeedTrayTwo.F= TRAYS(i).Inlet.F; |
---|
| 530 | |
---|
| 531 | "FeedTrayBottom Inlet Temperature" |
---|
| 532 | FeedTrayTwo.T = TRAYS(i).Inlet.T; |
---|
| 533 | |
---|
| 534 | "FeedTrayBottom Inlet Pressure" |
---|
| 535 | FeedTrayTwo.P = TRAYS(i).Inlet.P; |
---|
| 536 | |
---|
| 537 | "FeedTrayBottom Inlet Composition" |
---|
| 538 | FeedTrayTwo.z = TRAYS(i).Inlet.z; |
---|
| 539 | |
---|
| 540 | "FeedTrayBottom Inlet Vapour Fraction" |
---|
| 541 | FeedTrayTwo.v = TRAYS(i).Inlet.v; |
---|
| 542 | |
---|
| 543 | "FeedTrayBottom Inlet Enthalpy" |
---|
| 544 | FeedTrayTwo.h = TRAYS(i).Inlet.h; |
---|
| 545 | |
---|
| 546 | else |
---|
| 547 | |
---|
[832] | 548 | "Inlet Tray - Flow Sealed" |
---|
| 549 | 0*'mol/h'= TRAYS(i).Inlet.F; |
---|
| 550 | |
---|
[787] | 551 | "Inlet Tray - Temperature" |
---|
[832] | 552 | 300*'K' = TRAYS(i).Inlet.T; |
---|
| 553 | |
---|
[787] | 554 | "Inlet Tray - Pressure" |
---|
[832] | 555 | 1*'atm' = TRAYS(i).Inlet.P; |
---|
| 556 | |
---|
[787] | 557 | "Inlet Tray - Composition" |
---|
[832] | 558 | 0.1 = TRAYS(i).Inlet.z; |
---|
| 559 | |
---|
[787] | 560 | "Inlet Tray - Vapour Fraction" |
---|
[832] | 561 | 0 = TRAYS(i).Inlet.v; |
---|
| 562 | |
---|
[787] | 563 | "Inlet Tray - Enthalpy" |
---|
[832] | 564 | 0*'J/mol' = TRAYS(i).Inlet.h; |
---|
[787] | 565 | |
---|
| 566 | end |
---|
| 567 | |
---|
| 568 | end |
---|
| 569 | |
---|
| 570 | end |
---|
| 571 | |
---|
| 572 | end |
---|
| 573 | |
---|
| 574 | Model Section_Column3 as SectionColumnBasic |
---|
| 575 | |
---|
| 576 | ATTRIBUTES |
---|
| 577 | Pallete = true; |
---|
| 578 | Icon = "icon/SectionColumn"; |
---|
| 579 | Brief = "Model of a column section."; |
---|
| 580 | Info = |
---|
| 581 | "== Model of a column section containing == |
---|
| 582 | * NumberOfTrays TRAYS. |
---|
| 583 | * Three Feed Inlets. |
---|
| 584 | |
---|
| 585 | == Specify == |
---|
| 586 | * the feed stream of each tray (Inlet); |
---|
| 587 | * the Murphree eficiency for each tray Emv; |
---|
| 588 | * the InletLiquid stream of the top tray; |
---|
| 589 | * the InletVapour stream of the bottom tray. |
---|
| 590 | |
---|
| 591 | == Initial Conditions == |
---|
| 592 | * the TRAYS temperature (OutletLiquid.T); |
---|
| 593 | * the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F); |
---|
| 594 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray. |
---|
| 595 | "; |
---|
| 596 | |
---|
| 597 | SET |
---|
| 598 | NumberOfFeeds = 3; |
---|
| 599 | |
---|
| 600 | VARIABLES |
---|
| 601 | |
---|
| 602 | in FeedTrayOne as stream (Brief="Feed stream", PosX=0, PosY=0.45); |
---|
| 603 | in FeedTrayTwo as stream (Brief="Feed stream", PosX=0, PosY=0.65); |
---|
| 604 | in FeedTrayThree as stream (Brief="Feed stream", PosX=0, PosY=0.85); |
---|
| 605 | |
---|
| 606 | EQUATIONS |
---|
| 607 | |
---|
| 608 | for i in 1:NumberOfTrays do |
---|
| 609 | |
---|
| 610 | if i equal FeedTrayLocation(1) then |
---|
| 611 | |
---|
| 612 | "FeedTrayTop Inlet Flow" |
---|
| 613 | FeedTrayOne.F= TRAYS(i).Inlet.F; |
---|
| 614 | |
---|
| 615 | "FeedTrayTop Inlet Temperature" |
---|
| 616 | FeedTrayOne.T = TRAYS(i).Inlet.T; |
---|
| 617 | |
---|
| 618 | "FeedTrayTop Inlet Pressure" |
---|
| 619 | FeedTrayOne.P = TRAYS(i).Inlet.P; |
---|
| 620 | |
---|
| 621 | "FeedTrayTop Inlet Composition" |
---|
| 622 | FeedTrayOne.z = TRAYS(i).Inlet.z; |
---|
| 623 | |
---|
| 624 | "FeedTrayTop Inlet Vapour Fraction" |
---|
| 625 | FeedTrayOne.v = TRAYS(i).Inlet.v; |
---|
| 626 | |
---|
| 627 | "FeedTrayTop Inlet Enthalpy" |
---|
| 628 | FeedTrayOne.h = TRAYS(i).Inlet.h; |
---|
| 629 | |
---|
| 630 | else if i equal FeedTrayLocation(2) then |
---|
| 631 | |
---|
| 632 | "FeedTrayBottom Inlet Flow" |
---|
| 633 | FeedTrayTwo.F= TRAYS(i).Inlet.F; |
---|
| 634 | |
---|
| 635 | "FeedTrayBottom Inlet Temperature" |
---|
| 636 | FeedTrayTwo.T = TRAYS(i).Inlet.T; |
---|
| 637 | |
---|
| 638 | "FeedTrayBottom Inlet Pressure" |
---|
| 639 | FeedTrayTwo.P = TRAYS(i).Inlet.P; |
---|
| 640 | |
---|
| 641 | "FeedTrayBottom Inlet Composition" |
---|
| 642 | FeedTrayTwo.z = TRAYS(i).Inlet.z; |
---|
| 643 | |
---|
| 644 | "FeedTrayBottom Inlet Vapour Fraction" |
---|
| 645 | FeedTrayTwo.v = TRAYS(i).Inlet.v; |
---|
| 646 | |
---|
| 647 | "FeedTrayBottom Inlet Enthalpy" |
---|
| 648 | FeedTrayTwo.h = TRAYS(i).Inlet.h; |
---|
| 649 | |
---|
| 650 | else if i equal FeedTrayLocation(3) then |
---|
| 651 | |
---|
| 652 | "FeedTrayBottom Inlet Flow" |
---|
| 653 | FeedTrayThree.F= TRAYS(i).Inlet.F; |
---|
| 654 | |
---|
| 655 | "FeedTrayBottom Inlet Temperature" |
---|
| 656 | FeedTrayThree.T = TRAYS(i).Inlet.T; |
---|
| 657 | |
---|
| 658 | "FeedTrayBottom Inlet Pressure" |
---|
| 659 | FeedTrayThree.P = TRAYS(i).Inlet.P; |
---|
| 660 | |
---|
| 661 | "FeedTrayBottom Inlet Composition" |
---|
| 662 | FeedTrayThree.z = TRAYS(i).Inlet.z; |
---|
| 663 | |
---|
| 664 | "FeedTrayBottom Inlet Vapour Fraction" |
---|
| 665 | FeedTrayThree.v = TRAYS(i).Inlet.v; |
---|
| 666 | |
---|
| 667 | "FeedTrayBottom Inlet Enthalpy" |
---|
| 668 | FeedTrayThree.h = TRAYS(i).Inlet.h; |
---|
| 669 | |
---|
| 670 | else |
---|
| 671 | |
---|
[832] | 672 | "Inlet Tray - Flow Sealed" |
---|
| 673 | 0*'mol/h'= TRAYS(i).Inlet.F; |
---|
| 674 | |
---|
[787] | 675 | "Inlet Tray - Temperature" |
---|
[832] | 676 | 300*'K' = TRAYS(i).Inlet.T; |
---|
| 677 | |
---|
[787] | 678 | "Inlet Tray - Pressure" |
---|
[832] | 679 | 1*'atm' = TRAYS(i).Inlet.P; |
---|
| 680 | |
---|
[787] | 681 | "Inlet Tray - Composition" |
---|
[832] | 682 | 0.1 = TRAYS(i).Inlet.z; |
---|
| 683 | |
---|
[787] | 684 | "Inlet Tray - Vapour Fraction" |
---|
[832] | 685 | 0 = TRAYS(i).Inlet.v; |
---|
| 686 | |
---|
[787] | 687 | "Inlet Tray - Enthalpy" |
---|
[832] | 688 | 0*'J/mol' = TRAYS(i).Inlet.h; |
---|
[787] | 689 | |
---|
| 690 | end |
---|
| 691 | |
---|
| 692 | end |
---|
| 693 | |
---|
| 694 | end |
---|
| 695 | |
---|
| 696 | end |
---|
| 697 | |
---|
| 698 | end |
---|
| 699 | |
---|
[793] | 700 | |
---|
[806] | 701 | Model ColumnBasic |
---|
[787] | 702 | |
---|
| 703 | ATTRIBUTES |
---|
[806] | 704 | Pallete = false; |
---|
| 705 | Brief = "Model of a basic column."; |
---|
| 706 | Info = |
---|
| 707 | "Model of a basic column containing a vetor of TRAYS numbered from the top-down."; |
---|
| 708 | |
---|
| 709 | PARAMETERS |
---|
| 710 | outer PP as Plugin (Brief="External Physical Properties", Type="PP"); |
---|
| 711 | outer NComp as Integer (Brief="Number of components"); |
---|
| 712 | |
---|
| 713 | NumberOfFeeds as Integer (Brief="Number of Feed Trays",Default=3,Protected=true); |
---|
| 714 | FeedTrayLocation(NumberOfFeeds) as Integer (Brief="Feed tray Location",Default=2); |
---|
| 715 | NumberOfTrays as Integer (Brief="Number of trays", Default=8); |
---|
| 716 | LiqSideTrayIndex(NumberOfTrays) as Integer (Brief="Liquid Side Tray Index", Default=0,Hidden=true); |
---|
| 717 | VapSideTrayIndex(NumberOfTrays) as Integer (Brief="Vapour Side Tray Index", Default=0,Hidden=true); |
---|
| 718 | LiquidSideStreamLocation as Integer (Brief="Liquid Side Stream Location", Default=2); |
---|
| 719 | VapourSideStreamLocation as Integer (Brief="Vapour Side Stream Location", Default=2); |
---|
| 720 | g as acceleration (Brief="Gravity Acceleration",Default=9.81,Hidden=true); |
---|
| 721 | Mw(NComp) as molweight (Brief="Component Mol Weight",Hidden=true); |
---|
| 722 | |
---|
| 723 | VapourFlowModel as Switcher (Valid = ["Reepmeyer", "Feehery_Fv", "Roffel_Fv", "Klingberg", "Wang_Fv", "Elgue"], Default = "Reepmeyer"); |
---|
| 724 | LiquidFlowModel as Switcher (Valid = ["default", "Wang_Fl", "Olsen", "Feehery_Fl", "Roffel_Fl"], Default = "default"); |
---|
| 725 | |
---|
[838] | 726 | TrayDiameter as length (Brief="Tray Diameter",Default=1.600); |
---|
| 727 | TraySpacing as length (Brief="Tray Spacing",Default=0.600); |
---|
| 728 | Fraction_HoleArea as fraction (Brief="Fraction of the active area that is occupied by the holes with respect to the total tray area",Default=0.10); |
---|
| 729 | Fraction_DowncomerArea as fraction (Brief="Fraction of the downcomer area with respect to the total tray area",Default=0.20); |
---|
| 730 | WeirLength as length (Brief="Weir length", Default = 1); |
---|
| 731 | WeirHeight as length (Brief="Weir height", Default= 0.05); |
---|
| 732 | TrayLiquidPasses as positive (Brief="Number of liquid passes in the tray", Lower = 1,Default=1); |
---|
| 733 | HeatSupply as heat_rate (Brief="Rate of heat supply",Default = 0); |
---|
| 734 | AerationFraction as Real (Brief="Aeration fraction", Default = 1); |
---|
| 735 | DryPdropCoeff as Real (Brief="Dry pressure drop coefficient", Default= 0.60); |
---|
[806] | 736 | |
---|
[838] | 737 | PlateArea as area (Brief="Plate area = Atray - Adowncomer",Protected=true); |
---|
| 738 | TrayVolume as volume (Brief="Total Volume of the tray",Protected=true); |
---|
| 739 | HolesArea as area (Brief="Total holes area",Protected=true); |
---|
| 740 | |
---|
[806] | 741 | FeeheryCoeff as Real (Brief="Feeherys correlation coefficient", Unit='1/m^4', Default=1,Hidden=true); |
---|
| 742 | ElgueCoeff as Real (Brief="Elgues correlation coefficient", Unit='kg/m/mol^2', Default=1,Hidden=true); |
---|
| 743 | OlsenCoeff as Real (Brief="Olsens correlation coefficient", Default=1,Hidden=true); |
---|
| 744 | |
---|
[838] | 745 | Pi as constant (Brief="Pi Number",Default=3.14159265, Symbol = "\pi",Hidden=true); |
---|
| 746 | zero_flow as flow_mol (Brief = "Stream Flow closed",Default = 0, Hidden=true); |
---|
| 747 | low_flow as flow_mol (Brief = "Low stream Flow",Default = 1E-6, Hidden=true); |
---|
[806] | 748 | VapourFlow as Switcher (Valid = ["on", "off"], Default = "on",Hidden=true); |
---|
| 749 | LiquidFlow as Switcher (Valid = ["on", "off"], Default = "on",Hidden=true); |
---|
| 750 | |
---|
| 751 | SET |
---|
| 752 | VapSideTrayIndex(VapourSideStreamLocation) =1; |
---|
| 753 | LiqSideTrayIndex(LiquidSideStreamLocation) =1; |
---|
| 754 | Mw = PP.MolecularWeight(); |
---|
[838] | 755 | zero_flow = 0 * 'kmol/h'; |
---|
| 756 | low_flow = 1E-6 * 'kmol/h'; |
---|
[806] | 757 | |
---|
[838] | 758 | PlateArea = 0.25*Pi*(TrayDiameter^2)*(1-Fraction_DowncomerArea); |
---|
| 759 | TrayVolume = 0.25*Pi*(TrayDiameter^2)*TraySpacing; |
---|
| 760 | HolesArea = 0.25*Pi*(TrayDiameter^2)*Fraction_HoleArea; |
---|
[806] | 761 | |
---|
| 762 | VARIABLES |
---|
| 763 | |
---|
| 764 | INITIALIZATION as InitializeSection (Brief = "Column Model Initialization"); |
---|
| 765 | CONTROL as ControlSection (Brief = "Control"); |
---|
| 766 | |
---|
[809] | 767 | out TCI as control_signal (Brief="Temperature Indicator", Protected = true, PosX=1, PosY=0.55); |
---|
[806] | 768 | out PCI as control_signal (Brief="Pressure Indicator", Protected = true, PosX=0, PosY=0.23); |
---|
| 769 | |
---|
| 770 | TRAYS(NumberOfTrays) as tray (Brief="Number of trays in the Column Section"); |
---|
| 771 | VapourDrawOffFlow as flow_mol (Brief="Vapour Draw Off Stream Molar Flow Rate"); |
---|
| 772 | LiquidDrawOffFlow as flow_mol (Brief="Vapour Draw Off Stream Molar Flow Rate"); |
---|
| 773 | MurphreeEff as Real (Brief="Murphree efficiency for All Trays",Lower=0.01,Upper=1); |
---|
| 774 | |
---|
| 775 | out VapourDrawOff as vapour_stream (Brief="Vapour Outlet in the section", PosX=1, PosY=0.388,Protected = true); |
---|
| 776 | out LiquidDrawOff as liquid_stream (Brief="Liquid Outlet in the section", PosX=1, PosY=0.47,Protected = true); |
---|
| 777 | |
---|
| 778 | CONNECTIONS |
---|
| 779 | |
---|
| 780 | #Connecting Intermediate Trays |
---|
| 781 | TRAYS([2:NumberOfTrays]).OutletVapour to TRAYS([1:NumberOfTrays-1]).InletVapour; |
---|
| 782 | TRAYS([1:NumberOfTrays-1]).OutletLiquid to TRAYS([2:NumberOfTrays]).InletLiquid; |
---|
| 783 | |
---|
| 784 | INITIAL |
---|
| 785 | |
---|
| 786 | for i in 1:NumberOfTrays do |
---|
| 787 | |
---|
| 788 | "The initial temperature of the TRAYS" |
---|
| 789 | TRAYS(i).OutletLiquid.T = INITIALIZATION.TopTemperature+(INITIALIZATION.BottomTemperature-INITIALIZATION.TopTemperature)*((i-1)/(NumberOfTrays-1)); |
---|
| 790 | |
---|
| 791 | "The initial Level of the TRAYS" |
---|
[838] | 792 | TRAYS(i).Level = INITIALIZATION.LevelFraction*WeirHeight; |
---|
[806] | 793 | end |
---|
| 794 | |
---|
| 795 | for i in 1:NComp-1 do |
---|
| 796 | |
---|
| 797 | for j in 1:NumberOfTrays do |
---|
| 798 | |
---|
| 799 | "The initial composition of the TRAYS - Normalized" |
---|
| 800 | TRAYS(j).OutletLiquid.z(i) = INITIALIZATION.TopComposition(i)/sum(INITIALIZATION.TopComposition) +(INITIALIZATION.BottomComposition(i)/sum(INITIALIZATION.BottomComposition)-INITIALIZATION.TopComposition(i)/sum(INITIALIZATION.TopComposition) )*((j-1)/(NumberOfTrays-1)); |
---|
| 801 | end |
---|
| 802 | |
---|
| 803 | end |
---|
| 804 | |
---|
| 805 | EQUATIONS |
---|
| 806 | |
---|
| 807 | VapourDrawOff.F*VapSideTrayIndex= TRAYS.VapourSideStream.F; |
---|
| 808 | VapourDrawOff.T = TRAYS(VapourSideStreamLocation).VapourSideStream.T; |
---|
| 809 | VapourDrawOff.P = TRAYS(VapourSideStreamLocation).VapourSideStream.P; |
---|
| 810 | VapourDrawOff.z = TRAYS(VapourSideStreamLocation).VapourSideStream.z; |
---|
| 811 | |
---|
| 812 | LiquidDrawOff.F*LiqSideTrayIndex= TRAYS.LiquidSideStream.F; |
---|
| 813 | LiquidDrawOff.T = TRAYS(LiquidSideStreamLocation).LiquidSideStream.T; |
---|
| 814 | LiquidDrawOff.P = TRAYS(LiquidSideStreamLocation).LiquidSideStream.P; |
---|
| 815 | LiquidDrawOff.z = TRAYS(LiquidSideStreamLocation).LiquidSideStream.z; |
---|
| 816 | |
---|
| 817 | VapourDrawOffFlow = VapourDrawOff.F; |
---|
| 818 | LiquidDrawOffFlow = LiquidDrawOff.F; |
---|
| 819 | |
---|
| 820 | "Tray Temperature Indicator" |
---|
| 821 | #TCI*'K' = TRAYS(min([NumberOfTrays, CONTROL.Tindicator_TrayNumber])).OutletVapour.T; |
---|
| 822 | TCI*'K' = TRAYS(CONTROL.Tindicator_TrayNumber).OutletVapour.T; |
---|
| 823 | |
---|
| 824 | "Tray Pressure Indicator" |
---|
| 825 | #PCI*'atm' = TRAYS(min([NumberOfTrays, CONTROL.Pindicator_TrayNumber])).OutletVapour.P; |
---|
| 826 | PCI*'atm' = TRAYS(CONTROL.Pindicator_TrayNumber).OutletVapour.P; |
---|
| 827 | |
---|
| 828 | for i in [1:NumberOfTrays] do |
---|
| 829 | |
---|
| 830 | "Murphree Efficiency" |
---|
| 831 | TRAYS(i).OutletVapour.z = MurphreeEff * (TRAYS(i).yideal - TRAYS(i).InletVapour.z) + TRAYS(i).InletVapour.z; |
---|
| 832 | |
---|
| 833 | "Level of clear liquid over the weir" |
---|
[838] | 834 | TRAYS(i).Level = TRAYS(i).ML*TRAYS(i).vL/PlateArea; |
---|
[806] | 835 | |
---|
| 836 | "Geometry Constraint" |
---|
[838] | 837 | TrayVolume = TRAYS(i).ML* TRAYS(i).vL + TRAYS(i).MV*TRAYS(i).vV; |
---|
[806] | 838 | |
---|
| 839 | "Energy Holdup" |
---|
[838] | 840 | TRAYS(i).E = TRAYS(i).ML*TRAYS(i).OutletLiquid.h + TRAYS(i).MV*TRAYS(i).OutletVapour.h - TRAYS(i).OutletLiquid.P*TrayVolume; |
---|
[806] | 841 | |
---|
| 842 | "Energy Balance" |
---|
| 843 | diff(TRAYS(i).E) = ( TRAYS(i).Inlet.F*TRAYS(i).Inlet.h + TRAYS(i).InletLiquid.F*TRAYS(i).InletLiquid.h + TRAYS(i).InletVapour.F*TRAYS(i).InletVapour.h- TRAYS(i).OutletLiquid.F*TRAYS(i).OutletLiquid.h - TRAYS(i).OutletVapour.F*TRAYS(i).OutletVapour.h |
---|
[838] | 844 | -TRAYS(i).VapourSideStream.F*TRAYS(i).VapourSideStream.h - TRAYS(i).LiquidSideStream.F*TRAYS(i).LiquidSideStream.h + HeatSupply ); |
---|
[806] | 845 | |
---|
| 846 | switch LiquidFlow |
---|
| 847 | case "on": |
---|
| 848 | switch LiquidFlowModel |
---|
| 849 | case "default": |
---|
| 850 | "Francis Equation" |
---|
[838] | 851 | TRAYS(i).OutletLiquid.F*TRAYS(i).vL = 1.84*'1/s'*WeirLength*((TRAYS(i).Level-(AerationFraction*WeirHeight))/(AerationFraction))^2; |
---|
[806] | 852 | |
---|
| 853 | case "Wang_Fl": |
---|
[838] | 854 | TRAYS(i).OutletLiquid.F*TRAYS(i).vL = 1.84*'m^0.5/s'*WeirLength*((TRAYS(i).Level-(AerationFraction*WeirHeight))/(AerationFraction))^1.5; |
---|
[806] | 855 | |
---|
| 856 | case "Olsen": |
---|
[838] | 857 | TRAYS(i).OutletLiquid.F / 'mol/s'= WeirLength*TrayLiquidPasses*TRAYS(i).rhoL/sum(Mw*TRAYS(i).OutletVapour.z)/(0.665*OlsenCoeff)^1.5 * ((TRAYS(i).ML*sum(Mw*TRAYS(i).OutletLiquid.z)/TRAYS(i).rhoL/PlateArea)-WeirHeight)^1.5 * 'm^0.5/mol'; |
---|
[806] | 858 | |
---|
| 859 | case "Feehery_Fl": |
---|
[838] | 860 | TRAYS(i).OutletLiquid.F = WeirLength*TRAYS(i).rhoL/sum(Mw*TRAYS(i).OutletLiquid.z) * ((TRAYS(i).Level-WeirHeight)/750/'mm')^1.5 * 'm^2/s'; |
---|
[806] | 861 | |
---|
| 862 | case "Roffel_Fl": |
---|
[838] | 863 | TRAYS(i).OutletLiquid.F = 2/3*TRAYS(i).rhoL/sum(Mw*TRAYS(i).OutletLiquid.z)*WeirLength*(TRAYS(i).ML*sum(Mw*TRAYS(i).OutletLiquid.z)/(PlateArea*1.3)/TRAYS(i).rhoL)^1.5*sqrt(2*g/ |
---|
| 864 | (2*(1 - 0.3593/'Pa^0.0888545'*abs(TRAYS(i).OutletVapour.F*sum(Mw*TRAYS(i).OutletVapour.z)/(PlateArea*1.3)/sqrt(TRAYS(i).rhoV))^0.177709)-1)); #/'(kg/m)^0.0888545/s^0.177709'; |
---|
[806] | 865 | end |
---|
[838] | 866 | when TRAYS(i).Level < (AerationFraction *WeirHeight) switchto "off"; |
---|
[806] | 867 | |
---|
| 868 | case "off": |
---|
| 869 | "Low level" |
---|
| 870 | TRAYS(i).OutletLiquid.F = 0 * 'mol/h'; |
---|
[838] | 871 | when TRAYS(i).Level > (AerationFraction * WeirHeight) switchto "on"; |
---|
[806] | 872 | end |
---|
| 873 | |
---|
| 874 | switch VapourFlow |
---|
| 875 | case "on": |
---|
| 876 | switch VapourFlowModel |
---|
| 877 | case "Reepmeyer": |
---|
[838] | 878 | TRAYS(i).InletVapour.F*TRAYS(i).vV = sqrt((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)/(TRAYS(i).rhoV*DryPdropCoeff))*HolesArea; |
---|
[806] | 879 | |
---|
| 880 | case "Feehery_Fv": |
---|
[838] | 881 | TRAYS(i).InletVapour.F = TRAYS(i).rhoV/PlateArea/FeeheryCoeff/sum(Mw*TRAYS(i).OutletVapour.z) * sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-(TRAYS(i).rhoV*g*TRAYS(i).ML*TRAYS(i).vL/PlateArea))/TRAYS(i).rhoV); |
---|
[806] | 882 | |
---|
| 883 | case "Roffel_Fv": |
---|
[838] | 884 | TRAYS(i).InletVapour.F^1.08 * 0.0013 * 'kg/m/mol^1.08/s^0.92*1e5' = (TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)*1e5 - (AerationFraction*sum(TRAYS(i).M*Mw)/(PlateArea*1.3)*g*1e5) * (TRAYS(i).rhoV*HolesArea/sum(Mw*TRAYS(i).OutletVapour.z))^1.08 * 'm^1.08/mol^1.08'; |
---|
[806] | 885 | |
---|
| 886 | case "Klingberg": |
---|
[838] | 887 | TRAYS(i).InletVapour.F * TRAYS(i).vV = PlateArea * sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-TRAYS(i).rhoL*g*TRAYS(i).Level)/TRAYS(i).rhoV); |
---|
[806] | 888 | |
---|
| 889 | case "Wang_Fv": |
---|
[838] | 890 | TRAYS(i).InletVapour.F * TRAYS(i).vV = PlateArea * sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-TRAYS(i).rhoL*g*TRAYS(i).Level)/TRAYS(i).rhoV*DryPdropCoeff); |
---|
[806] | 891 | |
---|
| 892 | case "Elgue": |
---|
[838] | 893 | TRAYS(i).InletVapour.F = sqrt((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)/ElgueCoeff); |
---|
[806] | 894 | end |
---|
| 895 | when TRAYS(i).InletVapour.F < 1e-6 * 'kmol/h' switchto "off"; |
---|
| 896 | |
---|
| 897 | case "off": |
---|
| 898 | TRAYS(i).InletVapour.F = 0 * 'mol/s'; |
---|
[836] | 899 | when TRAYS(i).InletVapour.P > TRAYS(i).OutletVapour.P switchto "on"; |
---|
[806] | 900 | end |
---|
| 901 | |
---|
| 902 | end |
---|
| 903 | |
---|
| 904 | end |
---|
| 905 | |
---|
[810] | 906 | Model Distillation_kettle_cond as ColumnBasic |
---|
[806] | 907 | |
---|
| 908 | ATTRIBUTES |
---|
[787] | 909 | Pallete = true; |
---|
[715] | 910 | Icon = "icon/DistillationKettleCond"; |
---|
| 911 | Brief = "Model of a distillation column with dynamic condenser and dynamic reboiler."; |
---|
| 912 | Info = |
---|
| 913 | "== Specify == |
---|
| 914 | * the feed stream of each tray (Inlet); |
---|
| 915 | * the Murphree eficiency for each tray Emv; |
---|
| 916 | * the pump pressure difference; |
---|
| 917 | * the heat supllied in reboiler and condenser; |
---|
[735] | 918 | * the condenser vapor outlet flow (OutletVapour.F); |
---|
| 919 | * the reboiler liquid outlet flow (OutletLiquid.F); |
---|
[715] | 920 | * both splitter outlet flows OR one of the splitter outlet flows and the splitter frac. |
---|
| 921 | |
---|
| 922 | == Initial Conditions == |
---|
[735] | 923 | * the TRAYS temperature (OutletLiquid.T); |
---|
| 924 | * the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F); |
---|
| 925 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray; |
---|
[715] | 926 | |
---|
[735] | 927 | * the condenser temperature (OutletLiquid.T); |
---|
[715] | 928 | * the condenser liquid level (Level); |
---|
[735] | 929 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions; |
---|
[715] | 930 | |
---|
[735] | 931 | * the reboiler temperature (OutletLiquid.T); |
---|
[715] | 932 | * the reboiler liquid level (Level); |
---|
[735] | 933 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions. |
---|
[715] | 934 | "; |
---|
[806] | 935 | |
---|
| 936 | SET |
---|
| 937 | NumberOfFeeds = 1; |
---|
[715] | 938 | |
---|
[806] | 939 | VARIABLES |
---|
[715] | 940 | |
---|
[838] | 941 | CONDENSER as condenser; |
---|
| 942 | REBOILER as reboiler; |
---|
| 943 | SPLITTER as splitter_column (Brief="splitter which separate reflux and distillate",Protected=true); |
---|
| 944 | PUMP as pump; |
---|
| 945 | RefluxRatio as positive (Brief="Reflux Ratio",Default=10, Lower = 0.5); |
---|
| 946 | |
---|
[806] | 947 | out TI_reb as control_signal (Brief="Temperature Indicator of Reboiler ",Protected = true, PosX=0.515, PosY=1); |
---|
| 948 | out LI_reb as control_signal (Brief="Level Indicator of Reboiler", Protected = true, PosX=1, PosY=0.935); |
---|
| 949 | out PI_reb as control_signal (Brief="Pressure Indicator of Reboiler", Protected = true, PosX=0.41, PosY=1); |
---|
| 950 | |
---|
| 951 | out TI_cond as control_signal (Brief="Temperature Indicator of Condenser ", Protected = true, PosX=0.52, PosY=0); |
---|
| 952 | out LI_cond as control_signal (Brief="Level Indicator of Condenser", Protected = true, PosX=1, PosY=0.084); |
---|
| 953 | out PI_cond as control_signal (Brief="Pressure Indicator of Condenser", Protected = true, PosX=0.42, PosY=0); |
---|
| 954 | |
---|
| 955 | in FeedTray as stream (Brief="Feed stream", PosX=0, PosY=0.55); |
---|
| 956 | |
---|
| 957 | |
---|
[838] | 958 | |
---|
[806] | 959 | in HeatToReboiler as power (Brief="Heat supplied to Reboiler",Protected = true, PosX=1, PosY=0.855); |
---|
| 960 | in HeatToCondenser as power (Brief="Heat supplied to Condenser", Protected = true, PosX=1, PosY=0.01); |
---|
| 961 | |
---|
| 962 | |
---|
| 963 | ConnectorHeatReboiler as power (Brief="Heat supplied to Reboiler", Hidden=true); |
---|
| 964 | ConnectorHeatCondenser as power (Brief="Heat supplied to Condenser", Hidden=true); |
---|
| 965 | |
---|
[721] | 966 | in ConnectorCondenserVout as stream (Brief="Connector for Vapour outlet stream From Condenser", Hidden=true); |
---|
| 967 | in ConnectorReboilerLout as stream (Brief="Connector for Liquid outlet stream From Reboiler", Hidden=true); |
---|
| 968 | in ConnectorSplitterOut as stream (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true); |
---|
[715] | 969 | |
---|
[806] | 970 | out VapourDistillate as vapour_stream (Brief="Vapour outlet stream From Condenser", PosX=0.615, PosY=0); |
---|
| 971 | out LiquidDistillate as liquid_stream (Brief="Liquid outlet stream From Top Splitter", PosX=1, PosY=0.235); |
---|
| 972 | out BottomProduct as liquid_stream (Brief="Liquid outlet stream From Reboiler", PosX=0.62, PosY=1); |
---|
| 973 | |
---|
[715] | 974 | EQUATIONS |
---|
| 975 | |
---|
[838] | 976 | SPLITTER.RefluxRatio = RefluxRatio; |
---|
| 977 | |
---|
[806] | 978 | TI_reb = REBOILER.TI ; |
---|
| 979 | LI_reb = REBOILER.LI; |
---|
| 980 | PI_reb = REBOILER.PI; |
---|
| 981 | |
---|
| 982 | TI_cond = CONDENSER.TI; |
---|
| 983 | LI_cond = CONDENSER.LI; |
---|
| 984 | PI_cond = CONDENSER.PI; |
---|
| 985 | |
---|
| 986 | for i in 1:NumberOfTrays do |
---|
[715] | 987 | |
---|
[806] | 988 | if i equal FeedTrayLocation(1) then |
---|
[715] | 989 | |
---|
[806] | 990 | "FeedTrayTop Inlet Flow" |
---|
| 991 | FeedTray.F= TRAYS(i).Inlet.F; |
---|
[715] | 992 | |
---|
[806] | 993 | "FeedTrayTop Inlet Temperature" |
---|
| 994 | FeedTray.T = TRAYS(i).Inlet.T; |
---|
[715] | 995 | |
---|
[806] | 996 | "FeedTrayTop Inlet Pressure" |
---|
| 997 | FeedTray.P = TRAYS(i).Inlet.P; |
---|
| 998 | |
---|
| 999 | "FeedTrayTop Inlet Composition" |
---|
| 1000 | FeedTray.z = TRAYS(i).Inlet.z; |
---|
| 1001 | |
---|
| 1002 | "FeedTrayTop Inlet Vapour Fraction" |
---|
| 1003 | FeedTray.v = TRAYS(i).Inlet.v; |
---|
| 1004 | |
---|
| 1005 | "FeedTrayTop Inlet Enthalpy" |
---|
| 1006 | FeedTray.h = TRAYS(i).Inlet.h; |
---|
| 1007 | |
---|
| 1008 | else |
---|
| 1009 | |
---|
[832] | 1010 | "Inlet Tray - Flow Sealed" |
---|
| 1011 | 0*'mol/h'= TRAYS(i).Inlet.F; |
---|
| 1012 | |
---|
[806] | 1013 | "Inlet Tray - Temperature" |
---|
[832] | 1014 | 300*'K' = TRAYS(i).Inlet.T; |
---|
| 1015 | |
---|
[806] | 1016 | "Inlet Tray - Pressure" |
---|
[832] | 1017 | 1*'atm' = TRAYS(i).Inlet.P; |
---|
| 1018 | |
---|
[806] | 1019 | "Inlet Tray - Composition" |
---|
[832] | 1020 | 0.1 = TRAYS(i).Inlet.z; |
---|
| 1021 | |
---|
[806] | 1022 | "Inlet Tray - Vapour Fraction" |
---|
[832] | 1023 | 0 = TRAYS(i).Inlet.v; |
---|
| 1024 | |
---|
[806] | 1025 | "Inlet Tray - Enthalpy" |
---|
[832] | 1026 | 0*'J/mol' = TRAYS(i).Inlet.h; |
---|
[806] | 1027 | |
---|
| 1028 | end |
---|
| 1029 | |
---|
| 1030 | end |
---|
| 1031 | |
---|
[715] | 1032 | # Condenser Connector Equations |
---|
| 1033 | ConnectorCondenserVout.T = VapourDistillate.T; |
---|
| 1034 | ConnectorCondenserVout.P = VapourDistillate.P; |
---|
| 1035 | ConnectorCondenserVout.F = VapourDistillate.F; |
---|
| 1036 | ConnectorCondenserVout.z = VapourDistillate.z; |
---|
| 1037 | |
---|
| 1038 | # Splitter Connector Equations |
---|
| 1039 | ConnectorSplitterOut.T = LiquidDistillate.T; |
---|
| 1040 | ConnectorSplitterOut.P = LiquidDistillate.P; |
---|
| 1041 | ConnectorSplitterOut.F = LiquidDistillate.F; |
---|
| 1042 | ConnectorSplitterOut.z = LiquidDistillate.z; |
---|
| 1043 | |
---|
| 1044 | # Reboiler Connector Equations |
---|
| 1045 | ConnectorReboilerLout.T = BottomProduct.T; |
---|
| 1046 | ConnectorReboilerLout.P = BottomProduct.P; |
---|
| 1047 | ConnectorReboilerLout.F = BottomProduct.F; |
---|
| 1048 | ConnectorReboilerLout.z = BottomProduct.z; |
---|
| 1049 | |
---|
[719] | 1050 | HeatToReboiler = ConnectorHeatReboiler; |
---|
| 1051 | HeatToCondenser = ConnectorHeatCondenser; |
---|
| 1052 | |
---|
[715] | 1053 | CONNECTIONS |
---|
| 1054 | #vapor |
---|
[735] | 1055 | REBOILER.OutletVapour to TRAYS(NumberOfTrays).InletVapour; |
---|
| 1056 | TRAYS(1).OutletVapour to CONDENSER.InletVapour; |
---|
[715] | 1057 | |
---|
| 1058 | #liquid |
---|
[735] | 1059 | CONDENSER.OutletLiquid to SPLITTER.Inlet; |
---|
[838] | 1060 | SPLITTER.Reflux to PUMP.Inlet; |
---|
[806] | 1061 | PUMP.Outlet to TRAYS(1).InletLiquid; |
---|
[735] | 1062 | TRAYS(NumberOfTrays).OutletLiquid to REBOILER.InletLiquid; |
---|
[715] | 1063 | |
---|
| 1064 | #Connectors |
---|
[735] | 1065 | CONDENSER.OutletVapour to ConnectorCondenserVout; |
---|
[838] | 1066 | SPLITTER.Distillate to ConnectorSplitterOut; |
---|
[735] | 1067 | REBOILER.OutletLiquid to ConnectorReboilerLout; |
---|
[715] | 1068 | |
---|
[719] | 1069 | ConnectorHeatReboiler to REBOILER.InletQ; |
---|
| 1070 | ConnectorHeatCondenser to CONDENSER.InletQ; |
---|
[715] | 1071 | |
---|
| 1072 | end |
---|
| 1073 | |
---|
[808] | 1074 | Model Distillation_thermosyphon_subcooling as ColumnBasic |
---|
[715] | 1075 | ATTRIBUTES |
---|
| 1076 | Pallete = true; |
---|
| 1077 | Icon = "icon/DistillationThermosyphonSubcooling"; |
---|
| 1078 | Brief = "Model of a distillation column with steady condenser and steady reboiler."; |
---|
| 1079 | Info = |
---|
| 1080 | "== Specify == |
---|
| 1081 | * the feed stream of each tray (Inlet); |
---|
| 1082 | * the Murphree eficiency for each tray Emv; |
---|
| 1083 | * the pump head; |
---|
| 1084 | * the condenser pressure drop; |
---|
| 1085 | * the heat supllied in top and bottom tanks; |
---|
| 1086 | * the heat supllied in condenser and reboiler; |
---|
| 1087 | * the Outlet1 flow in the bottom splitter (spbottom.Outlet1.F) that corresponds to the bottom product; |
---|
| 1088 | * both top splitter outlet flows OR one of the splitter outlet flows and the splitter frac. |
---|
| 1089 | |
---|
| 1090 | == Initial Conditions == |
---|
[735] | 1091 | * the TRAYS temperature (OutletLiquid.T); |
---|
| 1092 | * the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F); |
---|
| 1093 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray; |
---|
[715] | 1094 | |
---|
[735] | 1095 | * the top tank temperature (OutletLiquid.T); |
---|
[715] | 1096 | * the top tank liquid level (Level); |
---|
[735] | 1097 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions; |
---|
[715] | 1098 | |
---|
[735] | 1099 | * the bottom tank temperature (OutletLiquid.T); |
---|
[715] | 1100 | * the bottom tank liquid level (Level); |
---|
[735] | 1101 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions. |
---|
[715] | 1102 | "; |
---|
| 1103 | |
---|
[808] | 1104 | SET |
---|
[715] | 1105 | |
---|
[808] | 1106 | NumberOfFeeds = 1; |
---|
[715] | 1107 | |
---|
| 1108 | VARIABLES |
---|
[721] | 1109 | CONDENSER as condenserSteady (Brief="steady state condenser with subcooling"); |
---|
[808] | 1110 | REFLUX_DRUM as TankL (Brief="Reflux Drum"); |
---|
[838] | 1111 | SPLITTER_TOP as splitter_column (Brief="splitter which separate reflux and distillate",Protected=true); |
---|
[721] | 1112 | PUMP as pump (Brief="pump in reflux stream"); |
---|
| 1113 | REBOILER as reboilerSteady (Brief="steady state reboiler (thermosyphon)"); |
---|
[808] | 1114 | SUMP as SumpTank (Brief="vessel in the bottom of column"); |
---|
| 1115 | SPLITTER_BOTTOM as splitter2 (Brief="splitter to separate the bottom product and the stream to reboiler"); |
---|
| 1116 | |
---|
| 1117 | in FeedTray as stream (Brief="Feed stream", PosX=0, PosY=0.55); |
---|
| 1118 | |
---|
[838] | 1119 | RefluxRatio as positive (Brief="Reflux Ratio",Default=10, Lower = 0.5); |
---|
| 1120 | |
---|
[808] | 1121 | in Q_cond as power (Brief="Heat supplied to Condenser",PosX=1, PosY=0, Protected=true); |
---|
| 1122 | in Q_reb as power (Brief="Heat supplied to Reboiler", PosX=1, PosY=0.96,Protected=true); |
---|
| 1123 | in Q_drum as power (Brief="Heat supplied to Reflux Drum", PosX=1, PosY=0.30,Protected=true); |
---|
| 1124 | |
---|
| 1125 | out TI_reb as control_signal (Brief="Temperature Indicator of Reboiler ",Protected = true, PosX=1, PosY=0.92); |
---|
[809] | 1126 | out PI_reb as control_signal (Brief="Pressure Indicator of Reboiler", Protected = true, PosX=1, PosY=0.875); |
---|
[808] | 1127 | |
---|
| 1128 | out TI_cond as control_signal (Brief="Temperature Indicator of Condenser ", Protected = true, PosX=0.52, PosY=0); |
---|
| 1129 | out PI_cond as control_signal (Brief="Pressure Indicator of Condenser", Protected = true, PosX=0.42, PosY=0); |
---|
| 1130 | |
---|
[809] | 1131 | out TI_drum as control_signal (Brief="Temperature Indicator of Reflux Drum ", Protected = true, PosX=1, PosY=0.20); |
---|
| 1132 | out PI_drum as control_signal (Brief="Pressure Indicator of Reflux Drum", Protected = true, PosX=1, PosY=0.155); |
---|
[808] | 1133 | out LI_drum as control_signal (Brief="Level Indicator of Reflux Drum", Protected = true, PosX=1, PosY=0.24); |
---|
| 1134 | |
---|
[809] | 1135 | out LI_sump as control_signal (Brief="Level Indicator of Column Sump", Protected = true, PosX=1, PosY=0.765); |
---|
| 1136 | out TI_sump as control_signal (Brief="Temperature Indicator of Column Sump ",Protected = true, PosX=1, PosY=0.71); |
---|
| 1137 | |
---|
| 1138 | out LiquidDistillate as liquid_stream (Brief="Liquid outlet stream From Top Splitter", PosX=1, PosY=0.33); |
---|
| 1139 | out BottomProduct as liquid_stream (Brief="Liquid outlet stream From Bottom Splitter", PosX=0.06, PosY=1); |
---|
[808] | 1140 | |
---|
[809] | 1141 | in ConnectorSplitterBottom as stream (Brief="Connector for Liquid outlet stream From Reboiler", Hidden=true); |
---|
| 1142 | in ConnectorSplitterTop as stream (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true); |
---|
[808] | 1143 | |
---|
| 1144 | ConnectorHeatReboiler as power (Brief="Connector for Heat supplied to Reboiler", Hidden=true); |
---|
| 1145 | ConnectorHeatCondenser as power (Brief="Connector for Heat supplied to Condenser", Hidden=true); |
---|
| 1146 | ConnectorHeatRefluxDrum as power (Brief="Connector for Heat supplied to Reflux Drum", Hidden=true); |
---|
[809] | 1147 | |
---|
[715] | 1148 | EQUATIONS |
---|
[781] | 1149 | |
---|
[838] | 1150 | SPLITTER_TOP.RefluxRatio = RefluxRatio; |
---|
[809] | 1151 | TI_reb = REBOILER.TI ; |
---|
| 1152 | PI_reb = REBOILER.PI; |
---|
| 1153 | |
---|
| 1154 | TI_cond = CONDENSER.TI; |
---|
| 1155 | PI_cond = CONDENSER.PI; |
---|
| 1156 | |
---|
| 1157 | TI_drum = REFLUX_DRUM.TI; |
---|
| 1158 | PI_drum = REFLUX_DRUM.PI; |
---|
| 1159 | LI_drum = REFLUX_DRUM.LI; |
---|
| 1160 | |
---|
| 1161 | TI_sump = SUMP.TI; |
---|
| 1162 | LI_sump = SUMP.LI; |
---|
| 1163 | |
---|
| 1164 | for i in 1:NumberOfTrays do |
---|
| 1165 | |
---|
| 1166 | if i equal FeedTrayLocation(1) then |
---|
| 1167 | |
---|
| 1168 | "FeedTrayTop Inlet Flow" |
---|
| 1169 | FeedTray.F= TRAYS(i).Inlet.F; |
---|
| 1170 | |
---|
| 1171 | "FeedTrayTop Inlet Temperature" |
---|
| 1172 | FeedTray.T = TRAYS(i).Inlet.T; |
---|
| 1173 | |
---|
| 1174 | "FeedTrayTop Inlet Pressure" |
---|
| 1175 | FeedTray.P = TRAYS(i).Inlet.P; |
---|
| 1176 | |
---|
| 1177 | "FeedTrayTop Inlet Composition" |
---|
| 1178 | FeedTray.z = TRAYS(i).Inlet.z; |
---|
| 1179 | |
---|
| 1180 | "FeedTrayTop Inlet Vapour Fraction" |
---|
| 1181 | FeedTray.v = TRAYS(i).Inlet.v; |
---|
| 1182 | |
---|
| 1183 | "FeedTrayTop Inlet Enthalpy" |
---|
| 1184 | FeedTray.h = TRAYS(i).Inlet.h; |
---|
| 1185 | |
---|
| 1186 | else |
---|
| 1187 | |
---|
[832] | 1188 | "Inlet Tray - Flow Sealed" |
---|
| 1189 | 0*'mol/h'= TRAYS(i).Inlet.F; |
---|
| 1190 | |
---|
[809] | 1191 | "Inlet Tray - Temperature" |
---|
[832] | 1192 | 300*'K' = TRAYS(i).Inlet.T; |
---|
| 1193 | |
---|
[809] | 1194 | "Inlet Tray - Pressure" |
---|
[832] | 1195 | 1*'atm' = TRAYS(i).Inlet.P; |
---|
| 1196 | |
---|
[809] | 1197 | "Inlet Tray - Composition" |
---|
[832] | 1198 | 0.1 = TRAYS(i).Inlet.z; |
---|
| 1199 | |
---|
[809] | 1200 | "Inlet Tray - Vapour Fraction" |
---|
[832] | 1201 | 0 = TRAYS(i).Inlet.v; |
---|
| 1202 | |
---|
[809] | 1203 | "Inlet Tray - Enthalpy" |
---|
[832] | 1204 | 0*'J/mol' = TRAYS(i).Inlet.h; |
---|
[809] | 1205 | |
---|
| 1206 | end |
---|
| 1207 | |
---|
| 1208 | end |
---|
| 1209 | |
---|
[721] | 1210 | # Heat Connector Equations |
---|
[808] | 1211 | Q_cond = ConnectorHeatCondenser; |
---|
| 1212 | Q_reb = ConnectorHeatReboiler; |
---|
| 1213 | Q_drum = ConnectorHeatRefluxDrum; |
---|
[715] | 1214 | |
---|
| 1215 | # Top Splitter Connector Equations |
---|
| 1216 | ConnectorSplitterTop.T = LiquidDistillate.T; |
---|
| 1217 | ConnectorSplitterTop.P = LiquidDistillate.P; |
---|
| 1218 | ConnectorSplitterTop.F = LiquidDistillate.F; |
---|
| 1219 | ConnectorSplitterTop.z = LiquidDistillate.z; |
---|
| 1220 | |
---|
| 1221 | # Bottom Splitter Connector Equations |
---|
| 1222 | ConnectorSplitterBottom.T = BottomProduct.T; |
---|
| 1223 | ConnectorSplitterBottom.P = BottomProduct.P; |
---|
| 1224 | ConnectorSplitterBottom.F = BottomProduct.F; |
---|
| 1225 | ConnectorSplitterBottom.z = BottomProduct.z; |
---|
[809] | 1226 | |
---|
[715] | 1227 | CONNECTIONS |
---|
| 1228 | #vapor |
---|
[808] | 1229 | REBOILER.OutletVapour to SUMP.InletVapour; |
---|
| 1230 | SUMP.OutletVapour to TRAYS(NumberOfTrays).InletVapour; |
---|
[735] | 1231 | TRAYS(1).OutletVapour to CONDENSER.InletVapour; |
---|
[715] | 1232 | |
---|
| 1233 | #liquid |
---|
[808] | 1234 | CONDENSER.OutletLiquid to REFLUX_DRUM.Inlet; |
---|
| 1235 | REFLUX_DRUM.OutletLiquid to SPLITTER_TOP.Inlet; |
---|
[838] | 1236 | SPLITTER_TOP.Reflux to PUMP.Inlet; |
---|
[735] | 1237 | PUMP.Outlet to TRAYS(1).InletLiquid; |
---|
[808] | 1238 | TRAYS(NumberOfTrays).OutletLiquid to SUMP.InletLiquid; |
---|
| 1239 | SUMP.OutletLiquid to SPLITTER_BOTTOM.Inlet; |
---|
[735] | 1240 | SPLITTER_BOTTOM.Outlet2 to REBOILER.InletLiquid; |
---|
[715] | 1241 | |
---|
| 1242 | #Connectors |
---|
[721] | 1243 | ConnectorHeatCondenser to CONDENSER.InletQ; |
---|
| 1244 | ConnectorHeatReboiler to REBOILER.InletQ; |
---|
[808] | 1245 | ConnectorHeatRefluxDrum to REFLUX_DRUM.InletQ; |
---|
[838] | 1246 | SPLITTER_TOP.Distillate to ConnectorSplitterTop; |
---|
[721] | 1247 | SPLITTER_BOTTOM.Outlet1 to ConnectorSplitterBottom; |
---|
[715] | 1248 | |
---|
| 1249 | end |
---|
| 1250 | |
---|
[810] | 1251 | Model Distillation_thermosyphon_cond as ColumnBasic |
---|
[715] | 1252 | ATTRIBUTES |
---|
| 1253 | Pallete = true; |
---|
| 1254 | Icon = "icon/DistillationThermosyphonCond"; |
---|
| 1255 | Brief = "Model of a distillation column with dynamic condenser and steady reboiler."; |
---|
| 1256 | Info = |
---|
| 1257 | "== Specify == |
---|
| 1258 | * the feed stream of each tray (Inlet); |
---|
| 1259 | * the Murphree eficiency for each tray Emv; |
---|
| 1260 | * the pump head; |
---|
[735] | 1261 | * the condenser vapor outlet flow (OutletVapour.F); |
---|
[715] | 1262 | * the heat supllied in bottom tank; |
---|
| 1263 | * the heat supllied in condenser and reboiler; |
---|
| 1264 | * the Outlet1 flow in the bottom splitter (spbottom.Outlet1.F) that corresponds to the bottom product; |
---|
| 1265 | |
---|
| 1266 | == Initial Conditions == |
---|
[735] | 1267 | * the TRAYS temperature (OutletLiquid.T); |
---|
| 1268 | * the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F); |
---|
| 1269 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray; |
---|
[715] | 1270 | |
---|
[735] | 1271 | * the condenser temperature (OutletLiquid.T); |
---|
[715] | 1272 | * the condenser liquid level (Level); |
---|
[735] | 1273 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions; |
---|
[715] | 1274 | |
---|
[735] | 1275 | * the bottom tank temperature (OutletLiquid.T); |
---|
[715] | 1276 | * the bottom tank liquid level (Level); |
---|
[735] | 1277 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions. |
---|
[715] | 1278 | "; |
---|
| 1279 | |
---|
[810] | 1280 | SET |
---|
[715] | 1281 | |
---|
[810] | 1282 | NumberOfFeeds = 1; |
---|
| 1283 | |
---|
[715] | 1284 | VARIABLES |
---|
[721] | 1285 | CONDENSER as condenser (Brief="dynamic condenser without subcooling"); |
---|
[838] | 1286 | SPLITTER_TOP as splitter_column (Brief="splitter which separate reflux and distillate",Protected=true); |
---|
[721] | 1287 | PUMP as pump (Brief="pump in reflux stream"); |
---|
[810] | 1288 | SUMP as SumpTank (Brief="vessel in the bottom of column"); |
---|
| 1289 | SPLITTER_BOTTOM as splitter2 (Brief="splitter who separate the bottom product and the stream to reboiler"); |
---|
[721] | 1290 | REBOILER as reboilerSteady (Brief="steady state reboiler (thermosyphon)"); |
---|
[838] | 1291 | RefluxRatio as positive (Brief="Reflux Ratio",Default=10, Lower = 0.5); |
---|
| 1292 | |
---|
[810] | 1293 | in FeedTray as stream (Brief="Feed stream", PosX=0, PosY=0.55); |
---|
| 1294 | |
---|
| 1295 | in Q_cond as power (Brief="Heat supplied to Condenser",PosX=1, PosY=0, Protected=true); |
---|
| 1296 | in Q_reb as power (Brief="Heat supplied to Reboiler", PosX=1, PosY=0.96,Protected=true); |
---|
| 1297 | |
---|
| 1298 | out TI_cond as control_signal (Brief="Temperature Indicator of Condenser ", Protected = true, PosX=0.48, PosY=0); |
---|
| 1299 | out LI_cond as control_signal (Brief="Level Indicator of Condenser", Protected = true, PosX=1, PosY=0.084); |
---|
| 1300 | out PI_cond as control_signal (Brief="Pressure Indicator of Condenser", Protected = true, PosX=0.40, PosY=0); |
---|
| 1301 | |
---|
| 1302 | out TI_reb as control_signal (Brief="Temperature Indicator of Reboiler ",Protected = true, PosX=1, PosY=0.92); |
---|
| 1303 | out PI_reb as control_signal (Brief="Pressure Indicator of Reboiler", Protected = true, PosX=1, PosY=0.878); |
---|
| 1304 | |
---|
| 1305 | out LI_sump as control_signal (Brief="Level Indicator of Column Sump", Protected = true, PosX=1, PosY=0.76); |
---|
| 1306 | out TI_sump as control_signal (Brief="Temperature Indicator of Column Sump ",Protected = true, PosX=1, PosY=0.71); |
---|
| 1307 | |
---|
| 1308 | out VapourDistillate as vapour_stream (Brief="Vapour outlet stream From Top Condenser", PosX=0.62, PosY=0); |
---|
| 1309 | out LiquidDistillate as liquid_stream (Brief="Liquid outlet stream From Top Splitter", PosX=1, PosY=0.185); |
---|
| 1310 | out BottomProduct as liquid_stream (Brief="Liquid outlet stream From Bottom Splitter", PosX=0.06, PosY=1); |
---|
| 1311 | |
---|
[721] | 1312 | in ConnectorSplitterBottom as stream (Brief="Connector for Liquid outlet stream From Reboiler", Hidden=true); |
---|
| 1313 | in ConnectorSplitterTop as stream (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true); |
---|
| 1314 | in ConnectorCondenserVout as stream (Brief="Connector for Vapour outlet stream From Top Condenser", Hidden=true); |
---|
[715] | 1315 | |
---|
[810] | 1316 | ConnectorHeatReboiler as power (Brief="Connector for Heat supplied to Reboiler", Hidden=true); |
---|
| 1317 | ConnectorHeatCondenser as power (Brief="Connector for Heat supplied to Condenser", Hidden=true); |
---|
[715] | 1318 | |
---|
| 1319 | EQUATIONS |
---|
[810] | 1320 | |
---|
| 1321 | for i in 1:NumberOfTrays do |
---|
| 1322 | |
---|
| 1323 | if i equal FeedTrayLocation(1) then |
---|
| 1324 | |
---|
| 1325 | "FeedTrayTop Inlet Flow" |
---|
| 1326 | FeedTray.F= TRAYS(i).Inlet.F; |
---|
| 1327 | |
---|
| 1328 | "FeedTrayTop Inlet Temperature" |
---|
| 1329 | FeedTray.T = TRAYS(i).Inlet.T; |
---|
| 1330 | |
---|
| 1331 | "FeedTrayTop Inlet Pressure" |
---|
| 1332 | FeedTray.P = TRAYS(i).Inlet.P; |
---|
| 1333 | |
---|
| 1334 | "FeedTrayTop Inlet Composition" |
---|
| 1335 | FeedTray.z = TRAYS(i).Inlet.z; |
---|
| 1336 | |
---|
| 1337 | "FeedTrayTop Inlet Vapour Fraction" |
---|
| 1338 | FeedTray.v = TRAYS(i).Inlet.v; |
---|
| 1339 | |
---|
| 1340 | "FeedTrayTop Inlet Enthalpy" |
---|
| 1341 | FeedTray.h = TRAYS(i).Inlet.h; |
---|
| 1342 | |
---|
| 1343 | else |
---|
| 1344 | |
---|
[832] | 1345 | "Inlet Tray - Flow Sealed" |
---|
| 1346 | 0*'mol/h'= TRAYS(i).Inlet.F; |
---|
| 1347 | |
---|
[810] | 1348 | "Inlet Tray - Temperature" |
---|
[832] | 1349 | 300*'K' = TRAYS(i).Inlet.T; |
---|
| 1350 | |
---|
[810] | 1351 | "Inlet Tray - Pressure" |
---|
[832] | 1352 | 1*'atm' = TRAYS(i).Inlet.P; |
---|
| 1353 | |
---|
[810] | 1354 | "Inlet Tray - Composition" |
---|
[832] | 1355 | 0.1 = TRAYS(i).Inlet.z; |
---|
| 1356 | |
---|
[810] | 1357 | "Inlet Tray - Vapour Fraction" |
---|
[832] | 1358 | 0 = TRAYS(i).Inlet.v; |
---|
| 1359 | |
---|
[810] | 1360 | "Inlet Tray - Enthalpy" |
---|
[832] | 1361 | 0*'J/mol' = TRAYS(i).Inlet.h; |
---|
[810] | 1362 | |
---|
| 1363 | end |
---|
| 1364 | |
---|
| 1365 | end |
---|
| 1366 | |
---|
[838] | 1367 | SPLITTER_TOP.RefluxRatio = RefluxRatio; |
---|
| 1368 | |
---|
[810] | 1369 | TI_reb = REBOILER.TI ; |
---|
| 1370 | PI_reb = REBOILER.PI; |
---|
| 1371 | |
---|
| 1372 | TI_cond = CONDENSER.TI; |
---|
| 1373 | PI_cond = CONDENSER.PI; |
---|
| 1374 | LI_cond = CONDENSER.LI; |
---|
| 1375 | |
---|
| 1376 | TI_sump = SUMP.TI; |
---|
| 1377 | LI_sump = SUMP.LI; |
---|
| 1378 | |
---|
[721] | 1379 | # Heat Connector Equations |
---|
[810] | 1380 | Q_cond = ConnectorHeatCondenser; |
---|
| 1381 | Q_reb = ConnectorHeatReboiler; |
---|
[715] | 1382 | |
---|
| 1383 | # Condenser Connector Equations |
---|
| 1384 | ConnectorCondenserVout.T = VapourDistillate.T; |
---|
| 1385 | ConnectorCondenserVout.P = VapourDistillate.P; |
---|
| 1386 | ConnectorCondenserVout.F = VapourDistillate.F; |
---|
| 1387 | ConnectorCondenserVout.z = VapourDistillate.z; |
---|
| 1388 | |
---|
| 1389 | # Top Splitter Connector Equations |
---|
| 1390 | ConnectorSplitterTop.T = LiquidDistillate.T; |
---|
| 1391 | ConnectorSplitterTop.P = LiquidDistillate.P; |
---|
| 1392 | ConnectorSplitterTop.F = LiquidDistillate.F; |
---|
| 1393 | ConnectorSplitterTop.z = LiquidDistillate.z; |
---|
| 1394 | |
---|
| 1395 | # Bottom Splitter Connector Equations |
---|
| 1396 | ConnectorSplitterBottom.T = BottomProduct.T; |
---|
| 1397 | ConnectorSplitterBottom.P = BottomProduct.P; |
---|
| 1398 | ConnectorSplitterBottom.F = BottomProduct.F; |
---|
| 1399 | ConnectorSplitterBottom.z = BottomProduct.z; |
---|
| 1400 | |
---|
| 1401 | CONNECTIONS |
---|
| 1402 | #vapor |
---|
[810] | 1403 | SUMP.OutletVapour to TRAYS(NumberOfTrays).InletVapour; |
---|
| 1404 | REBOILER.OutletVapour to SUMP.InletVapour; |
---|
| 1405 | TRAYS(1).OutletVapour to CONDENSER.InletVapour; |
---|
[715] | 1406 | |
---|
| 1407 | #liquid |
---|
[735] | 1408 | CONDENSER.OutletLiquid to SPLITTER_TOP.Inlet; |
---|
[838] | 1409 | SPLITTER_TOP.Reflux to PUMP.Inlet; |
---|
[735] | 1410 | PUMP.Outlet to TRAYS(1).InletLiquid; |
---|
[810] | 1411 | TRAYS(NumberOfTrays).OutletLiquid to SUMP.InletLiquid; |
---|
| 1412 | SUMP.OutletLiquid to SPLITTER_BOTTOM.Inlet; |
---|
[735] | 1413 | SPLITTER_BOTTOM.Outlet2 to REBOILER.InletLiquid; |
---|
[715] | 1414 | |
---|
| 1415 | #Connectors |
---|
[721] | 1416 | ConnectorHeatCondenser to CONDENSER.InletQ; |
---|
| 1417 | ConnectorHeatReboiler to REBOILER.InletQ; |
---|
[810] | 1418 | CONDENSER.OutletVapour to ConnectorCondenserVout; |
---|
[838] | 1419 | SPLITTER_TOP.Distillate to ConnectorSplitterTop; |
---|
[721] | 1420 | SPLITTER_BOTTOM.Outlet1 to ConnectorSplitterBottom; |
---|
[715] | 1421 | |
---|
| 1422 | end |
---|
| 1423 | |
---|
[811] | 1424 | Model Distillation_kettle_subcooling as ColumnBasic |
---|
[715] | 1425 | |
---|
| 1426 | ATTRIBUTES |
---|
| 1427 | Pallete = true; |
---|
| 1428 | Icon = "icon/DistillationKettleSubcooling"; |
---|
| 1429 | Brief = "Model of a distillation column with steady condenser and dynamic reboiler."; |
---|
| 1430 | Info = |
---|
| 1431 | "== Specify == |
---|
| 1432 | * the feed stream of each tray (Inlet); |
---|
| 1433 | * the Murphree eficiency for each tray (Emv); |
---|
| 1434 | * the pump pressure difference; |
---|
| 1435 | * the heat supllied in reboiler and condenser; |
---|
| 1436 | * the heat supllied in the top tank; |
---|
| 1437 | * the condenser pressure drop; |
---|
[735] | 1438 | * the reboiler liquid outlet flow (OutletLiquid.F); |
---|
[715] | 1439 | * both splitter outlet flows OR one of the splitter outlet flows and the splitter frac. |
---|
| 1440 | |
---|
| 1441 | == Initial Conditions == |
---|
[735] | 1442 | * the TRAYS temperature (OutletLiquid.T); |
---|
| 1443 | * the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F); |
---|
| 1444 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray; |
---|
[715] | 1445 | |
---|
[735] | 1446 | * the top tank temperature (OutletLiquid.T); |
---|
[715] | 1447 | * the top tank liquid level (Level); |
---|
[735] | 1448 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions; |
---|
[715] | 1449 | |
---|
[735] | 1450 | * the reboiler temperature (OutletLiquid.T); |
---|
[715] | 1451 | * the reboiler liquid level (Level); |
---|
[735] | 1452 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions. |
---|
[715] | 1453 | "; |
---|
| 1454 | |
---|
[811] | 1455 | SET |
---|
| 1456 | |
---|
| 1457 | NumberOfFeeds = 1; |
---|
| 1458 | |
---|
[715] | 1459 | VARIABLES |
---|
[721] | 1460 | CONDENSER as condenserSteady (Brief="steady state condenser with subcooling"); |
---|
[811] | 1461 | REFLUX_DRUM as TankL (Brief="reflux drum"); |
---|
[838] | 1462 | SPLITTER as splitter_column (Brief="splitter to separate reflux and distillate",Protected=true); |
---|
[721] | 1463 | PUMP as pump (Brief="pump in reflux stream"); |
---|
| 1464 | REBOILER as reboiler (Brief="kettle reboiler"); |
---|
[811] | 1465 | |
---|
[838] | 1466 | RefluxRatio as positive (Brief="Reflux Ratio",Default=10, Lower = 0.5); |
---|
| 1467 | |
---|
[811] | 1468 | in Q_cond as power (Brief="Heat supplied to Condenser", PosX=1, PosY=0.065, Protected=true); |
---|
| 1469 | in Q_reb as power (Brief="Heat supplied to Reboiler", PosX=1, PosY=0.905, Protected=true); |
---|
| 1470 | in Q_drum as power (Brief="Heat supplied to Top Vessel", PosX=1, PosY=0.30, Protected=true); |
---|
| 1471 | |
---|
| 1472 | in FeedTray as stream (Brief="Feed stream", PosX=0, PosY=0.55); |
---|
| 1473 | |
---|
| 1474 | out TI_reb as control_signal (Brief="Temperature Indicator of Reboiler ",Protected = true, PosX=0.515, PosY=1); |
---|
| 1475 | out LI_reb as control_signal (Brief="Level Indicator of Reboiler", Protected = true, PosX=1, PosY=0.935); |
---|
| 1476 | out PI_reb as control_signal (Brief="Pressure Indicator of Reboiler", Protected = true, PosX=0.41, PosY=1); |
---|
| 1477 | |
---|
| 1478 | out TI_cond as control_signal (Brief="Temperature Indicator of Condenser ", Protected = true, PosX=0.52, PosY=0); |
---|
| 1479 | out PI_cond as control_signal (Brief="Pressure Indicator of Condenser", Protected = true, PosX=0.42, PosY=0); |
---|
| 1480 | |
---|
| 1481 | out TI_drum as control_signal (Brief="Temperature Indicator of Reflux Drum ", Protected = true, PosX=1, PosY=0.20); |
---|
| 1482 | out PI_drum as control_signal (Brief="Pressure Indicator of Reflux Drum", Protected = true, PosX=1, PosY=0.155); |
---|
| 1483 | out LI_drum as control_signal (Brief="Level Indicator of Reflux Drum", Protected = true, PosX=1, PosY=0.24); |
---|
| 1484 | |
---|
| 1485 | |
---|
[721] | 1486 | in ConnectorSplitterOut as stream (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true); |
---|
| 1487 | in ConnectorReboilerLout as stream (Brief="Connector for Liquid outlet stream From Reboiler", Hidden=true); |
---|
| 1488 | ConnectorHeatReboiler as power (Brief="Connector for Heat supplied to Reboiler", Hidden=true); |
---|
| 1489 | ConnectorHeatCondenser as power (Brief="Connector for Heat supplied to Condenser", Hidden=true); |
---|
| 1490 | ConnectorHeatAccumulator as power (Brief="Connector for Heat supplied to TopVessel", Hidden=true); |
---|
[811] | 1491 | |
---|
| 1492 | out LiquidDistillate as liquid_stream (Brief="Liquid outlet stream From Top Splitter", PosX=1, PosY=0.33); |
---|
[721] | 1493 | out BottomProduct as liquid_stream (Brief="Liquid outlet stream From Reboiler", PosX=0.67, PosY=1); |
---|
[715] | 1494 | |
---|
[721] | 1495 | EQUATIONS |
---|
[715] | 1496 | |
---|
[838] | 1497 | SPLITTER.RefluxRatio = RefluxRatio; |
---|
| 1498 | |
---|
[811] | 1499 | TI_reb = REBOILER.TI ; |
---|
| 1500 | PI_reb = REBOILER.PI; |
---|
| 1501 | LI_reb = REBOILER.LI; |
---|
| 1502 | |
---|
| 1503 | TI_cond = CONDENSER.TI; |
---|
| 1504 | PI_cond = CONDENSER.PI; |
---|
| 1505 | |
---|
| 1506 | TI_drum = REFLUX_DRUM.TI; |
---|
| 1507 | PI_drum = REFLUX_DRUM.PI; |
---|
| 1508 | LI_drum = REFLUX_DRUM.LI; |
---|
| 1509 | |
---|
| 1510 | for i in 1:NumberOfTrays do |
---|
| 1511 | |
---|
| 1512 | if i equal FeedTrayLocation(1) then |
---|
| 1513 | |
---|
| 1514 | "FeedTrayTop Inlet Flow" |
---|
| 1515 | FeedTray.F= TRAYS(i).Inlet.F; |
---|
| 1516 | |
---|
| 1517 | "FeedTrayTop Inlet Temperature" |
---|
| 1518 | FeedTray.T = TRAYS(i).Inlet.T; |
---|
| 1519 | |
---|
| 1520 | "FeedTrayTop Inlet Pressure" |
---|
| 1521 | FeedTray.P = TRAYS(i).Inlet.P; |
---|
| 1522 | |
---|
| 1523 | "FeedTrayTop Inlet Composition" |
---|
| 1524 | FeedTray.z = TRAYS(i).Inlet.z; |
---|
| 1525 | |
---|
| 1526 | "FeedTrayTop Inlet Vapour Fraction" |
---|
| 1527 | FeedTray.v = TRAYS(i).Inlet.v; |
---|
| 1528 | |
---|
| 1529 | "FeedTrayTop Inlet Enthalpy" |
---|
| 1530 | FeedTray.h = TRAYS(i).Inlet.h; |
---|
| 1531 | |
---|
| 1532 | else |
---|
| 1533 | |
---|
[832] | 1534 | "Inlet Tray - Flow Sealed" |
---|
| 1535 | 0*'mol/h'= TRAYS(i).Inlet.F; |
---|
| 1536 | |
---|
[811] | 1537 | "Inlet Tray - Temperature" |
---|
[832] | 1538 | 300*'K' = TRAYS(i).Inlet.T; |
---|
| 1539 | |
---|
[811] | 1540 | "Inlet Tray - Pressure" |
---|
[832] | 1541 | 1*'atm' = TRAYS(i).Inlet.P; |
---|
| 1542 | |
---|
[811] | 1543 | "Inlet Tray - Composition" |
---|
[832] | 1544 | 0.1 = TRAYS(i).Inlet.z; |
---|
| 1545 | |
---|
[811] | 1546 | "Inlet Tray - Vapour Fraction" |
---|
[832] | 1547 | 0 = TRAYS(i).Inlet.v; |
---|
| 1548 | |
---|
[811] | 1549 | "Inlet Tray - Enthalpy" |
---|
[832] | 1550 | 0*'J/mol' = TRAYS(i).Inlet.h; |
---|
[811] | 1551 | |
---|
| 1552 | end |
---|
| 1553 | |
---|
| 1554 | end |
---|
| 1555 | |
---|
[721] | 1556 | #Heat Connectors |
---|
[811] | 1557 | Q_cond = ConnectorHeatCondenser; |
---|
| 1558 | Q_reb = ConnectorHeatReboiler; |
---|
| 1559 | Q_drum = ConnectorHeatAccumulator; |
---|
[715] | 1560 | |
---|
| 1561 | # Splitter Connector Equations |
---|
| 1562 | ConnectorSplitterOut.T = LiquidDistillate.T; |
---|
| 1563 | ConnectorSplitterOut.P = LiquidDistillate.P; |
---|
| 1564 | ConnectorSplitterOut.F = LiquidDistillate.F; |
---|
| 1565 | ConnectorSplitterOut.z = LiquidDistillate.z; |
---|
| 1566 | |
---|
| 1567 | # Reboiler Connector Equations |
---|
| 1568 | ConnectorReboilerLout.T = BottomProduct.T; |
---|
| 1569 | ConnectorReboilerLout.P = BottomProduct.P; |
---|
| 1570 | ConnectorReboilerLout.F = BottomProduct.F; |
---|
| 1571 | ConnectorReboilerLout.z = BottomProduct.z; |
---|
| 1572 | |
---|
| 1573 | CONNECTIONS |
---|
| 1574 | #vapor |
---|
[735] | 1575 | REBOILER.OutletVapour to TRAYS(NumberOfTrays).InletVapour; |
---|
| 1576 | TRAYS(1).OutletVapour to CONDENSER.InletVapour; |
---|
[715] | 1577 | |
---|
| 1578 | #liquid |
---|
[811] | 1579 | CONDENSER.OutletLiquid to REFLUX_DRUM.Inlet; |
---|
| 1580 | REFLUX_DRUM.OutletLiquid to SPLITTER.Inlet; |
---|
[838] | 1581 | SPLITTER.Reflux to PUMP.Inlet; |
---|
[811] | 1582 | PUMP.Outlet to TRAYS(1).InletLiquid; |
---|
[735] | 1583 | TRAYS(NumberOfTrays).OutletLiquid to REBOILER.InletLiquid; |
---|
[715] | 1584 | |
---|
| 1585 | #Connectors |
---|
[811] | 1586 | ConnectorHeatCondenser to CONDENSER.InletQ; |
---|
| 1587 | ConnectorHeatReboiler to REBOILER.InletQ; |
---|
| 1588 | ConnectorHeatAccumulator to REFLUX_DRUM.InletQ; |
---|
[715] | 1589 | |
---|
[838] | 1590 | SPLITTER.Distillate to ConnectorSplitterOut; |
---|
[811] | 1591 | REBOILER.OutletLiquid to ConnectorReboilerLout; |
---|
[715] | 1592 | |
---|
| 1593 | end |
---|
| 1594 | |
---|
[815] | 1595 | Model Rectifier as ColumnBasic |
---|
[715] | 1596 | |
---|
| 1597 | ATTRIBUTES |
---|
| 1598 | Pallete = true; |
---|
| 1599 | Icon = "icon/RefluxedCond"; |
---|
| 1600 | Brief = "Model of a rectifier column with dynamic condenser."; |
---|
| 1601 | Info = |
---|
| 1602 | "== Specify == |
---|
| 1603 | * the feed stream of each tray (Inlet); |
---|
| 1604 | * the Murphree eficiency for each tray Emv; |
---|
[735] | 1605 | * the InletVapour stream of the bottom tray unless its flow; |
---|
[715] | 1606 | * the pump pressure difference; |
---|
| 1607 | * the heat supllied in the condenser; |
---|
[735] | 1608 | * the condenser vapor outlet flow (OutletVapour.F); |
---|
[715] | 1609 | * both splitter outlet flows OR one of the splitter outlet flows and the splitter frac. |
---|
| 1610 | |
---|
| 1611 | == Initial Conditions == |
---|
[735] | 1612 | * the TRAYS temperature (OutletLiquid.T); |
---|
| 1613 | * the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F); |
---|
| 1614 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray; |
---|
[715] | 1615 | |
---|
[735] | 1616 | * the condenser temperature (OutletLiquid.T); |
---|
[715] | 1617 | * the condenser liquid level (Level); |
---|
[735] | 1618 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions; |
---|
[715] | 1619 | "; |
---|
| 1620 | |
---|
[815] | 1621 | SET |
---|
| 1622 | NumberOfFeeds = 1; |
---|
[715] | 1623 | |
---|
| 1624 | VARIABLES |
---|
[721] | 1625 | CONDENSER as condenser (Brief="dymamic condenser without subcooling"); |
---|
[815] | 1626 | SPLITTER as splitter2 (Brief="splitter which separate reflux and distillate"); |
---|
[721] | 1627 | PUMP as pump (Brief="pump in reflux stream"); |
---|
[815] | 1628 | |
---|
| 1629 | out TI_cond as control_signal (Brief="Temperature Indicator of Condenser ", Protected = true, PosX=0.52, PosY=0); |
---|
| 1630 | out LI_cond as control_signal (Brief="Level Indicator of Condenser", Protected = true, PosX=1, PosY=0.084); |
---|
| 1631 | out PI_cond as control_signal (Brief="Pressure Indicator of Condenser", Protected = true, PosX=0.42, PosY=0); |
---|
| 1632 | |
---|
| 1633 | in FeedTray as stream (Brief="Feed stream", PosX=0, PosY=0.55); |
---|
| 1634 | |
---|
[721] | 1635 | out VapourDistillate as vapour_stream (Brief="Vapour outlet stream From Top Condenser", PosX=0.66, PosY=0); |
---|
| 1636 | out LiquidDistillate as liquid_stream (Brief="Liquid outlet stream From Top Splitter", PosX=1, PosY=0.24); |
---|
[815] | 1637 | |
---|
| 1638 | in Q_cond as power (Brief="Heat supplied to Condenser",PosX=1, PosY=0.04, Protected = true); |
---|
| 1639 | |
---|
[721] | 1640 | in VapourInlet as stream (Brief="Vapour Inlet in the section", PosX=0.07, PosY=1); |
---|
| 1641 | out LiquidOutlet as liquid_stream (Brief="Liquid Outlet in the section", PosX=0.32, PosY=1); |
---|
[815] | 1642 | |
---|
[721] | 1643 | in ConnectorSplitterOut as stream (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true); |
---|
| 1644 | in ConnectorCondenserVout as stream (Brief="Connector for Vapour outlet stream From Top Condenser", Hidden=true); |
---|
| 1645 | VapourConnector as stream (Brief="Vapour connection at the middle TRAYS", PosX=0.55, PosY=0,Hidden=true); |
---|
| 1646 | ConnectorHeatCondenser as power (Brief="Connector for Heat supplied to Condenser", Hidden=true); |
---|
[715] | 1647 | |
---|
| 1648 | EQUATIONS |
---|
| 1649 | |
---|
[815] | 1650 | TI_cond = CONDENSER.TI; |
---|
| 1651 | LI_cond = CONDENSER.LI; |
---|
| 1652 | PI_cond = CONDENSER.PI; |
---|
| 1653 | |
---|
| 1654 | for i in 1:NumberOfTrays do |
---|
[715] | 1655 | |
---|
[815] | 1656 | if i equal FeedTrayLocation(1) then |
---|
[715] | 1657 | |
---|
[815] | 1658 | "FeedTrayTop Inlet Flow" |
---|
| 1659 | FeedTray.F= TRAYS(i).Inlet.F; |
---|
[715] | 1660 | |
---|
[815] | 1661 | "FeedTrayTop Inlet Temperature" |
---|
| 1662 | FeedTray.T = TRAYS(i).Inlet.T; |
---|
[721] | 1663 | |
---|
[815] | 1664 | "FeedTrayTop Inlet Pressure" |
---|
| 1665 | FeedTray.P = TRAYS(i).Inlet.P; |
---|
| 1666 | |
---|
| 1667 | "FeedTrayTop Inlet Composition" |
---|
| 1668 | FeedTray.z = TRAYS(i).Inlet.z; |
---|
| 1669 | |
---|
| 1670 | "FeedTrayTop Inlet Vapour Fraction" |
---|
| 1671 | FeedTray.v = TRAYS(i).Inlet.v; |
---|
| 1672 | |
---|
| 1673 | "FeedTrayTop Inlet Enthalpy" |
---|
| 1674 | FeedTray.h = TRAYS(i).Inlet.h; |
---|
| 1675 | |
---|
| 1676 | else |
---|
| 1677 | |
---|
[832] | 1678 | "Inlet Tray - Flow Sealed" |
---|
| 1679 | 0*'mol/h'= TRAYS(i).Inlet.F; |
---|
| 1680 | |
---|
[815] | 1681 | "Inlet Tray - Temperature" |
---|
[832] | 1682 | 300*'K' = TRAYS(i).Inlet.T; |
---|
| 1683 | |
---|
[815] | 1684 | "Inlet Tray - Pressure" |
---|
[832] | 1685 | 1*'atm' = TRAYS(i).Inlet.P; |
---|
| 1686 | |
---|
[815] | 1687 | "Inlet Tray - Composition" |
---|
[832] | 1688 | 0.1 = TRAYS(i).Inlet.z; |
---|
| 1689 | |
---|
[815] | 1690 | "Inlet Tray - Vapour Fraction" |
---|
[832] | 1691 | 0 = TRAYS(i).Inlet.v; |
---|
| 1692 | |
---|
[815] | 1693 | "Inlet Tray - Enthalpy" |
---|
[832] | 1694 | 0*'J/mol' = TRAYS(i).Inlet.h; |
---|
[815] | 1695 | |
---|
| 1696 | end |
---|
| 1697 | |
---|
| 1698 | end |
---|
| 1699 | |
---|
| 1700 | Q_cond = ConnectorHeatCondenser; |
---|
| 1701 | |
---|
[715] | 1702 | # Condenser Connector Equations |
---|
| 1703 | ConnectorCondenserVout.T = VapourDistillate.T; |
---|
| 1704 | ConnectorCondenserVout.P = VapourDistillate.P; |
---|
| 1705 | ConnectorCondenserVout.F = VapourDistillate.F; |
---|
| 1706 | ConnectorCondenserVout.z = VapourDistillate.z; |
---|
| 1707 | |
---|
| 1708 | # Splitter Connector Equations |
---|
| 1709 | ConnectorSplitterOut.T = LiquidDistillate.T; |
---|
| 1710 | ConnectorSplitterOut.P = LiquidDistillate.P; |
---|
| 1711 | ConnectorSplitterOut.F = LiquidDistillate.F; |
---|
| 1712 | ConnectorSplitterOut.z = LiquidDistillate.z; |
---|
| 1713 | |
---|
[735] | 1714 | LiquidOutlet.F= TRAYS(NumberOfTrays).OutletLiquid.F; |
---|
| 1715 | LiquidOutlet.T = TRAYS(NumberOfTrays).OutletLiquid.T; |
---|
| 1716 | LiquidOutlet.P = TRAYS(NumberOfTrays).OutletLiquid.P; |
---|
| 1717 | LiquidOutlet.z = TRAYS(NumberOfTrays).OutletLiquid.z; |
---|
[715] | 1718 | |
---|
| 1719 | VapourConnector.F= VapourInlet.F; |
---|
| 1720 | VapourConnector.T = VapourInlet.T; |
---|
| 1721 | VapourConnector.P = VapourInlet.P; |
---|
| 1722 | VapourConnector.z = VapourInlet.z; |
---|
| 1723 | VapourConnector.v = VapourInlet.v; |
---|
| 1724 | VapourConnector.h = VapourInlet.h; |
---|
| 1725 | |
---|
| 1726 | CONNECTIONS |
---|
| 1727 | #vapor |
---|
[735] | 1728 | TRAYS(1).OutletVapour to CONDENSER.InletVapour; |
---|
[715] | 1729 | |
---|
| 1730 | #liquid |
---|
[735] | 1731 | CONDENSER.OutletLiquid to SPLITTER.Inlet; |
---|
[721] | 1732 | SPLITTER.Outlet2 to PUMP.Inlet; |
---|
[735] | 1733 | PUMP.Outlet to TRAYS(1).InletLiquid; |
---|
[715] | 1734 | |
---|
| 1735 | #Connectors |
---|
[721] | 1736 | ConnectorHeatCondenser to CONDENSER.InletQ; |
---|
[735] | 1737 | VapourConnector to TRAYS(NumberOfTrays).InletVapour; |
---|
[721] | 1738 | SPLITTER.Outlet1 to ConnectorSplitterOut; |
---|
[815] | 1739 | CONDENSER.OutletVapour to ConnectorCondenserVout; |
---|
[715] | 1740 | |
---|
| 1741 | end |
---|
| 1742 | |
---|
[815] | 1743 | #* |
---|
[787] | 1744 | Model Rectifier_subcooling as Section_ColumnBasicX |
---|
[715] | 1745 | |
---|
| 1746 | ATTRIBUTES |
---|
| 1747 | Pallete = true; |
---|
| 1748 | Icon = "icon/RefluxedSubcooling"; |
---|
| 1749 | Brief = "Model of a rectifier column with steady condenser."; |
---|
| 1750 | Info = |
---|
| 1751 | "== Specify == |
---|
| 1752 | * the feed stream of each tray (Inlet); |
---|
| 1753 | * the Murphree eficiency for each tray Emv; |
---|
[735] | 1754 | * the InletVapour stream of the bottom tray unless its flow; |
---|
[715] | 1755 | * the pump head; |
---|
| 1756 | * the condenser pressure drop; |
---|
| 1757 | * the heat supllied in the top tank; |
---|
| 1758 | * the heat supllied in condenser; |
---|
| 1759 | * both top splitter outlet flows OR one of the splitter outlet flows and the splitter frac. |
---|
| 1760 | |
---|
| 1761 | == Initial Conditions == |
---|
[735] | 1762 | * the TRAYS temperature (OutletLiquid.T); |
---|
| 1763 | * the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F); |
---|
| 1764 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray; |
---|
[715] | 1765 | |
---|
[735] | 1766 | * the top tank temperature (OutletLiquid.T); |
---|
[715] | 1767 | * the top tank liquid level (Level); |
---|
[735] | 1768 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions; |
---|
[715] | 1769 | "; |
---|
| 1770 | |
---|
| 1771 | PARAMETERS |
---|
| 1772 | CondenserVapourFlow as Switcher (Valid = ["on", "off"], Default = "on",Hidden=true); |
---|
| 1773 | |
---|
| 1774 | VARIABLES |
---|
[721] | 1775 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1776 | # Column Vapour and Liquid Draw Sides |
---|
| 1777 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1778 | out VapourDrawOff as vapour_stream (Brief="Vapour Outlet in the section", PosX=1, PosY=0.485,Protected = true); |
---|
| 1779 | out LiquidDrawOff as liquid_stream (Brief="Liquid Outlet in the section", PosX=1, PosY=0.535,Protected = true); |
---|
| 1780 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1781 | # Column Devices |
---|
| 1782 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1783 | CONDENSER as condenserSteady (Brief="steady state condenser with subcooling"); |
---|
| 1784 | ACCUMULATOR as tank_cylindrical (Brief="vessel drum (layed cilinder)"); |
---|
[739] | 1785 | SPLITTER as splitter2 (Brief="plitter which separate reflux and distillate"); |
---|
[721] | 1786 | PUMP as pump (Brief="pump in reflux stream"); |
---|
| 1787 | alfaTopo as Real; |
---|
| 1788 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1789 | # Column Bottom outlets |
---|
| 1790 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1791 | in VapourInlet as stream (Brief="Vapour Inlet in the section", PosX=0.07, PosY=1); |
---|
| 1792 | out LiquidOutlet as liquid_stream (Brief="Liquid Outlet in the section", PosX=0.32, PosY=1); |
---|
| 1793 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1794 | # Heat Ports |
---|
| 1795 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1796 | in HeatToCondenser as power (Brief="Heat supplied to Condenser", PosX=1, PosY=0.070, Protected=true); |
---|
| 1797 | in HeatToAccumulator as power (Brief="Heat supplied to Top Accumulator", PosX=1, PosY=0.23, Protected=true); |
---|
| 1798 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1799 | # Column Top Product |
---|
| 1800 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1801 | out LiquidDistillate as liquid_stream (Brief="Liquid outlet stream From Top Splitter", PosX=1, PosY=0.285); |
---|
| 1802 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1803 | # Column Connectors |
---|
| 1804 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1805 | in ConnectorSplitterTop as stream (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true); |
---|
| 1806 | VapourConnector as stream (Brief="Vapour connection at the middle TRAYS", PosX=0.55, PosY=0,Hidden=true); |
---|
| 1807 | ConnectorHeatCondenser as power (Brief="Connector for Heat supplied to Condenser", Hidden=true); |
---|
| 1808 | ConnectorHeatAccumulator as power (Brief="Connector for Heat supplied to Top Accumulator", Hidden=true); |
---|
[715] | 1809 | |
---|
| 1810 | EQUATIONS |
---|
[721] | 1811 | #Heat Connectors |
---|
| 1812 | HeatToCondenser = ConnectorHeatCondenser; |
---|
| 1813 | HeatToAccumulator = ConnectorHeatAccumulator; |
---|
| 1814 | |
---|
[735] | 1815 | LiquidOutlet.F= TRAYS(NumberOfTrays).OutletLiquid.F; |
---|
| 1816 | LiquidOutlet.T = TRAYS(NumberOfTrays).OutletLiquid.T; |
---|
| 1817 | LiquidOutlet.P = TRAYS(NumberOfTrays).OutletLiquid.P; |
---|
| 1818 | LiquidOutlet.z = TRAYS(NumberOfTrays).OutletLiquid.z; |
---|
[715] | 1819 | |
---|
| 1820 | VapourConnector.F= VapourInlet.F; |
---|
| 1821 | VapourConnector.T = VapourInlet.T; |
---|
| 1822 | VapourConnector.P = VapourInlet.P; |
---|
| 1823 | VapourConnector.z = VapourInlet.z; |
---|
| 1824 | VapourConnector.v = VapourInlet.v; |
---|
| 1825 | VapourConnector.h = VapourInlet.h; |
---|
| 1826 | |
---|
| 1827 | # Splitter Connector Equations |
---|
| 1828 | ConnectorSplitterTop.T = LiquidDistillate.T; |
---|
| 1829 | ConnectorSplitterTop.P = LiquidDistillate.P; |
---|
| 1830 | ConnectorSplitterTop.F = LiquidDistillate.F; |
---|
| 1831 | ConnectorSplitterTop.z = LiquidDistillate.z; |
---|
| 1832 | |
---|
[729] | 1833 | VapourDrawOff.F*VapSideTrayIndex= TRAYS.VapourSideStream.F; |
---|
| 1834 | VapourDrawOff.T = TRAYS(VapourSideStreamLocation).VapourSideStream.T; |
---|
| 1835 | VapourDrawOff.P = TRAYS(VapourSideStreamLocation).VapourSideStream.P; |
---|
| 1836 | VapourDrawOff.z = TRAYS(VapourSideStreamLocation).VapourSideStream.z; |
---|
| 1837 | |
---|
| 1838 | LiquidDrawOff.F*LiqSideTrayIndex= TRAYS.LiquidSideStream.F; |
---|
| 1839 | LiquidDrawOff.T = TRAYS(LiquidSideStreamLocation).LiquidSideStream.T; |
---|
| 1840 | LiquidDrawOff.P = TRAYS(LiquidSideStreamLocation).LiquidSideStream.P; |
---|
| 1841 | LiquidDrawOff.z = TRAYS(LiquidSideStreamLocation).LiquidSideStream.z; |
---|
| 1842 | |
---|
| 1843 | VapourDrawOffFlow = VapourDrawOff.F; |
---|
| 1844 | LiquidDrawOffFlow = LiquidDrawOff.F; |
---|
| 1845 | |
---|
[715] | 1846 | switch CondenserVapourFlow |
---|
| 1847 | |
---|
| 1848 | case "on": |
---|
[735] | 1849 | CONDENSER.InletVapour.F*TRAYS(1).vV = alfaTopo *Ah * sqrt(2*(TRAYS(1).OutletVapour.P - |
---|
| 1850 | CONDENSER.OutletLiquid.P + 1e-8 * 'atm') / (alfa*TRAYS(1).rhoV)); |
---|
| 1851 | when CONDENSER.InletVapour.F < 1e-6 * 'kmol/h' switchto "off"; |
---|
[715] | 1852 | |
---|
| 1853 | case "off": |
---|
[735] | 1854 | CONDENSER.InletVapour.F = 0 * 'mol/s'; |
---|
| 1855 | when TRAYS(1).OutletVapour.P > CONDENSER.OutletLiquid.P + 1e-1 * 'atm' switchto "on"; |
---|
[715] | 1856 | |
---|
| 1857 | end |
---|
| 1858 | |
---|
| 1859 | CONNECTIONS |
---|
| 1860 | #vapor |
---|
[735] | 1861 | TRAYS(1).OutletVapour to CONDENSER.InletVapour; |
---|
[715] | 1862 | |
---|
| 1863 | #liquid |
---|
[735] | 1864 | CONDENSER.OutletLiquid to ACCUMULATOR.Inlet; |
---|
[721] | 1865 | ACCUMULATOR.Outlet to SPLITTER.Inlet; |
---|
| 1866 | SPLITTER.Outlet2 to PUMP.Inlet; |
---|
[735] | 1867 | PUMP.Outlet to TRAYS(1).InletLiquid; |
---|
[715] | 1868 | |
---|
| 1869 | #Connectors |
---|
[735] | 1870 | VapourConnector to TRAYS(NumberOfTrays).InletVapour; |
---|
[721] | 1871 | ConnectorHeatCondenser to CONDENSER.InletQ; |
---|
| 1872 | SPLITTER.Outlet1 to ConnectorSplitterTop; |
---|
| 1873 | ConnectorHeatAccumulator to ACCUMULATOR.InletQ; |
---|
[715] | 1874 | |
---|
| 1875 | end |
---|
| 1876 | |
---|
[787] | 1877 | Model Reboiled_Stripping_kettle as Section_ColumnBasicX |
---|
[715] | 1878 | |
---|
| 1879 | ATTRIBUTES |
---|
| 1880 | Pallete = true; |
---|
| 1881 | Icon = "icon/ReboiledKettle"; |
---|
| 1882 | Brief = "Model of a reboiled stripping column with dynamic reboiler."; |
---|
| 1883 | Info = |
---|
| 1884 | "== Specify == |
---|
| 1885 | * the feed stream of each tray (Inlet); |
---|
| 1886 | * the Murphree eficiency for each tray Emv; |
---|
| 1887 | * the vapour flow leaving the top of the column; |
---|
[735] | 1888 | * the InletLiquidiquid stream of the top tray; |
---|
[715] | 1889 | * the heat supllied in the reboiler; |
---|
[735] | 1890 | * the reboiler liquid outlet flow (OutletLiquid.F); |
---|
[715] | 1891 | |
---|
| 1892 | == Initial Conditions == |
---|
[735] | 1893 | * the TRAYS temperature (OutletLiquid.T); |
---|
| 1894 | * the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F); |
---|
| 1895 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray; |
---|
[715] | 1896 | |
---|
[735] | 1897 | * the reboiler temperature (OutletLiquid.T); |
---|
[715] | 1898 | * the reboiler liquid level (Level); |
---|
[735] | 1899 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions. |
---|
[715] | 1900 | "; |
---|
| 1901 | |
---|
| 1902 | VARIABLES |
---|
[721] | 1903 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1904 | # Column Vapour and Liquid Draw Sides |
---|
| 1905 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1906 | out VapourDrawOff as vapour_stream (Brief="Vapour Outlet in the section", PosX=1, PosY=0.28,Protected = true); |
---|
| 1907 | out LiquidDrawOff as liquid_stream (Brief="Liquid Outlet in the section", PosX=1, PosY=0.33,Protected = true); |
---|
| 1908 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1909 | # Column Devices |
---|
| 1910 | #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1911 | REBOILER as reboiler (Brief="Kettle Reboiler"); |
---|
| 1912 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1913 | # Heat Port |
---|
| 1914 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1915 | in HeatToReboiler as power (Brief="Heat supplied to Reboiler",PosX=1, PosY=0.865,Protected=true); |
---|
| 1916 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1917 | # Column Top Outlets |
---|
| 1918 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1919 | in LiquidInlet as stream (Brief="Liquid Inlet in the section", PosX=0.30, PosY=0); |
---|
| 1920 | out VapourOutlet as vapour_stream (Brief="Vapour Outlet in the section", PosX=0.07, PosY=0); |
---|
| 1921 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1922 | # Column Bottom Product |
---|
| 1923 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1924 | out BottomProduct as liquid_stream (Brief="Liquid outlet stream From Reboiler", PosX=0.68, PosY=1); |
---|
| 1925 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1926 | # Column Connectors |
---|
| 1927 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 1928 | in ConnectorReboilerLout as stream (Brief="Connector for Liquid outlet stream From Reboiler", Hidden=true); |
---|
| 1929 | LiquidConnector as stream (Brief="Liquid connection at the middle TRAYS", PosX=0.75, PosY=1,Hidden=true); |
---|
| 1930 | ConnectorHeatReboiler as power (Brief="Connector for Heat supplied to Reboiler", Hidden=true); |
---|
[715] | 1931 | |
---|
| 1932 | CONNECTIONS |
---|
| 1933 | #vapor |
---|
[735] | 1934 | REBOILER.OutletVapour to TRAYS(NumberOfTrays).InletVapour; |
---|
[715] | 1935 | |
---|
| 1936 | #liquid |
---|
[735] | 1937 | TRAYS(NumberOfTrays).OutletLiquid to REBOILER.InletLiquid; |
---|
[715] | 1938 | |
---|
| 1939 | #Connectors |
---|
[735] | 1940 | REBOILER.OutletLiquid to ConnectorReboilerLout; |
---|
| 1941 | LiquidConnector to TRAYS(1).InletLiquid; |
---|
[721] | 1942 | ConnectorHeatReboiler to REBOILER.InletQ; |
---|
[715] | 1943 | |
---|
| 1944 | EQUATIONS |
---|
[721] | 1945 | #Heat Connectors |
---|
| 1946 | HeatToReboiler = ConnectorHeatReboiler; |
---|
| 1947 | |
---|
[715] | 1948 | LiquidConnector.F= LiquidInlet.F; |
---|
| 1949 | LiquidConnector.T = LiquidInlet.T; |
---|
| 1950 | LiquidConnector.P = LiquidInlet.P; |
---|
| 1951 | LiquidConnector.z = LiquidInlet.z; |
---|
| 1952 | LiquidConnector.v = LiquidInlet.v; |
---|
| 1953 | LiquidConnector.h = LiquidInlet.h; |
---|
| 1954 | |
---|
[735] | 1955 | VapourOutlet.F= TRAYS(1).OutletVapour.F; |
---|
| 1956 | VapourOutlet.T = TRAYS(1).OutletVapour.T; |
---|
| 1957 | VapourOutlet.P = TRAYS(1).OutletVapour.P; |
---|
| 1958 | VapourOutlet.z = TRAYS(1).OutletVapour.z; |
---|
[715] | 1959 | |
---|
| 1960 | # Reboiler Connector Equations |
---|
| 1961 | ConnectorReboilerLout.T = BottomProduct.T; |
---|
| 1962 | ConnectorReboilerLout.P = BottomProduct.P; |
---|
| 1963 | ConnectorReboilerLout.F = BottomProduct.F; |
---|
| 1964 | ConnectorReboilerLout.z = BottomProduct.z; |
---|
| 1965 | |
---|
| 1966 | VapourDrawOff.F*VapSideTrayIndex= TRAYS.VapourSideStream.F; |
---|
| 1967 | VapourDrawOff.T = TRAYS(VapourSideStreamLocation).VapourSideStream.T; |
---|
| 1968 | VapourDrawOff.P = TRAYS(VapourSideStreamLocation).VapourSideStream.P; |
---|
| 1969 | VapourDrawOff.z = TRAYS(VapourSideStreamLocation).VapourSideStream.z; |
---|
| 1970 | |
---|
| 1971 | LiquidDrawOff.F*LiqSideTrayIndex= TRAYS.LiquidSideStream.F; |
---|
| 1972 | LiquidDrawOff.T = TRAYS(LiquidSideStreamLocation).LiquidSideStream.T; |
---|
| 1973 | LiquidDrawOff.P = TRAYS(LiquidSideStreamLocation).LiquidSideStream.P; |
---|
| 1974 | LiquidDrawOff.z = TRAYS(LiquidSideStreamLocation).LiquidSideStream.z; |
---|
| 1975 | |
---|
| 1976 | VapourDrawOffFlow = VapourDrawOff.F; |
---|
| 1977 | LiquidDrawOffFlow = LiquidDrawOff.F; |
---|
| 1978 | |
---|
| 1979 | end |
---|
| 1980 | |
---|
[787] | 1981 | Model Reboiled_Stripping_thermosyphon as Section_ColumnBasicX |
---|
[715] | 1982 | |
---|
| 1983 | ATTRIBUTES |
---|
| 1984 | Pallete = true; |
---|
| 1985 | Icon = "icon/ReboiledThermosyphon"; |
---|
| 1986 | Brief = "Model of a reboiled stripping column with steady reboiler."; |
---|
| 1987 | Info = |
---|
| 1988 | "== Specify == |
---|
| 1989 | * the feed stream of each tray (Inlet); |
---|
| 1990 | * the Murphree eficiency for each tray (Emv); |
---|
| 1991 | * the vapour flow leaving the top of the column; |
---|
[735] | 1992 | * the InletLiquidiquid stream of the top tray; |
---|
[715] | 1993 | * the heat supllied in bottom tank; |
---|
| 1994 | * the heat supllied in the reboiler; |
---|
| 1995 | * the Outlet1 flow in the bottom splitter (spbottom.Outlet1.F) that corresponds to the bottom product; |
---|
| 1996 | |
---|
| 1997 | == Initial Conditions == |
---|
[735] | 1998 | * the TRAYS temperature (OutletLiquid.T); |
---|
| 1999 | * the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F); |
---|
| 2000 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray; |
---|
[715] | 2001 | |
---|
[735] | 2002 | * the bottom tank temperature (OutletLiquid.T); |
---|
[715] | 2003 | * the bottom tank liquid level (Level); |
---|
[735] | 2004 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions. |
---|
[715] | 2005 | "; |
---|
[721] | 2006 | |
---|
[715] | 2007 | VARIABLES |
---|
[721] | 2008 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2009 | # Column Vapour and Liquid Draw Sides |
---|
| 2010 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2011 | out VapourDrawOff as vapour_stream (Brief="Vapour Outlet in the section", PosX=1, PosY=0.28,Protected = true); |
---|
| 2012 | out LiquidDrawOff as liquid_stream (Brief="Liquid Outlet in the section", PosX=1, PosY=0.33,Protected = true); |
---|
| 2013 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2014 | # Column Devices |
---|
| 2015 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2016 | REBOILER as reboilerSteady (Brief="steady state reboiler (thermosyphon)"); |
---|
[739] | 2017 | SPLITTER as splitter2 (Brief="splitter which separate the bottom product and the stream to reboiler"); |
---|
[721] | 2018 | ACCUMULATOR as tank (Brief="vessel in the bottom of column"); |
---|
| 2019 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2020 | # Heat Ports |
---|
| 2021 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2022 | in HeatToReboiler as power (Brief="Heat supplied to Reboiler",PosX=1, PosY=0.84,Protected=true); |
---|
| 2023 | in HeatToAccumulator as power (Brief="Heat supplied to Bottom Vessel",PosX=0, PosY=0.91,Protected=true); |
---|
| 2024 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2025 | # Column Top Outlets |
---|
| 2026 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2027 | in LiquidInlet as stream (Brief="Liquid Inlet in the section", PosX=0.32, PosY=0); |
---|
| 2028 | out VapourOutlet as vapour_stream (Brief="Vapour Outlet in the section", PosX=0.07, PosY=0); |
---|
| 2029 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2030 | # Column Bottom Product |
---|
| 2031 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2032 | out BottomProduct as liquid_stream (Brief="Liquid outlet stream From Reboiler", PosX=0.19, PosY=1); |
---|
| 2033 | LiquidConnector as stream (Brief="Liquid connection at the middle TRAYS", PosX=0.75, PosY=1,Hidden=true); |
---|
| 2034 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2035 | # Column Connectors |
---|
| 2036 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2037 | in ConnectorSplitterBottom as stream (Brief="Connector for Liquid outlet stream From Reboiler", Hidden=true); |
---|
| 2038 | ConnectorHeatReboiler as power (Brief="Connector for Heat supplied to Reboiler", Hidden=true); |
---|
| 2039 | ConnectorHeatAccumulator as power (Brief="Connector for Heat supplied to Bottom Accumulator", Hidden=true); |
---|
[715] | 2040 | |
---|
| 2041 | CONNECTIONS |
---|
| 2042 | #vapor |
---|
[735] | 2043 | REBOILER.OutletVapour to TRAYS(NumberOfTrays).InletVapour; |
---|
[715] | 2044 | |
---|
| 2045 | #liquid |
---|
[735] | 2046 | TRAYS(NumberOfTrays).OutletLiquid to ACCUMULATOR.Inlet; |
---|
[721] | 2047 | ACCUMULATOR.Outlet to SPLITTER.Inlet; |
---|
[735] | 2048 | SPLITTER.Outlet2 to REBOILER.InletLiquid; |
---|
[715] | 2049 | |
---|
| 2050 | #Connectors |
---|
[735] | 2051 | LiquidConnector to TRAYS(1).InletLiquid; |
---|
[721] | 2052 | ConnectorHeatReboiler to REBOILER.InletQ; |
---|
| 2053 | ConnectorHeatAccumulator to ACCUMULATOR.InletQ; |
---|
| 2054 | SPLITTER.Outlet1 to ConnectorSplitterBottom; |
---|
[715] | 2055 | |
---|
| 2056 | EQUATIONS |
---|
[721] | 2057 | #Heat Connectors |
---|
| 2058 | HeatToReboiler = ConnectorHeatReboiler; |
---|
| 2059 | HeatToAccumulator = ConnectorHeatAccumulator; |
---|
| 2060 | |
---|
[715] | 2061 | ConnectorSplitterBottom.T = BottomProduct.T; |
---|
| 2062 | ConnectorSplitterBottom.P = BottomProduct.P; |
---|
| 2063 | ConnectorSplitterBottom.F = BottomProduct.F; |
---|
| 2064 | ConnectorSplitterBottom.z = BottomProduct.z; |
---|
| 2065 | |
---|
| 2066 | LiquidConnector.F= LiquidInlet.F; |
---|
| 2067 | LiquidConnector.T = LiquidInlet.T; |
---|
| 2068 | LiquidConnector.P = LiquidInlet.P; |
---|
| 2069 | LiquidConnector.z = LiquidInlet.z; |
---|
| 2070 | LiquidConnector.v = LiquidInlet.v; |
---|
| 2071 | LiquidConnector.h = LiquidInlet.h; |
---|
| 2072 | |
---|
[735] | 2073 | VapourOutlet.F= TRAYS(1).OutletVapour.F; |
---|
| 2074 | VapourOutlet.T = TRAYS(1).OutletVapour.T; |
---|
| 2075 | VapourOutlet.P = TRAYS(1).OutletVapour.P; |
---|
| 2076 | VapourOutlet.z = TRAYS(1).OutletVapour.z; |
---|
[715] | 2077 | |
---|
| 2078 | VapourDrawOff.F*VapSideTrayIndex= TRAYS.VapourSideStream.F; |
---|
| 2079 | VapourDrawOff.T = TRAYS(VapourSideStreamLocation).VapourSideStream.T; |
---|
| 2080 | VapourDrawOff.P = TRAYS(VapourSideStreamLocation).VapourSideStream.P; |
---|
| 2081 | VapourDrawOff.z = TRAYS(VapourSideStreamLocation).VapourSideStream.z; |
---|
| 2082 | |
---|
| 2083 | LiquidDrawOff.F*LiqSideTrayIndex= TRAYS.LiquidSideStream.F; |
---|
| 2084 | LiquidDrawOff.T = TRAYS(LiquidSideStreamLocation).LiquidSideStream.T; |
---|
| 2085 | LiquidDrawOff.P = TRAYS(LiquidSideStreamLocation).LiquidSideStream.P; |
---|
| 2086 | LiquidDrawOff.z = TRAYS(LiquidSideStreamLocation).LiquidSideStream.z; |
---|
| 2087 | |
---|
| 2088 | VapourDrawOffFlow = VapourDrawOff.F; |
---|
| 2089 | LiquidDrawOffFlow = LiquidDrawOff.F; |
---|
[721] | 2090 | |
---|
[715] | 2091 | end |
---|
| 2092 | |
---|
| 2093 | Model ReactiveDistillation |
---|
| 2094 | |
---|
| 2095 | ATTRIBUTES |
---|
| 2096 | Pallete = true; |
---|
| 2097 | Icon = "icon/DistillationReac"; |
---|
| 2098 | Brief = "Model of a reactive distillation column with dynamic condenser and reboiler."; |
---|
| 2099 | Info = |
---|
| 2100 | "== Specify == |
---|
| 2101 | * the reaction related variables for each tray, condenser and reboiler; |
---|
| 2102 | * the feed stream of each tray (Inlet); |
---|
| 2103 | * the Murphree eficiency for each tray Emv; |
---|
| 2104 | * the pump pressure difference; |
---|
| 2105 | * the heat supllied in reboiler and condenser; |
---|
[735] | 2106 | * the condenser vapor outlet flow (OutletVapour.F); |
---|
| 2107 | * the reboiler liquid outlet flow (OutletLiquid.F); |
---|
[715] | 2108 | * both splitter outlet flows OR one of the splitter outlet flows and the splitter frac. |
---|
| 2109 | |
---|
| 2110 | == Initial Conditions == |
---|
[735] | 2111 | * the TRAYS temperature (OutletLiquid.T); |
---|
| 2112 | * the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F); |
---|
| 2113 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray; |
---|
[715] | 2114 | |
---|
[735] | 2115 | * the condenser temperature (OutletLiquid.T); |
---|
[715] | 2116 | * the condenser liquid level (Level); |
---|
[735] | 2117 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions; |
---|
[715] | 2118 | |
---|
[735] | 2119 | * the reboiler temperature (OutletLiquid.T); |
---|
[715] | 2120 | * the reboiler liquid level (Level); |
---|
[735] | 2121 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions. |
---|
[715] | 2122 | "; |
---|
| 2123 | |
---|
| 2124 | PARAMETERS |
---|
[721] | 2125 | outer PP as Plugin (Type="PP"); |
---|
| 2126 | outer NComp as Integer (Brief="Number of Components"); |
---|
| 2127 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2128 | # Section Column Feed Tray - Side Streams Location and Numbering |
---|
| 2129 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2130 | NumberOfTrays as Integer (Brief="Number of TRAYS", Default=2); |
---|
| 2131 | FeedTrayIndex(NumberOfTrays) as Integer (Brief="Feed Tray Index", Default=0,Hidden=true); |
---|
| 2132 | LiqSideTrayIndex(NumberOfTrays) as Integer (Brief="Liquid Side Tray Index", Default=0,Hidden=true); |
---|
| 2133 | VapSideTrayIndex(NumberOfTrays) as Integer (Brief="Vapour Side Tray Index", Default=0,Hidden=true); |
---|
| 2134 | FeedTrayLocation as Integer (Brief="Feed tray Location", Default=2); |
---|
| 2135 | LiquidSideStreamLocation as Integer (Brief="Liquid Side Stream Location", Default=2); |
---|
| 2136 | VapourSideStreamLocation as Integer (Brief="Vapour Side Stream Location", Default=2); |
---|
| 2137 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2138 | # Section Column Flow Model for Liquid and Vapour |
---|
| 2139 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2140 | VapourFlow as Switcher (Valid = ["on", "off"], Default = "off",Hidden=true); |
---|
[715] | 2141 | TrayVapourFlow as Switcher (Valid = ["on", "off"], Default = "off",Hidden=true); |
---|
[721] | 2142 | TrayLiquidFlow as Switcher (Valid = ["on", "off"], Default = "off",Hidden=true); |
---|
| 2143 | alfacond as Real; |
---|
| 2144 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2145 | # Section Column Tray Geometry and Auxiliar Parameters |
---|
| 2146 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2147 | V as volume (Brief="Total Volume of the tray",Hidden=true); |
---|
| 2148 | Q as power (Brief="Rate of heat supply",Hidden=true); |
---|
| 2149 | Ap as area (Brief="Plate area = Atray - Adowncomer",Hidden=true); |
---|
| 2150 | Ah as area (Brief="Total holes area",Hidden=true); |
---|
| 2151 | lw as length (Brief="Weir length",Hidden=true); |
---|
| 2152 | g as acceleration (Brief="Gravity Acceleration",Default=9.81,Hidden=true); |
---|
| 2153 | hw as length (Brief="Weir height",Hidden=true); |
---|
| 2154 | beta as fraction (Brief="Aeration fraction"); |
---|
| 2155 | alfa as fraction (Brief="Dry pressure drop coefficient"); |
---|
[715] | 2156 | |
---|
| 2157 | VolumeOfTray as volume (Brief="Total Volume of the tray"); |
---|
| 2158 | HeatSupply as heat_rate (Brief="Rate of heat supply"); |
---|
[721] | 2159 | PlateArea as area (Brief="Plate area = Atray - Adowncomer"); |
---|
| 2160 | HolesArea as area (Brief="Total holes area"); |
---|
| 2161 | WeirLength as length (Brief="Weir length"); |
---|
| 2162 | WeirHeight as length (Brief="Weir height"); |
---|
[715] | 2163 | |
---|
| 2164 | stoic(NComp) as Real(Brief="Stoichiometric matrix"); |
---|
| 2165 | Hr as energy_mol; |
---|
| 2166 | |
---|
| 2167 | SET |
---|
| 2168 | |
---|
[721] | 2169 | FeedTrayIndex(FeedTrayLocation) =1; |
---|
[715] | 2170 | VapSideTrayIndex(FeedTrayLocation) =1; |
---|
| 2171 | LiqSideTrayIndex(FeedTrayLocation) =1; |
---|
| 2172 | |
---|
| 2173 | V=VolumeOfTray; |
---|
| 2174 | Q=HeatSupply; |
---|
| 2175 | Ap=PlateArea; |
---|
| 2176 | Ah=HolesArea; |
---|
| 2177 | lw=WeirLength; |
---|
| 2178 | hw=WeirHeight ; |
---|
| 2179 | |
---|
| 2180 | VARIABLES |
---|
[721] | 2181 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2182 | # Tray Initialization For the Whole Column Section |
---|
| 2183 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2184 | INITIALIZATION as InitializeSection (Brief = "Column Model Initialization"); |
---|
| 2185 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2186 | # Column Feed - Side Stream Flow and Murphree Efficiency |
---|
| 2187 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2188 | in FeedTray as stream (Brief="Feed stream", PosX=0, PosY=0.55); |
---|
| 2189 | VapourDrawOffFlow as flow_mol (Brief = "Stream Molar Flow Rate"); |
---|
| 2190 | LiquidDrawOffFlow as flow_mol (Brief = "Stream Molar Flow Rate"); |
---|
| 2191 | MurphreeEff as Real (Brief = "Murphree efficiency"); |
---|
| 2192 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2193 | # Column Vapour and Liquid Draw Sides |
---|
| 2194 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2195 | out VapourDrawOff as vapour_stream (Brief="Vapour Outlet in the section", PosX=1, PosY=0.395,Protected = true); |
---|
| 2196 | out LiquidDrawOff as liquid_stream (Brief="Liquid Outlet in the section", PosX=1, PosY=0.44,Protected = true); |
---|
| 2197 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2198 | # Heat Ports |
---|
| 2199 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2200 | in HeatToReboiler as power (Brief="Heat supplied to Reboiler", PosX=1, PosY=0.89); |
---|
| 2201 | in HeatToCondenser as power (Brief="Heat supplied to Condenser", PosX=1, PosY=0.07); |
---|
| 2202 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2203 | # Column Devices |
---|
| 2204 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
[715] | 2205 | TRAYS(NumberOfTrays) as trayReac; |
---|
[721] | 2206 | CONDENSER as condenserReact; |
---|
| 2207 | REBOILER as reboilerReact; |
---|
[739] | 2208 | SPLITTER as splitter2; |
---|
[721] | 2209 | PUMP as pump; |
---|
| 2210 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2211 | # Column Outlets |
---|
| 2212 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2213 | out VapourDistillate as vapour_stream (Brief="Vapour outlet stream From Condenser", PosX=0.59, PosY=0); |
---|
| 2214 | out LiquidDistillate as liquid_stream (Brief="Liquid outlet stream From Splitter", PosX=1, PosY=0.23); |
---|
| 2215 | out BottomProduct as liquid_stream (Brief="Liquid outlet stream From Reboiler", PosX=0.67, PosY=1); |
---|
| 2216 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2217 | # Column Connectors |
---|
| 2218 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2219 | in ConnectorCondenserVout as stream (Brief="Connector for Vapour outlet stream From Top Condenser", Hidden=true); |
---|
| 2220 | in ConnectorSplitterOut as stream (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true); |
---|
| 2221 | in ConnectorReboilerLout as stream (Brief="Connector for Liquid outlet stream From Reboiler", Hidden=true); |
---|
| 2222 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2223 | # Column Control Ports |
---|
| 2224 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2225 | out TRI as control_signal (Brief="Temperature of Reboiler Indicator", PosX=1, PosY=0.83); |
---|
| 2226 | out LRI as control_signal (Brief="Level of Reboiler Indicator", PosX=1, PosY=0.855); |
---|
| 2227 | out TCI as control_signal (Brief="Temperature of Condenser Indicator", PosX=1, PosY=0.01); |
---|
| 2228 | out LCI as control_signal (Brief="Level of Condenser Indicator", PosX=1, PosY=0.04); |
---|
[715] | 2229 | |
---|
| 2230 | CONNECTIONS |
---|
| 2231 | #vapor |
---|
[735] | 2232 | REBOILER.OutletVapour to TRAYS(NumberOfTrays).InletVapour; |
---|
| 2233 | TRAYS(1).OutletVapour to CONDENSER.InletVapour; |
---|
| 2234 | TRAYS([2:NumberOfTrays]).OutletVapour to TRAYS([1:NumberOfTrays-1]).InletVapour; |
---|
[715] | 2235 | |
---|
| 2236 | #liquid |
---|
[735] | 2237 | CONDENSER.OutletLiquid to SPLITTER.Inlet; |
---|
[721] | 2238 | SPLITTER.Outlet2 to PUMP.Inlet; |
---|
[735] | 2239 | PUMP.Outlet to TRAYS(1).InletLiquid; |
---|
| 2240 | TRAYS([1:NumberOfTrays-1]).OutletLiquid to TRAYS([2:NumberOfTrays]).InletLiquid; |
---|
| 2241 | TRAYS(NumberOfTrays).OutletLiquid to REBOILER.InletLiquid; |
---|
[715] | 2242 | |
---|
| 2243 | #Connectors |
---|
[735] | 2244 | CONDENSER.OutletVapour to ConnectorCondenserVout; |
---|
[721] | 2245 | SPLITTER.Outlet1 to ConnectorSplitterOut; |
---|
[735] | 2246 | REBOILER.OutletLiquid to ConnectorReboilerLout; |
---|
[715] | 2247 | |
---|
[721] | 2248 | INITIAL |
---|
| 2249 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2250 | # Tray Initialization |
---|
| 2251 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2252 | for i in 1:NumberOfTrays do |
---|
[715] | 2253 | |
---|
[721] | 2254 | "The initial temperature of the TRAYS" |
---|
[735] | 2255 | TRAYS(i).OutletLiquid.T = INITIALIZATION.TopTemperature+(INITIALIZATION.BottomTemperature-INITIALIZATION.TopTemperature)*((i-1)/(NumberOfTrays-1)); |
---|
[721] | 2256 | |
---|
| 2257 | "The initial Level of the TRAYS" |
---|
| 2258 | TRAYS(i).Level = INITIALIZATION.LevelFraction*hw; |
---|
| 2259 | end |
---|
| 2260 | |
---|
| 2261 | for i in 1:NComp-1 do |
---|
| 2262 | |
---|
| 2263 | for j in 1:NumberOfTrays do |
---|
| 2264 | |
---|
| 2265 | "The initial composition of the TRAYS - Normalized" |
---|
[735] | 2266 | TRAYS(j).OutletLiquid.z(i) = INITIALIZATION.TopComposition(i)/sum(INITIALIZATION.TopComposition) +(INITIALIZATION.BottomComposition(i)/sum(INITIALIZATION.BottomComposition)-INITIALIZATION.TopComposition(i)/sum(INITIALIZATION.TopComposition) )*((j-1)/(NumberOfTrays-1)); |
---|
[721] | 2267 | end |
---|
| 2268 | |
---|
| 2269 | end |
---|
| 2270 | |
---|
[715] | 2271 | EQUATIONS |
---|
| 2272 | |
---|
| 2273 | for i in [1:NumberOfTrays] do |
---|
| 2274 | |
---|
| 2275 | "Murphree Efficiency" |
---|
[735] | 2276 | TRAYS(i).OutletVapour.z = MurphreeEff * (TRAYS(i).yideal - TRAYS(i).InletVapour.z) + TRAYS(i).InletVapour.z; |
---|
[715] | 2277 | |
---|
| 2278 | "Energy Holdup" |
---|
[735] | 2279 | TRAYS(i).E = TRAYS(i).ML*TRAYS(i).OutletLiquid.h + TRAYS(i).MV*TRAYS(i).OutletVapour.h - TRAYS(i).OutletLiquid.P*V; |
---|
[715] | 2280 | |
---|
| 2281 | "Component Molar Balance" |
---|
[735] | 2282 | diff(TRAYS(i).M)=TRAYS(i).Inlet.F*TRAYS(i).Inlet.z + TRAYS(i).InletLiquid.F*TRAYS(i).InletLiquid.z + TRAYS(i).InletVapour.F*TRAYS(i).InletVapour.z- TRAYS(i).OutletLiquid.F*TRAYS(i).OutletLiquid.z - TRAYS(i).OutletVapour.F*TRAYS(i).OutletVapour.z- |
---|
[715] | 2283 | TRAYS(i).LiquidSideStream.F*TRAYS(i).LiquidSideStream.z-TRAYS(i).VapourSideStream.F*TRAYS(i).VapourSideStream.z + stoic*TRAYS(i).r3*TRAYS(i).ML*TRAYS(i).vL; |
---|
| 2284 | |
---|
| 2285 | "Energy Balance" |
---|
[735] | 2286 | diff(TRAYS(i).E) = ( TRAYS(i).Inlet.F*TRAYS(i).Inlet.h + TRAYS(i).InletLiquid.F*TRAYS(i).InletLiquid.h + TRAYS(i).InletVapour.F*TRAYS(i).InletVapour.h- TRAYS(i).OutletLiquid.F*TRAYS(i).OutletLiquid.h - TRAYS(i).OutletVapour.F*TRAYS(i).OutletVapour.h |
---|
[715] | 2287 | -TRAYS(i).VapourSideStream.F*TRAYS(i).VapourSideStream.h - TRAYS(i).LiquidSideStream.F*TRAYS(i).LiquidSideStream.h + Q ) + Hr * TRAYS(i).r3 *TRAYS(i). vL*TRAYS(i).ML; |
---|
| 2288 | |
---|
| 2289 | "Level of clear liquid over the weir" |
---|
| 2290 | TRAYS(i).Level = TRAYS(i).ML*TRAYS(i).vL/Ap; |
---|
| 2291 | |
---|
| 2292 | "Geometry Constraint" |
---|
| 2293 | V = TRAYS(i).ML* TRAYS(i).vL + TRAYS(i).MV*TRAYS(i).vV; |
---|
| 2294 | |
---|
| 2295 | end |
---|
| 2296 | |
---|
| 2297 | switch VapourFlow |
---|
| 2298 | case "on": |
---|
| 2299 | "Pressure Drop through the condenser" |
---|
[735] | 2300 | CONDENSER.InletVapour.F*TRAYS(1).vV / 'm^2' = |
---|
| 2301 | sqrt((TRAYS(1).OutletVapour.P - CONDENSER.OutletLiquid.P + 1e-8 * 'atm')/(TRAYS(1).rhoV*alfacond)); |
---|
| 2302 | when TRAYS(1).OutletVapour.P < CONDENSER.OutletLiquid.P switchto "off"; |
---|
[715] | 2303 | |
---|
| 2304 | case "off": |
---|
| 2305 | "Prato selado" |
---|
[735] | 2306 | CONDENSER.InletVapour.F = 0.0 * 'mol/s'; |
---|
| 2307 | when TRAYS(1).OutletVapour.P > CONDENSER.OutletLiquid.P + 1e-3 * 'atm' switchto "on"; |
---|
[715] | 2308 | end |
---|
| 2309 | |
---|
| 2310 | for i in [1:NumberOfTrays] do |
---|
| 2311 | |
---|
| 2312 | switch TrayLiquidFlow |
---|
| 2313 | case "on": |
---|
| 2314 | "Francis Equation" |
---|
[735] | 2315 | TRAYS(i).OutletLiquid.F*TRAYS(i).vL = 1.84*'1/s'*lw*((TRAYS(i).Level-(beta*hw)+1e-6*'m')/(beta))^2; |
---|
[715] | 2316 | when TRAYS(i).Level < (beta * hw) switchto "off"; |
---|
| 2317 | |
---|
| 2318 | case "off": |
---|
| 2319 | "Low level" |
---|
[735] | 2320 | TRAYS(i).OutletLiquid.F = 0 * 'mol/h'; |
---|
[715] | 2321 | when TRAYS(i).Level > (beta * hw) + 1e-6*'m' switchto "on"; |
---|
| 2322 | end |
---|
| 2323 | |
---|
| 2324 | switch TrayVapourFlow |
---|
| 2325 | case "on": |
---|
[735] | 2326 | TRAYS(i).InletVapour.F*TRAYS(i).vV = sqrt((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P - TRAYS(i).Level*g*TRAYS(i).rhoL + 1e-8 * 'atm')/(TRAYS(i).rhoV*alfa))*Ah; |
---|
| 2327 | when TRAYS(i).InletVapour.P < TRAYS(i).OutletVapour.P + TRAYS(i).Level*g*TRAYS(i).rhoL switchto "off"; |
---|
[715] | 2328 | |
---|
| 2329 | case "off": |
---|
[735] | 2330 | TRAYS(i).InletVapour.F = 0 * 'mol/s'; |
---|
| 2331 | when TRAYS(i).InletVapour.P > TRAYS(i).OutletVapour.P + TRAYS(i).Level*g*TRAYS(i).rhoL + 3e-2 * 'atm' switchto "on"; |
---|
[715] | 2332 | end |
---|
| 2333 | |
---|
| 2334 | end |
---|
| 2335 | |
---|
| 2336 | # Connecting Trays |
---|
| 2337 | FeedTray.F*FeedTrayIndex= TRAYS.Inlet.F; |
---|
| 2338 | FeedTray.T = TRAYS.Inlet.T; |
---|
| 2339 | FeedTray.P = TRAYS.Inlet.P; |
---|
| 2340 | FeedTray.z = TRAYS.Inlet.z; |
---|
| 2341 | FeedTray.v = TRAYS.Inlet.v; |
---|
| 2342 | FeedTray.h = TRAYS.Inlet.h; |
---|
| 2343 | |
---|
[721] | 2344 | HeatToReboiler = REBOILER.InletQ; |
---|
| 2345 | HeatToCondenser = CONDENSER.InletQ; |
---|
[715] | 2346 | |
---|
| 2347 | "Reboiler Temperature indicator" |
---|
[735] | 2348 | TRI * 'K' = REBOILER.OutletLiquid.T; |
---|
[715] | 2349 | |
---|
| 2350 | "Reboiler Level indicator" |
---|
| 2351 | LRI*REBOILER.V = REBOILER.Level*REBOILER.Across; |
---|
| 2352 | |
---|
| 2353 | "Condenser Temperature indicator" |
---|
[735] | 2354 | TCI * 'K' = CONDENSER.OutletLiquid.T; |
---|
[715] | 2355 | |
---|
| 2356 | "Condenser Level indicator" |
---|
| 2357 | LCI*CONDENSER.V = CONDENSER.Level*CONDENSER.Across; |
---|
| 2358 | |
---|
| 2359 | # Condenser Connector Equations |
---|
| 2360 | ConnectorCondenserVout.T = VapourDistillate.T; |
---|
| 2361 | ConnectorCondenserVout.P = VapourDistillate.P; |
---|
| 2362 | ConnectorCondenserVout.F = VapourDistillate.F; |
---|
| 2363 | ConnectorCondenserVout.z = VapourDistillate.z; |
---|
| 2364 | |
---|
| 2365 | # Splitter Connector Equations |
---|
| 2366 | ConnectorSplitterOut.T = LiquidDistillate.T; |
---|
| 2367 | ConnectorSplitterOut.P = LiquidDistillate.P; |
---|
| 2368 | ConnectorSplitterOut.F = LiquidDistillate.F; |
---|
| 2369 | ConnectorSplitterOut.z = LiquidDistillate.z; |
---|
| 2370 | |
---|
| 2371 | # Reboiler Connector Equations |
---|
| 2372 | ConnectorReboilerLout.T = BottomProduct.T; |
---|
| 2373 | ConnectorReboilerLout.P = BottomProduct.P; |
---|
| 2374 | ConnectorReboilerLout.F = BottomProduct.F; |
---|
| 2375 | ConnectorReboilerLout.z = BottomProduct.z; |
---|
| 2376 | |
---|
| 2377 | VapourDrawOff.F*VapSideTrayIndex= TRAYS.VapourSideStream.F; |
---|
| 2378 | VapourDrawOff.T = TRAYS(VapourSideStreamLocation).VapourSideStream.T; |
---|
| 2379 | VapourDrawOff.P = TRAYS(VapourSideStreamLocation).VapourSideStream.P; |
---|
| 2380 | VapourDrawOff.z = TRAYS(VapourSideStreamLocation).VapourSideStream.z; |
---|
| 2381 | |
---|
| 2382 | LiquidDrawOff.F*LiqSideTrayIndex= TRAYS.LiquidSideStream.F; |
---|
| 2383 | LiquidDrawOff.T = TRAYS(LiquidSideStreamLocation).LiquidSideStream.T; |
---|
| 2384 | LiquidDrawOff.P = TRAYS(LiquidSideStreamLocation).LiquidSideStream.P; |
---|
| 2385 | LiquidDrawOff.z = TRAYS(LiquidSideStreamLocation).LiquidSideStream.z; |
---|
| 2386 | |
---|
| 2387 | VapourDrawOffFlow = VapourDrawOff.F; |
---|
| 2388 | LiquidDrawOffFlow = LiquidDrawOff.F; |
---|
| 2389 | |
---|
| 2390 | end |
---|
| 2391 | |
---|
[833] | 2392 | *# |
---|
| 2393 | |
---|
[715] | 2394 | Model Packed_Section_ColumnBasic |
---|
[833] | 2395 | |
---|
| 2396 | ATTRIBUTES |
---|
[715] | 2397 | Pallete = false; |
---|
| 2398 | Icon = "icon/PackedSectionColumn"; |
---|
| 2399 | Brief = "Model of a packed column section."; |
---|
| 2400 | Info = |
---|
| 2401 | "== Model of a packed column section containing == |
---|
| 2402 | * NStages theoretical stages. |
---|
| 2403 | |
---|
| 2404 | == Specify == |
---|
| 2405 | * the feed stream of each tray (Inlet); |
---|
[735] | 2406 | * the InletLiquid stream of the top tray; |
---|
| 2407 | * the InletVapour stream of the bottom tray; |
---|
[715] | 2408 | * the total pressure drop (dP) of the section. |
---|
| 2409 | |
---|
| 2410 | == Initial Conditions == |
---|
[735] | 2411 | * the stages temperature (OutletLiquid.T); |
---|
[715] | 2412 | * the stages liquid holdup; |
---|
[735] | 2413 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray. |
---|
[715] | 2414 | "; |
---|
| 2415 | |
---|
| 2416 | PARAMETERS |
---|
| 2417 | |
---|
[833] | 2418 | outer PP as Plugin (Brief = "External Physical Properties", Type="PP"); |
---|
| 2419 | outer NComp as Integer (Brief="Number of components"); |
---|
[715] | 2420 | |
---|
[833] | 2421 | NumberOfFeeds as Integer (Brief="Number of Feed Stages",Default=3,Protected=true); |
---|
| 2422 | NumberOfStages as Integer (Brief="Number of Stages", Default=3); |
---|
| 2423 | FeedStageLocation(NumberOfFeeds) as Integer (Brief="Feed Stage Location",Default=2); |
---|
[715] | 2424 | |
---|
[721] | 2425 | PackingHeight as length (Brief="Height of packing"); |
---|
| 2426 | HeatSupply as heat_rate (Brief="Rate of heat supply"); |
---|
[715] | 2427 | ColumnDiameter as length (Brief="Column diameter"); |
---|
[721] | 2428 | VoidFraction as Real (Brief="Void fraction of packing, m^3/m^3"); |
---|
| 2429 | ResistanceCoeff as positive (Brief="Resistance coefficient on the liquid load", Default=1); |
---|
| 2430 | AreaPerPackingVol as Real (Brief="surface area per packing volume", Unit='m^2/m^3'); |
---|
[715] | 2431 | |
---|
[833] | 2432 | V as volume (Brief="Total Volume of the tray",Protected=true); |
---|
| 2433 | Pi as constant (Brief="Pi Number",Default=3.14159265, Symbol = "\pi",Hidden=true); |
---|
| 2434 | Gconst as acceleration (Brief="Gravity Acceleration",Default=9.81,Hidden=true); |
---|
| 2435 | hs as length (Brief="Height of the packing stage",Protected=true); |
---|
| 2436 | ColumnArea as area (Brief="Column Sectional Cross Area",Protected=true); |
---|
[715] | 2437 | |
---|
| 2438 | VARIABLES |
---|
[833] | 2439 | |
---|
| 2440 | INITIALIZATION as InitializeStage (Brief = "Column Model Initialization"); |
---|
[715] | 2441 | |
---|
[833] | 2442 | in LiquidInlet as stream (Brief="Liquid Inlet in the section", PosX=0.70, PosY=0); |
---|
| 2443 | out VapourOutlet as vapour_stream (Brief="Vapour Outlet in the section", PosX=0.30, PosY=0); |
---|
| 2444 | |
---|
| 2445 | in VapourInlet as stream (Brief="Vapour Inlet in the section", PosX=0.30, PosY=1); |
---|
| 2446 | out LiquidOutlet as liquid_stream (Brief="Liquid Outlet in the section", PosX=0.70, PosY=1); |
---|
| 2447 | |
---|
| 2448 | LiquidConnector as stream (Brief="Liquid connection at the middle STAGES", PosX=0.75, PosY=1,Hidden=true); |
---|
| 2449 | VapourConnector as stream (Brief="Vapour connection at the middle STAGES", PosX=0.55, PosY=0,Hidden=true); |
---|
| 2450 | |
---|
[715] | 2451 | STAGES(NumberOfStages) as packedStage; |
---|
[721] | 2452 | PressureDrop as pressure; |
---|
[715] | 2453 | |
---|
| 2454 | SET |
---|
| 2455 | |
---|
[833] | 2456 | ColumnArea = 0.25*Pi*ColumnDiameter^2; |
---|
[715] | 2457 | hs =PackingHeight/NumberOfStages; |
---|
| 2458 | V = hs * ColumnArea; |
---|
| 2459 | |
---|
| 2460 | CONNECTIONS |
---|
| 2461 | |
---|
[735] | 2462 | STAGES([2:NumberOfStages]).OutletVapour to STAGES([1:NumberOfStages-1]).InletVapour; |
---|
| 2463 | STAGES([1:NumberOfStages-1]).OutletLiquid to STAGES([2:NumberOfStages]).InletLiquid; |
---|
[833] | 2464 | LiquidConnector to STAGES(1).InletLiquid; |
---|
| 2465 | VapourConnector to STAGES(NumberOfStages).InletVapour; |
---|
[715] | 2466 | |
---|
| 2467 | EQUATIONS |
---|
| 2468 | |
---|
[833] | 2469 | "PresureDrop" |
---|
[715] | 2470 | STAGES.deltaP = PressureDrop/NumberOfStages; |
---|
| 2471 | |
---|
[833] | 2472 | LiquidConnector.F= LiquidInlet.F; |
---|
| 2473 | LiquidConnector.T = LiquidInlet.T; |
---|
| 2474 | LiquidConnector.P = LiquidInlet.P; |
---|
| 2475 | LiquidConnector.z = LiquidInlet.z; |
---|
| 2476 | LiquidConnector.v = LiquidInlet.v; |
---|
| 2477 | LiquidConnector.h = LiquidInlet.h; |
---|
| 2478 | |
---|
| 2479 | VapourConnector.F= VapourInlet.F; |
---|
| 2480 | VapourConnector.T = VapourInlet.T; |
---|
| 2481 | VapourConnector.P = VapourInlet.P; |
---|
| 2482 | VapourConnector.z = VapourInlet.z; |
---|
| 2483 | VapourConnector.v = VapourInlet.v; |
---|
| 2484 | VapourConnector.h = VapourInlet.h; |
---|
| 2485 | |
---|
| 2486 | LiquidOutlet.F= STAGES(NumberOfStages).OutletLiquid.F; |
---|
| 2487 | LiquidOutlet.T = STAGES(NumberOfStages).OutletLiquid.T; |
---|
| 2488 | LiquidOutlet.P = STAGES(NumberOfStages).OutletLiquid.P; |
---|
| 2489 | LiquidOutlet.z = STAGES(NumberOfStages).OutletLiquid.z; |
---|
| 2490 | |
---|
| 2491 | VapourOutlet.F= STAGES(1).OutletVapour.F; |
---|
| 2492 | VapourOutlet.T = STAGES(1).OutletVapour.T; |
---|
| 2493 | VapourOutlet.P = STAGES(1).OutletVapour.P; |
---|
| 2494 | VapourOutlet.z = STAGES(1).OutletVapour.z; |
---|
| 2495 | |
---|
[715] | 2496 | for i in [1:NumberOfStages] do |
---|
| 2497 | |
---|
| 2498 | "Energy Balance" |
---|
[735] | 2499 | diff(STAGES(i).E) = ( STAGES(i).Inlet.F*STAGES(i).Inlet.h + STAGES(i).InletLiquid.F*STAGES(i).InletLiquid.h + STAGES(i).InletVapour.F*STAGES(i).InletVapour.h- STAGES(i).OutletLiquid.F*STAGES(i).OutletLiquid.h |
---|
| 2500 | - STAGES(i).OutletVapour.F*STAGES(i).OutletVapour.h + HeatSupply ); |
---|
[715] | 2501 | |
---|
| 2502 | "Energy Holdup" |
---|
[735] | 2503 | STAGES(i).E = STAGES(i).ML*STAGES(i).OutletLiquid.h + STAGES(i).MV*STAGES(i).OutletVapour.h - STAGES(i).OutletLiquid.P*V; |
---|
[715] | 2504 | |
---|
| 2505 | "Geometry Constraint" |
---|
| 2506 | V*VoidFraction= STAGES(i).ML*STAGES(i).vL + STAGES(i).MV*STAGES(i).vV; |
---|
| 2507 | |
---|
| 2508 | "Volume flow rate of vapor, m^3/m^2/s" |
---|
[833] | 2509 | #STAGES(i).uV * ((ColumnArea)*VoidFraction - STAGES(i).Al) = STAGES(i).OutletVapour.F * STAGES(i).vV; |
---|
| 2510 | STAGES(i).uV * ((ColumnArea)*VoidFraction - STAGES(i).Al) = STAGES(i).InletVapour.F * STAGES(i).vV; |
---|
[715] | 2511 | |
---|
| 2512 | "Liquid holdup" |
---|
| 2513 | STAGES(i).hl = STAGES(i).ML*STAGES(i).vL/V/VoidFraction; |
---|
| 2514 | |
---|
| 2515 | "Particle diameter" |
---|
| 2516 | STAGES(i).dp = 6 * (1-VoidFraction)/AreaPerPackingVol; |
---|
| 2517 | |
---|
| 2518 | "Wall Factor" |
---|
| 2519 | STAGES(i).invK = (1 + (2*STAGES(i).dp/(3*ColumnDiameter*(1-VoidFraction)))); |
---|
| 2520 | |
---|
| 2521 | "Reynolds number of the vapor stream" |
---|
| 2522 | STAGES(i).Rev*STAGES(i).invK = STAGES(i).dp*STAGES(i).uV*STAGES(i).rhoV / (STAGES(i).miV*(1-VoidFraction)); |
---|
| 2523 | |
---|
| 2524 | "Pressure drop and Vapor flow" |
---|
| 2525 | STAGES(i).deltaP/hs = ResistanceCoeff*AreaPerPackingVol*STAGES(i).uV^2*STAGES(i).rhoV*STAGES(i).invK / (2*(VoidFraction-STAGES(i).hl)^3); |
---|
| 2526 | |
---|
| 2527 | "Liquid holdup" |
---|
[833] | 2528 | STAGES(i).hl = (12*STAGES(i).miL*AreaPerPackingVol^2*STAGES(i).uL/STAGES(i).rhoL/Gconst)^1/3; |
---|
[715] | 2529 | |
---|
| 2530 | "Area occupied by the liquid" |
---|
| 2531 | STAGES(i).Al = STAGES(i).ML*STAGES(i).vL/hs; |
---|
| 2532 | |
---|
| 2533 | end |
---|
| 2534 | |
---|
| 2535 | INITIAL |
---|
| 2536 | |
---|
| 2537 | for i in 1:NumberOfStages do |
---|
| 2538 | |
---|
| 2539 | "The initial temperature of the STAGES" |
---|
[735] | 2540 | STAGES(i).OutletLiquid.T = INITIALIZATION.TopStageTemperature+(INITIALIZATION.BottomStageTemperature-INITIALIZATION.TopStageTemperature)*((i-1)/(NumberOfStages-1)); |
---|
[715] | 2541 | |
---|
| 2542 | "The Liquid Holdup of the STAGES" |
---|
[721] | 2543 | STAGES(i).ML = INITIALIZATION.LiquidMolarHoldup; |
---|
[715] | 2544 | |
---|
| 2545 | end |
---|
| 2546 | |
---|
| 2547 | for i in 1:NComp-1 do |
---|
| 2548 | |
---|
| 2549 | for j in 1:NumberOfStages do |
---|
| 2550 | |
---|
| 2551 | "The initial composition of the TRAYS" |
---|
[735] | 2552 | STAGES(j).OutletLiquid.z(i) = INITIALIZATION.TopStageComposition(i)/sum(INITIALIZATION.TopStageComposition) +(INITIALIZATION.BottomStageComposition(i)/sum(INITIALIZATION.BottomStageComposition)-INITIALIZATION.TopStageComposition(i)/sum(INITIALIZATION.TopStageComposition) )*((j-1)/(NumberOfStages-1)); |
---|
[715] | 2553 | |
---|
| 2554 | end |
---|
| 2555 | |
---|
| 2556 | end |
---|
| 2557 | |
---|
| 2558 | end |
---|
| 2559 | |
---|
| 2560 | Model Packed_Section_Column as Packed_Section_ColumnBasic |
---|
| 2561 | |
---|
| 2562 | ATTRIBUTES |
---|
| 2563 | Pallete = true; |
---|
| 2564 | Icon = "icon/PackedSectionColumn"; |
---|
| 2565 | Brief = "Model of a packed column section."; |
---|
| 2566 | Info = |
---|
| 2567 | "== Model of a packed column section containing == |
---|
| 2568 | * NStages theoretical stages. |
---|
| 2569 | |
---|
| 2570 | == Specify == |
---|
| 2571 | * the feed stream of each tray (Inlet); |
---|
[735] | 2572 | * the InletLiquid stream of the top tray; |
---|
| 2573 | * the InletVapour stream of the bottom tray; |
---|
[715] | 2574 | * the total pressure drop (dP) of the section. |
---|
| 2575 | |
---|
| 2576 | == Initial Conditions == |
---|
[735] | 2577 | * the stages temperature (OutletLiquid.T); |
---|
[715] | 2578 | * the stages liquid holdup; |
---|
[735] | 2579 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray. |
---|
[715] | 2580 | "; |
---|
| 2581 | |
---|
[833] | 2582 | SET |
---|
| 2583 | NumberOfFeeds = 1; |
---|
| 2584 | |
---|
[715] | 2585 | VARIABLES |
---|
| 2586 | |
---|
[833] | 2587 | in FeedStage as stream (Brief="Feed stream", PosX=0, PosY=0.55); |
---|
[715] | 2588 | |
---|
| 2589 | EQUATIONS |
---|
[833] | 2590 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2591 | # Equating Feed Tray Variables to Trays Variables |
---|
| 2592 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2593 | for i in 1:NumberOfStages do |
---|
[715] | 2594 | |
---|
[833] | 2595 | if i equal FeedStageLocation(1) then |
---|
| 2596 | |
---|
| 2597 | "FeedTrayTop Inlet Flow" |
---|
| 2598 | FeedStage.F= STAGES(i).Inlet.F; |
---|
| 2599 | |
---|
| 2600 | "FeedTrayTop Inlet Temperature" |
---|
| 2601 | FeedStage.T = STAGES(i).Inlet.T; |
---|
| 2602 | |
---|
| 2603 | "FeedTrayTop Inlet Pressure" |
---|
| 2604 | FeedStage.P = STAGES(i).Inlet.P; |
---|
| 2605 | |
---|
| 2606 | "FeedTrayTop Inlet Composition" |
---|
| 2607 | FeedStage.z = STAGES(i).Inlet.z; |
---|
| 2608 | |
---|
| 2609 | "FeedTrayTop Inlet Vapour Fraction" |
---|
| 2610 | FeedStage.v = STAGES(i).Inlet.v; |
---|
| 2611 | |
---|
| 2612 | "FeedTrayTop Inlet Enthalpy" |
---|
| 2613 | FeedStage.h = STAGES(i).Inlet.h; |
---|
| 2614 | |
---|
| 2615 | else |
---|
| 2616 | |
---|
| 2617 | "Inlet Tray - Flow Sealed" |
---|
| 2618 | 0*'mol/h'= STAGES(i).Inlet.F; |
---|
[715] | 2619 | |
---|
[833] | 2620 | "Inlet Tray - Temperature" |
---|
| 2621 | 300*'K' = STAGES(i).Inlet.T; |
---|
[715] | 2622 | |
---|
[833] | 2623 | "Inlet Tray - Pressure" |
---|
| 2624 | 1*'atm' = STAGES(i).Inlet.P; |
---|
[715] | 2625 | |
---|
[833] | 2626 | "Inlet Tray - Composition" |
---|
| 2627 | 0.1 = STAGES(i).Inlet.z; |
---|
| 2628 | |
---|
| 2629 | "Inlet Tray - Vapour Fraction" |
---|
| 2630 | 0 = STAGES(i).Inlet.v; |
---|
| 2631 | |
---|
| 2632 | "Inlet Tray - Enthalpy" |
---|
| 2633 | 0*'J/mol' = STAGES(i).Inlet.h; |
---|
[715] | 2634 | |
---|
| 2635 | end |
---|
| 2636 | |
---|
[833] | 2637 | end |
---|
| 2638 | |
---|
| 2639 | end |
---|
| 2640 | |
---|
[715] | 2641 | Model PackedDistillation_kettle_cond as Packed_Section_ColumnBasic |
---|
| 2642 | ATTRIBUTES |
---|
| 2643 | Pallete = true; |
---|
| 2644 | Icon = "icon/PackedDistillationKettleCond"; |
---|
| 2645 | Brief = "Model of a distillation column with dynamic condenser and dynamic reboiler."; |
---|
| 2646 | Info = |
---|
| 2647 | "== Specify == |
---|
| 2648 | * the feed stream of each tray (Inlet); |
---|
| 2649 | * the pump pressure difference; |
---|
| 2650 | * the total pressure drop (dP) of the packing; |
---|
| 2651 | * the heat supllied in reboiler and condenser; |
---|
[735] | 2652 | * the condenser vapor outlet flow (OutletVapour.F); |
---|
| 2653 | * the reboiler liquid outlet flow (OutletLiquid.F); |
---|
[715] | 2654 | * both splitter outlet flows OR one of the splitter outlet flows and the splitter frac. |
---|
| 2655 | |
---|
| 2656 | == Initial Conditions == |
---|
[735] | 2657 | * the stages temperature (OutletLiquid.T); |
---|
[715] | 2658 | * the stages initial molar holdup; |
---|
[735] | 2659 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each stage; |
---|
[715] | 2660 | |
---|
[735] | 2661 | * the condenser temperature (OutletLiquid.T); |
---|
[715] | 2662 | * the condenser liquid level (Level); |
---|
[735] | 2663 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions; |
---|
[715] | 2664 | |
---|
[735] | 2665 | * the reboiler temperature (OutletLiquid.T); |
---|
[715] | 2666 | * the reboiler liquid level (Level); |
---|
[735] | 2667 | * (NoComps - 1) OutletLiquid (OR OutletVapour) compositions. |
---|
[715] | 2668 | "; |
---|
| 2669 | |
---|
| 2670 | PARAMETERS |
---|
| 2671 | |
---|
[721] | 2672 | VapourFlow as Switcher (Valid = ["on", "off"], Hidden = true ,Default = "on"); |
---|
[715] | 2673 | |
---|
| 2674 | VARIABLES |
---|
[721] | 2675 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2676 | # Column Devices |
---|
| 2677 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2678 | CONDENSER as condenser (Brief="dymamic condenser"); |
---|
| 2679 | REBOILER as reboiler (Brief="kettle reboiler"); |
---|
[739] | 2680 | SPLITTER as splitter2 (Brief="splitter which separate reflux and distillate"); |
---|
[721] | 2681 | PUMP as pump (Brief="pump in reflux stream"); |
---|
| 2682 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2683 | # Heat Ports |
---|
| 2684 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2685 | in HeatToReboiler as power (Brief="Heat supplied to Reboiler",Protected = true, PosX=1, PosY=0.885); |
---|
| 2686 | in HeatToCondenser as power (Brief="Heat supplied to Condenser", Protected = true, PosX=1, PosY=0.034); |
---|
| 2687 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2688 | # Column Products |
---|
| 2689 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2690 | out VapourDistillate as vapour_stream (Brief="Vapour outlet stream From Top Condenser", PosX=0.66, PosY=0); |
---|
| 2691 | out LiquidDistillate as liquid_stream (Brief="Liquid outlet stream From Top Splitter", PosX=1, PosY=0.20); |
---|
| 2692 | out BottomProduct as liquid_stream (Brief="Liquid outlet stream From Reboiler", PosX=0.68, PosY=1); |
---|
| 2693 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2694 | # Column Connectors |
---|
| 2695 | #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ |
---|
| 2696 | in ConnectorCondenserVout as stream (Brief="Connector for Vapour outlet stream From Top Condenser", Hidden=true); |
---|
| 2697 | in ConnectorSplitterOut as stream (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true); |
---|
| 2698 | in ConnectorReboilerLout as stream (Brief="Connector for Liquid outlet stream From Reboiler", Hidden=true); |
---|
| 2699 | ConnectorHeatReboiler as power (Brief="Heat supplied to Reboiler", Hidden=true); |
---|
| 2700 | ConnectorHeatCondenser as power (Brief="Heat supplied to Condenser", Hidden=true); |
---|
[715] | 2701 | |
---|
| 2702 | CONNECTIONS |
---|
| 2703 | |
---|
| 2704 | #vapor |
---|
[735] | 2705 | REBOILER.OutletVapour to STAGES(NumberOfStages).InletVapour; |
---|
| 2706 | STAGES(1).OutletVapour to CONDENSER.InletVapour; |
---|
[715] | 2707 | |
---|
| 2708 | #liquid |
---|
[735] | 2709 | CONDENSER.OutletLiquid to SPLITTER.Inlet; |
---|
[721] | 2710 | SPLITTER.Outlet2 to PUMP.Inlet; |
---|
[735] | 2711 | PUMP.Outlet to STAGES(1).InletLiquid; |
---|
| 2712 | STAGES(NumberOfStages).OutletLiquid to REBOILER.InletLiquid; |
---|
[715] | 2713 | |
---|
| 2714 | #Connectors |
---|
[721] | 2715 | ConnectorHeatReboiler to REBOILER.InletQ; |
---|
| 2716 | ConnectorHeatCondenser to CONDENSER.InletQ; |
---|
[715] | 2717 | |
---|
[735] | 2718 | CONDENSER.OutletVapour to ConnectorCondenserVout; |
---|
[721] | 2719 | SPLITTER.Outlet1 to ConnectorSplitterOut; |
---|
[735] | 2720 | REBOILER.OutletLiquid to ConnectorReboilerLout; |
---|
[715] | 2721 | |
---|
| 2722 | EQUATIONS |
---|
[721] | 2723 | HeatToReboiler = ConnectorHeatReboiler; |
---|
| 2724 | HeatToCondenser = ConnectorHeatCondenser; |
---|
| 2725 | |
---|
[715] | 2726 | # Condenser Connector Equations |
---|
| 2727 | ConnectorCondenserVout.T = VapourDistillate.T; |
---|
| 2728 | ConnectorCondenserVout.P = VapourDistillate.P; |
---|
| 2729 | ConnectorCondenserVout.F = VapourDistillate.F; |
---|
| 2730 | ConnectorCondenserVout.z = VapourDistillate.z; |
---|
| 2731 | |
---|
| 2732 | # Splitter Connector Equations |
---|
| 2733 | ConnectorSplitterOut.T = LiquidDistillate.T; |
---|
| 2734 | ConnectorSplitterOut.P = LiquidDistillate.P; |
---|
| 2735 | ConnectorSplitterOut.F = LiquidDistillate.F; |
---|
| 2736 | ConnectorSplitterOut.z = LiquidDistillate.z; |
---|
| 2737 | |
---|
| 2738 | # Reboiler Connector Equations |
---|
| 2739 | ConnectorReboilerLout.T = BottomProduct.T; |
---|
| 2740 | ConnectorReboilerLout.P = BottomProduct.P; |
---|
| 2741 | ConnectorReboilerLout.F = BottomProduct.F; |
---|
| 2742 | ConnectorReboilerLout.z = BottomProduct.z; |
---|
| 2743 | |
---|
| 2744 | end |
---|
| 2745 | |
---|