- Timestamp:
- Sep 1, 2008, 7:45:47 PM (15 years ago)
- Location:
- branches/gui
- Files:
-
- 4 edited
-
eml/stage_separators/columnTeste.mso (modified) (25 diffs)
-
eml/stage_separators/tray.mso (modified) (1 diff)
-
eml/stage_separators/trayTeste.mso (modified) (2 diffs)
-
sample/stage_separators/sample_columnTeste.mso (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/gui/eml/stage_separators/columnTeste.mso
r612 r613 66 66 67 67 PARAMETERS 68 outer PP as Plugin (Brief = "External Physical Properties", Type="PP"); 69 outer NComp as Integer; 70 NTrays as Integer (Brief="Number of trays", Default=2); 71 FeedTray as Integer (Brief="Number of Feed tray", Default=2); 72 topdown as Integer (Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1,Hidden=true); 73 top as Integer (Brief="Number of top tray",Hidden=true); 74 bot as Integer (Brief="Number of bottom tray",Hidden=true); 75 g as acceleration (Default=9.81,Hidden=true); 68 outer PP as Plugin (Brief = "External Physical Properties", Type="PP"); 69 outer NComp as Integer; 70 NTrays as Integer (Brief="Number of trays", Default=2); 71 SpecialTrayIndex(NTrays) as Integer (Brief="Number of trays", Default=0,Hidden=true); 72 FeedTray as Integer (Brief="Number of Feed tray", Default=2); 73 topdown as Integer (Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1,Hidden=true); 74 top as Integer (Brief="Number of top tray",Hidden=true); 75 bot as Integer (Brief="Number of bottom tray",Hidden=true); 76 g as acceleration (Default=9.81,Hidden=true); 77 Mw(NComp) as molweight (Brief="Component Mol Weight"); 78 79 80 Ah as area (Brief="Total holes area"); 81 lw as length (Brief="Weir length"); 82 hw as length (Brief="Weir height"); 83 beta as fraction (Brief="Aeration fraction"); 84 alfa as fraction (Brief="Dry pressure drop coefficient"); 85 w as Real (Brief="Feeherys correlation coefficient", Unit='1/m^4', Default=1); 86 btray as Real (Brief="Elgues correlation coefficient", Unit='kg/m/mol^2', Default=1); 87 fw as Real (Brief="Olsens correlation coefficient", Default=1); 88 Np as Real (Brief="Number of liquid passes in the tray", Default=1); 76 89 77 90 VapourFlow as Switcher (Valid = ["on", "off"], Default = "on",Hidden=true); … … 83 96 top = (NTrays-1)*(1-topdown)/2+1; 84 97 bot = NTrays/top; 98 SpecialTrayIndex(FeedTray) =1; 99 Mw = PP.MolecularWeight(); 85 100 86 101 VARIABLES 87 zero as stream (Brief="Dummy stream",Hidden=true);88 102 in Inlet as stream (Brief="Feed stream", PosX=0, PosY=0.50); 89 103 trays(NTrays) as trayTeste; … … 94 108 out Lout as liquid_stream (Brief="Feed stream", PosX=0.75, PosY=1); 95 109 out Vout as vapour_stream (Brief="Feed stream", PosX=0.55, PosY=0); 110 111 LiquidConnector as stream (Brief="Feed stream", PosX=0.75, PosY=1,Hidden=true); 112 VapourConnector as stream (Brief="Feed stream", PosX=0.55, PosY=0,Hidden=true); 96 113 97 114 CONNECTIONS 115 98 116 trays([top+topdown:topdown:bot]).OutletV to trays([top:topdown:bot-topdown]).InletV; 99 117 trays([top:topdown:bot-topdown]).OutletL to trays([top+topdown:topdown:bot]).InletL; 100 118 119 LiquidConnector to trays(1).InletL; 120 VapourConnector to trays(NTrays).InletV; 121 101 122 EQUATIONS 102 123 103 # Zero stream definition 104 zero.F = 0 * 'kmol/h'; 105 zero.T = 300 * 'K'; 106 zero.P = 1 * 'atm'; 107 zero.z = 1/NComp; 108 zero.v = 0; 109 zero.h = 0 * 'J/mol'; 110 111 # Connecting Feed Tray 112 Inlet.F = trays(FeedTray).Inlet.F; 113 Inlet.T = trays(FeedTray).Inlet.T; 114 Inlet.P = trays(FeedTray).Inlet.P; 115 Inlet.z = trays(FeedTray).Inlet.z; 116 Inlet.v = trays(FeedTray).Inlet.v; 117 Inlet.h = trays(FeedTray).Inlet.h; 118 119 if FeedTray equal 1 120 121 then 122 123 for i in [2:NTrays] do 124 125 zero.F = trays(i).Inlet.F; 126 zero.T = trays(i).Inlet.T; 127 zero.P = trays(i).Inlet.P; 128 zero.z = trays(i).Inlet.z; 129 zero.v = trays(i).Inlet.v; 130 zero.h = trays(i).Inlet.h; 131 132 end 133 134 else 135 136 if FeedTray equal NTrays 137 138 then 139 140 for i in [1:NTrays-1] do 141 142 zero.F = trays(i).Inlet.F; 143 zero.T = trays(i).Inlet.T; 144 zero.P = trays(i).Inlet.P; 145 zero.z = trays(i).Inlet.z; 146 zero.v = trays(i).Inlet.v; 147 zero.h = trays(i).Inlet.h; 148 149 end 150 151 else 152 153 for i in [2:FeedTray-1] do 154 155 zero.F = trays(i).Inlet.F; 156 zero.T = trays(i).Inlet.T; 157 zero.P = trays(i).Inlet.P; 158 zero.z = trays(i).Inlet.z; 159 zero.v = trays(i).Inlet.v; 160 zero.h = trays(i).Inlet.h; 161 162 end 163 164 for i in [FeedTray+1: NTrays-1] do 165 166 zero.F = trays(i).Inlet.F; 167 zero.T = trays(i).Inlet.T; 168 zero.P = trays(i).Inlet.P; 169 zero.z = trays(i).Inlet.z; 170 zero.v = trays(i).Inlet.v; 171 zero.h = trays(i).Inlet.h; 172 173 end 174 175 end 176 177 end 124 LiquidConnector.F= Lin.F; 125 LiquidConnector.T = Lin.T; 126 LiquidConnector.P = Lin.P; 127 LiquidConnector.z = Lin.z; 128 LiquidConnector.v = Lin.v; 129 LiquidConnector.h = Lin.h; 130 131 VapourConnector.F= Vin.F; 132 VapourConnector.T = Vin.T; 133 VapourConnector.P = Vin.P; 134 VapourConnector.z = Vin.z; 135 VapourConnector.v = Vin.v; 136 VapourConnector.h = Vin.h; 137 138 Lout.F= trays(NTrays).OutletL.F; 139 Lout.T = trays(NTrays).OutletL.T; 140 Lout.P = trays(NTrays).OutletL.P; 141 Lout.z = trays(NTrays).OutletL.z; 142 #Lout.v = trays(NTrays).OutletL.v; 143 # Lout.h = trays(NTrays).OutletL.h; 144 145 Vout.F= trays(1).OutletV.F; 146 Vout.T = trays(1).OutletV.T; 147 Vout.P = trays(1).OutletV.P; 148 Vout.z = trays(1).OutletV.z; 149 #Vout.v = trays(1).OutletV.v; 150 #Vout.h = trays(1).OutletV.h; 151 152 # Connecting Trays 153 Inlet.F*SpecialTrayIndex= trays.Inlet.F; 154 Inlet.T = trays.Inlet.T; 155 Inlet.P = trays.Inlet.P; 156 Inlet.z = trays.Inlet.z; 157 Inlet.v = trays.Inlet.v; 158 Inlet.h = trays.Inlet.h; 178 159 179 160 for i in [1:NTrays] do … … 184 165 case "default": 185 166 "Francis Equation" 186 trays(i).OutletL.F*trays(i).vL = 1.84*'1/s'* trays(i).lw*((trays(i).Level-(trays(i).beta*trays(i).hw))/(trays(i).beta))^2;167 trays(i).OutletL.F*trays(i).vL = 1.84*'1/s'*lw*((trays(i).Level-(beta*hw))/(beta))^2; 187 168 188 169 case "Wang_Fl": 189 trays(i).OutletL.F*trays(i).vL = 1.84*'m^0.5/s'* trays(i).lw*((trays(i).Level-(trays(i).beta*trays(i).hw))/(trays(i).beta))^1.5;170 trays(i).OutletL.F*trays(i).vL = 1.84*'m^0.5/s'*lw*((trays(i).Level-(beta*hw))/(beta))^1.5; 190 171 191 172 case "Olsen": 192 trays(i).OutletL.F / 'mol/s'= trays(i).lw*trays(i).Np*trays(i).rhoL/sum(trays(i).Mw*trays(i).OutletV.z)/(0.665*trays(i).fw)^1.5 * ((trays(i).ML*sum(trays(i).Mw*trays(i).OutletL.z)/trays(i).rhoL/trays(i).Ap)-trays(i).hw)^1.5 * 'm^0.5/mol';173 trays(i).OutletL.F / 'mol/s'= lw*Np*trays(i).rhoL/sum(Mw*trays(i).OutletV.z)/(0.665*fw)^1.5 * ((trays(i).ML*sum(Mw*trays(i).OutletL.z)/trays(i).rhoL/trays(i).Ap)-hw)^1.5 * 'm^0.5/mol'; 193 174 194 175 case "Feehery_Fl": 195 trays(i).OutletL.F = trays(i).lw*trays(i).rhoL/sum(trays(i).Mw*trays(i).OutletL.z) * ((trays(i).Level-trays(i).hw)/750/'mm')^1.5 * 'm^2/s';176 trays(i).OutletL.F = lw*trays(i).rhoL/sum(Mw*trays(i).OutletL.z) * ((trays(i).Level-hw)/750/'mm')^1.5 * 'm^2/s'; 196 177 197 178 case "Roffel_Fl": 198 trays(i).OutletL.F = 2/3*sqrt(2*g)*trays(i).rhoL/sum( trays(i).Mw*trays(i).OutletL.z)*trays(i).lw*(2*trays(i).btemp-1)*(trays(i).ML*sum(trays(i).Mw*trays(i).OutletL.z)/(trays(i).Ap*1.3)/trays(i).rhoL/(2*trays(i).btemp-1))^1.5;179 trays(i).OutletL.F = 2/3*sqrt(2*g)*trays(i).rhoL/sum(Mw*trays(i).OutletL.z)*lw*(2*trays(i).btemp-1)*(trays(i).ML*sum(Mw*trays(i).OutletL.z)/(trays(i).Ap*1.3)/trays(i).rhoL/(2*trays(i).btemp-1))^1.5; 199 180 end 200 when trays(i).Level < ( trays(i).beta * trays(i).hw) switchto "off";181 when trays(i).Level < (beta *hw) switchto "off"; 201 182 202 183 case "off": 203 184 "Low level" 204 185 trays(i).OutletL.F = 0 * 'mol/h'; 205 when trays(i).Level > ( trays(i).beta * trays(i).hw) + 1e-6*'m' switchto "on";186 when trays(i).Level > (beta * hw) + 1e-6*'m' switchto "on"; 206 187 end 207 188 208 trays(i).btemp = 1 - 0.3593/'Pa^0.0888545'*abs(trays(i).OutletV.F*sum( trays(i).Mw*trays(i).OutletV.z)/(trays(i).Ap*1.3)/sqrt(trays(i).rhoV))^0.177709; #/'(kg/m)^0.0888545/s^0.177709';189 trays(i).btemp = 1 - 0.3593/'Pa^0.0888545'*abs(trays(i).OutletV.F*sum(Mw*trays(i).OutletV.z)/(trays(i).Ap*1.3)/sqrt(trays(i).rhoV))^0.177709; #/'(kg/m)^0.0888545/s^0.177709'; 209 190 210 191 switch VapourFlow … … 212 193 switch VapourFlowModel 213 194 case "Reepmeyer": 214 trays(i).InletV.F*trays(i).vV = sqrt((trays(i).InletV.P - trays(i).OutletV.P)/(trays(i).rhoV* trays(i).alfa))*trays(i).Ah;195 trays(i).InletV.F*trays(i).vV = sqrt((trays(i).InletV.P - trays(i).OutletV.P)/(trays(i).rhoV*alfa))*Ah; 215 196 216 197 case "Feehery_Fv": 217 trays(i).InletV.F = trays(i).rhoV/trays(i).Ap/ trays(i).w/sum(trays(i).Mw*trays(i).OutletV.z) * sqrt(((trays(i).InletV.P - trays(i).OutletV.P)-(trays(i).rhoV*g*trays(i).ML*trays(i).vL/trays(i).Ap))/trays(i).rhoV);198 trays(i).InletV.F = trays(i).rhoV/trays(i).Ap/w/sum(Mw*trays(i).OutletV.z) * sqrt(((trays(i).InletV.P - trays(i).OutletV.P)-(trays(i).rhoV*g*trays(i).ML*trays(i).vL/trays(i).Ap))/trays(i).rhoV); 218 199 219 200 case "Roffel_Fv": 220 trays(i).InletV.F^1.08 * 0.0013 * 'kg/m/mol^1.08/s^0.92*1e5' = (trays(i).InletV.P - trays(i).OutletV.P)*1e5 - ( trays(i).beta*sum(trays(i).M*trays(i).Mw)/(trays(i).Ap*1.3)*g*1e5) * (trays(i).rhoV*trays(i).Ah/sum(trays(i).Mw*trays(i).OutletV.z))^1.08 * 'm^1.08/mol^1.08';201 trays(i).InletV.F^1.08 * 0.0013 * 'kg/m/mol^1.08/s^0.92*1e5' = (trays(i).InletV.P - trays(i).OutletV.P)*1e5 - (beta*sum(trays(i).M*Mw)/(trays(i).Ap*1.3)*g*1e5) * (trays(i).rhoV*Ah/sum(Mw*trays(i).OutletV.z))^1.08 * 'm^1.08/mol^1.08'; 221 202 222 203 case "Klingberg": … … 224 205 225 206 case "Wang_Fv": 226 trays(i).InletV.F * trays(i).vV = trays(i).Ap * sqrt(((trays(i).InletV.P - trays(i).OutletV.P)-trays(i).rhoL*g*trays(i).Level)/trays(i).rhoV* trays(i).alfa);207 trays(i).InletV.F * trays(i).vV = trays(i).Ap * sqrt(((trays(i).InletV.P - trays(i).OutletV.P)-trays(i).rhoL*g*trays(i).Level)/trays(i).rhoV*alfa); 227 208 228 209 case "Elgue": 229 trays(i).InletV.F = sqrt((trays(i).InletV.P - trays(i).OutletV.P)/ trays(i).btray);210 trays(i).InletV.F = sqrt((trays(i).InletV.P - trays(i).OutletV.P)/btray); 230 211 end 231 212 when trays(i).InletV.F < 1e-6 * 'kmol/h' switchto "off"; … … 240 221 end 241 222 223 242 224 #*---------------------------------------------------------------------- 243 225 * Model of a distillation column containing: … … 247 229 * - a splitter which separate reflux and distillate; 248 230 * - a pump in reflux stream; 249 *--------------------------------------------------------------------- *#231 *--------------------------------------------------------------------- 250 232 Model Distillation_kettle_condTeste 251 233 ATTRIBUTES … … 339 321 * - a pump in reflux stream. 340 322 * 341 * ------------------------------------------------------------------ *#323 * ------------------------------------------------------------------ 342 324 Model Distillation_thermosyphon_subcoolingTeste 343 325 ATTRIBUTES … … 435 417 * - a pump in reflux stream. 436 418 * 437 * ------------------------------------------------------------------* #419 * ------------------------------------------------------------------* 438 420 Model Distillation_thermosyphon_condTeste 439 421 ATTRIBUTES … … 526 508 * - a pump in reflux stream. 527 509 * 528 * ------------------------------------------------------------------* #510 * ------------------------------------------------------------------* 529 511 Model Distillation_kettle_subcoolingTeste 530 512 ATTRIBUTES … … 613 595 * - a splitter which separate reflux and distillate; 614 596 * - a pump in reflux stream; 615 *---------------------------------------------------------------------* #597 *---------------------------------------------------------------------* 616 598 Model RectifierTeste 617 599 ATTRIBUTES … … 693 675 * - a pump in reflux stream. 694 676 * 695 * ------------------------------------------------------------------* #677 * ------------------------------------------------------------------* 696 678 Model Rectifier_subcoolingTeste 697 679 ATTRIBUTES … … 773 755 * - a splitter which separate reflux and distillate; 774 756 * - a pump in reflux stream; 775 *---------------------------------------------------------------------* #757 *---------------------------------------------------------------------* 776 758 Model Refluxed_StrippingTeste 777 759 ATTRIBUTES … … 853 835 * - a pump in reflux stream. 854 836 * 855 * ------------------------------------------------------------------* #837 * ------------------------------------------------------------------* 856 838 Model Refluxed_Stripping_subcoolingTeste 857 839 ATTRIBUTES … … 933 915 * - a splitter which separate reflux and distillate; 934 916 * - a pump in reflux stream; 935 *---------------------------------------------------------------------* #917 *---------------------------------------------------------------------* 936 918 Model Refluxed_AbsorptionTeste 937 919 ATTRIBUTES … … 1014 996 * - a pump in reflux stream. 1015 997 * 1016 * ------------------------------------------------------------------* #998 * ------------------------------------------------------------------* 1017 999 Model Refluxed_Absorption_subcoolingTeste 1018 1000 ATTRIBUTES … … 1094 1076 * - a kettle reboiler; 1095 1077 * 1096 * ------------------------------------------------------------------* #1078 * ------------------------------------------------------------------* 1097 1079 Model Reboiled_Stripping_kettleTeste 1098 1080 ATTRIBUTES … … 1154 1136 * - steady state reboiler (thermosyphon); 1155 1137 * 1156 * ------------------------------------------------------------------* #1138 * ------------------------------------------------------------------* 1157 1139 Model Reboiled_Stripping_thermosyphonTeste 1158 1140 ATTRIBUTES … … 1217 1199 * - a kettle reboiler; 1218 1200 * 1219 * ------------------------------------------------------------------* #1201 * ------------------------------------------------------------------* 1220 1202 Model Reboiled_Absorption_kettleTeste 1221 1203 ATTRIBUTES … … 1277 1259 * - steady state reboiler (thermosyphon); 1278 1260 * 1279 * ------------------------------------------------------------------* #1261 * ------------------------------------------------------------------* 1280 1262 Model Reboiled_Absorption_thermosyphonTeste 1281 1263 ATTRIBUTES … … 1336 1318 * Reactive Distillation Column 1337 1319 * 1338 * ------------------------------------------------------------------* #1320 * ------------------------------------------------------------------* 1339 1321 Model ReactiveDistillationTeste 1340 1322 ATTRIBUTES … … 1423 1405 * - NStages = theoretical number of equilibrium stages. 1424 1406 * 1425 *---------------------------------------------------------------------* #1407 *---------------------------------------------------------------------* 1426 1408 Model Packed_Section_ColumnTeste 1427 1409 ATTRIBUTES … … 1479 1461 * - a splitter which separate reflux and distillate; 1480 1462 * - a pump in reflux stream; 1481 *---------------------------------------------------------------------* #1463 *---------------------------------------------------------------------* 1482 1464 Model PackedDistillation_kettle_cond as Packed_Section_ColumnTeste 1483 1465 ATTRIBUTES … … 1542 1524 1543 1525 end 1526 -
branches/gui/eml/stage_separators/tray.mso
r522 r613 197 197 end 198 198 199 btemp = 1 - 0.3593/'Pa^0.0888545'* (OutletV.F*sum(Mw*OutletV.z)/(Ap*1.3)/sqrt(rhoV))^0.177709; #/'(kg/m)^0.0888545/s^0.177709';199 btemp = 1 - 0.3593/'Pa^0.0888545'*abs(OutletV.F*sum(Mw*OutletV.z)/(Ap*1.3)/sqrt(rhoV))^0.177709; #/'(kg/m)^0.0888545/s^0.177709'; 200 200 201 201 switch VapourFlow -
branches/gui/eml/stage_separators/trayTeste.mso
r612 r613 137 137 "; 138 138 139 PARAMETERS140 Ah as area (Brief="Total holes area");141 lw as length (Brief="Weir length");142 hw as length (Brief="Weir height");143 beta as fraction (Brief="Aeration fraction");144 alfa as fraction (Brief="Dry pressure drop coefficient");145 w as Real (Brief="Feeherys correlation coefficient", Unit='1/m^4', Default=1);146 btray as Real (Brief="Elgues correlation coefficient", Unit='kg/m/mol^2', Default=1);147 fw as Real (Brief="Olsens correlation coefficient", Default=1);148 Np as Real (Brief="Number of liquid passes in the tray", Default=1);149 Mw(NComp) as molweight (Brief="Component Mol Weight");150 151 # VapourFlow as Switcher (Valid = ["on", "off"], Default = "on",Hidden=true);152 # LiquidFlow as Switcher (Valid = ["on", "off"], Default = "on",Hidden=true);153 # VapourFlowModel as Switcher (Valid = ["Reepmeyer", "Feehery_Fv", "Roffel_Fv", "Klingberg", "Wang_Fv", "Elgue"], Default = "Reepmeyer");154 # LiquidFlowModel as Switcher (Valid = ["default", "Wang_Fl", "Olsen", "Feehery_Fl", "Roffel_Fl"], Default = "default");155 156 SET157 Mw = PP.MolecularWeight();158 159 139 VARIABLES 160 140 rhoL as dens_mass; … … 171 151 rhoV = PP.VapourDensity(InletV.T, InletV.P, InletV.z); 172 152 173 #* switch LiquidFlow174 case "on":175 switch LiquidFlowModel176 case "default":177 "Francis Equation"178 OutletL.F*vL = 1.84*'1/s'*lw*((Level-(beta*hw))/(beta))^2;179 180 case "Wang_Fl":181 OutletL.F*vL = 1.84*'m^0.5/s'*lw*((Level-(beta*hw))/(beta))^1.5;182 183 case "Olsen":184 OutletL.F / 'mol/s'= lw*Np*rhoL/sum(Mw*OutletV.z)/(0.665*fw)^1.5 * ((ML*sum(Mw*OutletL.z)/rhoL/Ap)-hw)^1.5 * 'm^0.5/mol';185 186 case "Feehery_Fl":187 OutletL.F = lw*rhoL/sum(Mw*OutletL.z) * ((Level-hw)/750/'mm')^1.5 * 'm^2/s';188 189 case "Roffel_Fl":190 OutletL.F = 2/3*sqrt(2*g)*rhoL/sum(Mw*OutletL.z)*lw*(2*btemp-1)*(ML*sum(Mw*OutletL.z)/(Ap*1.3)/rhoL/(2*btemp-1))^1.5;191 end192 when Level < (beta * hw) switchto "off";193 194 case "off":195 "Low level"196 OutletL.F = 0 * 'mol/h';197 when Level > (beta * hw) + 1e-6*'m' switchto "on";198 end199 200 btemp = 1 - 0.3593/'Pa^0.0888545'*abs(OutletV.F*sum(Mw*OutletV.z)/(Ap*1.3)/sqrt(rhoV))^0.177709; #/'(kg/m)^0.0888545/s^0.177709';201 202 switch VapourFlow203 case "on":204 switch VapourFlowModel205 case "Reepmeyer":206 InletV.F*vV = sqrt((InletV.P - OutletV.P)/(rhoV*alfa))*Ah;207 208 case "Feehery_Fv":209 InletV.F = rhoV/Ap/w/sum(Mw*OutletV.z) * sqrt(((InletV.P - OutletV.P)-(rhoV*g*ML*vL/Ap))/rhoV);210 211 case "Roffel_Fv":212 InletV.F^1.08 * 0.0013 * 'kg/m/mol^1.08/s^0.92*1e5' = (InletV.P - OutletV.P)*1e5 - (beta*sum(M*Mw)/(Ap*1.3)*g*1e5) * (rhoV*Ah/sum(Mw*OutletV.z))^1.08 * 'm^1.08/mol^1.08';213 214 case "Klingberg":215 InletV.F * vV = Ap * sqrt(((InletV.P - OutletV.P)-rhoL*g*Level)/rhoV);216 217 case "Wang_Fv":218 InletV.F * vV = Ap * sqrt(((InletV.P - OutletV.P)-rhoL*g*Level)/rhoV*alfa);219 220 case "Elgue":221 InletV.F = sqrt((InletV.P - OutletV.P)/btray);222 end223 when InletV.F < 1e-6 * 'kmol/h' switchto "off";224 225 case "off":226 InletV.F = 0 * 'mol/s';227 when InletV.P > OutletV.P + Level*g*rhoL + 1e-1 * 'atm' switchto "on";228 end229 *#230 153 end 231 154 -
branches/gui/sample/stage_separators/sample_columnTeste.mso
r612 r613 28 28 using "stage_separators/columnTeste"; 29 29 30 # column section with 2 trays31 30 FlowSheet SectionColumn_Test_with2tray 32 31 33 32 PARAMETERS 34 33 PP as Plugin(Brief="Physical Properties", … … 41 40 42 41 DEVICES 43 sec 44 feed 45 reb 46 cond 47 42 sec as Section_ColumnTeste; 43 feed as liquid_stream; 44 reb as vapour_stream; 45 cond as liquid_stream; 46 48 47 SET 49 48 NComp = PP.NumberOfComponents; 50 sec.FeedTray = 2; 49 sec.FeedTray=1; 50 sec.NTrays = 2; 51 #COLUMN 52 sec.trays.V = 4 * 'ft^3'; 53 sec.Ah = 0.394 * 'ft^2'; 54 sec.lw = 20.94 * 'in'; 55 sec.hw = 0.125 * 'ft'; 56 sec.trays.Q = 0 * 'kW'; 57 sec.beta = 0.6; 58 sec.alfa = 4; 59 sec.trays.Ap = 3.94 * 'ft^2'; 51 60 52 61 CONNECTIONS 53 feed to sec.Inlet;#Feed at stage 2 54 62 feed to sec.Inlet; 63 reb to sec.Vin; 64 cond to sec.Lin; 65 55 66 SPECIFY 56 67 feed.F = 113.4 * 'kmol/h'; … … 71 82 sec.trays(1).OutletV.F = 150 * 'kmol/h'; 72 83 73 SET74 sec.NTrays = 2;75 #COLUMN76 sec.trays.V = 4 * 'ft^3';77 sec.trays.Ah = 0.394 * 'ft^2';78 sec.trays.lw = 20.94 * 'in';79 sec.trays.hw = 0.125 * 'ft';80 sec.trays.Q = 0 * 'kW';81 sec.trays.beta = 0.6;82 sec.trays.alfa = 4;83 sec.trays.Ap = 3.94 * 'ft^2';84 85 84 INITIAL 86 85 sec.trays.OutletL.T = 270 *'K'; 87 sec.trays.Level = 0.9 * sec. trays.hw;86 sec.trays.Level = 0.9 * sec.hw; 88 87 sec.trays.OutletL.z(1) = 0.5; 89 88 90 OPTIONS89 OPTIONS 91 90 TimeStep = 10; 92 91 TimeEnd = 1000; 93 92 94 93 end 95
Note: See TracChangeset
for help on using the changeset viewer.
