This function allows calculating the compressibility factor, the fugacity coefficient and density of a pure compound for both liquid and vapour phases by means of the Peng-Robinson equation of state (PR EOS).
Temperature, pressure, molar weigth, and supercritical properties of the selected compound are needed to be specified as parameters.
It is a Beta version, full version will include multicomponent systems and other equations of state.
This MATLAB function provides an educational tool to predict fluid phase equilibrium and can be easily extended to a multicomponent system.
Requirements:
· MATLAB Release: R13
Download