ARTE TECNOLOGICA

L'equazione della vorticità

20.06.2014 21:49

L'equazione della vorticità descrive l'evoluzione nel tempo della vorticità, ed è quindi utile nella comprensione dei moti rotatori nei fluidi. Nel caso di fluidi barotropici non viscosi essa stabilisce che il vettore vorticità è solidale con la materia, si sposta con essa e la sua evoluzione è tale da garantire la conservazione del momento angolare.

L'equazione, nella sua variante per fluidi incompressibili non viscosi, è stata ricavata da Helmholtz ^[1].

La forma generale è data da:

$\frac{\operatorname d}{\operatorname dt} \frac{\bar \omega}{\rho} - \left(\frac{\bar \omega}{\rho} \cdot \nabla\right) \frac{\operatorname d \bar r}{\operatorname dt} = \frac{1}{\rho} \nabla \times \frac{\operatorname d^2 \bar r}{\operatorname dt^2}$

dove d/dt rappresenta la derivata totale rispetto al tempo, ω è la vorticità, $\bar r$ è la posizione. Quindi la variazione temporale della vorticità è determinata dal rotore delle forze in gioco.

Nel caso di fluidi baroclini non viscosi l'equazione diventa:

$\frac{\operatorname d}{\operatorname dt} \frac{\bar \omega}{\rho} - \left(\frac{\bar \omega}{\rho} \cdot \nabla\right) \frac{\operatorname d \bar r}{\operatorname dt} = \frac{\bar B}{\rho}$

dove il vettore B è la baroclinicità ^[2] dato da:

$\bar B = \frac{1}{\rho^2} \nabla \rho \times \nabla p$

Nel caso di fluidi barotropici non viscosi essa si riduce a:

$\frac{\operatorname d}{\operatorname dt} \frac{\bar \omega}{\rho} - \left(\frac{\bar \omega}{\rho} \cdot \nabla\right) \frac{\operatorname d \bar r}{\operatorname dt} = 0$

Indice

Derivazione formale
Significato fisico
- Conservazione del momento angolare
- Vorticità solidale con la materia
Equazione della vorticità per moti sinottici
Note
Bibliografia
Voci correlate

Derivazione formale
Modifica

Si parte calcolando il rotore della derivata totale della velocità:

$\frac{\operatorname d^2 \bar r}{\operatorname dt^2} = \frac{\partial}{\partial t} \frac{\operatorname d \bar r}{\operatorname dt} + \left(\frac{\operatorname d \bar r}{\operatorname dt} \cdot \nabla \right) \frac{\operatorname d \bar r}{\operatorname dt}$

La seguente identità, dovuta a Lagrange ^[3], agevola i calcoli

$\frac{\operatorname d^2 \bar r}{\operatorname dt^2} = \frac{\partial}{\partial t} \frac{\operatorname d \bar r}{\operatorname dt} + \bar \omega \times \frac{\operatorname d \bar r}{\operatorname dt} + \frac{1}{2} \nabla \left(\frac{\operatorname d \bar r}{\operatorname dt}\right)^ 2$

calcolando il rotore di questa eguaglianza, e utilizzando il fatto che

$\nabla \cdot \bar \omega = 0$

si ottiene:

$\nabla \times \frac{\operatorname d^2 \bar r}{\operatorname dt^2} = \frac{\operatorname d \bar \omega}{\operatorname dt} - \left(\bar \omega \cdot \nabla\right) \frac{\operatorname d \bar r}{\operatorname dt} + \bar \omega \left(\nabla \cdot \frac{\operatorname d \bar r}{\operatorname dt} \right)$

utilizzando l'equazione di continuità della massa

$\frac{1}{\rho} \frac{\operatorname d \rho}{\operatorname dt} = - \nabla \cdot \frac{\operatorname d \bar r}{\operatorname dt}$

si ottiene

$\nabla \times \frac{\operatorname d^2 \bar r}{\operatorname dt^2} = \frac{\operatorname d \bar \omega}{\operatorname dt} - \left(\bar \omega \cdot \nabla\right) \frac{\operatorname d \bar r}{\operatorname dt} - \frac{\bar \omega}{\rho} \frac{\operatorname d \rho}{\operatorname dt}$

dividendo per ρ e utilizzando il fatto che

$\frac{\operatorname d}{\operatorname dt} \frac{\bar \omega}{\rho} = \frac{1}{\rho} \frac{\operatorname d \bar \omega}{\operatorname dt} - \frac{\bar \omega}{\rho^2} \frac{\operatorname d\rho}{\operatorname dt}$

