On Capturing Laminar/Turbulent Regions Over a Wing Using WMLES

Cette étude démontre que pour simuler avec précision à la fois les régions laminaires et turbulentes sur une aile en utilisant la simulation des grandes échelles modélisée par paroi (WMLES), il est nécessaire de combiner un maillage adapté aux variations de l'épaisseur de la couche limite avec l'introduction de perturbations instationnaires en amont du point neutre pour déclencher correctement la transition.

L'essentiel

Imaginez que vous essayez de prédire comment l'air glisse sur l'aile d'un avion. C'est un peu comme essayer de filmer une rivière : parfois l'eau coule doucement et uniformément (c'est le régime laminaire), et parfois elle devient turbulente, avec des tourbillons et des remous chaotiques (c'est le régime turbulent).

Le défi pour les ingénieurs est de créer une simulation informatique capable de voir les deux états en même temps, car sur une vraie aile, l'air commence souvent par être calme avant de devenir turbulent.

Voici l'histoire de cette recherche, racontée simplement :

1. Le Problème : La Caméra trop petite ou trop grosse

Les chercheurs ont utilisé un outil puissant appelé WMLES. Pour faire simple, imaginez que c'est une caméra qui filme l'écoulement de l'air.

• Pour voir les détails fins de l'eau calme (le régime laminaire), il faut une caméra très haute définition, avec des pixels minuscules collés à la surface.
• Pour voir les gros tourbillons de l'eau agitée (le régime turbulent), la caméra peut être un peu moins précise, car les détails fins sont moins importants là-bas.

Le problème, c'est que sur une aile d'avion, l'air passe du calme au turbulent. Si vous utilisez une caméra avec des pixels fixes (une grille de calcul constante) :

• Si les pixels sont trop gros, vous voyez bien les gros tourbillons, mais vous manquez totalement les détails fins de l'eau calme au début.
• Si les pixels sont trop petits, vous voyez parfaitement l'eau calme, mais vous vous perdez dans les détails inutiles pour les tourbillons, et la simulation devient si lente qu'elle ne fonctionne plus bien.

C'est comme essayer de prendre une photo d'une foule : si vous zoomez trop sur un visage, vous ne voyez plus la foule entière. Si vous zoomez trop loin, vous ne voyez plus les visages.

2. L'Expérience : Essayer différentes "Caméras"

Les chercheurs ont testé plusieurs configurations sur une aile de type NACA0012 (une forme standard d'aile d'avion).

• Le premier essai (Tout turbulent) : Ils ont supposé que l'air était turbulent partout, dès le début. Avec une grille "standard" (pixels de taille moyenne), ça a plutôt bien fonctionné pour la partie turbulente, mais c'était un peu flou au tout début.
• Le deuxième essai (Laminaire puis turbulent) : Ils ont essayé de simuler la transition naturelle.
  • Quand ils ont mis des pixels très fins pour voir le début calme, la simulation a échoué à créer les tourbillons plus loin. C'est comme si la caméra était si proche du sol qu'elle ne voyait plus les nuages au-dessus.
  • Quand ils ont gardé des pixels plus gros, ils voyaient les tourbillons, mais l'air restait "calme" artificiellement trop longtemps.

3. La Solution : Une Caméra Intelligente et un "Coup de Pouce"

Pour résoudre ce casse-tête, les chercheurs ont eu deux idées de génie :

A. La grille qui s'adapte (La caméra intelligente)
Au lieu d'avoir des pixels de taille fixe, ils ont créé une grille qui change de taille selon l'épaisseur de la couche d'air.

• Là où l'air est fin (début de l'aile), les pixels sont minuscules.
• Là où l'air s'épaissit (vers la fin de l'aile), les pixels s'agrandissent.
  C'est comme un zoom automatique qui reste net sur le visage (le début calme) et qui s'adapte pour voir le paysage (la fin turbulente).

B. Le "Coup de Pouce" (Le déclencheur)
Même avec la bonne grille, la simulation avait du mal à passer du calme au turbulent naturellement. C'est comme si l'eau dans la rivière hésitait à se mettre à bouger.
Les chercheurs ont donc ajouté un petit "coup de pouce" numérique au début de l'aile : ils ont injecté de petites vibrations (des ondes) pour forcer l'air à devenir turbulent au bon moment. C'est comme donner un petit coup de pied à une balle pour qu'elle commence à rouler.

4. Le Résultat Final

En combinant une grille intelligente (qui s'adapte à la forme de l'écoulement) et un petit coup de pouce pour déclencher la turbulence, les chercheurs ont réussi à obtenir une simulation parfaite.

• Ils voient clairement l'air calme au début.
• Ils voient la transition précise.
• Ils voient les tourbillons à la fin.

En résumé

Ce papier nous dit que pour simuler l'air sur une aile d'avion avec précision, on ne peut pas utiliser la même "règle" partout. Il faut être flexible : utiliser des détails fins là où c'est nécessaire, et des détails plus larges là où ça va vite. Et parfois, il faut juste donner un petit coup de pouce à la simulation pour qu'elle se comporte comme la réalité. C'est une avancée majeure pour concevoir des avions plus efficaces et plus silencieux.

Résumé technique

Titre de l'étude

Capture des régions laminaire et turbulente au-dessus d'une aile par la simulation des grandes échelles avec modélisation de paroi (WMLES)

1. Problématique

L'objectif principal de cette étude est d'évaluer la capacité de la Simulation des Grandes Échelles avec modélisation de paroi (WMLES) à prédire avec précision les écoulements autour d'une aile présentant à la fois des zones laminaire et turbulente.

• Contexte : Les méthodes de résolution directe (DNS) et les LES résolues aux parois (WRLES) sont trop coûteuses en calcul pour les géométries complexes à haut nombre de Reynolds. La WMLES offre un compromis en modélisant la couche interne de la couche limite et en résolvant la partie externe.
• Défi spécifique : La difficulté réside dans la génération d'une grille adaptée. La couche limite laminaire est très mince, tandis que la couche limite turbulente est beaucoup plus épaisse (augmentant d'un ordre de grandeur lors de la transition).
  • Une grille grossière (adaptée à la zone turbulente) échoue à résoudre la région laminaire.
  • Une grille fine (adaptée à la région laminaire) place le premier point hors paroi dans la sous-couche tampon, ce qui empêche le modèle de paroi de fonctionner correctement dans la zone turbulente, retardant ainsi la transition.
