Two-way coupling of gravity waves and wind farm wakes: a reduced-order boundary-layer model

Ce papier présente un modèle d'ordre réduit computationnellement efficace qui capture avec succès le couplage bidirectionnel entre les ondes de gravité et les sillages des parcs éoliens en linéarisant les équations de Boussinesq non hydrostatiques et en couplant la dynamique de la couche limite et de l'atmosphère libre au moyen d'une inversion de capping, la validation par rapport à des simulations aux grandes échelles confirmant sa capacité à reproduire des caractéristiques d'écoulement clés telles que le blocage en amont et une récupération accélérée du sillage.

L'essentiel

Imaginez un immense parc éolien non pas comme un simple ensemble d'éoliennes en rotation, mais comme une main géante et invisible tendue vers le ciel, tentant de saisir une poignée de vent pour produire de l'électricité. Cet article porte sur la compréhension de ce qui se produit lorsque cette « main » devient si grande qu'elle ne se contente pas d'aspirer le vent ; elle repousse en réalité l'atmosphère elle-même, créant une danse complexe entre le vent et l'air qui le surmonte.

Voici l'histoire de cette danse, décomposée en concepts simples :

Le Problème : Un Parc Éolien Trop Grand pour ses Brises

Autrefois, les parcs éoliens étaient assez petits pour être comparés à quelques galets dans une rivière. L'eau (le vent) s'écoulait facilement autour d'eux, et la rivière ne s'en préoccupait guère. Mais aujourd'hui, les parcs éoliens sont immenses — parfois aussi hauts que toute la couche d'air dans laquelle nous vivons (la couche limite atmosphérique).

Lorsqu'un parc éolien aussi vaste tente de dérober de l'énergie au vent, il ralentit l'air. Comme l'air ne peut tout simplement pas disparaître, ce ralentissement force l'air à se déplacer vers le haut et vers le bas pour faire de la place. Pensez-y comme à un wagon de métro bondé : si tout le monde s'arrête soudainement de avancer, ils doivent se déplacer vers le haut ou vers le bas pour éviter de se cogner les uns aux autres.

L'Effet « Trampoline » (Ondes de Gravité)

L'atmosphère n'est pas un simple espace vide ; elle possède des couches. Juste au-dessus du parc éolien, il existe un « plafond » distinct appelé l'inversion de capping. Vous pouvez imaginer ce plafond comme un trampoline ou une couverture lourde tendue au-dessus du parc éolien.

Lorsque le parc éolien ralentit l'air, il pousse l'air vers le haut, ce qui creuse une bosse dans ce plafond-trampoline.

