# Example 9.1     
 P1   =   1.0 Mpa
 T1   =   300.0 degC
 V1   =   50.0 m/s
 P2   =   150.0 kPa
 V2   =   200.0 m/s
 h1   =   3051.2 kJ/kg
 s1   =   7.1228 kJ/kg/degK
 s2 = s1  ; kJ/kg/degK 
 sf    =   1.4355 kJ/kg/degK
 sfg   =   5.7897 kJ/kg/degK
 x2 = (s2 - sf ) / sfg   
 hf   =   467.1 kJ/kg
 hfg   =   2226.5 kJ/kg
 h2 = hf + x2 * hfg  ; kJ/kg 
 w = h1 + V1^2 / 2 - h2 - V2^2 / 2  ; kJ/kg 

# Example 9.2     
 P1   =   1.0 Mpa
 T1   =   300.0 degC
 V1   =   30.0 m/s
 h1   =   3051.2 kJ/kg
 s1   =   7.1228 kJ/kg/degK
 P2   =   0.3 Mpa
 T2   =   159.1 degC
 h2   =   2780.2 kJ/kg
 V2 = sqrt(2 * (h1 - h2 + V1^2 / 2))  ; m/s 

# Example 9.3     
 T1   =   -20.0 degC
 P2   =   1.0 Mpa
 T2   =   40.0 degC
 s2   =   1.7148 kJ/kg/degK
 s1   =   1.7395 kJ/kg/degK

# Example 9.4     
 P1   =   100.0 kPa
 T1   =   290.0 degK
 P2   =   1000.0 kPa
 k  =   1.4 
 T2 = T1 * (P2 / P1)^((k-1)/k)  ; degK 
 CP0   =   1.004 kJ/kg/degK
 w = CP0 * (T1 - T2)  ; kJ/kg 


# Example 9.4E     
 P1   =   14.7 psi
 T1   =   520.0 degR
 P2   =   147.0 psi
 k  =   1.4 
 T2 = T1 * (P2 / P1)^((k-1)/k)  ; degR 
 CP0   =   0.240 Btu/lbm/degR
 w = CP0* (T1 - T2)  ; Btu/lbm 


# Example 9.5     
 mdot1  =   2.0 kg/s
 P1   =   300.0 kPa
 T1   =   200.0 degC
 T2   =   20.0 degC
 P2   =   300.0 kPa
 h1   =   2865.5 kJ/kg
 s1   =   7.3115 kJ/kg/degK
 h2   =   83.9 kJ/kg
 s2   =   0.2966 kJ/kg/degK
 h3   =   2725.3 kJ/kg
 s3   =   6.9918 kJ/kg/degK
 mdot2 = mdot1 * (h1 - h3) / (h3 - h2) ; kg/s 
 mdot3 = mdot1 + mdot2  ; kg/s 
 Sdotgen = mdot3 * s3 - mdot1 * s1 - mdot2 * s2  ; kW/degK 


# Example 9.6     
 V0  =   40.0 liter
 P1   =   100.0 kPa
 T1   =   17.0 degC
 P2   =   1000.0 kPa
 sT1   =   6.8352 kJ/kg/degK
 R  =   0.2870 kJ/kg degK
 sT2 = sT1 + R * ln(P2 / P1)  ; kJ/kg/degK 
 T2   =   555.7 degK
 u2   =   401.4900 kJ/kg
 m1 = P1 * V0 / R / T1  ; kg 
 m2 = P2 * V0 / R / T2  ; kg 
 min0 = m2 - m1  ; kg 
 hin   =   290.4 kJ/kg
 u1   =   207.2 kJ/kg
 u2   =   401.5 kJ/kg
 W12 = min0 * hin + m1 * u1 - m2 * u2  ; kJ 


# Example 9.7     
 P1   =   100.0 kPa
 T1   =   30.0 degC
 P2   =   5.0 Mpa
 h  =   0.0010040 m^3/kg
 W1 = h * (P2 - P1)  ; kJ/kg 


# Example 9.8     
 P1   =   1.0 Mpa
 T1   =   300.0 degC
 P2   =   15.0 kPa
 wa   =   600.0 kJ/kg
 h1   =   3051.2 kJ/kg
 s1   =   7.1228 kJ/kg/degK
 sf    =   0.7548 kJ/kg/degK
 sfg   =   7.2536 kJ/kg/degK
 xes = (s1 - sf) / sfg   
 hf    =   225.9000 kJ/kg
 hfg   =   2373.1000 kJ/kg
 hes = hf + xes * hfg  ; kJ/kg 
 ws  = h1 - hes  ; kJ/kg 
 h = wa / ws  
 he = h1 - wa  ; kJ/kg 
 xe = (he - hf) / hfg  

# Example 9.9     
 T1   =   1600.0 degK
 P2   =   100.0 kPa
 T2   =   830.0 degK
 h   =   0.85 
 h1   =   1757.3 kJ/kg
 s1   =   8.6905 kJ/kg/degK
 he   =   855.30 kJ/kg
 w  = h1 - he  ; kJ/kg 
 ws  = w / h  ; kJ/kg 
 hes = h1 - ws  ; kJ/kg 
 Tes   =   683.7 degK
 ses   =   7.7148 kJ/kg/degK
 R  =   0.2870 kJ/kg degK
 P1 = P2 / exp((ses - s1) / R) ; kPa 

# Example 9.10     
 P1   =   100.0 kPa
 T1   =   300.0 degK
 P2   =   150.0 kPa
 h   =   0.70 
 CP0   =   1.0040 kJ/kg/degK
 k  =   0.29 
 T2 = T1 * (P2 / P1)^k  ; degK 
 hes = CP0 * (T1 - T2)  ; kJ/kg 
 ws = hes / h   ; kJ/kg 
 Te = T1 - ws / CP0  ; degK 


# Homework Problems:

TBD