Changeset 735 for branches/gui/eml/stage_separators/tray.mso
 Timestamp:
 Feb 26, 2009, 11:00:39 PM (13 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/gui/eml/stage_separators/tray.mso
r709 r735 47 47 VapourSideStream as vapour_stream (Brief="vapour Sidestream", Hidden=true, Symbol="_{outV}"); 48 48 49 in InletL as stream (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}");50 in InletV as stream (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}");51 out OutletL as liquid_stream (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}");52 out OutletV as vapour_stream (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}");49 in InletLiquid as stream (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}"); 50 in InletVapour as stream (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}"); 51 out OutletLiquid as liquid_stream (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}"); 52 out OutletVapour as vapour_stream (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}"); 53 53 54 54 … … 65 65 EQUATIONS 66 66 "Component Molar Balance" 67 diff(M)=Inlet.F*Inlet.z + InletL .F*InletL.z + InletV.F*InletV.z OutletL.F*OutletL.z  OutletV.F*OutletV.z67 diff(M)=Inlet.F*Inlet.z + InletLiquid.F*InletLiquid.z + InletVapour.F*InletVapour.z OutletLiquid.F*OutletLiquid.z  OutletVapour.F*OutletVapour.z 68 68 LiquidSideStream.F*LiquidSideStream.zVapourSideStream.F*VapourSideStream.z; 69 69 70 70 "Molar Holdup" 71 M = ML*OutletL .z + MV*OutletV.z;71 M = ML*OutletLiquid.z + MV*OutletVapour.z; 72 72 73 73 "Mol fraction normalisation" 74 sum(OutletL .z)= 1.0;75 sum(OutletL .z)= sum(OutletV.z);74 sum(OutletLiquid.z)= 1.0; 75 sum(OutletLiquid.z)= sum(OutletVapour.z); 76 76 77 77 "Liquid Volume" 78 vL = PP.LiquidVolume(OutletL .T, OutletL.P, OutletL.z);78 vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z); 79 79 80 80 "Vapour Volume" 81 vV = PP.VapourVolume(OutletV .T, OutletV.P, OutletV.z);81 vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z); 82 82 83 83 "Chemical Equilibrium" 84 PP.LiquidFugacityCoefficient(OutletL .T, OutletL.P, OutletL.z)*OutletL.z = PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, yideal)*yideal;84 PP.LiquidFugacityCoefficient(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z)*OutletLiquid.z = PP.VapourFugacityCoefficient(OutletVapour.T, OutletVapour.P, yideal)*yideal; 85 85 86 86 "Thermal Equilibrium" 87 OutletV .T = OutletL.T;87 OutletVapour.T = OutletLiquid.T; 88 88 89 89 "Mechanical Equilibrium" 90 OutletV .P = OutletL.P;90 OutletVapour.P = OutletLiquid.P; 91 91 92 92 "Thermal Equilibrium Vapour Side Stream" 93 OutletV .T = VapourSideStream.T;93 OutletVapour.T = VapourSideStream.T; 94 94 95 95 "Thermal Equilibrium Liquid Side Stream" 96 OutletL .T = LiquidSideStream.T;96 OutletLiquid.T = LiquidSideStream.T; 97 97 98 98 "Mechanical Equilibrium Vapour Side Stream" 99 OutletV .P= VapourSideStream.P;99 OutletVapour.P= VapourSideStream.P; 100 100 101 101 "Mechanical Equilibrium Liquid Side Stream" 102 OutletL .P = LiquidSideStream.P;102 OutletLiquid.P = LiquidSideStream.P; 103 103 104 104 "Composition Liquid Side Stream" 105 OutletL .z= LiquidSideStream.z;105 OutletLiquid.z= LiquidSideStream.z; 106 106 107 107 "Composition Vapour Side Stream" 108 OutletV .z= VapourSideStream.z;108 OutletVapour.z= VapourSideStream.z; 109 109 110 110 end … … 120 120 * the Liquid inlet stream 121 121 * the Vapour inlet stream 122 * the Vapour outlet flow (OutletV .F)122 * the Vapour outlet flow (OutletVapour.F) 123 123 124 124 == Initial == 125 * the plate temperature (OutletL .T)126 * the liquid height (Level) OR the liquid flow OutletL .F127 * (NoComps  1) OutletL compositions125 * the plate temperature (OutletLiquid.T) 126 * the liquid height (Level) OR the liquid flow OutletLiquid.F 127 * (NoComps  1) OutletLiquid compositions 128 128 129 129 == Options == … … 149 149 150 150 "Liquid Density" 151 rhoL = PP.LiquidDensity(OutletL .T, OutletL.P, OutletL.z);151 rhoL = PP.LiquidDensity(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z); 152 152 153 153 "Vapour Density" 154 rhoV = PP.VapourDensity(InletV .T, InletV.P, InletV.z);154 rhoV = PP.VapourDensity(InletVapour.T, InletVapour.P, InletVapour.z); 155 155 156 156 end … … 176 176 * the Liquid inlet stream; 177 177 * the Vapour inlet stream; 178 * the Vapour outlet flow (OutletV .F);178 * the Vapour outlet flow (OutletVapour.F); 179 179 * the reaction related variables. 180 180 181 181 == Initial == 182 * the plate temperature (OutletL .T)183 * the liquid height (Level) OR the liquid flow OutletL .F184 * (NoComps  1) OutletL compositions182 * the plate temperature (OutletLiquid.T) 183 * the liquid height (Level) OR the liquid flow OutletLiquid.F 184 * (NoComps  1) OutletLiquid compositions 185 185 "; 186 186 … … 197 197 198 198 199 in InletL as stream (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}");200 in InletV as stream (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}");201 out OutletL as liquid_stream (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}");202 out OutletV as vapour_stream (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}");199 in InletLiquid as stream (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}"); 200 in InletVapour as stream (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}"); 201 out OutletLiquid as liquid_stream (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}"); 202 out OutletVapour as vapour_stream (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}"); 203 203 204 204 yideal(NComp) as fraction; … … 221 221 222 222 "Molar Concentration" 223 OutletL .z = vL * C;223 OutletLiquid.z = vL * C; 224 224 225 225 "Reaction" 226 r3 = exp(7150*'K'/OutletL .T)*(4.85e4*C(1)*C(2)  1.23e4*C(3)*C(4))*'l/mol/s';226 r3 = exp(7150*'K'/OutletLiquid.T)*(4.85e4*C(1)*C(2)  1.23e4*C(3)*C(4))*'l/mol/s'; 227 227 228 228 "Molar Holdup" 229 M = ML*OutletL .z + MV*OutletV.z;229 M = ML*OutletLiquid.z + MV*OutletVapour.z; 230 230 231 231 "Thermal Equilibrium Vapour Side Stream" 232 OutletV .T = VapourSideStream.T;232 OutletVapour.T = VapourSideStream.T; 233 233 234 234 "Thermal Equilibrium Liquid Side Stream" 235 OutletL .T = LiquidSideStream.T;235 OutletLiquid.T = LiquidSideStream.T; 236 236 237 237 "Mechanical Equilibrium Vapour Side Stream" 238 OutletV .P= VapourSideStream.P;238 OutletVapour.P= VapourSideStream.P; 239 239 240 240 "Mechanical Equilibrium Liquid Side Stream" 241 OutletL .P = LiquidSideStream.P;241 OutletLiquid.P = LiquidSideStream.P; 242 242 243 243 "Composition Liquid Side Stream" 244 OutletL .z= LiquidSideStream.z;244 OutletLiquid.z= LiquidSideStream.z; 245 245 246 246 "Composition Vapour Side Stream" 247 OutletV .z= VapourSideStream.z;247 OutletVapour.z= VapourSideStream.z; 248 248 249 249 "Mol fraction normalisation" 250 sum(OutletL .z)= 1.0;250 sum(OutletLiquid.z)= 1.0; 251 251 252 252 "Liquid Volume" 253 vL = PP.LiquidVolume(OutletL .T, OutletL.P, OutletL.z);253 vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z); 254 254 255 255 "Vapour Volume" 256 vV = PP.VapourVolume(OutletV .T, OutletV.P, OutletV.z);256 vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z); 257 257 258 258 "Thermal Equilibrium" 259 OutletV .T = OutletL.T;259 OutletVapour.T = OutletLiquid.T; 260 260 261 261 "Mechanical Equilibrium" 262 OutletV .P = OutletL.P;262 OutletVapour.P = OutletLiquid.P; 263 263 264 264 Vol = ML*vL; 265 265 266 266 "Liquid Density" 267 rhoL = PP.LiquidDensity(OutletL .T, OutletL.P, OutletL.z);267 rhoL = PP.LiquidDensity(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z); 268 268 269 269 "Vapour Density" 270 rhoV = PP.VapourDensity(InletV .T, InletV.P, InletV.z);270 rhoV = PP.VapourDensity(InletVapour.T, InletVapour.P, InletVapour.z); 271 271 272 272 "Chemical Equilibrium" 273 PP.LiquidFugacityCoefficient(OutletL .T, OutletL.P, OutletL.z)*OutletL.z = PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, yideal)*yideal;274 275 sum(OutletL .z)= sum(OutletV.z);273 PP.LiquidFugacityCoefficient(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z)*OutletLiquid.z = PP.VapourFugacityCoefficient(OutletVapour.T, OutletVapour.P, yideal)*yideal; 274 275 sum(OutletLiquid.z)= sum(OutletVapour.z); 276 276 277 277 end … … 293 293 294 294 == Initial == 295 * the plate temperature (OutletL .T)295 * the plate temperature (OutletLiquid.T) 296 296 * the liquid molar holdup ML 297 * (NoComps  1) OutletL compositions297 * (NoComps  1) OutletLiquid compositions 298 298 "; 299 299 … … 313 313 314 314 Inlet as stream (Brief="Feed stream", PosX=0, PosY=0.4932, Symbol="_{in}"); 315 in InletL as stream (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}");316 in InletV as stream (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}");317 out OutletL as liquid_stream (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}");318 out OutletV as vapour_stream (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}");315 in InletLiquid as stream (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}"); 316 in InletVapour as stream (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}"); 317 out OutletLiquid as liquid_stream (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}"); 318 out OutletVapour as vapour_stream (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}"); 319 319 320 320 M(NComp) as mol (Brief="Molar Holdup in the tray", Default=0.01, Lower=0, Upper=100); … … 346 346 347 347 "Component Molar Balance" 348 diff(M)=Inlet.F*Inlet.z + InletL .F*InletL.z + InletV.F*InletV.z OutletL.F*OutletL.z  OutletV.F*OutletV.z;348 diff(M)=Inlet.F*Inlet.z + InletLiquid.F*InletLiquid.z + InletVapour.F*InletVapour.z OutletLiquid.F*OutletLiquid.z  OutletVapour.F*OutletVapour.z; 349 349 350 350 "Molar Holdup" 351 M = ML*OutletL .z + MV*OutletV.z;351 M = ML*OutletLiquid.z + MV*OutletVapour.z; 352 352 353 353 "Mol fraction normalisation" 354 sum(OutletL .z)= 1.0;354 sum(OutletLiquid.z)= 1.0; 355 355 356 356 "Liquid Volume" 357 vL = PP.LiquidVolume(OutletL .T, OutletL.P, OutletL.z);357 vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z); 358 358 359 359 "Vapour Volume" 360 vV = PP.VapourVolume(OutletV .T, OutletV.P, OutletV.z);360 vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z); 361 361 362 362 "Chemical Equilibrium" 363 PP.LiquidFugacityCoefficient(OutletL .T, OutletL.P, OutletL.z)*OutletL.z = PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, OutletV.z)*OutletV.z;363 PP.LiquidFugacityCoefficient(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z)*OutletLiquid.z = PP.VapourFugacityCoefficient(OutletVapour.T, OutletVapour.P, OutletVapour.z)*OutletVapour.z; 364 364 365 365 "Thermal Equilibrium" 366 OutletV .T = OutletL.T;366 OutletVapour.T = OutletLiquid.T; 367 367 368 368 "Mechanical Equilibrium" 369 OutletL .P = OutletV.P;369 OutletLiquid.P = OutletVapour.P; 370 370 371 371 "Liquid Density" 372 rhoL = PP.LiquidDensity(OutletL .T, OutletL.P, OutletL.z);372 rhoL = PP.LiquidDensity(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z); 373 373 374 374 "Vapour Density" 375 rhoV = PP.VapourDensity(InletV .T, InletV.P, InletV.z);375 rhoV = PP.VapourDensity(InletVapour.T, InletVapour.P, InletVapour.z); 376 376 377 377 "Liquid viscosity" 378 miL = PP.LiquidViscosity(OutletL .T, OutletL.P, OutletL.z);378 miL = PP.LiquidViscosity(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z); 379 379 380 380 "Vapour viscosity" 381 miV = PP.VapourViscosity(InletV .T, InletV.P, InletV.z);381 miV = PP.VapourViscosity(InletVapour.T, InletVapour.P, InletVapour.z); 382 382 383 383 "Volume flow rate of liquid, m^3/m^2/s" 384 uL * Al = OutletL .F * vL;385 386 deltaP = InletV .P  OutletV.P;384 uL * Al = OutletLiquid.F * vL; 385 386 deltaP = InletVapour.P  OutletVapour.P; 387 387 388 388 end … … 458 458 in Inlet as stream (Brief="Feed stream", PosX=0, PosY=0.4932, Symbol="_{in}"); 459 459 in InletFV as stream (Brief="Feed stream", PosX=0, PosY=0.4932, Symbol="_{in}"); 460 in InletL as stream (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}");461 in InletV as stream (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}");462 out OutletL as liquid_stream (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}");463 out OutletV as vapour_stream (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}");460 in InletLiquid as stream (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}"); 461 in InletVapour as stream (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}"); 462 out OutletLiquid as liquid_stream (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}"); 463 out OutletVapour as vapour_stream (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}"); 464 464 465 465 M_liq(NComp) as mol (Brief="Liquid Molar Holdup in the tray"); … … 479 479 EQUATIONS 480 480 "Component Molar Balance" 481 diff(M_liq)=Inlet.F*Inlet.z + InletL .F*InletL.z482  OutletL .F*OutletL.z + interf.NL;483 484 diff(M_vap)=InletFV.F*InletFV.z + InletV .F*InletV.z485  OutletV .F*OutletV.z  interf.NV;481 diff(M_liq)=Inlet.F*Inlet.z + InletLiquid.F*InletLiquid.z 482  OutletLiquid.F*OutletLiquid.z + interf.NL; 483 484 diff(M_vap)=InletFV.F*InletFV.z + InletVapour.F*InletVapour.z 485  OutletVapour.F*OutletVapour.z  interf.NV; 486 486 487 487 "Energy Balance" 488 diff(E_liq) = Inlet.F*Inlet.h + InletL .F*InletL.h489  OutletL .F*OutletL.h + Q + interf.E_liq;490 491 diff(E_vap) = InletFV.F*InletFV.h + InletV .F*InletV.h492  OutletV .F*OutletV.h  interf.E_vap;488 diff(E_liq) = Inlet.F*Inlet.h + InletLiquid.F*InletLiquid.h 489  OutletLiquid.F*OutletLiquid.h + Q + interf.E_liq; 490 491 diff(E_vap) = InletFV.F*InletFV.h + InletVapour.F*InletVapour.h 492  OutletVapour.F*OutletVapour.h  interf.E_vap; 493 493 494 494 "Molar Holdup" 495 M_liq = ML*OutletL .z;496 497 M_vap = MV*OutletV .z;495 M_liq = ML*OutletLiquid.z; 496 497 M_vap = MV*OutletVapour.z; 498 498 499 499 "Energy Holdup" 500 E_liq = ML*(OutletL .h  OutletL.P*vL);501 502 E_vap = MV*(OutletV .h  OutletV.P*vV);500 E_liq = ML*(OutletLiquid.h  OutletLiquid.P*vL); 501 502 E_vap = MV*(OutletVapour.h  OutletVapour.P*vV); 503 503 504 504 "Energy Rate through the interface" 505 interf.E_liq = interf.htL*interf.a*(interf.TOutletL .T)+sum(interf.NL)*interf.hL;506 507 interf.E_vap = interf.htV*interf.a*(OutletV .Tinterf.T)+sum(interf.NV)*interf.hV;505 interf.E_liq = interf.htL*interf.a*(interf.TOutletLiquid.T)+sum(interf.NL)*interf.hL; 506 507 interf.E_vap = interf.htV*interf.a*(OutletVapour.Tinterf.T)+sum(interf.NV)*interf.hV; 508 508 509 509 "Mass Conservation" … … 514 514 515 515 "Mol fraction normalisation" 516 sum(OutletL .z)= 1.0;517 sum(OutletL .z)= sum(OutletV.z);516 sum(OutletLiquid.z)= 1.0; 517 sum(OutletLiquid.z)= sum(OutletVapour.z); 518 518 sum(interf.x)=1.0; 519 519 sum(interf.x)=sum(interf.y); 520 520 521 521 "Liquid Volume" 522 vL = PP.LiquidVolume(OutletL .T, OutletL.P, OutletL.z);522 vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z); 523 523 "Vapour Volume" 524 vV = PP.VapourVolume(OutletV .T, OutletV.P, OutletV.z);524 vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z); 525 525 526 526 "Chemical Equilibrium" … … 535 535 536 536 "Total Mass Transfer Rates" 537 interf.NL(1:NC1)=interf.a*sumt(interf.kL*(interf.x(1:NC1)OutletL .z(1:NC1)))/vL+538 OutletL .z(1:NC1)*sum(interf.NL);537 interf.NL(1:NC1)=interf.a*sumt(interf.kL*(interf.x(1:NC1)OutletLiquid.z(1:NC1)))/vL+ 538 OutletLiquid.z(1:NC1)*sum(interf.NL); 539 539 540 540 # interf.NL(1:NC1)=0.01*'kmol/s'; 541 541 542 interf.NV(1:NC1)=interf.a*sumt(interf.kV*(OutletV .z(1:NC1)interf.y(1:NC1)))/vV+543 OutletV .z(1:NC1)*sum(interf.NV);542 interf.NV(1:NC1)=interf.a*sumt(interf.kV*(OutletVapour.z(1:NC1)interf.y(1:NC1)))/vV+ 543 OutletVapour.z(1:NC1)*sum(interf.NV); 544 544 545 545 "Mechanical Equilibrium" 546 OutletV .P = OutletL.P;547 interf.P=OutletL .P;546 OutletVapour.P = OutletLiquid.P; 547 interf.P=OutletLiquid.P; 548 548 end 549 549 … … 558 558 * the Liquid inlet stream 559 559 * the Vapour inlet stream 560 * the Vapour outlet flow (OutletV .F)560 * the Vapour outlet flow (OutletVapour.F) 561 561 562 562 == Initial == 563 * the plate temperature of both phases (OutletL .T and OutletV.T)563 * the plate temperature of both phases (OutletLiquid.T and OutletVapour.T) 564 564 * the liquid height (Level) OR the liquid flow holdup (ML) 565 565 * the vapor holdup (MV) 566 * (NoComps  1) OutletL compositions566 * (NoComps  1) OutletLiquid compositions 567 567 "; 568 568 … … 584 584 EQUATIONS 585 585 "Liquid Density" 586 rhoL = PP.LiquidDensity(OutletL .T, OutletL.P, OutletL.z);586 rhoL = PP.LiquidDensity(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z); 587 587 "Vapour Density" 588 rhoV = PP.VapourDensity(InletV .T, InletV.P, InletV.z);588 rhoV = PP.VapourDensity(InletVapour.T, InletVapour.P, InletVapour.z); 589 589 590 590 switch LiquidFlow 591 591 case "on": 592 592 "Francis Equation" 593 # OutletL .F*vL = 1.84*'m^0.5/s'*lw*((Level(beta*hw))/(beta))^1.5;594 OutletL .F*vL = 1.84*'1/s'*lw*((Level(beta*hw))/(beta))^2;593 # OutletLiquid.F*vL = 1.84*'m^0.5/s'*lw*((Level(beta*hw))/(beta))^1.5; 594 OutletLiquid.F*vL = 1.84*'1/s'*lw*((Level(beta*hw))/(beta))^2; 595 595 when Level < (beta * hw) switchto "off"; 596 596 597 597 case "off": 598 598 "Low level" 599 OutletL .F = 0 * 'mol/h';599 OutletLiquid.F = 0 * 'mol/h'; 600 600 when Level > (beta * hw) + 1e6*'m' switchto "on"; 601 601 end … … 603 603 switch VapourFlow 604 604 case "on": 605 InletV .F*vV = sqrt((InletV.P  OutletV.P)/(rhoV*alfa))*Ah;606 when InletV .F < 1e6 * 'kmol/h' switchto "off";605 InletVapour.F*vV = sqrt((InletVapour.P  OutletVapour.P)/(rhoV*alfa))*Ah; 606 when InletVapour.F < 1e6 * 'kmol/h' switchto "off"; 607 607 608 608 case "off": 609 InletV .F = 0 * 'mol/s';610 when InletV .P > OutletV.P + Level*g*rhoL + 1e1 * 'atm' switchto "on";609 InletVapour.F = 0 * 'mol/s'; 610 when InletVapour.P > OutletVapour.P + Level*g*rhoL + 1e1 * 'atm' switchto "on"; 611 611 end 612 612 end
Note: See TracChangeset
for help on using the changeset viewer.