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

Questo articolo presenta simulazioni numeriche utilizzando Nek5000 ed EllipSys per analizzare la transizione e la risposta del flusso a lungo termine di un profilo alare FFA-W3 a bassi numeri di Reynolds, convalidando la capacità del codice EllipSys di catturare le instabilità di Kelvin-Helmholtz e identificando una modulazione lenta e periodica del coefficiente di forza normale legata a oscillazioni a bassa frequenza e al potenziale scoppio della bolla.

L'essenziale

Immaginate una pala di una turbina eolica come un'ala gigante che ruota. Proprio come l'ala di un aeroplano, ha bisogno di un flusso d'aria fluido per funzionare in modo efficiente. Ma quando il vento colpisce l'ala con determinate angolazioni, l'aria può "incastrarsi" e separarsi dalla superficie, creando un groviglio caotico e vorticoso. Questo articolo è come un esperimento in una galleria del vento ad alta tecnologia, ma invece di usare un modello fisico, i ricercatori ne hanno costruito uno virtuale all'interno di un supercomputer per osservare esattamente come si comporta l'aria.

Ecco la storia della loro scoperta, suddivisa in concetti semplici:

1. La galleria del vento virtuale

I ricercatori volevano studiare una sezione specifica di una massiccia pala di una turbina eolica (da una turbina da 10 MW). Hanno utilizzato due diversi programmi per computer, NEK5000 ed ELLIPSYS, per simulare il flusso d'aria sopra questa pala.

Pensate a NEK5000 come a una telecamera di alta gamma, ultra-precisa, che cattura ogni minimo dettaglio ma è molto lenta e costosa da utilizzare. ELLIPSYS è come una telecamera leggermente più veloce ed efficiente. Il team ha prima dovuto dimostrare che la telecamera "più veloce" (ELLIPSYS) potesse vedere le stesse cose di quella "di alta gamma". Hanno scoperto che, sebbene ELLIPSYS perdesse alcune increspature minuscole e deboli nell'aria fluida (perché la rendeva un po' troppo levigata), era eccellente nel catturare i grandi vortici caotici che sono fondamentali per le prestazioni della pala.

2. Quanto deve essere larga la galleria?

Prima di avviare le simulazioni lunghe, dovevano capire quanto dovesse essere larga la loro "galleria del vento" virtuale. Se la galleria è troppo stretta, potrebbe comprimere l'aria e creare risultati falsi. Se è troppo larga, spreca potenza di calcolo.

