Physics-based modelling of wind-turbine wakes turbulence in neutral atmospheric boundary layers

Cet article présente un modèle physique basé sur les budgets d'énergie cinétique turbulente et de contrainte de Reynolds pour prédire l'intensité de turbulence ajoutée par les sillages d'éoliennes, offrant ainsi une approche plus cohérente et prédictive que les modèles empiriques existants.

L'essentiel

🌬️ Le Secret des Sillages des Éoliennes : Une Nouvelle Carte pour Mieux Prévoir le Vent

Imaginez que vous marchez dans un champ de blé par une journée venteuse. Si vous passez devant une personne, elle crée un "trou" dans le vent derrière elle. Les autres personnes qui passent derrière elle recevront moins de vent et seront plus secouées par des rafales désagréables.

C'est exactement ce qui se passe avec les éoliennes. Quand une éolienne tourne pour produire de l'électricité, elle freine le vent. Cela crée un "sillage" (une zone de vent plus faible et plus turbulent) derrière elle. Si vous placez une deuxième éolienne dans ce sillage, elle produira moins d'électricité et s'usera plus vite à cause des secousses.

Le problème, c'est que les ingénieurs ont du mal à prédire exactement à quoi ressemble ce sillage. Les modèles actuels sont souvent des "devinettes" basées sur des formules empiriques (des règles du pouce) ou supposent que le sillage est parfaitement rond et symétrique, comme une boule de fumée. Mais dans la réalité, le vent au-dessus du sol n'est jamais uniforme : il est plus fort en haut, plus faible en bas, et il y a des tourbillons partout.

C'est ici que cette nouvelle étude intervient.

🧱 L'Analogie de la Cuisine : La Recette du Turbulence

Les auteurs de ce papier (des chercheurs de l'IFP Energies nouvelles et de l'Université de Durham) disent : "Arrêtons de deviner, commençons à comprendre la physique."

Imaginez que vous essayez de prédire comment une goutte d'encre se diffuse dans un verre d'eau agitée.

• Les anciennes méthodes disaient : "Regardez comment ça s'est passé la dernière fois, on va faire la même chose." (C'est ce qu'on appelle les modèles empiriques).
• La nouvelle méthode dit : "Analysons la vitesse de l'eau, la taille des tourbillons et la façon dont l'encre se déplace, puis calculons la trajectoire exacte."

Les chercheurs ont créé un modèle basé sur la physique pure. Au lieu de simplement regarder la vitesse du vent (qui manque), ils se sont concentrés sur l'énergie du chaos, ce qu'on appelle la turbulence.

🔍 Comment ils ont fait ? (Le Laboratoire Virtuel)

Pour créer leur modèle, ils n'ont pas seulement utilisé des formules sur un tableau noir. Ils ont construit un laboratoire virtuel géant grâce à des supercalculateurs.

