Équilibre de phases

Fichier:Binary phase diagram of NiTI (phase).JPG

Diagramme de phases du couple titane (Ti) - nickel (Ni). Les courbes délimitant les divers domaines représentent les conditions de transition de phase. Dans certains domaines, plusieurs phases peuvent coexister dans un équilibre stable.

En physique, et plus particulièrement en thermodynamique, un équilibre de phases est un état de la matière dans lequel deux phases ou plus coexistent sans évolution.

Un diagramme de phases est un diagramme représentant les domaines d'existence et de coexistence des phases d'un corps pur ou d'un mélange en fonction des conditions opératoires de pression, température et composition. La règle des phases permet de déterminer le nombre maximal de phases pouvant coexister selon les conditions opératoires, ainsi que le nombre de conditions qu'un opérateur peut faire évoluer indépendamment l'une de l'autre sans provoquer de transition de phase.

Deux phases sont en équilibre si elles ont atteint simultanément :

• un équilibre mécanique : elles ont la même pression ;
• un équilibre thermique : elles ont la même température ;
• un équilibre diffusif : le potentiel chimique de chaque espèce chimique est homogène entre les deux phases.

Cette dernière condition est équivalente, pour chaque espèce, à l'homogénéité de sa fugacité entre les phases.

Sous l'hypothèse simplificatrice de la solution idéale, l'équilibre liquide-vapeur peut être calculé par les loi de Raoult et loi de Henry. L'équation de Schröder-van Laar permet de calculer les équilibres liquide-solide. Les lois permettant de calculer les propriétés colligatives donnent l'évolution des pression et température d'équilibre de phases d'un solvant liquide en présence d'un soluté très dilué.

Définitions

États et phases

États de la matière et changements d'état

• 
  Diagramme de phase de l'eau montrant l'état de l'eau pure en fonction de la température et de la pression. Au-delà du point critique, le liquide et la vapeur ne se distinguent plus. Au point triple les trois états coexistent.
• Principaux changements d'état de la matière. Un changement d'état implique une enthalpie de changement d'état.
  Principaux changements d'état de la matière. Un changement d'état implique une enthalpie de changement d'état.
• Diagramme de phase du mélange bismuth-antimoine montrant, en fonction de la température et de la composition, les domaines d'existence de l'état liquide (Liq), de l'état solide (S) et d'un domaine de coexistence des deux états liquide et solide (Liq+S).
  Diagramme de phase du mélange bismuth-antimoine montrant, en fonction de la température et de la composition, les domaines d'existence de l'état liquide (Liq), de l'état solide (S) et d'un domaine de coexistence des deux états liquide et solide (Liq+S).

Selon les conditions de pression, température et composition (conditions opératoires), la matière peut se présenter sous plusieurs formes appelées états de la matière, dont les principaux, étudiés en thermodynamique, sont le gaz, le liquide et le solide. L'eau, par exemple, peut se présenter sous forme de glace, qui est son état solide, de liquide, ou de vapeur, qui est son état gazeux. Aux hautes pressions et températures, au-delà du point critique, les états gazeux et liquide ne se distinguent plus : cet état est appelé état supercritique. Le point critique de l'eau se situe à Modèle:Unité et Modèle:Unité.

Le passage d'un état à un autre est appelé changement d'état. Les principaux changements d'état sont : la vaporisation, de l'état liquide à l'état gazeux, la liquéfaction ou condensation liquide, de l'état gazeux à l'état liquide, la fusion, de l'état solide à l'état liquide, la solidification, de l'état liquide à l'état solide, la sublimation, de l'état solide à l'état gazeux, et la condensation solide, de l'état gazeux à l'état solide.

Un changement d'état d'un corps pur s'effectue à température et pression constantes. À température constante, un changement d'état d'un mélange s'effectue à pression variable ; à pression constante il s'effectue à température variable. Il existe des cas particuliers dans lesquels le changement d'état d'un mélange s'effectue à pression et température constantes, le mélange se comportant alors comme un corps purs : c'est notamment le cas des azéotropes pour les changements liquide-vapeur, et des points de fusion congruents pour les changements liquide-solide. Dans les hétéroazéotropes, un gaz se condense en deux phases liquides distinctes, de même que dans des eutectiques, un liquide se condense en deux phases solides distinctes.

Phases et transitions de phases

• Diagramme de phases de l'eau : en bleu l'état solide, en vert l'état liquide, en orange l'état gazeux. Les chiffres romains correspondent aux diverses phases de glace (état solide).
  Diagramme de phases de l'eau : en bleu l'état solide, en vert l'état liquide, en orange l'état gazeux. Les chiffres romains correspondent aux diverses phases de glace (état solide).
• Le diamant (à gauche) et le graphite (à droite) sont deux formes solides du carbone, constituant autant de phases solides distinctes. Dans les conditions de pression et température atmosphériques ordinaires, le graphite est stable, le diamant métastable.
  Le diamant (à gauche) et le graphite (à droite) sont deux formes solides du carbone, constituant autant de phases solides distinctes. Dans les conditions de pression et température atmosphériques ordinaires, le graphite est stable, le diamant métastable.
• Diagramme de phases cuivre-argent. α est une phase solide de cuivre contenant de l'argent, β une phase solide d'argent contenant du cuivre. Les deux phases coexistent dans le domaine α+β. Le point E est un eutectique.
  Diagramme de phases cuivre-argent. $\alpha$ est une phase solide de cuivre contenant de l'argent, $\beta$ une phase solide d'argent contenant du cuivre. Les deux phases coexistent dans le domaine $\alpha +\beta$. Le point E est un eutectique.

Selon les conditions de pression, température et composition (conditions opératoires), un même état peut se présenter sous plusieurs formes appelées phases. Pour un corps pur (composé d'une seule espèce chimique) comme pour un mélange (composé de plusieurs espèces chimiques), il n'existe qu'une seule phase gazeuse : la diffusion de la matière induit naturellement que deux gaz mis en contact, quels qu'ils soient, se mélangent intimement pour former un système homogène.

Les corps purs ne présentent qu'une seule phase liquide, quelles que soient les conditions opératoires, à l'exception notable de certains isotopes de l'hélium présentant, en plus du liquide ordinaire, une phase liquide superfluide aux très basses températures (quelques kelvins). Les mélanges peuvent présenter plusieurs phases liquides simultanément. Par exemple, un mélange d'huile et d'eau constitue deux phases liquides distinctes : une phase liquide organique et une phase liquide aqueuse non miscibles (dans les faits, chacune des deux phases contient une infime quantité de l'autre phase dissoute).

Les corps purs et les mélanges peuvent présenter plusieurs phases solides. L'eau pure solide (glace), par exemple, présente, selon les conditions opératoires, plus de onze formes cristallines qui constituent autant de phases solides distinctes. Les mélanges (alliages) de métaux présentent souvent, selon les conditions opératoires, de multiples formes cristallines qui sont autant de phases solides. La faculté d'un corps simple de se présenter sous plusieurs formes solides est appelée allotropie (cas du carbone avec, entre autres, le graphite et le diamant)Modèle:,, pour un corps composé elle est appelée polymorphisme (cas du carbonate de calcium avec l'aragonite, la calcite et la vatérite).

Lorsque les conditions opératoires sont modifiées, un matériau peut changer de phase. Ce changement est appelé transition de phase. Un changement d'état est une transition de phase, mais toutes les transitions de phase ne s'accompagnent pas nécessairement d'un changement d'état. Par exemple, la transition d'une forme de glace de l'eau à une autre est une transition de phase dans laquelle l'état de l'eau reste l'état solide.

Continuité, homogénéité et dispersion

• Erreur lors de la création de la vignette :
  Homogénéité
  Cette hypothèse est simplificatrice.
  Dans la colonne d'un baromètre à mercure, la pression est quasi nulle au sommet et égale à la pression atmosphérique au pied.
• Dispersion Une phase n'est pas nécessairement d'un seul tenant. Un précipité est une phase solide dispersée dans une phase liquide en d'innombrables cristaux.
  Dispersion
  Une phase n'est pas nécessairement d'un seul tenant.
  Un précipité est une phase solide dispersée dans une phase liquide en d'innombrables cristaux.
• Continuité Une discontinuité de propriété lorsque l'on modifie une condition opératoire correspond à une transition de phase. Masse volumique de l'eau en fonction de la température. À Modèle:Unité la masse volumique augmente brusquement : la glace (solide) fond et devient liquide. L'eau est un rare exemple de liquide plus dense que le solide (la glace flotte).
  Continuité
  Une discontinuité de propriété lorsque l'on modifie une condition opératoire correspond à une transition de phase.
  Masse volumique de l'eau en fonction de la température. À Modèle:Unité la masse volumique augmente brusquement : la glace (solide) fond et devient liquide. L'eau est un rare exemple de liquide plus dense que le solide (la glace flotte).

Une phase peut être définie comme une Modèle:Citation Dans cette définition, les propriétés de la phase ne dépendent pas de la position dans la phase, et toute propriété a donc la même valeur en tout point du système. Ce n'est pas le cas dans la réalité. Par exemple, par effet de la gravité, la pression d'une colonne de fluide n'est pas la même en tête et en pied de colonne, ce qu'énonce la loi de l'hydrostatique. Cette différence de pression peut entrainer, surtout pour les fluides compressibles comme les gaz, une importante variation de masse volumique au long de la colonne. Dans une colonne d'eau liquide de Modèle:Unité, par exemple, la différence de pression est d'environ Modèle:Unité entre le pied et la tête de la colonne ; si la pression en tête de colonne est la pression atmosphérique, une bulle d'air présente dans cette colonne a une masse volumique deux fois plus importante en pied qu'en tête. L'hypothèse de l'homogénéité des propriétés est raisonnable sur de petits volumes de phase dans lesquels, par exemple, l'effet de la loi de l'hydrostatique peut être négligée : dans un verre d'eau, la pression en fond de verre peut être considérée comme égale à la pression de surface. Cette hypothèse est généralement retenue en thermodynamique, étant donné que l'équilibre de phases est étudié dans les parages immédiats des interfaces entre ces phases ; de plus, dans le cadre d'une modélisation de la phase, elle permet de simplifier la mise en équation, en évitant de faire dépendre les propriétés du fluide de son emplacement dans l'espace. Toutefois, le système ne peut être trop petit, car à une certaine échelle l'effet de la tension superficielle ne peut plus être négligé et doit être pris en compte.

Une phase peut cependant être dispersée dans une autre phase. Une phase dispersée se présente comme un ensemble de sous-domaines séparés les uns des autres par une autre phase et ayant tous les mêmes propriétés (état, pression, température, composition, forme cristalline, masse volumiqueModèle:Etc.). L'ensemble de ces sous-domaines est considéré comme une phase unique. Par exemple, un précipité ou une suspension est une phase solide divisée en petits éléments dans une phase liquide ; une émulsion est une phase liquide divisée en gouttelettes dans une autre phase liquide ; un aérosol est une phase liquide divisée en gouttelettes dans une phase gazeuse.

Une phase peut aussi être définie comme un Modèle:Citation La continuité signifie qu'il n'y a pas de « saut » de valeur d'une propriété de ce système si l'on modifie les conditions opératoires telles que la pression, la température ou la composition, aussi petite cette modification soit-elle. Un saut de valeur est caractéristique d'une transition de phase (voir la section Classification des transitions de phase). Par exemple, lorsque l'on comprime un gaz à température constante, sa masse volumique augmente continument. Lorsque la pression dépasse la pression de vapeur saturante du gaz, la masse volumique du système augmente fortement et brusquement : la phase n'est plus la même, le système est devenu liquide, de densité beaucoup plus importante que le gaz. De même, un matériau ferromagnétique ou ferrimagnétique dont on augmente la température perd brusquement son aimantation permanente à la température de Curie et devient paramagnétique.

Équilibre et stabilité

Équilibre de phases

• Équilibre liquide-vapeur du brome.
  Équilibre liquide-vapeur du brome.
• Deux phases liquides en équilibre : de l'huile (jaune, en haut) sur de l'eau (transparente, en bas).
  Deux phases liquides en équilibre : de l'huile (jaune, en haut) sur de l'eau (transparente, en bas).

Selon les conditions opératoires, un corps pur ou un mélange peut se présenter dans une phase ou une autre, voire sous plusieurs phases simultanément. La coexistence de plusieurs phases est appelée équilibre de phases. Les phases en équilibre peuvent être dans le même état ou dans des états différents. Lorsqu'une phase liquide est en équilibre avec une phase gazeuse, cette dernière est généralement appelée vapeur, l'équilibre est appelé équilibre liquide-vapeur plutôt que équilibre liquide-gaz. À pression atmosphérique normale (Modèle:Unité), les phases solide (glace) et liquide de l'eau coexistent à Modèle:Unité, ses phases gazeuse (vapeur d'eau) et liquide à Modèle:Unité. Pour un corps pur les trois états coexistent de façon stable au point triple, un unique couple de pression et de température : pour l'eau ce point se situe à Modèle:Unité et Modèle:Unité. Plus généralement, le terme de point triple désigne tout point de coexistence de trois phases, quels que soient leurs états. Le diagramme de phase de l'eau montre plusieurs points triples en plus du point solide-liquide-vapeur : plusieurs points avec le liquide et deux formes de glace, et plusieurs points avec trois formes de glace.

Dans les équilibres de phases d'un mélange, les diverses phases ont le plus souvent des compositions différentes. Par exemple, un mélange d'huile et d'eau présente simultanément une phase liquide d'huile saturée en eau et une phase liquide aqueuse saturée en huile ; de même, un mélange d'eau et d'éthanol présente le plus souvent un équilibre liquide-vapeur dans lequel les deux phases ont des compositions différentes. Il existe des cas particuliers d'équilibre de phases d'un mélange dans lesquels les phases ont la même composition : c'est notamment le cas des azéotropes pour les équilibres liquide-vapeur et celui des points de fusion congruents pour les équilibres liquide-solide.

Stabilité et métastabilité

• Exemple d'état métastable. La perturbation appliquée à cette eau liquide à température négative, en état de surfusion, déclenche son passage à l'état solide.
  Exemple d'état métastable. La perturbation appliquée à cette eau liquide à température négative, en état de surfusion, déclenche son passage à l'état solide.
• Diagramme de phases du carbone montrant les domaines de stabilité et de métastabilité du graphite et du diamant. La graphite est la forme stable dans les conditions atmosphériques.
  Diagramme de phases du carbone montrant les domaines de stabilité et de métastabilité du graphite et du diamant. La graphite est la forme stable dans les conditions atmosphériques.