si ottiene l'asserto. Nel caso di moti baroclini non viscosi si ha, dall'equazione del momento di Navier-Stokes:

$\frac{\operatorname d^2 \bar r}{\operatorname dt^2} = - \frac{1}{\rho} \nabla p - \nabla \phi$

dove p è la pressione, Φ è il geopotenziale. Dato che il rotore di un gradiente è nullo la parte della forza dovuta a ϕ non contribuisce a modificare la vorticità. Calcolando il rotore del termine dovuto alla pressione si ottiene:

$- \nabla \times \frac{\nabla p}{\rho} = \frac{1}{\rho^2} \nabla \rho \times \nabla p$

quindi otteniamo la versione dell'equazione per fluidi baroclini non viscosi.

Nei fluidi barotropici la densità è funzione della sola pressione, quindi i gradienti delle due grandezze sono paralleli. Pertanto la baroclinicità si annulla, e si ottiene la versione dell'equazione per fluidi barotropici non viscosi.

Significato fisico
Modifica

Il significato fisico dell'equazione della vorticità può essere meglio compreso considerando il caso semplificato di un fluido incompressibile non viscoso, non sottoposto a forze esterne. In questo caso il moto è inerziale, il momento angolare si conserva, e il vettore vorticità è solidale con la materia, cioè, se in un dato istante la vorticità è allineata a una linea materiale di lunghezza infinitesima $\delta$ , essa rimane sempre così allineata.

Conservazione del momento angolare

Con un semplice esempio si verifica che l'equazione della vorticità implica la conservazione del momento angolare del fluido, in assenza di forze con rotore non nullo quali la viscosità.

Si consideri un fluido ideale, con densità costante e privo di viscosità, in moto in un tubo ideale perfettamente liscio posto lungo l'asse x. Il moto è di tipo stazionario: il fluido si sposta con velocità u lungo l'asse x, e ruota su se stesso sul piano yz. Per comodità consideriamo che, per ogni sezione verticale di spessore infinitesimo dx, il moto rotatorio avvenga con velocità angolare ω uguale su tutta la sezione. Quindi ogni sezione verticale si sposta verso destra con velocità u(x), e ruota su se stessa con velocità angolare ω(x).

A un certo punto il tubo si restringe, la superficie della sezione verticale passa da $\Sigma_A$ a $\Sigma_B$ . È facile prevedere come varia il moto: u si modifica in modo da garantire la costanza della portata volumetrica, quindi

$u_A \Sigma_A = u_B \Sigma_B$

La velocità angolare ω varia in modo da conservare il momento angolare. Il momento angolare di ogni sezione è dato dal prodotto tra la velocità angolare e il momento di inerzia:

$L = I\bar \omega = \frac{1}{2} m r^2$

quindi, tenendo conto del fatto che la superficie è proporzionale al quadrato del raggio si ottiene:

$\omega_A \Sigma_A = \omega_B \Sigma_B$

Verifichiamo che l'equazione della vorticità porta esattamente allo stesso risultato. In ogni sezione del tubo la vorticità è diretta lungo l'asse x, ed è data da:

$\omega_x = \frac{\partial v}{\partial z} - \frac{\partial w}{\partial y}$

dove v e w sono rispettivamente la componente y e z della velocità. L'equazione della vorticità si riduce a:

$\frac{\operatorname d \omega_x}{\operatorname dt} = \frac{\partial \omega_x}{\partial t} + u \frac{\partial \omega_x}{\partial x} = \omega_x \frac{\partial u}{\partial x}$

Siccome il moto è stazionario, per ogni x si ha $\frac{\partial \omega_x}{\partial t} = 0$ . Quindi si ricava subito:

$\omega_x \propto u$

Utilizzando il fatto che $u_A \Sigma_A = u_B \Sigma_B$ si ottiene

$\omega_{xA} \Sigma_A = \omega_{xB} \Sigma_B$

\operatorname dato che la velocità angolare è uniforme su tutta la sezione, si può facilmente verificare che:

$\omega_x = 2 \bar \omega$

Quindi si ottiene la relazione ottenuta partendo dalla conservazione del momento angolare.

Si può proseguire il ragionamento notando che, siccome il fluido è incomprimibile, la superficie della sezione e il suo spessore sono inversamente proporzionali. Si arriva dunque alla relazione:

$\frac{\omega_{xA}}{\delta_{xA}} = \frac{\omega_{xB}}{\delta_{xB}}$