• Objectif : Documenter les exigences de résolution de la grille pour capturer les deux régimes et proposer des pratiques améliorées pour la génération de grilles et la modélisation dans le cadre WMLES.

2. Méthodologie

Les simulations ont été réalisées sur un profil NACA0012 à 0° d'incidence et un nombre de Reynolds basé sur la corde de $3 \times 10^6$ (Mach 0.3).

• Code et Modèles : Utilisation du code FUN3D (NASA LaRC) basé sur un volume fini non structuré. Le modèle de paroi utilisé est une loi de paroi d'équilibre (loi de Spalding) avec le premier point hors paroi comme lieu d'échange.
• Scénarios simulés :
  1. Cas entièrement turbulent : Toute la surface est supposée turbulente.
  2. Cas laminaire-turbulent : La transition est prédite via des calculs de stabilité linéaire (méthode N-Factor, seuil N=9), situant le début de la transition à $x \approx 0.36$.
• Stratégies de grille et de forçage :
  • Grille structurée C-type : Grille uniforme en direction normale à la paroi (premier point à $\Delta z^+ \approx 3 \times 10^{-4}$).
  • Grille non structurée adaptative : Grille générée avec un nombre fixe de points par épaisseur de couche limite (estimée via une simulation RANS précurseur), permettant d'ajuster la densité de points selon l'épaisseur locale.
  • Forçage des perturbations : Introduction de perturbations harmoniques (fréquence la plus amplifiée, $F_0 = 6.5 \times 10^{-5}$) en amont du point neutre pour déclencher la transition de manière contrôlée.