1. Le Coup de Pouce : Le parc éolien pousse l'air vers le haut.
2. Le Rebond : Le « trampoline » (l'air stable situé au-dessus) veut se rétracter vers le bas. Ce rebond crée des ondulations, connues sous le nom d'ondes de gravité.
3. La Rétroaction : Ces ondulations ne restent pas simplement là ; elles repoussent le parc éolien vers le bas. C'est comme si le trampoline repoussait vos pieds. Cela crée des changements de pression qui peuvent soit bloquer le vent avant qu'il n'atteigne les turbines (les rendant moins efficaces), soit aider à accélérer le vent derrière le parc (aidant à la récupération du sillage).

L'Ancienne Méthode vs. La Nouvelle Méthode

L'Ancienne Méthode (Le Marteau Lourd) :
Les scientifiques utilisaient auparavant des simulations informatiques ultra-complexes appelées « Simulations des Grandes Échelles » (LES) pour étudier ce phénomène. Imaginez essayer de simuler chaque molécule d'air et chaque minuscule ondulation du trampoline. C'est incroyablement précis, mais cela demande tellement de puissance informatique que c'est comme essayer de compter chaque grain de sable sur une plage juste pour voir comment la marée se déplace. C'est trop lent pour planifier de nouveaux parcs éoliens ou les optimiser en temps réel.

La Nouvelle Méthode (Le Croquis Intelligent) :
Les auteurs de cet article ont créé un « modèle d'ordre réduit ». Imaginez cela comme un croquis intelligent plutôt qu'une peinture photoréaliste.

• Ils ont simplifié les mathématiques en se concentrant uniquement sur les parties les plus importantes : le mouvement vertical de l'air et les ondulations sur le « trampoline ».
• Ils ont traité le parc éolien comme une force continue plutôt que de simuler chaque pale de turbine individuellement.
• Ils ont utilisé un astucieux tour de mathématiques (mélangeant des méthodes spectrales et des différences finies) pour résoudre les équations rapidement.

Ce Qu'ils Ont Découvert

Ils ont testé leur « croquis intelligent » contre le « marteau lourd » (les simulations ultra-complexes) et des données réelles. Voici ce qu'ils ont découvert :

1. Le Blocage : Lorsque le parc éolien se trouve dans une atmosphère stable (comme une journée calme avec un « plafond » clair), les ondes de gravité créent un « vent de face » avant même que le parc ne commence. C'est comme essayer de courir contre un fort vent de face qui se forme avant même que vous n'atteigniez l'obstacle. Cela ralentit considérablement le vent avant qu'il ne frappe les turbines.
2. La Récupération : Derrière le parc éolien, le « trampoline » se rétracte vers le bas, créant un « vent arrière » qui pousse l'air vers l'avant. Cela aide la vitesse du vent à se rétablir beaucoup plus rapidement qu'elle ne le ferait par une journée calme et neutre.
3. Précision : Leur modèle simplifié correspondait presque parfaitement aux résultats des simulations ultra-complexes, mais il s'exécutait des milliers de fois plus vite.

La Conclusion

Cet article offre aux ingénieurs un outil rapide et fiable pour prédire comment les immenses parcs éoliens interagiront avec le ciel. Au lieu d'attendre des jours qu'un supercalculateur leur dise comment un parc se comportera, ils peuvent désormais utiliser ce modèle pour voir en quelques secondes comment l'effet « trampoline » de l'atmosphère aidera ou entravera le parc éolien. Il comble le fossé entre les suppositions simples et les super-simulations impossibles à exécuter, nous aidant à concevoir de meilleurs parcs éoliens qui fonctionnent avec l'atmosphère, et non simplement contre elle.

Résumé technique

Résumé technique : Couplage bidirectionnel des ondes de gravité et des sillages de parcs éoliens : un modèle de couche limite à ordre réduit

Énoncé du problème
L'expansion rapide des parcs éoliens à l'échelle des services publics a donné naissance à des installations dont la taille est comparable à l'épaisseur de la couche limite atmosphérique (CLA). Ce changement d'échelle crée une interaction bidirectionnelle complexe entre le parc éolien et l'atmosphère, particulièrement sous stratification thermique. Dans une couche limite conventionnellement neutre (CLCN), caractérisée par une inversion de couverture et une atmosphère libre stratifiée, l'extraction d'énergie par les turbines induit des mouvements verticaux pour satisfaire la continuité. Ces mouvements déplacent l'inversion de couverture, générant des ondes de gravité à grande échelle (à la fois interfaciales et internes). Ces ondes, à leur tour, induisent des gradients de pression qui rétroagissent sur l'écoulement de la couche limite, provoquant un blocage significatif en amont et modifiant les taux de récupération du sillage en aval.

Bien que la Simulation des Grandes Échelles (SGE) puisse capturer ces phénomènes, elle est coûteuse en calcul en raison des domaines spatiaux massifs requis pour la propagation des ondes de gravité et de la sensibilité de ces ondes aux conditions aux limites. Les modèles à ordre réduits précédents ont tenté de résoudre ce problème mais souffrent de limitations : certains supposent des champs d'écoulement verticalement uniformes, tandis que d'autres couplent des modèles d'ondes linéarisés avec des modèles de sillage d'ingénierie, aboutissant à des systèmes hautement paramétrés qui dépendent de choix empiriques pour la diffusion turbulente, les taux de récupération du sillage et les méthodes de superposition. Il subsiste un besoin pour un cadre unifié qui résout la structure verticale continue de la CLA tout en maintenant l'efficacité computationnelle et en réduisant la dépendance aux paramètres.

Méthodologie
Les auteurs développent un modèle numérique linéarisé à ordre réduit qui résout la structure verticale continue de la CLA et capture le couplage bidirectionnel entre les écoulements de sillage et les ondes de gravité. L'approche comprend les étapes suivantes :

1. Équations gouvernantes : Le modèle linéarise les équations de Boussinesq bidimensionnelles, stationnaires et non hydrostatiques autour d'un écoulement de fond avec une vitesse horizontale $U(z)$. Les équations sont séparées en deux régions : la CLA (où existent la turbulence et le forçage du parc) et l'atmosphère libre (où la stratification existe mais où la turbulence et le forçage sont négligés).
2. Formulation de la couche limite : Au sein de la CLA, le modèle dérive une équation gouvernante pour la perturbation de vitesse verticale ($w$) en éliminant la vitesse horizontale et la pression. La contrainte de cisaillement turbulente est modélisée à l'aide de l'hypothèse de Boussinesq avec une approche de longueur de mélange.
3. Formulation de l'atmosphère libre : Dans l'atmosphère libre, les équations se réduisent à une équation de Helmholtz pour la vitesse verticale, résolue via une transformée de Fourier pour obtenir des solutions pour les ondes internes et interfaciales.
4. Mécanisme de couplage : Les deux régions sont couplées à l'inversion de couverture ($z=H$) par des conditions aux limites cinématiques et dynamiques.
   • Cinématique : La couche d'inversion se déplace avec le fluide, reliant la vitesse verticale au déplacement de l'inversion ($\eta$).
   • Dynamique : Une condition de continuité de pression est dérivée en comparant les pressions avant et après le déplacement de l'inversion. Cela produit une condition aux limites clé (Équation 2.19) qui relie le gradient de vitesse verticale au sommet de la CLA à la vitesse verticale elle-même. Cette condition englobe les effets des ondes interfaciales (via la gravité réduite $g_r$) et des ondes internes (via la fréquence de flottabilité $N$).
5. Discrétisation numérique : Le système est résolu en utilisant une méthode mixte spectrale aux différences finies. Le domaine vertical est discrétisé par différences finies, tandis que la direction horizontale est traitée par transformées de Fourier. Cette approche est efficace en calcul et accueille naturellement des domaines périodiques.
6. Forçage du parc éolien : Le parc est modélisé comme un réseau semi-infini de disques actionneurs. Le terme de forçage est calculé itérativement pour tenir compte du déficit de sillage local et du rapport entre les effets de sillage locaux et moyennés latéralement ($\beta_n$).

Résultats clés
Le modèle a été validé contre trois ensembles de données SGE distincts représentant différents régimes atmosphériques (stratification et nombres de Froude variables) :

• Comportement qualitatif : Les comparaisons entre une couche limite véritablement neutre (CLVN) et une CLCN démontrent que les ondes de gravité modifient considérablement la dynamique de l'écoulement. Dans la CLCN, le modèle reproduit :
  • Blocage en amont : Un gradient de pression adverse se développe en amont du parc, provoquant de forts effets de blocage global absents dans le cas CLVN.
  • Récupération du sillage : Un gradient de pression favorable se forme à l'intérieur et en aval du parc, accélérant la récupération du sillage par rapport au cas CLVN.
  • Déplacement vertical : Le modèle capture le déplacement de l'inversion de couverture, montrant que pour des écoulements subcritiques (stratification forte), le déplacement maximal se déplace vers le début du parc, tandis que pour des écoulements supercritiques, il se produit près de la fin.
• Accord quantitatif :
  • Le modèle montre un bon accord avec les données SGE pour les perturbations de vitesse dans le sens du vent à hauteur de moyeu et les vitesses moyennées verticalement.
  • Il reproduit avec succès l'amplitude du blocage en amont et l'accélération de la récupération du sillage dans les cas subcritiques.
  • Dans le cas subcritique fortement stratifié (Q00), le modèle surestime légèrement la récupération du sillage (entraînant une petite perturbation de vitesse positive), ce que les auteurs attribuent à la simplification 2D qui pourrait surestimer le gradient de pression favorable induit par les ondes de gravité.
• Efficacité computationnelle : Le modèle est efficace en calcul, chaque itération prenant environ 3 secondes sur un PC standard sans parallélisation, ce qui le rend adapté aux balayages de paramètres et aux tâches d'optimisation où la SGE est prohibitive.

Portée et affirmations
L'article prétend fournir un cadre unifié qui comble le fossé entre la SGE coûteuse en calcul et les modèles d'ingénierie hautement paramétrés. Ses contributions principales sont :

1. Fidélité physique : Il résout la structure verticale continue de la CLA sans s'appuyer sur la « discrétisation en couches » (par exemple, diviser la CLA en couches distinctes de parc et supérieures) utilisée dans les approches à ordre réduit précédentes.
2. Réduction de la paramétrisation : En dérivant le couplage dynamiquement à partir de premiers principes (physique de l'inversion de couverture) plutôt qu'en couplant des modèles d'ondes linéaires avec des modèles de sillage empiriques, le cadre réduit le nombre de paramètres d'ingénierie arbitraires (tels que les taux de récupération du sillage ou les méthodes de superposition).
3. Couplage bidirectionnel : Le modèle capture naturellement la rétroaction des ondes de gravité sur l'écoulement de la couche limite grâce à la condition aux limites dérivée, reproduisant avec précision à la fois le blocage en amont et l'accélération du sillage en aval.

Les auteurs positionnent ce cadre comme une fondation pour les futures paramétrisations des effets des parcs éoliens dans les modèles météorologiques et climatiques, et comme un outil pour l'évaluation efficace des performances des parcs éoliens et des interactions entre parcs. Ils notent que les travaux futurs étendront le cadre à trois dimensions et intégreront des physiques supplémentaires, y compris les effets non linéaires, le forçage de Coriolis et des représentations améliorées de la turbulence.