Visualizza wikitesto di Vorticità

{{F|fisica|agosto 2010}}
{{C|Alcuni passaggi non sono chiari.|fisica|gennaio 2008}}
In [[fluidodinamica]] si definisce '''vorticità''' la [[grandezza vettoriale]]:

:<math>\bar \omega = \vec\nabla \times \vec{v}</math>

ovverosia il [[rotore (matematica)|rotore]] della [[velocità]]; rappresenta inoltre la [[densità superficiale]] di [[circolazione (fluidodinamica)|circolazione]]:

:<math>\bar \omega = \frac {\partial \Gamma}{\partial r^2}</math>  
:dove:   
:*<math>\Gamma = \oint_{\partial S}     \vec{v}\cdot\vec{t} dS</math>
:*<math>\vec{t} </math>è il vettore tangente alla traiettoria

La vorticità è poi collegata alla [[velocità angolare|velocità di rotazione]] di un elemento di fluido e in particolare si dimostra studiandone la deformazione che:

:<math>\omega = 2 \Omega </math>

dove <math>\Omega</math> è la velocità istantanea di rotazione rigida dell'elemento di fluido.

In generale la dinamica della vorticità è descritta dall'[[equazione della vorticità]].

==Fluidi incomprimibili==
In un [[fluido incomprimibile]] la vorticità nasce solo in presenza di contorni solidi e deriva dalla ''condizione di aderenza'', ossia che il flusso non scorra sulla parete. La velocità relativa tra fluido e contorno solido è in altre parole considerata nulla e le particelle a diretto contatto con la parete vi aderiranno, mentre quelle appena più distanti avranno una certa velocità relativa non nulla.

Consideriamo una corrente uniforme ed una lastra piana di lunghezza ''Δy'' e di spessore trascurabile immersa nello stesso verso della corrente presente nel campo. Consideriamo una superficie bidimensionale infinitesima di fluido ''dx Δy'', dove con ''x'' si è indicato l'asse parallelo al flusso e con ''y'' l'asse normale alla lastra, a contatto con la faccia superiore della lastra stessa. La circolazione al contorno sarà perciò in base al [[teorema del rotore]]:

:<math>\operatorname d\Gamma=\frac {\operatorname d (x,y)}{\operatorname dt} \cdot \operatorname d(x,y) = \nabla \times \frac {\operatorname d (x,y)}{\operatorname dt} \cdot \operatorname dx \Delta y = \omega \operatorname dx \Delta y </math>

L'elemento avrà sulla faccia superiore una velocità lungo l'asse ''x'', mentre sulla faccia inferiore, a contatto con la lastra, dovrà soddisfare la [[condizione al contorno]] di aderenza, perciò l'equazione si riduce a:

:<math>\omega \, \operatorname dx \, \Delta y = - \frac {\operatorname dx}{\operatorname dt} \operatorname dx</math>

da cui si deduce che la quantità di vorticità totale presente ad una certa coordinata ''x'' della lastra sarà:

:<math>\omega \Delta y = - \frac {\operatorname dx}{\operatorname dt} </math>

la vorticità sulla faccia superiore della piastra sarà di conseguenza:

:<math>\omega_- = -\frac {1}{\Delta y} \frac {\operatorname dx}{\operatorname dt}</math>

Con analoghe considerazioni si arriva a dire che la vorticità totale presente sulla faccia inferiore sarà:

:<math>\omega_+ = \frac {1}{\Delta y} \frac {\operatorname dx}{\operatorname dt}</math>

Eseguendo la somma di questi due termini si può notare che la vorticità globale del flusso è nulla, come lo era prima di incontrare la lastra.

Nel caso di corpi con uno spessore, come per esempio un [[profilo alare]] simmetrico ad incidenza nulla, qualitativamente la dinamica è la stessa, ma ora il corpo sposterà le particelle di fluido durante il suo moto anche in direzione normale alla corrente indisturbata per soddisfare la ''condizione di impermeabilità'' (il corpo non è penetrabile dal fluido). Si formerà quindi un campo di moto intorno al corpo che sarà soluzione dell'[[equazione di Laplace]] che darà come risultato un [[flusso potenziale]] caratterizzato da un [[punto di ristagno]] anteriore ''r'' e da un punto di ristagno posteriore ''r + Δr'' ed in particolare sarà nota la velocità alla parete. Indicando con ''r'' l'ascissa curvilinea lungo il profilo, la vorticità complessiva presente in un punto della superficie del profilo sarà:

