Fully Turbulent Wakes at Low Reynolds Numbers: the Case of the Thin Flat Plate

Cet article démontre, par simulation numérique directe et comparaison expérimentale, que l'écoulement de sillage derrière une plaque plane bidimensionnelle mince devient pleinement turbulent à un nombre de Reynolds relativement bas de 400, présentant des caractéristiques statistiques et spectrales indiscernables de celles des sillages turbulents à nombres de Reynolds plus élevés, un chemin de transition qui diffère significativement de celui des cylindres circulaires ou carrés canoniques.

L'essentiel

Imaginez que vous tenez une fine plaque de carton (comme une carte à jouer) face à un vent fort. Alors que le vent frappe la carte, il crée derrière elle une traînée d'air désordonnée et agitée appelée « sillage ». Pendant longtemps, les scientifiques ont cru que pour que ce sillage devienne véritablement chaotique et « turbulent », il faudrait que le vent souffle très fort, ou que la carte ait une forme spécifique comme un tuyau rond ou un bloc carré.

Ce document raconte une histoire différente. Les chercheurs ont découvert que si vous utilisez une plaque plate et mince, l'air derrière elle devient pleinement chaotique et turbulent à une vitesse de vent bien plus faible que ce que l'on imaginait. En fait, cela se produit à une vitesse où, pour d'autres formes, l'air est encore relativement calme et ordonné.

Voici une décomposition de leurs découvertes utilisant des analogies simples :

1. La surprise du « Seuil de Turbulence »

Imaginez la turbulence comme une piste de danse bondée.

• L'ancienne croyance (Cylindres ronds) : Si vous avez un poteau rond dans le vent, l'air derrière lui commence par une danse calme et rythmée (un balancement d'avant en arrière). Il faut beaucoup d'énergie (vitesse élevée) avant que les danseurs ne commencent à se cogner, à tourner follement et à créer un désordre chaotique (turbulence). Cette transition se fait lentement sur une large gamme de vitesses.
• La nouvelle découverte (Plaque plate mince) : Les chercheurs ont découvert que pour une plaque plate mince, la « piste de danse » passe d'un état calme à un mosh pit sauvage presque instantanément. Même à une vitesse de vent relativement faible (nombre de Reynolds 400), l'air derrière la plaque est déjà pleinement chaotique. Elle ne passe pas par les étapes rythmiques et lentes que traversent les poteaux ronds ; elle saute directement à la fête.

2. Comment ils l'ont prouvé

Pour être sûrs qu'ils n'imaginaient pas les choses, l'équipe a agi comme des détectives comparant des scènes de crime.

• La simulation (Le laboratoire virtuel) : Ils ont utilisé des superordinateurs pour simuler l'impact du vent sur la plaque à basse vitesse (Re 150 et Re 400).
• Le test en conditions réelles (La soufflerie) : Ils ont également observé de réelles expériences où le vent soufflait beaucoup plus fort (Re 12 500 et Re 19 700).
• La correspondance : Lorsqu'ils ont comparé la simulation à basse vitesse (Re 400) avec les expériences réelles à haute vitesse, les motifs correspondaient parfaitement. Les « empreintes digitales » de la turbulence — la façon dont l'air se déplaçait, la quantité d'énergie qu'il possédait et la manière dont il tourbillonnait — étaient identiques.
• Le groupe témoin : Lorsqu'ils ont examiné la simulation à une vitesse encore plus basse (Re 150), les motifs étaient totalement différents. L'air était encore dans la phase « calme », pas encore chaotique. Cela a prouvé que la transition vers le chaos se produit quelque part entre 150 et 400.

3. L'« Empreinte digitale » de la turbulence

Comment savoir si un écoulement est véritablement turbulent ? Le document recherche des « signes de vie » spécifiques dans les données :

• Le spectre d'énergie (Le son du vent) : Dans un écoulement calme, l'énergie est concentrée sur quelques notes spécifiques (comme une flûte jouant une note unique). Dans un écoulement turbulent, cela ressemble à un bruit blanc ou à de la friture, avec de l'énergie répartie sur une immense gamme de fréquences. Les chercheurs ont constaté qu'à Re 400, le « son » du vent derrière la plaque était déjà rempli de cette friture chaotique, tout comme dans les expériences à haute vitesse.
• L'« Intermittence » (Le cri occasionnel) : Dans un écoulement véritablement turbulent, l'air ne se contente pas de tourbillonner doucement ; il présente des poussées soudaines et intenses de vitesse et de rotation. Les chercheurs ont trouvé ces « cris » dans les données à Re 400, mais ils étaient absents à Re 150.

4. Pourquoi est-ce différent ?

Le document suggère que la raison de ce saut soudain est la forme de l'objet.

• Objets ronds ou carrés : Lorsque le vent frappe un objet rond ou carré, la partie arrière de l'objet agit comme un bouclier, stabilisant l'écoulement de l'air derrière lui. Il faut beaucoup d'énergie pour briser cette stabilité.
• La plaque mince : Parce que la plaque est très mince, il n'y a pas d'« arrière » pour protéger l'air. Les fluctuations de pression (la poussée et la traction de l'air) sont directement liées aux vortex tourbillonnants dès le départ. C'est comme essayer de faire tenir un crayon en équilibre sur sa pointe plutôt que de faire tenir une boule de bowling ; le crayon (la plaque mince) est intrinsèquement instable et bascule dans le chaos beaucoup plus vite.

