Changeset 599 for branches/gui/eml/pressure_changers
- Timestamp:
- Aug 15, 2008, 10:31:59 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/gui/eml/pressure_changers/compressor.mso
r598 r599 21 21 22 22 Model centrifugal_compressor 23 24 # Do Not Forget to Update The Losses by Mechanical Efficiency !!!! 23 25 24 26 ATTRIBUTES … … 38 40 VARIABLES 39 41 40 PolyCoeff as positive (Brief = "Polytropic Coefficient", Lower= 0);41 IseCoeff as positive (Brief = "Isentropic Coefficient", Lower=1 e-3);42 Cp as cp_mol (Brief = "Heat Capacity ");43 Cv as cv_mol (Brief = "Heat Capacity ");42 PolyCoeff as positive (Brief = "Polytropic Coefficient", Lower=1E-6); 43 IseCoeff as positive (Brief = "Isentropic Coefficient", Lower=1E-6); 44 Cp as cp_mol (Brief = "Heat Capacity at constant Pressure"); 45 Cv as cv_mol (Brief = "Heat Capacity at constant Volume"); 44 46 Pratio as positive (Brief = "Pressure Ratio", Symbol ="P_{ratio}"); 45 47 Pdrop as press_delta (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P"); 48 Pincrease as press_delta (Brief = "Pressure Increase", DisplayUnit = 'kPa', Symbol ="P_{incr}"); 46 49 47 Head as energy_mass (Brief = " IsentropicHead");50 Head as energy_mass (Brief = "Head"); 48 51 Tiso as temperature (Brief = "Isentropic Temperature"); 49 52 … … 86 89 Outlet.P = Inlet.P - Pdrop; 87 90 91 "Pressure Increase" 92 Outlet.P = Inlet.P + Pincrease; 93 88 94 "Mass Density" 89 95 rho = PP.VapourDensity(Inlet.T, Inlet.P, Inlet.z); … … 104 110 IseCoeff * Cv = Cp; 105 111 106 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++107 112 switch CompressorType 108 113 … … 110 115 111 116 "Isentropic Head" 112 Head = (0.5*Zfac_in+0.5*Zfac_out)*(1/Mwm)*(IseCoeff/(IseCoeff-1 ))*Rgas*Inlet.T*((Outlet.P/Inlet.P)^((IseCoeff-1)/IseCoeff) - 1);117 Head = (0.5*Zfac_in+0.5*Zfac_out)*(1/Mwm)*(IseCoeff/(IseCoeff-1.001))*Rgas*Inlet.T*((Outlet.P/Inlet.P)^((IseCoeff-1.001)/IseCoeff) - 1); 113 118 114 119 "Fluid Power" 115 120 FluidPower*IsentropicEff = Head*sum(Mw*Inlet.z)*Inlet.F; 116 121 122 "Discharge Temperature" 123 Outlet.T = Inlet.T*((Outlet.P/Inlet.P)^((IseCoeff-1.001)/IseCoeff)); 124 117 125 case "Polytropic Operation": 118 126 119 127 "Polytropic Head" 120 Head = (0.5*Zfac_in+0.5*Zfac_out)*(1/Mwm)*(PolyCoeff/(PolyCoeff-1 ))*Rgas*Inlet.T*((Outlet.P/Inlet.P)^((PolyCoeff-1)/PolyCoeff) - 1);128 Head = (0.5*Zfac_in+0.5*Zfac_out)*(1/Mwm)*(PolyCoeff/(PolyCoeff-1.001))*Rgas*Inlet.T*((Outlet.P/Inlet.P)^((PolyCoeff-1.001)/PolyCoeff) - 1); 121 129 122 130 "Fluid Power" 123 131 FluidPower*PolytropicEff = Head*sum(Mw*Inlet.z)*Inlet.F; 132 133 "Discharge Temperature" 134 Outlet.T = Inlet.T*((Outlet.P/Inlet.P)^((PolyCoeff-1.001)/PolyCoeff)); 124 135 125 136 end … … 128 139 PP.VapourEntropy(Tiso, Outlet.P, Outlet.z) = PP.VapourEntropy(Inlet.T, Inlet.P, Inlet.z); 129 140 130 "Discharge Temperature"131 IsentropicEff * (Outlet.T- Inlet.T) = (Tiso - Inlet.T);141 #"Calculate Politropic Coefficient" 142 #PolyCoeff*(ln(Outlet.T/Inlet.T)) = (PolyCoeff-1)*(ln(Outlet.P/Inlet.P)); 132 143 133 "Calculate Politropic Coefficient" 134 PolyCoeff*(ln(Outlet.T/Inlet.T)) = (PolyCoeff-1)*(ln(Outlet.P/Inlet.P)); 135 136 #"Polytropic Efficiency" 137 #PolytropicEff * (PolyCoeff-1) * IseCoeff = PolyCoeff * (IseCoeff-1); 144 "Polytropic Efficiency" 145 PolytropicEff * (PolyCoeff-1) * IseCoeff = PolyCoeff * (IseCoeff-1); 138 146 139 147 "Fluid Power"
Note: See TracChangeset
for help on using the changeset viewer.