# Example 6.1 d0 = 0.2 m A0 = pi() * d0^2 / 4 ; m^2 v0 = 0.1 m/s T0 = 25.0 degC P0 = 150.0 kPa R = 0.2870 kJ/kg degK n = R * T0 / P0 ; m^3/kg mdot = v0 * A0 / n ; kg/s # Example 6.2 P = 600.0 kPa m0 = 1.0 kg n = 0.001 m^3/kg W0 = P * m0 * n ; kJ # Example 6.3 Pi0 = 1.0 Mpa Ti = 60.0 degC mdotr = 0.2 kg/s Pe0 = 0.950 Mpa Te = 35.0 degC hi_r = 441.89 kJ/kg he_r = 249.10 kJ/kg hi_w = 42.00 kJ/kg he_w = 83.95 kJ/kg mdotw = mdotr * (hi_r - he_r) / (he_w - hi_w) ; kg/s QdotCV = mdotr * (he_r - hi_r) ; kW QdotCV = mdotw * (he_w - hi_r) ; kW # Example 6.4 Pii = 0.6 Mpa Ti = 200.0 degC vi = 50.0 m/s Pe = 0.150 Mpa ve = 600.0 m/s hi = 2850.10 kJ/kg he = hi + (vi^2 - ve^2) / 2 ; kJ/kg hf = 467.10 kJ/kg hfg = 2226.50 kJ/kg x = (he - hf) / hfg # Example 6.4E Pii = 100.0 psi Ti = 400.0 degF vi = 200.0 ft/s Pe = 20.800 psi ve = 2000.0 ft/s hi = 1227.50 Btu/lbm he = hi + (vi^2 - ve^2) / 2 ; Btu/lbm hf = 198.31 Btu/lbm hfg = 958.81 Btu/lbm x = (he - hf) / hfg # Example 6.5 Pii = 1.5 Mpa Ti = 35.0 degC Pe = 291.0 kPa hi = 346.800 kJ/kg he = hi ; kJ/kg hf = 134.40 kJ/kg hfg = 1296.40 kJ/kg x = (he - hf) / hfg # Example 6.6 mdot = 1.5 kg/s Pi0 = 2.0 Mpa Ti = 350.0 degC hi = 3137.000 kJ/kg Zi = 6.000 m Vi = 50.000 m/s Pe = 0.100 Mpa xe = 1.000 he = 2675.500 kJ/kg Ze = 3.000 m Ve = 100.000 m/s QdotCV = -8.50 kW WdotCV = QdotCV + mdot * (hi + Vi^2 / 2 + grav() * Zi) \ - mdot * (he + Ve^2 / 2 + grav() *Ze) WdotCV / (1 kg / s); kJ/kg # Example 6.7 q12 = 0.0 kJ/kg Wdot = -50.0 kW P1 = 100.0 kPa T1 = 280.0 degK h1 = 198.000 kJ/kg V1 = 0.000 m/s P2 = 1100.000 kPa h2 = 401.500 kJ/kg V2 = 25.000 m/s h3 = 257.900 kJ/kg V3 = 25.000 m/s w = q12 + h1 - h2 + (V1^2 - V2^2) / 2 ; kJ/kg mdot = Wdot / w ; kg/s q23 = h3 - h2 + (V3^2 - V2^2) / 2 ; kJ/kg Qdotcool = -mdot * q23 ; kW # Example 6.8 mdot = 1.5 kg/s Pi0 = 90.0 kPa Ti = 10.0 degC Zi = 0.000 m d = 0.040 m Pe = 400.000 kPa Te = 10.0 degC Ze = -15.000 m n = 0.001001 m^3/kg Pcorr = 101.30 kPa Wdot = mdot * (grav() * (Ze - Zi) - (Pe + Pcorr - Pi0) * n) ; W # Example 6.9 w45 = -4.0 kJ/kg P1 = 2.0 Mpa T1 = 300.0 degC h1 = 3023.5 kJ/kg P2 = 1.900 Mpa T2 = 290.0 degC h2 = 3002.5 kJ/kg P3 = 15.000 kPa x3 = 0.9 hf = 226.0 kJ/kg hfg = 2373.1 kJ/kg h3 = hf + x3 * hfg ; kJ/kg P4 = 14.000 kPa T4 = 45.0 degC h4 = 188.5 kJ/kg q12 = h2 - h1 ; kJ/kg w23 = h2 - h3 ; kJ/kg q34 = h4 - h3 ; kJ/kg h5 = h4 - w45 ; kJ/kg q15 = h1 - h5 ; kJ/kg # Example 6.10 W0 = 5.0 kW mdot = 0.1 kg/s P1 = 100.0 kPa T1 = -20.0 degC h1 = 387.2 kJ/kg P2 = 800.000 kPa T2 = 50.0 degC h2 = 435.1 kJ/kg T3 = 30.000 degC x3 = 0.0 h3 = 241.8 kJ/kg T4 = -25.000 degC hf = 167.4 kJ/kg hfg = 215.6 kJ/kg h4 = h3 ; kJ/kg x4 = (h4 - hf) / hfg Qdot = mdot * (h1 - h4) ; kW Qdot = mdot * (h2 - h1) - W0 ; kW # Example 6.11 P = 1.4 Mpa T1 = 300.0 degC u2 = 3040.4 kJ/kg h1 = u2 T2 = 452.0 degC # Example 6.12 P1 = 350.0 kPa T1 = 300.0 degC V0 = 0.4 m^3 P2 = 1.4 Mpa n1 = 0.524 m^3/kg m1 = V0 / n1 ; kg h1 = 3040.4 kJ/kg u1 = 2548.9 kJ/kg T2 = 300.0 degC n2 = 0.1823 m^3/kg u2 = 2785.2 kJ/kg e2 = V0 * (h1 - u2) / n2 - m1 * (h1 - u1) ; kJ T3 = 350.0 degC n3 = 0.2 m^3/kg u3 = 2869.1 kJ/kg e3 = V0 * (h1 - u3) / n3 - m1 * (h1 - u1) ; kJ T4 = (abs(e3) * T2 + abs(e2) * T3) / (abs(e2) + abs(e3)) ; degC n4 = 0.197 m^3/kg m2 = V0 / n4 ; kg mi0 = m2 - m1 ; kg ## Example 6.13 #BEGIN # V0 = 2.0 m^3 # T1 = 40.0 degC # fliq = 0.5 # fvap = 0.5 # T2 = 10.0 degC # hg1 = 1470.2 kJ/kg # hg2 = 1452.0 kJ/kg # hg_ave = 0.5 * (hg1 + hg2) ; kJ/kg # nf1 = 0.00172500 m^3/kg # nf2 = 0.00160000 m^3/kg # uf1 = 368.7 kJ/kg # uf2 = 226.0 kJ/kg # ng1 = 0.0831300 m^3/kg # nfg2 = 0.2038100 m^3/kg # ug1 = 1341.0 kJ/kg # ufg2 = 1099.7 kJ/kg # mf1 = fliq * V0 / nf1 ; kg # mg1 = fvap * V0 / ng1 ; kg # m1 = mf1 + mg1 ; kg # Dmh = m1 * hg_ave - mf1 * uf1 - mg1 * ug1 ; kJ # x2 = 0.05 # e = V0 * (hg_ave - uf2 - x2 * ufg2) / (nf2 + nfg2 * x2) - Dmh ; kJ # n2 = nf2 + x2 * nfg2 ; m^3/kg # m2 = V0 / n2 ; kg # me = m1 - m2 ; kg #END #solve(e, x2) ## Example 6.14 #BEGIN # V0 = 50.0 ft^3 # T1 = 100.0 degF # fliq = 0.5 # fvap = 0.5 # T2 = 50.0 degF # hg1 = 631.8 Btu/lbm # hg2 = 624.3 Btu/lbm # hg_ave = 0.5 * (hg1 + hg2) ; Btu/lbm # nf1 = 0.02747000 ft^3/lbm # nf2 = 0.02564000 ft^3/lbm # uf1 = 153.9 Btu/lbm # uf2 = 97.2 Btu/lbm # ng1 = 1.4168000 ft^3/lbm # nfg2 = 3.2647000 ft^3/lbm # ug1 = 576.2 Btu/lbm # ufg2 = 472.8 Btu/lbm # mf1 = fliq * V0 / nf1 ; lbm # mg1 = fvap * V0 / ng1 ; lbm # m1 = mf1 + mg1 ; lbm # Dmh = m1 * hg_ave - mf1 * uf1 - mg1 * ug1 ; Btu # x2 = 0.05 # e = V0 * (hg_ave - uf2 - x2 * ufg2) / (nf2 + nfg2 * x2) - Dmh ; kJ # n2 = nf2 + x2 * nfg2 ; ft^3/lbm # m2 = V0 / n2 ; lbm # me = m1 - m2 ; lbm #END #solve(e, x2) # Homework Problems: TBD