pcsaft.pcsaft_hres

pcsaft.pcsaft_hres()

Calculate the residual enthalpy for one phase of the system.

Parameters:
  • t (float) – Temperature (K)
  • rho (float) – Molar density (mol m-3)
  • x (ndarray, shape (n,)) – Mole fractions of each component. It has a length of n, where n is the number of components in the system.
  • params (dict) –

    A dictionary containing PC-SAFT parameters that can be passed for use in PC-SAFT:

    m : ndarray, shape (n,)
    Segment number for each component.
    s : ndarray, shape (n,)
    Segment diameter for each component. For ions this is the diameter of the hydrated ion. Units of Angstrom.
    e : ndarray, shape (n,)
    Dispersion energy of each component. For ions this is the dispersion energy of the hydrated ion. Units of K.
    k_ij : ndarray, shape (n,n)
    Binary interaction parameters between components in the mixture. (dimensions: ncomp x ncomp)
    e_assoc : ndarray, shape (n,)
    Association energy of the associating components. For non associating compounds this is set to 0. Units of K.
    vol_a : ndarray, shape (n,)
    Effective association volume of the associating components. For non associating compounds this is set to 0.
    dipm : ndarray, shape (n,)
    Dipole moment of the polar components. For components where the dipole term is not used this is set to 0. Units of Debye.
    dip_num : ndarray, shape (n,)
    The effective number of dipole functional groups on each component molecule. Generally this is set to 1, but some implementations use this as an adjustable parameter that is fit to data.
    z : ndarray, shape (n,)
    Charge number of the ions
    dielc : float
    Dielectric constant of the medium to be used for electrolyte calculations.
    assoc_scheme : list, shape (n,)
    The types of association sites for each component. Use None for molecules without association sites. If a molecule has multiple association sites, use a nested list for that component to specify the association scheme for each site. The accepted association schemes are those given by Huang and Radosz (1990): 1, 2A, 2B, 3A, 3B, 4A, 4B, 4C. If e_assoc and vol_a are given but assoc_scheme is not, the 2B association scheme is assumed (which would, for example, correspond to one hydroxyl functional group).
Returns:

hres – Residual enthalpy (J mol-1)

Return type:

float