Changeset 46
- Timestamp:
- Nov 7, 2006, 2:47:55 PM (16 years ago)
- Location:
- mso
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
mso/eml/stage_separators/column.mso
r42 r46 167 167 "Pressure Drop through the tray" 168 168 trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - ttop.Outlet.P)/"Pa"))/2 * 169 trays(top).Ah/trays(top).vV * sqrt(2*(trays(top).OutletV.P -169 trays(top).Ah/trays(top).vV /2* sqrt(2*(trays(top).OutletV.P - 170 170 ttop.Outlet.P + 1e-8 * "atm") / (trays(top).alfa*trays(top).rhoV)); 171 171 … … 173 173 (trays([top+topdown:topdown:bot]).OutletV.P - 174 174 trays([top+topdown:topdown:bot]).InletL.P)/"Pa"))/2 * 175 trays([top+topdown:topdown:bot]).Ah/trays([top+topdown:topdown:bot]).vV *175 trays([top+topdown:topdown:bot]).Ah/trays([top+topdown:topdown:bot]).vV /2* 176 176 sqrt(2*(trays([top+topdown:topdown:bot]).OutletV.P - 177 177 trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * "atm") / -
mso/eml/stage_separators/reboiler.mso
r38 r46 139 139 "Vapourisation Fraction" 140 140 OutletV.v = 1.0; 141 end 142 143 Model reboilerSteady_fakeH 144 PARAMETERS 145 ext PP as CalcObject; 146 ext NComp as Integer; 147 DP as press_delta (Brief="Pressure Drop in the reboiler"); 148 k as Real (Brief = "Flow Constant", Unit="mol/J"); 149 150 VARIABLES 151 in InletL as stream; #(Brief="Liquid inlet stream"); 152 out OutletV as stream; #(Brief="Vapour outlet stream"); 153 in Q as heat_rate (Brief="Heat supplied"); 154 155 EQUATIONS 156 "Molar Balance" 157 InletL.F = OutletV.F; 158 InletL.z = OutletV.z; 159 160 "Energy Balance" 161 InletL.F*InletL.h + Q = OutletV.F*OutletV.h; 162 163 "Pressure" 164 DP = InletL.P - OutletV.P; 165 166 "Fake Vapourisation Fraction" 167 OutletV.v = 1.0; 168 169 "Fake output temperature" 170 OutletV.T = 300*"K"; 171 172 "Pressure Drop through the reboiler" 173 OutletV.F = k*Q; 141 174 end 142 175 -
mso/sample/controllers/sample_tanks_pid.mso
r29 r46 127 127 128 128 OPTIONS 129 NLASolver = "sundials"; 129 130 time = [0:0.1:30] * "h" ; 130 131 end -
mso/sample/miscellaneous/sample_pend.mso
r30 r46 64 64 time = [0:0.1:30]; 65 65 integration = "original"; # original, index0 or index1 66 accuracyType = "scalar"; # scalar or vector67 66 relativeAccuracy = 1e-8; 68 67 absoluteAccuracy = 1e-9; -
mso/sample/reactors/sample_pfr.mso
r43 r46 78 78 #mode = "steady"; 79 79 time = [0:0.05:3] * "s"; 80 NLASolver = "sundials"; 81 #DAESolver = "dassl"; 80 82 end -
mso/sample/stage_separators/sample_column.mso
r33 r46 73 73 OPTIONS 74 74 #relativeAccuracy = 1e-3; 75 NLASolver = "sundials"; 75 76 time = [0:10:1000]; 76 77 end -
mso/sample/stage_separators/sample_tank.mso
r42 r46 11 11 NComp = PP.NumberOfComponents; 12 12 13 VARIABLES 14 Qtank as heat_rate (Brief="Heat rate supplied to tank"); 15 13 16 DEVICES 14 17 t as tank; … … 17 20 CONNECTIONS 18 21 s to t.Inlet; 22 Qtank to t.Q; 19 23 20 24 SPECIFY … … 25 29 s.v = 0.698; 26 30 27 t.Q= 0 * "J/s";31 Qtank = 0 * "J/s"; 28 32 t.Outlet.F = 179 * "kmol/h"; 29 33 … … 50 54 PP.VapourModel = "PR"; 51 55 NComp = PP.NumberOfComponents; 56 57 VARIABLES 58 Qtank as heat_rate (Brief="Heat rate supplied to tank"); 52 59 53 60 DEVICES … … 57 64 CONNECTIONS 58 65 s to t.Inlet; 66 Qtank to t.Q; 59 67 60 68 SPECIFY … … 65 73 s.v = 0.368; 66 74 67 t.Q= 0 * "J/s";75 Qtank = 0 * "J/s"; 68 76 t.Outlet.F = 490 * "kmol/h"; 69 77
Note: See TracChangeset
for help on using the changeset viewer.