FluidodinamicaProcesso e acqua

Fluidodinamica computazionale

Risolvere un’equazione differenziale alle derivate parziali su una griglia — un campo di velocità nello spazio e nel tempo. Un risolutore Lattice-Boltzmann autonomo riproduce il risultato più iconico della fluidodinamica: la scia di vortici di von Kármán generata dal flusso attorno a un cilindro, validata su flusso in canale, numero di Strouhal, resistenza e portanza, e una cavità trascinata dal coperchio.

DjiniousLabModellazione e simulazione

Un campo di vorticità che mostra vortici rossi e blu alternati che si staccano da un cilindro in una scia di von Kármán

errore del flusso in canale contro soluzione analitica: 0.27%errore del flusso in canale contro soluzione analitica
numero di Strouhal (frequenza di distacco): 0.20numero di Strouhal (frequenza di distacco)
coefficiente di resistenza del cilindro: 0.94coefficiente di resistenza del cilindro
la resistenza oscilla al doppio della frequenza della portanza: 2×la resistenza oscilla al doppio della frequenza della portanza
notebook di progettazione: 10notebook di progettazione
requisiti PASS: 5 / 5requisiti PASS

Un intero campo di velocità, non una manciata di stati.

Ogni altro modello di questo portafoglio traccia poche variabili nel tempo. Il flusso di un fluido è diverso: è un campo — velocità e pressione in ogni punto dello spazio, accoppiate da un’equazione differenziale alle derivate parziali, che evolvono insieme. Risolverlo su una griglia è la fluidodinamica computazionale, il dominio di COMSOL e ANSYS Fluent. Questo programma costruisce da zero un risolutore autonomo e lo usa per riprodurre il singolo risultato più riconoscibile del campo — la scia di vortici che si stacca da un cilindro in una corrente fluente — poi lo valida contro benchmark da manuale.

Ottieni la fisica giusta su un problema che puoi verificare.

Il risolutore è un metodo Lattice-Boltzmann — una via cinetica, a distribuzione di particelle, alle equazioni di Navier-Stokes incomprimibili, esplicita, locale e robusta dove una risoluzione diretta della pressione è delicata. Prima di fidarsene su qualcosa di complesso, il programma lo valida sul flusso in canale guidato da pressione, dove la risposta esatta è una parabola: il profilo simulato corrisponde a quello analitico entro lo 0,27%. Un risolutore che centra il caso che puoi verificare sulla carta è uno di cui puoi fidarti sul caso che non puoi.

Un profilo di velocità del flusso in canale simulato sovrapposto alla parabola analitica — L’ancora di correttezza: il flusso in canale guidato da pressione (di Poiseuille) ha un profilo di velocità parabolico esatto, e il risolutore Lattice-Boltzmann lo riproduce con un errore L2 dello 0,27%. Valida prima contro la forma chiusa, poi lascia il risolutore libero sui flussi che non ne hanno una.

Il campo di vorticità di una scia di vortici di von Kármán — vortici rossi e blu alternati che si staccano da un cilindro — Il protagonista: il flusso attorno a un cilindro al di sopra di un numero di Reynolds critico stacca vortici alternati — la scia di vortici di von Kármán. Il campo di vorticità mostra il pattern caratteristico di vortici controrotanti (rossi e blu) che si staccano dal cilindro e vengono convetti a valle. Il risultato più riconoscibile della fluidodinamica, risolto su una griglia.

La scia non stacca sempre vortici — finché non lo fa.

A basso numero di Reynolds il flusso dietro il cilindro è una coppia stazionaria e simmetrica di bolle di ricircolo. Alza il numero di Reynolds oltre un valore critico e quella simmetria si rompe: la scia diventa instabile e inizia a staccare vortici a una frequenza periodica netta — una biforcazione di Hopf nel flusso stesso. Il programma esplora il numero di Reynolds e cattura la transizione, con l’ampiezza di oscillazione della scia che salta da sostanzialmente zero a un distacco pienamente sviluppato.

