Phase equilibrium of a pure substance¶

[1]:

from feos.si import BAR, KELVIN
from feos.eos import EquationOfState, PhaseEquilibrium
from feos.pcsaft import PcSaftParameters

parameters = PcSaftParameters.from_json(
    substances=['methanol'],
    pure_path='../parameters/pcsaft/gross2002.json'
)
eos = EquationOfState.pcsaft(parameters)

[2]:

vle_t = PhaseEquilibrium.pure(eos, 350 * KELVIN)
vle_t

[2]:

	temperature	density
phase 1	350.00000 K	62.68366 mol/m³
phase 2	350.00000 K	23.13883 kmol/m³

[3]:

vle_p = PhaseEquilibrium.pure(eos, 15 * BAR)
vle_p

[3]:

	temperature	density
phase 1	426.13707 K	575.28700 mol/m³
phase 2	426.13707 K	20.09406 kmol/m³