Trempez un anneau dans de l’eau savonneuse et soufflez doucement : un mince film liquide s’étire, se pince et se détend en quelques millisecondes pour former une sphère irisée parfaite flottant dans l’air. Observez ensuite une mousse de savon à la loupe : ce qui semble être un enchevêtrement aléatoire obéit en réalité à une géométrie implacable. Dans toute la mousse, les films liquides ne se rejoignent jamais par quatre ou cinq, mais exclusivement par trois le long de canaux liquides appelés bords de Plateau à un angle exact de $120^\circ$ ($2\pi/3$), et ces bords se croisent exclusivement par quatre aux sommets selon l’angle tétraédrique exact de $\arccos(-1/3) \approx 109,4712^\circ$.

Pourquoi une bulle de savon libre est-elle toujours sphérique ? Pourquoi, lorsqu’on relie une petite bulle et une grosse bulle par un tube, est-ce la petite bulle qui se vide dans la grosse — et comment nos 300 millions d’alvéoles pulmonaires échappent-elles à cet effondrement à chaque expiration ? Pourquoi l’eau pure ou l’alcool pur ne peuvent-ils jamais former de bulles durables, alors qu’une trace de tensioactif confère au film une élasticité auto-cicatrisante ? Et comment la conjecture centenaire de Lord Kelvin sur le pavage de l’espace en cellules de même volume et d’aire minimale a-t-elle été battue en 1994 par la structure de Weaire-Phelan ?

Cet article mobilise la thermodynamique, la géométrie différentielle, le calcul des variations, la théorie géométrique de la mesure et l’hydrodynamique interfaciale pour dévoiler toute la physique des bulles et des mousses.

1. Thermodynamique de la tension superficielle et asymétrie moléculaire

Au cœur (bulk) de l’eau liquide, chaque molécule $\text{H}_2\text{O}$ est entourée de façon isotrope par ses voisines et liée par des liaisons hydrogène et des forces de van der Waals. La résultante moyenne des forces y est nulle et l’énergie potentielle cohésive est minimale (coordinence $z_b \approx 12$, énergie $-u_0$).

À l’interface liquide-vapeur, cette symétrie est brutalement rompue : la densité moléculaire de la phase gazeuse est environ mille fois plus faible que celle du liquide. Une molécule située à la surface perd près de la moitié de ses voisines attractives ($z_s < z_b$) et possède donc un excès d’énergie potentielle. Pour amener des molécules du volume vers la surface et accroître l’aire interfaciale de $dA$, il faut fournir un travail réversible positif contre les forces de cohésion.

La différentielle totale de l’enthalpie libre de Gibbs $G(T, P, A, \{N_i\})$ s’écrit :

d’où découle la définition thermodynamique rigoureuse de la tension superficielle $\gamma$ :

Elle possède deux interprétations physiques équivalentes :

1. Énergie libre en excès par unité d’aire ($\text{J/m}^2$) : À volume $V$ fixé, un système minimise son énergie libre $\gamma A$ en minimisant son aire $A$. En vertu de l’inégalité isopérimétrique tridimensionnelle $A^3 \ge 36\pi V^2$, la sphère est l’unique forme minimisant la surface pour un volume donné.
2. Force de traction tangentielle par unité de longueur ($\text{N/m}$) : Selon la formule de Kirkwood-Buff (1949) en mécanique statistique, $\gamma$ est l’intégrale de l’anisotropie du tenseur des pressions entre la composante normale $P_N$ et la composante tangentielle $P_T(z)$ à travers l’interface : $$\gamma = \int_{-\infty}^{+\infty} \bigl[ P_N - P_T(z) \bigr]\,dz$$

2. Géométrie différentielle des surfaces courbes et équation de Young-Laplace

Lorsqu’une interface est courbe, la tension tangentielle $\gamma$ engendre une force normale nette dirigée vers le côté concave, imposant une surpression $\Delta P = P_{\text{int}} - P_{\text{ext}} > 0$. Soient $R_1, R_2$ les deux rayons de courbure principaux ($\kappa_1 = 1/R_1, \kappa_2 = 1/R_2$) d’une surface régulière $\Sigma \subset \mathbb{R}^3$. La courbure moyenne $H$ et la courbure de Gauss $K$ sont définies par :

