2.1.4 Vyjádření izotermní kompresibility z tlakového viriálního rozvoje

Ze stavové rovnice z = 1+Bp/(RT), kde B = B(T), odvoďte vztah pro koeficient izotermní stlačitelnosti $..$ .

↓

κ_T = RT/[R T p+B p²].

Nejprve definujeme použitou stavovou rovnici

> z:=1+B*p/(R*T);

a vyjádříme z ní virální koeficient

> Vm:=solve(z=p*Vm/(R*T),Vm);

Derivací pak zjistíme vztah pro stačitelnost

> kappa:=-diff(ln(Vm),p);

Tento vztah pak můžeme dále zjednodušit

> kappa:=factor(kappa);

>