Changeset 109
 Timestamp:
 Jan 11, 2007, 8:03:27 PM (16 years ago)
 Location:
 mso
 Files:

 2 edited
Legend:
 Unmodified
 Added
 Removed

mso/eml/pressure_changers/valve.mso
r108 r109 27 27 * * Isentalpic 28 28 * * Liquid 29 * 29 * 30 30 * 31 31 * Author: Estefane Horn, Núbia do Carmo Ferreira … … 54 54 55 55 SET 56 rho60F = 99 .022 * "kg/m^3";56 rho60F = 999.2 * "kg/m^3"; 57 57 58 58 EQUATIONS … … 69 69 Outlet.z = Inlet.z; 70 70 71 "Valve Equation" 72 #if Pdiff >= 0 then 73 # "Flow" 71 if Pdiff > 0 then 72 "Valve Equation  Flow" 74 73 Qv = fc*cv*sqrt(Pdiff/Gf); 75 #else76 # "Closed"77 # Pdiff < 0 * "kmol/h";78 #end74 else 75 "Valve Equation  Closed" 76 Qv = 0 * "m^3/h"; 77 end 79 78 80 79 "Calculate Gf" … … 92 91 end 93 92 94 #*95 * Model of a valve_linear96 **#97 #*  Specify:98 * * Opening fraction (x)99 *100 **#101 102 103 93 Model valve_linear as valve_basic 104 94 … … 113 103 end 114 104 115 #*116 * Model of a valve_parabolic117 **#118 #*  Specify:119 * * Opening fraction (x)120 *121 **#122 123 105 Model valve_parabolic as valve_basic 124 106 … … 136 118 end 137 119 138 #* 139 * Model of a valve_equal 140 **# 141 #*  Specify: 142 * * Opening fraction (x) 143 * 144 **# 145 146 Model valve_equal as valve_basic 120 Model valve_equal as valve_basic 147 121 148 122 PARAMETERS … … 158 132 159 133 end 160 161 #*162 * Model of a valve_quick163 **#164 #*  Specify:165 * * Opening fraction (x)166 *167 **#168 134 169 135 Model valve_quick as valve_basic … … 245 211 (1Inlet.v)*PP.LiquidVolume((Inlet.T+Outlet.T)/2, (Inlet.P+Outlet.P)/2, Outlet.z); 246 212 247 #if Inlet.P > Outlet.P then248 #"Flow"249 Outlet.F * v = cv*x*sqrt((Inlet.P  Outlet.P)*rho_ref / rho ) ;250 #else251 #"Closed"252 #Outlet.F = 0 * "kmol/h";253 #end254 end 213 if Inlet.P > Outlet.P then 214 "Flow" 215 Outlet.F * v = k*x*sqrt((Inlet.P  Outlet.P)*rho_ref / rho ) ; 216 else 217 "Closed" 218 Outlet.F = 0 * "kmol/h"; 219 end 220 end 
mso/sample/pressure_changers/sample_valve.mso
r81 r109 59 59 OPTIONS 60 60 relativeAccuracy = 1e6; 61 time = [0:20:1000];62 61 mode = "steady"; 63 64 62 end 65 63 … … 97 95 OPTIONS 98 96 relativeAccuracy = 1e6; 99 time = [0:20:1000];100 97 mode = "steady"; 101 102 98 end 103 99 … … 128 124 V1.Qv = 113.56 * "m^3/h"; 129 125 V1.x = 0.5; 130 V1.cv = 3.997 * "m^3/h/kPa^0.5"; 126 V1.cv = 10 * "m^3/h/kPa^0.5"; 127 128 OPTIONS 129 relativeAccuracy = 1e6; 130 mode = "steady"; 131 131 end
Note: See TracChangeset
for help on using the changeset viewer.