La formulazione funzione di corrente–vorticità semplifica i problemi di flusso laminare a basso numero di Reynolds sostituendo un’equazione del quarto ordine con due equazioni accoppiate del secondo ordine. Per un flusso stazionario, incomprimibile e piano di Stokes, si definisce la funzione di corrente in modo che la continuità sia soddisfatta automaticamente, quindi la si mette in relazione con la vorticità tramite (\nabla^2\psi=-\omega_z). Il problema governante diventa un’equazione di Poisson per (\psi) accoppiata a un’equazione di Laplace per (\omega_z), generalmente più facile da risolvere analiticamente e numericamente.
Concetto fondamentale: Nel flusso di Stokes, utilizzare (\nabla^2\psi=-\omega_z) e (\nabla^2\omega_z=0) invece di risolvere direttamente l’equazione biharmonica (\nabla^4\psi=0). Ciò riduce l’ordine differenziale, impone naturalmente la conservazione della massa e fornisce un quadro pratico per geometrie di operazioni unitarie che coinvolgono pareti, particelle, canali e mezzi porosi.
Perché il flusso laminare a basso numero di Reynolds è particolarmente adatto
Dominano le forze viscose
Il flusso laminare a basso numero di Reynolds si verifica quando il numero di Reynolds è molto inferiore a uno:
[ \mathrm{Re}=\frac{\rho U L}{\mu}\ll 1 ]
In questa condizione, i termini inerziali sono trascurabili rispetto alle forze viscose e di pressione. L’equazione della quantità di moto si riduce all’equazione di Stokes:
[ -\nabla p+\mu\nabla^2\mathbf{v}+\mathbf{f}=0 ]
dove (\mathbf{f}) rappresenta le forze di volume, quando pertinenti.
Le equazioni governanti diventano lineari
Poiché l’accelerazione convettiva non lineare è assente, velocità e pressione dipendono linearmente dal moto delle pareti, dai gradienti di pressione e dalle forze di volume. Questa linearità rende disponibile il principio di sovrapposizione per le soluzioni analitiche e migliora la robustezza dei solutori numerici.
Il flusso risponde senza ritardo inerziale
Nell’approssimazione ideale del flusso di Stokes, il campo di velocità si adegua immediatamente alle variazioni di pressione o al moto delle pareti. Questo comportamento pseudostazionario è utile nell’analisi di apparecchiature azionate lentamente, fluidi altamente viscosi e flussi su scala microscopica.
Come la funzione di corrente elimina un vincolo
La continuità è soddisfatta automaticamente
Per un flusso incomprimibile bidimensionale, si definisce la funzione di corrente (\psi) mediante
[ u=\frac{\partial\psi}{\partial y}, \qquad v=-\frac{\partial\psi}{\partial x} ]
L’equazione di continuità è
[ \frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0 ]
La sostituzione delle definizioni della funzione di corrente soddisfa identicamente questa equazione. Pertanto, non è necessario risolvere separatamente la continuità a ogni passaggio.
Le linee di corrente derivano direttamente da (\psi)
Le linee a (\psi) costante sono linee di corrente. La differenza tra i valori della funzione di corrente di due linee di corrente rappresenta la portata volumetrica per unità di profondità tra di esse, in base alla convenzione di segno adottata.
Ciò conferisce alla formulazione un valore sia matematico sia fisico: la variabile numerica descrive direttamente il modello di flusso e la velocità di trasporto.
La vorticità collega la rotazione alla curvatura della velocità
Per un flusso piano, la vorticità fuori dal piano è
[ \omega_z=\frac{\partial v}{\partial x} -\frac{\partial u}{\partial y} ]
Utilizzando le definizioni della funzione di corrente si ottiene
[ \omega_z=-\nabla^2\psi ]
o, in modo equivalente,
[ \nabla^2\psi=-\omega_z ]
Il segno cambia se si utilizza una convenzione alternativa per la funzione di corrente; pertanto, la convenzione deve rimanere coerente per tutta la derivazione e nel codice.
Come viene disaccoppiata l’equazione del quarto ordine
Formulazione diretta con la funzione di corrente
Calcolando il rotore dell’equazione della quantità di moto di Stokes si elimina la pressione. Per viscosità costante e in assenza di un contributo di forza di volume rotazionale, l’equazione risultante per la funzione di corrente è
[ \nabla^4\psi=0 ]
Questa è l’equazione biharmonica. In coordinate cilindriche o sferiche, l’operatore corrispondente viene modificato dalla geometria, spesso indicato schematicamente come
[ E^4\psi=0 ]
Formulazione accoppiata del secondo ordine
Si introduce la vorticità come campo indipendente. Il problema del quarto ordine diventa
[ \nabla^2\psi=-\omega_z ]
e
[ \nabla^2\omega_z=0 ]
La prima equazione è un’equazione di Poisson per la funzione di corrente. La seconda è un’equazione di Laplace per la vorticità.
Questa è la semplificazione centrale: due equazioni del secondo ordine sono generalmente più facili da discretizzare, risolvere e associare alle condizioni al contorno rispetto a un’unica equazione del quarto ordine.
Perché un ordine inferiore è vantaggioso
Le discretizzazioni del secondo ordine sono ampiamente supportate dai metodi alle differenze finite, ai volumi finiti e agli elementi finiti. Inoltre, richiedono informazioni al contorno in una forma più vicina al problema fisico, come la velocità della parete, l’impermeabilità o un comportamento prescritto della vorticità.
Un’equazione del quarto ordine richiede condizioni al contorno aggiuntive sulle derivate, che possono essere difficili da specificare accuratamente in corrispondenza di pareti solide, angoli, interfacce e superfici delle particelle.
Vantaggi numerici per le operazioni unitarie
La pressione viene eliminata dalla soluzione principale
La formulazione funzione di corrente–vorticità determina il campo di velocità senza risolvere la pressione come incognita principale. La pressione può essere ricostruita successivamente se sono richiesti la perdita di pressione, la forza o il consumo di potenza.
Ciò è utile nei problemi in cui inizialmente la struttura del flusso è più importante del campo di pressione dettagliato, come la circolazione di sospensioni o il ricircolo attorno alle particelle.
La formulazione è adatta alle geometrie strutturate
Per canali piani, interstizi, cavità, film e passaggi porosi idealizzati, le due equazioni possono essere discretizzate su griglie strutturate. I sistemi lineari risultanti possono essere risolti con metodi standard per le equazioni ellittiche.
Le fasi numeriche tipiche sono:
- Applicare le condizioni al contorno per (\psi) e per la velocità.
- Risolvere l’equazione della vorticità.
- Risolvere l’equazione di Poisson per (\psi).
- Ricavare (u) e (v) dalle derivate di (\psi).
- Iterare fino alla convergenza dei campi accoppiati.
Per il flusso puro di Stokes, le equazioni sono lineari. L’iterazione è principalmente associata all’accoppiamento, al trattamento della vorticità al contorno, alla geometria o alle estensioni che includono proprietà variabili e fisica aggiuntiva.
La conservazione della massa è integrata nel recupero della velocità
Le velocità ricavate da (\psi) soddisfano automaticamente l’incomprimibilità, a parte gli errori di discretizzazione. Ciò riduce il rischio di ottenere un campo di velocità che violi la continuità, una preoccupazione comune nelle formulazioni pressione–velocità.
Le condizioni al contorno diventano la principale difficoltà numerica
La condizione di non scorrimento prescrive entrambe le componenti della velocità su una parete solida stazionaria. Nella forma con funzione di corrente, l’impermeabilità determina comunemente un valore costante di (\psi) lungo la parete, mentre la velocità tangenziale fornisce informazioni sulle derivate.
La vorticità sulla parete spesso non viene prescritta direttamente. Deve essere dedotta dall’equazione della funzione di corrente e dalla condizione sulla velocità in prossimità della parete, utilizzando una discretizzazione appropriata del contorno.
Vantaggi analitici nei problemi di ingegneria chimica
La sovrapposizione diventa semplice
Poiché il flusso di Stokes è lineare, le soluzioni possono essere sommate. Un moto complesso può essere scomposto in componenti più semplici, come traslazione, rotazione, flusso indotto dalla pressione e moto delle pareti.
Ad esempio, il flusso attorno a una particella in movimento può essere costruito combinando un problema di flusso uniforme con un campo perturbato che soddisfi le condizioni al contorno della particella.
La simmetria riduce il dominio
Le equazioni della funzione di corrente e della vorticità preservano la simmetria geometrica. Se un canale, una disposizione di particelle o una cavità è simmetrica, può essere necessario analizzare solo una parte del dominio.
Ciò riduce la complessità algebrica e facilita la verifica dei casi limite.
La separazione delle variabili rimane utile
In geometrie cartesiane, cilindriche e sferiche semplici, si possono assumere forme come
[ \psi(x,y)=X(x)Y(y) ]
oppure utilizzare sviluppi di Fourier, funzioni proprie o multipoli. Le equazioni del secondo ordine producono equazioni differenziali ordinarie più gestibili rispetto alla manipolazione diretta dell’equazione biharmonica della funzione di corrente.
I risultati classici delle operazioni unitarie vengono recuperati naturalmente
Diversi risultati familiari possono essere interpretati attraverso lo stesso quadro:
- Il flusso completamente sviluppato in un canale presenta un profilo semplice della funzione di corrente e una vorticità determinata dal gradiente di pressione.
- Il flusso in interstizi stretti è governato principalmente dalla resistenza viscosa e può essere approssimato utilizzando le ipotesi di lubrificazione.
- Il flusso a basso numero di Reynolds attorno alle particelle fornisce i campi di resistenza e perturbazione utilizzati nell’analisi della sedimentazione e delle sospensioni.
- Le formulazioni in coordinate cilindriche descrivono flussi idealizzati in capillari, condotti anulari e canali porosi.
La forma delle coordinate del laplaciano e la relazione della vorticità devono essere modificate opportunamente per le geometrie cilindriche o sferiche; la formula cartesiana non può essere riutilizzata semplicemente senza modifica.
Applicazioni nelle operazioni unitarie
Lavorazione dei polimeri e fluidi altamente viscosi
I fusi polimerici e i liquidi concentrati presentano spesso viscosità elevate e numeri di Reynolds caratteristici bassi. La formulazione aiuta a calcolare profili di velocità, sforzo alla parete, ricircolo e comportamento del tempo di residenza in passaggi azionati lentamente.
Per i fluidi non newtoniani, tuttavia, la semplice equazione armonica della vorticità generalmente non si applica più nella sua forma originaria. La viscosità dipende dal tasso di deformazione, quindi le equazioni della quantità di moto e della vorticità diventano accoppiate attraverso il modello costitutivo.
Microcanali e passaggi stretti
Le piccole scale dimensionali riducono il numero di Reynolds, rendendo comuni le ipotesi di flusso a basso numero di Reynolds nei sistemi microfluidici e di trasporto su scala microscopica. I metodi basati sulla funzione di corrente possono risolvere il flusso attorno a ostacoli, restringimenti, espansioni e contorni mobili.
Sono particolarmente utili quando l’obiettivo principale è il modello locale del flusso, la struttura della velocità indipendente dalla pressione o il trasporto di particelle.
Sospensioni e trasporto di particelle
Il flusso attorno a colloidi, gocce, bolle e particelle sospese di maggiori dimensioni determina la resistenza idrodinamica e la migrazione. La soluzione per (\psi) e (\omega) fornisce il campo perturbato necessario per stimare le forze e le interazioni tra particelle.
Il metodo rende inoltre evidente la reversibilità del flusso: invertendo il moto imposto alla particella o alla parete, si invertono il campo di velocità e le linee di corrente nel limite ideale del flusso lineare di Stokes.
Idealizzazioni dei mezzi porosi
Alla scala dei pori, le soluzioni del flusso a basso numero di Reynolds aiutano a spiegare la resistenza viscosa e il trasporto indotto dalla pressione attraverso passaggi semplificati. Possono inoltre fungere da elementi costitutivi locali per modelli di permeabilità su scala maggiore.
Un mezzo poroso reale richiede generalmente un metodo numerico che risolva la geometria o un modello macroscopico come il flusso di Darcy. La formulazione con funzione di corrente è direttamente applicabile soprattutto a geometrie dei pori bidimensionali o opportunamente ridotte.
Comprendere i compromessi
Le condizioni al contorno della vorticità non sono sempre semplici
L’introduzione della vorticità riduce l’ordine differenziale, ma aggiunge un campo i cui valori al contorno possono essere indiretti. Approssimazioni inadeguate della vorticità sulla parete possono generare sforzi di taglio imprecisi o comportamenti spuri in prossimità della parete.
Questo compromesso è particolarmente importante quando la grandezza di interesse è lo sforzo alla parete, la resistenza o la perdita di pressione, anziché soltanto il campo di velocità globale.
Angoli e singolarità richiedono attenzione
Angoli acuti, espansioni improvvise e contatti ravvicinati tra particelle e pareti possono generare forti gradienti di velocità o singolarità matematiche. Una griglia grossolana o un trattamento del contorno di basso ordine può produrre una vorticità locale fuorviante anche quando il flusso globale appare ragionevole.
Sono necessari un raffinamento locale della griglia e un confronto accurato con soluzioni analitiche limite.
La pressione non è più immediatamente disponibile
La pressione scompare quando l’equazione della quantità di moto viene trasformata nella forma in termini di vorticità. Se sono richieste la perdita di pressione, la potenza di pompaggio o il bilancio delle forze, la pressione deve essere ricavata dalle equazioni originali di Stokes o da una procedura equivalente di post-elaborazione.
La formulazione semplifica quindi la soluzione della velocità, ma non elimina la pressione dall’analisi ingegneristica complessiva.
Le ipotesi possono non essere valide
Le equazioni di base assumono incomprimibilità, comportamento newtoniano, viscosità costante e inerzia trascurabile. Possono richiedere modifiche per:
- Fusi polimerici non newtoniani.
- Fluidi a viscosità variabile.
- Flussi comprimibili o fortemente dipendenti dalla temperatura.
- Geometrie tridimensionali.
- Effetti inerziali quando il numero di Reynolds non è sufficientemente basso.
- Forze di volume o effetti multifisici che contribuiscono alla generazione di vorticità.
I metodi basati sulla funzione di corrente sono principalmente bidimensionali
Una singola funzione di corrente scalare è particolarmente conveniente per il flusso piano incomprimibile. I flussi completamente tridimensionali richiedono generalmente potenziali vettoriali, formulazioni velocità–pressione o altri metodi specializzati.
Il metodo è quindi efficace per geometrie ridotte, ma non dovrebbe essere applicato a operazioni unitarie intrinsecamente tridimensionali.
Scegliere l’approccio giusto in base al proprio obiettivo
La formulazione è più efficace quando la geometria e le ipotesi costitutive supportano un modello bidimensionale, incomprimibile e a basso numero di Reynolds.
- Se l’obiettivo principale sono i profili di velocità o le linee di corrente: risolvere (\nabla^2\psi=-\omega_z) insieme a (\nabla^2\omega_z=0), poiché la continuità viene imposta automaticamente e il modello del flusso deriva direttamente da (\psi).
- Se l’obiettivo principale è la robustezza numerica: utilizzare discretizzazioni accoppiate del secondo ordine, trattando con attenzione la vorticità sulle pareti e le regioni angolari.
- Se l’obiettivo principale è la derivazione analitica: sfruttare linearità, simmetria, separazione delle variabili e sovrapposizione prima di tentare una soluzione biharmonica generale.
- Se l’obiettivo principale è la perdita di pressione o la potenza di pompaggio: utilizzare la soluzione funzione di corrente–vorticità per la velocità, quindi ricostruire la pressione dall’equazione della quantità di moto di Stokes.
- Se l’obiettivo principale è la lavorazione dei polimeri o la reologia delle sospensioni: mantenere la formulazione come punto di partenza concettuale, ma includere l’appropriata relazione costitutiva non newtoniana prima di applicarla quantitativamente.
- Se l’obiettivo principale è un’operazione unitaria tridimensionale: preferire un metodo tridimensionale velocità–pressione o basato su potenziali vettoriali, a meno che la geometria consenta una riduzione bidimensionale giustificata.
La relazione funzione di corrente–vorticità trasforma l’analisi dei flussi a basso numero di Reynolds in un problema di ordine inferiore e coerente con la conservazione, preservando al contempo l’intuizione fisica necessaria per progettare e interpretare le operazioni unitarie dell’ingegneria chimica.
Tabella riepilogativa:
| Aspetto chiave | Indicazione pratica |
|---|---|
| Equazioni fondamentali | Risolvere (\nabla^2\psi=-\omega_z) e (\nabla^2\omega_z=0) invece dell’equazione biharmonica del quarto ordine. |
| Vantaggio principale | Due equazioni accoppiate del secondo ordine sono più facili da discretizzare e da associare alle condizioni al contorno. |
| Conservazione | La definizione della velocità tramite la funzione di corrente soddisfa automaticamente la continuità incomprimibile. |
| Significato fisico | Le curve della funzione di corrente mostrano le linee di corrente e la portata; la vorticità rappresenta la rotazione locale del fluido. |
| Usi numerici | Adatta a canali, cavità, interstizi stretti, particelle, microflussi e passaggi porosi idealizzati. |
| Strumenti analitici | La linearità consente la sovrapposizione, la riduzione per simmetria, la separazione delle variabili e le soluzioni classiche. |
| Principali difficoltà | Il trattamento della vorticità sulle pareti, gli angoli acuti, le singolarità, il recupero della pressione e il raffinamento della griglia richiedono attenzione. |
| Limiti di applicabilità | È più adatta a flussi bidimensionali, incomprimibili, newtoniani, a viscosità costante e a basso numero di Reynolds. |
| Applicazioni estese | I flussi non newtoniani, a viscosità variabile, tridimensionali o con effetti inerziali richiedono formulazioni modificate. |
| Risultati ingegneristici | Utilizzare la soluzione della velocità per valutare profili, linee di corrente, sforzi di taglio, resistenza, comportamento del tempo di residenza e trasporto. |
Avete bisogno di apparecchiature affidabili per la didattica, la ricerca o lo sviluppo di processo? Contattate LABPARK per discutere le vostre esigenze. LABPARK fornisce impianti pilota per operazioni unitarie didattiche e professionali destinati all’ingegneria chimica, ai bioprocessi e alle biotecnologie, nonché alle applicazioni ambientali e al trattamento delle acque. I nostri sistemi aiutano università, istituti di ricerca e aziende a dimostrare i principi dei flussi a basso numero di Reynolds, studiare il trasporto indotto dalla pressione, analizzare il comportamento di particelle e sospensioni e collegare i modelli teorici ai dati operativi pratici. Possiamo aiutarvi a selezionare una configurazione di impianto pilota in linea con gli obiettivi del vostro corso, la capacità del laboratorio, le esigenze di ricerca e gli obiettivi di sviluppo del processo.