Generazione di profili del vento casuali e correlati spazialmente per simulazioni CFD del vento

1. Introduzione

Una valutazione accurata delle azioni del vento sulle strutture dipende non solo dal modello di carico applicato, ma fondamentalmente da come viene rappresentato lo strato limite atmosferico (ABL) in arrivo. Nel design ingegneristico pratico, il vento è un fenomeno altamente instabile, turbolento e spazialmente correlato. Tuttavia, la maggior parte degli approcci basati sui codici, incluso il modello di vento standard Eurocode, riducono questa complessità a quantità statiche equivalenti per garantire semplicità, robustezza e margini di sicurezza conservativi.

Nelle simulazioni di vento basate su Computational Fluid Dynamics (CFD), la definizione della condizione al contorno dell'inlet è un fattore decisivo per l'accuratezza delle pressioni superficiali previste, dei carichi di picco e degli effetti dinamici.

Il metodo documentato qui implementa un generatore di afflusso turbolento casuale e spazialmente correlato su un piano di ingresso verticale (piano YZ) ed è destinato a:

• Simulazioni transitorie di vento CFD

• Valutazione delle pressioni e delle forze rilevanti per il carico

• Validazione rispetto a condizioni di afflusso simili a quelle di una galleria del vento

• Applicazioni avanzate al di là dell'ambito conservativo dell'Eurocode

Questo approccio va significativamente oltre il modello di vento standard dell'Eurocode (EC), che è progettato principalmente per la derivazione di carichi statici di design, non per la fisica del flusso risolta nel tempo.

2. Panoramica del Metodo Implementato

La generazione dell'afflusso consiste in cinque componenti accoppiate:

2.1 Piano di Ingresso Discretizzato Spazialmente

• Ingresso definito su un piano YZ a x costante = x₀

• Griglia regolare: n_y × n_z punti (es. 20 × 100 → 2.000 sonde)

• Ogni punto della griglia rappresenta un segnale di velocità dipendente dal tempo

• Non viene applicata alcuna media spanwise (importante per preservare le fluttuazioni)

Questo permette di rappresentare fisicamente gli effetti di coerenza verticale e laterale.

2.2 Profilo del Vento Medio (Legge di Potenza, Applicazione Fisica)

Il profilo del vento medio è definito utilizzando una formulazione a legge di potenza, un'approssimazione ingegneristica empirica che rappresenta l'aumento della velocità del vento con l'altezza all'interno dello strato limite atmosferico. Il modello collega la velocità a una condizione di vento di riferimento e un esponente dipendente dal terreno. Sebbene più semplice della formulazione del logaritmo dello strato limite derivata dalla teoria della similarità, l'approccio a legge di potenza fornisce una definizione di afflusso pratica e ampiamente utilizzata per le applicazioni ingegneristiche. La formulazione garantisce una distribuzione verticale della velocità coerente riferita a un'altezza specificata.

2.3 Modello di Turbolenza: Spettro di Kaimal

La turbolenza in ciascun punto dell'inlet è generata utilizzando lo spettro di Kaimal, che rappresenta il comportamento dello strato limite atmosferico e definisce come viene distribuita l'energia turbolenta tra le frequenze. Questo approccio cattura le caratteristiche realistiche dipendenti dalla frequenza delle fluttuazioni del vento piuttosto che assumere un comportamento costante. Il modello riflette l'intensità della turbolenza e le scale di lunghezza caratteristiche per descrivere la struttura delle variazioni di velocità. Come risultato, la turbolenza generata include contenuti energetici fisicamente coerenti che variano con la frequenza. In ciascun punto dell'inlet, la turbolenza è generata utilizzando l'autospettro di Kaimal, coerente con la teoria dello strato limite atmosferico [1]:

Questo produce contenuti energetici dipendenti dalla frequenza, assenti nell'approccio statico dell'Eurocode.

2.4 Coerenza Spaziale e Matrice Spettrale di Croce

Un modello di coerenza spaziale è applicato per garantire una correlazione fisicamente realistica tra i punti dell'inlet descrivendo come la similarità della turbolenza diminuisca con la distanza e la frequenza. Questa funzione di coerenza è utilizzata per definire la relazione tra segnali in posizioni diverse, riflettendo fluttuazioni del vento spazialmente correlate. Basandosi su questa coerenza, viene costruita una matrice densità di spettri di croce completa, combinando spettri individuali con la loro correlazione spaziale. La matrice risultante rappresenta sia la distribuzione dell'energia sia il collegamento spaziale delle fluttuazioni di velocità attraverso l'inlet. Per garantire una correlazione fisicamente realistica tra i punti dell'inlet, viene applicato un modello di coerenza dipendente dalla distanza:

Questo risulta in una matrice densità di spettri di croce completa:

Caratteristiche principali:

• Turbolenza completamente correlata spazialmente

• Decadimento della correlazione dipendente dalla frequenza

• Cattura strutture realistiche delle raffiche attraverso l'inlet

2.5 Decomposizione di Cholesky Robusta (Stabilità Numerica)