Ampiezza di oscillazione della scia in funzione del numero di Reynolds — quasi zero a basso Re, ampia sopra il valore critico — La transizione di distacco: sotto un numero di Reynolds critico (~47) la scia è stazionaria e l’ampiezza di oscillazione è sostanzialmente zero; sopra di esso la scia stacca vortici e l’ampiezza salta. Una biforcazione nel flusso, catturata esplorando un singolo numero adimensionale.

Il distacco lascia un’impronta sulle forze.

La scia di vortici non è solo bella — è misurabile. Sonda la scia e la sua velocità oscilla a una frequenza che, resa adimensionale, è il numero di Strouhal ≈ 0,2, il valore che ogni manuale di fluidodinamica cita. E il distacco scuote il cilindro: la forza di portanza oscilla alla frequenza di distacco mentre la resistenza ondeggia esattamente al doppio — la firma di von Kármán, e il motivo per cui ciminiere, cavi e impalcati di ponti devono essere progettati contro le vibrazioni indotte dai vortici.

Resistenza e portanza sul cilindro nel tempo — la portanza che oscilla alla frequenza di distacco, la resistenza al doppio — La firma delle forze: i vortici alternati spingono il cilindro lateralmente, quindi la portanza oscilla alla frequenza di distacco e la resistenza ondeggia esattamente al doppio (coefficiente di resistenza ≈ 0,94). Questa è la vibrazione indotta dai vortici — il motivo per cui le strutture snelle in un flusso sono progettate esattamente contro questo.

L’altro test canonico, superato anch’esso.

La cavità trascinata dal coperchio — un fluido in una scatola trascinato da un coperchio superiore in movimento — è l’altro caso di validazione CFD standard, con decenni di soluzioni di riferimento pubblicate. Il risolutore produce il vortice di ricircolo primario atteso, con il suo centro che cade esattamente dove lo mette la letteratura (Ghia et al.). Due benchmark canonici superati, da un unico risolutore autonomo scritto nel notebook.

Un flusso di cavità trascinata dal coperchio — un vortice di ricircolo primario con frecce di velocità che circolano nella scatola — La cavità trascinata dal coperchio: un coperchio superiore in movimento trascina il fluido in un unico vortice di ricircolo primario, con il suo centro che cade dove lo mettono le soluzioni di riferimento. L’altro benchmark classico del campo CFD — stesso risolutore, una scatola chiusa invece di una scia aperta.

Ogni numero riderivato all’approvazione finale.

Il notebook di V&V ricostruisce da zero il risolutore e rideriva ogni requisito rispetto a teoria e benchmark, stampando un quadro PASS/FAIL.

Risultato

Profilo di Poiseuille contro soluzione analitica: 0,27%
Ampiezza di distacco dei vortici / U₀: 0,87
Numero di Strouhal: 0,20
Coefficiente di resistenza / firma: 0,94; resistenza = 2× portanza
Centro del vortice primario di cavità: (0,61, 0,72)

Requisito

Profilo di Poiseuille contro soluzione analitica: R-01 < 2%
Ampiezza di distacco dei vortici / U₀: R-02
Numero di Strouhal: R-03 ≈ 0,15-0,2
Coefficiente di resistenza / firma: R-04
Centro del vortice primario di cavità: R-05

Un vero risolutore a campo, a fedeltà didattica.

Il deliverable è dieci notebook e il dossier — un campo PDE discretizzato, non una rete concentrata acausale, quindi non c’è blocco personalizzato né canvas. Il risolutore è Lattice-Boltzmann 2D incomprimibile, a basso Mach, con collisione a tempo di rilassamento singolo; la risoluzione della griglia è modesta, il cilindro usa un contorno a rimbalzo a gradini, i coefficienti di resistenza sono indicativi anziché validati da pubblicazione, e non c’è modello di turbolenza. Ciò che il programma dimostra è la capacità che mancava di più al portafoglio — risolvere un’equazione differenziale alle derivate parziali come campo spaziale su una griglia — con un risolutore autonomo (nessuna dipendenza CFD esterna) che riproduce la scia di vortici, il numero di Strouhal, la firma delle forze, e due benchmark canonici, ogni numero rieseguibile.