Introduction
La couche limite est un concept fondamental en mécanique des fluides, désignant la zone d'interface entre un corps et le fluide environnant lors d'un mouvement relatif entre les deux, et ce, en raison de la viscosité du fluide. Lorsqu'un fluide réel s'écoule le long d'une paroi, les vitesses sont nulles sur la paroi et atteignent la vitesse de l'écoulement non perturbé à une certaine distance de l'obstacle. Dans cette zone, la vitesse varie donc entre zéro et un maximum. La couche limite peut être laminaire ou turbulente, selon les conditions d'écoulement. Cet article se concentre sur la couche limite turbulente, ses caractéristiques, sa modélisation et son importance dans divers domaines d'ingénierie. La compréhension et la modélisation des équations de la couche limite sont peut-être l'une des plus importantes avancées de la dynamique des fluides.
Définition et Formation de la Couche Limite
Lorsqu'un fluide se déplace autour d'un obstacle, les champs de vitesse et de pression sont perturbés. La couche limite est la zone de fluide contiguë à la paroi, dans laquelle se produit l'accommodation entre la vitesse nulle en paroi et la pleine vitesse locale du fluide libre. De nombreuses propriétés de l'écoulement (frottement, transferts thermiques et présence éventuelle de décollements) dépendent de cette couche limite.
La notion de couche limite est due à Ludwig Prandtl, en 1904. Elle a permis d'interpréter un bon nombre de phénomènes que la théorie des fluides parfaits ne suffisait pas à expliquer. C'est la couche limite qui conditionne directement la résistance de frottement du corps.
Transition Laminaire-Turbulent
L'écoulement dans la couche limite peut être laminaire ou turbulent. Initialement, près du bord d'attaque d'une surface, la couche limite est laminaire, caractérisée par un écoulement ordonné où les vecteurs de vitesse sont parallèles à un même plan. Cependant, lorsque le nombre de Reynolds dépasse une certaine valeur critique, une transition vers la turbulence se produit. Cette transition est influencée par divers facteurs :
- Nombre de Reynolds: Lorsque le nombre de Reynolds dépasse une certaine valeur critique, on observe une transition de la couche limite laminaire vers la turbulence.
- Courbure: L'effet de la courbure de paroi sur le point de transition est faible, et le nombre de Reynolds de transition est pratiquement inchangé par rapport au cas de la paroi plane.
- État de la surface (rugosité): L'effet des rugosités de paroi est d'avancer la transition, c'est-à-dire de la provoquer à un nombre de Reynolds plus faible. Toutefois, il existe une taille critique de rugosité en dessous de laquelle la position de la transition reste inchangée.
La couche limite turbulente, qui s’épaissit plus rapidement que la couche limite laminaire, résiste par ailleurs beaucoup mieux aux gradients de pression adverses.
Lire aussi: Couches bébé : guide d'achat
Caractéristiques de la Couche Limite Turbulente
La couche limite turbulente est caractérisée par un écoulement désordonné et chaotique, avec des fluctuations de vitesse importantes. Contrairement à la couche limite laminaire, où l'écoulement se fait en couches parallèles, la couche limite turbulente présente un mélange intense dû aux tourbillons de différentes tailles. Cette turbulence a plusieurs conséquences :
- Épaississement: La couche limite turbulente s'épaissit plus rapidement que la couche limite laminaire. L'épaisseur d'une couche limite varie suivant l'état de la surface (plus la surface est lisse, plus la couche est mince) mais aussi suivant la vitesse du fluide. L’accélération de l’écoulement hors de la couche limite contribue à amincir la couche limite.
- Résistance aux gradients de pression: La couche limite turbulente est plus résistante aux gradients de pression défavorables que la couche limite laminaire.
- Transfert de quantité de mouvement, de chaleur et de masse: La turbulence favorise le mélange et augmente les taux de transfert de quantité de mouvement, de chaleur et de masse.
Modélisation de la Couche Limite Turbulente
Le traitement de la couche limite turbulente est plus complexe en raison de la dépendance du fluide à la variable temps. Une des techniques les plus courantes, lorsque 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.
Décomposition de Reynolds et Contraintes de Reynolds
Dans la décomposition de Reynolds, les propriétés instantanées du fluide (comme la vitesse et la pression) sont décomposées en une composante moyenne et une composante fluctuante. Par exemple, la vitesse instantanée u est exprimée comme u = <u> + u', où <u> est la vitesse moyenne et u' est la fluctuation de vitesse.
En appliquant cette décomposition aux équations de Navier-Stokes, on obtient les équations de Reynolds-moyennées (RANS). Ces équations contiennent des termes supplémentaires, appelés contraintes de Reynolds, qui représentent l'effet des fluctuations turbulentes sur l'écoulement moyen. La contrainte de Reynolds partagé dans la couche limite turbulente est connu sous le nom de contrainte partagé de Reynolds.
Modèles de Turbulence
La solution aux équations de la couche limite turbulente nécessite l'utilisation de modèles de turbulence, dont le but est d'exprimer la contrainte de Reynolds partagé dans des termes connu de variable du fluide et ses dérivés. Puisque les contraintes de Reynolds introduisent plus d'inconnues que d'équations, il est nécessaire de recourir à des modèles de turbulence pour fermer le système d'équations. Ces modèles visent à exprimer les contraintes de Reynolds en fonction des variables moyennes de l'écoulement et de leurs dérivées. Il existe de nombreux modèles de turbulence, allant des modèles algébriques simples aux modèles plus complexes basés sur des équations de transport pour l'énergie cinétique turbulente et d'autres quantités.
Lire aussi: Causes et symptômes de la fausse couche
- Modèles algébriques (de longueur de mélange): Ces modèles sont les plus simples et utilisent une longueur de mélange pour estimer les contraintes de Reynolds.
- Modèles à une équation (par exemple, Spalart-Allmaras): Ces modèles résolvent une équation de transport pour une variable liée à la turbulence, comme la viscosité turbulente.
- Modèles à deux équations (par exemple, k-ε, k-ω): Ces modèles résolvent deux équations de transport pour l'énergie cinétique turbulente (k) et une autre variable, comme le taux de dissipation de l'énergie cinétique turbulente (ε) ou la fréquence spécifique de dissipation (ω).
- Modèles de contraintes de Reynolds (RSM): Ces modèles résolvent directement des équations de transport pour les contraintes de Reynolds, offrant une modélisation plus précise des écoulements complexes.
- Large Eddy Simulation (LES): Cette approche résout les équations de Navier-Stokes filtrées, en résolvant directement les grandes échelles de turbulence et en modélisant les petites échelles.
- Direct Numerical Simulation (DNS): Cette approche résout directement les équations de Navier-Stokes sans aucune modélisation de la turbulence, mais elle est extrêmement coûteuse en termes de calcul et ne peut être appliquée qu'à des écoulements simples à faible nombre de Reynolds.
Analyse d'échelle et Simplification des Équations de Navier-Stokes
En utilisant 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.
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 ou 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. 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 u0 est la vitesse du fluide en dehors de la couche limite, où u et u0 sont parallèles. Ces approximations sont utilisées dans un grand nombre de problèmes scientifiques et d'ingénierie. L'analyse précédente est pour toute couche limite (laminaire ou turbulente), mais les équations sont principalement utilisées pour étudier la couche limite laminaire.
Importance et Applications
La compréhension et le contrôle de la couche limite turbulente sont cruciaux dans de nombreux domaines de l'ingénierie, notamment :
- Aéronautique: La couche limite sur les ailes d'avion influence la portance, la traînée et la stabilité. Le contrôle de la couche limite permet de réduire la traînée, d'augmenter la portance et d'améliorer l'efficacité énergétique.
- Automobile: La couche limite sur la carrosserie des voitures affecte la traînée aérodynamique et la consommation de carburant.
- Énergie éolienne: La couche limite sur les pales d'éolienne influence la production d'énergie.
- Transfert de chaleur: La couche limite thermique joue un rôle important dans les processus de transfert de chaleur, par exemple dans les échangeurs de chaleur et les systèmes de refroidissement.
- Génie chimique: La notion de couche limite est très importante en génie des procédés, car c'est dans cette couche que se trouvent généralement les plus importantes résistances aux transferts de matière et de chaleur.
Contrôle de la Couche Limite
Le décollement de la couche limite, générateur de perte de portance et d’augmentation de traînée, nuit considérablement aux performances des profils d’ailes, il faut donc chercher à le réduire. Pour contrôler le décollement, il est apparu naturel d’intervenir au niveau de l’écoulement de proche paroi, car c’est en agissant au niveau de la couche limite qu’il est possible d’influer sur les caractéristiques aérodynamiques globales d’un profil. volets à fentes au bord de fuite des ailes. Ces hypersustentateurs sont entièrement déployés à l’atterrissage. Pour des raisons de sécurité, l’atterrissage doit s’effectuer à la vitesse la plus faible possible. Pour maintenir une portance suffisante à faible vitesse, il est nécessaire d’augmenter considérablement l’angle d’incidence. Les becs à fente de bord d’attaque ré-injectent de l’air à vitesse élevée sur l’extrados et retardent ainsi le décollement de la couche limite.
Avec le développement de nouvelles technologies et l’amélioration des connaissances relatives à la dynamique interne de la couche limite turbulente, plusieurs techniques ont été explorées telles que la modification de géométrie pour obtenir un gradient de pression favorable plus long, l’utilisation de l’aspiration pour stabiliser la couche limite, ou encore le refroidissement de la surface. Le contrôle de la couche limite ou le contrôle de l'écoulement laminaire permet de réduire considérablement la résistance au frottement, et d'améliorer la qualité aérodynamique de l'écoulement aux grands angles d'attaque.
Lire aussi: Couche-Tard : Stratégies et développement
Couche limite sur une plaque plane
Pour une plaque plane de longueur (L), on dispose des corrélations qui suivent.
Couche limite laminaire
L'épaisseur de la couche limite est donnée par la relation :
(\frac{\delta }{x}=4,64\cdot R{{e}_{x}}^{-1/2})
Le coefficient de frottement (qui n'est autre que le rapport de la contrainte à la paroi par le terme d'énergie cinétique (\frac{\rho \cdot {{u}{\infty }}^{2}}{2})), est le suivant (avec (R{{e}{x}}=\frac{\rho \cdot u\cdot x}{\mu }) ; (Re=\frac{\rho \cdot u\cdot L}{\mu })) :
localement : ({{{}^{f}/{}{2}}{x}}=0,353\cdot R{{e}{x}}^{-1/2}) ; globalement : ({}^{f}/{}{2}=0,646\cdot R{{e}^{-1/2}})
Ceci permet de calculer la contrainte au niveau de la paroi.
Couche limite modérément turbulente
Pour (10^5 < Re < 10^7), l'épaisseur de la couche limite est donnée par la relation :
(\frac{\delta }{x}=0,376\cdot R{{e}_{x}}^{-0,2})
L'épaisseur de la sous-couche limite laminaire est la suivante :
(\frac{\delta }{x}=72,4\cdot R{{e}_{x}}^{-0,9})
Les coefficients de frottement local et global sont les suivants :
({{{}^{f}/{}{2}}{x}}=0,0293\cdot R{{e}{x}}^{-0,2}) ; ({}^{f}/{}{2}=0,0366\cdot R{{e}^{-0,2}})
Couche limite turbulente
Pour (10^7 < Re ), l'épaisseur de la couche limite et le coefficient de frottement global sont donnés par les expressions de Schoenhen :
(\left{ \begin{array}{r l}& \frac{\delta }{x}=0,54\cdot \sqrt{{}^{f}/{}{2}} \& \frac{0,171}{\sqrt{{}^{f}/{}{2}}}=\log \ \left( 2\cdot Re\cdot {}^{f}/{}_{2} \right) \\end{array} \right.)
Couche limite fortement turbulente
Pour (10^7 < Re < 10^9,) on donne l'approximation suivante, sous forme explicite (Van Karman & Schliching) :
(\left{ \begin{array} {r l}& \frac{\delta }{x}=0,22\cdot R{{e}{x}}^{1/6} \& \& {}^{f}/{}{2}=0,456\cdot {{\left( \log \ Re \right)}^{-2,58}} \\end{array} \right.)
Couche limite dans une conduite
À l'intérieur d'une conduite, on observe également le développement d'une couche limite. Toutefois, dans cette géométrie confinée, au bout d'une distance suffisamment longue par rapport à l'entrée dans la conduite, la couche limite envahit toute la conduite.
Ainsi, si la couche limite envahit toute la conduite alors que la couche limite n'était encore que laminaire, l'écoulement dans la conduite sera laminaire (figure de gauche), on comprend que ceci se produira en particulier si la conduite a un faible diamètre. Si au contraire, la couche limite turbulente a pu se développer avant que la couche limite n'envahisse toute la conduite, l'écoulement sera turbulent (figure de droite).
L'épaisseur de la sous-couche limite laminaire dans une conduite circulaire est globalement proportionnelle à l'inverse de la vitesse moyenne, et est donnée par l'expression :
(\frac{\delta }{D}=60\cdot R{{e}^{-7/8}})
La longueur d'établissement ({L}_{e}) du régime d'écoulement correspond à la distance entre l'entrée dans la conduite et la position à partir de laquelle la couche limite a envahi toute la conduite.
tags: #couche #limite #turbulente #définition