Mécanique des fluides

La mécanique des fluides est un domaine de la physique dédié à l’étude du comportement des fluides (liquides, gaz et plasmas) et des forces internes associées. C’est une branche de la mécanique des milieux continus qui modélise la matière à l’aide de particules assez petites pour relever de l’analyse mathématique mais assez grandes par rapport aux molécules pour être décrites par des fonctions continues.

Elle se divise en deux parties, la statique des fluides qui est l’étude des fluides au repos et la dynamique des fluides, qui est l’étude des fluides en mouvement.

Aujourd’hui, la dynamique des fluides est un domaine actif de la recherche avec de nombreux problèmes non résolus ou partiellement résolus. Elle utilise systématiquement des méthodes numériques appelées « mécanique des fluides numérique » (MFN), ou en anglais computational fluid dynamics (CFD).

Historique

Jusqu'au XVIII^e siècle

Avant qu'elle ne soit étudiée, la mécanique des fluides a été largement employée pour des applications quotidiennes comme l'irrigation en agriculture, les canaux, les fontaines, etc. La sédentarisation des humains a entraîné la nécessaire invention de moyens de maîtrise de l'eau. L'irrigation à petite échelle serait née vers 6500 av. J.-C. à la fin du Néolithique. On commence à trouver de grands ouvrages hydrauliques (canaux, irrigation gravitaire) vers 3000 av. J.-C. Vers cette époque des instruments ont déjà été inventés pour mesurer le niveau des crues, des zones de marécages sont drainées et asséchées, barrages et digues pour se protéger des crues sont construits sur le Nil, le Fleuve Jaune et l'Euphrate. Il est possible que les plus vieux aqueducs aient été construits en Crète au II^e millénaire av. J.-C. et en Palestine au XI^e siècle av. J.-C.

L’étude de l'eau et de son comportement mécanique ne passe des applications concrètes à la théorie que tardivement. À Alexandrie au III^e siècle av. J.-C., Archimède étudie avec les disciples d'Euclide et, en revenant à Syracuse, formule des principes qui sont à l’origine de la statique des fluides notamment avec son principe éponyme. Héron d'Alexandrie au I^er siècle a poursuivi le travail de statique des fluides en découvrant le principe de la pression et surtout du débit.

Durant l'antiquité tardive les grands travaux hydrauliques se poursuivent et se raffinent avec des aqueducs, des systèmes de distribution et d'assainissement de l'eau, mais aussi les fontaines et les bains. Ces travaux sont décrits par Frontin. Comme la plupart des sciences, l'hydrostatique et l'hydraulique disparaissent en partie de l'Europe pendant le moyen âge, la migration du savoir se faisant de l'ancien empire gréco-romain vers l'empire arabe. L'âge d'or islamique voit d'abord la traduction des œuvres d'Archimède, d'Euclide, et la publication du Livre des mécanismes ingénieuxLivre des mécanismes ingénieux ou Kitāb al-Ḥiyal, ouvrage traitant de l'hydraulique et de l'hydrostatique d'Archimède.

Du point de vue des édifices hydrauliques, si le Moyen Âge voit la disparition du système d'irrigation de la Mésopotamie à cause des invasions mongoles provoquant l'effondrement de la population locale, au VII^e siècle sous la dynastie Sui s'achève la première étape des travaux du Grand Canal qui relie Nord et Sud de la Chine.

La mécanique des fluides n'est étudiée à nouveau en Europe qu'avec les études de Léonard de Vinci au XV^e siècle qui décrit à la fois les multiples types d'écoulements et formule le principe de conservation de la masse ou principe de continuité, prenant ainsi la suite de Héron. C'est lui qui jette les fondements de la discipline et introduit de nombreuses notions d'hydrodynamiques dont les lignes de courant. Comprenant intrinsèquement la problématique de résistance à l'écoulement, il conçoit le parachute, l'anémomètre et la pompe centrifuge.

Époque moderne

Personnalités de l'histoire de la mécanique des fluides (liste non exhaustive).

Il faut attendre l'inclusion des mathématiques à la physique pour que la mécanique des fluides gagne en profondeur. En 1738 Daniel Bernoulli établit des lois applicables aux fluides non visqueux en utilisant le principe de conservation de l'énergie mécanique. La naissance du calcul différentiel permet à Jean le Rond D'Alembert en 1749 d'exposer, en 137 pages, les bases de l'hydrodynamique en présentant le principe de la pression interne d'un fluide, du champ de vitesse et des dérivées partielles appliquées aux fluides. Leonhard Euler complète plus tard l'analyse de D'Alembert sur la pression interne et les équations de dynamique des fluides incompressibles.

En 1755 Euler publie ainsi un traité qui donne les équations à dérivées partielles décrivant les fluides parfaits incompressibles. Un peu avant, en 1752, D'Alembert relève le paradoxe à son nom qui montre que les équations contredisent la pratique : un corps plongé dans un fluide se mouvrait sans résistance d'après la théorie, ce que l'observation contredit directement. L'introduction par Henri Navier en 1820 de la notion de frottement sous forme d'un nouveau terme dans les équations mathématiques de mécanique des fluides. George Gabriel Stokes aboutit en 1845 à une équation permettant de décrire un écoulement de fluide visqueux. Les équations de Navier-Stokes marqueront toute la suite de l'histoire de la mécanique des fluides.