1. La Simulation (LES) : Ils ont simulé le vent soufflant sur une éolienne (représentée comme un disque poreux) dans des conditions très réalistes (vent plus fort en haut, plus faible en bas, sol rugueux ou lisse). C'est comme filmer un film ultra-réaliste de l'écoulement de l'air, pixel par pixel.
2. L'Analyse des Bilans : Ils ont regardé ce qui se passait dans ce film. Ils ont calculé l'équation de la "turbulence" (l'énergie du désordre). Ils ont vu que certaines forces créent du chaos (comme les bords des pales) et d'autres le calment.
3. La Simplification : Ils ont pris ces équations complexes et les ont simplifiées pour qu'elles deviennent une formule mathématique simple et rapide à calculer, tout en restant fidèle à la réalité physique.

🎯 La Grande Révélation : Le Sillage n'est pas une Boule !

Leur découverte principale, c'est que le sillage d'une éolienne n'est pas une boule symétrique.

• Près du sol : Le vent est freiné par le frottement du sol, donc le sillage est étouffé.
• Au niveau des pales : C'est là que le chaos est le plus grand.
• En haut : Le sillage s'étale différemment.

Leur nouveau modèle capture cette asymétrie. Il sait que le vent ne se comporte pas de la même façon en haut et en bas. C'est comme si, au lieu de dessiner un cercle parfait pour le sillage, ils dessinaient une forme de fer à cheval ou de goutte d'eau déformée, beaucoup plus réaliste.

🛠️ Pourquoi est-ce utile pour nous ?

Ce modèle est un outil de précision pour les ingénieurs qui conçoivent les parcs éoliens.

1. Mieux placer les éoliennes : En sachant exactement où le vent va être perturbé, on peut espacer les éoliennes de manière optimale pour maximiser la production d'électricité.
2. Économiser de l'argent : Moins d'erreurs de prédiction signifie moins de pertes d'énergie et moins de pannes dues aux secousses imprévues.
3. Calculer plus vite : Contrairement aux simulations complexes qui prennent des jours, cette nouvelle formule est rapide. On peut l'utiliser pour tester des milliers de configurations de parcs éoliens en quelques secondes.

🏁 En Résumé

Cette étude est comme passer d'une carte dessinée à main levée (les anciens modèles) à un GPS haute précision (le nouveau modèle).

En comprenant la physique profonde de la turbulence dans l'atmosphère, les chercheurs ont créé un outil qui permet de mieux prédire comment le vent interagit avec les éoliennes. C'est une étape clé pour rendre l'énergie éolienne plus efficace, moins chère et plus fiable pour l'avenir.

Le message clé : Ne sous-estimez pas le chaos du vent. En le modélisant avec justesse, on peut mieux l'apprivoiser.

Résumé technique

1. Problématique et Contexte

La conception des parcs éoliens et le contrôle en temps réel reposent souvent sur des modèles de flux simplifiés (analytiques) pour estimer les pertes de charge et la production d'énergie. Ces modèles se composent généralement de deux sous-modèles : un pour le déficit de vitesse et un pour la turbulence ajoutée par la sillage (wake-added turbulence).

• Limites actuelles : Bien que le déficit de vitesse soit bien modélisé, les modèles de turbulence ajoutée sont souvent empiriques ou supposent une symétrie axiale. Or, les éoliennes opèrent dans des couches limites atmosphériques (ABL) cisaillées et turbulentes, ce qui crée des sillages asymétriques (notamment verticalement).
• Besoins : Une estimation précise de l'intensité turbulente locale est cruciale pour prédire les interactions entre les turbines et le rendement énergétique. Les modèles existants ne capturent pas correctement la structure tridimensionnelle de la turbulence ni son évolution dans les couches limites neutres.

2. Méthodologie

Les auteurs proposent un nouveau modèle analytique tridimensionnel basé sur la physique, développé en plusieurs étapes :

A. Génération de données de référence (LES)

• Outil : Utilisation du solveur waLBerla-wind, basé sur la méthode de Boltzmann sur réseau (LBM), pour effectuer des simulations aux grandes échelles (LES).
• Configuration : Disques poreux modélisant des éoliennes dans des couches limites neutres.
• Paramètres : Étude de différentes rugosités de sol (de lisse à rugueux) et de coefficients de poussée ($C_T$). Les simulations couvrent des conditions de bordure de couche atmosphérique (ABL) avec des intensités turbulentes allant de 5 % à 13 %.
• Validation : Les résultats LES sont validés par rapport à des données expérimentales en soufflerie (Bastankhah & Porté-Agel, Stein & Kaltenbach).

B. Analyse budgétaire et simplification des équations

• Équation de transport de l'Énergie Cinétique Turbulente (TKE) : Les auteurs analysent le budget de l'équation de transport de la TKE ajoutée au sillage ($k_w$).
• Simplification : Ils identifient les termes dominants (production, advection, transport turbulent, dissipation) et négligent les termes mineurs (diffusion visqueuse, transport par pression).
• Modélisation des termes :
  • Utilisation de l'hypothèse de diffusion de gradient pour le transport turbulent.
  • Modélisation de la dissipation via une échelle de temps turbulente ($\tau_{k,w}$) plutôt qu'une longueur de mélange, car les données LES montrent que $\tau_{k,w}$ a une forme 3D plus simple.
  • Introduction d'une viscosité turbulente dépendante de l'écoulement pour éviter une dissipation excessive dans le sillage proche.

C. Développement du modèle analytique (Zone lointaine)

• Hypothèses de sillage lointain : Pour obtenir une forme explicite (algébrique) et éviter les schémas numériques coûteux, des hypothèses de similitude sont appliquées à l'advection et au transport dans la zone lointaine ($x > 4D$).
• Formulation finale : Le modèle aboutit à une équation algébrique fermée pour la TKE ajoutée ($k_w$) et la contrainte normale de Reynolds ($u'u'$), reliant ces grandeurs aux gradients de vitesse moyenne, à la viscosité turbulente de l'écoulement incident et à des paramètres calibrés.
• Extension à $u'u'$ : La même méthodologie est appliquée à la contrainte normale de Reynolds, en intégrant un terme de corrélation pression-déformation (modélisé via un terme rapide dominant).

3. Contributions Clés

1. Modèle Physique 3D : C'est le premier modèle analytique basé sur la physique pour la turbulence ajoutée qui prend explicitement en compte l'asymétrie induite par le cisaillement de la couche limite, abandonnant l'hypothèse de symétrie axiale.
2. Approche Hybride : Combinaison rigoureuse de l'analyse budgétaire LES, d'approximations RANS classiques et d'hypothèses de sillage lointain pour dériver une forme analytique simple.
3. Modélisation de la Dissipation : Remplacement des modèles de longueur de mélange par une échelle de temps turbulente ($\tau$) qui s'adapte mieux aux écoulements hors équilibre dans le sillage.
4. Généralité : Le modèle est calibré sur des données LES mais validé sur des données expérimentales "invisibles" (non utilisées pour l'ajustement), démontrant sa robustesse.

4. Résultats

• Comparaison avec LES : Le modèle analytique reproduit avec précision la structure 3D de la TKE et de $u'u'$ dans la zone lointaine. Il capture la transition d'une distribution bimodale (près du rotor) à une distribution plus uniforme (loin du rotor).
  • Note : Une surestimation est observée dans le sillage proche ($x/D < 4$), ce qui est attendu car les hypothèses de similitude ne s'appliquent pleinement qu'en zone lointaine.
• Validation Expérimentale :
  • Cas Bastankhah & Porté-Agel : Le modèle suit très bien les profils latéraux et verticaux de l'intensité turbulente ajoutée, surpassant souvent les modèles empiriques existants (comme le modèle I&Q) qui surestiment la turbulence sous le moyeu.
  • Cas Stein & Kaltenbach : Le modèle reproduit correctement l'effet de la rugosité du sol (cas lisse vs rugueux) sur la contrainte normale de Reynolds, là où d'autres modèles échouent à capturer la magnitude correcte dans les conditions rugueuses.
• Indépendance au Nombre de Reynolds : Des simulations à échelle utilitaire ($D=150$m) confirment que les résultats sont indépendants du nombre de Reynolds dans la plage étudiée.

5. Signification et Perspectives

• Outil Prédictif : Ce modèle offre un outil plus cohérent physiquement et prédictif pour les solveurs de flux de parcs éoliens, permettant une estimation plus précise de la production d'énergie et des charges de fatigue.
• Efficacité Computationnelle : Contrairement aux modèles basés sur l'intégration numérique (comme celui de Bastankhah et al. 2019), ce modèle est analytique et explicite, ce qui le rend idéal pour les tâches intensives en calcul comme l'optimisation de la disposition des parcs éoliens.
• Travaux Futurs : Les auteurs prévoient d'étendre ce modèle pour inclure les effets de la stabilité atmosphérique (stratification thermique) et d'appliquer la méthode à l'échelle de grands parcs éoliens.

En résumé, cet article comble un vide important dans la modélisation éolienne en fournissant une solution analytique tridimensionnelle, physiquement fondée et validée expérimentalement, pour la turbulence ajoutée par les sillages, dépassant les limitations des approches purement empiriques ou symétriques.