Numerical simulations of transition and long-term response of a wind turbine airfoil

Cet article présente des simulations numériques utilisant Nek5000 et EllipSys pour analyser la transition et la réponse d'écoulement à long terme d'un profil FFA-W3 à bas nombres de Reynolds, validant la capacité du code EllipSys à capturer les instabilités de Kelvin-Helmholtz et identifiant une modulation lente et périodique du coefficient de force normale liée à des oscillations de basse fréquence et à un potentiel éclatement de bulles.

L'essentiel

Imaginez une pale d'éolienne comme une immense aile tournante. Tout comme l'aile d'un avion, elle a besoin d'un flux d'air fluide pour fonctionner efficacement. Mais lorsque le vent frappe la pale selon certains angles, l'air peut s'« accrocher » et se détacher de la surface, créant un désordre chaotique et tourbillonnant. Ce document est comme une expérience de soufflerie de haute technologie, mais au lieu d'utiliser un modèle physique, les chercheurs en ont construit un virtuel à l'intérieur d'un superordinateur pour observer exactement comment cet air se comporte.

Voici l'histoire de leur découverte, décomposée en concepts simples :

1. La soufflerie virtuelle

Les chercheurs voulaient étudier une tranche spécifique d'une pale d'éolienne massive (d'une turbine de 10 MW). Ils ont utilisé deux programmes informatiques différents, NEK5000 et ELLIPSYS, pour simuler l'écoulement de l'air sur cette pale.

Considérez NEK5000 comme une caméra haut de gamme, ultra-précise, qui capture chaque petit détail mais qui est très lente et coûteuse à faire fonctionner. ELLIPSYS est comme une caméra légèrement plus rapide et plus efficace. L'équipe a d'abord dû prouver que la « caméra plus rapide » (ELLIPSYS) pouvait voir les mêmes choses que la « haute gamme ». Ils ont découvert que, bien qu'ELLIPSYS manque quelques ondulations infimes et ténues dans l'air fluide (car elle lisse un peu trop les choses), elle est excellente pour capturer les grands tourbillons chaotiques qui comptent réellement pour la performance de la pale.

2. Quelle largeur doit avoir le tunnel ?

Avant de lancer les simulations de longue durée, ils ont dû déterminer la largeur nécessaire de leur « soufflerie virtuelle ». Si le tunnel est trop étroit, il pourrait comprimer l'air et créer des résultats faussés. S'il est trop large, il gaspille la puissance de calcul.

Ils ont testé un tunnel « étroit » (10 % de la largeur de l'aile) contre un tunnel « large » (20 % de la largeur).

• L'analogie : Imaginez que vous observez la course d'une rivière. Si vous ne regardez qu'une étroite bande de la rivière, manquez-vous les grandes vagues ?
• Le résultat : Ils ont trouvé que le tunnel étroit était en fait suffisant. Les grandes vagues et les tourbillons se formaient parfaitement dans l'espace étroit. Cela signifiait qu'ils pouvaient économiser beaucoup de temps de calcul en utilisant la boîte de simulation plus petite et plus étroite.

3. La « bulle » et le « battement »

La partie la plus intéressante de l'étude s'est produite sur le dessus de l'aile (le côté aspiration).

• La bulle de séparation : À mesure que l'air s'écoule sur l'aile, il se détache momentanément, créant une petite poche d'air en recirculation appelée « bulle de séparation laminaire » (LSB). Voyez cela comme un minuscule tourbillon temporaire à la surface de l'aile.
• L'instabilité : À l'intérieur de cette bulle, l'air ne reste pas immobile ; il vibre et s'enroule en vagues (comme des rides sur un étang). Les chercheurs ont observé la croissance de ces vagues. Ils ont découvert que le principal « roulement » dans cette bulle est un type d'instabilité appelé mode Kelvin-Helmholtz.
• La découverte : Ils ont confirmé que le programme informatique « plus rapide » (ELLIPSYS) pouvait prédire avec précision la croissance de ces vagues et le comportement de la bulle, correspondant aux résultats du programme ultra-précis.

4. Le pouls lent et rythmique (La grande surprise)

Après avoir validé leurs outils, ils ont laissé la simulation tourner pendant très longtemps (l'équivalent de 50 passages de l'air devant l'aile). C'est là qu'ils ont découvert quelque chose de spécial.

Pendant que l'air s'agitait avec des mouvements rapides et chaotiques, ils ont remarqué un pouls très lent et rythmique dans la force qui pousse sur l'aile.

• L'analogie : Imaginez un battement de tambour. Le mouvement rapide et chaotique de l'air est comme un roulement de tambour rapide et aigu. Le pouls lent qu'ils ont trouvé est comme un battement de cœur profond et lent qui se produit une fois toutes les 48 secondes (en temps de simulation).
• L'effet : Ce battement de cœur lent a fait osciller la force sur l'aile de haut en bas d'environ 10,5 %.
• Le lien avec les éoliennes réelles : Lorsqu'ils ont traduit cela pour une véritable éolienne en rotation, ils ont réalisé que ce pouls lent se produit une fois tous les 7,7 rotations complètes de la pale.

5. Pourquoi cela arrive-t-il ?

Les chercheurs pensent que ce pouls lent est causé par un cycle où l'air « décroche » (se bloque) et « décroche » (se relâche) sur l'aile.

• Le cycle : L'air se bloque, créant une grande bulle. Ensuite, quelque chose déclenche l'éclatement de la bulle, et l'air se réattache de manière fluide. Puis, la pression remonte, la bulle se forme, et le cycle recommence.
• Le déclencheur : Ils soupçonnent que cela se produit parce que l'air tourbillonne vers l'arrière si fortement sur l'aile qu'il crée un état d'« instabilité absolue » — une façon sophistiquée de dire que l'air est si turbulent qu'il ne peut s'empêcher d'osciller de lui-même.

6. L'essentiel

Cette étude est une réussite pour la modélisation informatique. Ils ont prouvé qu'un programme informatique plus rapide et plus efficace (ELLIPSYS) peut être utilisé pour étudier la physique complexe des éoliennes, à condition de le vérifier d'abord par rapport à la « référence absolue ».

Ils ont découvert que même sur une pale d'éolienne épaisse, il existe une « respiration » très lente et rythmique du flux d'air qui se produit environ toutes les 8 rotations de la pale. Cette respiration provoque des variations importantes de la portance (la force qui fait tourner la turbine). Comprendre ce rythme lent est crucial car, bien qu'il ne brise pas la turbine immédiatement, ces oscillations de force lentes et importantes pourraient finir par fatiguer les matériaux après de nombreuses années de fonctionnement.

En bref : Ils ont construit une soufflerie virtuelle, prouvé qu'un ordinateur plus rapide pouvait faire le travail, et découvert que les pales d'éoliennes possèdent un « battement de cœur » lent et rythmé causé par des bulles d'air qui se forment et éclatent sur leur surface.

Résumé technique

Résumé technique : Simulations numériques de la transition et de la réponse à long terme d'un profil d'éolienne

Énoncé du problème
Le décollement de l'écoulement sur les pales d'éoliennes, provoqué par des angles d'incidence élevés et des géométries de profils épais, impacte de manière significative les performances et les charges structurelles. Des études antérieures ont identifié plusieurs fréquences naturelles dans les écoulements décollés, notamment le détachement de sillage, les modes de Kelvin-Helmholtz (KH) et les oscillations à basse fréquence (LFO). Les LFO sont caractérisées par de fortes variations de portance et sont souvent attribuées au décrochage et au rattrapage périodiques de l'écoulement, potentiellement liés à une instabilité absolue au sein des bulles de séparation laminaire (LSB). Bien que les LFO aient été observées dans diverses études sur les profils aérodynamiques, leur occurrence à de faibles angles d'attaque (AoA) sur des profils d'éoliennes épais, ainsi que les mécanismes spécifiques qui les pilotent dans ces conditions, nécessitent des investigations plus approfondies. Cette étude se concentre sur la transition et la réponse à long terme d'une section du profil de l'éolienne de référence DTU 10-MW (série FFA-W3) à un nombre de Reynolds de corde ($Re_c$) de $1 \times 10^5$.

Méthodologie
La recherche utilise des simulations de grande turbulence (LES) avec résolution de paroi, sans modèles de sous-maille, pour résoudre les équations de Navier-Stokes incompressibles. Deux solveurs sont utilisés :

1. NEK5000 : Un code de la méthode des éléments spectraux (SEM) connu pour sa dissipation numérique minimale et sa haute précision.
2. ELLIPSYS : Un code de volumes finis développé à la DTU/Risø, traditionnellement utilisé pour la RANS, mais ici appliqué à l'LES implicite.

L'étude comprend une validation systématique d'ELLIPSYS par rapport à NEK5000. Les étapes méthodologiques clés incluent :

• Sensibilité du domaine : Évaluation de la largeur du domaine en ordonnée ($L_z$) en comparant les simulations avec $L_z = 0,1c$ et $L_z = 0,2c$ pour assurer la résolution des structures d'écoulement pertinentes.
• Validation du solveur : Comparaison des champs d'écoulement moyens, de l'évolution des perturbations et des caractéristiques spectrales entre les deux solveurs.
• Analyse de stabilité : Utilisation des équations de stabilité paraboliques (PSE) et de la décomposition propre orthogonale spectrale (SPOD) pour analyser la croissance des instabilités (modes de Tollmien-Schlichting et KH) et la structure spatiale des modes cohérents.
• Intégration à long terme : Utilisation du solveur ELLIPSYS, plus efficace sur le plan computationnel, pour simuler l'évolution de l'écoulement sur 50 passages de corde ($T = 50c/U_\infty$) afin de capturer les phénomènes de basse fréquence.

Résultats clés

• Largeur du domaine : Une largeur en ordonnée de 10 % de la corde ($L_z = 0,1c$) a été jugée suffisante pour capturer l'évolution des principales perturbations et reproduire les statistiques d'écoulement moyennées, correspondant aux résultats d'un domaine plus large ($L_z = 0,2c$).
• Validation du solveur : ELLIPSYS a démontré une étroite concordance avec NEK5000 pour les champs d'écoulement moyens et l'évolution des perturbations les plus amplifiées. Cependant, ELLIPSYS a sous-estimé l'amplitude des ondes de Tollmien-Schlichting (TS) dans la couche limite attachée en raison d'une dissipation numérique plus élevée. À l'inverse, ELLIPSYS a prédit avec précision l'évolution du mode KH au sein de la bulle de séparation laminaire (LSB), s'alignant bien avec les prédictions de la PSE.
• Mécanisme de transition : La transition sur l'extrados est pilotée par l'instabilité KH se développant sur la couche de cisaillement séparée. Le mode KH le plus amplifié présente une fréquence de $f \approx 15$, et son évolution en aval est bien capturée par les deux solveurs et la PSE.
• Oscillations à basse fréquence (LFO) : Les simulations à long terme ont révélé une modulation lente du coefficient de force normale ($C_N$) avec une amplitude de 10,5 % et une période de 48 passages de corde. Cela correspond à une fréquence $f = 0,021$ et un nombre de Strouhal $St = 0,0012$.
  • Cette fréquence est environ 720 fois plus faible que la fréquence primaire de l'instabilité KH.
  • Dans le contexte de l'éolienne DTU 10-MW, cette période équivaut à environ 7,7 rotations de pale.
  • Le $St$ observé est inférieur aux valeurs typiquement rapportées dans la littérature pour d'autres profils, survenant à un angle d'attaque effectif relativement faible de $3,1^\circ - 3,3^\circ$.

Signification et affirmations
L'article établit qu'ELLIPSYS est un outil viable pour l'LES avec résolution de paroi des écoulements transitoires, offrant une alternative plus efficace sur le plan computationnel aux codes d'éléments spectraux pour les intégrations de longue durée. L'étude confirme que le décrochage et le rattrapage périodiques de l'écoulement, potentiellement déclenchés par une instabilité absolue due à un flux de retour élevé (jusqu'à -17 % sur l'extrados), sont responsables des LFO observées.

Une découverte significative est l'identification de LFO sur un profil épais (épaisseur de 24 %) à un faible angle d'attaque, une condition distincte des décrochages de profils minces typiquement associés à ces phénomènes dans la littérature. Les auteurs notent que le flux de retour sur les deux côtés du profil est suffisant pour permettre une instabilité absolue, ce qui peut piloter l'éclatement périodique de la LSB en un écoulement totalement décroché. L'étude contribue à la compréhension de la manière dont la géométrie épaisse des profils influence l'apparition et les caractéristiques des oscillations à basse fréquence, soulignant que ces phénomènes peuvent survenir à des angles d'incidence et des nombres de Strouhal plus faibles que ceux documentés précédemment pour les profils standards.