:<math>\omega = - \frac {1}{\Delta r}\frac {\operatorname d \bar r}{\operatorname dt} </math>

dove ''Δr'' è lo spessore locale dello strato di fluido presente a diretto contatto con la superficie del corpo interessato dalla presenza di vorticità.

La vorticità complessiva sulla faccia superiore sarà:

:<math>-{\int_{r}^{r+\Delta r}\frac {\operatorname d r}{\operatorname dt}_{-} \, \operatorname dr}</math>

mentre quella sulla faccia inferiore sarà:

:<math>{\int_{r}^{r+\Delta r} \frac {\operatorname dr}{\operatorname dt}_{+} \, \operatorname dr}</math>

La somma delle due vorticità è ancora nulla.

Questo non è casuale perché per i fluidi incomprimibili sussiste un teorema di conservazione della vorticità. Per i problemi tridimensionali questo può essere facilmente verificato partendo dal fatto che in questo caso vale comunque:

:<math>\nabla \cdot (\frac {\operatorname d \bar r}{\operatorname dt})=0</math>
 
Per i problemi bidimensionali la rigorosità di questo teorema è stata dimostrata rigorosamente solo di recente da J.C.Wu. e si è visto che la conservazione continua a valere se si considera oltre alla vorticità presente nel fluido anche quella presente all'interno dei corpi immersi dentro al fluido definendola come il doppio della loro velocità angolare, coerentemente anche al significato fisico della vorticità. Si può quindi scrivere in questo caso:

:<math>\frac {\operatorname d}{\operatorname dt}{\int_{R}\omega \, dR} = 0 \, </math>

dove il dominio ''R'' va scelto come spiegato considerando anche i corpi presenti dentro il campo di moto.

==Dinamica della vorticità in un fluido incomprimibile==
L'equazione della dinamica della vorticità può essere scritta facendo il rotore di tutti i termini dell'equazione di [[Legge di conservazione della quantità di moto|bilancio della quantità di moto]] scritta per un fluido incomprimibile nella forma:

:<math>{\partial \over\partial t} \frac {\operatorname d \bar r}{\operatorname dt}+\bar \omega \times \frac {\operatorname d \bar r}{\operatorname dt} =-\nabla\left({p\over \rho} + \frac 12 \left(\frac {\operatorname d \bar r}{\operatorname dt}\right)^2+ 2 \omega\right)+d_{\frac {\operatorname d \bar r}{\operatorname dt}}\nabla^2 \frac {\operatorname d \bar r}{\operatorname dt}</math>

dove <math>\bar f = - \nabla \omega </math>

Considerando le seguenti relazioni vettoriali:

:<math>\nabla^2 \frac {\operatorname d \bar r}{\operatorname dt}=\nabla\left(\nabla \cdot \frac {\operatorname d \bar r}{\operatorname dt}\right)-\nabla \times \left(\nabla \times\frac {\operatorname d \bar r}{\operatorname dt}\right)</math>

:<math>\nabla \times \left(\bar \omega \times \frac {\operatorname d \bar r}{\operatorname dt}\right)=\bar \omega \nabla \cdot \frac {\operatorname d \bar r}{\operatorname dt}-\frac {\operatorname d \bar r}{\operatorname dt} \nabla \cdot \bar \omega + \frac {\operatorname d \bar r}{\operatorname dt}\cdot \nabla \bar \omega-\bar \omega \cdot \nabla \frac {\operatorname d \bar r}{\operatorname dt}</math>

:<math>\nabla \cdot \bar \omega=\nabla \cdot \left(\nabla \times \frac {\operatorname d \bar r}{\operatorname dt}\right)=0</math>

:<math>\nabla \times \left(\nabla (\cdot)\right)=0</math>

:<math>\nabla \times {\partial \over\partial t} (\frac {\operatorname d \bar r}{\operatorname dt})={\partial \over\partial t} \left(\nabla \times \frac {\operatorname d \bar r}{\operatorname dt}\right)</math>

si ottiene, definendo <math>d_{\frac {\operatorname d \bar r}{\operatorname dt}}</math> la [[diffusività cinematica]]:

