A generalized model for complex distillation processes

which taking into account both multiphase (two liquid phases) and reaction on the trays is presented in this paper. The algorithm is worked out with globally convergent Newton homotopy continuation method for solving highly nonlinear systems. The algorithm for the calculation of liquid -liquid equilibria is featured by convergence

rapidity and reliability.