Couche limite

    La couche limite est la zone d'interface entre un corps et le fluide environnant lors d'un mouvement relatif entre les deux, conséquence de sa viscosité.
   Elle est un élément important en mécanique des fluides, (aérodynamique, hydrodynamique), en météorologie, en océanographie, etc.
  Lorsqu'un fluide réel s'écoule le long d'une paroi supposée fixe, les vitesses sur la paroi sont nulles alors qu'à l'infini (c'est-à-dire loin de l'obstacle) elles sont égales à la vitesse de l'écoulement non perturbé. 
  Sur une normale à la paroi, la vitesse doit donc dans tous les cas varier entre 0 et un maximum. La loi de variation dépend de la viscosité du fluide qui induit un frottement entre les couches voisines : la couche la plus lente tend à freiner la couche la plus rapide qui, en retour, tend à l'accélérer.
Dans ces conditions, une forte viscosité égalise au maximum les vitesses. 
  Au contraire, si le fluide est peu visqueux, les différentes couches sont beaucoup plus indépendantes : la vitesse à l'infini se maintient jusqu'à une courte distance de l'obstacle et il y a une variation plus forte des vitesses dans la petite épaisseur de la couche limite.
Dans le premier cas, il faut utiliser les équations générales du fluide visqueux. Dans le second, on peut utiliser dans la couche limite des équations simplifiées complétées par des résultats expérimentaux. Les équations, également plus simples, du fluide parfait appliquées au-delà de la paroi « engraissée » par la couche limite fournissent les conditions aux limites pour le calcul.
En fait, ce n'est pas la viscosité elle-même qui intervient.
   Comme toujours en mécanique des fluides, c'est un nombre sans dimension qui caractérise le phénomène : le nombre de Reynolds. Celui-ci décrit le rapport des forces liées à la vitesse aux forces de frottement.
   Ainsi, au lieu d'augmenter la viscosité, on peut obtenir un phénomène semblable en diminuant la vitesse ou les dimensions de l'obstacle.

Théorie


