Physique
Un fluide n'a pas de forme propre : sous la moindre contrainte tangentielle, il se déforme indéfiniment. Il ne reste donc qu'un champ de vitesse, auquel on applique les trois conservations — masse, quantité de mouvement, énergie. Le terme d'inertie (u·∇)u, non linéaire, fait tout le reste : la beauté des tourbillons et la difficulté du domaine. Les cinq planches simulent, calculent ou résolvent exactement ; aucune n'est un dessin.
Une vraie simulation, ici, dans le navigateur : méthode de Boltzmann sur réseau, 240 × 80 cellules, neuf directions de propagation par cellule. Rien n'est préenregistré — montez le Reynolds et regardez l'écoulement changer de nature.
Ce qu'il faut guetter. En dessous de Re ≈ 47, le sillage est stationnaire : deux tourbillons restent collés derrière l'obstacle, l'image se fige. Au-delà, cette symétrie devient instable et les tourbillons se détachent alternativement — l'allée de von Kármán. Le simulateur retrouve ce seuil tout seul : il n'est écrit nulle part dans le code, il sort des équations. Le nombre de Strouhal mesuré, St = f·D/U, tourne autour de 0,2, la valeur observée dans la nature sur cinq décades de Reynolds — c'est elle qui donne sa note au vent dans les fils électriques. Il est ici un peu élevé parce que le canal est étroit (D/H = 0,2), ce qui resserre l'écoulement et accélère le lâcher. Au-delà de quelques centaines, la grille ne résout plus les petites échelles : le calcul reste stable, mais il faut le lire comme une caricature de turbulence, pas comme la turbulence.
Re = ρUL/μ compare l'inertie à la viscosité. Deux écoulements de même Reynolds sont dynamiquement semblables, quelles que soient leur taille et leur vitesse : c'est ce qui autorise la maquette en soufflerie.
Deux asymptotes et une chute. À gauche règne la loi de Stokes, CD = 24/Re : la traînée est proportionnelle à la vitesse, l'écoulement est réversible — filmé à l'envers, il reste plausible. À droite, CD se stabilise vers 0,4–0,5 : la traînée devient proportionnelle au carré de la vitesse, car ce qui coûte n'est plus le frottement mais la quantité de mouvement jetée dans le sillage. Entre les deux, aucune formule simple. Et vers Re ≈ 3·10⁵ survient la crise de traînée : la couche limite devient turbulente, colle mieux à la paroi, le sillage se rétrécit et la traînée chute brutalement — c'est exactement à quoi servent les alvéoles d'une balle de golf. Curiosité : en deux dimensions strictes, la loi de Stokes n'existe pas pour un cylindre — c'est le paradoxe de Stokes, aucune solution stationnaire ne satisfait les conditions à l'infini.
Même très faible, la viscosité ne s'ignore jamais contre une paroi : il faut bien que la vitesse y tombe à zéro. Toute la traînée, et tout le décrochage, se jouent dans cette pellicule.
L'équation qui se résout ici même. À droite, le navigateur intègre par tir l'équation de Falkner–Skan, f‴ + ((m+1)/2)·f·f″ + m(1 − f′²) = 0, en cherchant la pente f″(0) qui raccorde le profil à l'écoulement extérieur. Pour β = 0 elle redevient l'équation de Blasius et rend f″(0) = 0,33206, valeur tabulée depuis 1908. Faites glisser β vers les négatifs : le gradient de pression devient adverse, le fluide près de la paroi doit remonter la pente de pression, le profil se creuse, un point d'inflexion apparaît. À β = −0,1988 exactement, f″(0) s'annule : le frottement pariétal disparaît, l'écoulement décolle. C'est ce point qui fait décrocher une aile trop cabrée, et c'est pour le retarder qu'on met des générateurs de vortex sur les ailes.
Non, l'air ne « se rejoint pas » au bord de fuite. La portance vient d'une circulation autour du profil, et sa valeur est fixée par une seule condition physique : que l'écoulement quitte proprement le bord de fuite.
Le rôle du curseur Γ. L'écoulement potentiel autour d'un profil admet une infinité de solutions, une par valeur de la circulation — toutes mathématiquement valables. Décochez la condition de Kutta et promenez Γ : le point d'arrêt arrière se déplace sur l'extrados, et pour toute valeur sauf une, la vitesse devient infinie au bord de fuite. C'est cette absurdité que la viscosité interdit dans le monde réel. Kutta et Joukowski en tirent la seule valeur admissible, Γ = 4πUa·sin(α + β₀), et avec elle le théorème de portance L′ = ρUΓ, exact et sans approximation. La carte de pression montre le reste : dessus l'écoulement accélère donc la pression tombe, dessous elle monte — Bernoulli, mais Bernoulli conséquence de la circulation, non cause de la portance.
élément signature · turbulence
La turbulence n'a pas d'échelle privilégiée : les grands tourbillons en nourrissent de plus petits, qui en nourrissent de plus petits encore, jusqu'à ce que la viscosité les efface. C'est cette absence d'échelle qui la rend, aujourd'hui encore, hors de portée du calcul.
Le compte qui condamne le calcul direct. Kolmogorov, en 1941, tire de deux hypothèses de similitude le spectre E(k) ∝ ε2/3k−5/3 et l'échelle où la viscosité reprend la main, η = (ν³/ε)1/4, soit η/L = Re−3/4. Pour tout simuler sans rien modéliser, il faut donc (L/η)³ = Re9/4 points de grille — et autant de pas de temps. Poussez le curseur jusqu'au Reynolds d'une aile d'avion et lisez le nombre : il dépasse de plusieurs ordres de grandeur tout ce que l'humanité a jamais fait tourner. C'est la raison d'être des modèles de turbulence, et la raison pour laquelle l'existence même de solutions régulières aux équations de Navier–Stokes en trois dimensions reste l'un des sept problèmes du prix du millénaire.