Un déplacement normal virtuel $\delta\xi$ d’un élément de surface $dA = R_1 R_2\,d\theta_1\,d\theta_2$ induit une variation d’aire $\delta(dA) = dA\left(\frac{1}{R_1} + \frac{1}{R_2}\right)\delta\xi$ pour un volume balayé $\delta(dV) = dA\,\delta\xi$. La condition d’équilibre $\delta F = \gamma\,\delta(dA) - \Delta P\,\delta(dV) = 0$ conduit à l’équation de Young-Laplace (1805–1806) :

Alors qu’une goutte d’eau ne possède qu’une seule interface ($\Delta P_{\text{goutte}} = 2\gamma/R$), une bulle de savon dans l’air est une mince coquille liquide d’épaisseur $h \ll R$ possédant deux interfaces gaz-liquide (interne de rayon $R$ et externe de rayon $R+h$), d’où le facteur $4$ :

3. Surfaces minimales ($H = 0$), catastrophe de la caténoïde et problème de Plateau

Un film de savon ouvert tendu sur un cadre métallique non planaire subit la même pression atmosphérique sur ses deux faces ($\Delta P = 0$). L’équation de Young-Laplace impose alors que sa courbure moyenne soit identiquement nulle en tout point :

Ces surfaces à courbure moyenne nulle sont appelées surfaces minimales. Sous forme de graphe $z = u(x,y)$, la minimisation de la fonctionnelle d’aire $A[u] = \iint \sqrt{1+u_x^2+u_y^2}\,dx\,dy$ donne l’équation aux dérivées partielles d’Euler-Lagrange établie par Lagrange en 1760 :

Tout point non plan d’un film de savon ouvert est localement une selle hyperbolique ($K < 0$). La surface minimale de révolution tendue entre deux anneaux coaxiaux de rayon $R$ distants de $2h$ est la caténoïde découverte par Euler en 1744 :

Lorsque l’on écarte les deux anneaux au-delà du rapport critique ($\coth u_c = u_c \approx 1,19968$) :

la solution réelle disparaît par bifurcation nœud-col : le col du film s’effondre brusquement et se rompt pour former deux disques plats séparés (solution discontinue de Goldschmidt). Le problème de Plateau (existence d’une surface minimale bordée par toute courbe de Jordan fermée) fut résolu en 1930 par Jesse Douglas (première médaille Fields en 1936) et Tibor Radó.

4. Paradoxe des deux bulles, mûrissement d’Ostwald et surfactant pulmonaire (DPPC)

Si l’on relie par un tube muni d’un robinet une petite bulle de savon ($R_{\text{petite}}$) et une grosse bulle ($R_{\text{grosse}}$), l’ouverture du robinet provoque le dégonflement rapide de la petite bulle qui gonfle la grosse bulle ! En effet, d’après $\Delta P = 4\gamma/R$, la surpression est inversement proportionnelle au rayon :

  flowchart TD
    A["Deux bulles reliées par un tube à vanne :<br/>Petite bulle (R_petite) et Grosse bulle (R_grosse)"] --> B["Équation de Young-Laplace : ΔP = 4γ / R<br/>La surpression est inversement proportionnelle au rayon R"]
    B --> C["Pression petite bulle : ΔP_petite = 4γ / R_petite<br/>Pression grosse bulle : ΔP_grosse = 4γ / R_grosse"]
    C --> D["Déséquilibre : ΔP_petite > ΔP_grosse<br/>(Une plus forte courbure engendre une pression supérieure !)"]
    D --> E["Ouverture de la vanne : le gaz s'écoule de la petite bulle<br/>vers la grosse bulle"]
    E --> F["Rétroaction positive :<br/>La petite bulle rétrécit, sa pression grimpe et elle s'effondre"]
    D --> G["Crise physiologique dans les poumons :<br/>300 millions d'alvéoles de rayons variés (50–200 μm)"]
    G --> H["Les pneumocytes de type II sécrètent le surfactant (DPPC)"]
    H --> I["Quand l'alvéole rétrécit, la monocouche de DPPC se comprime,<br/>faisant chuter γ jusqu'à 0,5–2 mN/m"]
    I --> J["Égalisation des pressions (2γ_petite / R_petite ≈ 2γ_grosse / R_grosse)<br/>empêchant le collapsus alvéolaire (atélectasie) !"]