De façon générale, un système thermodynamique présentant une ou plusieurs phases est dans un état d'équilibre dit stable s'il revient à cet équilibre après avoir subit une perturbation (un apport d'énergie relativement faible)Modèle:,. À pression atmosphérique et Modèle:Unité, si l'on agite brièvement de l'eau à son équilibre solide-liquide, cet équilibre n'est pas modifié (une agitation prolongée peut apporter suffisamment de chaleur au système pour provoquer la fonte de la glace).

Certains états observés se maintiennent dans le temps jusqu'à ce qu'une perturbation déclenche une transformation souvent violente : un tel état est dit métastable. La transformation observée consiste en l'établissement de la ou des phases stables dans les conditions opératoires données. Ainsi, l'eau peut être maintenue liquide en dessous de Modèle:Unité à pression atmosphérique : cet état est appelé état de surfusion, il suffit alors de l'introduction d'une poussière ou d'une fluctuation de température pour que l'eau se solidifie brutalement ; cet état s'observe de façon naturelle lors des pluies verglaçantes. L'état métastable de l'eau non vaporisée au-dessus de Modèle:Unité est appelé surébullition ou retard à l'ébullition ; il est observé dans une tasse de café chauffé au microonde qui entre en ébullition brutale lorsque l'on y ajoute du sucre. Dans les conditions de pression et température atmosphériques, le graphite est la forme solide stable du carbone, le diamant est métastable.

Diagramme de phases

Types de diagramme

• Diagramme de phases du mélange méthane (Modèle:Fchim)-éthane (Modèle:Fchim). Pression de changement d'état en fonction de la température, à compositions fixes.
  Diagramme de phases du mélange méthane (Modèle:Fchim)-éthane (Modèle:Fchim). Pression de changement d'état en fonction de la température, à compositions fixes.
• Diagramme binaire isobare du mélange niobium (Nb)-antimoine (Sb). Température de changement d'état en fonction de la composition, à pression atmosphérique.
  Diagramme binaire isobare du mélange niobium (Nb)-antimoine (Sb). Température de changement d'état en fonction de la composition, à pression atmosphérique.
• Diagramme isotherme du mélange chloroforme-éthanol. Pression de changement d'état en fonction de la composition à Modèle:Unité.
  Diagramme isotherme du mélange chloroforme-éthanol. Pression de changement d'état en fonction de la composition à Modèle:Unité.
• Diagramme du mélange éthanol-acétate d'éthyle. Fraction molaire de l'éthanol en phase vapeur en fonction de la fraction molaire de l'éthanol en phase liquide à l'équilibre, à pression atmosphérique.
  Diagramme du mélange éthanol-acétate d'éthyle. Fraction molaire de l'éthanol en phase vapeur en fonction de la fraction molaire de l'éthanol en phase liquide à l'équilibre, à pression atmosphérique.
• Diagramme ternaire du mélange plomb (Pb)-bismuth (Bi)-étain (Sn).
  Diagramme ternaire du mélange plomb (Pb)-bismuth (Bi)-étain (Sn).

Un diagramme de phases représente les domaines d'existence des phases d'un matériau donné (corps pur ou mélange) en fonction des conditions opératoires qui déterminent l'existence de ces phases : la pression, la température et la composition. Les courbes qui délimitent ces domaines représentent les transitions de phases.

Pour un corps pur, un diagramme de phases représente la pression en fonction de la température. Un corps pur ne présente qu'une seule phase gazeuse et une seule phase liquide. Par contre, il peut présenter plusieurs phases solides.

Pour un mélange, plusieurs types de diagrammes sont souvent utilisés :

• un diagramme représentant la pression en fonction de la température à composition constante ;
• un diagramme isobare, représentant, à pression donnée constante, l'évolution de la température en fonction de la composition ;
• un diagramme isotherme, représentant, à température donnée constante, l'évolution de la pression en fonction de la composition ;
• un diagramme représentant, pour deux phases en équilibre, la composition de l'une des phases en fonction de la composition de l'autre phase ; dans ce diagramme la pression ou la température est fixée.

Les pentes des courbes d'un diagramme représentant la pression en fonction de la température sont données théoriquement par la formule de Clapeyron ou par les formules d'Ehrenfest.

Construction d'un diagramme isobare ou isotherme

Fichier:Diagramme de phase isotherme - Px - principe.jpg

Diagramme de phases isotherme. Principes.

L'exemple pris ici concerne un équilibre liquide-vapeur (gaz), mais peut être transposé à tout autre équilibre de phases. Le diagramme isotherme ci-contre a été construit pour deux corps (espèces chimiques), notés $1$ et $2$, de pressions de vapeur saturante respectives $P_1$ et $P_2$ à la température du diagramme. Avec $P_1<P_2$, le corps $2$ est plus léger que le corps $1$ : il se liquéfie à une pression plus élevée pour une même température.

Un diagramme isotherme représente, à température donnée constante, la pression en fonction de la composition d'un mélange. Un diagramme isobare représente, à pression donnée constante, la température en fonction de la composition d'un mélange. Ces deux types de diagramme sont utilisés pour des mélanges binaires (constitués de deux espèces chimiques). Les fractions molaires de ces deux espèces répondent donc à la contrainte $x_1+x_2=1$. Une courbe de pression ou de température peut ainsi n'être représentée qu'en fonction d'une seule de ces deux fractions : à température donnée $P=P\left(x_2\right)$, à pression donnée $T=T\left(x_2\right)$. Les courbes reportées sur ces diagrammes représentent les conditions opératoires dans lesquelles débutent les transitions de phase, elles délimitent les domaines d'existence des différentes phases du mélange. Un domaine du diagramme dans lequel n'existe qu'une seule phase est délimité par des courbes fonctions de la composition de cette phase. Par exemple, dans le diagramme ci-contre, le domaine d'existence de la phase liquide seule (aux hautes pressions) est délimité par la courbe de bulle (bleue) $P=P\left(x_2^{\text{l}}\right)$ tracée en fonction de la composition du liquide ; cette courbe représente les conditions de transition du liquide vers la vapeur dans lesquelles apparait la première bulle de gaz. Le domaine d'existence de la phase gaz seule (aux basses pressions) est délimité par la courbe de rosée (orange) $P=P\left(x_2^{\text{g}}\right)$ tracée en fonction de la composition du gaz ; cette courbe représente les conditions de transition de la vapeur vers le liquide dans lesquelles apparait la première goutte de liquide.

Entre deux domaines représentant des phases seules, il peut exister un domaine dans lequel les deux phases coexistent en équilibre. Par exemple, dans le diagramme ci-contre, il existe, entre les domaines du liquide seul et du gaz seul, un domaine dans lequel les deux phases coexistent. Dans ce type de domaine, les compositions des phases en équilibre sont données par les intersections des courbes limitrophes et de la droite horizontale isobare ou isotherme (selon le type de diagramme). Dans le diagramme ci-contre, au point $𝖢$ dans le domaine de coexistence des phases liquide et gaz, la composition de la phase liquide est donnée par le point ${𝖢}^{\text{l}}$ sur la courbe de bulle et la composition de la phase vapeur par le point ${𝖢}^{\text{g}}$ sur la courbe de rosée, les trois points étant à la même pression $P_{𝖢}$. Au point $𝖢$ le corps $2$ a une fraction molaire $x_2^{\text{l}}$ en phase liquide et $x_2^{\text{g}}$ en phase gaz. Le point $𝖢$, lui, se situe à une fraction molaire $x_2$ comprise entre $x_2^{\text{l}}$ et $x_2^{\text{g}}$.

On suppose que le mélange est constitué d'une quantité de matière totale $n$, dont une quantité $n_2$ de l'espèce $2$. On a par définition :

$x_2=\frac{n_2}{n}$

Lorsque le mélange est biphasique, la quantité de matière se répartit entre une quantité $n^{\text{l}}$ en phase liquide, dont une quantité $n_2^{\text{l}}$ de l'espèce $2$, et une quantité $n^{\text{g}}$ en phase vapeur, dont une quantité $n_2^{\text{g}}$ de l'espèce $2$. Les bilans de matière donnent :

$n=n^{\text{l}}+n^{\text{g}}$

$n_2=n_2^{\text{l}}+n_2^{\text{g}}$

Les fractions molaires sont définies par :

$x_2^{\text{l}}=\frac{n_2^{\text{l}}}{n^{\text{l}}}$

$x_2^{\text{g}}=\frac{n_2^{\text{g}}}{n^{\text{g}}}$

En considérant les définitions des fractions molaires et les bilans de matière, on montre la règle des moments ou règle du levierModèle:,Modèle:,Modèle:, :

Règle des moments ou règle du levier : $\frac{\left[𝖢{𝖢}^{\text{l}}\right]}{\left[{𝖢}^{\text{g}}{𝖢}^{\text{l}}\right]}=\frac{x_2-x_2^{\text{l}}}{x_2^{\text{g}}-x_2^{\text{l}}}=\frac{n^{\text{g}}}{n^{\text{l}}+n^{\text{g}}}=\frac{n^{\text{g}}}{n}$

Autrement dit, au point $𝖢$, le rapport des longueurs des segments de droite $\left[𝖢{𝖢}^{\text{l}}\right]$ et $\left[{𝖢}^{\text{g}}{𝖢}^{\text{l}}\right]$ permet de déterminer la proportion $n^{\text{g}}/n$ de gaz dans le mélange. Ce paramètre est appelé taux de vaporisation ou fraction de vapeur :

Taux de vaporisation ou fraction de vapeur : ${\tau }^{\text{g}}=\frac{n^{\text{g}}}{n^{\text{l}}+n^{\text{g}}}=\frac{n^{\text{g}}}{n}$

On peut de même définir un taux de liquéfaction ou fraction de liquide : ${\tau }^{\text{l}}=\frac{\left[𝖢{𝖢}^{\text{g}}\right]}{\left[{𝖢}^{\text{l}}{𝖢}^{\text{g}}\right]}=\frac{n^{\text{l}}}{n}=1-{\tau }^{\text{g}}$.

Thermodynamique des équilibres de phases

Conditions d'équilibre

Conditions générales

Soit un système isolé (qui n'échange ni énergie ni matière avec l'extérieur) composé deux de phases, notées $a$ et $b$. Ce système subit une transformation afin d'atteindre l'équilibre. On écrit pour chacune des deux phases la variation d'énergie interne $U$ :

$\mathrm{d}{U}^a=-P^a\mathrm{d}{V}^a+T^a\mathrm{d}{S}^a+\sum _{i=1}^N{\mu }_i^a\mathrm{d}{n}_i^a$

$\mathrm{d}{U}^b=-P^b\mathrm{d}{V}^b+T^b\mathrm{d}{S}^b+\sum _{i=1}^N{\mu }_i^b\mathrm{d}{n}_i^b$

avec

$X^p$

la grandeur

$X$

relative à la phase

$p$

. Soit, pour le système isolé global, les fonctions globales énergie interne globale

$U=U^a+U^b$

, entropie globale

$S=S^a+S^b$

, volume global

$V=V^a+V^b$

et quantité globale

$n_i=n_i^a+n_i^b$

de chacun des

$N$

corps

$i$

. On a les relations sur le système global isolé, en l'absence de réaction chimique :

• $\mathrm{d}U=\mathrm{d}{U}^a+\mathrm{d}{U}^b=0$ (premier principe de la thermodynamique) ;
• $\mathrm{d}S=\mathrm{d}{S}^a+\mathrm{d}{S}^b\ge 0$ (deuxième principe de la thermodynamique) ;
• $\mathrm{d}V=\mathrm{d}{V}^a+\mathrm{d}{V}^b=0$ (absence de travail des forces de pression) ;
• $\mathrm{d}{n}_i=\mathrm{d}{n}_i^a+\mathrm{d}{n}_i^b=0$ pour tout $i\in [1,\cdots ,N]$ (conservation de la matière).

Les grandeurs globales du système global isolé (de l'ensemble constitué par les deux phases) sont constantes (hormis l'entropie), mais les grandeurs propres à chacun des deux phases, elles, peuvent varier au cours de la transformation, leurs variations étant opposées. Par exemple, le volume global du système isolé, $V=V^a+V^b$, est constant, mais les volumes des phases $V^a$ et $V^b$ peuvent varier et $\mathrm{d}{V}^b=-\mathrm{d}{V}^a$.

On obtient les relations globales d'évolution du système isolé :

Évolution $\mathrm{d}U=-(P^a-P^b)\mathrm{d}{V}^a+T^a\mathrm{d}{S}^a+T^b\mathrm{d}{S}^b+\sum _{i=1}^N({\mu }_i^a-{\mu }_i^b)\mathrm{d}{n}_i^a=0$ $\mathrm{d}S=\mathrm{d}{S}^a+\mathrm{d}{S}^b\ge 0$

L'équilibre est atteint lorsque l'entropie globale du système atteint un maximum. Les variations d'entropie des deux phases sont alors liées par : $\mathrm{d}S=\mathrm{d}{S}^a+\mathrm{d}{S}^b=0$. On peut donc écrire pour le système global à l'équilibreModèle:, :

À l'équilibre $\mathrm{d}U=-(P^a-P^b)\mathrm{d}{V}^a+(T^a-T^b)\mathrm{d}{S}^a+\sum _{i=1}^N({\mu }_i^a-{\mu }_i^b)\mathrm{d}{n}_i^a=0$

Quand les deux phases sont à l'équilibre, cela ne signifie pas qu'il n'y a pas d'échanges entre elles. Les deux phases continuent à échanger de la matière et de l'énergie, mais globalement les échanges se font dans les deux sens et se compensent : l'état d'équilibre obtenu est qualifié de dynamique ou stationnaire. Les variations $\mathrm{d}{V}^a$, $\mathrm{d}{S}^a$ et $\mathrm{d}{n}_i^a$ pour tout $i\in [1,\cdots ,N]$ sont donc non nulles à l'équilibre. La relation précédente implique par conséquent l'homogénéité des potentiels mécanique (pression), thermique (température) et chimiques (potentiels chimiques) entre les deux phases à l'équilibre :

• $P^a=P^b$ ;
• $T^a=T^b$ ;
• ${\mu }_i^a={\mu }_i^b$ pour tout $i\in [1,\cdots ,N]$.

Les conditions d'équilibre sont donc données parModèle:,Modèle:, :

Conditions d'équilibre de deux phases $a$ et $b$

Équilibre mécanique	homogénéité de la pression $P$
Équilibre thermique	homogénéité de la température $T$
Équilibre diffusif	homogénéité des potentiels chimiques pour tout corps $i$ : ${\mu }_i^a={\mu }_i^b$

Si le système contient une troisième phase $c$, celle-ci doit être à l'équilibre avec les phases $a$ et $b$ : elle sera par conséquent à l'équilibre aux mêmes pression, température et potentiels chimiques. Ces conditions sont donc généralisables à plus de deux phases en équilibre (équilibres gaz-liquide-liquide ou gaz-liquide-solide par exemple, ou multiples phases liquides ou solides).

Note sur l'homogénéité des potentiels chimiques

L'homogénéité des potentiels chimiques signifie que pour le corps 1 ${\mu }_1^a={\mu }_1^b$, pour le corps 2 ${\mu }_2^a={\mu }_2^b$, … , pour le corps $N$ ${\mu }_N^a={\mu }_N^b$, mais pas que ces potentiels sont égaux entre corps différents : ${\mu }_1^a\ne {\mu }_2^a\ne \cdots \ne {\mu }_N^a$.

Homogénéité des fugacités

Dans n'importe quelle phase $p$, la fugacité $f_i^p$ du corps $i$ est définie par la relation entre potentiels chimiques :

Fugacité ${\mu }_i^p-{\mu }_i^{\bullet ,*}=RT\ln \left(\frac{f_i^p}{P}\right)$

avec : Modèle:Colonnes

La fugacité $f_i^p$ lie les propriétés d'un corps $i$, pur ou en mélange, dans une phase $p$ quelconque réelle, aux propriétés de ce même corps à l'état de gaz parfait pur aux mêmes pression et température que la phase réelle. On peut écrire pour deux phases $a$ et $b$ en équilibre, en introduisant les fugacités respectives $f_i^a$ et $f_i^b$ de tout corps $i$ :

${\mu }_i^a={\mu }_i^{\bullet ,*}+RT\ln \left(\frac{f_i^a}{P}\right)$

${\mu }_i^b={\mu }_i^{\bullet ,*}+RT\ln \left(\frac{f_i^b}{P}\right)$

À l'équilibre, les deux phases $a$ et $b$ sont aux mêmes pression et température, en conséquence le potentiel ${\mu }_i^{\bullet ,*}$ est le même dans ces deux relations, et l'on peut écrire, en considérant que ${\mu }_i^a={\mu }_i^b$Modèle:,Modèle:, :

Condition d'équilibre des phases : $f_i^a=f_i^b$

En termes de condition d'équilibre des phases, l'égalité des fugacités est donc équivalente à l'égalité des potentiels chimiques : si l'une est vérifiée, alors l'autre est aussi vérifiée, et réciproquement. Néanmoins, les potentiels chimiques ne pouvant être déterminés qu'à une constante additive près, l'égalité des potentiels chimiques est difficilement utilisable, tandis que les fugacités peuvent être déterminées de façon absolue, ce qui les rend particulièrement utiles.

Règle des phases

Erreur lors de la création de la vignette :

Diagramme de phase de l'eau.
Les courbes représentent les équilibres de deux phases : si la température est fixée, la pression est subie, et réciproquement. Dans les domaines délimités par les courbes une seule phase est stable : on peut y faire varier la pression indépendamment de la température. Le point critique et le point triple ne peuvent être modifiés.

Soit un mélange quelconque composé de $N$ espèces chimiques réparties en $p$ phases. Les conditions d'équilibre stable d'un système hétérogène définies plus haut impliquent que :

• les $p$ pressions sont liées par les $p-1$ égalités :

$P^a=P^b=\cdots =P^p$

• les $p$ températures sont liées par les $p-1$ égalités :

$T^a=T^b=\cdots =T^p$

• pour chacun des $N$ constituants, les $p$ potentiels chimiques sont liés par les $p-1$ égalités :

$\{\begin{array}{c}{\mu }_1^a={\mu }_1^b=\cdots ={\mu }_1^p\\ {\mu }_2^a={\mu }_2^b=\cdots ={\mu }_2^p\\ \cdots \\ {\mu }_N^a={\mu }_N^b=\cdots ={\mu }_N^p\end{array}$

Pour chaque phase, les diverses variables sont également liées par la relation de Gibbs-Duhem, soit $p$ équations supplémentaires :

$\{\begin{array}{c}{V}^a\mathrm{d}{P}^a-S^a\mathrm{d}{T}^a=\sum _{i=1}^N{n}_i^a\mathrm{d}{\mu }_i^a\\ V^b\mathrm{d}{P}^b-S^b\mathrm{d}{T}^b=\sum _{i=1}^N{n}_i^b\mathrm{d}{\mu }_i^b\\ \cdots \\ V^p\mathrm{d}{P}^p-S^p\mathrm{d}{T}^p=\sum _{i=1}^N{n}_i^p\mathrm{d}{\mu }_i^p\end{array}$

On a donc un total de $X=(2+N)p$ variables et $Y=(2+N)(p-1)+p$ équations. La variance $v=X-Y$ indique le nombre de variables, c'est-à-dire le nombre de conditions opératoires (parmi la pression, la température et les compositions des phases), que l'opérateur peut fixer indépendamment l'une de l'autre en restant dans le domaine d'équilibre stable des phases, d'où la règle des phasesModèle:, :

Règle des phases : $v=N+2-p$

Par exemple, pour un corps pur, soit $N=1$ :

• en présence d'une seule phase (par exemple un liquide seul), soit $p=1$, on a $v=2$ : l'opérateur peut imposer la pression et la température indépendamment l'une de l'autre ;
• en présence de deux phases (par exemple un équilibre liquide-gaz), soit $p=2$, on a $v=1$ : l'opérateur peut imposer la pression mais la température d'équilibre est subie (et inversement) ;
• en présence de trois phases (par exemple le point triple de l'équilibre solide-liquide-gaz), soit $p=3$, on a $v=0$ : l'opérateur ne peut imposer aucune condition de l'équilibre.

Un mélange de $N$ constituants peut au plus se présenter sous $p=N+2$ phases à l'équilibre stable. Il ne peut y avoir au plus qu'une seule phase gaz, mais il peut y avoir plusieurs phases liquides et solides à l'équilibre.

Théorème de Gibbs-Konovalov

• Erreur lors de la création de la vignette :
  Diagramme de phases isobare d'un azéotrope négatif. Au maximum des courbes de température, les compositions des deux phases sont identiques.
• Diagramme de phases du mélange chlorure d'hydrogène (Modèle:Fchim) - eau (Modèle:Fchim) montrant quatre points de fusion congruents à Modèle:Unité, Modèle:Unité, Modèle:Unité et Modèle:Unité.
  Diagramme de phases du mélange chlorure d'hydrogène (Modèle:Fchim) - eau (Modèle:Fchim) montrant quatre points de fusion congruents à Modèle:Unité, Modèle:Unité, Modèle:Unité et Modèle:Unité.

Dans un diagramme de phases isobare ou isotherme d'un mélange, les courbes de transition de phase peuvent présenter des extrémums (maximum ou minimum), respectivement de pression ou de température, en fonction de la composition. Le théorème de Gibbs-Konovalov établit queModèle:, :

Modèle:Début citation blocAux points d'extrémum de pression ou de température, les phases en équilibre ont la même composition.Modèle:Fin citation bloc

Dans le diagramme de phases d'un équilibre liquide-vapeur, par exemple, si la courbe de bulle atteint un extrémum, alors la courbe de rosée atteint un extrémum à la même composition. Pour la pression et la température considérées, les phases liquide et vapeur en équilibre ont la même composition, et le mélange se comporte comme un corps pur : il s'agit d'un azéotrope. Dans un équilibre liquide-solide, ce cas de figure est appelé point de fusion congruent.

Stabilité d'un équilibre

Le deuxième principe de la thermodynamique induit qu'un système à l'équilibre n'est stable que si l'entropie $S$ est maximale. Ceci induit que, pour un système à pression et température données, l'enthalpie libre $G$ est minimale. Ce principe implique également que l'entropie est une fonction concave de ses variables, et l'enthalpie libre, à pression et température données, une fonction convexe de ses variables.

L'entropie dans un changement d'état

On apporte, à pression $P$ et quantité de matière $n$ constantes, de la chaleur $Q$ à un corps pur dans un récipient fermé de volume $V$ pouvant varier. L'énergie interne $U$ du corps pur, en présence du travail des seules forces de pression, évolue selon le premier principe de la thermodynamique :

$\mathrm{d}U=-P\mathrm{d}V+\delta Q_{P,n}$

En introduisant l'enthalpie $H=U+PV$, on a, à pression $P$ constante :

$\delta Q_{P,n}=\mathrm{d}U+P\mathrm{d}V=\mathrm{d}[U+PV]=\mathrm{d}H$

La transformation est réversible, on introduit l'entropie $S$ :

$\delta Q_{P,n}=T\mathrm{d}S$

La chaleur apportée au corps pur est donc liée à l'entropie et à l'enthalpie selon :

$\delta Q_{P,n}=T\mathrm{d}S=\mathrm{d}H$

Puisque la température thermodynamique $T$ est positive, on obtient :

Modèle:Equarefa : ${\left(\frac{\partial S}{\partial H}\right)}_{P,n}=\frac{1}{T}>0$

La transformation ayant lieu à pression et quantité de matière constantes, on a par définition de la capacité thermique isobare $C_P$ :

$C_P=T{\left(\frac{\partial S}{\partial T}\right)}_{P,n}={\left(\frac{\partial H}{\partial S}\right)}_{P,n}{\left(\frac{\partial S}{\partial T}\right)}_{P,n}={\left(\frac{\partial H}{\partial T}\right)}_{P,n}$

On peut donc écrire la dérivée seconde :

${\left(\frac{{\partial }^{2}S}{{\partial H}^2}\right)}_{P,n}={\left(\frac{\partial }{\partial H}{\left(\frac{\partial S}{\partial H}\right)}_{P,n}\right)}_{P,n}={\left(\frac{\partial \frac{1}{T}}{\partial H}\right)}_{P,n}=-\frac{1}{T^2}{\left(\frac{\partial T}{\partial H}\right)}_{P,n}=-\frac{1}{C_P{T}^2}$

Les conditions de stabilité thermodynamique donnent $C_P\ge 0$. On obtient :

Modèle:Equarefa : ${\left(\frac{{\partial }^{2}S}{{\partial H}^2}\right)}_{P,n}=-\frac{1}{C_P{T}^2}\le 0$

Les équations Modèle:Equarefl et Modèle:Equarefl montrent que, à pression et composition constantes, l'entropie $S$ est une fonction croissante concave de l'enthalpie $H$.

Fichier:Entropie d'un corps pur - changement d'état.jpg

Changement d'état d'un corps pur. La courbe bleue ABC et la courbe orange CDE sont concaves. Cependant entre B et D l'entropie ABCDE est convexe. Pour conserver la concavité de l'entropie, le corps pur change d'état et l'entropie passe par le segment gris BD.

Sur la figure ci-contre, la courbe ABC bleue représente l'évolution, à pression et quantité de matière constantes, de l'entropie d'un corps pur dans une Modèle:Nobr, la courbe CDE orange celle du même corps pur dans une Modèle:Nobr. Les deux courbes sont concaves : si l'on joint deux points bleus, ou deux points orange, par un segment de droite, celui-ci est en-dessous de la courbe considérée. Cependant, la courbe globale ABCDE n'est pas dérivable au point C. De plus, elle est localement convexe entre B et D : si l'on joint un point de BC à un point de CD par un segment de droite, celui-ci est au-dessus de la courbe. Entre les points B et D, le corps pur ne peut être en équilibre stable s'il suit la courbe BCD. Pour conserver ses propriétés de continuité, dérivabilité et concavité, l'entropie passe par le segment de droite BD gris. Au point B le segment gris est tangent à la courbe bleue, au point D il est tangent à la courbe orange. La courbe ABDE est globalement dérivable et concave, contrairement à la courbe ABCDE. Le segment BD est la seule trajectoire qui rend l'entropie globalement dérivable et concave de A à E, tout autre point que B ou D, toute autre courbe qu'un segment de droite entre ces deux points engendre une courbe d'entropie localement non dérivable ou convexe.

Aux basses enthalpies, en amont du point B, le corps pur est exclusivement en Modèle:Nobr. En aval du point D, aux hautes enthalpies, le corps pur est exclusivement en Modèle:Nobr. Entre B et D, pour une enthalpie donnée, le segment gris représente une entropie supérieure aux entropies des Modèle:Nobr et $b$ seules représentées respectivement par les courbes bleue et orange : cette entropie maximale dans ces conditions de pression et température représente le point de stabilité du système, mais ne représente pas une phase seule. Entre B et D, le corps pur subit un changement d'état et les deux Modèle:Nobr et $b$, à enthalpie donnée, coexistent en un équilibre stable.

La relation $\mathrm{d}H=T\mathrm{d}S$ induit que sur le segment BD la température est constante : à la pression donnée, les deux phases du corps pur ne peuvent coexister qu'à une seule température, conformément à la règle des phases. Ainsi, les grandeurs $S^{𝖡}$, $H^{𝖡}$, $S^{𝖣}$ et $H^{𝖣}$ sont établies aux mêmes pression $P$, température $T_{a\to b}$ et quantité de matière $n$. La température de changement d'état est donnée par le rapport :

Température de changement d'état : $T_{a\to b}=\frac{H^{𝖣}-H^{𝖡}}{S^{𝖣}-S^{𝖡}}$

L'écart $S^{𝖣}-S^{𝖡}$ rapporté à la quantité de matière totale $n$ mise en jeu est l'entropie de changement d'état du corps pur. L'écart $H^{𝖣}-H^{𝖡}$ rapporté à la quantité de matière totale $n$ mise en jeu est l'enthalpie de changement d'état du corps pur. L'enthalpie de changement d'état est la chaleur qu'il faut apporter à une mole de corps pur pour passer entièrement d'un état à l'autre :

Entropie de changement d'état : ${\mathrm{\Delta }}_{a\to b}S=\frac{S^{𝖣}-S^{𝖡}}{n}$ Enthalpie de changement d'état : ${\mathrm{\Delta }}_{a\to b}H=\frac{H^{𝖣}-H^{𝖡}}{n}=T_{a\to b}{\mathrm{\Delta }}_{a\to b}S$

L'équation du segment BD gris s'écrit :

$S=(1-{\tau }^b)\times S^{𝖡}+{\tau }^b\times S^{𝖣}$

$H=(1-{\tau }^b)\times H^{𝖡}+{\tau }^b\times H^{𝖣}$

Les points B et D représentant la même quantité $n$ de corps pur, on peut écrire, en introduisant les grandeurs molaires respectives :

$S=(1-{\tau }^b)n\times {\overline{S}}^{𝖡}+{\tau }^{b}n\times {\overline{S}}^{𝖣}$

$H=(1-{\tau }^b)n\times {\overline{H}}^{𝖡}+{\tau }^{b}n\times {\overline{H}}^{𝖣}$

Pour ${\tau }^b=0$ le corps pur est au point B, entièrement en Modèle:Nobr ; pour ${\tau }^b=1$ il est au point D, entièrement en Modèle:Nobr. Pour $0<{\tau }^b<1$ le corps pur est biphasique, entre B et D : une quantité $n^a=(1-{\tau }^b)n$ de Modèle:Nobr coexiste en un équilibre stable avec une quantité $n^b={\tau }^{b}n$ de Modèle:Nobr. Le paramètre ${\tau }^b$ est le taux de changement d'état, ou fraction de Modèle:Nobr, c'est-à-dire la quantité de Modèle:Nobr rapportée à la quantité de matière totale $n$ mise en jeu lors de la transformation :

Taux de changement d'état ou fraction de phase : ${\tau }^b=\frac{n^b}{n^a+n^b}=\frac{n^b}{n}$

Autrement dit, si l'on dispose, à la pression $P$ et la température $T_{a\to b}$, d'une part d'une quantité $n$ de corps pur en Modèle:Nobr et que l'on en prélève une fraction $1-{\tau }^b$, et d'autre part d'une quantité $n$ de corps pur en Modèle:Nobr et que l'on en prélève une fraction ${\tau }^b$, si l'on met en contact ces deux fractions, alors celles-ci sont d'emblée en équilibre stable, sans nécessité d'une transformation pour atteindre cet équilibre.

Les branches BC bleue et CD orange sont des branches de métastabilité des phases correspondantes. Si le corps pur évolue sur ces branches (y compris au-delà de C pour la Modèle:Nobr et en deçà de C pour la Modèle:Nobr), la moindre perturbation peut provoquer un brusque changement d'état, amenant l'entropie sur la courbe ABDE stable, à enthalpie constante. L'état métastable d'un liquide non solidifié en dessous de sa température de fusion est appelé surfusion, il s'observe de façon naturelle pour l'eau lors des pluies verglaçantes (voir également la vidéo présentée dans la section Stabilité et métastabilité). L'état métastable d'un liquide non vaporisé au-dessus de sa température d'ébullition est appelé surébullition ou retard à l'ébullition, il s'observe dans une tasse de café chauffé au microonde qui entre en ébullition brutale lorsque l'on y ajoute du sucre.

Le même raisonnement peut être tenu pour les changements d'état des mélanges ; il doit cependant être tenu compte de la variation des compositions des phases en cours de transition.

Équilibre de phases selon la composition

On considère un mélange constitué de deux espèces chimiques notées $1$ et $2$. Si l'on fait varier la pression, la température et la composition du mélange, l'énergie interne du système varie selon :

$\mathrm{d}U=-P\mathrm{d}V+T\mathrm{d}S+{\mu }_1\mathrm{d}{n}_1+{\mu }_2\mathrm{d}{n}_2$

avec :

• $n_1$ et $n_2$ les quantités totales respectives de $1$ et $2$ dans le mélange ;
• ${\mu }_1$ et ${\mu }_2$ les potentiels chimiques respectifs des deux corps.

Par la suite, on considère une transformation à pression et température données, seule la variation de la composition est étudiée. À pression et température constantes, on introduit l'enthalpie libre $G=U+PV-TS$ et la relation précédente devient :

$\mathrm{d}U+P\mathrm{d}V-T\mathrm{d}S=\mathrm{d}[U+PV-TS]=\mathrm{d}G={\mu }_1\mathrm{d}{n}_1+{\mu }_2\mathrm{d}{n}_2$

Le théorème d'Euler permet d'écrire, en introduisant l'enthalpie libre molaire $\overline{G}$ et les fractions molaires $x_1$ et $x_2$ :

$G=n_1{\mu }_1+n_2{\mu }_2$

$\overline{G}=x_1{\mu }_1+x_2{\mu }_2$

À pression et température constantes, le deuxième principe de la thermodynamique induit que l'enthalpie libre ne peut que décroître lors d'une transformation spontanée ; à l'équilibre, l'enthalpie libre atteint un minimum. Lorsque l'on mélange deux corps purs, l'enthalpie libre du mélange résultant ne peut donc être, à l'équilibre, qu'inférieure à la somme des enthalpies libres des corps purs considérés séparément, soit :

$\overline{G}\le x_1{\mu }_1^{*}+x_2{\mu }_2^{*}$

avec ${\mu }_1^{*}$ et ${\mu }_2^{*}$ les potentiels chimiques des deux corps purs. Un mélange tel que $\overline{G}>x_1{\mu }_1^{*}+x_2{\mu }_2^{*}$ est instable. Une grandeur molaire n'est fonction que de la pression, de la température et des fractions molaires : $\overline{G}=\overline{G}(P,T,x_1,x_2)$. Pour un mélange binaire, les deux fractions étant liées par la relation $x_1+x_2=1$, on peut écrire la fonction comme ne dépendant, arbitrairement, que de $x_2$ pour la composition : $\overline{G}=\overline{G}(P,T,x_2)$. Ceci permet de tracer des diagrammes à pression et température constantes comme celui de la figure ci-contre. À pression et température données, l'enthalpie libre est une fonction continue, dérivable et convexe selon la composition pour tout $0\le x_2\le 1$.

Fichier:Enthalpie libre mélange binaire 2 phases.jpg

Enthalpie libre d'un mélange binaire à pression et température constantes. Stabilité et équilibre de phases en fonction de la composition. Une phase est stable si son enthalpie libre est minimale. Entre B et D, l'enthalpie libre est minimale pour un équilibre des deux phases.

On suppose que le mélange binaire peut se présenter, aux pression et température données, sous une phase $a$ ou une phase $b$ selon la composition du mélange. Soient ${\overline{G}}^a$ et ${\overline{G}}^b$ leurs enthalpies libres molaires respectives :

${\overline{G}}^a=x_1{\mu }_1^a+x_2{\mu }_2^a$

${\overline{G}}^b=x_1{\mu }_1^b+x_2{\mu }_2^b$

La figure ci-contre donne un exemple de diagramme représentant simultanément deux courbes ${\overline{G}}^a$ et ${\overline{G}}^b$ représentatives de la plupart des fonctions de ce type. Ces deux courbes sont continues, dérivables et convexes pour toutes les compositions.

Pour une composition donnée (une valeur de $x_2$), la phase stable est la phase ayant une enthalpie libre minimale. Pour les faibles valeurs de $x_2$, entre les points A et B, on a ${\overline{G}}^a<{\overline{G}}^b$, c'est donc la phase $a$ qui est stable et $\overline{G}={\overline{G}}^a$. Pour les fortes valeurs de $x_2$, entre les points D et E, on a ${\overline{G}}^b<{\overline{G}}^a$, c'est donc la phase $b$ qui est stable et $\overline{G}={\overline{G}}^b$.

Entre les points B et D, si l'on suit le même raisonnement, la fonction $\overline{G}$ devrait passer par le point C et ne serait plus dérivable ni convexe. La fonction passe par le segment gris BD, tangent aux deux courbes en ces points : ce segment est la seule fonction qui rend l'enthalpie libre continue, dérivable et convexe sur toutes les compositions. Entre B et D les deux phases $a$ et $b$ sont en équilibreModèle:,. Par exemple, pour la composition $x_2$ donnée sur le graphe, on a ${\overline{G}}^b<{\overline{G}}^a$, mais la phase $b$ est métastable : le mélange n'est stable que sous la forme d'un équilibre des deux phases représenté par le point P sur le segment gris, auquel $\overline{G}<{\overline{G}}^b$. À cet équilibre, la phase $a$ a la composition $x_2^{𝖡}$ du point B et la phase $b$ la composition $x_2^{𝖣}$ du point D. Avec la composition $x_2$, la règle des moments donne la fraction de phaseModèle:, :

Fraction de phase $b$ : ${\tau }^b=\frac{n^b}{n^a+n^b}=\frac{n^b}{n_1+n_2}=\frac{x_2-x_2^{𝖡}}{x_2^{𝖣}-x_2^{𝖡}}$

avec $n^a$ et $n^b$ les quantités respectives des phases $a$ et $b$. Entre les points B et D on a par conséquentModèle:, :

$\overline{G}=(1-{\tau }^b)\times {\overline{G}}^{𝖡}+{\tau }^b\times {\overline{G}}^{𝖣}$

avec les enthalpies libres respectives aux points B et D :

${\overline{G}}^{𝖡}={\overline{G}}^a\left(x_2^{𝖡}\right)=(1-x_2^{𝖡}){\mu }_1^a\left(x_2^{𝖡}\right)+x_2^{𝖡}{\mu }_2^a\left(x_2^{𝖡}\right)$

${\overline{G}}^{𝖣}={\overline{G}}^b\left(x_2^{𝖣}\right)=(1-x_2^{𝖣}){\mu }_1^b\left(x_2^{𝖣}\right)+x_2^{𝖣}{\mu }_2^b\left(x_2^{𝖣}\right)$

La continuité et la dérivabilité de la fonction $\overline{G}$ aux points B et D induisent qu'à l'équilibre de deux phases les potentiels chimiques d'une espèce chimique dans les deux phases sont égaux :

${\mu }_1^a\left(x_2^{𝖡}\right)={\mu }_1^b\left(x_2^{𝖣}\right)$ ${\mu }_2^a\left(x_2^{𝖡}\right)={\mu }_2^b\left(x_2^{𝖣}\right)$

Par définition, le potentiel chimique est l'enthalpie libre molaire partielle. Pour l'espèce $1$ dans le mélange binaire, on a :

${\mu }_1={\left(\frac{\partial G}{\partial n_1}\right)}_{P,T,n_2}={\left(\frac{\partial n\overline{G}}{\partial n_1}\right)}_{P,T,n_2}=\overline{G}+n{\left(\frac{\partial \overline{G}}{\partial n_1}\right)}_{P,T,n_2}$

avec $n=n_1+n_2$ la quantité de matière totale du mélange, $x_1=n_1/n$ et $x_2=n_2/n$. La relation $x_1+x_2=1$ permettant de ne retenir, arbitrairement, que la seule variable de composition $x_2$, on a $\overline{G}=\overline{G}(P,T,x_2)$. La dérivation à pression et température constantes donne :

${\left(\frac{\partial \overline{G}}{\partial n_1}\right)}_{P,T,n_2}={\left(\frac{\partial \overline{G}}{\partial x_2}\right)}_{P,T}{\left(\frac{\partial x_2}{\partial n_1}\right)}_{n_2}={\left(\frac{\partial \overline{G}}{\partial x_2}\right)}_{P,T}(-\frac{x_2}{n})$

${\mu }_1=\overline{G}-{\left(\frac{\partial \overline{G}}{\partial x_2}\right)}_{P,T}{x}_2$

On a pour l'espèce $2$ :

${\mu }_2={\left(\frac{\partial G}{\partial n_2}\right)}_{P,T,n_1}=\overline{G}+{\left(\frac{\partial \overline{G}}{\partial x_2}\right)}_{P,T}(1-x_2)$

On obtient la dérivée de la fonction $\overline{G}$ :

${\left(\frac{\partial \overline{G}}{\partial x_2}\right)}_{P,T}={\mu }_2-{\mu }_1$

Au point B, la fonction ${\overline{G}}^a$ donne :

${\left(\frac{\partial {\overline{G}}^a}{\partial x_2}\right)}_{P,T}^{𝖡}={\mu }_2^a\left(x_2^{𝖡}\right)-{\mu }_1^a\left(x_2^{𝖡}\right)$

Au point D, la fonction ${\overline{G}}^b$ donne :

${\left(\frac{\partial {\overline{G}}^b}{\partial x_2}\right)}_{P,T}^{𝖣}={\mu }_2^b\left(x_2^{𝖣}\right)-{\mu }_1^b\left(x_2^{𝖣}\right)$

La fonction $\overline{G}$ étant dérivable pour tout $0\le x_2\le 1$, aux points B et D les deux dérivées précédentes sont égales à la pente du segment BD :

${\left(\frac{\partial {\overline{G}}^a}{\partial x_2}\right)}_{P,T}^{𝖡}={\left(\frac{\partial {\overline{G}}^b}{\partial x_2}\right)}_{P,T}^{𝖣}=\frac{{\overline{G}}^{𝖣}-{\overline{G}}^{𝖡}}{x_2^{𝖣}-x_2^{𝖡}}$

On pose :

$\mathrm{\Gamma }={\mu }_2^a\left(x_2^{𝖡}\right)-{\mu }_1^a\left(x_2^{𝖡}\right)={\mu }_2^b\left(x_2^{𝖣}\right)-{\mu }_1^b\left(x_2^{𝖣}\right)=\frac{[(1-x_2^{𝖣}){\mu }_1^b\left(x_2^{𝖣}\right)+x_2^{𝖣}{\mu }_2^b\left(x_2^{𝖣}\right)]-[(1-x_2^{𝖡}){\mu }_1^a\left(x_2^{𝖡}\right)+x_2^{𝖡}{\mu }_2^a\left(x_2^{𝖡}\right)]}{x_2^{𝖣}-x_2^{𝖡}}$

On a, en réarrangeant le dernier terme :

$\begin{array}{cc}\mathrm{\Gamma }& =\frac{{\mu }_1^b\left(x_2^{𝖣}\right)-{\mu }_1^a\left(x_2^{𝖡}\right)}{x_2^{𝖣}-x_2^{𝖡}}+\frac{x_2^{𝖣}[{\mu }_2^b\left(x_2^{𝖣}\right)-{\mu }_1^b\left(x_2^{𝖣}\right)]-x_2^{𝖡}[{\mu }_2^a\left(x_2^{𝖡}\right)-{\mu }_1^a\left(x_2^{𝖡}\right)]}{x_2^{𝖣}-x_2^{𝖡}}\\ & =\frac{{\mu }_1^b\left(x_2^{𝖣}\right)-{\mu }_1^a\left(x_2^{𝖡}\right)}{x_2^{𝖣}-x_2^{𝖡}}+\frac{x_2^{𝖣}\mathrm{\Gamma }-x_2^{𝖡}\mathrm{\Gamma }}{x_2^{𝖣}-x_2^{𝖡}}\\ & =\frac{{\mu }_1^b\left(x_2^{𝖣}\right)-{\mu }_1^a\left(x_2^{𝖡}\right)}{x_2^{𝖣}-x_2^{𝖡}}+\mathrm{\Gamma }\end{array}$

$\frac{{\mu }_1^b\left(x_2^{𝖣}\right)-{\mu }_1^a\left(x_2^{𝖡}\right)}{x_2^{𝖣}-x_2^{𝖡}}=0$

On en déduit que :

${\mu }_1^a\left(x_2^{𝖡}\right)={\mu }_1^b\left(x_2^{𝖣}\right)$

et, puisque $\mathrm{\Gamma }={\mu }_2^a\left(x_2^{𝖡}\right)-{\mu }_1^a\left(x_2^{𝖡}\right)={\mu }_2^b\left(x_2^{𝖣}\right)-{\mu }_1^b\left(x_2^{𝖣}\right)$ :

${\mu }_2^a\left(x_2^{𝖡}\right)={\mu }_2^b\left(x_2^{𝖣}\right)$

On considère $x_2$ comme la seule variable de composition de l'enthalpie libre molaire $\overline{G}=\overline{G}(P,T,x_2)$, puisque l'on a $x_1+x_2=1$. On a les dérivées :

${\left(\frac{\partial \overline{G}}{\partial x_2}\right)}_{P,T}={\mu }_2-{\mu }_1$

${\left(\frac{{\partial }^2\overline{G}}{{\partial x_2}^2}\right)}_{P,T}={\left(\frac{\partial {\mu }_2}{\partial x_2}\right)}_{P,T}-{\left(\frac{\partial {\mu }_1}{\partial x_2}\right)}_{P,T}$

Le théorème de dérivation des fonctions composées donne pour ${\mu }_2$ (qui, comme $\overline{G}$, est considéré comme ne dépendant que de $x_2$ en composition) :

${\left(\frac{\partial {\mu }_2}{\partial n_2}\right)}_{P,T,n_1}={\left(\frac{\partial {\mu }_2}{\partial x_2}\right)}_{P,T}{\left(\frac{\partial x_2}{\partial n_2}\right)}_{n_1}={\left(\frac{\partial {\mu }_2}{\partial x_2}\right)}_{P,T}(\frac{1}{n}-\frac{x_2}{n})$

${\left(\frac{\partial {\mu }_2}{\partial x_2}\right)}_{P,T}=\frac{n}{x_1}{\left(\frac{\partial {\mu }_2}{\partial n_2}\right)}_{P,T,n_1}$

On obtient de la même façon pour ${\mu }_1$ :

$-{\left(\frac{\partial {\mu }_1}{\partial x_2}\right)}_{P,T}=\frac{n}{x_2}{\left(\frac{\partial {\mu }_1}{\partial n_1}\right)}_{P,T,n_2}$

La dérivée seconde de $\overline{G}$ se réécrit :

${\left(\frac{{\partial }^2\overline{G}}{{\partial x_2}^2}\right)}_{P,T}=\frac{n}{x_1}{\left(\frac{\partial {\mu }_2}{\partial n_2}\right)}_{P,T,n_1}+\frac{n}{x_2}{\left(\frac{\partial {\mu }_1}{\partial n_1}\right)}_{P,T,n_2}$

Les conditions de stabilité thermodynamique imposent queModèle:, :

${\left(\frac{\partial {\mu }_1}{\partial n_1}\right)}_{P,T,n_2}\ge 0$

${\left(\frac{\partial {\mu }_2}{\partial n_2}\right)}_{P,T,n_1}\ge 0$

On a donc :

${\left(\frac{{\partial }^2\overline{G}}{{\partial x_2}^2}\right)}_{P,T}\ge 0$

L'enthalpie libre molaire $\overline{G}=\overline{G}(P,T,x_2)$, à pression et température constantes, est convexe par rapport à $x_2$.

Exemple de construction d'un diagramme de phases

Soit un mélange de deux corps, notés $1$ et $2$, à une température $T$ donnée. Seules les variations de sa composition (des fractions molaires $x_1$ et $x_2$) et de la pression $P$ sont étudiées. On suppose qu'en fonction de la pression et de la composition le mélange peut être liquide ou gazeux. Soient les enthalpies libres molaires des deux phases :

pour le liquide : $\frac{{\overline{G}}^{\text{l}}}{RT}=x_1\ln x_1+x_2\ln x_2+x_1{x}_2(A_{2,1}{x}_1+A_{1,2}{x}_2)$ pour le gaz : $\frac{{\overline{G}}^{\bullet }}{RT}=x_1\ln \left(\frac{x_{1}P}{P_1^{\text{sat}}}\right)+x_2\ln \left(\frac{x_{2}P}{P_2^{\text{sat}}}\right)$

Soient ${\mu }_1^{\text{l,*}}$ et ${\mu }_2^{\text{l,*}}$ les potentiels chimiques des deux corps purs à l'état liquide à $P$. L'enthalpie libre molaire ${\overline{G}}^{\text{l}}$ du mélange liquide à $P$ est calculée selon :

$\frac{{\overline{G}}^{\text{l}}}{RT}=\frac{x_1{\mu }_1^{\text{l,*}}+x_2{\mu }_2^{\text{l,*}}}{RT}+x_1\ln x_1+x_2\ln x_2+x_1{x}_2(A_{2,1}{x}_1+A_{1,2}{x}_2)$

L'enthalpie libre de mélange molaire est calculée, pour les besoins de l'exemple, en employant le modèle de Margules à deux paramètres, $A_{1,2}$ et $A_{2,1}$, pour le terme d'excès. Ce modèle ne dépend pas de la pression.

La phase gazeuse est considérée comme un mélange de gaz parfaits. Soient ${\mu }_1^{\bullet ,*}$ et ${\mu }_2^{\bullet ,*}$ les potentiels chimiques des deux corps purs à l'état de gaz parfait à $P$. L'enthalpie libre molaire ${\overline{G}}^{\bullet }$ du mélange gazeux à $P$ est calculée selon :

$\frac{{\overline{G}}^{\bullet }}{RT}=\frac{x_1{\mu }_1^{\bullet ,*}+x_2{\mu }_2^{\bullet ,*}}{RT}+x_1\ln x_1+x_2\ln x_2$

Soit ${\mu }_i^{\bullet ,*,\text{sat}}$ le potentiel chimique du Modèle:Nobr à sa pression de vapeur saturante $P_i^{\text{sat}}$ à la température $T$. On a la relation :

${\mu }_i^{\bullet ,*}={\mu }_i^{\bullet ,*,\text{sat}}+RT\ln \frac{P}{P_i^{\text{sat}}}$

À la température $T$ et la pression $P_i^{\text{sat}}$, le Modèle:Nobr pur est à l'équilibre liquide-vapeur. Son potentiel en phase gaz (considérée comme un gaz parfait) est donc égal à son potentiel en phase liquide ${\mu }_i^{\text{l,*,sat}}$ dans ces conditions :

${\mu }_i^{\bullet ,*,\text{sat}}={\mu }_i^{\text{l,*,sat}}$

Un liquide étant peu compressible, son potentiel chimique peut être considéré comme indépendant de la pression, d'où :

${\mu }_i^{\text{l,*}}\approx {\mu }_i^{\text{l,*,sat}}$

On a par conséquent :

$\frac{{\overline{G}}^{\bullet }}{RT}\approx \frac{x_1{\mu }_1^{\text{l,*}}+x_2{\mu }_2^{\text{l,*}}}{RT}+x_1\ln \left(\frac{x_{1}P}{P_1^{\text{sat}}}\right)+x_2\ln \left(\frac{x_{2}P}{P_2^{\text{sat}}}\right)$

Les deux potentiels ${\mu }_1^{\text{l,*}}\approx {\mu }_1^{\text{l,*,sat}}$ et ${\mu }_2^{\text{l,*}}\approx {\mu }_2^{\text{l,*,sat}}$ peuvent être considérés comme des constantes, puisque l'on travaille à température $T$ constante. Arbitrairement, on pose Modèle:Nobr.

La figure 1 donne les tracés de ces deux fonctions avec, à la température de Modèle:Unité :

• Modèle:Nobr : butan-1-ol, Modèle:Nobr ;
• Modèle:Nobr : eau, Modèle:Nobr ;
• $A_{1,2}=\ln \left(21\right)$ et $A_{2,1}=\ln \left(4\right)$.

${\overline{G}}^{\text{l}}$ ne dépendant pas de la pression, la figure 1 ne donne qu'une seule courbe d'enthalpie libre pour le liquide. Les points 1 et 2 représentent respectivement le butan-1-ol pur à sa pression de vapeur saturante $P_1^{\text{sat}}$ et l'eau à $P_2^{\text{sat}}$. Par hypothèse, ces points représentent également ces corps liquides à des pressions plus élevées. Plusieurs tracés de ${\overline{G}}^{\bullet }$ sont donnés en fonction de la pression.

• Figure 1 - Enthalpies libres du liquide et du gaz du mélange butan-1-ol - eau calculées avec le modèle d'activité de Margules[42].
  Figure 1 - Enthalpies libres du liquide et du gaz du mélange butan-1-ol - eau calculées avec le modèle d'activité de Margules.
• Figure 2 - Diagramme de phases du mélange butan-1-ol - eau construit avec le modèle de Margules[42]. Le point J est un hétéroazéotrope.
  Figure 2 - Diagramme de phases du mélange butan-1-ol - eau construit avec le modèle de Margules. Le point J est un hétéroazéotrope.

À Modèle:Unité, on a ${\overline{G}}^{\text{l}}<{\overline{G}}^{\bullet }$ à toute composition : le mélange est exclusivement liquide. Le mélange est dans une phase liquide riche en butan-1-ol et pauvre en eau (Modèle:Nobr) jusqu'au point I, dans une phase riche en eau et pauvre en butan-1-ol (Modèle:Nobr) au-delà du point K. Ces deux phases liquides sont en équilibre entre les points I et K : il s'agit d'une démixtion. Ceci est valable à partir de Modèle:Unité.

À Modèle:Unité, on a ${\overline{G}}^{\text{l}}<{\overline{G}}^{\bullet }$ aux faibles fractions $x_2$, en amont du point I : le mélange est en Modèle:Nobr. On a ${\overline{G}}^{\text{l}}<{\overline{G}}^{\bullet }$ aux fortes fractions $x_2$, en aval du point K : le mélange est en Modèle:Nobr. Au point J la courbe ${\overline{G}}^{\bullet }$ est tangente au segment IK : le mélange est triphasique, le gaz est en équilibre avec les Modèle:Nobr et Modèle:Nobr, il s'agit d'un hétéroazéotrope. Entre les points I et J le mélange est biphasique : la Modèle:Nobr est en équilibre avec le gaz. Entre les points J et K le mélange est biphasique : le gaz est en équilibre avec la Modèle:Nobr.

À Modèle:Unité, on a ${\overline{G}}^{\text{l}}<{\overline{G}}^{\bullet }$ aux faibles fractions $x_2$, en amont du point E : le mélange est en Modèle:Nobr. On a ${\overline{G}}^{\bullet }<{\overline{G}}^{\text{l}}$ entre les points F et G : le mélange est gazeux. On a ${\overline{G}}^{\text{l}}<{\overline{G}}^{\bullet }$ aux fortes fractions $x_2$, en aval du point H : le mélange est en Modèle:Nobr. Entre les points E et F le mélange est biphasique : la Modèle:Nobr est en équilibre avec le gaz. Entre les points G et H le mélange est biphasique : le gaz est en équilibre avec la Modèle:Nobr.

À Modèle:Unité, on a ${\overline{G}}^{\text{l}}<{\overline{G}}^{\bullet }$ aux faibles fractions $x_2$, en amont du point C : le mélange est en Modèle:Nobr. On a ${\overline{G}}^{\bullet }<{\overline{G}}^{\text{l}}$ entre les points D et 2 : le mélange est gazeux. On a ${\overline{G}}^{\text{l}}={\overline{G}}^{\bullet }$ pour $x_2=1$ : l'eau est pure en équilibre liquide-vapeur au point 2. Entre les points C et D le mélange est biphasique : la Modèle:Nobr est en équilibre avec le gaz.

À Modèle:Unité, on a ${\overline{G}}^{\text{l}}<{\overline{G}}^{\bullet }$ aux faibles fractions $x_2$, en amont du point A : le mélange est en Modèle:Nobr. On a ${\overline{G}}^{\bullet }<{\overline{G}}^{\text{l}}$ aux fortes fractions $x_2$, en aval du point B : le mélange est gazeux. Entre les points A et B le mélange est biphasique : la Modèle:Nobr est en équilibre avec le gaz.

À Modèle:Unité, on a ${\overline{G}}^{\text{l}}={\overline{G}}^{\bullet }$ pour $x_2=0$ : le butan-1-ol est pur en équilibre liquide-vapeur au point 1. On a ${\overline{G}}^{\bullet }<{\overline{G}}^{\text{l}}$ aux autres fractions $x_2$, en aval du point 1 : le mélange est gazeux.

À Modèle:Unité, on a ${\overline{G}}^{\bullet }<{\overline{G}}^{\text{l}}$ à toute composition : le mélange est exclusivement gazeux.

La figure 2 donne le diagramme de phases issu de cette analyse.

Classification des transitions de phase

Classification d'Ehrenfest

Fichier:Phase transition first and second order de.svg

Classification d'Ehrenfest des transitions de phase. Modèle:1re ligne : absence de transition. Modèle:2e ligne : transition d'Modèle:Nobr. Modèle:3e ligne : transition d'Modèle:Nobr.

La classification d'Ehrenfest des transitions de phase repose sur le principe suivantModèle:,Modèle:,Modèle:,Modèle:, :

Modèle:Début citation blocUne transition de phase est d'Modèle:Nobr si la fonction enthalpie libre et ses dérivées jusqu'à l'Modèle:Nobr sont continues, tandis qu'une de ses dérivées d'Modèle:Nobr au moins est discontinue.Modèle:Fin citation bloc

La classification d'Ehrenfest est aujourd'hui abandonnée au profit d'une classification qui en reprend partiellement le principe ainsi que d'autres énoncés par Landau. Elle permet néanmoins de définir les principales transitions de phase étudiées en thermodynamique, notamment les changements d'état.

Invariance de l'enthalpie libreModèle:,Modèle:,Modèle:,

Dans ces conditions, pour toute transition d'une Modèle:Nobr à une Modèle:Nobr, l'enthalpie libre $G$ est continue, c'est-à-dire invariante lors de cette transition :

$G^b=G^a$

Ceci revient à l'égalité des potentiels chimiques des deux phases :

${\mu }^b=\frac{G^b}{n}=\frac{G^a}{n}={\mu }^a$

Transition d'Modèle:NobrModèle:,Modèle:,Modèle:,

Les dérivées partielles d'Modèle:Nobr de l'enthalpie libre sont :

${\left(\frac{\partial G}{\partial P}\right)}_T=V$

${\left(\frac{\partial G}{\partial T}\right)}_P=-S$

Dans une transition d'Modèle:Nobr les dérivées d'Modèle:Nobr sont discontinues, d'où :

$V^b\ne V^a$

$S^b\ne S^a$

avec :

• $V^a$ et $V^b$ les volumes des Modèle:Nobr et $b$ ;
• $S^a$ et $S^b$ les entropies des Modèle:Nobr et $b$.

Dans une transition d'Modèle:Nobr le système étudié change de volume. De plus, il y a une entropie de changement d'état ${\mathrm{\Delta }}_{a\to b}S$ et donc une enthalpie de changement d'état ${\mathrm{\Delta }}_{a\to b}H$ :

${\mathrm{\Delta }}_{a\to b}S=\frac{S^b-S^a}{n}$

${\mathrm{\Delta }}_{a\to b}H=T{\mathrm{\Delta }}_{a\to b}S$

Les transitions d'Modèle:Nobr correspondent aux changements d'état, qui nécessitent un transfert thermique (échange de chaleur) entre la matière étudiée et l'extérieur. Ce type de transition s'effectue par une succession continue d'équilibres des deux phases, l'une disparaissant progressivement au profit de l'autre. Si l'apport de chaleur est interrompu au cours de la transformation, alors celle-ci s'interrompt et les deux phases présentes subsistent sans autre modificationModèle:,Modèle:,Modèle:,Modèle:,Modèle:,.

Par exemple, si l'on place de l'eau liquide dans un récipient fermé pouvant changer de volume et que l'on chauffe à pression atmosphérique, l'eau se met à bouillir à Modèle:Unité et se vaporise progressivement. Si l'on coupe la chauffe au cours de cette transition et que l'on maintient le récipient à pression atmosphérique et Modèle:Unité, alors on observe deux états (vapeur et liquide) coexistant de façon stable tant que l'on ne rétablit pas la chauffe. Lorsque la transition est totalement accomplie, le volume résultant de gaz est bien supérieur au volume initial de liquide.

Une transition de phase d'ordre supérieur s'effectue sans enthalpie de changement d'état. Toutes les transitions de phase sans changement d'état sont donc des transitions de phase d'ordre supérieur à 1Modèle:,Modèle:,Modèle:,Modèle:,.

Transition d'Modèle:NobrModèle:,Modèle:,Modèle:,

Dans une transition d'Modèle:Nobr les dérivées d'Modèle:Nobr sont continues :

$V^b=V^a$

$S^b=S^a$

Les dérivées partielles d'Modèle:Nobr de l'enthalpie libre sont :

${\left(\frac{{\partial }^{2}G}{\partial T^2}\right)}_P=-{\left(\frac{\partial S}{\partial T}\right)}_P=-\frac{C_P}{T}$

$\left(\frac{{\partial }^{2}G}{\partial P\partial T}\right)=-{\left(\frac{\partial S}{\partial P}\right)}_T=\left(\frac{{\partial }^{2}G}{\partial T\partial P}\right)={\left(\frac{\partial V}{\partial T}\right)}_P=V\alpha$

${\left(\frac{{\partial }^{2}G}{\partial P^2}\right)}_T={\left(\frac{\partial V}{\partial P}\right)}_T=-V{\chi }_T$

La température $T$ et le volume $V=V^b=V^a$ sont communs aux deux phases. Dans une transition d'Modèle:Nobr les dérivées d'Modèle:Nobr sont discontinues, d'où :

$C_P^b\ne C_P^a$

${\alpha }^b\ne {\alpha }^a$

${\chi }_T^b\ne {\chi }_T^a$

avec :

• ${\alpha }^a$ et ${\alpha }^b$ les coefficients de dilatation isobare des Modèle:Nobr et $b$ ; $\alpha =\frac{1}{V}{\left(\frac{\partial V}{\partial T}\right)}_P$ ;
• ${\chi }_T^a$ et ${\chi }_T^b$ les coefficients de compressibilité isotherme des Modèle:Nobr et $b$ ; ${\chi }_T=-\frac{1}{V}{\left(\frac{\partial V}{\partial P}\right)}_T$ ;
• $C_P^a$ et $C_P^b$ les capacités thermiques isobares des Modèle:Nobr et $b$ ; $C_P=T{\left(\frac{\partial S}{\partial T}\right)}_P$.

Dans une transition de phase d'Modèle:Nobr, le corps pur ne change pas de volume, la transformation ne nécessite pas un apport de chaleur et les deux phases ne coexistent pas. Un exemple de transition d'Modèle:Nobr est la transition conducteur-supraconducteur d'un métalModèle:,Modèle:,. La transition liquide-gaz au point critique est également une transition d'Modèle:NobrModèle:,.

Transitions d'ordre supérieur

La condensation de Bose-Einstein d'un gaz parfait formé de bosons est un exemple de transition d'Modèle:Nobr. La transition Berezinski-Kosterlitz-Thouless est un cas extrême de transition d'ordre infini.

Formules de Clapeyron et d'Ehrenfest

Soit une transition de phase d'une phase notée $a$ à une autre notée $b$ s'effectuant à la température $T$ donnée et à la pression constante $P_{a\to b}$. La variation de la pression $P_{a\to b}$ en fonction de la température $T$ peut être déterminée théoriquement. Les formules suivantes donnent les pentes des courbes délimitant les domaines représentés dans les diagrammes de phases représentant la pression en fonction de la température.

Pour un changement d'état (une transition de phase d'Modèle:Nobr selon la classification d'Ehrenfest des transitions de phase), la pression de changement d'état $P_{a\to b}$ d'un corps pur varie en fonction de la température $T$ selon la formule de ClapeyronModèle:, :

Formule de Clapeyron : $\left(\frac{\mathrm{d}{P}_{a\to b}}{\mathrm{d}T}\right)=\frac{{\mathrm{\Delta }}_{a\to b}H}{T\cdot ({\overline{V}}^b-{\overline{V}}^a)}$ pour un corps pur.

La formule de Clapeyron étendue aux mélanges, pour un mélange de $N$ constituants notés $i$, pour une Modèle:Nobr de composition constante $x^a$ en équilibre avec une Modèle:Nobr de composition variable $x^b$, s'écritModèle:,Modèle:, :

Formule de Clapeyron étendue aux mélanges : ${\left(\frac{\partial P_{a\to b}}{\partial T}\right)}_{x^a}=\frac{{\overline{H}}^b-\sum _{i=1}^N{x}_i^b{\overline{H}}_i^a}{T\cdot ({\overline{V}}^b-\sum _{i=1}^N{x}_i^b{\overline{V}}_i^a)}$

Pour une transition de phase d'Modèle:Nobr selon la classification d'Ehrenfest, la pression de transition de phase $P_{a\to b}$ varie en fonction de la température $T$ selon les formules d'EhrenfestModèle:, :

Formules d'Ehrenfest : $\left(\frac{\mathrm{d}{P}_{a\to b}}{\mathrm{d}T}\right)=\frac{{\alpha }^b-{\alpha }^a}{{\chi }_T^b-{\chi }_T^a}$ $\left(\frac{\mathrm{d}{P}_{a\to b}}{\mathrm{d}T}\right)=\frac{{\overline{C}}_P^b-{\overline{C}}_P^a}{T\overline{V}\cdot ({\alpha }^b-{\alpha }^a)}$

avec : Modèle:Colonnes

Calcul d'un équilibre liquide-vapeur

Le calcul d'équilibre liquide-vapeur est le plus documenté dans la littérature. Son principe peut être décliné à tout autre équilibre de phases, notamment aux équilibres vapeur-liquide-liquide et vapeur-liquide-solide.

Énoncé du problème

Mise en équations

Soit, à pression $P$ et température $T$ données, un mélange de $N$ constituants, chaque constituant $i$ étant représenté par la quantité de matière $n_i$. On cherche à déterminer la composition des phases liquide et vapeur en présence. Il s'agit 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 ;

qui doivent répondre aux équations :

• $n_i=n_i^{\text{l}}+n_i^{\text{g}}$ du bilan de matière ;
• $f_i^{\text{g}}=f_i^{\text{l}}$ de l'équilibre liquide-vapeur.

On a donc $2N$ inconnues et $2N$ équations.

Reformulation

On définit les variables intermédiaires :

• $n^{\text{g}}=\sum _{i=1}^N{n}_i^{\text{g}}$ la quantité de matière totale en phase vapeur ;
• $n^{\text{l}}=\sum _{i=1}^N{n}_i^{\text{l}}$ la quantité de matière totale en phase liquide ;
• $n=\sum _{i=1}^N{n}_i$ la quantité de matière totale considérée ;
• $x_i^{\text{g}}=\frac{n_i^{\text{g}}}{n^{\text{g}}}$ la fraction molaire du constituant $i$ en phase vapeur ;
• $x_i^{\text{l}}=\frac{n_i^{\text{l}}}{n^{\text{l}}}$ la fraction molaire du constituant $i$ en phase liquide.

Soit pour chacun des $N$ constituants la fraction molaire globale, notée $x_i$ :

Fraction molaire globale du corps $i$ : $x_i=\frac{n_i}{n}$

On définit enfin ${\tau }^{\text{g}}$ la fraction molaire globale de la phase vapeur, ou taux de vaporisation (inconnue) :

Taux de vaporisation ou fraction de vapeur : ${\tau }^{\text{g}}=\frac{n^{\text{g}}}{n^{\text{g}}+n^{\text{l}}}$

Les fractions $x_i$ sont calculées directement à partir des données du problème. Les variables $n_i^{\text{g}}$, $n^{\text{g}}$, $x_i^{\text{g}}$, $n_i^{\text{l}}$, $n^{\text{l}}$, $x_i^{\text{l}}$ et ${\tau }^{\text{g}}$ sont des inconnues. Le problème est reformulé de façon à n'avoir à déterminer que les fractions $x_i^{\text{g}}$ et $x_i^{\text{l}}$, et le taux de vaporisation ${\tau }^{\text{g}}$, soit $2N+1$ inconnues. Les $2N+1$ équations résolues sontModèle:, :

• les $N$ bilans de matière $x_i=x_i^{\text{l}}(1-{\tau }^{\text{g}})+x_i^{\text{g}}{\tau }^{\text{g}}$, reformulés selon les fractions molaires ;
• les $N$ équilibres $x_i^{\text{g}}=K_i{x}_i^{\text{l}}$, reformulés à l'aide des coefficients de partage ;
• la contrainte $\sum _{i=1}^N{x}_i^{\text{g}}-\sum _{i=1}^N{x}_i^{\text{l}}=0$.

Le bilan de matière de chacun des $N$ constituants donne :

$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}}(n^{\text{l}}+n^{\text{g}}-n^{\text{g}})+x_i^{\text{g}}{n}^{\text{g}}$

En divisant par la quantité de matière totale $n=n^{\text{g}}+n^{\text{l}}$ on obtient :

$x_i=x_i^{\text{l}}(1-{\tau }^{\text{g}})+x_i^{\text{g}}{\tau }^{\text{g}}$

Les $N$ bilans de matière sont ainsi réécrits selon les fractions molaires. En sommant sur l'ensemble des constituants on a :

$\sum _{i=1}^N{x}_i=1=\left(\sum _{i=1}^N{x}_i^{\text{l}}\right)(1-{\tau }^{\text{g}})+\left(\sum _{i=1}^N{x}_i^{\text{g}}\right){\tau }^{\text{g}}$

Si l'on a $\sum _{i=1}^N{x}_i^{\text{g}}=\sum _{i=1}^N{x}_i^{\text{l}}$, alors :

$1=\left(\sum _{i=1}^N{x}_i^{\text{l}}\right)(1-{\tau }^{\text{g}})+\left(\sum _{i=1}^N{x}_i^{\text{l}}\right){\tau }^{\text{g}}=\sum _{i=1}^N{x}_i^{\text{l}}$

On vérifie donc $\sum _{i=1}^N{x}_i^{\text{g}}=\sum _{i=1}^N{x}_i^{\text{l}}=1$.

Les quantités des phases et des constituants peuvent être déduites à postériori selon :

$n^{\text{g}}={\tau }^{\text{g}}n$

$n^{\text{l}}=(1-{\tau }^{\text{g}})n$

$n_i^{\text{g}}=x_i^{\text{g}}{n}^{\text{g}}$

$n_i^{\text{l}}=x_i^{\text{l}}{n}^{\text{l}}$

Fugacités et coefficient de partage

Les fugacités sont des fonctions de la pression $P$ et de la température $T$ du mélange à l'équilibre, et de la composition $x_i$ de la phase considérée :

$f_i^{\text{g}}=f_i^{\text{g}}(P,T,x_1^{\text{g}},\cdots ,x_N^{\text{g}})$

$f_i^{\text{l}}=f_i^{\text{l}}(P,T,x_1^{\text{l}},\cdots ,x_N^{\text{l}})$

Les fugacités sont le plus souvent calculées selon deux approches : celle par coefficient de fugacité ${\varphi }_i$, nécessitant une équation d'état, et celle par coefficient d'activité ${\gamma }_i$, nécessitant un modèle de coefficient d'activité. Ces deux approches sont en théorie équivalentes, mais en pratique leur usage dépend de la pression considérée. L'approche par coefficient de fugacité est employée à toutes les pressions pour la phase gaz et à des pressions supérieures à Modèle:Unité pour la phase liquide. L'usage d'une même équation d'état pour représenter les deux phases (par exemple une équation d'état cubique comme celle de van der Waals) rend le calcul cohérent à l'approche du point critique, où les deux phases doivent se rejoindre. Par contre, aux basses pressions, les équations d'état, spécialement développées pour les gaz, représentent souvent assez mal les liquides. Les modèles de coefficient d'activité sont spécialement développés pour les liquides aux basses pressions (par exemple les modèles de Margules, UNIQUAC, UNIFAC, COSMOSPACE). Ils ne dépendent généralement pas de la pression, mais seulement de la température et de la composition.

L'approche par coefficient de fugacité ${\varphi }_i$ donneModèle:, :

Pour les espèces en phase gaz : $f_i^{\text{g}}=x_i^{\text{g}}{\varphi }_i^{\text{g}}P$ quelle que soit la pression. Pour les espèces en phase liquide : $f_i^{\text{l}}=x_i^{\text{l}}{\varphi }_i^{\text{l}}P$ aux hautes pressions.

L'approche par coefficient d'activité ${\gamma }_i$ donne :

Pour les espèces en phase liquide : $f_i^{\text{l}}=x_i^{\text{l}}{\gamma }_i^{\text{l}}{f}_i^{\text{l},*}$ aux basses pressions.

Cette approche nécessite de connaitre la fugacité $f_i^{\text{l},*}$ de l'espèce chimique $i$ à l'état de corps pur liquide dans les conditions de pression $P$ et température $T$ du mélange à l'équilibre.

Pour une espèce liquide dans les conditions du mélange, pour laquelle il existe une pression de vapeur saturante $P_i^{\text{sat}}$ à la température $T$ telle que $P>P_i^{\text{sat}}$, cette fugacité se calcule selon :

Pour les espèces liquides : $f_i^{\text{l},*}={\varphi }_i^{\text{g,*,sat}}{P}_i^{\text{sat}}{𝒫}_i$

avec :

• $P_i^{\text{sat}}$ la pression de vapeur saturante du corps $i$ pur à $T$ ;
• ${\varphi }_i^{\text{g,*,sat}}$ le coefficient de fugacité du corps $i$ pur à saturation en phase vapeur, à $P_i^{\text{sat}}$ et $T$ ;
• ${\overline{V}}_i^{\text{l},*}$ le volume molaire liquide du corps $i$ pur à $T$ ; dépend également de la pression ;
• ${𝒫}_i=\exp \left(\frac{\underset{P_i^{\text{sat}}}{\overset{P}{\int }}{\overline{V}}_i^{\text{l},*}\mathrm{d}P}{RT}\right)$ la correction de Poynting.

Pour une espèce supercritique dans les conditions du mélange, pour laquelle il n'existe pas de pression de vapeur saturante à la température du mélange (c'est-à-dire $T>T_{\text{c},i}$, la température critique de l'espèce), ou gazeuse, pour laquelle il existe une pression de vapeur saturante telle que $P<P_i^{\text{sat}}$, la fugacité $f_i^{\text{l},*}$ est fictive et peut être calculée par extrapolation à l'aide d'une constante de Henry $k_{\text{H},i}^{\text{l}}$ (déterminée expérimentalement) selon la relation :

Pour les espèces supercritiques ou gazeuses : $f_i^{\text{l},*}=\frac{k_{\text{H},i}^{\text{l}}}{{\gamma }_i^{\text{l},\mathrm{\infty }}}$

avec le coefficient d'activité à dilution infinie :

${\gamma }_i^{\text{l},\mathrm{\infty }}=\underset{x_i^{\text{l}}\to 0}{\lim }{\gamma }_i^{\text{l}}$

Pour une espèce gazeuse dans les conditions de pression et température étudiées, soit $P<P_i^{\text{sat}}$, mais pour laquelle la pression est proche de la pression de vapeur saturante, soit $P\approx P_i^{\text{sat}}$, l'approche pour les liquides peut être appliquée par extrapolation de la correction de Poynting en considérant, puisque les liquides sont peu compressibles, le volume molaire ${\overline{V}}_i^{\text{l},*}$ comme indépendant de la pression : on a alors ${𝒫}_i=\exp ({\overline{V}}_i^{\text{l},*}(P-P_i^{\text{sat}})/\left(RT\right))$. L'approche par constante de Henry est employée dans le cas $P\ll P_i^{\text{sat}}$.

Les approches par coefficient de fugacité et par coefficient d'activité étant théoriquement équivalentes à toutes les pressions, on a, pour les approches aux basses pressions :

${\varphi }_i^{\text{l}}={\gamma }_i^{\text{l}}\frac{f_i^{\text{l},*}}{P}$

Pour tout corps $i$, quelle que soit la pression, on définit le coefficient de partage $K_i$ à partir de l'égalité des fugacités $f_i^{\text{g}}=x_i^{\text{g}}{\varphi }_i^{\text{g}}P=f_i^{\text{l}}=x_i^{\text{l}}{\varphi }_i^{\text{l}}P$ :

Coefficient de partage : $K_i=\frac{x_i^{\text{g}}}{x_i^{\text{l}}}=\frac{{\varphi }_i^{\text{l}}}{{\varphi }_i^{\text{g}}}$

Les coefficients de partage $K_i$ sont des inconnues, ils dépendent de la pression, de la température et des compositions des deux phases :

$K_i=K_i(P,T,x_1^{\text{g}},\cdots ,x_N^{\text{g}},x_1^{\text{l}},\cdots ,x_N^{\text{l}})$

Dans le système d'équations décrivant l'équilibre de phases, les égalités des fugacités $f_i^{\text{g}}=f_i^{\text{l}}$ sont remplacées par les équations équivalentes $x_i^{\text{g}}=K_i{x}_i^{\text{l}}$.

Calcul d'un point de bulle

Un point de bulle est un couple de pression et température auquel apparait la première bulle de vapeur dans un liquide lors d'une vaporisation.

Calcul de la pression de bulle

Dans ce cas, la température $T$ est connue et le mélange est entièrement liquide, soit ${\tau }^{\text{g}}=0$. Le problème consiste à déterminer la pression $P$ à laquelle apparait la première bulle de gaz et la composition de celui-ci, les fractions $x_i^{\text{g}}$, soit $N+1$ inconnues. On a directement :

$n^{\text{l}}=n$

$n^{\text{g}}=0$

$x_i^{\text{l}}=x_i$

$n_i^{\text{l}}=n_i$

$n_i^{\text{g}}=0$

Les $N+1$ équations sont :

• les $N$ équilibres $x_i^{\text{g}}=K_i{x}_i^{\text{l}}$ ;
• la contrainte $\sum _{i=1}^N{x}_i^{\text{g}}=1$.

L'intégration des équilibres dans la contrainte donne l'équation à résoudre :

Équation du point de bulle : $\sum _{i=1}^N{K}_i{x}_i=1$

L'égalité des fugacités $f_i^{\text{g}}=f_i^{\text{l}}$ donne $x_i^{\text{g}}{\varphi }_i^{\text{g}}P=x_i^{\text{l}}{\varphi }_i^{\text{l}}P$, d'où $x_i^{\text{g}}P=K_i{x}_i^{\text{l}}P=K_i{x}_{i}P$ et, par sommation sur les $N$ constituants :

$\left(\sum _{i=1}^N{x}_i^{\text{g}}\right)P=P=\left(\sum _{i=1}^N{K}_i{x}_i\right)P$

Cette équation peut être utilisée pour résoudre le problème par la méthode itérative du point fixe. Connaissant la pression $P_k$ à l'itération $k$, la pression $P_{k+1}$ à l'itération $k+1$ est calculée selon :

$P_{k+1}={\left(\sum _{i=1}^N{K}_i{x}_i\right)}_k{P}_k$

À chaque itération, on calcule la somme :

${\mathrm{\Sigma }}^{\text{g}}=\sum _{i=1}^N{K}_i{x}_i$

Cette somme est censée valoir ${\mathrm{\Sigma }}^{\text{g}}=\sum _{i=1}^N{x}_i^{\text{g}}=1$, ce qui n'est vrai que lorsque le calcul a convergé. Pour pouvoir calculer les coefficients $K_i$, il est nécessaire de normer les fractions molaires au cours des itérations selon :

${\tilde{x}}_i^{\text{g}}=\frac{K_i{x}_i}{{\mathrm{\Sigma }}^{\text{g}}}$

Lorsque ${\mathrm{\Sigma }}^{\text{g}}=1$, on a bien ${\tilde{x}}_i^{\text{g}}=x_i^{\text{g}}$.

Le calcul est initialisé par le choix d'une pression $P$ et de fractions molaires $x_i^{\text{g}}$, et peut être effectué selon l'algorithme suivant ($k$ est l'itération en cours).

Algorithme de calcul de la pression de bulle

1. Calculer les coefficients de partage $K_i$ par les coefficients de fugacité et d'activité.
2. Calculer la somme ${\mathrm{\Sigma }}^{\text{g}}$.
3. Calculer les fractions molaires normées ${\tilde{x}}_i^{\text{g}}$.
4. Calculer une nouvelle pression $P_{k+1}={\mathrm{\Sigma }}^{\text{g}}{P}_k$.
5. Si la valeur absolue $|P_{k+1}-P_k|>ϵ$, une erreur acceptable, recommencer en 1 avec la nouvelle pression. Sinon, la pression de bulle est trouvée. FIN.

Calcul de la température de bulle

Dans ce cas, la pression $P$ est connue. Le problème consiste à déterminer la température $T$ à laquelle apparait la première bulle de gaz et la composition de celui-ci, soit les $N$ fractions $x_i^{\text{g}}$. Les équations sont les mêmes que pour le calcul de la pression de bulle. On peut effectuer le calcul en faisant une hypothèse de température, en calculant la pression de bulle correspondante selon l'algorithme exposé dans la section précédente, puis en modifiant la température jusqu'à ce que la pression calculée corresponde à la pression $P$ donnée.

Calcul d'un point de rosée

Un point de rosée est un couple de pression et température auquel apparait la première goutte de liquide dans un gaz lors d'une liquéfaction.

Calcul de la pression de rosée

Dans ce cas, la température $T$ est connue et le mélange est entièrement gazeux, soit ${\tau }^{\text{g}}=1$. Le problème consiste à déterminer la pression $P$ à laquelle apparait la première goutte de liquide et la composition de celui-ci, les fractions $x_i^{\text{l}}$, soit $N+1$ inconnues. On a directement :

$n^{\text{g}}=n$

$n^{\text{l}}=0$

$x_i^{\text{g}}=x_i$

$n_i^{\text{g}}=n_i$

$n_i^{\text{l}}=0$

Les $N+1$ équations sont :

• les $N$ équilibres $x_i^{\text{g}}=K_i{x}_i^{\text{l}}$ ;
• la contrainte $\sum _{i=1}^N{x}_i^{\text{l}}=1$.

L'intégration des équilibres dans la contrainte donne l'équation à résoudre :

Équation du point de rosée : $\sum _{i=1}^N\frac{x_i}{K_i}=1$

L'égalité des fugacités $f_i^{\text{g}}=f_i^{\text{l}}$ donne $x_i^{\text{g}}{\varphi }_i^{\text{g}}P=x_i^{\text{l}}{\varphi }_i^{\text{l}}P$, d'où $\frac{x_i^{\text{g}}}{K_i}P=\frac{x_i}{K_i}P=x_i^{\text{l}}P$ et, par sommation sur les $N$ constituants :

$\left(\sum _{i=1}^N\frac{x_i}{K_i}\right)P=\left(\sum _{i=1}^N{x}_i^{\text{l}}\right)P=P$

Cette équation peut être utilisée pour résoudre le problème par la méthode itérative du point fixe. Connaissant la pression $P_k$ à l'itération $k$, la pression $P_{k+1}$ à l'itération $k+1$ est calculée selon :

$P_{k+1}={\left(\sum _{i=1}^N\frac{x_i}{K_i}\right)}_k{P}_k$

À chaque itération, on calcule la somme :

${\mathrm{\Sigma }}^{\text{l}}=\sum _{i=1}^N\frac{x_i}{K_i}$

Cette somme est censée valoir ${\mathrm{\Sigma }}^{\text{l}}=\sum _{i=1}^N{x}_i^{\text{l}}=1$, ce qui n'est vrai que lorsque le calcul a convergé. Pour pouvoir calculer les coefficients $K_i$, il est nécessaire de normer les fractions molaires au cours des itérations selon :

${\tilde{x}}_i^{\text{l}}=\frac{x_i}{K_i{\mathrm{\Sigma }}^{\text{l}}}$

Lorsque ${\mathrm{\Sigma }}^{\text{l}}=1$, on a bien ${\tilde{x}}_i^{\text{l}}=x_i^{\text{l}}$.

Le calcul est initialisé par le choix d'une pression $P$ et de fractions molaires $x_i^{\text{l}}$, et peut être effectué selon l'algorithme suivant ($k$ est l'itération en cours).

Algorithme de calcul de la pression de rosée

1. Calculer les coefficients de partage $K_i$ par les coefficients de fugacité et d'activité.
2. Calculer la somme ${\mathrm{\Sigma }}^{\text{l}}$.
3. Calculer les fractions molaires normées ${\tilde{x}}_i^{\text{l}}$.
4. Calculer une nouvelle pression $P_{k+1}={\mathrm{\Sigma }}^{\text{l}}{P}_k$.
5. Si la valeur absolue $|P_{k+1}-P_k|>ϵ$, une erreur acceptable, recommencer en 1 avec la nouvelle pression. Sinon, la pression de rosée est trouvée. FIN.

Calcul de la température de rosée

Dans ce cas, la pression $P$ est connue. Le problème consiste à déterminer la température $T$ à laquelle apparait la première goutte de liquide et la composition de celui-ci, soit les $N$ fractions $x_i^{\text{l}}$. Les équations sont les mêmes que pour le calcul de la pression de rosée. On peut effectuer le calcul en faisant une hypothèse de température, en calculant la pression de rosée correspondante selon l'algorithme exposé dans la section précédente, puis en modifiant la température jusqu'à ce que la pression calculée corresponde à la pression $P$ donnée.

Calcul d'une vaporisation partielle

On connait la pression $P$, la température $T$ et les fractions molaires globales $x_i$. Le taux de vaporisation ${\tau }^{\text{g}}$ et les fractions dans les phases $x_i^{\text{g}}$ et $x_i^{\text{l}}$ sont inconnues, soit $2N+1$ inconnues.

En substituant l'équilibre dans le bilan de matière de chacun des $N$ constituants, on obtient :

$x_i=x_i^{\text{l}}(1-{\tau }^{\text{g}})+x_i^{\text{g}}{\tau }^{\text{g}}=(1+(K_i-1){\tau }^{\text{g}})x_i^{\text{l}}$

Si $1+(K_i-1){\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 obtient $2N$ équations combinant les bilans de matière et les équilibres :

• $x_i^{\text{l}}=\frac{x_i}{1+(K_i-1){\tau }^{\text{g}}}$ ;
• $x_i^{\text{g}}=\frac{K_i{x}_i}{1+(K_i-1){\tau }^{\text{g}}}$.

Les contraintes combinées sur les fractions molaires $x_i^{\text{g}}$ et $x_i^{\text{l}}$ donnent la dernière équation :

$\sum _{i=1}^N{x}_i^{\text{g}}-\sum _{i=1}^N{x}_i^{\text{l}}=\sum _{i=1}^N\frac{(K_i-1)x_i}{1+(K_i-1){\tau }^{\text{g}}}=0$

Le problème est réduit à trouver la solution de cette équation, due à Rachford et Rice (1952)Modèle:,Modèle:,Modèle:, :

Équation de Rachford-Rice : $f\left({\tau }^{\text{g}}\right)=\sum _{i=1}^N\frac{(K_i-1)x_i}{1+(K_i-1){\tau }^{\text{g}}}=0$

La dérivée $f^{\prime }$ de cette fonction s'écrit :

$f^{\prime }\left({\tau }^{\text{g}}\right)=-\sum _{i=1}^N\frac{{(K_i-1)}^2{x}_i}{{(1+(K_i-1){\tau }^{\text{g}})}^2}$

On a $f^{\prime }\left({\tau }^{\text{g}}\right)<0$ : cette fonction est monotone strictement décroissante en ${\tau }^{\text{g}}$. Si l'on trouve sur la plage $0<{\tau }^{\text{g}}<1$ une solution de l'équation $f\left({\tau }^{\text{g}}\right)=0$, cette solution est unique.

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}}=\frac{n^{\text{l}}}{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\frac{(K_i-1)x_i}{K_i-(K_i-1){\tau }^{\text{l}}}=0$

Cette fonction est monotone strictement croissante en ${\tau }^{\text{l}}$.

On peut résoudre l'équation de Rachford-Rice sur l'intervalle $0\le {\tau }^{\text{g}}\le 1$ par la méthode de Newton selon ($j$ est l'itération en cours) :

${\tau }_{j+1}^{\text{g}}={\tau }_j^{\text{g}}-\frac{f\left({\tau }_j^{\text{g}}\right)}{f^{\prime }\left({\tau }_j^{\text{g}}\right)}$

La méthode de dichotomie peut également être employée.

À chaque itération les sommes suivantes sont calculées :

${\mathrm{\Sigma }}^{\text{g}}=\sum _{i=1}^N\frac{K_i{x}_i}{1+(K_i-1){\tau }^{\text{g}}}$

${\mathrm{\Sigma }}^{\text{l}}=\sum _{i=1}^N\frac{x_i}{1+(K_i-1){\tau }^{\text{g}}}$

Ces sommes sont censées valoir ${\mathrm{\Sigma }}^{\text{g}}=\sum _{i=1}^N{x}_i^{\text{g}}=1$ et ${\mathrm{\Sigma }}^{\text{l}}=\sum _{i=1}^N{x}_i^{\text{l}}=1$, ce qui n'est vrai que lorsque le calcul a convergé. Pour pouvoir calculer les coefficients $K_i$, il est nécessaire de normer les fractions molaires au cours des itérations selon :

${\tilde{x}}_i^{\text{g}}=\frac{K_i{x}_i}{(1+(K_i-1){\tau }^{\text{g}})\cdot {\mathrm{\Sigma }}^{\text{g}}}$

${\tilde{x}}_i^{\text{l}}=\frac{x_i}{(1+(K_i-1){\tau }^{\text{g}})\cdot {\mathrm{\Sigma }}^{\text{l}}}$

Lorsque ${\mathrm{\Sigma }}^{\text{g}}={\mathrm{\Sigma }}^{\text{l}}=1$, on a bien ${\tilde{x}}_i^{\text{g}}=x_i^{\text{g}}$, ${\tilde{x}}_i^{\text{l}}=x_i^{\text{l}}$ et $f\left({\tau }^{\text{g}}\right)={\mathrm{\Sigma }}^{\text{g}}-{\mathrm{\Sigma }}^{\text{l}}=0$.

Le calcul est initialisé par le choix d'un taux de vaporisation $0<{\tau }^{\text{g}}<1$, de fractions molaires $x_i^{\text{g}}$ et $x_i^{\text{l}}$, et de coefficients de partage $K_i$. L'algorithme de résolution comporte deux boucles imbriquéesModèle:, : une sur les coefficients $K_i$, l'autre sur ${\tau }^{\text{g}}$.

Algorithme de calcul d'une vaporisation partielle

1. Calculer les sommes ${\mathrm{\Sigma }}^{\text{g}}$ et ${\mathrm{\Sigma }}^{\text{l}}$.
2. Calculer les fractions molaires normées ${\tilde{x}}_i^{\text{g}}$ et ${\tilde{x}}_i^{\text{l}}$.
3. Calculer les coefficients de partage $K_i$ par les coefficients de fugacité et d'activité.
4. Si les équilibres ne sont pas vérifiés ($|{\tilde{x}}_i^{\text{g}}-K_i{\tilde{x}}_i^{\text{l}}|>ϵ$ pour au moins un constituant), reprendre en 1 avec les nouvelles fractions et les nouveaux coefficients de partage. Sinon, passer en 5.
5. Résoudre l'équation de Rachford-Rice en ${\tau }^{\text{g}}$. On obtient ${\tau }_{k+1}^{\text{g}}$ ($k$ est l'itération en cours).
6. Si $|{\tau }_{k+1}^{\text{g}}-{\tau }_k^{\text{g}}|>ϵ$, alors reprendre en 1 avec le nouveau taux. Sinon, FIN.

Notes

• Le passage de l'étape 4 à l'étape 5 s'effectue même si ${\mathrm{\Sigma }}^{\text{g}}\ne 1$ ou ${\mathrm{\Sigma }}^{\text{l}}\ne 1$.
• À l'issue du calcul, si ${\tau }^{\text{g}}\ge 1$, le mélange est entièrement gazeux ; si ${\tau }^{\text{g}}\le 0$, le mélange est entièrement liquide. Des calculs préalables de pression de bulle et de rosée peuvent être effectués afin de vérifier ces possibilités et éviter de résoudre le problème de vaporisation partielle.

Équilibres idéaux

Solution idéale

Fichier:Fugacité - Henry - Lewis et Randall.jpg

Évolution de la fugacité en fonction de la fraction molaire à pression et température constantesModèle:,.

Le terme solution est pris au sens large de mélange. La notion de solution idéale s'applique à toutes les phases, quel que soit leur état. Dans une solution idéale, la fugacité idéale $f_i^{\text{id}}$ est supposée évoluer, à pression et température données, selon la relation :

Solution idéale : $\frac{f_i^{\text{id}}}{x_i}=\frac{f_i^{\circ }}{x_i^{\circ }}$

dans laquelle :

• $f_i^{\text{id}}$ est la fugacité idéale à la fraction molaire $x_i$ ;
• $f_i^{\circ }$ est la fugacité à la fraction molaire de référence $x_i^{\circ }$ (fixe).

Les deux fugacités sont prises aux mêmes pression et température, et pour la même phase. À pression et température données, la fugacité $f_i^{\circ }$ est une constante, et la fugacité $f_i^{\text{id}}$ varie linéairement en fonction de $x_i$. Les fugacités sont parfois exprimées en fonction de la fraction massique, de la concentration molaire ou de la concentration massique.

Dans les modèles idéaux, les coefficients de fugacité ${\varphi }_i$ et d'activité ${\gamma }_i$ sont égaux à 1. Pour un gaz, la solution idéale est le mélange de gaz parfaits. Pour les espèces chimiques en phase liquide, plusieurs états de référence sont souvent employés : le corps pur ($x_i^{\circ }=1$), le corps infiniment dilué ($x_i^{\circ }=0$), le corps à une concentration molaire de Modèle:Unité. Pour les espèces chimiques en phase solide, l'état de référence est le plus souvent le corps pur.

Équilibres liquide-vapeur idéaux

Loi de Raoult et loi de Henry

Pour un équilibre liquide-vapeur, l'égalité des fugacités s'écrit :

Équilibre liquide-vapeur : $f_i^{\text{g}}=f_i^{\text{l}}$

Pour la phase gazeuse, le modèle idéal est celui du gaz parfait. La fugacité de tout corps $i$ est égale à la pression partielle de ce corps :

$f_i^{\text{g}}=f_i^{\bullet }=x_i^{\text{g}}P$

Pour un corps $i$ en phase liquide, lorsque l'état de référence choisi est celui du corps pur, $x_i^{\text{l},\circ }=1$, alors $f_i^{\text{l},\circ }=f_i^{\text{l},*}$ est la fugacité du corps pur liquide. La fugacité suit la loi de Lewis et Randall : $f_i^{\text{l}}=x_i^{\text{l}}{f}_i^{\text{l},*}$. La fugacité du corps pur liquide se calcule de façon générale selon (voir la section Fugacités et coefficient de partage) :

$f_i^{\text{l},*}={\varphi }_i^{\text{g,*,sat}}{P}_i^{\text{sat}}{𝒫}_i^{\text{l}}$

Dans les équilibres idéaux ${\varphi }_i^{\text{g,*,sat}}=1$, et aux basses pressions ${𝒫}_i^{\text{l}}\approx 1$. On obtient la loi de Raoult :

Loi de Raoult : $x_i^{\text{g}}P=x_i^{\text{l}}{P}_i^{\text{sat}}$

Pour un corps $i$ en phase liquide, lorsque l'état de référence choisi est celui du corps à dilution infinie, $x_i^{\text{l},\circ }=0$, alors $f_i^{\text{l},\circ }=0$ et $\frac{f_i^{\text{l},\circ }}{x_i^{\text{l},\circ }}=k_{\text{H},i}$ la constante de Henry. On obtient la loi de Henry :

Loi de Henry : $x_i^{\text{g}}P=x_i^{\text{l}}{k}_{\text{H},i}$

Les lois de Raoult et de Henry considèrent toutes deux la phase gaz comme un mélange de gaz parfaits, ce qui restreint leur domaine d'application à des pressions de moins de Modèle:Unité. La loi de Raoult s'applique pour des fractions molaires en phase liquide plutôt grandes, $x_i^{\text{l}}\approx 1$, et la loi de Henry pour des fractions molaires plutôt petites, $x_i^{\text{l}}\approx 0$.

Exemple de calcul d'un équilibre idéal

On considère un équilibre liquide-vapeur dont tous les constituants suivent la loi de Raoult. Les données du problème sont la pression $P$, la température $T$ et les fractions molaires globales $x_i$. À température donnée, le coefficient de partage de chacun des $N$ constituants est une fonction de la pression et ne dépend pas de la composition :

$K_i\left(P\right)=\frac{P_i^{\text{sat}}{𝒫}_i^{\text{l}}\left(P\right)}{P}$

Les coefficients de partage sont donc invariants au cours d'un calcul par itérations.

Avant de calculer un équilibre liquide-vapeur, il est nécessaire de vérifier que le mélange n'est pas monophasique, c'est-à-dire uniquement liquide ou gazeux. 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\frac{x_i}{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$.

On effectue par conséquent le calcul selon l'algorithme suivant.

Algorithme de calcul

1. À pression et température de l'énoncé, calculer le coefficient de partage $K_i$ pour chacun des $N$ constituants.
2. 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. Sinon, passer en 3.
3. Si $\sum _{i=1}^N\frac{x_i}{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, passer en 4.
4. Sinon, le mélange est biphasique et présente un équilibre liquide-vapeur. Résoudre l'équation de Rachford-Rice : $\sum _{i=1}^N\frac{(K_i-1)x_i}{1+(K_i-1){\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}}=\frac{1-K_{\text{B}}}{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}}=\frac{x_{\text{A}}-x_{\text{A}}^{\text{l}}}{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 à Modèle:Unité et Modèle:Unité - 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 Modèle:Nobr et Modèle:Nobr. Le mélange est donc biphasique. Les résultats sont reportés dans le tableau suivant.

Mélange d'hydrocarbures à Modèle:Unité et Modèle:Unité - résultats.

Taux de vaporisation ${\tau }^{\text{g}}$	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

Équilibres liquide-solide idéaux

L'équation de Schröder-van Laar permet de déterminer, à pression constante, l'équilibre liquide-solide de deux phases idéales, dans lequel seule l'espèce $\text{A}$ est présente dans la phase solide. Elle s'écrit en fonction de la composition de la phase liquideModèle:, :

Équation de Schröder-van Laar $\ln x_{\text{A}}^{\text{l}}=\frac{{\mathrm{\Delta }}_{\text{fus}}{H}_{\text{A}}}{R{T}_{\text{fus,A}}}[1-\frac{T_{\text{fus,A}}}{T}]-\frac{{\mathrm{\Delta }}_{\text{fus}}{C}_{P,\text{A}}}{R}[1-\frac{T_{\text{fus,A}}}{T}+\ln \left(\frac{T_{\text{fus,A}}}{T}\right)]$

avec :

• $P$ la pression de l'équilibre liquide-solide ;
• $R$ la constante universelle des gaz parfaits ;
• $T$ la température de l'équilibre liquide-solide sous la pression $P$ ;
• $T_{\text{fus,A}}$ la température du fusion du corps $\text{A}$ pur sous la pression $P$ ;
• $x_{\text{A}}^{\text{l}}$ la fraction molaire du corps $\text{A}$ dans la phase liquide ;
• ${\mathrm{\Delta }}_{\text{fus}}{C}_{P,\text{A}}$ l'écart des capacités thermiques isobares molaires du corps $\text{A}$ pur liquide et solide à $T_{\text{fus,A}}$ :
  • ${\mathrm{\Delta }}_{\text{fus}}{C}_{P,\text{A}}={\overline{C}}_{P,\text{A}}^{\text{l,*}}\left(T_{\text{fus,A}}\right)-{\overline{C}}_{P,\text{A}}^{\text{s,*}}\left(T_{\text{fus,A}}\right)$ ;
  • ${\overline{C}}_{P,\text{A}}^{\text{l,*}}$ la capacité thermique isobare molaire du corps $\text{A}$ liquide pur ;
  • ${\overline{C}}_{P,\text{A}}^{\text{s,*}}$ la capacité thermique isobare molaire du corps $\text{A}$ solide pur ;
• ${\mathrm{\Delta }}_{\text{fus}}{H}_{\text{A}}$ l'enthalpie de fusion du corps $\text{A}$ pur à $T_{\text{fus,A}}$.

Propriétés colligatives

• Diagramme de phases d'un solvant pur (courbes pleines) et du même solvant en présence d'un soluté (pointillés). Les propriétés colligatives se traduisent par un déplacement des courbes d'équilibre solide-liquide et gaz-liquide.
  Diagramme de phases d'un solvant pur (courbes pleines) et du même solvant en présence d'un soluté (pointillés). Les propriétés colligatives se traduisent par un déplacement des courbes d'équilibre solide-liquide et gaz-liquide.
• À l'équilibre, la pression de part et d'autre de la membrane semi-perméable est plus élevée dans le compartiment de la solution que dans celui du solvant pur. La différence des pressions est appelée « pression osmotique ».
  À l'équilibre, la pression de part et d'autre de la membrane semi-perméable est plus élevée dans le compartiment de la solution que dans celui du solvant pur. La différence des pressions est appelée « pression osmotique ».

Dans une solution liquide, le « solvant » est l'espèce chimique majoritaire et le « soluté » l'espèce minoritaire. La comparaison des conditions opératoires de transition de phase du solvant pur avec celles de la solution montre que la présence du soluté modifie ces conditions : à température identique, la pression de transition est modifiée, de même qu'à pression identique la température de transition est modifiée. Dans le cas particulier où la phase en équilibre avec la solution ne contient que le solvant et où le soluté n'est présent qu'en très faible quantité dans la solution, les écarts constatés sont appelés « propriétés colligative ». Dans ces conditions, ces écarts peuvent être calculés, selon l'hypothèse de la solution idéale, à l'aide de lois ne dépendant que des propriétés du solvant et linéaires selon la concentration du soluté.

À température $T$ donnée, la pression d'équilibre liquide-vapeur de la solution est plus basse que celle du solvant pur. L'abaissement de la pression de vapeur saturante se calcule selon la loi de la tonométrie :

Loi de la tonométrie : $\frac{\mathrm{\Delta }{P}^{\text{sat}}}{P^{\text{sat}}}=x_{\sigma }^{\text{l}}$

avec :

• $\mathrm{\Delta }{P}^{\text{sat}}$ l'abaissement de la pression de vapeur saturante du solvant en présence du soluté (en Modèle:Unité) ; la pression d'équilibre liquide-vapeur de la solution vaut $P^{\text{sat}}-\mathrm{\Delta }{P}^{\text{sat}}$ ;
• $P^{\text{sat}}$ la pression de vapeur saturante du solvant pur à la température $T$ (en Modèle:Unité) ;
• $x_{\sigma }^{\text{l}}$ la fraction molaire du soluté en solution.

À pression $P$ donnée, la température d'équilibre liquide-vapeur de la solution est plus élevée que celle du solvant pur. L'élévation de la température d'ébullition se calcule selon la loi de l'ébulliométrie :

Loi de l'ébulliométrie : $\mathrm{\Delta }{T}_{\text{vap}}=K_{\text{éb}}\cdot x_{\sigma }^{\text{l}}$ Constante ébullioscopique : $K_{\text{éb}}=\frac{R{T_{\text{vap}}}^2}{{\mathrm{\Delta }}_{\text{vap}}H}$

avec :

• $\mathrm{\Delta }{T}_{\text{vap}}$ l'élévation de la température d'ébullition du solvant (en Modèle:Unité) ; la température d'ébullition de la solution vaut $T_{\text{vap}}+\mathrm{\Delta }{T}_{\text{vap}}$ ;
• $T_{\text{vap}}$ la température d'ébullition du solvant pur à la pression $P$ (en Modèle:Unité) ;
• ${\mathrm{\Delta }}_{\text{vap}}H$ l'enthalpie de vaporisation du solvant pur à $T_{\text{vap}}$ (en Modèle:Unité) ;
• $R$ la constante universelle des gaz parfaits (en Modèle:Unité).

À pression $P$ donnée, la température d'équilibre liquide-solide de la solution est plus basse que celle du solvant pur. L'abaissement de la température de fusion (ou température de solidification ou température de congélation) se calcule selon la loi de la cryométrie :

Loi de la cryométrie : $\mathrm{\Delta }{T}_{\text{fus}}=K_{\text{cry}}\cdot x_{\sigma }^{\text{l}}$ Constante cryoscopique : $K_{\text{cry}}=\frac{R{T_{\text{fus}}}^2}{{\mathrm{\Delta }}_{\text{fus}}H}$

avec :

• $\mathrm{\Delta }{T}_{\text{fus}}$ l'abaissement de la température de fusion du solvant (en Modèle:Unité) ; la température de fusion de la solution vaut $T_{\text{fus}}-\mathrm{\Delta }{T}_{\text{fus}}$ ;
• $T_{\text{fus}}$ la température de fusion du solvant pur à la pression $P$ (en Modèle:Unité) ;
• ${\mathrm{\Delta }}_{\text{fus}}H$ l'enthalpie de fusion du solvant pur à $T_{\text{fus}}$ (en Modèle:Unité).

À température $T$ donnée, si l'on place la solution et le solvant pur de part et d'autre d'une membrane semi-perméable (ne laissant passer que le solvant), la pression d'équilibre est plus élevée dans le compartiment de la solution que dans celui du solvant. La différence de pression entre les deux compartiments, appelée « pression osmotique », se calcule selon la loi de l'osmométrie :

Loi de van 't Hoff, ou loi de l'osmométrie : $\mathrm{\Pi }=c_{\sigma }^{\text{l}}RT$

avec :

• $\mathrm{\Pi }$ la pression osmotique (en Modèle:Unité) ; si au niveau de la membrane la pression du compartiment du solvant pur vaut $P$, celle du côté de la solution vaut $P+\mathrm{\Pi }$ ;
• $c_{\sigma }^{\text{l}}$ la concentration molaire du soluté en solution (en Modèle:Unité) ;
• $T$ la température (en Modèle:Unité).

Notations

Alphabet latin

Alphabet grec

Notes et références

Bibliographie

• Modèle:Ouvrage.
• Modèle:Ouvrage.
• Modèle:Ouvrage.
• Modèle:Ouvrage.
• Modèle:Ouvrage.
• Modèle:Ouvrage.
• Modèle:Ouvrage.
• Modèle:Ouvrage.
• Modèle:Ouvrage.