Hanno testato una galleria "stretta" (il 10% della larghezza dell'ala) contro una galleria "larga" (il 20% della larghezza).

• L'analogia: Immaginate di osservare il corso di un fiume. Se guardate solo una striscia stretta del fiume, perderete le grandi onde?
• Il risultato: Hanno scoperto che la galleria stretta era in realtà sufficiente. Le grandi onde e i vortici si formavano perfettamente nello spazio ristretto. Ciò significava che potevano risparmiare molto tempo di calcolo utilizzando la scatola di simulazione più piccola e stretta.

3. La "bolla" e lo "flap"

La parte più interessante dello studio è avvenuta sulla parte superiore dell'ala (il lato di aspirazione).

• La bolla di separazione: Mentre l'aria scorre sull'ala, si stacca per un momento, creando una piccola sacca di aria che ricircola, chiamata "bolla di separazione laminare" (Laminar Separation Bubble - LSB). Pensate a questo come a un piccolo, temporaneo vortice sulla superficie dell'ala.
• L'instabilità: All'interno di questa bolla, l'aria non sta ferma; vibra e si arrotola in onde (come increspature su uno stagno). I ricercatori hanno osservato queste onde crescere. Hanno scoperto che il principale "rotolamento" in questa bolla è un tipo di instabilità chiamata modo Kelvin-Helmholtz.
• La scoperta: Hanno confermato che il programma per computer "più veloce" (ELLIPSYS) poteva prevedere accuratamente come queste onde crescessero e come si comportasse la bolla, corrispondendo ai risultati del programma ultra-preciso.

4. Il impulso lento e ritmico (La grande sorpresa)

Dopo aver validato i loro strumenti, hanno lasciato che la simulazione girasse per un tempo molto lungo (equivalente a 50 volte il passaggio dell'aria accanto all'ala). È qui che hanno scoperto qualcosa di speciale.

Mentre l'aria era in tumulto con movimenti rapidi e caotici, hanno notato un impulso ritmico molto lento nella forza che spingeva sull'ala.

• L'analogia: Immaginate un rullante. Il tumulto rapido dell'aria è come un rullo veloce e acuto. L'impulso lento che hanno trovato è come un battito cardiaco profondo e lento che avviene una volta ogni 48 secondi (nel tempo della simulazione).
• L'effetto: Questo battito cardiaco lento ha fatto sì che la forza sull'ala oscillasse su e giù di circa il 10,5%.
• La connessione con le turbine reali: Quando hanno tradotto questo dato su una vera turbina eolica rotante, hanno capito che questo impulso lento avviene una volta ogni 7,7 rotazioni complete della pala.

5. Perché succede questo?

I ricercatori ritengono che questo impulso lento sia causato da un ciclo in cui l'aria "stalla" (si blocca) e "destalla" (si rilascia) sull'ala.

• Il ciclo: L'aria si blocca, creando una grande bolla. Poi, qualcosa innesca lo scoppio della bolla, e l'aria si riattacca fluidamente. Successivamente, la pressione aumenta di nuovo, la bolla si forma e il ciclo si ripete.
• L'innesco: Sospettano che ciò accada perché l'aria sta ruotando all'indietro con tale forza sulla superficie dell'ala da creare uno stato di "instabilità assoluta", un modo elaborato per dire che l'aria è così turbolenta da non poter fare a meno di oscillare autonomamente.

6. In sintesi

Questo studio è un successo per la modellazione al computer. Hanno dimostato che un programma per computer più veloce ed efficiente (ELLIPSYS) può essere affidabile per studiare la complessa fisica delle turbine eoliche, a patto di verificarlo prima rispetto allo "standard di riferimento".

Hanno scoperto che, anche su una pala di turbina eolica spessa, esiste un "respiro" lento e ritmico del flusso d'aria che avviene circa ogni 8 rotazioni della pala. Questo respiro fa sì che la portanza (la forza che fa girare la turbina) salga e scenda significativamente. Comprendere questo ritmo lento è fondamentale perché, sebbene non possa rompere la turbina immediatamente, queste ampie oscillazioni di forza potrebbero alla fine affaticare i materiali nel corso di molti anni di operatività.

In breve: Hanno costruito una galleria del vento virtuale, hanno dimostrato che un computer veloce può svolgere il lavoro e hanno scoperto che le pale delle turbine eoliche hanno un "battito cardiaco" lento e ritmico causato dalla formazione e dallo scoppio di bolle d'aria sulla loro superficie.

Sommario tecnico

Sintesi Tecnica: Simulazioni Numeriche della Transizione e della Risposta a Lungo Termine di un Profilo Alare per Turbine Eoliche

Definizione del Problema
Il distacco del flusso sulle pale delle turbine eoliche, causato da elevati angoli di incidenza del flusso e geometrie di profili alari spessi, impatta significativamente le prestazioni e i carichi strutturali. Studi precedenti hanno identificato molteplici frequenze naturali nei flussi separati, tra cui lo shedding della scia, i modi di Kelvin-Helmholtz (KH) e le oscillazioni a bassa frequenza (LFO). Le LFO sono caratterizzate da forti variazioni di portanza e sono spesso attribuite allo stalling e unstalling periodico del flusso, potenzialmente legati all'instabilità assoluta all'interno delle bolle di separazione laminare (LSB). Sebbene le LFO siano state osservate in vari studi su profili alari, la loro occorrenza a bassi angoli di attacco (AoA) su profili alari spessi per turbine eoliche, e i meccanismi specifici che le guidano in tali condizioni, richiedono ulteriori indagini. Questo studio si concentra sulla transizione e sulla risposta a lungo termine di una sezione del profilo alare della turbina di riferimento DTU 10-MW (serie FFA-W3) a un numero di Reynolds di corda ($Re_c$) pari a $1 \times 10^5$.

Metodologia
La ricerca impiega simulazioni Large Eddy (LES) con risoluzione della parete (wall-resolved) senza modelli di sottogriglia per risolvere le equazioni di Navier-Stokes incomprimibili. Vengono utilizzati due solver:

1. NEK5000: Un codice a elementi spettrali (SEM) noto per la minima dissipazione numerica e l'alta accuratezza.
2. ELLIPSYS: Un codice a volumi finiti sviluppato presso DTU/Risø, tradizionalmente utilizzato per RANS, ma qui applicato alla LES implicita.

Lo studio prevede una validazione sistematica di ELLIPSYS rispetto a NEK5000. Le fasi metodologiche chiave includono:

• Sensibilità del Dominio: Valutazione della larghezza del dominio spanwise ($L_z$) confrontando le simulazioni con $L_z = 0.1c$ e $L_z = 0.2c$ per garantire la risoluzione delle strutture di flusso rilevanti.
• Validazione del Solver: Confronto dei campi di flusso medi, dell'evoluzione delle perturbazioni e delle caratteristiche spettrali tra i due solver.
• Analisi di Stabilità: Utilizzo delle Equazioni di Stabilità Parabolizzate (PSE) e della Decomposizione Spettrale delle Proprietà Ortogonali (SPOD) per analizzare la crescita delle instabilità (modi di Tollmien-Schlichting e KH) e la struttura spaziale dei modi coerenti.
• Integrazione a Lungo Termine: Utilizzo del solver ELLIPSYS, computazionalmente più efficiente, per simulare l'evoluzione del flusso su 50 passaggi di corda ($T = 50c/U_\infty$) per catturare i fenomeni a bassa frequenza.

Risultati Chiave

• Larghezza del Dominio: Una larghezza spanwise pari al 10% della corda ($L_z = 0.1c$) è stata ritenuta sufficiente per catturare l'evoluzione delle principali perturbazioni e riprodurre le statistiche del flusso mediato nel tempo, concordando con i risultati di un dominio più ampio ($L_z = 0.2c$).
• Validazione del Solver: ELLIPSYS ha dimostrato una stretta concordanza con NEK5000 per i campi di flusso medi e l'evoluzione delle perturbazioni più amplificate. Tuttavia, ELLIPSYS ha sottostimato l'ampiezza delle onde di Tollmien-Schlichting (TS) nello strato limite attaccato a causa di una maggiore dissipazione numerica. Al contrario, ELLIPSYS ha predetto accuratamente l'evoluzione del modo KH all'interno della bolla di separazione laminare (LSB), allineandosi bene con le previsioni PSE.
• Meccanismo di Transizione: La transizione sul lato ventrale (suction side) è guidata dall'instabilità KH che si sviluppa sullo strato di taglio separato. Il modo KH più amplificato ha una frequenza di $f \approx 15$, e la sua evoluzione a valle è ben catturata da entrambi i solver e dalla PSE.
• Oscillazioni a Bassa Frequenza (LFO): Le simulazioni a lungo termine hanno rivelato una lenta modulazione del coefficiente di forza normale ($C_N$) con un'ampiezza del 10.5% e un periodo di 48 passaggi di corda. Ciò corrisponde a una frequenza $f = 0.021$ e a un numero di Strouhal $St = 0.0012$.
  • Questa frequenza è circa 720 volte inferiore alla frequenza primaria dell'instabilità KH.
  • Nel contesto della turbina DTU 10-MW, questo periodo equivale a circa 7.7 rotazioni della pala.
  • Lo $St$ osservato è inferiore ai valori tipicamente riportati in letteratura per altri profili alari, verificandosi a un angolo di attacco efficace relativamente piccolo ($3.1^\circ - 3.3^\circ$).

Significatività e Rivendicazioni
Il documento stabilisce che ELLIPSYS è uno strumento valido per la LES wall-resolved di flussi transizionali, offrendo un'alternativa computazionalmente efficiente ai codici a elementi spettrali per integrazioni a lungo termine. Lo studio conferma che lo stalling e l'unstalling periodico del flusso, potenzialmente innescati dall'instabilità assoluta dovuta all'elevato flusso inverso (fino al -17% sul lato ventrale), è responsabile delle LFO osservate.

Un risultato significativo è l'identificazione delle LFO su un profilo alare spesso (spessore 24%) a un basso angolo di attacco, una condizione distinta dagli stall di profili sottili tipicamente associati a tali fenomeni nella letteratura scientifica. Gli autori osservano che il flusso inverso su entrambi i lati del profilo è sufficiente a permettere l'instabilità assoluta, che può guidare il bursting periodico della LSB in un flusso completamente stallato. Lo studio contribuisce alla comprensione di come la geometria del profilo alare influenzi l'insorgenza e le caratteristiche delle oscillazioni a bassa frequenza, evidenziando come tali fenomeni possano verificarsi ad angoli di incidenza e numeri di Strouhal inferiori rispetto a quanto documentato precedentemente per profili standard.