# Loi de Raoult

En physique, et plus particulièrement en thermodynamique, la loi de Raoult énonce que :
   « Dans une solution idéale, à température constante, la pression partielle en phase vapeur d'un constituant est proportionnelle à sa fraction molaire en phase liquide. »
Cette loi a été établie empiriquement par le physicien français François-Marie Raoult en 1882, elle est dérivée de sa loi de la tonométrie. Elle est utilisée dans de nombreux domaines de la chimie, de la physique et de la météorologie. Elle permet notamment de calculer les équilibres liquide-vapeur des solutions idéales liquides dont la phase vapeur est un mélange de gaz parfaits.

## Énoncé et démonstration
### Énoncé mathématique

Étant donné que pour le corps pur la pression est égale à sa pression de vapeur saturante, le coefficient de proportionnalité est égal à cette pression à la température du mélange. La forme mathématique de la loi de Raoult s'écrit, pour une espèce chimique $i$ :
* Loi de Raoult pression partielle de $i$ : $P_{i}=x_{i}^{\text{g}}P=x_{i}^{\text{l}}P_{i}^{\text{sat}}$
avec les notations :
* $P$ la pression totale du mélange ;
* $P_{i}$ la pression partielle du composant $i$, par définition $P_{i}=x_{i}^{\text{g}}P$ ;
* $P_{i}^{\text{sat}}$ la pression de vapeur saturante du composant $i$ pur à la température $T$ du mélange ;
* $x_{i}^{\text{g}}$ la fraction molaire du constituant $i$ dans la phase vapeur ;
* $x_{i}^{\text{l}}$ la fraction molaire du constituant $i$ dans la phase liquide.
Pour un corps pur, $x_{i}^{\text{g}}=x_{i}^{\text{l}}=1$ et par conséquent $P_{i}=P=P_{i}^{\text{sat}}$.
Dans un mélange, la relation de Duhem-Margules impose que si l'un des composants suit la loi de Raoult alors les autres composants suivent également la loi de Raoult, ou, pour les gaz dissouts, la loi de Henry.

### Démonstration

La loi de Raoult a été établie empiriquement par son auteur, sur base de sa loi de la tonométrie pour les liquides et de la loi de Dalton pour les gaz. La démonstration rigoureuse qui suit est postérieure au développement du formalisme mathématique de la thermodynamique par Willard Gibbs, qui introduisit notamment la notion de potentiel chimique, et par Gilbert Lewis, qui introduisit les notions de fugacité et d'activité chimique.
Dans une solution à l'équilibre, pour chaque constituant $i$ les potentiels chimiques dans la phase liquide $\mu _{i}^{\text{l}}$ et dans la phase vapeur $\mu _{i}^{\text{g}}$ sont égaux :
    (1) $\mu _{i}^{\text{l}}\!\left(P,T,x^{\text{l}}\right)=\mu _{i}^{\text{g}}\!\left(P,T,x^{\text{g}}\right)$
Pour la phase liquide et la phase vapeur, considérées comme des solutions idéales, les potentiels chimiques sont développés selon :
    (2) $\mu _{i}^{\text{l}}\!\left(P,T,x^{\text{l}}\right)=\mu _{i}^{\text{l,*}}\!\left(P,T\right)+RT\ln x_{i}^{\text{l}}$
    (3) $\mu _{i}^{\text{g}}\!\left(P,T,x^{\text{g}}\right)=\mu _{i}^{\text{g,*}}\!\left(P,T\right)+RT\ln x_{i}^{\text{g}}$
avec :
* $P$ la pression totale du mélange ;
* $R$ la constante universelle des gaz parfaits ;
* $T$ la température.
La variation isotherme du potentiel chimique en fonction de la pression permet d'écrire pour le composant $i$ pur en phase liquide :
    $\mu _{i}^{\text{l,*}}\!\left(P,T\right)=\mu _{i}^{\text{l,*}}\!\left(P_{i}^{\text{sat}},T\right)+\int _{P_{i}^{\text{sat}}}^{P}{\bar {V}}_{i}^{\text{l,*}}\!\left(P,T\right)\,\mathrm {d} P$
avec ${\bar {V}}_{i}^{\text{l,*}}$ le volume molaire du composant $i$ à l'état de liquide pur. Le dernier terme est généralement négligé, soit :
    $\mu _{i}^{\text{l,*}}\!\left(P,T\right)\approx \mu _{i}^{\text{l,*}}\!\left(P_{i}^{\text{sat}},T\right)$
Ceci traduit le fait que les liquides sont généralement peu sensibles à la pression et sont considérés comme incompressibles : leurs propriétés sont considérées, à température donnée, comme étant celles du liquide à saturation à cette même température, quelle que soit la pression. En introduisant cette relation dans (2), il vient :
    (4) $\mu _{i}^{\text{l}}\!\left(P,T,x^{\text{l}}\right)\approx \mu _{i}^{\text{l,*}}\!\left(P_{i}^{\text{sat}},T\right)+RT\ln x_{i}^{\text{l}}$
La variation isotherme du potentiel chimique en fonction de la pression permet d'écrire pour le composant $i$ pur en phase vapeur :
    $\mu _{i}^{\text{g,*}}\!\left(P,T\right)=\mu _{i}^{\text{g,*}}\!\left(P_{i}^{\text{sat}},T\right)+\int _{P_{i}^{\text{sat}}}^{P}{\bar {V}}_{i}^{\text{g,*}}\!\left(P,T\right)\,\mathrm {d} P$
avec ${\bar {V}}_{i}^{\text{g,*}}$ le volume molaire du composant $i$ à l'état de gaz pur. Si ce gaz est considéré comme parfait alors :
    ${\bar {V}}_{i}^{\text{g,*}}\!\left(P,T\right)={\frac {RT}{P}}$
    $\mu _{i}^{\text{g,*}}\!\left(P,T\right)=\mu _{i}^{\text{g,*}}\!\left(P_{i}^{\text{sat}},T\right)+RT\int _{P_{i}^{\text{sat}}}^{P}{\mathrm {d} P \over P}=\mu _{i}^{\text{g,*}}\!\left(P_{i}^{\text{sat}},T\right)+RT\ln \!\left({\frac {P}{P_{i}^{\text{sat}}}}\right)$
En introduisant cette relation dans (3), il vient :
    (5) $\mu _{i}^{\text{g}}\!\left(P,T,x^{\text{g}}\right)=\mu _{i}^{\text{g,*}}\!\left(P_{i}^{\text{sat}},T\right)+RT\ln \!\left({\frac {P}{P_{i}^{\text{sat}}}}\right)+RT\ln x_{i}^{\text{g}}$
Les relations (4) et (5) permettent donc de réécrire la relation (1) selon :
    (6) $\mu _{i}^{\text{l,*}}\!\left(P_{i}^{\text{sat}},T\right)+RT\ln x_{i}^{\text{l}}=\mu _{i}^{\text{g,*}}\!\left(P_{i}^{\text{sat}},T\right)+RT\ln \!\left({\frac {P}{P_{i}^{\text{sat}}}}\right)+RT\ln x_{i}^{\text{g}}$
La relation (1) étant aussi valable pour un corps pur à l'équilibre liquide-vapeur (saturation), celui-ci à la température $T$ étant à sa pression de vapeur saturante $P_{i}^{\text{sat}}$, il vient :
    $\mu _{i}^{\text{l,*}}\!\left(P_{i}^{\text{sat}},T\right)=\mu _{i}^{\text{g,*}}\!\left(P_{i}^{\text{sat}},T\right)$
La relation (6) produit donc :
    $x_{i}^{\text{l}}P_{i}^{\text{sat}}=x_{i}^{\text{g}}P=P_{i}$