Au sein d’une mousse, ce transfert diffusif du gaz des petites bulles vers les grosses constitue le mûrissement d’Ostwald (coarsening). En 1952, John von Neumann a démontré grâce au théorème de Gauss-Bonnet que dans une mousse sèche 2D, l’aire $A_n$ d’une bulle polygonale à $n$ côtés évolue selon la loi de von Neumann :

Les bulles ayant $n < 6$ côtés rétrécissent et disparaissent, les hexagones ($n = 6$) sont stationnaires et les bulles à $n > 6$ côtés grossissent (la généralisation 3D a été établie par MacPherson et Srolovitz en 2007).

Dans les poumons humains, quelque 300 millions d’alvéoles ($R \approx 50\text{–}200\,\mu\text{m}$) communiquent entre elles. Si la tension superficielle $\gamma$ du liquide alvéolaire était constante, les petites alvéoles se videraient dans les grandes (atélectasie). Les pneumocytes de type II évitent cette catastrophe en sécrétant le surfactant pulmonaire, riche en dipalmitoylphosphatidylcholine (DPPC) : lors de l’expiration, la compression de la monocouche insoluble de DPPC abaisse $\gamma(R)$ vers $0,5\text{–}2\,\text{mN/m}$, stabilisant l’ensemble des alvéoles.

5. Lois de Joseph Plateau (1873), théorème de Jean Taylor (1976) et théorème de la bulle double

Formulées expérimentalement par le physicien belge aveugle Joseph Plateau en 1873 et démontrées rigoureusement par la mathématicienne Jean E. Taylor en 1976 via la théorie géométrique de la mesure, les lois d’équilibre des mousses stipulent :

