On first-order thermodynamic equilibrium conditions for fluid-fluid and solid-phase interfaces

Résumé Technique : Sur les conditions d'équilibre thermodynamique du premier ordre pour les interfaces fluide-fluide et solide-phase

Énoncé du Problème
La détermination des propriétés interfaciales dans les systèmes solide-fluide est cruciale pour la compréhension de phénomènes en science des matériaux, en métallurgie et en géophysique, tels que la nucléation et la solidification. Bien que la théorie thermodynamique des systèmes hétérogènes repose sur des fondements solides grâce aux travaux de Gibbs (pour les interfaces fluides) et aux développements ultérieurs de Larché et Cahn (pour les systèmes solide-fluide), les cadres existants traitent souvent la thermodynamique de la masse, de l'interface et de la configuration comme des cadres théoriques distincts. Plus précisément, la formulation de Larché-Cahn, bien qu'elle étende rigoureusement le travail de Gibbs aux solides via la contrainte de « conservation du réseau », n'inclut pas explicitement les contributions d'excès interfaciaux dans le fonctionnel contraint. Par conséquent, elle ne peut pas décrire pleinement les systèmes où les effets de surface sont significatifs (par exemple, à l'échelle nanométrique) ou lorsque les interfaces sont non planes et anisotropes. De plus, les relations classiques telles que l'équation de Young-Laplace, l'équation de Shuttleworth et le vecteur de Cahn-Hoffman sont souvent dérivées à partir d'hypothèses ad hoc plutôt que d'un principe variationnel unifié.

Méthodologie
Les auteurs développent un cadre thermodynamique unifié basé sur une formulation variationnelle contrainte unique. En partant de l'approche de Larché-Cahn pour les solides en contact avec des fluides, les auteurs étendent le problème de minimisation contrainte pour inclure explicitement les quantités d'excès interfaciaux (énergies de surface et de ligne) en tant que sous-systèmes thermodynamiques distincts.

La méthodologie repose sur la distinction de trois classes de variations appliquées à un fonctionnel thermodynamique (énergie interne sous contraintes de l'entropie totale, du nombre de moles total et du volume total) :

1. Variations Extérieures : Changements locaux des champs thermodynamiques (densité d'entropie, densité molaire) à des points spatiaux fixes.
2. Variations Intérieures : Déformations géométriques de la configuration spatiale actuelle (déplacement de points matériels), qui modifient le volume, la surface et l'état de déformation sans changer l'identité de la phase.
3. Variations Configurationnelles : Variations de l'identité de phase des points matériels (mouvement de la frontière de phase), qui changent la configuration matérielle. Pour les solides cristallins, cela inclut une condition de « conservation forte des étiquettes de réseau » où la structure du réseau est préservée lors des déplacements virtuels.

L'analyse est menée d'abord pour les systèmes fluide-fluide afin d'établir la structure mathématique, puis étendue aux interfaces solide-fluide. La densité d'énergie interne est supposée dépendre des variables d'état locales (entropie, composition) et, pour les solides et les interfaces, des caractéristiques géométriques telles que le gradient de déformation ($F$) et la normale de l'interface ($n$).

Contributions Clés et Résultats

• Dérivation Unifiée des Relations Classiques : Le document démontre que les relations d'équilibre classiques émergent naturellement comme des conditions de stationnarité d'un seul fonctionnel thermodynamique sous différentes classes de variations :

  • Les variations extérieures permettent de retrouver les conditions d'équilibre de la masse (égalité des températures et des potentiels chimiques).
  • Les variations intérieures dans la masse produisent l'équilibre mécanique ($\nabla \cdot \mathbf{S} = 0$).
  • Les variations intérieures à l'interface produisent l'équation généralisée de Young-Laplace, la relation de Shuttleworth (liant la contrainte de surface et l'énergie libre de surface) et le vecteur de Cahn-Hoffman (pour les interfaces anisotropes).
  • Les variations configurationnelles produisent le tenseur énergie-moment de Eshelby et la force motrice du mouvement de la frontière de phase.

• Contrainte de Surface Généralisée et Équilibre :

  • Pour les interfaces fluide-fluide, le cadre retrouve l'équation standard de Young-Laplace ($2\gamma H + P_\alpha - P_\beta = 0$).
  • Pour les interfaces solide-fluide, les auteurs dérivent un tenseur de contrainte de surface généralisé $\boldsymbol{\sigma}_\Sigma$ qui tient compte de la dépendance à la fois vis-à-vis de l'orientation de l'interface ($n$) et du gradient de déformation ($F$). Ce tenseur est montré comme étant symétrique sous l'indifférence de cadre.
  • La condition d'équilibre à l'interface est exprimée par un équilibre des tractions : $\nabla_\Sigma \cdot \boldsymbol{\sigma}_\Sigma - (\mathbf{S}_s - \mathbf{S}_l)\mathbf{n} = 0$, où $\mathbf{S}_s$ et $\mathbf{S}_l$ sont les contraintes de Cauchy dans le solide et le fluide, respectivement.
  • Une nouvelle condition est dérivée : $\frac{\partial \omega_\Sigma}{\partial \mathbf{F}} \mathbf{F}^T \mathbf{n} = 0$, qui stipule que la contrainte énergétique interfaciale dépendante de la déformation ne produit aucune traction dans la direction normale à l'interface.

• Équilibre Configurationnel et Conservation du Réseau :

  • Le document distingue la conservation de « l'étiquette faible » (Larché-Cahn), où le réseau est une continuation triviale, de la conservation de « l'étiquette forte », où le transport des étiquettes de réseau est explicitement suivi.
  • Sous une conservation forte, la force motrice du mouvement de l'interface est régie par la composante normale du tenseur d'Eshelby ($\mathbf{N} \cdot \mathbf{E} \mathbf{N}$), plutôt que par la seule densité de grand potentiel scalaire.
  • La condition classique de Larché-Cahn ($\omega_s = \omega_l$) est récupérée comme un cas particulier où l'état de déformation est maintenu fixe pendant l'avancement de l'interface.

• Distinction Entre Pressions Mécanique et Thermodynamique :

  • Les auteurs clarifient que dans les solides déformables, la densité du grand potentiel de la masse ($\omega_s$) et la pression thermodynamique ($P_s$) sont deux quantités distinctes en raison de la résistance élastique à la déformation. Cela conduit à une équation de Young-Laplace généralisée où le saut de pression est équilibré par les termes de contrainte de surface et la courbure, plutôt que par la seule énergie libre de surface.

Signification et Revendications
Le document revendique de fournir une « base théorique commune » pour la thermodynamique de la masse, de l'interface et de la configuration. En formulant l'équilibre des systèmes hétérogènes comme un problème de minimisation contrainte unique, les auteurs montrent que les théories historiquement séparées (Gibbs, Larché-Cahn, Shuttleworth, Cahn-Hoffman, Eshelby) sont des conditions de stationnarité complémentaires d'un même fonctionnel.

Le travail clarifie les hypothèses selon lesquelles des résultats spécifiques sont valides, particulièrement concernant la distinction entre l'équilibre mécanique (équilibre des forces) et l'équilibre configurationnel (forces motrices de transformation de phase). Les auteurs notent que leur formulation est limitée aux conditions nécessaires du premier ordre (équilibre) et ne vérifie pas la convexité locale (stabilité). Ils limitent également la portée aux phases cristallines uniques en contact avec des fluides, laissant les interfaces solide-solide cohérentes aux travaux futurs.

Le cadre est présenté comme un outil pour interpréter les simulations atomistiques et développer des méthodes de calcul rigoureuses pour les propriétés thermodynamiques interfaciales, particulièrement dans les contextes où les effets de surface sont significatifs et où les approximations classiques des fluides échouent. Le document stipule explicitement que, bien qu'il récupère la voie mécanique de Kirkwood-Buff comme un cas limite, il explique pourquoi cette voie échoue pour les interfaces cristallines où la contrainte de surface et l'énergie libre de surface divergent en raison de la déformation et de la dépendance à l'orientation.