:<math>{\partial \over\partial t} (\nabla \times \frac {\operatorname d \bar r}{\operatorname dt})+\bar \omega \nabla \cdot \frac {\operatorname d \bar r}{\operatorname dt}-\frac {\operatorname d \bar r}{\operatorname dt} \nabla \cdot \bar \omega + \frac {\operatorname d \bar r}{\operatorname dt}\cdot \nabla \bar \omega-\bar \omega \cdot \nabla \frac {\operatorname d \bar r}{\operatorname dt}=\nabla \times \left(\nabla \left({p\over \rho} + \frac 12 \left(\frac {\operatorname d \bar r}{\operatorname dt}\right)^2 + 2 \omega\right)\right)+d_{\frac {\operatorname d \bar r}{\operatorname dt}} \nabla \times \left(\nabla \left(\nabla \cdot \frac {\operatorname d \bar r}{\operatorname dt}\right)-\nabla \times \left(\nabla \times \frac {\operatorname d \bar r}{\operatorname dt}\right)\right)</math>

e quindi, poiché valgono:

:<math>\nabla \frac {\operatorname d \bar r}{\operatorname dt}=0</math>

e:

:<math>\omega = \nabla \times \frac {\operatorname d \bar r}{\operatorname dt}</math>

si ottiene:

:<math>{\partial \over\partial t} \left(\nabla \times \frac {\operatorname d \bar r}{\operatorname dt}\right)+\frac {\operatorname d \bar r}{\operatorname dt}\cdot \nabla \bar \omega -\bar \omega \cdot \nabla \frac {\operatorname d \bar r}{\operatorname dt}=-d_{\frac {\operatorname d \bar r}{\operatorname dt}} \nabla \times(\nabla \times(\nabla \times \frac {\operatorname d \bar r}{\operatorname dt}))</math>

riordinando i termini:

:<math>{\partial \bar \omega \over\partial t}+\frac {\operatorname d \bar r}{\operatorname dt}\cdot \nabla \bar \omega=\bar \omega \cdot \nabla \frac {\operatorname d \bar r}{\operatorname dt}+d_{\frac {\operatorname d \bar r}{\operatorname dt}} \nabla^2 \bar \omega</math>

Nel caso in cui la vorticità sia nulla diremo che il flusso trattato sarà irrotazionale. Quest'ultimo concetto è piuttosto rilevante in quanto una corrente irrotazionale ammette l'esistenza di una funzione ''f'', definita a meno di una costante, tale per cui <math>\overrightarrow V=\nabla (f)</math>, tale funzione sarà la funzione potenziale del campo. Inoltre essendo per il [[teorema del rotore]]:

:<math>\iint_{S} (\nabla \times \frac {\operatorname d \bar r}{\operatorname dt}) \cdot \operatorname d \bar r^2 = \oint_C \frac {\operatorname d \bar r}{\operatorname dt} \cdot \operatorname d \bar r</math>

dove <math>\hat N</math> è il [[versore]] normale alla superficie ''S'' e <math>d \bar s</math> è il tratto infinitesimo della curva presa nel campo, la [[circuitazione]] di <math>\frac {\operatorname d \bar r}{\operatorname dt}</math> sarà nulla per qualsiasi circuito chiuso preso nel campo.

In un [[fluido incomprimibile]] la vorticità rappresenta la risultante delle azioni [[viscosità|viscose]] agenti su una particella fluida. Se questa grandezza è zero, l'equazione che ne regola il moto è identica a quella di un fluido non viscoso. Attenzione, questo non vuol dire che le azioni viscose siano nulle, ma solo la loro risultante. La vorticità non dice nulla sull'entità delle forze viscose, ma solo sulla loro rilevanza dinamica.

Nella teoria dello [[strato limite]] di Prandtl, per esempio, si divide il dominio in due parti, lo strato limite, dove la vorticità è decisamente non nulla, e la parte esterna, in cui il flusso è irrotazionale, e dove è quindi possibile definire una funzione potenziale della velocità;

:<math>\nabla \times \frac {\operatorname d \bar r}{\operatorname dt} = 0</math>

questo semplifica moltissimo il problema matematico che diventa a questo punto facilmente risolvibile con l'aiuto di un [[calcolatore]].

==Note==
<references/>

