Surface tension of a pure substance¶

[1]:

from feos.si import *
from feos.dft import HelmholtzEnergyFunctional, PhaseEquilibrium, PlanarInterface
from feos.pcsaft import PcSaftParameters

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

[2]:

vle = PhaseEquilibrium.pure(functional, 350 * 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 = (PlanarInterface
    .from_pdgt(vle, 1024)
    .solve()
    .surface_tension)
surface_tension_pdgt

[3]:

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

[4]:

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

[4]:

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