Cette suite prend corps dans la seconde moitié du XVIII^e siècle et la première du XX^e siècle :

• développements dans les domaine incompressible ou compressible avec la création du concept de couche limite par Ludwig Prandtl qui s'avèrera très fructueux, particulièrement pour l'aérodynamique et l'hydrodynamique navale,
• étude du domaine nouveau que constitue le supersonique,
• étude des écoulements en milieu poreux par Henry Darcy et des interfaces eau-air par Moritz Weber,
• l'étude des instabilités et de la turbulence, un chapitre bien loin d'être clos aujourd'hui. Ce domaine voit l'apparition d'écoles fondées par un précurseur : Prandtl à Göttingen ou l'école russe de mathématiciens par Kolmogorov.

Au cours de cette période un nouveau chapitre est ouvert par Ludwig Boltzmann avec la description statistique des gaz au niveau microscopique. Ce domaine sera développé par Martin Knudsen pour le domaine inaccessible à une description relevant de l'hypothèse du continu. David Enskog et Sydney Chapman montreront comment passer pour les gaz du niveau moléculaire au continu, permettant ainsi le calcul les coefficients de transport (diffusion, viscosité, conduction) à partir du potentiel d'interaction moléculaire.

Toutes les travaux théoriques s'appuient sur les travaux fondamentaux antérieurs de mathématiciens comme Leonhard Euler, Augustin Louis Cauchy ou Bernhard Riemann.

Par ailleurs le développement de nombreuses installations d'essai et de moyens de mesure permet d'obtenir de nombreux résultats. Tous ne sont pas explicables par la théorie et on voit apparaître un grand nombre de nombres adimensionnels permettant une explication et une justification d'essais effectués sur maquette en soufflerie ou bassin de carène. Deux mondes scientifiques se côtoient et très souvent s'ignorent jusqu'à la fin du XIX^e siècle. Ce gap disparaîtra sous l'impulsion de gens comme Theodore von Kármán ou Ludwig Prandtl au début du XX^e siècle.

Époque récente

Le calcul numérique naît dans la seconde moitié du XX^e siècle. Il va permettre l'éclosion d'une nouvelle branche de la mécanique des fluides, la mécanique des fluides numérique. Elle est basée sur l'avènement de calculateurs toujours plus puissants mais aussi de méthodes mathématiques permettant le calcul numérique. La puissance de calcul permet la réalisation d'« expériences numériques » qui concurrencent les moyens d'essai ou permettent l'interprétation plus aisée de ceux-ci. Ce type d'approche est couramment utilisée dans l'étude de la turbulence.

Le second fait d'importance dans cette période est l'augmentation considérable du nombres de personnes impliquées dans la recherche et développement. Les découvertes sont devenues plutôt le fait d'équipes que d'individus.

Statique des fluides

L'hydrostatique, ou statique des fluides, est l'étude des fluides immobiles. Ce domaine a de nombreuses applications comme la mesure de pression et de masse volumique.

Échelles et nature du problème hydrodynamique

Niveau microscopique

Au niveau le plus bas de la modélisation on décrit le milieu par position et vitesse de chaque particule constitutive et le potentiel d'interaction entre elles. Cette approche est bien sûr limitée par la quantité d'information qu'elle suppose. Elle est utilisée :

• en pratique dans les méthodes de dynamique moléculaire où elle constitue une véritable expérience numérique possible pour un liquide comme pour un gaz,
• en théorie pour des tentatives de construction ab initio d'un système formel de description macroscopique du milieu. Ce type d'approche est extrêmement difficile et peu de résultats ont été obtenus depuis les travaux de Jean Leray. En particulier l'existence de solutions régulières des équations de Navier-Stokes fait l'objet du prix Clay.

Pour les gaz et à un niveau moins détaillé on se contente de décrire la distribution statistique des vitesses et éventuellement de tous les autres degrés de liberté (énergie interne, rotation et vibration dans le cas de molécules). Ludwig Boltzmann a ainsi réussi à écrire l'équation cinétique qui porte son nom. Cette fonction du temps, de la position et de la vitesse peut être calculée à partir d'outils comme la simulation directe Monte Carlo ou la méthode de gaz sur réseau particulièrement bien adaptée aux milieu poreux. Il s'agit de calculs coûteux en raison de la dimension 7 du problème. Pour cette raison on utilise généralement un potentiel d'interaction peu réaliste physiquement mais conduisant à des résultats acceptables.

Niveau mésoscopique

Par ce vocable on entend la description de phénomènes descriptibles à une échelle grande devant la précédente mais petite devant l'échelle du continu.

Concept de particule élémentaire du fluide

La particule fluide décrit un fluide à l'échelle mésoscopique : c'est un volume de dimension suffisamment petite pour que les propriétés du fluide ne varient pas spatialement dans la particule et suffisamment grand pour qu'une quantité importante de molécules soient comprises dedans de manière à moyenner les fluctuations statistiques.

