Equazioni di Navier-Stokes

Da Wikipedia, l'enciclopedia libera.
Se riscontri problemi nella visualizzazione dei caratteri, clicca qui.
MECCANICA CLASSICA
Meccanica del continuo

Le equazioni di Navier-Stokes sono un sistema di equazioni differenziali alle derivate parziali che descrivono il comportamento di un fluido dal punto di vista macroscopico. L'ipotesi di base è che il fluido possa essere modellato come un continuo deformabile. Esse presuppongono perciò la continuità del fluido in esame, ovverosia il sistema perde di validità nello studio di un gas rarefatto.

Le equazioni debbono il loro nome a Claude-Louis Navier e a George Gabriel Stokes che le formalizzarono e la loro soluzione analitica generale rappresenta attualmente uno dei problemi irrisolti della matematica moderna (i cosiddetti 7 problemi per il millennio) per il quale vale il premio Clay; soluzioni analitiche particolari si hanno in casi estremamente semplificati mentre soluzioni approssimate si ottengono tipicamente ricorrendo a metodi propri dell'analisi numerica e all'uso congiunto del calcolatore.

Queste equazioni rappresentano l'approssimazione di Chapman-Enskog del prim'ordine delle equazioni di bilancio canoniche (sono quindi più generali delle equazioni di Eulero, che costituiscono l'approssimazione precedente, e più particolari delle equazioni di Burnett che costituiscono quella successiva). Le relazioni costitutive che contengono sono anch'esse lineari: la legge di Newton e la legge di Fourier.

Il modello matematico[modifica | modifica wikitesto]

La maggiore efficienza predittiva di tali equazioni rispetto a quello di Eulero viene pagata in termini di difficoltà di soluzione. Nel caso generale coinvolgono infatti cinque equazioni scalari differenziali alle derivate parziali e 20 variabili. Il bilancio tra equazioni e incognite avviene (come vedremo più avanti) con la definizione delle proprietà del fluido considerato, delle eventuali forze di campo in gioco e con considerazioni matematiche. Inoltre, a causa della loro non linearità, le equazioni di Navier-Stokes non ammettono quasi mai una soluzione analitica (ovvero una soluzione esatta), ma esclusivamente numerica (una soluzione approssimata con un metodo numerico).

Le equazioni di Navier-Stokes sono in grado di descrivere completamente qualsiasi flusso fluido, anche turbolento. In particolare per un flusso turbolento, dove cioè le traiettorie delle particelle di flusso non sono più costanti nel tempo, un approccio numerico di calcolo è chiamato generalmente simulazione numerica diretta. A causa del fatto che le risorse di calcolo necessarie alla loro risoluzione crescono con il numero di Reynolds (quasi con Re³) e che tale numero può avere valori dell'ordine di 106-109, tale approccio resta tecnicamente impossibile. In alternativa alla simulazione numerica è possibile adottare sistemi meno onerosi quali la formulazione LES o le equazioni mediate.

Le equazioni vengono completate dalle condizioni al contorno e dalle condizioni iniziali (condizioni imposte all'inizio temporale del fenomeno da studiare). Possono inoltre essere integrate dall'equazione di stato dei gas perfetti e dalle equazioni di conservazione delle singole specie gassose nel caso di una miscela di gas.

La soluzione delle equazioni fornisce il campo delle velocità del fluido. Da questo sarà poi possibile risalire a tutte le altre grandezze che caratterizzano il flusso.

Ipotesi del modello matematico[modifica | modifica wikitesto]

Il modello matematico che permette l'analisi della dinamica dei continui deformabili si basa sulle seguenti caratteristiche:

  • fluido continuo;
  • fluido chimicamente omogeneo e non reagente;
  • fluido privo di cariche elettriche.

Ipotesi di fluido continuo[modifica | modifica wikitesto]

Viene trascurata la natura discontinua della materia, in questo modo sarà possibile far tendere a zero un volume di fluido, senza che questo possa restare privo di materia.

Un parametro fondamentale che caratterizza il mezzo dal punto di vista della continuità è il numero di Knudsen, definito come il rapporto tra il cammino libero medio di una particella costituente il fluido e una lunghezza caratteristica del flusso:

\mathrm{Kn} = \frac {l}{L}

Se il numero di Knudsen è molto minore di uno, allora è possibile considerare il fluido continuo. Altrimenti sarà necessario studiare il comportamento del gas unicamente su base statistica, mediante la teoria cinetica dei gas, la quale analizza statisticamente la distribuzione delle velocità molecolari e da questa ricava tutte le proprietà del gas.

Ipotesi di fluido chimicamente omogeneo e non reagente[modifica | modifica wikitesto]

Verranno trascurate le perturbazioni dovute alla non omogeneità del flusso ed alle reazioni chimiche. Ciò non sarà del tutto possibile per flussi reagenti quali quelli all'interno di una camera di combustione ad esempio.

Ipotesi di fluido privo di cariche elettriche[modifica | modifica wikitesto]

Verranno trascurate le perturbazioni dovute al campo elettromagnetico. L'interazione di flussi con campi elettromagnetici è studiata dalla magnetofluidodinamica.

Derivazione delle equazioni di Navier-Stokes[modifica | modifica wikitesto]

Le equazioni di Navier-Stokes sono la formalizzazione matematica di tre principi fisici ai quali i fluidi rispondono, imposta la condizione di continuo deformabile:

Per questo motivo sono spesso nominate anche equazioni di bilancio.

Nei successivi paragrafi indicheremo sempre il vettore velocità del fluido con la notazione  \vec u , mentre p e ρ indicheranno rispettivamente la pressione statica e la densità del fluido stesso. Il simbolo  \vec a rappresenterà il vettore delle accelerazioni di campo.

Descrizione del moto lagrangiana ed euleriana[modifica | modifica wikitesto]

È possibile descrivere temporalmente il moto di un fluido attraverso due punti di vista.

Il primo, detto lagrangiano o materiale, segue la traiettoria di ogni particella di fluido, identificata solitamente dalle sue coordinate iniziali, analizzando le variazioni delle sue proprietà fisiche (come ad esempio la densità o la temperatura). Indicando le coordinate iniziali come (a; b; c; 0), le coordinate della particella saranno variabili dipendenti:

\begin{cases} x = x (a; b; c; t) \\ y = y (a; b; c; t) \\ z = z (a; b; c; t) \end{cases}

così come le altre proprietà specifiche.

Il punto di vista euleriano, viceversa, osserva le variazioni delle proprietà fisiche per ogni data posizione spaziale (x; y; z). Le coordinate spaziali (assieme alla variabile temporale) saranno perciò variabili indipendenti. Le variabili dipendenti sono perciò funzione di quelle spaziali e temporali. Ad esempio, per la velocità:

\vec u = \vec u (x; y; z; t).

Teorema del trasporto di Reynolds[modifica | modifica wikitesto]

Exquisite-kfind.png Per approfondire, vedi teorema del trasporto di Reynolds e teorema della divergenza.

Per comodità di trattazione riportiamo il teorema del trasporto di Reynolds, che, per una proprietà:

A = A (x;y;z;t)

contenuta in un volume arbitrario V, che si muova con il fluido, ed abbia superficie S, è indicato come:

 \frac{d}{dt} \int_{V(t)} A \, dV = \int_{V(t)} \frac{\partial A}{\partial t} \, dV + \int_{S(t)} A \vec u \cdot \hat n \, dS .

Ricordando il teorema della divergenza è possibile esprimere la precedente anche come:

 \frac{d}{dt} \int_{V(t)} A \, dV = \int_{V(t)} \left( \frac{\partial A}{\partial t} + \nabla \cdot \left( A \vec u \right) \right) \, dV

e ricordando che:

 \nabla \cdot \left( A \vec u \right) = \left( \vec u \cdot \nabla \right) A + A \nabla \cdot \vec u

nonché la definizione di derivata totale, è possibile esprimere il teorema in una forma molto utile:

 \frac{d}{dt} \int_{V(t)} A \, dV = \int_{V(t)} \left( \frac{DA}{Dt} + A \nabla \cdot \vec u \right) \, dV.

Equazione di continuità[modifica | modifica wikitesto]

Exquisite-kfind.png Per approfondire, vedi legge della conservazione della massa (fisica).

Punto di vista lagrangiano[modifica | modifica wikitesto]

Il principio di conservazione della massa, nel caso di moto di un fluido, può essere espresso dal punto di vista lagrangiano affermando che:

« resta invariata nel tempo la massa contenuta in un volume (deformabile) che si muove con il fluido. »

In questo caso dunque, in termini matematici:

 \frac{dm}{dt} = \frac{d}{dt} \int_{V(t)} \rho \, dV = 0 \, .

Applicando il teorema del trasporto di Reynolds alla densità ρ (massa per unità di volume), otteniamo l'equazione di continuità in forma di divergenza:

 \frac {\partial \rho} {\partial t} + \nabla \cdot \left( \rho \vec u \right) = 0

che può essere riportata anche in forma indiciale:

 \frac {\partial \rho} {\partial t} + \frac{\partial \left( \rho u_i \right)}{\partial x_i} = 0

che può essere riportata anche in forma estesa:

 \frac {\partial \rho} {\partial t} + \frac{\partial \left( \rho u \right)}{\partial x} + \frac{\partial \left( \rho v \right)}{\partial y} + \frac{\partial \left( \rho w \right)}{\partial z} = 0

o in termini della derivata totale:

 \frac {D \rho} {D t} + \rho \nabla \cdot \vec u = 0 \, .

Punto di vista euleriano[modifica | modifica wikitesto]

Lo stesso principio di conservazione, dal punto di vista euleriano, può essere così espresso:

« la variazione di massa contenuta in un volume fisso eguaglia la differenza tra i flussi di massa entranti ed i flussi di massa uscenti (opposto del flusso netto). »

Un generico flusso di massa per unità di superficie, che passi attraverso una coppia di facce P e Q di un volume, è considerato come il prodotto tra la densità ρ del fluido, la componente della velocità in direzione perpendicolare alla faccia considerata e l'area della faccia medesima.

Considerando l'ipotesi di elemento infinitesimo possiamo approssimare il valore del flusso nel punto centrale di ogni faccia con il suo valore medio e calcolare il valore del flusso su una faccia a partire dal valore assunto sulla faccia precedente tramite una serie di Taylor troncata al primo ordine:

 \Phi_{P_x} = \rho u \, dy dz \qquad \Phi_{Q_x} = \rho u \, dy dz +\frac {\partial \rho u} {\partial x}  \, dx dy dz

dove con Px e Qx si sono indicate le facce normali (cioè perpendicolari) alla direzione x. Seguendo l'enunciato del principio, ovverosia calcolando la differenza dei flussi, otteniamo:

 \Phi_x = \Phi_{Q_x} - \Phi_{P_x} = \rho u \, dy dz + \frac {\partial \rho u} {\partial x} \, dx dy dz -\rho u \, dy dz = \frac {\partial \rho u} {\partial x} \, dx dy dz

Estendendo il ragionamento alle altre direzioni spaziali otteniamo che il flusso netto sarà uguale a:

 \Phi = \left( \frac {\partial \rho u} {\partial x} + \frac {\partial \rho v} {\partial y} + \frac {\partial \rho w} {\partial z} \right) \, dx dy dz = \left[ \nabla \cdot \left( \rho \vec u \right) \right] \, dx dy dz.

Ponendo ora la variazione di massa nel tempo uguale all'opposto del flusso netto:

 \frac {\partial \rho} {\partial t} \, dx dy dz = \left[ \nabla \cdot \left( \rho \vec u \right) \right] \, dx dy dz

ed infine, per unità di volume:

 \frac {\partial \rho} {\partial t} + \nabla \cdot \left( \rho \vec u \right) = 0

riottenendo l'espressione precedentemente mostrata.

Equazione di bilancio della quantità di moto[modifica | modifica wikitesto]

Exquisite-kfind.png Per approfondire, vedi legge di conservazione della quantità di moto.

Punto di vista lagrangiano[modifica | modifica wikitesto]

La conservazione della quantità di moto (definita come prodotto della massa per la velocità o, per unità di volume, della densità per la velocità) si esprime affermando che:

« la variazione temporale della quantità di moto di un sistema coincide con la risultante delle forze esterne al sistema »

e matematicamente:

 \frac{d \vec Q}{dt} = \vec F_e

dove, appunto, con Fe si è indicata la somma delle forze esterne, di massa (come ad esempio la forza di gravità) e di superficie (quali ad esempio le forze viscose).

Introducendo questa differenziazione nelle forze ed una formulazione integrale:

\frac{d}{dt} \int_{V(t)} \rho \vec u \, dV = \int_{V(t)} \rho \vec F_V \, dV + \int_{S(t)} \vec F_S \, dS

Il primo membro può essere trasformato in forma più conveniente mediante il teorema del trasporto di Reynolds:

 \frac{d}{dt} \int_{V(t)} \rho \vec u \, dV = \int_{V(t)} \left( \frac{D \rho \vec u}{Dt} + \rho \vec u \left( \nabla \cdot \vec u \right) \right) \, dV

che può essere ridotta nella forma:

 \frac{d}{dt} \int_{V(t)} \rho \vec u \, dV = \int_{V(t)} \rho \frac{D \vec u}{Dt} \, dV + \int_{V(t)} \vec u \left( \frac{D \rho}{Dt} + \rho \nabla \cdot \vec u \right) \, dV

dove l'ultimo integrale coincide con l'equazione di continuità ed è perciò nullo.

Se si applica il teorema della divergenza all'ultimo integrale dell'equazione della quantità di moto, sarà possibile scriverlo come integrale di volume. L'equazione si trasforma quindi come segue:

\int_{V(t)} \left( \rho \frac{D \vec u}{Dt} - \rho \vec F_V - \nabla \cdot \underline{\underline T } \right) \, dV = 0

dove con T con doppia sottolineatura si è indicato il tensore delle tensioni. Dato che l'uguaglianza precedente deve valere per qualsiasi arbitrario volume di integrazione, dovrà essere nullo l'integrando:

 \rho \frac{D \vec u}{Dt} = \rho \vec F_V + \nabla \cdot \underline{\underline T }

che esprime l'equazione della quantità di moto (per unità di volume).

Punto di vista euleriano[modifica | modifica wikitesto]

Il secondo principio della dinamica esprime la conservazione della quantità di moto e, per un elemento del fluido può essere enunciato come segue:

« la variazione, nel tempo, della quantità di moto del fluido contenuto nel volume di controllo τ, sommata al flusso netto di quantità di moto attraverso la superficie σ, uguaglia la risultante delle forze esterne agenti sul fluido contenuto nel volume stesso. »

Quindi, con formulazione integrale:

 \frac{d}{dt} \int_V \rho \vec u \, dV + \oint_S \left( \rho \vec u \right) \vec u \cdot \hat n \, dS = \int_V \vec F_V \, dV + \oint_S \underline{\underline T } \, dS

dove il volume (come la superficie S che lo racchiude) non è funzione del tempo.

Il tensore delle tensioni per un fluido[modifica | modifica wikitesto]

Exquisite-kfind.png Per approfondire, vedi relazioni costitutive e tensione interna.
Generica forza agente su di una superficie (orientata).

Il tensore delle tensioni o tensore degli sforzi è un tensore bidimensionale, caratterizzato da nove componenti Tik che rappresentano le tre componenti degli sforzi nelle tre direzioni spaziali di un certo sistema di riferimento cartesiano. In forma matematica:

 \underline{\underline T } = \begin{Bmatrix} T_{11} & T_{12} & T_{13} \\ T_{21} & T_{22} & T_{23} \\ T_{31} & T_{32} & T_{33} \\ \end{Bmatrix} .

Se si considera una generica superficie, orientata secondo il versore n come nella figura, e la risultante delle forze elementari che le molecole di fluido in prossimità della faccia positiva esercitano su quelle in prossimità della faccia negativa si ha che lo sforzo relativo alla superficie S ed al versore n è:

\vec T_n = \frac{\vec F_n}{S}

Il pedice n indica che la forza dipende dalle dimensioni e dall'orientamento della superficie, mentre lo sforzo dal solo orientamento.

Nel caso statico, cioè di un fluido in quiete, a differenza della meccanica dei solidi, l'unica forza sarà puramente normale e lo sforzo risultante sarà chiamato pressione. In un fluido in moto invece, l'attrito tra strati adiacenti e che si muovano a differente velocità darà luogo a forze risultanti oblique.

Al contrario dei solidi elastici, per i quali gli sforzi dipendono essenzialmente dalla deformazione attuale delle particelle, per i fluidi gli sforzi dipendono dalla velocità di deformazione.

Fluido non micropolare[modifica | modifica wikitesto]

Nel caso in cui il fluido sia non micropolare, cioè con tensore degli sforzi simmetrico, quindi Tik = Tki, cosicché le nove componenti sono ridotte a sei quantità indipendenti. Ciò è dovuto al fatto che i momenti meccanici agenti sulle facce di un certo volume (per esempio in forma di parallelepipedo rettangolo), rispetto ad un certo asse (per esempio z), che passi per il centro del volume, sono:

 \left| \vec M_z \right| = T'_{12} \, bc \, \frac{a}{2} - T''_{12} \, bc \left( - \frac{a}{2} \right) + T'_{21} \, ac \left( - \frac{b}{2} \right) - T'_{21} \, ac \, \frac{b}{2}

mentre l'equazione del momento meccanico di un moto attorno ad un baricentro di un parallelepipedo è:

 \left| \vec M_z \right| = I_z \, \omega_z = \rho \, abc \, \frac{a^2 + b^2}{12} \, \omega_z

dove con Iz si è indicato il momento d'inerzia attorno all'asse z e con ωz la velocità angolare. Eguagliando le precedenti espressioni si ottiene:

 \rho \, \frac{a^2 + b^2}{12} \, \omega_z = \frac{T'_{12} + T''_{12}}{2} - \frac{T'_{21} + T''_{21}}{2} .

Al tendere del volume a 0, le lunghezze a, b e c tenderanno a 0, mentre gli sforzi sulle facce opposte tenderanno ad un valore comune. Resta quindi:

 T_{12} - T_{21} = 0

che vale anche per gli altri assi.

Relazioni tra sforzi e velocità di deformazione: fluidi newtoniani isotropi[modifica | modifica wikitesto]

Exquisite-kfind.png Per approfondire, vedi fluido newtoniano.

Un fluido si definisce newtoniano quando la sua viscosità non varia con la velocità e, per questo motivo, la relazione matematica che lega il tensore degli sforzi alle componenti del tensore della velocità di deformazione è lineare.

Desiderando trovare le relazioni che legano sforzi e velocità di deformazione, analizziamo i casi più semplici per poi sommarne gli effetti (grazie alla linearità del problema), ricavando il caso generale.

Il caso più semplice in assoluto sarà il caso statico: come già osservato gli sforzi saranno puramente normali, mentre il tensore delle velocità di deformazione (che indicheremo con \underline{\underline \varepsilon }) è nullo per ipotesi. In termini matematici:

\begin{cases}
T_{ik} = - p \qquad &\mathrm{per \ } i = k \\
T_{ik} = 0 \qquad &\mathrm{per \ } i \ne k .\\
\end{cases}

Consideriamo adesso un flusso in moto, dove però, per un particolare sistema di riferimento cartesiano, gli sforzi siano puramente normali alle superfici di un elemento di forma parellelepipeda (sistema di riferimento degli assi principali di deformazione). Per esempio supponiamo che sia:

\begin{cases}
T_{ik} = 0 \qquad &\mathrm{per \ } i \ne k \\
T_{11} > T_{22} = T_{33} \, .\\
\end{cases}

Gli effetti del sistema di sforzi precedente su di un fluido, sono differenti nel caso di fluido isotropo (come ad esempio acqua ed aria) oppure anisotropo (come ad esempio il sangue, le cui molecole conferiscono al fluido proprietà diverse nelle differenti direzioni). L'esperienza fisica dimostra che i fluidi che interessano l'aerodinamica e l'idrodinamica sono fluidi newtoniani ed isotropi, altrimenti detti fluidi stokesiani. Analizzeremo perciò un fluido isotropo, dove cioè dovrà essere ε12 = 0:

\begin{cases}
\varepsilon_{ik} = 0 \qquad &\mathrm{per \ } i \ne k \\
\varepsilon_{11} > 0 \\
\varepsilon_{22} = \varepsilon_{33} < \varepsilon_{11} \, .\\
\end{cases}

Resta infine da considerare il caso più generale, dove cioè tutte le componenti degli sforzi saranno diversi da zero:

T_{ik} \ne 0 \quad \forall i;k \qquad \Rightarrow \qquad \varepsilon_{ik} \ne 0 \quad \forall i;k.

Ogni componente del tensore degli sforzi sarà una certa funzione, lineare per fluidi newtoniani, delle componenti del tensore delle velocità di deformazione. Sviluppando tale funzione in serie di Taylor (arrestata al primo ordine per la sua proprietà di linearità), si ottiene:

T_{ik} = f_0 + f_1(\varepsilon_{mn}).

Resta ora da ricavare tali funzioni lineari: trattando il problema in un sistema di riferimento particolare quale quello degli assi principali di deformazione, si ha:

\begin{cases}
T_{11} = f_0 + a_1 \varepsilon_{11} + b_1 \varepsilon_{22} + b_1 \varepsilon_{33}\\
T_{22} = f_0 + a_2 \varepsilon_{22} + b_2 \varepsilon_{11} + b_2 \varepsilon_{33}\\
T_{33} = f_0 + a_3 \varepsilon_{33} + b_3 \varepsilon_{11} + b_3 \varepsilon_{22}.\\
\end{cases}

Nel primo caso analizzato sarà quindi:

f_0 = -p.\

A causa del fatto che studiamo un fluido stokesiano, vi è inoltre completa equivalenza di comportamento tra le tre direzioni principali di deformazione x1, x2, x3 e quindi:

\begin{cases}
a_1 = a_2 = a_3 = a\\
b_1 = b_2 = b_3 = b\\
\end{cases}

e dunque il sistema iniziale si potrà scrivere come:

\begin{cases}
T_{11} = f_0 + (a-b) \varepsilon_{11} + b \sum_{i=1}^3 \varepsilon_{ii}\\
T_{22} = f_0 + (a-b) \varepsilon_{22} + b \sum_{i=1}^3 \varepsilon_{ii}\\
T_{33} = f_0 + (a-b) \varepsilon_{33} + b \sum_{i=1}^3 \varepsilon_{ii}.\\
\end{cases}

Infine, tenendo conto che

\sum_{i=1}^3 \varepsilon_{ii} = \nabla \cdot \vec u

e ponendo per comodità

\begin{cases}
a - b = 2\mu\\
b = \lambda\\
\end{cases}

si ottiene:

\begin{cases}
T_{11} = -p + \lambda \nabla \cdot \vec u + 2\mu \, \varepsilon_{11}\\
T_{22} = -p + \lambda \nabla \cdot \vec u + 2\mu \, \varepsilon_{22}\\
T_{33} = -p + \lambda \nabla \cdot \vec u + 2\mu \, \varepsilon_{33}\\
\end{cases}

dove il secondo termine a secondo membro descrive l'effetto della viscosità dovuto alla variazione di volume di una particella di fluido.

Non resta ora che generalizzare il sistema di equazioni precedente al caso di una terna di riferimento qualsiasi:

\begin{cases}
T_{kk} = -p + \lambda \nabla \cdot \vec u + 2\mu \, \varepsilon_{kk}\\
T_{ik} = 2\mu \, \varepsilon_{ik} \qquad i \ne k.\\
\end{cases}

La prima equazione del sistema precedente evidenzia il fatto che, nel caso generale, i tre sforzi normali sono differenti tra loro. La loro media è:

\begin{align}
\frac{T_{11} + T_{22} + T_{33}}{3} &= - p + \lambda \nabla \cdot \vec u + \frac{2}{3} \; \mu \left( \varepsilon_{11} + \varepsilon_{22} + \varepsilon_{33} \right) = \\
&= - p + \lambda \nabla \cdot \vec u + \frac{2}{3} \; \mu \nabla \cdot \vec u \\
&= - p + \mu' \nabla \cdot \vec u
\end{align}

dove con μ' si è indicata la viscosità di volume (o in terminologia anglosassone bulk viscosity), la quale descrive la differenza tra lo sforzo normale medio e la pressione di un fluido, dovuta alla viscosità. Il valore della viscosità di volume in genere è trascurabile per i gas, in particolare per quelli monoatomici.

Conservazione dell'energia[modifica | modifica wikitesto]

Exquisite-kfind.png Per approfondire, vedi legge di conservazione dell'energia.

Il primo principio della termodinamica, ovvero il principio di conservazione dell'energia può essere espresso dicendo che la variazione nell'unità di tempo dell'energia totale del fluido contenuto nel volume di controllo sommata al flusso netto di energia totale attraverso le facce del volume di controllo uguaglia la somma della potenza delle forze agenti sull'elemento di fluido e del flusso netto di energia termica trasmessa all'elemento di fluido per conduzione.

Come si nota in questa formulazione viene trascurata l'energia trasmessa all'elemento per irraggiamento. Formalizzando matematicamente questo principio si sfrutterà il concetto di energia totale per unità di massa E che è uno scalare definito come:

E = e + \frac {1} {2} V^2

cioè la somma tra l'energia interna delle molecole e l'energia meccanica degli elementini di fluido.

Nell'enunciato si parla di flusso netto di energia totale: come per la quantità di moto si indicherà questo flusso come il prodotto tra il flusso di massa e l'energia totale per unità di massa trasportata in ogni direzione:

 \Phi_E=\frac {\partial E \rho u} {\partial x} + \frac {\partial E \rho v} {\partial y} + \frac {\partial E \rho w} {\partial z}

La potenza degli sforzi agenti sull'elementino di fluido considerato comprende sia la potenza sviluppata dagli sforzi viscosi del tensore  S sia gli sforzi associati alla pressione.

Ricorrendo alla definizione di potenza come prodotto di una forza per una velocità, si potrà scrivere:

P_S=\frac {\partial (S_{xx}u+S_{yx}v+S_{zx}w)} {\partial x}+\frac {\partial (S_{xy}u+S_{yy}v+S_{zy}w)} {\partial y}+\frac {\partial (S_{xz}u+S_{yz}v+S_{zz}w)} {\partial z}

per quanto riguarda gli sforzi viscosi, mentre per la pressione sarà:

 P_p=-\left(\frac {\partial p u} {\partial x} + \frac {\partial p v} {\partial y} + \frac {\partial p w} {\partial z}\right)

La potenza delle forze di campo si definisce come:

P_c=\rho a_xu+\rho a_yv+\rho a_zw

Per quanto riguarda la potenza termica trasmessa per conduzione attraverso le facce dell'elementino è necessaria la definizione di un vettore \vec q=[q_x, q_y, q_z]^T flusso termico. Sarà possibile scrivere:

-\left(\frac {\partial q_x} {\partial x} + \frac {\partial q_y} {\partial y} + \frac {\partial q_z} {\partial z}\right)

L'equazione completa che formalizza il primo principio della termodinamica per i fluidi in movimento sarà quindi:

 \frac {\partial \rho E} {\partial t} + \frac {\partial E \rho u} {\partial x} + \frac {\partial E \rho v} {\partial y} + \frac {\partial E \rho w} {\partial z} = -\left(\frac {\partial p u} {\partial x} + \frac {\partial p v} {\partial y} + \frac {\partial p w} {\partial z}\right) +
+ \frac {\partial (S_{xx}u+S_{yx}v+S_{zx}w)} {\partial x}+\frac {\partial (S_{xy}u+S_{yy}v+S_{zy}w)} {\partial y}+\frac {\partial (S_{xz}u+S_{yz}v+S_{zz}w)} {\partial z} +
+ \rho a_xu+\rho a_yv + \rho a_zw - \left(\frac {\partial q_x} {\partial x} + \frac {\partial q_y} {\partial y} + \frac {\partial q_z} {\partial z}\right)

Osservazioni e chiusura del problema[modifica | modifica wikitesto]

Le 3 equazioni (due equazioni scalari ed un'equazione vettoriale) appena derivate sono insufficienti, da sole, alla chiusura del problema della determinazione del campo di moto del fluido. Infatti le equazioni contengono 20 incognite:

Queste equazioni sono del tutto generali e per la loro applicazione è necessaria una sorta di specializzazione delle stesse alla situazione di lavoro.

Per la chiusura del problema è quindi necessario definire le proprietà termofisiche del fluido in esame (che permettono di definire la conducibilità termica, la densità, l'energia interna e una o più equazioni di stato in grado di determinare anche temperatura e pressione) e il campo di forze in cui si muove (determinando il vettore di accelerazioni di campo). Inoltre si osserva che il tensore degli sforzi viscosi S è simmetrico, con la conseguenza che le incognite effettivamente contenute sono 6 e non 9 e sono determinabili sperimentalmente o teoricamente specificando il tipo di fluido. Saranno successivamente necessarie le condizioni iniziali e le condizioni al contorno, trattandosi di equazioni differenziali (problema di Cauchy o problema di von Neumann).

Equazioni di Eulero[modifica | modifica wikitesto]

Exquisite-kfind.png Per approfondire, vedi equazioni di Eulero.

Relazioni di salto[modifica | modifica wikitesto]

Discontinuità di contatto[modifica | modifica wikitesto]

Onda d'urto[modifica | modifica wikitesto]

Exquisite-kfind.png Per approfondire, vedi onda d'urto (fluidodinamica).

Le equazioni in forma adimensionale[modifica | modifica wikitesto]

Le equazioni scritte nei paragrafi precedenti sono in forma dimensionale, nel senso che ogni termine possiede dimensioni fisiche della grandezza considerata:

  •  \mathrm{ \left[ \frac{kg}{m^3 \, s} \right] \qquad } nella prima equazione;
  • \mathrm{ \left[ \frac{kg}{m^2 \, s^2} \right] \qquad } nelle tre equazioni della quantità di moto;
  • \mathrm{ \left[ \frac{kg}{m \, s^3} \right] \qquad } nell'ultima equazione.

Di conseguenza, volendo confrontare tra loro i numerosi coefficienti per sapere quali di essi sia il più preponderante nei vari casi in esame, bisognerebbe calcolare il valore di ogni singolo termine. Un metodo pratico per ovviare a questa necessità è quello di dividere ogni coefficiente per una certa grandezza omogenea di riferimento, in tal modo i coefficienti risulteranno adimensionali. Queste grandezze di riferimento saranno scelte in base alle condizioni al contorno ed alle condizioni iniziali del particolare problema fluidodinamico che si vuole esaminare. Qui sono indicate con il pedice 0 (zero):

\rho^* = \frac{\rho}{\rho_0} \qquad {u_i}^* = \frac{u_i}{U_0} \qquad p^* = \frac{p}{p_0} \qquad t^* = \frac{t}{t_0} \qquad x^* = \frac{x_i}{L_0} \qquad T^* = \frac{T}{T_0}

L'equazione di conservazione della massa[modifica | modifica wikitesto]

L'equazione di conservazione della massa scritta nella forma:

\frac{D \rho}{D t} + \rho \vec{\nabla} \cdot \vec{u} = \frac{\partial \rho}{\partial t} + (\vec{u} \cdot \vec{\nabla}) \rho + \rho \vec{\nabla} \cdot \vec{u} = 0

può essere resa adimensionale esprimendola nella forma:

\mathrm{St} \; \frac{\partial \rho^*}{\partial t^*} + (\vec{u}^* \cdot \vec{\nabla}^*) \rho^* + \rho^* \vec{\nabla}^* \cdot \vec{u}^* = 0

dove con il simbolo St si è indicato il gruppo adimensionale, detto numero di Strouhal:

  • \mathrm{St} = \frac{L_0}{U_0 \; t_0}.

Le equazioni di conservazione della quantità di moto[modifica | modifica wikitesto]

Le equazioni di conservazione della quantità di moto possono essere adimensionalizzate nella forma:

\mathrm{St} \left( \rho^* \frac{\partial \vec {u}^*}{\partial t^*} \right) + \rho^* \left( \vec{u}^* \cdot \vec{\nabla}^* \right) \vec{u}^* = - \frac{1}{\mathrm{Ru}} \vec{\nabla}^* p^* - \frac{1}{\mathrm{Fr}^2} \rho^* \vec k + \frac{1}{\mathrm{Re}} {\nabla^*}^2 \vec{u}^* + \frac{1}{3 \, \mathrm{Re}} \vec{\nabla}^* \left( \vec{\nabla}^* \cdot \vec{u}^* \right)

dove i simboli indicano i seguenti gruppi adimensionali:

Nel caso in cui la viscosità dinamica \mu non sia costante, si troverà un valore di riferimento \mu_0 e si utilizzerà all'interno dell'equazione il valore adimensionale \mu^* .

L'equazione di conservazione dell'energia termica[modifica | modifica wikitesto]

L'equazione di conservazione dell'energia termica, dato che quella dell'energia meccanica condurrebbe a gruppi adimensionali già visti per le equazioni della quantità di moto, viene espressa in funzione di termini adimensionali:

\mathrm{St} \left( \rho^* \frac{\partial T^*}{\partial t^*} \right) + \rho^* \left( \vec{u}^* \cdot \vec{\nabla}^* \right) T^* = \mathrm{St} \,  \frac{\mathrm{Ec}}{\mathrm{Ru}} \, \frac{\partial p^*}{\partial t^*} + \frac{\mathrm{Ec}}{\mathrm{Ru}} \left( \vec{u}^* \cdot \vec{\nabla}^* \right) p^* + \frac{\mathrm{Ec}}{\mathrm{Re}} \Phi^* + \frac{1}{\mathrm{Pr} \, \mathrm{Re}} {\nabla^*}^2 T^* + \frac{\mathrm{Nu}}{\mathrm{Pr} \, \mathrm{Re}} \rho^* \dot{q}^*

dove i simboli indicano i seguenti gruppi adimensionali:

Bibliografia[modifica | modifica wikitesto]

  • (EN) R. Byron Bird, Warren E. Stewart; Edwin N. Lightfoot, Transport Phenomena, 2ª ed., New York, Wiley, 2005, ISBN 0-470-11539-4.

Voci correlate[modifica | modifica wikitesto]