Los cálculos de equilibrio líquido-vapor con ecuaciones de estado son ampliamente utilizados en la industria petrolera y consumen gran cantidad de tiempo en las simulaciones numéricas. Por esta razón, en el presente trabajo se estudian los principales modelos de cálculo flash basados en sustituciones sucesivas y las diferentes partes que conforman estos algoritmos. Los análisis son realizados sobre los diferentes métodos para determinar la fracción de vapor, correlaciones para estimar las constantes de equilibrio y las diferentes propuestas para actualizar las relaciones de equilibrio en el proceso iterativo. Los resultados de la evaluación muestran que el método de Leibovici- Neoschil con la rutina de Newton-Raphson, es la opción más eficiente para determinar la fracción de vapor, ya que requiere un menor número de iteraciones y, por ende, un menor tiempo de cómputo. Para las sensibilidades realizadas con respecto al tipo de correlación utilizada para estimar las relaciones de equilibrio iniciales, se encontró que la propuesta de Wilson es la más efectiva, ya que se obtiene la convergencia de un mayor número de casos. Del análisis de los métodos acelerados de Risnes-Jessen, Mehra y Crowe-Nishio para la actualización de las relaciones de equilibrio se pudo concluir que el algoritmo de Mehra tipo I es el más rápido de todas las opciones, sin llegar a converger en ningún caso a la solución trivial. Por último, se propone un nuevo algoritmo que combina los métodos de Crowe-Nishio, Mehra tipo I y Leibovici, y con el cual se obtiene convergencia de un 100 % de los casos y el menor tiempo de cómputo con respecto a las metodologías bajo estudio.

The calculations of vapor-liquid equilibrium with equations of state are widely used in the oil industry, and they consume a lot time in numerical simulations. Therefore, in this paper several flash calculation models based on successive substitutions and the different parts that defines these algorithms are studied. Basically, the evaluations are performed with different methods to determine the fraction of vapor, correlations to estimate the equilibrium constants and the various proposals for upgrading equilibrium ratios in the iterative process. The results of the evaluation show that the method of Leibovici-Neoschil with Newton-Raphson routine is the most efficient way to determine the fraction of vapor since it requires fewer iterations and time. Sensitivity analyses are performed considering the correlation used to estimate initial equilibrium relationships. It has been found that Wilson’s proposal is the most effective since the convergence has been achieved in a greater number of cases. The analysis of the accelerated methods of Risnes-Jessen, Mehra et al. and Crowe Nishio of updating the equilibrium relationships could conclude that the type I algorithm of Mehra is the fastest of all the options, without having to converge in any case to the trivial solution. Finally, a new algorithm which combines the methods from Crowe Nishio, Mehra type I and Leibovici has been proposed. With this model the convergence is achieved in 100% of the cases. Additionally, the computation time is lower than any of the other methods under study.