On peut effectuer dans cette particule un bilan de masse, de quantité de mouvement et d'énergie en utilisant les flux correspondants sur les limites du domaine. Cette approche conduit à l'écriture des équations de conservation correspondantes et, par passage à la limite, aux équations descriptives du phénomène. Cette méthode est aussi la base de la description numérique, le volume élémentaire étant alors la maille élémentaire du calcul.

Suppression des détails de taille intermédiaire

La géométrie étudiée peut comprendre des détails dont la prise en compte explicite va rendre le problème coûteux, par exemple une rugosité de la surface ou le détail de la géométrie d'un milieu poreux. Dans ce dernier cas les méthodes bien connues de la prise de moyenne volumique ou de l'homogénéisation permettent le calcul de quantités intervenant sous forme de coefficients comme le coefficient de diffusion dans l'équation de Darcy. Dans le cas d'une rugosité l'homogénéisation aboutit à l'écriture d'une relation de saut à la paroi, c'est-à-dire une relation liant toute valeur à sa dérivée spatiale.

On peut faire également entrer dans cette catégorie les phénomènes de raréfaction dans un choc ou une couche pariétale. Dans ces régions d'espace les équations du continu sont invalides sur une distance de quelques libres parcours moyens. On peut généralement les ignorer. Lorsque ce n'est pas les cas leur modélisation aboutit comme précédemment à des équations de saut. Les relations de Rankine-Hugoniot en sont un exemple.

Enfin, et ce n'est pas le moindre problème, on peut faire disparaître toutes les fluctuations d'un écoulement turbulent par des méthodes de moyennage très diverses, pouvant ramener le problème à une simple diffusion équivalente. Là aussi le but est de simplifier le calcul, possible par la simulation directe, mais coûteux.

Niveau macroscopique

Le niveau macroscopique résulte donc d'une simplification drastique de tous les détails du problème, lesquels sont tout de même présents au travers des coefficients qui interviennent dans les équations descriptives, des conditions aux limites et de l'équation d'état du milieu.

Compressible et incompressible

Ces notions qui séparent nettement deux types d'écoulements ont une origine microscopique :

• le caractère compressible généralement associé à un gaz est lié au fait qu'un tel milieu est formé d'objets très espacée ayant des interactions rares, caractérisées par un potentiel particule-particule. Ceci est vrai même dans le cas de milieux contenant des espèces chargées en faible proportion, où les électrons ne sont pas totalement libres et accompagnent (statistiquement) les ions (diffusion ambipolaire). La connaissance de ces potentiels, aujourd'hui d'origine spectroscopique, est suffisante pour permettre le calcul de toutes les propriétés de transport du milieu : coefficients de diffusion binaire et thermique (équations de Stefan-Maxwell), viscosités dynamique et volumique, conductivité. Ce caractère de milieu peu dense n'est pas affecté par un changement de pression donc une variation du libre parcours moyen entre deux collisions.
• le caractère d'incompressibilité associé aux liquides est lié aux liaisons que voit une particule dans un tel milieu. Elle est en effet liée à plusieurs voisins, même si ces liaisons ne sont pas aussi strictes que dans un solide. Ce caractère interdit une approche formelle comme dans les gaz : les propriétés de transport sont mesurées, la théorie ne permettant que d'expliquer les variations avec la température par exemple. Ce caractère d'incompressibilité n'est pas insurmontable : une pression très élevée de quelques centaines de GPa tel que rencontré dans le noyau terrestre met en évidence une variation de masse volumique des composants liquides.

Équations de la mécanique des fluides

On considère ici les équations de Navier-Stokes pour un fluide simple (newtonien), qui sont la pierre angulaire du domaine et à partir desquelles on déduit de nombreuses autres lois.

Ces équations sont écrites dans un repère fixe, mais il existe deux façons d'exprimer les différentes grandeurs en fonction de la position : soit en fonction des coordonnées actuelles dans le repère (description eulérienne), soit en fonction des coordonnées occupées à un certain instant initial (description lagrangienne). Dans le premier cas le vecteur ${\vec {v}}(x,y,z,t)$ représente la vitesse à l'instant t et au point de coordonnées ($x,y,z$) (mais à différents instants il ne s'agira pas de la même portion de matière), dans le second cas ${\vec {v}}(x_{0},y_{0},z_{0},t)$ représente la vitesse à l'instant t de la matière qui à l'instant initial occupait la position $(x_{0},y_{0},z_{0})$ (et qui à l'instant t se trouve en un point différent $(x(x_{0},y_{0},z_{0},t),y(x_{0},y_{0},z_{0},t),z(x_{0},y_{0},z_{0},t))$). On utilise le plus souvent la description eulérienne.

Équations de base

On peut obtenir ces équations par au moins deux voies :

• à partir des relations de conservation de la masse, de la quantité de mouvement et de l'énergie,
• à partir de l'équation de Boltzmann décrivant l'évolution moléculaire par la méthode de Chapman-Enskog. Cette méthode n'est utilisable pour les gaz en raison de la simplicité relative des interactions au niveau microscopique dans ce cas.

Dans la première méthode apparaissent le tenseur des contraintes (ou tenseur de pression, incluant contraintes visqueuses et pression) et le flux de chaleur. Pour ces deux quantités on fait l'hypothèse qu'elles sont liées à un gradient :