==Bibliografia==
*Clancy, L.J. (1975), ''Aerodynamics'', Pitman Publishing Limited, London ISBN 0-273-01120-0
* "''[https://web.archive.org/web/20060427140337/http://www.weather.com/glossary/v.html Weather Glossary]''"' The Weather Channel Interactive, Inc.. 2004.
* "''[http://www.tpub.com/content/aerographer/14010/css/14010_18.htm Vorticity]''". Integrated Publishing.
* {{Cita pubblicazione | nome=G. K. | cognome=Batchelor | linkautore=George Batchelor | titolo=An Introduction to Fluid Dynamics | editore=Cambridge University Press | anno=2000 |annooriginale=1967 | isbn=0-521-66396-2}}
* Ohkitani, K., "''Elementary Account Of Vorticity And Related Equations''". Cambridge University Press. January 30, 2005. ISBN 0-521-81984-9
* [[Alexandre Chorin|Chorin,  Alexandre J.]], "''Vorticity and Turbulence''". Applied Mathematical Sciences, Vol 103, Springer-Verlag. March 1, 1994. ISBN 0-387-94197-5
* [[Andrew Majda|Majda, Andrew J.]], Andrea L. Bertozzi, "''Vorticity and Incompressible Flow''". Cambridge University Press; 2002. ISBN 0-521-63948-4
* Tritton, D. J., "''Physical Fluid Dynamics''". Van Nostrand Reinhold, New York. 1977. ISBN 0-19-854493-6
* Arfken, G., "''Mathematical Methods for Physicists''", 3rd ed. Academic Press, Orlando, FL. 1985. ISBN 0-12-059820-5

==Voci correlate==
*[[Equazione della vorticità]]
* [[Circolazione (fluidodinamica)]]
*[[Vorticità potenziale]]
*[[Funzione di corrente]]
* [[Flusso potenziale]]
*[[Flusso potenziale incomprimibile]]
* [[Legge di Biot-Savart]]
* [[Teorema di Kutta-Žukovskij]]
*[[Equazione di Bernoulli]]
*[[Equazioni di Navier-Stokes]]
*[[Enstrofia]]

== Altri progetti ==
{{interprogetto}}

== Collegamenti esterni ==
* Weisstein, Eric W., "''[http://scienceworld.wolfram.com/physics/Vorticity.html Vorticity]''". Scienceworld.wolfram.com.
* Doswell III, Charles A.,  "''[http://www.cimms.ou.edu/~doswell/vorticity/vorticity_primer.html A Primer on Vorticity for Application in Supercells and Tornadoes]''". Cooperative Institute for Mesoscale Meteorological Studies, Norman, Oklahoma. 
* Cramer, M. S., "''Navier-Stokes Equations -- [https://web.archive.org/web/20040723231025/http://www.navier-stokes.net/nsvint.htm Vorticity Transport Theorems]: Introduction''". Foundations of Fluid Mechanics.
* Parker, Douglas, "''ENVI 2210 - Atmosphere and Ocean Dynamics, [https://web.archive.org/web/20040504110849/http://www.env.leeds.ac.uk/envi2210/lectures/lect9.html 9: Vorticity]''". School of the Environment, University of Leeds. September 2001. 
* [[James R. Graham|Graham, James R.]],  "''Astronomy 202: Astrophysical Gas Dynamics''". Astronomy Department, [[UC Berkeley]]. 
** "''[https://web.archive.org/web/20040702114551/http://astron.berkeley.edu/~jrg/ay202/node92.html The vorticity equation: incompressible and barotropic fluids]''". 
** "''[https://web.archive.org/web/20040702055351/http://astron.berkeley.edu/~jrg/ay202/node93.html Interpretation of the vorticity equation]''". 
** "''[https://web.archive.org/web/20040702130555/http://astron.berkeley.edu/~jrg/ay202/node94.html Kelvin's vorticity theorem for incompressible or barotropic flow]''". 
* "''[http://www.scd.ucar.edu/css/software/spherepack/ Spherepack 3.1] {{Webarchive|url=https://web.archive.org/web/20040622182810/http://www.scd.ucar.edu/css/software/spherepack/ |date=22 giugno 2004 }}''". (includes a collection of FORTRAN vorticity program) 
* "''{{collegamento interrotto|1=[http://132.206.43.151:5080/realtime/main_page.html Mesoscale Compressible Community (MC2)] |data=aprile 2018 |bot=InternetArchiveBot }} Real-Time Model Predictions''". (Potential vorticity analysis)

{{Portale|Aviazione|Fisica|Ingegneria}}

[[Categoria:Fluidodinamica]]