Three phase flash calculations using free energy minimisation

Abstract The distillation of systems exhibiting partial miscibility over part of the composition range is growing in interest and the modeling of three phase distillation columns is highly desirable. This paper presents work which has been produced as part of a project directed towards such objective. An algorithm for flash calculations in which partially miscible systems are involved is described. The algorithm is based on the principle of minimisation of the system Gibbs free energy which is becoming increasingly popular allowing as it does for the simultaneous satisfaction of the constraints of equilibria and mass balance requirements. The minimisation procedure is based on a Newton-Raphson algorithm and so in order to avoid trivial solutions it is necessary to have an adequate choice of initial conditions. In the present paper this is taken into due account and a policy for selecting initial compositions and flowrates is presented. Numerical examples are presented for the ternary systems ethanol/ethylacetate/water and methanol/ethylacetate/water.