>
restart;
Použijeme zadanou Antoineovu rovici
>
rov:=log[10](p)=5.92617-1101.47/(t+215.15);
jejímž řešením dostáváme normální teplotu varu
>
tnbv:=solve(subs(p=101.325,rov),t);
>
Tnbv:=tnbv+273.15;
Vyjadřime derivaci logaritmu tlaku podle teploty
>
p:=solve(rov,p);
>
t:=T-273.15;
>
der:=factor(diff(ln(p),T));
a) Tuto derivaci dosadíme do Clausiovy-Clapeyronovy rovnice (s aproximaci idealního plynu)
>
rov1:=der=dH/(R*T^2);
a vypočteme výparné teplo za teploty normálního bodu varu
>
R:=8.314;
>
solve(subs(T=Tnbv,rov1),dH);
b) Nejprve vypočteme molární objem kapaliny a její kompresibilitní faktor
>
M:=0.072;
>
rhol:=800;
>
Vl:=M/rhol;
>
p:=101325;
>
zl:=p*Vl/(R*Tnbv);
Dosadíme kompresibilitní faktor plynu
>
zg:=0.95;
Dříve vypočtenou derivaci dosadíme do Clausiovy-Clapeyronovy rovnice (se započtením neideálního
chování plynu a objemu kapaliny)
>
rov2:=der=dH/((zg-zl)*R*T^2);
a opět vypočteme výparné teplo za teploty normálního bodu varu
>
solve(subs(T=Tnbv,rov2),dH);
>