Changeset 918 for branches/gui/eml/stage_separators
- Timestamp:
- Feb 25, 2010, 1:45:10 PM (14 years ago)
- Location:
- branches/gui/eml/stage_separators
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/gui/eml/stage_separators/flash.mso
r879 r918 19 19 using "tank"; 20 20 21 Model flash 21 Model flash as VesselVolume 22 22 23 23 ATTRIBUTES … … 61 61 62 62 Mw(NComp) as molweight (Brief="Mol Weight", Hidden=true); 63 Gconst as acceleration (Brief="Gravity Acceleration",Default=9.81,Hidden=true);64 63 65 64 Levelpercent_Initial as positive (Brief="Initial liquid height in Percent", Default = 0.70); … … 71 70 Mw=PP.MolecularWeight(); 72 71 73 Gconst = 9.81 * 'm/(s^2)';74 75 72 VARIABLES 76 77 Geometry as VesselVolume (Brief="Vessel Geometry", Symbol=" ");78 73 79 74 in Inlet as stream (Brief="Feed Stream", PosX=0, PosY=0.48, Symbol="_{in}"); … … 123 118 124 119 "Energy Holdup" 125 E = LiquidHoldup*OutletLiquid.h + VapourHoldup*OutletVapour.h - OutletLiquid.P* Geometry.Vtotal;120 E = LiquidHoldup*OutletLiquid.h + VapourHoldup*OutletVapour.h - OutletLiquid.P*Vtotal; 126 121 127 122 "Mol fraction normalisation" … … 151 146 152 147 "Static Head" 153 Pstatic = PP.LiquidDensity(OutletLiquid.T, Peq, OutletLiquid.z) * Gconst * Geometry.Level;148 Pstatic = PP.LiquidDensity(OutletLiquid.T, Peq, OutletLiquid.z) * Gconst * Level; 154 149 155 150 "Mechanical Equilibrium for the Liquid Phase" … … 163 158 164 159 "Geometry Constraint" 165 Geometry.Vtotal = LiquidHoldup * vL + VapourHoldup * vV;160 Vtotal = LiquidHoldup * vL + VapourHoldup * vV; 166 161 167 162 "Temperature indicator" … … 172 167 173 168 "Level indicator" 174 LI* Geometry.Vtotal= Geometry.Vfilled;169 LI*Vtotal= Vfilled; 175 170 176 171 "Liquid Level" 177 LiquidHoldup * vL = Geometry.Vfilled;172 LiquidHoldup * vL = Vfilled; 178 173 179 174 end -
branches/gui/eml/stage_separators/tank.mso
r913 r918 37 37 38 38 PARAMETERS 39 P Ias positive (Brief="Pi value", Default=3.141593,Hidden=true, Symbol="\pi");39 Pi as positive (Brief="Pi value", Default=3.141593,Hidden=true, Symbol="\pi"); 40 40 Gconst as acceleration (Brief="Gravity Acceleration",Default=9.81,Hidden=true); 41 41 … … 53 53 54 54 Gconst = 9.81 * 'm/(s^2)'; 55 Vhead_elliptical = (P I*Diameter^3)/12;56 Vhead_hemispherical = (P I*Diameter^3)/6;57 Vcylinder = 0.25*(P I*Diameter^2)*Lenght;55 Vhead_elliptical = (Pi*Diameter^3)/12; 56 Vhead_hemispherical = (Pi*Diameter^3)/6; 57 Vcylinder = 0.25*(Pi*Diameter^2)*Lenght; 58 58 radius = 0.5*Diameter; 59 59 … … 72 72 73 73 "Vessel Cross Section Area" 74 Across = 0.25*(P I*Diameter^2);74 Across = 0.25*(Pi*Diameter^2); 75 75 76 76 switch Heads … … 84 84 85 85 "Vessel Filled Volume" 86 Vfilled = 0.25*P I*(((Diameter*Level)/(0.25*Diameter))^2)*(0.25*Diameter-Level/3);86 Vfilled = 0.25*Pi*(((Diameter*Level)/(0.25*Diameter))^2)*(0.25*Diameter-Level/3); 87 87 88 88 else 89 89 90 90 "Vessel Filled Volume" 91 Vfilled = 0.25*P I*(Diameter^2)*(Level - 0.25*Diameter/3);91 Vfilled = 0.25*Pi*(Diameter^2)*(Level - 0.25*Diameter/3); 92 92 93 93 end … … 101 101 102 102 "Vessel Filled Volume" 103 Vfilled = 0.25*P I*(Level^2)*(2*Diameter-4*Level/3);103 Vfilled = 0.25*Pi*(Level^2)*(2*Diameter-4*Level/3); 104 104 105 105 else 106 106 107 107 "Vessel Filled Volume" 108 Vfilled = 0.25*P I*((2/3)*((0.5*Diameter)^3) - (0.25*(Diameter)^3) + Level*Diameter^2);108 Vfilled = 0.25*Pi*((2/3)*((0.5*Diameter)^3) - (0.25*(Diameter)^3) + Level*Diameter^2); 109 109 110 110 end
Note: See TracChangeset
for help on using the changeset viewer.