• le flux de chaleur est proportionnel au gradient de température (loi de Fourier),
• le tenseur des contraintes est proportionnel au tenseur des vitesses de déformations (hypothèse de Stokes). En une dimension d'espace cette expression s'exprime en disant que la contrainte visqueuse est proportionnelle au taux de cisaillement. Ceci définit un écoulement « newtonien ».

Le mécanisme sous-jacent dans les deux cas n'est pas très apparent : on se doute que cette proportionnalité est liée à une linéarisation des équations qui décrivent le problème exact sous-jacent. C'est là un processus général en physique mathématique.

La méthode partant du microscopique permet d'éclairer cet aspect. Les équations de Navier-Stokes sont l'expression d'une petite perturbation de la fonction de distribution microscopique des vitesses et, éventuellement, des énergies internes (statistique de Maxwell-Boltzmann). A contrario les équations d'Euler décrivent le cas correspondant à l'équilibre thermodynamique local.

Il faut alors donner les coefficients qui interviennent : pression, viscosité et conductivité. La pression est définie par l'équation d'état. Les propriétés de transport, viscosités, conductivité peut résulter dans le cas du gaz d'un calcul effectué à partir du niveau microscopique (du potentiel interatomique). Pour les liquides ces quantités relèvent de l'expérience.

Exemple : fluide incompressible

• étape 1 : équations de conservation

équations de conservation pour la masse et la quantité de mouvement

$\mathbf {\nabla } \cdot \mathbf {V} =0$

${\begin{array}{rcl}\rho {\dfrac {\partial \mathbf {V} }{\partial t}}+\rho \mathbf {\nabla } \cdot \left(\mathbf {V} \mathbf {V} \right)&=&\mathbf {\nabla } \cdot {\mathsf {P}}\end{array}}$

• étape 2 : loi de comportement

${\mathsf {P}}=\underbrace {\,\mu \,\mathbf {\nabla } ^{2}\mathbf {V} \,} _{cisaillement}-\underbrace {\,p\,{\mathsf {I}}\,} _{pression}$

• étape 3 : équation d'état et expression de la viscosité

$p=...\,,\;\;\;\mu =...$

avec

ρ	masse volumique
V	vitesse
t	temps
P	tenseur de pression (contraintes)
I	tenseur unité
p	pression
μ	viscosité dynamique

 

Similitude

La similitude est la mise en évidence de nombres sans dimensions permettant de réduire le nombre de paramètres intervenant dans les équations afin de simplifier son analyse, éventuellement de définir des expériences à l'échelle du laboratoire. Elle est basée sur l'invariance d'échelle qui assure la covariance des équations : celles-ci sont valides dans tout référentiel galiléen.

On peut alors par un changement de variable faire apparaître des nombres adimensionnels et diminuer ainsi le nombre de variables d'un problème.

Exemple : nombre de Reynolds

Reprenons l'exemple précédent. On définit :

• une longueur de référence L^* caractéristique du problème traité,
• une vitesse de référence V^*, par exemple issue de la condition aux limites.

À partir de ces valeurs on déduit les variables réduites :

- espace	${\tilde {\mathbf {x} }}={\frac {\mathbf {x} }{L^{*}}}$
- temps	${\tilde {t}}={\frac {V^{*}}{L^{*}}}\,t$
- vitesse	${\tilde {\mathbf {V} }}={\frac {\mathbf {V} }{V^{*}}}$
- pression	${\tilde {p}}={\frac {p}{\rho {V^{*}}^{2}}}$

Le système en variables réduites s'écrit :

$\mathbf {\nabla } _{\tilde {x}}\cdot {\tilde {\mathbf {V} }}=0$

${\frac {\partial {\tilde {\mathbf {V} }}}{\partial {\tilde {t}}}}+\mathbf {\nabla } _{\tilde {x}}\cdot \left({\tilde {\mathbf {V} }}{\tilde {\mathbf {V} }}\right)=-\mathbf {\nabla } {\tilde {p}}+{\frac {1}{Re}}\nabla _{\tilde {x}}^{2}\,{\tilde {\mathbf {V} }}$

$\mathbf {\nabla } _{\tilde {x}}=L^{*}\mathbf {\nabla }$ est l'opérateur nabla adimensionné et $Re={\frac {\rho V^{*}L^{*}}{\mu }}$ le nombre de Reynolds.

Le problème ne dépend plus explicitement des dimensions physiques : l'équation ci-dessus décrit une famille de problèmes (et donc de solutions) se déduisant l'un de l'autre par transformation de l'espace et du temps.

 

Instabilités et turbulence

Instabilités

Instabilité de Kelvin-Helmholtz matérialisée par des nuages.

L'instabilité des solutions des équations est due au terme non-linéaire de transport de quantité de mouvement V ⋅ ∇V. Elles correspondent à une bifurcation de la solution obtenue pour une certaine valeur du nombre de Reynolds. On rencontre divers types d'instabilités :

