This paper presents a relatively simple procedure for determining azeotrope coordinates for an isothermal or isobaric binary system, at low pressure and with a single or double homogeneous azeotrope, based on an expression for g E . Previous investigations, some more comprehensive, have used rather complex procedures for binary and multicomponent systems. have treated fi nding the temperatures and compositions of all azeotropes for isobaric multicomponent systems as a root-fi nding problem solved by a homotopy method. Similarly, such systems have been treated by as a global optimization problem. Both groups used the Wilson (1964) equation for solution of examples, but comparisons were not made with experimental results. For an isobaric fi vecomponent system used by both groups, the two methods agreed in the number and nature of azeotropes, but differed somewhat in values of the coordinates, particularly for compositions of the A relatively simple procedure is described to calculate the coordinates, temperature (T az ), pressure (P az ), and composition (x 1az ), of an azeotrope (for either an isothermal or isobaric binary system at low pressure and with a single or double homogeneous azeotrope) from an expression for the liquid-phase excess molar Gibbs function (g E ). General results are based on (∂g E /∂x 1 ) az = RT az ln[p 2 *(T az )/p 1 *(T az )], and P az = (γ 1 p 1 *) az = (γ 2 p 2 *) az , where R is the gas constant, p* is saturation vapour pressure, and γ is liquid-phase activity coeffi cient. Specifi c results are given for the Redlich-Kister, van Laar, Wilson, and NRTL equations. Numerical examples are provided for both an isothermal and an isobaric system. The procedure provides a means to obtain azeotrope coordinates that are consistent with a g E expression obtained either from experimental data or from a model. It is applicable to polyazeotropy, whereas the criteria given previously are generally not applicable to polyazeotropy.
On décrit une procédure relativement simple pour calculer les coordonnées, la température (T az ), la pression (P az ) et la composition (x 1az ) d'un azéotrope (pour un système isotherme ou isobare binaire à faible pression et avec un azéotrope simple ou double homogène) à partir d'une expression pour la fonction de Gibbs d'excès molaire de la phase liquide (g E ). Les résultats généraux reposent sur (∂g E /∂x 1 ) az = RT az ln[p 2 *(T az )/ p 1 *(T az )], et P az = (γ 1 p 1 *) az = (γ 2 p 2 *) az , où R est la constante des gaz, p* la pression de vapeur saturante et γ le coeffi cient d'activité de la phase liquide. Des résultats spécifi ques sont donnés pour les équations de Redlich-Kister, van Laar, Wilson et NRTL. Des exemples numériques sont présentés à la fois pour un système isotherme et un système isobare. Cette méthode offre un moyen d'obtenir des paramètres d'azéotrope qui soient consistants avec une expression de gE obtenue à partir de données expérimentales ou d'un modèle. Elle est applicable à la polyazéotropie, ce qui n'est généralement pas le cas des critères donnés précédemment.