vl,vv = vdw_vsat(T,a,b)

Returns the saturation liquid and vapour volumes of the vdW cubic equation of state.

Inputs:

• T: Temperature (Kelvin)
• a: Atraction parameter [m^6*Pa/mol]
• b: Covolume [m^3/mol]

Outputs:

• vl : Saturation Liquid Volume [m^3]. Returns NaN if the value is outside the range of the ancillary.
• vv : Saturation Vapour Volume [m^3]. Returns NaN if the value is outside the range of the ancillary.

vl = vdw_vlsat(T,a,b)

Returns the saturation liquid volume of the vdW cubic equation of state.

Inputs:

• T: Temperature (Kelvin)
• a: Atraction parameter [m^6*Pa/mol]
• b: Covolume [m^3/mol]

Outputs:

• vl : Saturation Liquid Volume [m^3]. Returns NaN if the value is outside the range of the ancillary.

vv = vdw_vvsat(T,a,b)

Returns the saturation vapour volume of the vdW cubic equation of state.

Inputs:

• T: Temperature (Kelvin)
• a: Atraction parameter [m^6*Pa/mol]
• b: Covolume [m^3/mol]

Outputs:

• vv : Saturation Vapour Volume [m^3]. Returns NaN if the value is outside the range of the ancillary.

p = vdw_psat(T,a,b)

Returns the saturation pressure of the vdW cubic equation of state.

Inputs:

• T: Temperature (Kelvin)
• a: Atraction parameter [m^6*Pa/mol]
• b: Covolume [m^3/mol]

Outputs:

• p : Saturation pressure [Pa]. Returns NaN if the value is outside the range of the ancillary.