dove $\delta_x$ è lo spessore della sezione. La grandezza $\frac{\omega_{x}}{\delta_x}$ rappresenta la vorticità potenziale della sezione, che dunque si mantiene costante come conseguenza della conservazione del momento angolare.

Per approfondire, vedi vorticità potenziale.

Vorticità solidale con la materia

Si può verificare che la variazione della distanza $\delta$ tra due elementi di fluido contigui A e B è data da:

$\frac{\operatorname d}{\operatorname dt} \delta = \left(\delta \cdot \nabla\right) \frac{\operatorname d \bar r}{\operatorname dt}$

infatti, considerando per semplicità la sola componente x, abbiamo:

$\frac{\operatorname d}{\operatorname dt} \delta_x = \frac{\operatorname d \bar r}{\operatorname dt}'_A - \frac{\operatorname d \bar r}{\operatorname dt}'_B = \delta_x \frac{\partial}{\partial x} \frac{\operatorname d \bar x}{\operatorname dt}$

dove $\frac{\operatorname d \bar x}{\operatorname dt}$ è la componente x della velocità.

Quindi la legge con cui varia la vorticità in fluidi incompressibili, non viscosi e non sottoposti a forze esterne, è esattamentela stessa con cui variano le distanze tra elementi contigui di fluido. Ne risulta che se in un dato istante la vorticità è allineata alla linea materiale $\delta$ , essa rimane sempre così allineata.

Equazione della vorticità per moti sinottici
Modifica

Il concetto di vorticità ha grande importanza in meteorologia e in climatologia dato che è alla base della dinamica dei fluidi rotanti. Dalla conservazione di grandezze ad essa collegate, quali la vorticità potenziale, si deducono comportamenti fondamentali come le onde di Rossby.

Il modo più rapido per ricavare l'equazione della vorticità per moti su scala sinottica è calcolare il rotore delle equazioni del moto approssimato ^[4]:

$\frac{\operatorname d^2x}{\operatorname dt^2} - f \frac{\operatorname dy}{\operatorname dt} = - \frac{1}{\rho}\frac{\partial p}{\partial x}$

$\frac{\operatorname d^2y}{\operatorname dt^2} + f \frac{\operatorname dz}{\operatorname dt} = - \frac{1}{\rho}\frac{\partial p}{\partial y}$

dove compaione le componenti zonale (cioè rivolta verso Est) e meridionale (cioè rivolta verso Nord) della velocità, f è lavorticità planetaria, cioè la vorticità intrinseca del moto rotatorio della superficie terrestre. Naturalmente, dato che nei moti sinottici si considerano solo le componenti orizzontali della velocità, la vorticità ha la sola componente verticale non nulla. Il risultato è ^[5]:

$\frac{\operatorname d}{\operatorname dt} (\bar \omega + f) = - (\bar \omega + f) \left(\frac{\partial}{\partial x} \frac{\operatorname dx}{\operatorname dt}+ \frac{\partial}{\partial y} \frac{\operatorname dy}{\operatorname dt}\right) + \frac{1}{\rho^2} \left(\frac{\partial \rho}{\partial x}\frac{\partial p}{\partial y} - \frac{\partial \rho}{\partial y}\frac{\partial p}{\partial x}\right)$

dove ω è la vorticità relativa, cioè la vorticità del fluido rispetto alla superficie terrestre. L'ultimo termine rappresenta la baroclinicità. Pertanto nell'approssimazione barotropica l'equazione si riduce a:

Note
Modifica

^ Feynman, vol 2, 40-10
^ Gill, p238
^ Gill, p227
^ Holton, p40
^ Holton, p101

Bibliografia
Modifica

Adrian Gill, Atmosphere-Ocean dynamics, ISBN 0-12-283522-0
James R Holton, An introduction to dynamic meteorology, ISBN 978-0-12-354015-7, 4th edition
Richard Feynman, Lectures on physics, ISBN 0-8053-9049-9, The definitive edition

Voci correlate
Modifica

Indietro

Contatti

Italydee via L. Ori 1/B
Civitanova Marche
62012 (MC) 0733 70336 italydee@gmail.com cell. 334 6733765

L'equazione della vorticità

Indice

Derivazione formaleModifica

Significato fisicoModifica

Conservazione del momento angolare

Vorticità solidale con la materia

Equazione della vorticità per moti sinotticiModifica

NoteModifica

BibliografiaModifica

Voci correlateModifica