### Correction de Poynting

Si l'on ne néglige pas le terme $\int _{P_{i}^{\text{sat}}}^{P}{\bar {V}}_{i}^{\text{l,*}}\!\left(P,T\right)\,\mathrm {d} P$ dans l'expression du potentiel chimique en phase liquide, celui-ci dans la relation (4) s'exprime exactement selon :
    $\mu _{i}^{\text{l}}\!\left(P,T,x^{\text{l}}\right)=\mu _{i}^{\text{l,*}}\!\left(P_{i}^{\text{sat}},T\right)+\int _{P_{i}^{\text{sat}}}^{P}{\bar {V}}_{i}^{\text{l,*}}\!\left(P,T\right)\,\mathrm {d} P+RT\ln x_{i}^{\text{l}}$
    $\mu _{i}^{\text{l}}\!\left(P,T,x^{\text{l}}\right)=\mu _{i}^{\text{l,*}}\!\left(P_{i}^{\text{sat}},T\right)+RT\ln \!\left(x_{i}^{\text{l}}\exp \!\left({\frac {\int _{P_{i}^{\text{sat}}}^{P}{\bar {V}}_{i}^{\text{l,*}}\!\left(P,T\right)\,\mathrm {d} P}{RT}}\right)\right)$
avec ${\bar {V}}_{i}^{\text{l,*}}\!\left(P,T\right)$ le volume molaire du composant $i$ liquide pur à $P$ et $T$. Le facteur exponentiel est appelé correction de Poynting ou facteur de Poynting :
${\mathcal {P}}_{i}\!\left(P,T\right)=\exp \!\left({\frac {\int _{P_{i}^{\text{sat}}}^{P}{\bar {V}}_{i}^{\text{l,*}}\!\left(P,T\right)\,\mathrm {d} P}{RT}}\right)$
La loi de Raoult avec correction de Poynting devient :
* Loi de Raoult et correction de Poynting pression partielle de $i$ : $P_{i}=x_{i}^{\text{g}}P=x_{i}^{\text{l}}P_{i}^{\text{sat}}{\mathcal {P}}_{i}$
Le liquide pur $i$ n'existe de façon stable à la température $T$ que pour des pressions $P\geq P_{i}^{\text{sat}}\!\left(T\right)$. Pour des pressions $P<P_{i}^{\text{sat}}\!\left(T\right)$ le liquide $i$ pur soit n'existe pas, soit existe sous une forme thermodynamiquement instable : on peut alors ne pas trouver de valeur pour ${\bar {V}}_{i}^{\text{l,*}}\!\left(P,T\right)$. Toutefois, les liquides étant peu compressibles, le volume molaire liquide ${\bar {V}}_{i}^{\text{l,*}}$ peut être considéré comme indépendant de la pression ; il est souvent donné à saturation du liquide : ${\bar {V}}_{i}^{\text{l,*}}\!\left(P,T\right)\approx {\bar {V}}_{i}^{\text{l,*}}\!\left(P_{i}^{\text{sat}}\!\left(T\right),T\right)\approx {\bar {V}}_{i}^{\text{l,*,sat}}\!\left(T\right)$, auquel cas :
    ${\mathcal {P}}_{i}\!\left(P,T\right)\approx \exp \!\left({\frac {{\bar {V}}_{i}^{\text{l,*,sat}}\!\left(T\right)\cdot \left(P-P_{i}^{\text{sat}}\right)}{RT}}\right)$
