## Peng-Robinson superancillaries

EoSSuperancillaries.pr_vsatFunction
vl,vv = pr_vsat(T,a,b)

Returns the saturation liquid and vapour volumes of the Peng-Robinson 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.
EoSSuperancillaries.pr_vlsatFunction
vl = pr_vlsat(T,a,b)

Returns the saturation liquid volume of the Peng-Robinson 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.
EoSSuperancillaries.pr_vvsatFunction
vv = pr_vvsat(T,a,b)

Returns the saturation vapour volume of the Peng-Robinson 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.
EoSSuperancillaries.pr_psatFunction
p = pr_psat(T,a,b)

Returns the saturation pressure of the Peng-Robinson 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.