L'essentiel

Ce document modifie notre compréhension de la façon dont l'air circule autour des objets plats. Il prouve que les plaques plates minces créent des sillages pleinement turbulents à des vitesses étonnamment basses, bien plus basses que pour les objets ronds ou carrés. La transition vers le chaos n'est pas un processus lent et graduel pour ces formes ; c'est un changement soudain et fondamental qui se produit très tôt dans la plage de vitesse.

Les chercheurs n'ont pas discuté de la manière dont cela s'applique à la construction de ponts, à la conception de voitures ou de dispositifs médicaux. Ils se sont strictement concentrés sur la preuve que ce phénomène existe et sur la manière dont la physique de l'écoulement de l'air diffère de ce que nous pensions auparavant.

Résumé technique

Résumé technique : Sillage pleinement turbulent à bas nombres de Reynolds : le cas de la plaque plane mince

Énoncé du problème
La simulation numérique de l'écoulement nominalement bidimensionnel (2D) contournant une plaque plane mince perpendiculaire à un flux uniforme s'est historiquement révélée complexe, les études de mécanique des fluides numérique (CFD) produisant des résultats contradictoires concernant la dynamique du sillage à bas nombres de Reynolds ($Re < 1000$). Alors que les données expérimentales suggèrent que les paramètres intégraux (tels que les coefficients de traînée et les nombres de Strouhal) dépendent faiblement de $Re$ dans le régime 3D, les simulations ont rapporté des comportements de sillage très différents. Une question centrale demeure : le sillage d'une plaque plane mince passe-t-il à un état pleinement turbulent immédiatement après l'apparition des instabilités 3D, à l'instar des sillages de cylindres canoniques, ou suit-il un chemin de transition multi-étapes distinct ? Les auteurs examinent si le sillage derrière une plaque plane mince est déjà pleinement turbulent à un nombre de Reynolds relativement bas de $Re = 400$, un régime où les sillages de cylindres circulaires et carrés canoniques présentent généralement une forte dépendance à $Re$ et des caractéristiques de fluctuations non turbulentes.

Méthodologie
L'étude emploie une approche multidimensionnelle combinant des simulations numériques directes (DNS) et une validation expérimentale :

1. Simulations numériques : Deux codes DNS indépendants ont été utilisés pour garantir la robustesse :

   • Méthode d'élément Spectral/hp (SEM) : Implémentée dans Nektar++, utilisée pour la DNS principale de haute précision à $Re = 400$ (DNS400) et $Re = 150$ (DNS150).
   • Méthode des volumes finis (FVM) : Implémentée dans OpenFOAM, utilisée pour les études de convergence de maillage et la validation croisée (DNS400V).
   • Le domaine de calcul a été tronqué et soumis à des conditions de non-glissement sur la plaque, des conditions d'entrée de flux libre, et des conditions de sortie spécifiques proposées par Dong et al. La périodicité a été imposée dans la direction de l'empilement (spanwise).
   • La convergence du maillage a été rigoureusement testée, garantissant la résolution de l'échelle de longueur de Kolmogorov ($\eta$) avec des valeurs de $\Delta x/\eta$ comparables ou supérieures à la littérature existante pour des géométries similaires.

2. Données expérimentales : Deux ensembles de données expérimentales ont été acquis via la vélocimétrie par images de particules (PIV) stéréoscopique résolue en temps dans une soufflerie à l'Université de Calgary :

   • EX12K : $Re = 12\,500$
   • EX20K : $Re = 19\,700$
     Ces ensembles de données représentent des régimes pleinement turbulents et servent de référence pour les simulations à bas $Re$.