Ceci permet de calculer la correction de Poynting à toute pression, y compris celles pour lesquelles, à température donnée, le corps $i$ pur n'existe pas sous forme de liquide stable.
Aux basses pressions (moins de 10 atm, domaine d'application de la loi des gaz parfaits), la correction de Poynting est négligeable :
    ${\mathcal {P}}_{i}\!\left(P,T\right)\approx 1$
Pour un mélange, la relation de Duhem-Margules induit que si l'un des composants suit la loi de Raoult avec correction de Poynting alors les autres composants suivent également la loi de Raoult avec correction de Poynting ou la loi de Henry avec correction de Poynting.
Avec la loi de Raoult sans correction de Poynting, connaissant pour un mélange liquide sa composition (fractions molaires $x_{i}^{\text{l}}$) et la pression de vapeur saturante de chacun des composants $P_{i}^{\text{sat}}$ à la température $T$, la pression se calcule explicitement selon :
    $P=\sum _{i=1}^{N}x_{i}^{\text{l}}P_{i}^{\text{sat}}$
La loi de Raoult avec correction de Poynting nécessite un calcul itératif, puisque l'expression est implicite en $P$ :
    $P=\sum _{i=1}^{N}x_{i}^{\text{l}}P_{i}^{\text{sat}}{\mathcal {P}}_{i}\!\left(P,T\right)$

## Limites et extension de la loi de Raoult
### Limites de la loi idéale
#### Écarts des équilibres liquide-vapeur réels

La loi de Raoult s’applique à des systèmes idéaux dans lesquels les interactions entre molécules de même espèce ou d'espèces différentes sont toutes identiques. Les écarts à la loi de Raoult observés dans les équilibres liquide-vapeur réels peuvent être positifs ou négatifs selon l'attraction électronique des composants du mélange.

##### Écart positif

    Lorsque dans un mélange composé de deux constituants A et B les molécules du même constituant (interactions A-A et B-B) s'attirent plus fortement que les molécules d'espèces différentes entre elles (interaction A-B), l’ajout par exemple de molécules B dans le liquide A pur a pour effet de diminuer les forces d’attraction entre les molécules A et la volatilité de A est augmentée. Il y a alors une déviation positive par rapport à la loi de Raoult : la pression totale de la vapeur en équilibre au-dessus de la solution liquide est supérieure à celle calculée.

##### Écart négatif

    Si, par contre, les molécules de A et de B sont plus attirées entre elles (interaction A-B) que par celles du même constituant (interactions A-A et B-B), l'ajout de molécules B au liquide A pur augmente les forces d’attraction auxquelles sont soumises les molécules A et celles-ci deviennent moins volatiles. Il y a alors une déviation négative par rapport à la loi de Raoult : la pression totale de la vapeur en équilibre au-dessus de la solution liquide est inférieure à celle calculée.
Le plus souvent les pressions d'équilibre d'un mélange binaire réel restent comprises, pour une température donnée, entre les pressions de vapeur saturante des deux composants à cette même température.

#### Les azéotropes

Lorsque le comportement non idéal d'un mélange binaire devient suffisamment important, les pressions d'équilibre du mélange peuvent être soit plus grandes que la plus grande des pressions de vapeur saturante des composants, soit plus petites que la plus petite des pressions de vapeur saturante des composants, ce qui est représenté sur la figure ci-dessus. Dans les deux cas, les courbes d'équilibre liquide-vapeur présentent nécessairement un extrémum (respectivement un maximum et un minimum) qui correspond, en vertu du théorème de Gibbs-Konovalov, à un azéotrope. L'azéotrope est dit positif si sa pression est supérieure à la pression prévue par la loi de Raoult et négatif dans le cas contraire.
La loi de Raoult est fondamentalement incapable de représenter des azéotropes, c'est-à-dire des mélanges pour lesquels, dans certaines conditions de pression, température et composition, le liquide et la vapeur en équilibre ont la même composition. Dans ces conditions la pression d'ébullition et la pression de rosée (voir plus bas) sont égales : un mélange azéotropique se comporte donc comme un corps pur, néanmoins il y a bien plusieurs espèces chimiques en présence. Par exemple, le mélange eau - éthanol contenant 96 % d'éthanol présente un azéotrope à 78,1 °C sous une pression de 1 atm : la loi de Raoult ne peut pas représenter ce mélange.
Un azéotrope est une preuve de non-idéalité du mélange. Dans le cas du mélange eau - éthanol, cette non-idéalité se traduit également par le fait qu'un mélange d'un litre d'eau et d'un litre d'éthanol a un volume de 1,92 litre : les forces d'attraction entre molécules des deux espèces produisent une contraction du volume. La loi de Raoult étant basée sur l'hypothèse de l'idéalité des phases liquide et vapeur, ce type de comportement est donc en dehors de son champ d'application.
Soit un mélange binaire contenant les espèces 1 et 2. On note, à une température $T$ donnée :
* $x_{1}^{\text{g}}$ et $x_{2}^{\text{g}}$ les fractions molaires en phase vapeur ($x_{1}^{\text{g}}+x_{2}^{\text{g}}=1$) ;
* $x_{1}^{\text{l}}$ et $x_{2}^{\text{l}}$ les fractions molaires en phase liquide ($x_{1}^{\text{l}}+x_{2}^{\text{l}}=1$) ;
* $P_{1}^{\text{sat}}$ et $P_{2}^{\text{sat}}$ les pressions de vapeur saturante à la température $T$.
Selon la loi de Raoult, l'équilibre liquide-vapeur est donné par les deux équations :
    $x_{1}^{\text{g}}P=x_{1}^{\text{l}}P_{1}^{\text{sat}}$
    $x_{2}^{\text{g}}P=x_{2}^{\text{l}}P_{2}^{\text{sat}}$
À la température $T$ donnée, on cherche la pression et la composition d'un éventuel azéotrope, point auquel les deux phases ont la même composition, soit $x_{1}=x_{1}^{\text{g}}=x_{1}^{\text{l}}$ et $x_{2}=x_{2}^{\text{g}}=x_{2}^{\text{l}}$. Les équations précédentes donnent :
    $x_{1}P=x_{1}P_{1}^{\text{sat}}$
    $x_{2}P=x_{2}P_{2}^{\text{sat}}$
Ces équations admettent deux solutions triviales :
* $P=P_{1}^{\text{sat}}$, avec $x_{1}=1$ et $x_{2}=0$, le corps 1 est pur ;
* $P=P_{2}^{\text{sat}}$, avec $x_{1}=0$ et $x_{2}=1$, le corps 2 est pur.
Il n'existe pas de solution avec $x_{1}\neq 0$ et $x_{2}\neq 0$ car on obtient $P=P_{1}^{\text{sat}}=P_{2}^{\text{sat}}$, or ces deux pressions ne sont égales qu'à un éventuel point de Bancroft. La loi de Raoult est donc incapable de représenter un azéotrope.

#### Cas des gaz dissouts, loi de Henry

La loi de Raoult s'applique à des composants qui présentent une pression de vapeur saturante $P_{i}^{\text{sat}}$ à la température $T$ du mélange. Pour un gaz dissout, c'est-à-dire un fluide supercritique dont la température critique est inférieure à la température du mélange, soit $T_{{\text{c}},i}<T$, la loi de Raoult n'est pas applicable, puisque le gaz ne possède pas de pression de vapeur saturante à cette température. Pour les gaz dissouts, le pendant de la loi de Raoult est la loi de Henry.
La loi de Henry établit la relation suivante entre la pression partielle $P_{i}$ d'un constituant gazeux et sa fraction molaire $x_{i}^{\text{l}}$ dans un solvant liquide $s$ :
$P_{i}=x_{i}^{\text{l}}k_{{\text{H}},i,s}$
où $k_{{\text{H}},i,s}$ est appelé constante de Henry, bien qu'elle dépende de la température et de la pression. La constante de Henry est spécifique du gaz $i$ dans le solvant $s$, c'est-à-dire qu'elle doit être déterminée pour chaque couple « gaz $i$ - solvant $s$ » et n'est pas valable si le gaz $i$ est considéré avec un solvant autre que le solvant $s$ pour lequel elle a été déterminée.

### Extension aux équilibres liquide-vapeur réels

La loi de Raoult n'est valable que si le liquide et sa vapeur sont tous deux des solutions idéales, la phase gaz étant de surcroît un mélange de gaz parfaits. L'application de la loi de Raoult est donc restreinte à des cas particuliers, pour certains mélanges de produits et certaines plages de température et de pression. Elle constitue néanmoins une base pour calculer les équilibres liquide-vapeur des mélanges réels (y compris les azéotropes) à l'aide de facteurs correctifs tels que les coefficients de fugacité pour la phase gaz et les coefficients d'activité pour la phase liquide :
$x_{i}^{\text{g}}\phi _{i}^{\text{g}}P=x_{i}^{\text{l}}\gamma _{i}^{\text{l}}\phi _{i}^{\text{g,*,sat}}P_{i}^{\text{sat}}{\mathcal {P}}_{i}$
avec :
* $\phi _{i}^{\text{g}}$ le coefficient de fugacité du corps $i$ en phase gaz, à $P$, $T$ et composition du mélange gazeux ;
* $\phi _{i}^{\text{g,*,sat}}$ le coefficient de fugacité du corps $i$ pur à saturation en phase gaz, à $P_{i}^{\text{sat}}$ et $T$ ;
* $\gamma _{i}^{\text{l}}$ le coefficient d'activité du corps $i$ en phase liquide, à $P$, $T$ et composition du mélange liquide ;
* ${\mathcal {P}}_{i}$ la correction de Poynting.

## Applications
### Calcul d'un point d'ébullition
#### Pression d'ébullition

La pression d'ébullition d'un mélange, à température donnée, est la pression au-delà de laquelle (pour des pressions plus fortes) le mélange est entièrement liquide : à la pression d'ébullition apparaît la première bulle de vapeur. En deçà de la pression d'ébullition (pour des pressions plus faibles), le mélange présente un équilibre liquide-vapeur, ou est entièrement gazeux (en deçà de la pression de rosée). Soit un mélange liquide de $N$ composants, chaque composant $i$ étant représenté par sa fraction molaire $x_{i}^{\text{l}}$, à la température $T$.
Nous connaissons la température $T$ du mélange, nous pouvons donc calculer la pression de vapeur saturante $P_{i}^{\text{sat}}$ de chacun des composants. La composition de la phase liquide est connue. On suppose que tous les composants du mélange suivent la loi de Raoult. La loi de Raoult sans correction de Poynting donne, en sommant sur l'ensemble des constituants :
    $\sum _{i=1}^{N}x_{i}^{\text{l}}P_{i}^{\text{sat}}=\sum _{i=1}^{N}x_{i}^{\text{g}}P=\left(\sum _{i=1}^{N}x_{i}^{\text{g}}\right)P=P$
La pression d'ébullition $P_{\text{b}}$ se calcule donc selon$^,$ :
* Pression d'ébullition : $P_{\text{b}}=\sum _{i=1}^{N}x_{i}^{\text{l}}P_{i}^{\text{sat}}$
Toutefois avec la correction de Poynting ${\mathcal {P}}_{i}\!\left(P\right)=\exp \!\left({\frac {\int _{P_{i}^{\text{sat}}}^{P}{\bar {V}}_{i}^{\text{l,*}}\!\left(P,T\right)\,\mathrm {d} P}{RT}}\right)$ le calcul n'est pas explicite. En effet, l'équation à résoudre en $P$ devient :
    $P_{\text{b}}=\sum _{i=1}^{N}x_{i}^{\text{l}}P_{i}^{\text{sat}}{\mathcal {P}}_{i}\!\left(P_{\text{b}}\right)$
On pourra employer la méthode du point fixe. Soit la fonction $f$ de $P$ :
    $f\!\left(P\right)=\sum _{i=1}^{N}x_{i}^{\text{l}}P_{i}^{\text{sat}}{\mathcal {P}}_{i}\!\left(P\right)$
On initialise le calcul par la solution obtenue par la loi de Raoult sans correction de Poynting :
    $P_{0}=\sum _{i=1}^{N}x_{i}^{\text{l}}P_{i}^{\text{sat}}$
puis on calcule la pression $P_{n+1}$ à partir de la pression $P_{n}$ obtenue à l'itération précédente par :
    $P_{n+1}=f\!\left(P_{n}\right)$
jusqu'à une erreur $\epsilon$ acceptable :
    $\left|P_{n+1}-P_{n}\right|<\epsilon$
Exemple
    On considère un mélange d'hydrocarbures aromatiques dans les conditions données par le tableau suivant.
Mélange d'hydrocarbures aromatiques à 50 °C - données. :
| Constituant $i$ | Fraction molaire liquide $x_{i}^{\text{l}}$ | Pression de vapeur saturante $P_{i}^{\text{sat}}$ en mmHg |
| Benzène | 0,4 | 271,29 |
| Toluène | 0,4 | 92,11 |
| Éthylbenzène | 0,2 | 35,16 |
    La pression de bulle est de $P_{\text{b}}=0{,}4\times 271{,}29+0{,}4\times 92{,}11+0{,}2\times 35{,}16=152{,}39\,{\mathsf {mmHg}}$.

#### Température d'ébullition

La température d'ébullition d'un mélange, à pression donnée, est la température en deçà de laquelle (pour des températures plus basses) le mélange est entièrement liquide : à la température d'ébullition apparaît la première bulle de vapeur. Au-delà de la température d'ébullition (pour des températures plus élevées), le mélange présente un équilibre liquide-vapeur, ou est entièrement gazeux (au-delà de la température de rosée). Soit un mélange liquide de $N$ composants, chaque composant $i$ étant représenté par sa fraction molaire $x_{i}^{\text{l}}$, à la pression $P$.
La solution n'est pas explicite, car l'expression de la loi de Raoult (avec ou sans correction de Poynting) n'est pas explicite en température. La pression étant donnée et la température étant l'inconnue, les pressions de vapeur saturante $P_{i}^{\text{sat}}=P_{i}^{\text{sat}}\!\left(T\right)$ et les corrections de Poynting ${\mathcal {P}}_{i}\!\left(T\right)=\exp \!\left({\frac {\int _{P_{i}^{\text{sat}}}^{P}{\bar {V}}_{i}^{\text{l,*}}\!\left(P,T\right)\,\mathrm {d} P}{RT}}\right)$ sont des fonctions de $T$. Si l'on pose la fonction $f$ de la température $T$ :
    $f\!\left(T\right)=\sum _{i=1}^{N}x_{i}^{\text{l}}P_{i}^{\text{sat}}\!\left(T\right){\mathcal {P}}_{i}\!\left(T\right)$
la température d'ébullition $T_{\text{b}}$ est la solution de l'équation : $f\!\left(T\right)=P$, la pression $P$ étant une donnée du problème. On pourra travailler par itérations imbriquées, en faisant varier la température $T$ et en calculant la pression d'ébullition correspondante jusqu'à ce que cette pression soit la pression de l'énoncé$^,$.

### Calcul d'un point de rosée
#### Pression de rosée

La pression de rosée d'un mélange, à température donnée, est la pression en deçà de laquelle (pour des pressions plus faibles) le mélange est entièrement gazeux : à la pression de rosée apparaît la première goutte de liquide. Au-delà de la pression de rosée (pour des pressions plus fortes), le mélange présente un équilibre liquide-vapeur, ou est entièrement liquide (au-delà de la pression d'ébullition). Soit un mélange gazeux de $N$ composants, chaque composant $i$ étant représenté par sa fraction molaire $x_{i}^{\text{g}}$, à la température $T$.
Nous connaissons la température $T$ du mélange, nous pouvons donc calculer la pression de vapeur saturante $P_{i}^{\text{sat}}$ de chacun des composants. La composition de la phase vapeur est connue. On suppose que tous les composants du mélange suivent la loi de Raoult. La loi de Raoult sans correction de Poynting donne pour tout constituant $i$ :
    ${x_{i}^{\text{g}} \over P_{i}^{\text{sat}}}={x_{i}^{\text{l}} \over P}$