Profil de vitesse au travers d’une couche limite au dessus d’une plaque
   La définition même de la couche limite réside dans le fait qu'elle représente la région de l'écoulement où les effets visqueux sont aussi importants que les effets inertiels (en termes d'ordre de grandeur).
   Ce n'est en effet pas le cas loin de la paroi, où l'écoulement est alors dit « d'Euler », et où les effets visqueux ne se font pratiquement pas ressentir. 
   Un fluide parfait est par définition non conducteur et a ses coefficients de Lamé nuls (c'est-à-dire pas de viscosité).
On définit en général l'épaisseur de la couche limite telle que : u(\Delta(x))=0.99\cdot U_e~
avec :
  • U_e~, la vitesse uniforme de l'écoulement sans obstacle
  • \Delta(x)~, l'épaisseur de la couche limite en fonction de x.
  Le profil de vitesse au sein de la couche limite dépend de sa nature laminaire ou turbulente. L'endroit où se produit la transition entre l'écoulement laminaire et l'écoulement turbulent dépend du nombre de Reynoldsprofil).
   Dans une couche où Re est suffisamment élevé, la couche est turbulente. Si Re < 1, l’équation est linéaire car les phénomènes de diffusion dominent et l'écoulement est laminaire. (Re) et de la géométrie de l'écoulement (

Équations de la couche limite

    La compréhension et la modélisation des équations de la couche limite sont peut-être une des plus importantes avancées de la dynamique des fluides. 
   En utilisant l'analyse d'échelle, les équations de Navier-Stokes peuvent être écrites sous forme simplifiée. En effet, les équations de Navier-Stokes originales sont elliptiques alors que les équations simplifiées sont paraboliques. 
   Cela simplifie grandement la résolution des équations. La simplification repose sur la division en deux de l'espace dans lequel s'écoule le fluide : la couche limite et le reste de l'espace (le reste étant facile à résoudre par de nombreuses méthodes).
     La couche limite est alors gouvernée par des équations différentielles partielles faciles à résoudre. 
  Les équations de Navier-Stokes et de continuité, pour un écoulement bidimensionnel incompressible en coordonnées cartésiennes, sont :
 {\partial u\over\partial x}+{\partial v\over\partial y}=0
 u{\partial u \over \partial x}+v{\partial u \over \partial y}=-{1\over \rho} {\partial p \over \partial x}+{\nu}\left({\partial^2 u\over \partial x^2}+{\partial^2 u\over \partial y^2}\right)
 u{\partial v \over \partial x}+v{\partial v \over \partial y}=-{1\over \rho} {\partial p \over \partial y}+{\nu}\left({\partial^2 v\over \partial x^2}+{\partial^2 v\over \partial y^2}\right)
où u et v sont les composantes de la vitesse, \rho est la masse volumique, p la pression, et \nu est la viscosité cinématique du fluide en un point.
Un écoulement qui a un nombre de Reynolds élevé peut être simplifié.
    La simplification consiste à diviser l'espace en deux régions. La première est la région où l'écoulement du fluide n'est pas affecté par la viscosité — la majorité de l'espace — , l'autre région — proche des surfaces du domaine — est la région où la viscosité joue un rôle important (couche limite). Alors  u et  v sont respectivement la vitesse sur la ligne de courant et la vitesse normale à la ligne de courant à l'intérieur de la couche limite. En utilisant l'analyse d'échelle, les équations de mouvement pour la couche limite se simplifient et deviennent :
 {\partial u\over\partial x}+{\partial v\over\partial y}=0
 u{\partial u \over \partial x}+v{\partial u \over \partial y}=-{1\over \rho} {\partial p \over \partial x}+{\nu}{\partial^2 u\over \partial y^2}
et si le fluide est incompressible, ce qui est le cas d'un liquide dans les conditions standards :
 {1\over \rho} {\partial p \over \partial y}=0
  Une analyse asymptotique montre que v, la vitesse normale, est petite comparée à u, la vitesse sur une ligne de courant, et que les propriétés de ses variations dans la direction de la ligne de courant sont généralement moins importantes que dans la direction normale.
La pression statique p est indépendante de y, alors la pression au bord de la couche limite est la pression de la ligne de courant. 
  La pression externe peut être calculée en appliquant le théorème de Bernoulli. Alors  u_0  u et  u_0 sont parallèles. En remplaçant p, les équations deviennent : est la vitesse du fluide en dehors de la couche limite, où
 u{\partial u \over \partial x}+v{\partial u \over \partial y}=u_0{\partial u_0 \over \partial x}+{\nu}{\partial^2 u\over \partial y^2}
avec les conditions limites
 {\partial u\over\partial x}+{\partial v\over\partial y}=0
   Pour un fluide dans lequel la pression statique  p ne dépend pas de la direction d'écoulement du fluide :
 {\partial p\over\partial x}=0
donc  u_0 reste constant.
Les équations de mouvement simplifiées sont :
 u{\partial u \over \partial x}+v{\partial u \over \partial y}={\nu}{\partial^2 u\over \partial y^2}
On utilise ces approximations dans un grand nombre de problèmes scientifiques et d'ingénierie. L'analyse précédente concerne toute couche limite (laminaire ou turbulente), mais les équations sont principalement utilisées pour étudier la couche limite laminaire.
    En effet, la moyenne de la vitesse correspond à la vitesse instantanée, car les fluctuations de vitesses sont absentes.

Couche limite turbulente

Le traitement de la couche limite turbulente est rendu difficile par la dépendance du fluide à la variable temps. 
   Une des techniques les plus courantes, quand le fluide est considéré comme turbulent, est d'appliquer la décomposition de Reynolds. 
   Dans ce cas, les propriétés instantanées du fluide sont décomposées entre la moyenne et les fluctuations à la moyenne. En appliquant cette technique, les équations de la couche limite donnent une couche limite pleinement turbulente ; elles sont peu abordées dans les publications scientifiques et techniques
 {\partial \overline{u}\over\partial x}+{\partial \overline{v}\over\partial y}=0
 \overline{u}{\partial \overline{u} \over \partial x}+\overline{v}{\partial \overline{u} \over \partial y}=-{1\over \rho} {\partial \overline{p} \over \partial x}+ \nu \left({\partial^2 \overline{u}\over \partial x^2}+{\partial^2 \overline{u}\over \partial y^2}\right)-\frac{\partial}{\partial y}(\overline{u'v'})-\frac{\partial}{\partial x}(\overline{u'^2})
 \overline{u}{\partial \overline{v} \over \partial x}+\overline{v}{\partial \overline{v} \over \partial y}=-{1\over \rho} {\partial \overline{p} \over \partial y}+\nu \left({\partial^2 \overline{v}\over \partial x^2}+{\partial^2 \overline{v}\over \partial y^2}\right)-\frac{\partial}{\partial x}(\overline{u'v'})-\frac{\partial}{\partial y}(\overline{v'^2})
En utilisant la même technique de calcul pour l'équation instantanée, les équations deviennent dans sa forme classique :
 {\partial \overline{u}\over\partial x}+{\partial \overline{v}\over\partial y}=0
 \overline{u}{\partial \overline{u} \over \partial x}+\overline{v}{\partial \overline{u} \over \partial y}=-{1\over \rho} {\partial \overline{p} \over \partial x}+{\nu}{\partial^2 \overline{u}\over \partial y^2}-\frac{\partial}{\partial y}(\overline{u'v'})
 {\partial \overline{p} \over \partial y}=0
   Le terme additionnel \overline{u'v'}dans la couche limite turbulente est connu sous le nom de « contrainte partagée de Reynolds ». La solution aux équations de la couche limite turbulente nécessite l'utilisation de modèles de turbulences, dont le but est d'exprimer la contrainte de Reynolds partagée dans des termes connus de variables du fluide et ses dérivés. 
   Le manque de précision et la non généralisation de ces modèles constituent l'obstacle majeur dans le succès des prédictions des propriétés des écoulements turbulents, dans la science de la dynamique des fluides.

Paramètres caractéristiques de la couche limites


Paramètres caractéristiques de la couche limites.
  1. Épaisseur de déplacement
    Les lignes de courants en écoulement visqueux sont déplacées par rapport à leurs position en fluide non-visqueux, ce déplacement est exploité pour définir une épaisseur tel que les aires A et A' soient égales.
     {A=\int_0^\infty (U_e-u)\, \mathrm dy}=A'=\delta_1 U_e\,\!
    d’où :  \delta_1=\int_0^\infty (1-\frac{u}{U_e})\, \mathrm dy.\!
  2. Épaisseur de la quantité de mouvement
    De la même manière, on définit l'épaisseur de la quantité de mouvement \delta_2\,\! :
     \delta_2=\int_0^\infty \frac{u}{U_e}\ (1-\frac{u}{U_e}) \,\mathrm dy.\!
  3. Épaisseur en énergie
    Similairement,on définit l'épaisseur en énergie \delta_3\,\! :
     \delta_3=\int_0^\infty \frac{u}{U_e}\ (1-\frac{u^2}{U_e^2}) \,\mathrm dy.\!

Domaines d'application

Anatomie

Article détaillé : Chair de poule.

Aéronautique

Aérodynamique


Générateurs de tourbillon installées au centre de la dérive d'un Panavia Tornado. Les six petites lames provoquent l'apparition de petites couches limites turbulentes.
La couche limite joue un rôle majeur dans les performances d'une surface portante : par exemple, le décollement de la couche limite sur une aile d'avion provoque une chute de la portance et une augmentation de la traînée de l'aile, ce qui correspond à une baisse notable des performances aérodynamiques de l'avion.
    Le décollement de la couche limite survient lorsque l'angle d'incidence de l'aile devient trop important, ce qui correspond pratiquement à une assiette cabrée de l'avion (à l'atterrissage par exemple). 
   Si cet angle est trop important, il se produit le phénomène de décrochage : la couche limite est fortement décollée et la portance peut chuter de façon très importante, plus ou moins brutalement. Ce phénomène est à l'origine de nombreux accidents aériens, la perte de portance pouvant entraîner la perte de contrôle de l'appareil.
Sur certains avions on trouve de petites lames, placées soit sur les ailes soit à l'arrière du fuselage, qui permettent de produire une couche limite turbulente qui résiste au décollement. Ces lames sont appelées « générateurs de tourbillons » (vortex generator en anglais).

Effets indésirables


Entrée d'air d'un SAAB 37 Viggen
La couche limite peut sérieusement perturber le fonctionnement d'un moteur à réaction, d'une part à cause des turbulences dans le flux d'air ingéré par le moteur, et d'autre part en réduisant son efficacité à cause de la faible vitesse de l'air au niveau de la couche. Ce problème ne se pose pas lorsque l'entrée d'air est frontale (dans le nez de l'avion) ou que le réacteur est contenu dans une nacelle fixée sous les ailes (cas de la grande majorité des avions civils).
Par contre, lorsque l'entrée d'air est située le long du fuselage (cas des avions militaires surtout), elle est le plus souvent légèrement écartée de celui-ci pour être placée hors de la couche limite. Une plaque métallique est parfois ajoutée juste avant l'entrée d'air pour maintenir la couche limite contre le fuselage : on parle alors de « piège à couche limite ».

Météorologie


Représentation de la couche limite atmosphérique.
En météorologie, on appelle couche limite planétaire la zone de l'atmosphère entre la surface (terre ou mer), où le frottement ralentit le déplacement de l'air, et l'atmosphère libre où cette dernière devient négligeable. Elle varie entre 0,5 et 3 km d'épaisseur selon la stabilité de l'air et la rugosité de la surface, la moyenne étant de 1 500 mètres[1]. L'étude théorique de cette tranche d'atmosphère divise en fait la couche limite planétaire comme la superposition de deux couches dont les épaisseurs sont très inégales[2] :
  • La couche d'Ekman (d'après le nom du physicien suédois Vagn Walfrid Ekman) dans laquelle le vent est causé par un équilibre entre le gradient de pression, la force de Coriolis, due à la rotation quotidienne de la Terre, et une portion de la friction diminuant graduellement jusqu'à l'atmosphère libre. La vitesse et la direction au sommet de cette couche est approximativement celle du vent géostrophique alors qu'elle diminue graduellement et tourne vers la plus basse pression à mesure qu'on descend vers le sol
  • La couche de surface ou couche limite de turbulence atmosphériqueimmédiatement au contact du sol et dont l'épaisseur ne dépasse pas le dixième de celle de l'ensemble de la couche limite. La vitesse de l'air y est causée par la convection due aux différences de températures et par les effets dynamiques du relief. Le flux y est turbulent.
On parle également d'une sous-couche rugueuse tout près de la surface, qui varie de quelques centimètres à quelques dizaines de mètres selon les aspérités du relief. La vitesse y tend vers zéro. Une conséquence pratique de cette couche est que la mesure des vents par un anémomètre se fait à 10 mètres du sol pour éviter autant que possible son effet de friction[2]. La mesure représente ainsi la vitesse moyenne du vent dans la couche de surface.
Les échanges de matière, d'énergie et de mouvement se produisant au sein de la couche limite planétaire sont primordiaux en météorologique. On y retrouve la plupart des éléments à méso-échelle qui mène au déclenchement de la convection profonde et une bonne partie des éléments qui mènent aux systèmes à l'échelle synoptique. La paramétrisation de la couche limite est donc primordiale dans la mise au point des modèles de prévision numérique du temps. Le rotor des éoliennes sont également placés le plus haut possible au-dessus du sol pour profiter du vent maximal, leur puissance étant environ proportionnelle au cube de la vitesse du vent[2].