• instabilité de cisaillement bidimensionnelle pour des profils de vitesse perpendiculaire à l'écoulement ayant un point d'inflexion (instabilité de Kelvin-Helmholtz). Le tourbillon généré est dans le plan de l'écoulement ;
• instabilités centrifuges de type Taylor-Couette qui se crée lorsque le moment cinétique r V( r ) décroit lorsque l'on séloigne du centre de courbure. Le tourbillon généré est perpendiculaire à l'écoulement, conduisant par exemple aux tourbillons de Görtler. Il existe nombre d'autres instabilités de type inertiel telles les instabilité elliptique et instabilité de Crow rencontrées en aéronautique ou en géophysique.

De plus les interfaces soumises à une accélération ou à un champ de gravité peuvent être le siège d'instabilités : Rayleigh-Taylor, Richtmyer-Meshkov, etc.

Transition vers la turbulence

Strioscopie de l'écoulement transitionnel créé par une chandelle.

Le passage de l'état laminaire d'un écoulement vers un état totalement turbulent peut emprunter plusieurs voies :

• transition naturelle : une perturbation quelconque est amplifiée comme le montre une analyse de stabilité (équation de Orr-Sommerfeld). Les perturbations d'abord régulières (ondes de Tollmien-Schlichting) se déforment, créent des tourbillons longitudinaux qui sont eux-mêmes déformés et qui finissent par créer des régions (« spots ») turbulentes, lesquelles finissent par occuper tout l'espace.
• transition « by-pass » : ce phénomène, présent dans les couches limites, désigne une transition forcée par la contamination par une turbulence externe. Les premières étapes de la transition naturelle sont contournées, d'où le nom.
• transition par rugosité : les irrégularités de la paroi sont un moteur puissant de déstabilisation de la couche limite.
• transition par décollement de l'écoulement moyen créant une couche de cisaillement instable.

Il n'existe pas de modèle universel de transition. Ceci est aisément compréhensible dans le cas de la transition naturelle où la source de l'instabilité peut être diverse et où de plus son amplitude joue un rôle. De même on ne maîtrise pas forcément une turbulence externe. En pratique on utilise des critères expérimentaux valides sur telle ou telle configuration.

Turbulence

Visualisation par fluorescence induite par laser d'un jet turbulent.

La turbulence est un phénomène étudié depuis Léonard de Vinci mais encore mal compris. Il n'existe pas de théorie permettant de décrire le phénomène à partir des équations de Navier-Stokes. La cascade turbulente se manifeste par un transfert d'énergie des grandes structures créées par les gradients de vitesse - encore le terme V ⋅ ∇V - vers les petits tourbillons détruits par dissipation visqueuse. Un résultat majeur obtenu par Kolmogorov est la description des échelles intermédiaires où se produit une diffusion de l'énergie cinétique par mélange et étirement/repli des tourbillons. Cette région possède une propriété d'auto-similitude : les transferts se produisent identiquement à toutes les échelles. Ce résultat illustre la capacité explicative de l'approche physique statistique et systèmes dynamiques.

Il faut mentionner l'existence d'une turbulence quasi-bidimensionnelle obtenue lorsque l'une des dimensions du problème est limitée : c'est le cas de l'atmosphère où les grands tourbillons excèdent largement la « hauteur utile » où peut se développer une troisième dimension. Il se produit alors une double cascade d'énergie.

En pratique l'approche physique statistique ne permet pas un calcul global. De même la résolution directe des équations est beaucoup trop coûteuse et ne sert qu'à générer des expériences numériques servant de test à une théorie. En pratique la mécanique des fluides numérique utilise une méthode où les moments des corrélations statistiques des variables issus d'une prise de moyenne sont modélisés par une hypothèse physique raisonnable. Il existe plusieurs modèles, chacun étant plus ou moins adaptée à une situation donnée.

Les effets de la turbulence sur l'écoulement sont importants. Directement ils favorisent les échanges de masse, quantité de mouvement et énergie. Ce phénomène augmente également le bruit acoustique. Il a aussi un effet indirect en modifiant la structure globale d'une région, par exemple la région décollée d'une couche limite ou un jet.

Lois de comportement

La loi de comportement d'un milieu solide ou fluide (voire intermédiaire) relie les contraintes σ_ij exercées dans le milieu aux déformations ε_ij du milieu et/ou à leurs dérivées par rapport au temps.

Fluide newtonien

Pour beaucoup de fluides le tenseur des contraintes peut s'écrire comme la somme d'un terme isotrope (la pression p) et d'un déviateur (le cisaillement):

$\varepsilon _{ij}=\underbrace {-p\,\delta _{ij}} _{pression}+\underbrace {\mu \left({\frac {\partial V_{i}}{\partial x_{j}}}+{\frac {\partial V_{j}}{\partial x_{i}}}\right)} _{cisaillement}$

δ_ij est le symbole de Kronecker, μ la viscosité dynamique et V la vitesse.

En réalité il existe toujours un terme de viscosité volumique  μ' div V δ_ij  correspondant à une variation isotrope de volume et dû à des interactions moléculaires inélastiques. Ce terme est généralement négligé quoique mesurable et, dans le cas des gaz, calculable. Très petit, il est supposé nul dans l'hypothèse de Stokes.

