ex16.m

function ex16()
% Clausius-Rankine
addpath('/pfad/zu/water95');  

% gegebene Werte
T1 = celsiusToKelvin(500);
p1 = 60e5;
p2 = 0.1e5;

% Punkt 3 im Nassdampfbereich, ggb. p3 = p2, x3 = 0
p3 = p2;
T3 = saturationTemperature(p3);
[~, rhop3, rhopp3] = saturationPressure(T3);
rho3 = rhop3;      % x3 = 0
h3 = enthalpy(rho3, T3)
s3 = entropy(rho3, T3)

% Punkt 4 im Wasserbereich, ggb. p4 = p1, s4 = s3
p4 = p1;
s4 = s3;
T4 = temperaturePS(p4, s4);
rho4 = density(p4, T4);
h4 = enthalpy(rho4, T4)

% Punkt 1 im Heißdampfbereich, ggb. p1, T1
rho1 = density(p1, T1);
h1 = enthalpy(rho1, T1)
s1 = entropy(rho1, T1)

% Punkt 2 im Nassdampfbereich, ggb. p2, s2 = s1
s2 = s1;
rhop2 = rhop3;
rhopp2 = rhopp3;
T2 = T3;
sp2 = entropy(rhop2, T2)
spp2 = entropy(rhopp2, T2)
x2 = (s2 - sp2)/(spp2 - sp2)

hp2 = enthalpy(rhop2, T2)
hpp2 = enthalpy(rhopp2, T2)
h2 = hp2 + x2*(hpp2 - hp2)

eta = 1 + (h3 - h2)/(h1 - h4)
rW = 1 - (h4 - h3)/(h1 - h2)