et, en sommant sur l'ensemble des constituants :
    $\sum _{i=1}^{N}{x_{i}^{\text{g}} \over P_{i}^{\text{sat}}}={\sum _{i=1}^{N}x_{i}^{\text{l}} \over P}={1 \over P}$
La pression de rosée $P_{\text{r}}$ se calcule donc selon$^,$ :
* Pression de rosée : ${1 \over P_{\text{r}}}=\sum _{i=1}^{N}{x_{i}^{\text{g}} \over P_{i}^{\text{sat}}}$
Toutefois avec la correction de Poynting ${\mathcal {P}}_{i}\!\left(P\right)=\exp \!\left({\frac {\int _{P_{i}^{\text{sat}}}^{P}{\bar {V}}_{i}^{\text{l,*}}\!\left(P,T\right)\,\mathrm {d} P}{RT}}\right)$ le calcul n'est pas explicite. En effet, l'équation à résoudre en $P$ devient :
    ${1 \over P_{\text{r}}}=\sum _{i=1}^{N}{x_{i}^{\text{g}} \over P_{i}^{\text{sat}}{\mathcal {P}}_{i}\!\left(P_{\text{r}}\right)}$
On pourra employer la méthode du point fixe. Soit la fonction $f$ de $P$ :
    $f\!\left(P\right)=\sum _{i=1}^{N}{x_{i}^{\text{g}} \over P_{i}^{\text{sat}}{\mathcal {P}}_{i}\!\left(P\right)}$
On initialise le calcul par la solution obtenue par la loi de Raoult sans correction de Poynting :
    ${1 \over P_{0}}=\sum _{i=1}^{N}{x_{i}^{\text{g}} \over P_{i}^{\text{sat}}}$
puis on calcule la pression $P_{n+1}$ à partir de la pression $P_{n}$ obtenue à l'itération précédente par :
    ${1 \over P_{n+1}}=f\!\left(P_{n}\right)$
jusqu'à une erreur $\epsilon$ acceptable :
    $\left|P_{n+1}-P_{n}\right|<\epsilon$
Exemple
    On considère un mélange d'hydrocarbures aromatiques dans les conditions données par le tableau suivant.
Mélange d'hydrocarbures aromatiques à 50 °C - données. :
| Constituant $i$ | Fraction molaire vapeur $x_{i}^{\text{g}}$ | Pression de vapeur saturante $P_{i}^{\text{sat}}$ en mmHg |
| Benzène | 0,4 | 271,29 |
| Toluène | 0,4 | 92,11 |
| Éthylbenzène | 0,2 | 35,16 |
    La pression de rosée est de $P_{\text{r}}={1 \over {0{,}4 \over 271{,}29}+{0{,}4 \over 92{,}11}+{0{,}2 \over 35{,}16}}=86{,}92\,{\mathsf {mmHg}}$.

#### Température de rosée