Certains matériaux comme les verres ont un comportement qui passe continûment de l'état solide à l'état liquide. C'est vraisemblablement le cas du verre commun si l'on en croit les mesures de viscosité dans la plage où celles-ci sont faisables en un temps raisonnable ou celle du Silly Putty.

Fluides non newtoniens

Versement vertical d'un shampooing qui crée un jet latéral par effet Kaye.

De nombreux fluides ont des comportements différents, particulièrement en cisaillement. Ce comportement est lié à leur composition : phase solide en suspension, polymère, etc. Leur étude relève de la rhéologie. On présente généralement leur comportement sous un cisaillement simple pour lequel la viscosité est la pente de la courbe contrainte-déformation :

• le fluide de Bingham (boue de forage, dentifrice) qui a un comportement visqueux newtonien passé un seuil de déformation correspondant à la dislocation de sa structure au repos.
• fluides rhéofluidifiants ou pseudo-plastiques (sang, peintures, pâte à papier, etc.) dont la viscosité apparente diminue avec la contrainte appliquée, le phénomène étant lié à une diminution des liaisons internes, affectées par l'écoulement. Certains fluides de Bingham ont un comportement pseudo-plastique.
• à l'inverse certains fluides sont rhéoépaississant ou dilatants comme les suspensions concentrées.

La relation containtes-déformation n'est pas suffisante pour caractériser certains fluides dont le comportement est plus complexe :

• dépendance en fonction du temps comme les fluides thixotropes comme le ketchup dont la viscosité apparente diminue avec le temps sous contrainte constante. Ce phénomène est lié à une déstructuration plus ou moins rapide du milieu. Plus rarement on rencontre des fluides antithixotropes comme le latex.
• Un autre caractère possible est l'élasticité de fluides comme certaines résines polyacrilamides capables d'aligner leurs chaînes macromoléculaires dans le sens de l'écoulement.

Ces caractéristiques peuvent donner naissance à des comportement remarquables comme :

• l'effet Weissenberg pour les fluides pseudo-plastiques ;
• l'effet de siphon ouvert pour les fluides élastiques ;
• l'effet Kaye pour les fluides thixotropes.

Les comportements peuvent être décrits par des modèles rhéologiques obtenus en ordonnant de manière plus ou moins complexes des éléments de base : ressort pour l'élasticité, amortisseur pour le comportement visqueux, patin pour la pseudo-plasticité. On obtient ainsi le modèle de Kelvin-Voigt ou le modèle de Maxwell pour décrire la viscoélasticité.

Les caractéristiques sont mesurées à l'aide de rhéomètres ou, dans le cas des polymères, peuvent être prédites.

Types d'écoulement (milieu homogène)

Stationnarité, instationnarité

Allée de tourbillons de Karman.

Un écoulement peut être stationnaire ou instationnaire ou les deux à la fois. Prenons l'exemple de l'écoulement autour d'un cylindre infini :

• à bas nombre de Reynolds basé sur le diamètre l'écoulement est laminaire et stationnaire ;
• lorsque le nombre de Reynolds augmente ( Re ≈ 10 ), on voit apparaître une région de recirculation à l'arrière : cet écoulement est stationnaire ;
• cette recirculation entraîne une instabilité de type Kelvin-Helmholtz ( Re ≈ 100 ) et l'apparition d'allées de Kármán par appariement des tourbillons produits de part et d'autre ;
• à plus grand nombre de Reynolds ( Re ≈ 1000 ) l'écoulement de sillage devient turbulent donc instationnaire mais stationnaire en moyenne statistique.

Vorticité

Tourbillon de bout d'aile.

Les tourbillons peuvent naître dans une région décollée comme la recirculation dans l'exemple précédent. Il s'agit alors d'un phénomène entretenu d'origine visqueuse.

Ils peuvent également avoir pour origine une dissymétrie des conditions aux limites : c'est le cas des extrémités d'une aile d'avion. Dans ce cas il s'agit d'un phénomène inertiel non entretenu (en un point de l'espace donné). Les tourbillons ainsi créés sont de grande taille et peu affectés par la viscosité, ce leur confère une grande durée de vie.

Mathématiquement le tourbillon (ou vorticité) se définit comme le rotationnel de la vitesse ou la moitié de cette valeur. On sait écrire une équation de transport pour cette quantité qui est à la base des études sur la turbulence vue sous l'angle mécanique des fluides et non sous l'angle statistique comme dans l'étude de la cascade turbulente.

Compressibilité

Balle de fusil en vol transsonique.

Tous les fluides sont visqueux jusqu'à un certain degré. La compressibilité de l'eau par exemple vaut environ 5 × 10^-10 m² N^-1, ce qui suppose des pressions de l'ordre du kbar pour obtenir un effet mesurable. Cette faible valeur permet dans le cas général de faire l'approximation de masse volumique constante. Les écoulements dans lesquels cette approximation est valide sont généralement tels que la température y est sensiblement constante et où l'on peut par suite supposer la viscosité constante. L'équation de conservation de l'énergie est découplée et les équations de Navier-Stokes réduites à une forme plus simple. Si de plus on suppose le nombre de Reynolds petit ( Re ≈ 1 ) on aboutit à l'équation de Stokes. Dans le cas d'un écoulement irrotationnel on montre que la vitesse découle d'un potentiel : on parle d'écoulement potentiel.