3. Analyse : L'étude compare les quantités intégrales globales, les champs de vitesse moyens, les contraintes de Reynolds, les budgets d'énergie cinétique turbulente (TKE), les spectres d'énergie et les propriétés statistiques des gradients de vitesse (spécifiquement la non-gaussianité) à travers les différents cas de $Re$.

Contributions clés et résultats

• Apparition précoce de la turbulence : Le résultat principal est que le sillage d'une plaque plane mince devient pleinement turbulent à $Re = 400$. C'est nettement inférieur à la transition observée pour les sillages de cylindres circulaires ($Re \approx 10^5$) ou carrés ($Re \approx 6000$) canoniques.
• Accord quantitatif avec les expériences à haut $Re$ : Les résultats DNS à $Re = 400$ (DNS400) montrent un accord remarquable avec les données expérimentales à $Re = 12\,500$ et $19\,700$ (EX12K, EX20K) concernant :
  • Quantités globales : Le coefficient de traînée moyen ($C_D$) et le nombre de Strouhal ($St$) sont des fonctions à variation lente de $Re$ pour $Re \ge 400$.
  • Structure de l'écoulement : Les lignes de courant moyennes, les profils de vitesse et les distributions de contraintes de Reynolds ($u'^2, v'^2, w'^2$) sont statistiquement indiscernables entre $Re=400$ et les expériences à haut $Re$ dans l'incertitude de mesure.
  • Budgets d'énergie : La distribution spatiale de la production de TKE et des termes d'advection dans l'équation de transport de $k$ est cohérente entre le cas $Re=400$ et les cas à haut $Re$.
• Signatures turbulentes à bas $Re$ : Le cas $Re=400$ présente des caractéristiques définitives de la turbulence absentes dans le cas $Re=150$ :
  • Spectres d'énergie : Les composantes de vitesse fluctuantes affichent des spectres à large bande avec une mise à l'échelle de la plage inertielle approximativement de type $f^{-5/3}$ (loi de Kolmogorov), indiquant une cascade d'énergie. En revanche, les spectres de $Re=150$ montrent des transferts d'énergie entre seulement quelques modes.
  • Non-gaussianité : Les fonctions de densité de probabilité (PDF) des taux de déformation ($\Sigma$) et de rotation ($\Omega$) présentent des queues lourdes, signifiant l'intermittence et les événements rares de haute intensité caractéristiques de la turbulence. Le cas $Re=150$ présente des distributions gaussiennes.
• Chemin de transition distinct : La transition à $Re=150$ implique des instabilités 3D mais manque de la structure cohérente et pleinement turbulente observée à $Re=400$. Le sillage de $Re=150$ présente un reflux plus faible, une récupération de sillage plus lente et une évolution non monotone de la demi-largeur du sillage, différant fondamentalement du comportement des cylindres canoniques où les instabilités 3D (Modes A, B, C) évoluent sur une large plage de $Re$ avant la turbulence complète.

Signification et affirmations
L'article soutient que le chemin vers la transition à la turbulence pour une plaque plane mince est fondamentalement différent de celui des cylindres canoniques. Alors que les sillages de cylindres subissent une séquence prolongée d'instabilités (Modes A, B, C) sur une large plage de nombres de Reynolds avant d'atteindre un état pleinement turbulent, le sillage de la plaque plane mince passe à un état pleinement turbulent presque immédiatement après l'apparition des instabilités 3D (entre $Re=150$ et $Re=400$).

Les auteurs suggèrent que cette différence peut être attribuée à l'absence d'un corps arrière longitudinal de la part de la géométrie. Pour les cylindres, le corps arrière stabilise l'écoulement et découple les fluctuations de pression du point de séparation. En revanche, pour les plaques planes minces, les fluctuations de pression sont directement couplées à la région de formation de vortex, ce qui pourrait favoriser une transition vers la turbulence plus précoce et plus rapide. Ces conclusions aident à réconcilier les divergences de la littérature précédente concernant la sensibilité des sillages de plaques minces aux rapports d'aspect et aux nombres de Reynolds, établissant que le régime « pleinement turbulent » pour cette géométrie commence à un $Re$ beaucoup plus bas que ce qui était précédemment supposé pour les corps épais canoniques.