The VLLE flash is important in water and hydrocarbons mixtures, hydrocarbon and CO 2 rich mixtures, and hydrocarbon methane rich mixtures that are encountered in reservoir performance and recovery studies. A robust VLLE flash algorithm is proposed. The equilibrium and mass balance equations are solved as a constrained minimization problem. An inverse barrier function is used to handle the inequality constrains to solve for the phase fractions. It warrants always arriving to the solution. The challenging cases analyzed showed that the initialization procedure proposed, together with successive substitution iteration in the outer loop, is a good method for a stable VLLE flash algorithm, even near critical points. Whenever the result is in the region outside the three-phase physical domain, the solution suggests that the system has fewer phases. In one of the cases analyzed, a region with three liquid phases was encountered and the algorithm found two different solutions with positive phase fractions. V C 2015 American Institute of Chemical Engineers AIChE J, 61: [3081][3082][3083][3084][3085][3086][3087][3088][3089][3090][3091][3092][3093] 2015

No abstract