Toutefois la compressibilité d'un liquide n'est jamais nulle et il est possible d'y propager une onde de choc, laquelle suppose une discontinuité des toutes les variables comme indiqué par les relations de Rankine-Hugoniot. Celles-ci sont relatives aux équations d'Euler, donc à un milieu sans viscosité. Cette discontinuité n'existe qu'au point de vue macroscopique puisque la théorie cinétique montre pour les gaz une variation rapide sans discontinuité sur une distance de quelques libres parcours moyens.

L'onde de choc résulte d'une propriété remarquable des équations d'Euler : leur caractère hyperbolique. L'information dans le milieu est transportée par les caractéristiques. Ceci a donné lieu par le passé à des méthodes de résolution par construction géométrique dans des cas assez simples comme une tuyère ou l'onde accompagnant un objet en vol supersonique . Cette propriété est aujourd'hui à la base des méthodes de résolution numérique par volumes finis : les solveurs de Riemann.

Écoulements visqueux et non visqueux, couche limite

Hors problème de turbulence les effets dits visqueux, en fait tous les effets liés au transport de masse (diffusion), de quantité de mouvement (cisaillement) et d'énergie (conduction), sont généralement confinés à des régions particulières, généralement une paroi et dans ce cas on parle de couche limite. Un immense progrès dans la compréhension de ce phénomène a été faite au début du XX^e siècle. Elle a permis l'avènement de l'aérodynamique moderne grâce à l'analyse que permet son caractère parabolique : l'information ne remonte pas l'écoulement. En outre la relative simplicité des équations autorise la mise en évidence de solutions approchées.

Écoulement en milieu inhomogène

Écoulements polyphasiques

Cavitation générée par une hélice.

Ce domaine de la mécanique des fluides s'intéresse ce qui se passe lorsque l’on a affaire à plusieurs phases qui s’écoulent ensemble. Dans la majorité des cas il s'agit d'un milieu diphasique où une phase mineure en volume est dispersée dans la phase majeure. On peut distinguer en fonction du milieu majoritaire :

