Modelling the Vapor–Liquid Equilibrium of Polymer Solutions Using a Cubic Equation of State

This work presents the vapor–liquid equilibrium calculations in an isothermal flash, applied to polymer solutions, using the Peng–Robinson cubic equation of state modified by Stryjek–Vera, and the mixing rule introduced by Wong–Sandler. This rule allows combining the rigid lattice thermodynamic model of Flory–Huggins to the Peng–Robinson–Stryjek–Vera equation of state. As the Gibbs free energy must be minimum in the equilibrium state, a stochastic optimization method, the simulated annealing algorithm, was used to find out the extreme of this thermodynamic potential.