La couche limite turbulente est un concept fondamental en mécanique des fluides, particulièrement pertinent dans des domaines tels que le génie civil, l'aérodynamique et les transferts thermiques. Elle décrit une région mince de fluide adjacente à une surface solide où les effets de la viscosité sont significatifs et où l'écoulement devient turbulent. Comprendre cette couche est crucial pour la conception et l'analyse de nombreux systèmes d'ingénierie.
Profils de Vitesse et Couche Limite
Lorsqu'un fluide s'écoule le long d'une surface solide, un gradient de vitesse se forme. À la paroi, la vitesse du fluide est nulle, une condition connue sous le nom d'adhérence. La vitesse augmente progressivement à mesure que l'on s'éloigne de la paroi, jusqu'à atteindre une vitesse maximale au centre de la conduite. Cette variation de vitesse est due à la viscosité du fluide.
La notion de couche limite est essentielle car c'est dans cette zone que se concentrent les plus importantes résistances aux transferts de matière et de chaleur. Considérons un fluide s'écoulant le long d'une surface plane. La vitesse du fluide est nulle à la paroi, puis croît progressivement jusqu'à atteindre la vitesse (u{\infty}) loin de la paroi. Cet effet est confiné dans une épaisseur appelée couche limite, notée (\delta), où (u{\delta} = 0,99 \cdot u_{\infty}).
Transition Laminaire-Turbulent
Dans le cas d'une plaque plane, l'écoulement dans la couche limite est initialement laminaire près du bord d'attaque. Cependant, à une certaine distance (xc), l'écoulement transite vers un régime turbulent. Cette transition est caractérisée par le nombre de Reynolds critique (Rec), tel que (3 \times 10^5 < Rec = \frac{\rho \cdot u{\infty} \cdot x_c}{\mu} < 3 \times 10^6), où (\rho) est la masse volumique du fluide et (\mu) sa viscosité dynamique.
Entre la paroi et la couche limite turbulente, une sous-couche limite laminaire persiste, bien que d'épaisseur réduite. Cette théorie est cruciale pour l'étude des pompes centrifuges, des coques de bateau et des fuselages. En génie des procédés, la couche limite joue un rôle essentiel dans les transferts de matière et de chaleur.
Lire aussi: Laminaire vs. Turbulente
Couche Limite sur une Plaque Plane
Pour une plaque plane de longueur (L), différentes corrélations permettent d'évaluer l'épaisseur de la couche limite et le coefficient de frottement.
Couche Limite Laminaire
L'épaisseur de la couche limite laminaire est donnée par :[\frac{\delta}{x} = 4,64 \cdot Rex^{-1/2}]où (Rex = \frac{\rho \cdot u \cdot x}{\mu}).
Le coefficient de frottement local ({{{}^{f}/{}{2}}{x}}) et global ({}^{f}/{}{2}) sont respectivement :[{{{}^{f}/{}{2}}{x}} = 0,353 \cdot Rex^{-1/2}][{}^{f}/{}_{2} = 0,646 \cdot Re^{-1/2}]
Couche Limite Modérément Turbulente
Pour (10^5 < Re < 10^7), l'épaisseur de la couche limite est :[\frac{\delta}{x} = 0,376 \cdot Re_x^{-0,2}]
L'épaisseur de la sous-couche limite laminaire est :[\frac{\delta}{x} = 72,4 \cdot Re_x^{-0,9}]
Lire aussi: Choisir les meilleures couches de piscine pour bébés
Les coefficients de frottement local et global sont :[{{{}^{f}/{}{2}}{x}} = 0,0293 \cdot Rex^{-0,2}][{}^{f}/{}{2} = 0,0366 \cdot Re^{-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 :[\begin{cases}\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{cases}]
Couche Limite Fortement Turbulente
Pour (10^7 < Re < 10^9), l'approximation suivante (Van Karman & Schliching) est utilisée :[\begin{cases}\frac{\delta}{x} = 0,22 \cdot Rex^{1/6} \{}^{f}/{}{2} = 0,456 \cdot \left( \log \ Re \right)^{-2,58}\end{cases}]
Couche Limite dans une Conduite
À l'intérieur d'une conduite, une couche limite se développe également. Après une certaine distance, cette couche limite envahit toute la conduite. Si la couche limite est laminaire lorsqu'elle envahit la conduite, l'écoulement reste laminaire. Si elle est turbulente, l'écoulement devient turbulent.
L'épaisseur de la sous-couche limite laminaire dans une conduite circulaire est proportionnelle à l'inverse de la vitesse moyenne :[\frac{\delta}{D} = 60 \cdot Re^{-7/8}]où (D) est le diamètre de la conduite.
Lire aussi: Comprendre les pleurs nocturnes de bébé
La longueur d'établissement (L_e) est la distance nécessaire pour que la couche limite envahisse toute la conduite.
Modélisation Numérique de la Turbulence
Le logiciel COMSOL Multiphysics® propose différentes formulations pour résoudre les problèmes d'écoulement turbulent, notamment les modèles de turbulence L-VEL, yPlus algébrique, Spalart-Allmaras, k-ε, k-ω, k-ε bas Reynolds, SST et v2-f.
Modèles de Turbulence RANS
Les modèles RANS (Reynolds-averaged Navier-Stokes) sont basés sur l'observation que l'écoulement (u) contient de petites oscillations locales (u') et peut être traité comme une moyenne temporelle (U).
Modèles Algébriques
Les modèles algébriques utilisent des équations algébriques dépendant du champ de vitesse et de la distance aux parois pour décrire l'intensité de la turbulence. Une viscosité turbulente est calculée et ajoutée à la viscosité moléculaire du fluide.
Modèles Bas Reynolds
Ces modèles reproduisent le comportement de l'écoulement lorsque la distance à la paroi approche zéro. Ils donnent une description précise de la couche limite. La plupart des modèles utilisant ω sont des modèles bas Reynolds.
Régions de l'Écoulement Turbulent Proche d'une Paroi
L'écoulement turbulent proche d'une paroi plane peut être divisé en quatre régions :
- Sous-couche visqueuse (ou laminaire) : La vitesse du fluide évolue linéairement avec la distance à la paroi.
- Couche tampon : Les contraintes turbulentes commencent à dominer les contraintes visqueuses.
- Région de loi logarithmique : La vitesse moyenne évolue avec le logarithme de la distance à la paroi.
- Région d'écoulement libre : L'écoulement est complètement turbulent.
Lois de Paroi
Les lois de paroi ne résolvent pas l'écoulement dans la couche tampon, mais calculent analytiquement une vitesse non nulle à la paroi. Elles nécessitent moins de ressources de calcul.
Traitement Automatique de la Paroi
Cette fonctionnalité combine les avantages des lois de parois et des modèles bas Reynolds. Le logiciel mélange les deux formulations dans les éléments de frontière.
Différents Modèles de Turbulence RANS
Les huit modèles de turbulence RANS diffèrent par leur modélisation de l'écoulement en proche paroi, le nombre de variables additionnelles qu'ils résolvent et ce que ces variables représentent.
- Modèles L-VEL et yPlus algébrique : Calculent la viscosité turbulente à l'aide d'équations algébriques basées sur la vitesse locale et la distance à la paroi. Ils ne résolvent aucune équation de transport supplémentaire.
- Modèle Spalart-Allmaras : Ajoute une seule variable supplémentaire pour la viscosité cinématique turbulente non amortie. C'est un modèle bas Reynolds, robuste et avec un coût de résolution modéré.
- Modèle k-ε : Résout deux variables : k (énergie cinétique turbulente) et ε (taux de dissipation de l'énergie cinétique turbulente). Des lois de paroi sont utilisées. Il est populaire pour les applications industrielles en raison de son bon taux de convergence et de ses faibles besoins en mémoire.
- Modèle k-ω : Similaire au modèle k-ε, mais résout ω (taux de dissipation spécifique de l'énergie turbulente). C'est un modèle bas Reynolds, utile pour les écoulements internes, les écoulements avec de fortes courbures, les écoulements séparés et les jets.
- Modèle k-ε bas Reynolds : Similaire au modèle k-ε, mais ne nécessite pas de lois de paroi.
- Modèle SST : Combine le modèle k-ε dans la région d'écoulement libre et le modèle k-ω en proche paroi. C'est un modèle bas Reynolds, souvent privilégié pour les applications industrielles.
- Modèle v2-f : Décrit l'anisotropie de l'intensité de la turbulence dans la couche limite turbulente à l'aide de deux nouvelles équations, s'ajoutant aux deux équations de l'énergie cinétique turbulente (k) et du taux de dissipation (ε).
Vérification de la Précision des Simulations
Après avoir utilisé un modèle de turbulence, il est important de vérifier la précision de la solution.
- Pour les formulations avec loi de paroi, vérifiez la résolution à la paroi en unités visqueuses et le décalage de paroi.
- Pour les simulations utilisant un traitement bas Reynolds, vérifiez la distance adimensionnée au centre de la cellule.
Influence de la Rugosité de la Paroi
La rugosité de la paroi (ou l'état de surface) joue un rôle important dans la couche limite. À partir d'une certaine vitesse, et en fonction du profil de la paroi et de sa rugosité, la couche limite s'épaissit en passant en régime turbulent. L'état de surface influence l'épaisseur de la couche limite.
tags: #couche #limite #turbulente #définition