La température de rosée d'un mélange, à pression donnée, est la température au-delà de laquelle (pour des températures plus élevées) le mélange est entièrement gazeux : à la température de rosée apparaît la première goutte de liquide. En deçà de la température de rosée (pour des températures plus basses), le mélange présente un équilibre liquide-vapeur, ou est entièrement liquide (en deçà de la température d'ébullition). Soit un mélange gazeux de $N$ composants, chaque composant $i$ étant représenté par sa fraction molaire $x_{i}^{\text{g}}$, à la pression $P$.
La solution n'est pas explicite, car l'expression de la loi de Raoult (avec ou sans correction de Poynting) n'est pas explicite en température. La pression étant donnée et la température étant l'inconnue, les pressions de vapeur saturante $P_{i}^{\text{sat}}=P_{i}^{\text{sat}}\!\left(T\right)$ et les corrections de Poynting ${\mathcal {P}}_{i}\!\left(T\right)=\exp \!\left({\frac {\int _{P_{i}^{\text{sat}}}^{P}{\bar {V}}_{i}^{\text{l,*}}\!\left(P,T\right)\,\mathrm {d} P}{RT}}\right)$ sont des fonctions de $T$. Si l'on pose la fonction $f$ de la température $T$ :
    $f\!\left(T\right)=\sum _{i=1}^{N}{x_{i}^{\text{g}} \over P_{i}^{\text{sat}}\!\left(T\right){\mathcal {P}}_{i}\!\left(T\right)}$
la température de rosée $T_{\text{r}}$ est la solution de l'équation : $f\!\left(T\right)={1 \over P}$, la pression $P$ étant une donnée du problème. On pourra travailler par itérations imbriquées, en faisant varier la température $T$ et en calculant la pression de rosée correspondante jusqu'à ce que cette pression soit la pression de l'énoncé$^,$.

### Calcul d'un équilibre liquide-vapeur
#### Énoncé du problème

Soit, à pression $P$ et température $T$ données, un mélange de $N$ constituants, chaque constituant $i$ étant représenté par la quantité $n_{i}$. On cherche à déterminer la composition des phases vapeur et liquide en présence.
Nous avons $2\cdot N$ inconnues, il s'agit en effet de calculer pour chacun des $N$ constituants :
* $n_{i}^{\text{g}}$ la quantité du constituant $i$ en phase vapeur ;
* $n_{i}^{\text{l}}$ la quantité du constituant $i$ en phase liquide.
Soient :
* $n^{\text{g}}=\sum _{i=1}^{N}n_{i}^{\text{g}}$ la quantité de matière totale en phase vapeur (inconnue) ;
* $n^{\text{l}}=\sum _{i=1}^{N}n_{i}^{\text{l}}$ la quantité de matière totale en phase liquide (inconnue) ;
* $x_{i}^{\text{g}}={n_{i}^{\text{g}} \over n^{\text{g}}}$ la fraction molaire du constituant $i$ en phase vapeur (inconnue) ;
* $x_{i}^{\text{l}}={n_{i}^{\text{l}} \over n^{\text{l}}}$ la fraction molaire du constituant $i$ en phase liquide (inconnue).
On suppose que tous les composants du mélange suivent la loi de Raoult. La pression et la température étant données, les pressions de vapeur saturante $P_{i}^{\text{sat}}=P_{i}^{\text{sat}}\!\left(T\right)$ et les corrections de Poynting ${\mathcal {P}}_{i}\!\left(P,T\right)=\exp \!\left({\frac {\int _{P_{i}^{\text{sat}}}^{P}{\bar {V}}_{i}^{\text{l,*}}\!\left(P,T\right)\,\mathrm {d} P}{RT}}\right)$ sont des constantes. On définit pour chacun des $N$ constituants le coefficient de partage $K_{i}$, qui est une constante en appliquant la loi de Raoult à pression et température données :
$K_{i}={x_{i}^{\text{g}} \over x_{i}^{\text{l}}}={P_{i}^{\text{sat}}{\mathcal {P}}_{i} \over P}$
Nous avons les $2\cdot N$ équations :
* équilibres selon la loi de Raoult : $x_{i}^{\text{g}}=K_{i}x_{i}^{\text{l}}$ ;
* bilans de matière : $n_{i}=n_{i}^{\text{l}}+n_{i}^{\text{g}}$.

#### Réduction du problème

Soit $n$ la quantité de matière totale considérée, il s'agit d'une valeur connue du problème. On a, par bilan de matière global :
    $n=\sum _{i=1}^{N}n_{i}=n^{\text{g}}+n^{\text{l}}$
Soit pour chacun des $N$ constituants la fraction molaire globale, notée $x_{i}$ :
$x_{i}={n_{i} \over n}$
Puisque l'on connaît la quantité de matière $n_{i}$ de chacun des $N$ constituants, les fractions molaires $x_{i}$ sont connues également.
On définit d'autre part $\tau ^{\text{g}}$ la fraction molaire globale de la phase vapeur, ou taux de vaporisation :
$\tau ^{\text{g}}={n^{\text{g}} \over n^{\text{g}}+n^{\text{l}}}={\sum _{i=1}^{N}n_{i}^{\text{g}} \over n}$
Cette grandeur est une inconnue du problème. Le système de $2\cdot N$ équations à $2\cdot N$ inconnues est réduit à une équation à une inconnue, $\tau ^{\text{g}}$, due à Rachford et Rice (1952)$^,$$^,$$^,$ :
* Équation de Rachford-Rice : $f\!\left(\tau ^{\text{g}}\right)=\sum _{i=1}^{N}{\left(K_{i}-1\right)x_{i} \over 1+\left(K_{i}-1\right)\tau ^{\text{g}}}=0$
On a par définition des fractions molaires dans les deux phases :
    $n_{i}^{\text{g}}=x_{i}^{\text{g}}n^{\text{g}}$
    $n_{i}^{\text{l}}=x_{i}^{\text{l}}n^{\text{l}}$
En introduisant ces relations dans le bilan de matière de chacun des $N$ constituants :
    $n_{i}=n_{i}^{\text{l}}+n_{i}^{\text{g}}=x_{i}^{\text{l}}n^{\text{l}}+x_{i}^{\text{g}}n^{\text{g}}=x_{i}^{\text{l}}\left(n^{\text{l}}+n^{\text{g}}-n^{\text{g}}\right)+x_{i}^{\text{g}}n^{\text{g}}$
et en divisant par $n^{\text{g}}+n^{\text{l}}$ on obtient pour chacun des $N$ constituants :
    $x_{i}=x_{i}^{\text{l}}\left(1-\tau ^{\text{g}}\right)+x_{i}^{\text{g}}\tau ^{\text{g}}$
En substituant dans cette dernière relation l'équation d'équilibre chimique du constituant $i$ :
    $x_{i}^{\text{g}}=K_{i}x_{i}^{\text{l}}$
    $x_{i}=x_{i}^{\text{l}}\left(1-\tau ^{\text{g}}\right)+K_{i}\tau ^{\text{g}}x_{i}^{\text{l}}=\left(1+\left(K_{i}-1\right)\tau ^{\text{g}}\right)x_{i}^{\text{l}}$
On remarque que si $1+\left(K_{i}-1\right)\tau ^{\text{g}}=0$, alors $x_{i}=0$ : le constituant $i$ n'est pas présent dans le mélange ; ce cas est donc exclu. On peut donc écrire :
    $x_{i}^{\text{l}}={x_{i} \over 1+\left(K_{i}-1\right)\tau ^{\text{g}}}$
d'où également, toujours d'après la relation de l'équilibre chimique :
    $x_{i}^{\text{g}}={K_{i}x_{i} \over 1+\left(K_{i}-1\right)\tau ^{\text{g}}}$
