Surface tension 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'
)
functional = feos.HelmholtzEnergyFunctional.pcsaft(parameters)

[2]:

vle = feos.PhaseEquilibrium.pure(functional, 350 * si.KELVIN)
vle

[2]:

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

[3]:

surface_tension_pdgt = (feos.PlanarInterface
    .from_pdgt(vle, 1024)
    .solve()
    .surface_tension)
surface_tension_pdgt

[3]:

$20.191\,\mathrm{\frac{mN}{m}}$

[4]:

surface_tension_tanh = (feos.PlanarInterface
    .from_tanh(vle, 1024, 100 * si.ANGSTROM, 500 * si.KELVIN)
    .solve()
    .surface_tension
)
surface_tension_tanh

[4]:

$20.191\,\mathrm{\frac{mN}{m}}$