• liquide contenant un gaz sous forme de bulles (ébullition, échangeur de chaleur), de liquide (nombreuses applications industrielles, en particulier l'industrie pétrolière), de solide (suspensions diluées) ou de vide (cavitation),
• gaz contenant un liquide sous forme de gouttes (sprays).

Cette systématisation des phénomènes peut faire illusion : cela cache des problèmes de natures très différentes. Par exemple les bulles et leur interaction avec leur environnement constituent à elles seules un vrai problème physique que l'on doit aborder avant même de s'intéresser au problème diphasique.

Pour le traitement théorique et numérique du problème on distingue les méthodes cinétique où l'on suit chaque élément de la phase diluée en lui appliquant les lois d'interaction ad hoc (par exemple dans l'équation de Mason-Weaver) et méthodes bifluides où des équations de Navier-Stokes couplées sont écrites pour chaque phase, moyennant certaines hypothèses sur le moyennage des phases (exemple des méthode « volume de fluide (en) ». Cette méthode est plus économique mais pose souvent des problèmes de conditions aux limites où les hypothèses ne sont pas respectées.

Il faut noter que les systèmes diphasiques sont susceptibles de montrer des instabilités spécifiques, un exemple remarquable étant le geyser.

En taille et fraction suffisante les éléments dispersés peuvent affecter la turbulence.

Écoulements en milieu poreux

Microtomographie d'un composite Ti₂AlC / Al.

Les écoulements en milieu poreux sont présents dans de nombreux domaines comme l'hydrologie, les protections thermiques, etc. Il s'agit souvent de fluides homogènes mais on rencontre des cas hétérogènes comme dans l'extraction pétrolière. Ce sont par nature des écoulements de fluide à faible vitesse, généralement décrits par l'équation de Stokes à l'échelle du pore. La loi de Darcy établie expérimentalement est démontrable par prise de moyenne volumique ou homogénéisation sous cette condition. L'extension à des écoulements plus rapides (loi de Darcy-Forchheimer) se fait en introduisant un nombre de Reynolds. Pour les gaz on sait également traiter tous les régimes d'écoulement depuis le moléculaire jusqu'au continu (équation de Darcy-Klinkenberg).

La quantité importante dans le domaine est la perméabilité. Celle-ci est mesurable. Elle a longtemps été évaluée théoriquement par des modèles utilisant des porosités de forme simple, respectant la porosité (par exemple la loi de Kozeny-Carman). Ces méthodes ont une prédictabilité limitée aux variations et non aux valeurs absolues. Ceci a changé avec l'avènement de la microtomographie qui permet une simulation numérique directe du phénomène à l'échelle du pore.

Mécanique des fluides numérique

La mécanique des fluides numérique consiste à étudier les mouvements d'un fluide, ou leurs effets, par la résolution numérique des équations régissant le fluide. En fonction des approximations choisies, qui sont en général le résultat d'un compromis en termes de besoins de représentation physique par rapport aux ressources de calcul ou de modélisation disponibles, les équations résolues peuvent être les équations d'Euler, les équations de Navier-Stokes, etc.

La mécanique des fluides numérique a grandi d'une curiosité mathématique pour devenir un outil essentiel dans pratiquement toutes les branches de la dynamique des fluides, de la propulsion aérospatiale aux prédictions météorologiques en passant par le dessin des coques de bateaux. Dans le domaine de la recherche, cette approche est l'objet d'un effort important, car elle permet l'accès à toutes les informations instantanées (vitesse, pression, concentration) pour chaque point du domaine de calcul, pour un coût global généralement modique par rapport aux expériences correspondantes.

Domaines d’application

La mécanique des fluides au sens strict a de nombreuses applications dans divers domaines comme l'ingénierie navale, l'aéronautique, l'étude de l'écoulement du sang (hémodynamique), la météorologie, la climatologie ou encore l'océanographie.

Il existe également un grand nombre de domaines plus spécialisés qui peuvent s’écarter de la définition restrictive comme l’électro-fluidodynamique, la microfluidique ou l’étude des écoulements polyphasiques. Elle est actuellement étendue à des écoulements tels que ceux des glaciers ou du manteau terrestre.

Références

Ouvrages de référence

• Historique
  • (en) Olivier Darrigol, Worlds of Flow. A History of Hydrodynamics from the Berboullis to Prandtl., Oxford University Press, 2005 (ISBN 978-0-19-856843-8, lire en ligne)
• Ouvrages généraux
  • Patrick Chassaing, Mécanique des fluides : Élément d'un premier parcours, Cépaduès, 2010 (ISBN 2-854-28929-3)
  • (en) Lev Landau et Evgueni Lifchits, Mécanique des fluides, Ellipses Marketing, 1989 (ISBN 2-729-89463-2)
  • (en) George K. Batchelor, An Introduction to Fluid Dynamics, Cambridge University Press, 2000 (ISBN 0-521-66396-2)
  • (en) Hermann Schlichting et Klaus Gernsten, Boundary Layer Theory, Springer, 1955 (ISBN 978-3-662-52917-1, lire en ligne)
• Passage du microscopique au continu
  • (en) Sydney Chapman et Thomas George Cowling, The Mathematical Theory of Non-uniform Gases, Cambridge University Press, 1991 (ISBN 0-521-40844-X)
  • (en) Joseph Oakland Hirschfelder, Charles Francis Curtiss et Robert Byron Bird, Molecular Theory of Gases and Liquids, John Wiley and Sons, 1966 (ISBN 978-0-471-40065-3)
• Turbulence
  • Christophe Bailly et Geneviève Compte-Bellot, Turbulence, CNRS Editions, 2003 (ISBN 2-271-06008-7)
  • (en) Stephen B. Pope, Turbulent Flows, Cambridge University Press, 2000 (ISBN 978-0-521-59886-6)
• Loi de comportement
  • (en) Rutherford Aris, Vectors, Tensors, and the Basic Equations of Fluid Mechanics., Dover Publications, 1962 (ISBN 0-486-66110-5)
  • Guy Couarraze, Jean-Louis Grossiord et Nicolas Huang, Initiation à la rhéologie, Éditions Lavoisier, 2014 (ISBN 2-743-01568-3)
  • (en) R. Byron Bird, Robert C. Amstrong et Ole Hassager, Dynamics of Polymeric Liquids : Fluid Dynamics, Wiley Intersciences, 1987 (ISBN 0-471-80245-X)
• Milieux poreux
  • (en) Kambiz Vafai (Ed.), Handbook of Porous Media., CRC Press, 2015 (ISBN 978-1-4398-8557-4, lire en ligne)
  • (en) Stephen Whitaker, The Method of Volume Averaging, Kluwer Academic Publishers, 2010 (ISBN 978-3-642-05194-4)
• Diphasique
  • (en) George Yadigaroglu and Geoffrey F. Hewitt (Eds.), Introduction to Multiphase Flow, Springer, 2018 (ISBN 978-3-319-58717-2)
  • (en) Christopher Earls Brennen, Cavitation and Bubble Dynamics, Oxford University Press, 1995 (ISBN 0-195-09409-3)
  • (en) Dimitri Gidaspow, Multiphase Flow and Fluidization, Academic Press, 2012 (ISBN 978-0-1228-2470-8)

Vulgarisation

• Étienne Guyon, Jean-Pierre Hulin et Luc Petit, Ce que disent les fluides : la science des écoulements en images, Belin, 2005 (ISBN 2-70-113557-5)
• Romain Rioboo, « Introduction à la mécanique des fluides », sur YouTube
• Isabelle Gallagher, « Kolmogorov, le spectre de la turbulence », sur Bibliothèque nationale de France

Articles connexes

• Formules de mécanique des fluides
  • Application à un fluide en rotation : Équations primitives atmosphériques
• Hydrostatique
• Dynamique des fluides
• Théorème de Kutta-Jukowski
• Oléohydraulique
• Réseaux hydrauliques
• Thermohydraulique
• Bélier hydraulique

Centres d'études en mécanique des fluides

• Office national d'études et de recherches aérospatiales