• Comparaisons : Les résultats WMLES sont comparés à des simulations RANS (modèle SA), des WRLES (résolution complète) et des données expérimentales.

3. Contributions Clés

L'article apporte plusieurs contributions techniques majeures :

1. Analyse des exigences de grille : Démonstration qu'une grille à espacement constant ne peut pas simultanément résoudre correctement la région laminaire (nécessite des points fins) et la région turbulente (nécessite que le premier point soit dans la couche logarithmique).
2. Stratégie de grille adaptative : Validation de l'utilisation d'une grille non structurée basée sur l'épaisseur de la couche limite locale (environ 20 points par épaisseur de couche limite) pour capturer la physique laminaire.
3. Nécessité du forçage explicite : Mise en évidence que, sans forçage externe (trip), la WMLES sur des grilles adaptées à la couche laminaire ne déclenche pas la transition à l'endroit prévu, retardant la turbulence. L'introduction de perturbations à la fréquence la plus amplifiée est cruciale pour obtenir une transition réaliste.
4. Validation hybride : Démonstration qu'il est possible d'obtenir des résultats précis pour les deux régimes en combinant une grille adaptative et un forçage de transition, là où les approches précédentes échouaient souvent sur l'un des deux régimes.

4. Résultats

• Cas entièrement turbulent : La WMLES capture bien la friction de peau et les profils de vitesse dans la zone turbulente, en accord avec les résultats RANS et WRLES, bien que la friction soit légèrement sous-estimée près du bord d'attaque en raison de la grille grossière.
• Cas laminaire-turbulent (Grille uniforme) :
  • Avec une grille grossière ($\Delta z = 3 \times 10^{-4}$) : La région turbulente est bien capturée, mais la région laminaire est mal résolue (peu de points dans la couche limite).
  • Avec une grille fine ($\Delta z = 3 \times 10^{-5}$) : La région laminaire est bien résolue, mais la transition est retardée ($x=0.6$ au lieu de $0.4$) car le premier point se trouve dans la sous-couche tampon, empêchant le modèle de paroi de fonctionner.
• Cas laminaire-turbulent (Grille adaptative + Forçage) :
  • L'utilisation d'une grille avec un nombre fixe de points par épaisseur de couche limite (20 points) permet de résoudre correctement la région laminaire.
  • L'ajout de perturbations (forçage) à la fréquence la plus amplifiée permet de déclencher la transition à l'endroit prédit par la stabilité linéaire ($x \approx 0.4$).
  • Les coefficients de pression et de friction de peau moyens correspondent très bien aux résultats RANS et aux attentes physiques dans les deux régimes.
  • Les profils de vitesse et les contraintes de Reynolds sont bien capturés.

5. Signification et Conclusion

Cette étude souligne que la prédiction précise des écoulements avec transition laminaire-turbulente en WMLES est un défi complexe lié à la génération de grilles.

• Conclusion principale : Il n'existe pas de grille universelle à espacement constant pour ces écoulements. La stratégie la plus prometteuse consiste à :
  1. Générer une grille dont la distribution normale à la paroi varie selon l'épaisseur estimée de la couche limite (via RANS précurseur).
  2. Inclure un mécanisme de déclenchement explicite (forçage de perturbations) pour assurer que la transition se produit à l'endroit correct, car les grilles WMLES sont trop grossières pour capturer la croissance naturelle des instabilités.
• Impact : Ces résultats ouvrent la voie à l'application de la WMLES sur des configurations aéronautiques plus complexes (3D, haute portance) où la prédiction précise de la transition est critique pour le calcul de la traînée et de la portance, tout en maintenant un coût de calcul raisonnable par rapport aux méthodes DNS/WRLES.