Poser le pied sur une flaque gelée par un matin d’hiver, admirer un patineur de vitesse négociant un virage à $60\,\text{km/h}$ ou observer une pierre de curling en granit de $20\,\text{kg}$ glisser sur quarante mètres sans presque ralentir, c’est être témoin de l’une des anomalies les plus fascinantes de la physique de la matière condensée et de la tribologie (la science du frottement, de l’usure et de la lubrification). Si notre intuition quotidienne considère comme évident que « la glace glisse », la question « Pourquoi la glace cristalline est-elle si exceptionnellement glissante ? » demeure une énigme scientifique profonde qui a suscité des débats passionnés parmi les physiciens depuis plus de 170 ans.

Lorsque deux solides cristallins secs ordinaires — tels que l’acier sur le verre, le granit sur le quartz ou le cuivre sur l’alumine — glissent l’un contre l’autre sans lubrifiant, leur coefficient de frottement cinétique est élevé, typiquement compris entre $\mu_k \approx 0.4$ et $0.8$. À l’échelle microscopique, les aspérités de surface forment des jonctions adhésives qui doivent être continuellement cisaillées avec une forte dissipation d’énergie. Pourtant, lorsqu’une lame de patin en acier ou une semelle de ski en polyéthylène glisse sur de la glace hexagonale ordinaire (Glace $I_h$) entre $-5^\circ\text{C}$ et $-10^\circ\text{C}$ — bien en dessous du point de fusion de $0^\circ\text{C}$ —, le coefficient de frottement chute brutalement à $\mu_k \approx 0.005\text{–}0.02$, soit un à deux ordres de grandeur plus bas que n’importe quel autre solide cristallin, surpassant même le polytétrafluoroéthylène (PTFE / Téflon, $\mu_k \approx 0.04\text{–}0.10$).

La plupart des manuels scolaires et des encyclopédies populaires proposent depuis plus d’un siècle l’explication suivante : « Comme l’eau se dilate en gelant, l’application d’une pression sur la glace abaisse son point de fusion. Lorsqu’un patineur se tient sur une lame très fine, son poids se concentre sur une surface minuscule, générant une pression énorme qui fait fondre la glace sous la lame pour former un film d’eau lubrifiant. »

C’est la célèbre théorie de la fusion sous pression (Pressure-Melting Theory), formulée thermodynamiquement par James Thomson en 1849 et appliquée au patinage par John Joly en 1886. Or, une vérification quantitative rigoureuse révèle un constat sans appel : bien que l’équation de Clapeyron soit parfaitement exacte, l’hypothèse de la fusion sous pression pour expliquer le patinage, le ski ou les glissades sur un trottoir verglacé est quantitativement fausse de plusieurs ordres de grandeur !

Dans cet article, nous allons d’abord dériver l’équation de Clausius-Clapeyron à partir de l’équilibre du potentiel chimique et démontrer mathématiquement pourquoi la pression d’un patin n’abaisse le point de fusion de la glace que d’environ un dixième de degré Celsius. Nous analyserons ensuite l’expérience classique du regel de Faraday et du fil lesté traversant un bloc de glace par la conduction thermique de Fourier, expliquant pourquoi un fil de cuivre traverse la glace alors qu’un fil de nylon sous la même pression échoue. Enfin, nous exposerons les trois véritables mécanismes physiques qui gouvernent la glissance de la glace :