Una decomposizione di Cholesky robusta è applicata per mantenere la stabilità numerica quando la matrice spettrale di croce diventa mal condizionata, specialmente vicino al suolo. La matrice è prima simmetrizzata per garantire coerenza, seguita da un caricamento diagonale adattivo per migliorare la stabilità. Una strategia di riprova automatica con regolarizzazione crescente è usata se necessario, assicurando una definitezza positiva. Questo approccio permette una generazione stabile e affidabile di dati di serie temporali correlati. Poiché la matrice spettrale di croce può diventare mal condizionata (specialmente vicino al suolo), viene usata una fattorizzazione robusta di Cholesky:

• Applicazione dell'Hermitian:

• Caricamento diagonale adattivo

• Strategia di riprova automatica con regolarizzazione crescente

Questo garantisce una definitezza positiva e una generazione stabile delle serie temporali.

2.6 Costruzione delle Serie Temporali (FFT Inversa)

Fasi casuali sono introdotte per rappresentare la turbolenza all'interno del segnale, mentre lo spettro è scalato e strutturato per garantire un risultato reale fisicamente realistico. Un'inversa FFT viene quindi applicata per trasformare l'informazione nel dominio della frequenza in un segnale di velocità dipendente dal tempo. Il risultato finale è un campo di velocità composto da un flusso medio combinato con componenti fluttuanti.

• Fasi casuali applicate per frequenza

• Spettro monolato scalato da sqrt(2Δ𝑓)

• Simmetria Hermitian applicata

• FFT Inversa → segnale di velocità risolto nel tempo in ogni punto dell'inlet

Risultato:

3. Visualizzazione e Diagnostica

Il metodo consente:

• Profili verticali istantanei 𝑈(𝑧,𝑡_o)

• Confronti di più linee di sonde (nessun artefatto di mediazione)

• Animazioni complete di contorno del piano di ingresso 𝑈(𝑦,𝑧,𝑡)

Questo rende osservabile direttamente la fisica del flusso, a differenza dei profili statici basati sull'EC.

1 Distribuzione dei punti di ingresso sul piano YZ con griglia spaziale utilizzata per la generazione di un flusso di turbolenza casuale

2 Confronto del profilo del vento medio e dei profili di fluttuazione istantanea in più posizioni con flessioni

3 Campo di vento all’ingresso verticale istantaneo che mostra la distribuzione della velocita turbolenta correlata nello spazio

4. Confronto con il Modello di Vento Standard Eurocode

La tabella 1 riassume le differenze concettuali e metodologiche fondamentali tra il modello di vento standard Eurocode (EC) e il corrente metodo di generazione dell'inflow CFD casuale. Mentre entrambi gli approcci mirano a rappresentare gli effetti del vento sulle strutture, essi si basano su filosofie di modellazione differenti e servono scopi ingegneristici diversi.

Tabella 1: Confronto con il Modello di Vento Standard Eurocode

5. Differenze Chiave nell'Interpretazione

L'approccio Eurocode è principalmente destinato a garantire carichi di progetto sicuri incorporando effetti di turbolenza tramite fattori di sicurezza parziali, ma non fornisce approfondimenti dettagliati sul timing o sui meccanismi dello sviluppo del carico. Al contrario, il metodo attuale permette una comprensione fisica più profonda risolvendo quando, dove e perché si verificano picchi di pressione, inclusi penetrazione delle raffiche, picchi di pressione indotti dalla separazione e effetti di carico spazialmente correlati. Questo lo rende particolarmente prezioso per strutture leggere, elementi di facciata e rivestimento, sistemi a membrana, apparecchiature montate sul tetto e studi che coinvolgono una risposta strutturale dinamica:

📘Eurocodice

• Progettato per produrre carichi di progetto sicuri

• Effetti di turbolenza incorporati nei fattori di sicurezza parziali

• Nessuna informazione su quando o come si verificano i carichi

🌬️Metodo Corrente

• Risolve quando, dove e perché si formano picchi di pressione
• Penetrazione delle raffiche
• Picchi di pressione indotti dalla separazione
• Caricamento correlato su grande scala

⭐Particolarmente importante per:

• Strutture leggere

• Elementi di rivestimento e facciata

• Membrane e attrezzature sul tetto

• Studi di risposta dinamica

6. Implicazioni Ingegneristiche

Il metodo proposto non mira a sostituire l'Eurocode per il progetto normativo; invece, dovrebbe essere inteso come un approccio complementare che fornisce approfondimenti fisici aggiuntivi oltre alle assunzioni tipicamente incorporate nelle procedure Eurocode. È particolarmente adatto per studi di sensibilità, scopi di ricerca e validazione, processi di ottimizzazione e per spiegare discrepanze che possono sorgere tra i risultati basati sull'Eurocode e le simulazioni CFD. Nelle applicazioni ingegneristiche pratiche, questo approccio spesso aiuta a chiarire perché i valori di pressione derivati da CFD differiscono dalle previsioni Eurocode, specialmente riguardo ai picchi di pressione locali e agli effetti di flusso transitorio, senza necessariamente indicare un errore in una metodologia o nell'altra.

7. Riassunto

Il metodo di generazione dell'inflow presentato introduce la modellazione della turbolenza fisicamente coerente combinata con la coerenza spaziale e la rappresentazione del campo di vento risolto nel tempo, supportata da un'implementazione numerica robusta. Rispetto all'approccio standard Eurocode, il metodo sposta la prospettiva analitica dagli involucri di carico conservativi verso una comprensione guidata dalla fisica del comportamento del flusso, abilitando in tal modo approfondimenti più profondi e analisi ingegneristiche del vento più avanzate.