Phase equilibrium of a pure substance¶

[1]:

import si_units as si
import feos

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

[2]:

vle_t = feos.PhaseEquilibrium.pure(eos, 350 * si.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 = feos.PhaseEquilibrium.pure(eos, 15 * si.BAR)
vle_p

[3]:

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