1. La préfusion de surface (Surface Premelting) et la couche quasi-liquide nanométrique (Quasi-Liquid Layer, QLL) (Michael Faraday 1850, Woldemar Weyl 1951, Neville Fletcher 1968, J. G. Dash) ;
2. L’échauffement frictionnel (Frictional Heating) et la fusion dans l’épaisseur de peau thermique (Frank Philip Bowden & Thomas P. Hughes 1939) ;
3. La nanotribologie moderne du XXIe siècle : le roulement thermiquement activé des molécules d’eau de surface faiblement liées agissant comme des « roulements à billes moléculaires » (Weber & Bonn et al., 2018) et la nanorhéologie viscoélastique complexe ($G' \approx G''$) du film interfacial présentant une viscosité effective 50 à 100 fois supérieure à celle de l’eau liquide ordinaire (Canale & Bocquet et al., 2019).

1. L’anomalie tribologique de la glace et l’architecture multi-échelle

Comparons le coefficient de frottement cinétique d’Amontons-Coulomb $\mu_k = F_f / F_N$ pour différentes interfaces solides :

À $-100^\circ\text{C}$ ($173\,\text{K}$), la glace n’est absolument pas glissante ($\mu_k \approx 0.5$). Sa super-lubrification n’apparaît qu’au-dessus d’environ $-35^\circ\text{C}$ selon l’architecture suivante :

  flowchart TD
    A["Question centrale : Pourquoi la glace est-elle si glissante sous 0°C ?<br/>(Coefficient de frottement cinétique μ_k ≈ 0.005–0.02)"] --> B["Hypothèse 1 : Fusion sous pression<br/>(J. Thomson 1849 / J. Joly 1886)"]
    A --> C["Mécanisme 2 : Préfusion de surface & QLL<br/>(M. Faraday 1850 / W. Weyl 1951 / N. Fletcher 1968)"]
    A --> D["Mécanisme 3 : Échauffement frictionnel & Peau thermique<br/>(F. P. Bowden & T. P. Hughes 1939)"]
    A --> E["Mécanisme 4 : Nanotribologie & Viscoélasticité du XXIe siècle<br/>(Weber & Bonn 2018 / Canale & Bocquet 2019)"]

    B --> B1["Équation de Clausius-Clapeyron :<br/>dT_m / dP ≈ -0.0744 K/MPa"]
    B1 --> B2["La pression d'un patin P ≈ 1.6 MPa n'abaisse T_m que de ΔT_m ≈ -0.12°C !<br/>Même à la limite d'élasticité plastique (20 MPa), ΔT_m ≈ -1.5°C max<br/>[CONCLUSION : Réfuté quantitativement]"]

    C --> C1["Molécules H2O de surface sous-coordonnées (liaisons H manquantes)<br/>et paramètre d'étalement positif Δγ = γ_sv - γ_sl - γ_lv > 0"]
    C1 --> C2["Une couche quasi-liquide (QLL) d'épaisseur d(T) ≈ 1–30 nm<br/>existe spontanément à l'équilibre thermodynamique jusqu'à ~ -35°C"]

    D --> D1["Le flux de puissance frictionnelle q = μ P v est piégé dans<br/>une fine épaisseur de peau thermique δ_th = √(π α t_c) ≈ 10–20 μm"]
    D1 --> D2["Échauffement local des aspérités ΔT_fric ~ 5–20 K,<br/>portant l'interface à 0°C à grande vitesse"]

    E --> E1["Mobilité de surface d'Arrhenius (E_a ≈ 11.5 kJ/mol) :<br/>les molécules faiblement liées roulent comme des billes"]
    E1 --> E2["Nanorhéologie viscoélastique complexe (G' ≈ G'') :<br/>viscosité 50–100x celle de l'eau empêchant l'expulsion du film"]

    C2 --> F["Synthèse tribologique multi-échelle moderne :<br/>QLL & roulements moléculaires dominent au repos/basse vitesse ;<br/>chauffage frictionnel & EHL viscoélastique dominent à haute vitesse"]
    D2 --> F
    E2 --> F

2. Réfutation rigoureuse de la fusion sous pression par l’équation de Clapeyron

2.1 Structure cristalline de la glace hexagonale ($I_h$) et anomalie volumique

Dans le réseau cristallin de la glace hexagonale ordinaire (Glace $I_h$), chaque molécule $\text{H}_2\text{O}$ obéit aux règles de Bernal-Fowler et forme quatre liaisons hydrogène directionnelles ($\text{O}-\text{H}\cdots\text{O}$ de longueur $\approx 2.76\,\text{\AA}$) selon une géométrie tétraédrique parfaite ($109.5^\circ$). Ce réseau ouvert comporte de vastes cages interstitielles vides. Lors de la fusion à $T_m = 273.15\,\text{K}$ ($0^\circ\text{C}$), l’effondrement partiel du réseau tétraédrique augmente la masse volumique $\rho$ et diminue le volume massique $v = 1/\rho$ :

• Glace solide $I_h$ à $0^\circ\text{C}$ : $\rho_{\text{ice}} = 916.7\,\text{kg/m}^3 \implies v_{\text{ice}} = 1.0909 \times 10^{-3}\,\text{m}^3/\text{kg}$
• Eau liquide à $0^\circ\text{C}$ : $\rho_{\text{liquid}} = 999.84\,\text{kg/m}^3 \implies v_{\text{liquid}} = 1.0002 \times 10^{-3}\,\text{m}^3/\text{kg}$
• Variation de volume massique à la fusion : $$\Delta v_{\text{fus}} = v_{\text{liquid}} - v_{\text{ice}} = (1.0002 - 1.0909) \times 10^{-3}\,\text{m}^3/\text{kg} \approx -9.1 \times 10^{-5}\,\text{m}^3/\text{kg} < 0$$

2.2 Dérivation de l’équation de Clausius-Clapeyron et calcul numérique

Le long de la courbe de coexistence solide-liquide $P(T_m)$, l’égalité des enthalpies libres massiques $g_{\text{ice}}(T_m, P) = g_{\text{liquid}}(T_m, P)$ et la relation fondamentale $dg = -s\,dT + v\,dP$ imposent :

En introduisant l’entropie massique de fusion $\Delta s_{\text{fus}} = s_{\text{liquid}} - s_{\text{ice}} = \frac{L_f}{T_m}$ avec la chaleur latente de fusion $L_f = 3.34 \times 10^5\,\text{J/kg}$, on obtient l’équation de Clausius-Clapeyron :

En inversant cette relation, la pente d’abaissement du point de fusion par unité de pression vaut :

Il faut donc exercer une pression gigantesque de $13.44\,\text{MPa}$ (environ $133\,\text{atmosphères}$) pour abaisser le point de fusion de la glace d’un seul degré Celsius !

2.3 Pression exercée par un patineur et limite de fluage plastique de la glace

Considérons un patineur de masse $m = 75\,\text{kg}$ ($F_N = m g = 735.75\,\text{N}$) en appui sur une seule lame de longueur de contact $L = 30\,\text{cm} = 0.30\,\text{m}$ et d’épaisseur $w = 1.5\,\text{mm} = 1.5 \times 10^{-3}\,\text{m}$ ($A = L w = 4.5 \times 10^{-4}\,\text{m}^2$) :

Le point de fusion ne s’abaisse que de $0.12^\circ\text{C}$ ! Même en tenant compte de la concentration des contraintes sur les micro-aspérités de surface, la pression locale $P_{\text{local}}$ ne peut jamais dépasser la limite d’écrasement plastique (dureté à l’indentation) de la glace, qui vaut $\sigma_{y,\text{ice}} \approx 10\text{–}20\,\text{MPa}$. Même à cette limite physique absolue de $20\,\text{MPa}$, l’abaissement du point de fusion atteint au maximum $\Delta T_{m,\max} \approx -1.49^\circ\text{C}$. Cela ne peut en aucun cas expliquer pourquoi les hockeyeurs et skieurs glissent à $-15^\circ\text{C}$, ni pourquoi un piéton en bottes plates ($P \approx 34\,\text{kPa}$, $\Delta T_m \approx -0.0026^\circ\text{C}$) ou un palet de hockey glisse sans effort !

3. Le regel de Faraday et l’expérience du fil traversant la glace

Dans l’expérience classique de Bottomley (1872) où un fil lesté traverse un bloc de glace à $0^\circ\text{C}$ sans le couper en deux, la fusion sous le fil exige l’absorption continue de la chaleur latente $L_f = 3.34 \times 10^5\,\text{J/kg}$. Cette chaleur provient du regel de l’eau au-dessus du fil (à $T_0 = 0^\circ\text{C}$) et doit être conduite thermiquement à travers l’épaisseur du fil vers le bas (à $T_0 - |\Delta T_m|$) selon la loi de Fourier (Nye 1967) :

Un fil de cuivre ($\kappa_{\text{Cu}} \approx 400\,\text{W/(m}\cdot\text{K)}$) conduit rapidement cette chaleur latente et traverse le bloc à $v_{\text{wire}} \sim 10^{-5}\text{–}10^{-4}\,\text{m/s}$ (ce qui reste 5 à 6 ordres de grandeur trop lent pour un patin filant à $10\,\text{m/s}$), tandis qu’un fil de nylon de même diamètre ($\kappa_{\text{nylon}} \approx 0.25\,\text{W/(m}\cdot\text{K)}$, $1600$ fois moins conducteur) échoue totalement à traverser la glace malgré une pression identique ! Rappelons d’ailleurs que dès 1850, Michael Faraday avait rejeté l’explication par la pression de James Thomson et proposé avec une intuition géniale que la glace possède naturellement un film superficiel quasi-liquide.

4. Préfusion de surface et couche quasi-liquide nanométrique (QLL)

1. Sous-coordination des molécules de surface et critère de Lindemann : À l’interface cristal-vapeur, les molécules d’eau perdent leurs voisines supérieures et ne forment plus que 2 ou 3 liaisons hydrogène au lieu de 4. Leur constante de rappel $k_{\text{surf}}$ étant réduite de moitié, l’amplitude quadratique moyenne de leurs vibrations thermiques $\sqrt{\langle u_{\text{surf}}^2 \rangle} \propto \sqrt{k_B T / k_{\text{surf}}}$ franchit le seuil d’instabilité de Lindemann bien en dessous de $0^\circ\text{C}$ (dès $-35^\circ\text{C}$), engendrant une couche quasi-liquide (QLL) d’épaisseur $d(T) \sim 1\text{–}30\,\text{nm}$.
2. Modèle thermodynamique de mouillage interfacial (Weyl, Fletcher, Dash) : Comme le coefficient d’étalement $\Delta \gamma = \gamma_{sv} - \gamma_{sl} - \gamma_{lv} \approx 109 - 30 - 75.6 = +3.4\,\text{mJ/m}^2 > 0$ est positif, l’énergie libre interfaciale par unité de surface : $$\Delta G_{\text{surf}}(d) = -\left(\gamma_{sv} - \gamma_{sl} - \gamma_{lv}\right)\left(1 - e^{-d/\xi}\right) + \rho_l L_f \left(\frac{T_m - T}{T_m}\right) d$$ admet un minimum à l’équilibre ($\frac{\partial \Delta G_{\text{surf}}}{\partial d} = 0$) donnant l’épaisseur de la QLL : $$d(T) = \xi \ln\!\left(\frac{\Delta \gamma \, T_m}{\xi \rho_l L_f (T_m - T)}\right) \quad \text{ou} \quad d(T) \propto \left(T_m - T\right)^{-1/3}$$

5. Échauffement frictionnel et couche limite thermique (Bowden & Hughes, 1939)

À grande vitesse ($v \sim 5\text{–}20\,\text{m/s}$), la rugosité des lames et des skis ($R_q \sim 0.1\text{–}1\,\mu\text{m}$) dépasse largement l’épaisseur nanométrique de la QLL statique. En 1939, au laboratoire glaciaire du Jungfraujoch, Frank Philip Bowden et Thomas P. Hughes ont démontré que des skis en laiton ($\kappa \approx 110\,\text{W/(m}\cdot\text{K)}$) évacuent la chaleur de frottement et glissent très mal, tandis que des skis isolants en bois, ébonite ou PTFE ($\kappa \approx 0.2\,\text{W/(m}\cdot\text{K)}$) confinent la chaleur à la surface de la glace et réduisent drastiquement le frottement.

Pendant le temps de contact $t_c = L/v$, la puissance frictionnelle $q_{\text{fric}} = \mu_k P v$ reste confinée dans une mince épaisseur de peau thermique $\delta_{\text{th}} = \sqrt{\pi \alpha_{\text{ice}} t_c} = \sqrt{\pi \alpha_{\text{ice}} L / v} \approx 10\text{–}20\,\mu\text{m}$, produisant une élévation de température :

Cela explique à la fois la génération dynamique d’un film d’eau de plusieurs dizaines de nanomètres à $-10^\circ\text{C}$, le comportement abrasif de la neige semblable à du sable sec en dessous de $-35^\circ\text{C}$ en Antarctique (l’expédition tragique de Scott en 1912), et la remontée du frottement à $0^\circ\text{C}$ par succion capillaire des ménisques d’eau ($\Delta P_{\text{cap}} = -2\gamma_{lv}\cos\theta / h$).

6. Nanotribologie du XXIe siècle : Roulements à billes moléculaires et rhéologie viscoélastique (2018–2024)

1. Roulements à billes moléculaires et loi d’Arrhenius (Weber & Bonn et al., 2018) : Le coefficient de frottement de la glace entre $-100^\circ\text{C}$ et $0^\circ\text{C}$ suit une loi d’activation thermique d’Arrhenius $\mu(T) \propto \exp(E_a / k_B T)$ avec une énergie d’activation $E_a \approx 11.5\,\text{kJ/mol}$, correspondant à la rotation et au saut des molécules d’eau superficielles faiblement liées qui roulent sur le réseau cristallin comme des roulements à billes moléculaires.
2. Nanorhéologie viscoélastique complexe du film interfacial (Canale, Siria & Bocquet et al., Phys. Rev. X 2019) : L’eau liquide pure ($\eta_0 \approx 1.8\,\text{mPa}\cdot\text{s}$) est un très mauvais lubrifiant qui serait instantanément expulsé (squeeze-out) sous la pression mégapascal d’une lame de patin. Grâce à un microscope à force atomique à diapason de quartz (Stroke-Probe AFM), l’équipe de Lydéric Bocquet à l’ENS Paris a découvert que le film nanométrique cisaillé sur la glace présente une viscosité effective 50 à 100 fois supérieure à celle de l’eau ($\eta_{\text{eff}} \sim 0.05\text{–}0.2\,\text{Pa}\cdot\text{s}$, comparable à une huile moteur) ainsi qu’un fort module élastique de conservation $G' \approx G'' \sim 10^4\text{–}10^5\,\text{Pa}$ ! Cette « graisse viscoélastique auto-cicatrisante » génère une portance élastohydrodynamique qui empêche l’expulsion du film sous charge tout en offrant une résistance au cisaillement ultra-faible.