Dans les expressions des fractions $x_{i}^{\text{g}}$ et $x_{i}^{\text{l}}$ ainsi établies, seul le taux de vaporisation $\tau ^{\text{g}}$ est inconnu. Nous pouvons donc réduire le problème à la recherche de cette seule inconnue. Soit la fonction $f$ de $\tau ^{\text{g}}$ définie par :
    $f\!\left(\tau ^{\text{g}}\right)=\sum _{i=1}^{N}x_{i}^{\text{g}}-\sum _{i=1}^{N}x_{i}^{\text{l}}=\sum _{i=1}^{N}{\left(K_{i}-1\right)x_{i} \over 1+\left(K_{i}-1\right)\tau ^{\text{g}}}$
Cette fonction est monotone strictement décroissante en $\tau ^{\text{g}}$ pour $0<\tau ^{\text{g}}<1$. Si l'on trouve sur cette plage une solution à l'équation $f\!\left(\tau ^{\text{g}}\right)=0$, elle est unique et $\sum _{i=1}^{N}x_{i}^{\text{g}}=\sum _{i=1}^{N}x_{i}^{\text{l}}$.
Il reste à prouver que cette solution respecte les contraintes sur les fractions molaires, soit $\sum _{i=1}^{N}x_{i}^{\text{g}}=1$ et $\sum _{i=1}^{N}x_{i}^{\text{l}}=1$. Revenons à la relation pour chacun des constituants :
    $x_{i}=x_{i}^{\text{l}}\left(1-\tau ^{\text{g}}\right)+x_{i}^{\text{g}}\tau ^{\text{g}}$
En sommant sur l'ensemble des constituants :
    $\sum _{i=1}^{N}x_{i}=\sum _{i=1}^{N}x_{i}^{\text{l}}\left(1-\tau ^{\text{g}}\right)+\sum _{i=1}^{N}x_{i}^{\text{g}}\tau ^{\text{g}}$
or puisque $\sum _{i=1}^{N}x_{i}^{\text{g}}=\sum _{i=1}^{N}x_{i}^{\text{l}}$, alors :
    $\sum _{i=1}^{N}x_{i}=\sum _{i=1}^{N}x_{i}^{\text{g}}\left(1-\tau ^{\text{g}}\right)+\sum _{i=1}^{N}x_{i}^{\text{g}}\tau ^{\text{g}}=\sum _{i=1}^{N}x_{i}^{\text{g}}\left(1-\tau ^{\text{g}}+\tau ^{\text{g}}\right)=\sum _{i=1}^{N}x_{i}^{\text{g}}$
Par construction $\sum _{i=1}^{N}x_{i}=1$, d'où $\sum _{i=1}^{N}x_{i}^{\text{g}}=1$ et $\sum _{i=1}^{N}x_{i}^{\text{l}}=1$.
On peut résoudre cette équation sur l'intervalle $0\leq \tau ^{\text{g}}\leq 1$ par la méthode de Newton, en considérant la dérivée $f'$ de la fonction $f$ :
    $f'\!\left(\tau ^{\text{g}}\right)=-\sum _{i=1}^{N}{\left(K_{i}-1\right)^{2}x_{i} \over \left(1+\left(K_{i}-1\right)\tau ^{\text{g}}\right)^{2}}$