1. 1re loi de Plateau (3 films à $120^\circ$) : Les films de savon se rejoignent exclusivement par trois le long des bords de Plateau à des angles égaux de $120^\circ$ ($2\pi/3$) ($\sum_{i=1}^3 \hat{\mathbf{t}}_i = \mathbf{0} \implies \cos\theta = -1/2$). Une jonction quadruple à $90^\circ$ ($L_{\text{croix}} = 4\sqrt{2} \approx 5,657$) est instable et se scinde immédiatement en deux jonctions triples à $120^\circ$ selon un arbre de Steiner ($L_{\text{Steiner}} = 2(1+\sqrt{3}) \approx 5,464$, transition T1).
2. 2e loi de Plateau (4 bords à $109,47^\circ$) : Les bords de Plateau se croisent exclusivement par quatre aux sommets selon l’angle tétraédrique de Maraldi : $$\theta_{\text{Maraldi}} = \arccos\!\left(-\frac{1}{3}\right) \approx 109,47122^\circ \quad \left(109^\circ 28' 16''\right)$$
3. Théorème de la bulle double (Hutchings, Morgan, Ritoré, Ros 2000) : Lorsque deux bulles de rayons $R_1 > R_2$ fusionnent, la paroi séparatrice (septum) se courbe vers l’intérieur de la grande bulle avec un rayon de courbure $R_{\text{septum}}$ donné par : $$\frac{1}{R_{\text{septum}}} = \frac{1}{R_2} - \frac{1}{R_1} \quad \Longleftrightarrow \quad R_{\text{septum}} = \frac{R_1 R_2}{R_1 - R_2}$$

6. Conjecture de Kelvin (1887) contre structure de Weaire-Phelan (1994)

En 1887, Lord Kelvin chercha la partition de l’espace $\mathbb{R}^3$ en cellules de volume égal ($V=1$) minimisant l’aire interfaciale totale. En raison de la frustration géométrique (les dodécaèdres réguliers ne pavent pas l’espace), Kelvin proposa l’octaèdre tronqué (tétrakaidécaèdre à 6 carrés et 8 hexagones, cellule de Voronoï du réseau BCC) légèrement incurvé, d’aire normalisée $A_{\text{Kelvin}} \approx 5,30628$.

Après 107 ans de règne, les physiciens Denis Weaire et Robert Phelan (Trinity College Dublin) ont réfuté la conjecture de Kelvin en 1994 grâce à la structure de Weaire-Phelan (inspirée des clathrates hydrates de type I), dont la maille cubique réunit 8 polyèdres de même volume et de deux formes : 2 dodécaèdres pentagonaux irréguliers ($5^{12}$) et 6 tétrakaidécaèdres ($5^{12}6^2$ à 2 hexagones et 12 pentagones). Avec $88,9\%$ de faces pentagonales (angle interne de $108^\circ$ très proche de $109,47^\circ$), l’aire descend à $A_{\text{WP}} \approx 5,28825$ ($-0,34\%$ sous Kelvin), structure rendue célèbre par le Centre national de natation de Pékin (« Water Cube »).

7. Élasticité de Gibbs-Marangoni, interférences et pression de disjonction DLVO

Une faible tension superficielle ne suffit pas à stabiliser une bulle (l’éthanol pur éclate instantanément). La survie d’un film de savon repose sur l’élasticité de Gibbs-Marangoni :

  flowchart TD
    A["Perturbation externe (courant d'air, drainage, onde thermique)<br/>étirant et amincissant localement le film de savon"] --> B["L'aire locale ΔA augmente ;<br/>la concentration surfacique en tensioactif Γ chute localement"]
    B --> C["La baisse de Γ provoque une hausse locale de la tension<br/>superficielle γ_mince au-dessus de γ_épais environnant !"]
    C --> D["Apparition d'un gradient tangentiel ∇_s γ<br/>(contrainte de cisaillement de Marangoni τ_M = ∇_s γ vers la zone amincie)"]
    D --> E["Convection de Marangoni :<br/>La zone à fort γ tire la monocouche et entraîne l'eau sous-jacente par viscosité"]
    E --> F["L'épaisseur h s'auto-répare en quelques millisecondes,<br/>empêchant la rupture capillaire !"]
    F --> G["Le drainage gravitationnel amincit le film aux longueurs d'onde visibles (100–700 nm)"]
    G --> H["Les interférences en couche mince produisent des franges irisées"]
    H --> I["Pour h < 30 nm, la pression de disjonction DLVO Π(h) équilibre van der Waals ;<br/>le déphasage π crée le 'film noir de Newton' (h ≈ 4 nm)"]

Lorsqu’un film s’amincit localement, la densité surfacique de tensioactif $\Gamma$ diminue et la tension superficielle locale $\gamma_{\text{mince}}$ augmente. La contrainte de Marangoni $\boldsymbol{\tau}_M = \nabla_s \gamma$ et le module d’élasticité de Gibbs

ramènent immédiatement du liquide vers la zone amincie (le même effet explique les « larmes du vin »). Lorsque l’épaisseur $h$ atteint les longueurs d’onde visibles, les réflexions sur les deux faces produisent des interférences en couche mince ($2n_{\text{film}}h\cos\theta_t + \lambda/2 = m\lambda$). En dessous de $h \ll \lambda/4 \approx 30\,\text{nm}$, le déphasage de $\pi$ à la réflexion avant entraîne une interférence destructive pour toutes les longueurs d’onde : le film devient noir (film noir de Newton, $h \approx 4\,\text{nm}$), maintenu par la pression de disjonction DLVO $\Pi(h) = \Pi_{\text{vdW}}(h) + \Pi_{\text{el}}(h) + \Pi_{\text{sr}}(h)$.