On remarque que $f'<0$, ce qui indique que $f$ est décroissante.
On initialise le calcul avec $\tau _{0}^{\text{g}}=0,5$, puis on procède par itérations successives selon :
    $\tau _{n+1}^{\text{g}}=\tau _{n}^{\text{g}}-{f\!\left(\tau _{n}^{\text{g}}\right) \over f'\!\left(\tau _{n}^{\text{g}}\right)}$
jusqu'à une erreur $\epsilon$ acceptable :
    $\left|\tau _{n+1}^{\text{g}}-\tau _{n}^{\text{g}}\right|<\epsilon$
Les fractions molaires sont ensuite calculées par :
    $x_{i}^{\text{g}}={K_{i}x_{i} \over 1+\left(K_{i}-1\right)\tau ^{\text{g}}}$
    $x_{i}^{\text{l}}={x_{i} \over 1+\left(K_{i}-1\right)\tau ^{\text{g}}}$
et les quantités globales respectives des deux phases par :
    $n^{\text{g}}=\tau ^{\text{g}}n$
    $n^{\text{l}}=\left(1-\tau ^{\text{g}}\right)n$
d'où l'on obtient finalement les quantités dans chaque phase pour chacun des $N$ constituants :
    $n_{i}^{\text{g}}=x_{i}^{\text{g}}n^{\text{g}}$
    $n_{i}^{\text{l}}=x_{i}^{\text{l}}n^{\text{l}}$

##### Note - autre forme de l'équation de Rachford-Rice

    Plutôt que de rechercher le taux de vaporisation $\tau ^{\text{g}}$, on peut rechercher le taux de liquéfaction $\tau ^{\text{l}}$ défini par :
       $\tau ^{\text{l}}={n^{\text{l}} \over n^{\text{l}}+n^{\text{g}}}=1-\tau ^{\text{g}}$
    L'équation de Rachford-Rice se présente alors sous la forme :
       $f\left(\tau ^{\text{l}}\right)=\sum _{i=1}^{N}{\left(K_{i}-1\right)x_{i} \over K_{i}-\left(K_{i}-1\right)\tau ^{\text{l}}}=0$
    Cette fonction est monotone strictement croissante en $\tau ^{\text{l}}$ pour $0\leq \tau ^{\text{l}}\leq 1$.

#### Résolution du problème

Avant de calculer un équilibre liquide-vapeur, il est nécessaire de vérifier que le mélange n'est pas monophasique, soit uniquement liquide ou gazeux.
À température donnée, le coefficient de partage de chacun des $N$ constituants est une fonction de la pression :
    $K_{i}\!\left(P\right)={P_{i}^{\text{sat}}{\mathcal {P}}_{i}\!\left(P\right) \over P}$
Lorsque la pression est supérieure à la pression d'ébullition, soit $P>P_{\text{b}}$, le mélange est entièrement liquide. Pour une composition donnée, la fonction $f\!\left(P\right)=\sum _{i=1}^{N}x_{i}K_{i}\!\left(P\right)$ est une fonction décroissante de la pression. Lorsque $f\!\left(P\right)=1$ la pression est la pression d'ébullition $P_{\text{b}}$ aux température et composition données. Lorsque $P>P_{\text{b}}$ on a $f\!\left(P\right)<1$.
Lorsque la pression est inférieure à la pression de rosée, soit $P<P_{\text{r}}$, le mélange est entièrement gazeux. Pour une composition donnée, la fonction $g\!\left(P\right)=\sum _{i=1}^{N}{x_{i} \over K_{i}\!\left(P\right)}$ est une fonction croissante de la pression. Lorsque $g\!\left(P\right)=1$ la pression est la pression de rosée $P_{\text{r}}$ aux température et composition données. Lorsque $P<P_{\text{r}}$ on a $g\!\left(P\right)<1$.

##### Algorithme de calcul

On effectuera par conséquent le calcul selon l'algorithme suivant :
* à pression, température et composition de l'énoncé calculer la fraction molaire globale $x_{i}$ et le coefficient de partage $K_{i}$ pour chacun des $N$ constituants ;
* si $\sum _{i=1}^{N}x_{i}K_{i}<1$ le mélange est entièrement liquide, $\tau ^{\text{g}}=0$ et pour tout constituant $x_{i}^{\text{l}}=x_{i}$, fin du calcul ;
* si $\sum _{i=1}^{N}{x_{i} \over K_{i}}<1$ le mélange est entièrement gazeux, $\tau ^{\text{g}}=1$ et pour tout constituant $x_{i}^{\text{g}}=x_{i}$, fin du calcul ;
* sinon, le mélange est biphasique et présente un équilibre liquide-vapeur, résoudre l'équation de Rachford-Rice : $\sum _{i=1}^{N}{\left(K_{i}-1\right)x_{i} \over 1+\left(K_{i}-1\right)\tau ^{\text{g}}}=0$.
Dans le cas d'un mélange binaire, impliquant les constituants ${\text{A}}$ et ${\text{B}}$, après avoir vérifié la présence de deux phases comme précédemment, on calcule plus rapidement les compositions des phases selon :
* $x_{\text{A}}^{\text{l}}={1-K_{\text{B}} \over K_{\text{A}}-K_{\text{B}}}$ ;
* $x_{\text{B}}^{\text{l}}=1-x_{\text{A}}^{\text{l}}$ ;
* $x_{\text{A}}^{\text{g}}=K_{\text{A}}x_{\text{A}}^{\text{l}}$ ;
* $x_{\text{B}}^{\text{g}}=1-x_{\text{A}}^{\text{g}}$ ;
* $\tau ^{\text{g}}={x_{\text{A}}-x_{\text{A}}^{\text{l}} \over x_{\text{A}}^{\text{g}}-x_{\text{A}}^{\text{l}}}$.
Exemple
    On considère un mélange d'hydrocarbures dans les conditions données par le tableau suivant.
Mélange d'hydrocarbures à 3,5 bar et 300 K - données. :
| Constituant $i$ | Fraction molaire globale $x_{i}$ | Coefficient de partage $K_{i}$ |
| Propane | 0,2 | 2,525 |
| n-Butane | 0,3 | 0,7708 |
| Isobutane | 0,4 | 1,066 |
| n-Pentane | 0,05 | 0,2401 |
| Isopentane | 0,05 | 0,3140 |
    On a $\sum _{i=1}^{N}x_{i}K_{i}=1{,}1903$ et $\sum _{i=1}^{N}{x_{i} \over K_{i}}=1{,}2113$. Le mélange est donc biphasique. Les résultats sont reportés dans le tableau suivant.
Mélange d'hydrocarbures à 3,5 bar et 300 K - résultats. :
| Taux de vaporisation $\tau ^{\text{g}}$ | 0,524 | 0,524 |
| Constituant $i$ | Fraction molaire liquide $x_{i}^{\text{l}}$ | Fraction molaire vapeur $x_{i}^{\text{g}}$ |
| Propane | 0,1112 | 0,2807 |
| n-Butane | 0,3409 | 0,2628 |
| Isobutane | 0,3867 | 0,4121 |
| n-Pentane | 0,0831 | 0,0199 |
| Isopentane | 0,0781 | 0,0245 |

### Diagrammes de phase binaires idéaux
#### Principe

La relation de Duhem-Margules induit que dans un mélange binaire si l'un des corps suit la loi de Raoult alors le deuxième suit également la loi de Raoult (ou la loi de Henry). On considère par conséquent deux corps ${\text{A}}$ et ${\text{B}}$ suivant la loi de Raoult à l'équilibre liquide-vapeur. On a six inconnues :
* la température $T$ ;
* la pression $P$ ;
* la fraction molaire du corps ${\text{A}}$ en phase vapeur $x_{\text{A}}^{\text{g}}$ ;
* la fraction molaire du corps ${\text{A}}$ en phase liquide $x_{\text{A}}^{\text{l}}$ ;
* la fraction molaire du corps ${\text{B}}$ en phase vapeur $x_{\text{B}}^{\text{g}}$ ;
* la fraction molaire du corps ${\text{B}}$ en phase liquide $x_{\text{B}}^{\text{l}}$ ;
et quatre équations :
* la loi de Raoult pour le corps ${\text{A}}$ : $x_{\text{A}}^{\text{g}}P=x_{\text{A}}^{\text{l}}P_{\text{A}}^{\text{sat}}$ ;
* la loi de Raoult pour le corps ${\text{B}}$ : $x_{\text{B}}^{\text{g}}P=x_{\text{B}}^{\text{l}}P_{\text{B}}^{\text{sat}}$ ;
* la contrainte sur les fractions molaires de la phase vapeur : $x_{\text{A}}^{\text{g}}+x_{\text{B}}^{\text{g}}=1$ ;
* la contrainte sur les fractions molaires de la phase liquide : $x_{\text{A}}^{\text{l}}+x_{\text{B}}^{\text{l}}=1$.
Les pressions de vapeur saturante $P_{\text{A}}^{\text{sat}}$ et $P_{\text{B}}^{\text{sat}}$ sont des fonctions de la température.
On a donc deux degrés de liberté, conformément à la règle des phases : on peut étudier l'équilibre liquide-vapeur de ce système en fixant deux paramètres parmi les six inconnues définies plus haut. Par exemple, en fixant la température (diagrammes isothermes), on peut :
* fixer la composition du liquide et calculer la pression : on obtient le point d'ébullition ;
* fixer la composition du gaz et calculer la pression : on obtient le point de rosée ;
* fixer la composition du liquide et calculer la composition du gaz (et inversement).
On peut de façon similaire fixer la pression (diagrammes isobares) et la composition de l'une des phases pour calculer la température d'équilibre ou la composition de l'autre phase.

#### Diagrammes isothermes

La température $T$ est fixée. On veut tracer l'évolution de la pression $P$ en fonction de la composition.
Puisque la température est fixée, les pressions de vapeur saturante $P_{\text{A}}^{\text{sat}}$ et $P_{\text{B}}^{\text{sat}}$ sont constantes.
En considérant la contrainte sur la composition de la vapeur et la loi de Raoult :
    $1=x_{\text{A}}^{\text{g}}+x_{\text{B}}^{\text{g}}=x_{\text{A}}^{\text{l}}{P_{\text{A}}^{\text{sat}} \over P}+x_{\text{B}}^{\text{l}}{P_{\text{B}}^{\text{sat}} \over P}$
    $P=x_{\text{A}}^{\text{l}}P_{\text{A}}^{\text{sat}}+x_{\text{B}}^{\text{l}}P_{\text{B}}^{\text{sat}}$
Si l'on choisit, arbitrairement, le corps ${\text{B}}$ comme référence, avec la contrainte sur la composition liquide on a :
    $P=\left(1-x_{\text{B}}^{\text{l}}\right)P_{\text{A}}^{\text{sat}}+x_{\text{B}}^{\text{l}}P_{\text{B}}^{\text{sat}}$
On obtient la formule de la courbe d'ébullition :
$P=\left(P_{\text{B}}^{\text{sat}}-P_{\text{A}}^{\text{sat}}\right)x_{\text{B}}^{\text{l}}+P_{\text{A}}^{\text{sat}}$
En considérant la contrainte sur la composition du liquide et la loi de Raoult :
    $1=x_{\text{A}}^{\text{l}}+x_{\text{B}}^{\text{l}}=x_{\text{A}}^{\text{g}}{P \over P_{\text{A}}^{\text{sat}}}+x_{\text{B}}^{\text{g}}{P \over P_{\text{B}}^{\text{sat}}}$
    ${1 \over P}={x_{\text{A}}^{\text{g}} \over P_{\text{A}}^{\text{sat}}}+{x_{\text{B}}^{\text{g}} \over P_{\text{B}}^{\text{sat}}}$
Si l'on choisit, arbitrairement, le corps ${\text{B}}$ comme référence, avec la contrainte sur la composition vapeur on a :
    ${1 \over P}={1-x_{\text{B}}^{\text{g}} \over P_{\text{A}}^{\text{sat}}}+{x_{\text{B}}^{\text{g}} \over P_{\text{B}}^{\text{sat}}}=\left({1 \over P_{\text{B}}^{\text{sat}}}-{1 \over P_{\text{A}}^{\text{sat}}}\right)x_{\text{B}}^{\text{g}}+{1 \over P_{\text{A}}^{\text{sat}}}$
On obtient la formule de la courbe de rosée :
$P={P_{\text{A}}^{\text{sat}}P_{\text{B}}^{\text{sat}} \over \left(P_{\text{A}}^{\text{sat}}-P_{\text{B}}^{\text{sat}}\right)x_{\text{B}}^{\text{g}}+P_{\text{B}}^{\text{sat}}}$
La courbe d'ébullition est un segment de droite. La courbe de rosée est un arc d'hyperbole. Les deux courbes ne peuvent donc s'intercepter qu'en deux points, ces points correspondent aux points d'ébullition des deux corps purs, soit $\left(x_{\text{B}}^{\text{l}}=x_{\text{B}}^{\text{g}}=0,P=P_{\text{A}}^{\text{sat}}\right)$ et $\left(x_{\text{B}}^{\text{l}}=x_{\text{B}}^{\text{g}}=1,P=P_{\text{B}}^{\text{sat}}\right)$. La loi de Raoult est ainsi incapable de représenter un azéotrope, point auquel le liquide et la vapeur d'un mélange binaire ont la même composition.
Il est également possible de tracer la courbe de composition d'une phase en fonction de la composition de l'autre phase, par exemple de la vapeur en fonction du liquide : $x_{\text{B}}^{\text{g}}=x_{\text{B}}^{\text{g}}\!\left(x_{\text{B}}^{\text{l}}\right)$. En effet, à la température $T$ fixée la pression est donnée par les deux relations :
    $P=\left(P_{\text{B}}^{\text{sat}}-P_{\text{A}}^{\text{sat}}\right)x_{\text{B}}^{\text{l}}+P_{\text{A}}^{\text{sat}}={P_{\text{A}}^{\text{sat}}P_{\text{B}}^{\text{sat}} \over \left(P_{\text{A}}^{\text{sat}}-P_{\text{B}}^{\text{sat}}\right)x_{\text{B}}^{\text{g}}+P_{\text{B}}^{\text{sat}}}$
On obtient par conséquent :
$x_{\text{B}}^{\text{g}}={x_{\text{B}}^{\text{l}}P_{\text{B}}^{\text{sat}} \over \left(P_{\text{B}}^{\text{sat}}-P_{\text{A}}^{\text{sat}}\right)x_{\text{B}}^{\text{l}}+P_{\text{A}}^{\text{sat}}}$

#### Diagrammes isobares

La pression $P$ est fixée. On veut tracer l'évolution de la température $T$ en fonction de la composition.
Les pressions de vapeur saturante $P_{\text{A}}^{\text{sat}}$ et $P_{\text{B}}^{\text{sat}}$ sont des fonctions de la température. Il est donc nécessaire de disposer de corrélations permettant cette évaluation. En intégrant la formule de Clausius-Clapeyron, la pression de vapeur saturante d'un corps pur $i$ peut être écrite, par exemple, selon la formule de Rankine :
    $\ln \!\left({P_{i}^{\text{sat}} \over P}\right)=-{\Delta _{\text{vap}}H_{i} \over R}\left[{1 \over T}-{1 \over T_{{\text{éb}},i}}\right]$
avec :
* $T_{{\text{éb}},i}$ la température d'ébullition du corps pur $i$ sous la pression $P$, il s'agit d'une constante ;
* $P_{i}^{\text{sat}}$ la pression de vapeur saturante du corps pur $i$ à la température $T$ ; il s'agit donc d'une fonction telle que $P_{i}^{\text{sat}}=P_{i}^{\text{sat}}\!\left(T\right)$, et en particulier $P=P_{i}^{\text{sat}}\!\left(T_{{\text{éb}},i}\right)$ ;
* $\Delta _{\text{vap}}H_{i}$ l'enthalpie de vaporisation du corps pur $i$ à la température $T_{{\text{éb}},i}$ (donc sous la pression $P$), supposée constante dans l'établissement de la formule de Rankine (l'enthalpie de vaporisation dépend normalement de la température).
On pose les coefficients de partage :
    $K_{\text{A}}\!\left(T\right)={P_{\text{A}}^{\text{sat}}\!\left(T\right) \over P}={x_{\text{A}}^{\text{g}} \over x_{\text{A}}^{\text{l}}}$
    $K_{\text{B}}\!\left(T\right)={P_{\text{B}}^{\text{sat}}\!\left(T\right) \over P}={x_{\text{B}}^{\text{g}} \over x_{\text{B}}^{\text{l}}}$
En considérant la contrainte sur la composition de la vapeur et la loi de Raoult :
    $1=x_{\text{A}}^{\text{g}}+x_{\text{B}}^{\text{g}}=x_{\text{A}}^{\text{l}}K_{\text{A}}+x_{\text{B}}^{\text{l}}K_{\text{B}}$
Avec la contrainte sur la composition du liquide on a :
    $1=x_{\text{A}}^{\text{l}}K_{\text{A}}+\left(1-x_{\text{A}}^{\text{l}}\right)K_{\text{B}}=\left(1-x_{\text{B}}^{\text{l}}\right)K_{\text{A}}+x_{\text{B}}^{\text{l}}K_{\text{B}}$
d'où :
$x_{\text{A}}^{\text{l}}={1-K_{\text{B}} \over K_{\text{A}}-K_{\text{B}}}$
$x_{\text{B}}^{\text{l}}={1-K_{\text{A}} \over K_{\text{B}}-K_{\text{A}}}$
Avec la loi de Raoult on obtient :
$x_{\text{A}}^{\text{g}}=K_{\text{A}}{1-K_{\text{B}} \over K_{\text{A}}-K_{\text{B}}}$
$x_{\text{B}}^{\text{g}}=K_{\text{B}}{1-K_{\text{A}} \over K_{\text{B}}-K_{\text{A}}}$
Ainsi, en imposant la température $T$, la pression $P$ étant donnée, on peut calculer $K_{\text{A}}$ et $K_{\text{B}}$, puis $x_{\text{A}}^{\text{l}}$, $x_{\text{B}}^{\text{l}}$, $x_{\text{A}}^{\text{g}}$ et $x_{\text{B}}^{\text{g}}$. Si l'on choisit, arbitrairement, le corps ${\text{B}}$ comme référence, en faisant varier la température $T$ entre $T_{{\text{éb}},{\text{A}}}$ et $T_{{\text{éb}},{\text{B}}}$, la courbe $T=T\!\left(x_{\text{B}}^{\text{l}}\right)$ est la courbe d'ébullition, la courbe $T=T\!\left(x_{\text{B}}^{\text{g}}\right)$ est la courbe de rosée.
Il est également possible de tracer la courbe de composition d'une phase en fonction de la composition de l'autre phase, par exemple de la vapeur en fonction du liquide : $x_{\text{B}}^{\text{g}}=x_{\text{B}}^{\text{g}}\!\left(x_{\text{B}}^{\text{l}}\right)$. Il suffit, à la pression $P$ fixée :
* de poser $x_{\text{B}}^{\text{l}}$,
* de trouver $T$ entre $T_{{\text{éb}},{\text{A}}}$ et $T_{{\text{éb}},{\text{B}}}$ telle que $x_{\text{B}}^{\text{l}}={1-K_{\text{A}} \over K_{\text{B}}-K_{\text{A}}}$,
* de calculer $x_{\text{B}}^{\text{g}}=K_{\text{B}}x_{\text{B}}^{\text{l}}$.