Strato limite di quantità di moto

In fluidodinamica, lo strato limite, più precisamente detto strato limite di quantità di moto, è una zona ideale di un flusso dove gli effetti causati dalla prossimità di una parete sono molto più marcati rispetto ad una zona esterna. Solitamente tale zona è individuata in uno strato di fluido nelle immediate vicinanze di una superficie solida.

Zona interessata da uno strato limite laminare, dove cioè la componente orizzontale delle velocità è molto più piccola della componente orizzontale all'esterno dell'area di influenza della parete.

Teoria dello strato limite

La teoria dello strato limite si basa sull'intuizione dell'aerodinamico tedesco Ludwig Prandtl, che dimostrò poi con evidenze sperimentali agli inizi del '900, la validità delle ipotesi fatte. Prandtl indicò come strato limite la zona al cui interno la velocità differisce per più dell'1% dalla velocità della porzione esterna di fluido, assumendo per l'esterno il modello di moto a potenziale. Prandtl ipotizzò inoltre che le dimensioni dello strato limite fossero piccole rispetto alle dimensioni del campo di moto esterno, assumendo che lo strato inizi in corrispondenza del bordo di ingresso della lastra. Di conseguenza, per le ipotesi fatte, le linee di corrente rimangono indisturbate fin tanto che non intersecano il bordo dello strato limite, iniziando da questo punto in poi a deflettere leggermente a causa del rallentamento della corrente.

Per correnti ad alti numeri di Reynolds è possibile confinare gli effetti della viscosità all'interno di una zona limitata adiacente alla parete; il flusso nella zona esterna si può considerare potenziale. Questo tipo di approssimazione semplifica notevolmente le equazioni di Navier-Stokes, permettendo di ottenere soluzioni approssimate. Si ricorda che le equazioni di Navier-Stokes sono equazioni alle derivate parziali che, allo stato attuale degli studi matematici, non sono risolubili per via analitica (se non in casi molto semplificati), ma in generale solo per via numerica.

Lo strato limite può essere di tipo laminare, dove i filetti fluidi sono appunto lamine che seguono il contorno del corpo, oppure turbolento, dove i filetti di fluido seguono linee intricate e non stazionarie. Lo strato limite può transitare da laminare a turbolento attraverso una zona detta regione di transizione, ma non potrà mai tornare spontaneamente laminare. La transizione dallo strato laminare a quello turbolento avviene in una zona di transizione che presenta caratteristiche di non stazionarietà. La transizione viene infine favorita dall'aumentare della velocità (o più correttamente dal numero di Reynolds) e dalla rugosità superficiale del corpo e dalla forma meno affusolata del corpo.

Trattazione dello strato limite

Come già accennato precedentemente la trattazione dello strato limite dipende in buona parte dalle condizioni iniziali a seconda della velocità della corrente che lambisce la parete interessata e di conseguenza dal rispettivo numero di Reynolds della corrente. Ciò nonostante si può pervenire ad una trattazione generale dello strato limite a prescindere dalla sua natura laminare o turbolenta applicando in entrambi casi le ipotesi che il moto sia di tipo:

• permanente;
• bidimensionale;
• trascurabile l'effetto della forza peso;
• unidirezionale.

Strato limite laminare

Nel caso dello strato limite laminare il problema può essere analizzato ricorrendo all'utilizzo delle equazioni di bilancio di massa per un volume di controllo e utilizzando successivamente il bilancio di quantità di moto. Applicando quindi il bilancio di massa,

$\oint _{A}\rho \left({\vec {u}}\cdot d{\vec {A}}\right)=0$

dove con ${\vec {u}}$ viene indicata la velocità del fluido che lambisce la parete, con ${\vec {A}}$ una superficie orientata e con δ lo spessore dello strato limite laminare. Esplicitando i singoli termini, indicando con y la direzione perpendicolare alla parete e con x la direzione parallela sia alla parete che alla corrente:

$-\int _{0}^{\delta }\rho \,{\vec {u}}\,dy+\int _{0}^{\delta }\rho \,{\vec {u}}\,dy+{\frac {\partial }{\partial x}}\left(\int _{0}^{\delta }\rho \,{\vec {u}}\,dy\right)dx-q_{l}\,dx=0$

dove q_l rappresenta il flusso di massa entrante dal bordo superiore. Dopo aver semplificato si ricava la seguente espressione:

${\frac {\partial }{\partial x}}\left(\int _{0}^{\delta }\rho {\vec {u}}dy\right)=q_{l}\,$

A questo punto, si passa ad applicare il bilancio di quantità di moto:

$\sum _{k=1}^{n}{\vec {F}}_{e}=\oint _{A}\rho {\vec {u}}\left({\vec {u}}\cdot d{\vec {A}}\right)$

per il termine a sinistra dell'equazione del bilancio di quantità di moto, che rappresenta la somma delle forze esterne agenti sul volume di controllo (volume di profondità unitaria), si ottiene

$\sum _{k=1}^{n}{\vec {F}}_{e}=p\delta -\left(p+{\frac {\partial p}{\partial x}}dx\right)\left(\delta +{\frac {\partial \delta }{\partial x}}dx\right)+{\frac {1}{2}}\left(p+p+{\frac {\partial p}{\partial x}}dx\right){\frac {\partial \delta }{\partial x}}dx-\tau dx$

dove con τ si è indicato lo sforzo tangenziale sulla parete. Si noti inoltre che le pressioni agenti sulle facce si suppongano essere costanti, nonostante, in condizioni di corrente laminare, l'andamento delle pressioni ortogonalmente alle linee di corrente, per un fluido incomprimibile, sia di tipo idrostatico. Questa semplificazione nel caso di fluidi incomprimibili risulta comunque accettabile in virtù del fatto che lo strato limite si suppone essere di dimensioni molto ridotte. A questo punto semplificando e trascurando tutti gli infinitesimi di ordine superiore si ottiene infine

$\sum _{k=1}^{n}{\vec {F}}_{e}=-{\frac {\partial p}{\partial x}}dx\delta -\tau dx.$

Per quanto riguarda invece il termine a destra nell'equazione del bilancio di quantità di moto si ottiene

$\oint _{A}\rho {\vec {u}}\left({\vec {u}}\cdot d{\vec {A}}\right)=-\int _{0}^{\delta }\rho {\vec {u}}^{2}dy+\int _{0}^{\delta }\rho {\vec {u}}^{2}dy+{\frac {\partial }{\partial x}}\left(\int _{0}^{\delta }\rho {\vec {u}}^{2}dy\right)dx-{\vec {U}}_{0}q_{l}dx.$

Ricordando quindi l'espressione ricavata in precedenza per il termine $q_{l}$ ed inserendolo nell'equazione si ricava

$\oint _{A}\rho {\vec {u}}({\vec {u}}\cdot d{\vec {A}})=+{\frac {\partial }{\partial x}}\left(\int _{0}^{\delta }\rho {\vec {u}}^{2}dy\right)dx-{\vec {U}}_{0}{\frac {\partial }{\partial x}}\left(\int _{0}^{\delta }\rho {\vec {u}}dy\right)dx.$

Uguagliando i due termini ricavati si ottiene la seguente equazione

$-{\frac {\partial p}{\partial x}}\delta -\tau ={\frac {\partial }{\partial x}}\left(\int _{0}^{\delta }\rho {\vec {u}}^{2}dy\right)-{\vec {U}}_{0}{\frac {\partial }{\partial x}}\left(\int _{0}^{\delta }\rho {\vec {u}}dy\right).$

anche nota come equazione di von Kármán che, nel caso di moto a potenziale nello strato esterno allo strato limite laminare, supponendo che

${\frac {\partial p}{\partial x}}=\mathrm {costante}$

e

${\frac {\partial {\vec {U}}_{0}}{\partial x}}=\mathrm {costante}$

si semplifica nella seguente espressione

$\tau =-{\frac {\partial }{\partial x}}\left(\int _{0}^{\delta }\rho {\vec {u}}^{2}dy\right)+{\vec {U}}_{0}{\frac {\partial }{\partial x}}\left(\int _{0}^{\delta }\rho {\vec {u}}dy\right).$

L'equazione può inoltre scriversi anche in forma adimensionale, raccogliendo alcuni termini nel seguente modo

${\begin{aligned}{\frac {y}{\delta }}&=\eta \\{\frac {\vec {u}}{U}}_{0}&=u_{1}\end{aligned}}$

ottenendo quindi l'equazione di von Kármán in forma adimensionale

${\frac {\tau }{\rho U_{0}}}=-{\frac {\partial }{\partial x}}\delta \left(\int _{0}^{1}{\vec {u}}_{1}(1-{\vec {u}}_{1})d\eta \right){.}$

A questo punto l'equazione può essere integrata facilmente non appena sia nota la distribuzione delle velocità. Per fare ciò è quindi necessario conoscere la distribuzione delle velocità. Nel caso di strato limite laminare il problema è stato risolto da Blasius nel 1908 secondo cui la soluzione è espandibile in serie di potenze. Tuttavia la soluzione di Blasius risulta laboriosa per cui sono state proposte diverse altre soluzioni. Tra quelle che approssimano meglio la distribuzione delle velocità va citata sicuramente la distribuzione di Pohlhausen proposta nel 1921 che approssima la funzione con un polinomio di 4 grado del tipo

$\sum _{i=0}^{4}a_{i}\eta _{i}=u_{1}.$

Teoria di Blasius

Per ricavare la distribuzione delle velocità secondo il modello di Blasius è necessario un approccio diverso da quello precedente. In un problema fluidodinamico la velocità del fluido immediatamente adiacente ad una superficie solida è considerata uguale alla velocità della parete stessa. A causa del fatto che la velocità varia con continuità dalla parete verso l'esterno e che il suo valore tende asintoticamente al valore esterno, lo spessore dello strato limite deve essere arbitrariamente definito, per convenzione, come lo strato di fluido all'interno del quale la velocità varia da zero fino al 99% (o 90%, o 95% a seconda della particolare convenzione che si preferisce adottare) della velocità della corrente fluida indisturbata (detta appunto velocità asintotica).

Infatti gli sforzi tangenziali diminuiscono con continuità allontanandosi dalla parete. Un osservatore invece percepisce una variazione marcata degli sforzi a causa del fatto che lo spessore dello strato limite è molto piccolo rispetto alle dimensioni del corpo.

Il modello base di questa teoria prende in considerazione l'equazione di stato dei gas perfetti:

$p=\rho RT\!\$

dove p è la pressione, ρ la densità, R la costante specifica dei gas perfetti ed infine T la temperatura in kelvin.

E l'equazione di continuità della massa:

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

dove con il vettore V si è indicato il vettore velocità fluidodinamica, della quantità di moto:

${\frac {\partial \left(\rho u_{i}\right)}{\partial t}}+{\frac {\partial }{\partial x_{i}}}\left(\rho u_{i}u_{j}+p\delta _{ij}\right)={\frac {\partial \tau _{ij}}{\partial x_{i}}}+\rho f_{i}$

e dell'energia:

${\frac {\partial }{\partial t}}\left[\rho \left(e+{\frac {1}{2}}V^{2}\right)\right]+{\frac {\partial }{\partial x_{i}}}\left[\rho \left(e+{\frac {1}{2}}V^{2}\right)u_{i}\right]=\rho f_{i}u_{i}+{\frac {\partial \left(t_{ij}u_{i}\right)}{\partial x_{j}}}+{\frac {\partial }{\partial x_{i}}}\left(k{\frac {\partial T}{\partial x_{i}}}\right){.}$

Il sistema delle sei equazioni (tre scalari ed una vettoriale) viene semplificato con le seguenti ipotesi:

• il fluido è continuo (numero di Knudsen molto minore di 1), altrimenti non sono valide le equazioni precedenti;
• all'interno dello strato limite le variazioni lungo la coordinata parallela al flusso sono molto minori delle variazioni lungo la coordinata normale al flusso;
• flusso stazionario;
• flusso incomprimibile;
• flusso viscoso.

Infatti, indicando con δ lo spessore dello strato limite (funzione in generale di x, coordinata parallela alla parete ed al flusso) e con L una dimensione caratteristica del flusso (la lunghezza della parete o il diametro di un condotto, per esempio), si avrà certamente:

$\delta \ll L$

e perciò si può ipotizzare che le derivate rispetto ad x siano trascurabili rispetto a quelle rispetto ad y. Considerando fluidi incomprimibili (cioè a densità costante) e solo le direzioni x ed y (quindi un flusso bidimensionale) si semplifica dapprima l'equazione di continuità della massa, che per comodità scriviamo in forma vettoriale:

${\frac {D\rho }{Dt}}+\rho \nabla \cdot {\vec {V}}=0\quad \Rightarrow \quad \nabla \cdot {\vec {V}}=0$

ovvero, indicando con i pedici le derivate parziali:

$u_{x}+v_{y}=0.$

Analogamente viene semplificata anche l'equazione di conservazione della quantità di moto, in forma vettoriale:

$\rho {\frac {D{\vec {V}}}{Dt}}+\nabla p=\mu \nabla ^{2}{\vec {V}}+\left(\mu +\lambda \right)\nabla \left(\nabla \cdot {\vec {V}}\right)+\rho {\vec {f}}$

che diventa:

$\rho \left({\vec {V}}\cdot \nabla \right){\vec {V}}+\nabla p=\mu \nabla ^{2}{\vec {V}}$

ed in forma estesa:

${\begin{cases}uu_{x}+vu_{y}+{\frac {1}{\rho }}p_{x}={\mathcal {D}}_{p}u_{xx}+{\mathcal {D}}_{p}u_{yy}\\uv_{x}+vv_{y}+{\frac {1}{\rho }}p_{y}={\mathcal {D}}_{p}v_{xx}+{\mathcal {D}}_{p}v_{yy}\end{cases}}$

dove D_p rappresenta la diffusività cinematica.

Da queste equazioni si trascurano le derivate della componente verticale della velocità e la derivata seconda rispetto ad x della componente parallela al flusso della velocità. Ciò a causa del fatto che la velocità v è molto piccola rispetto ad u e le sue variazioni quindi sono minime. La mancanza del termine $u_{xx}$ è dovuta al fatto che le variazioni lungo l'asse x siano notevolmente più piccole rispetto a quelle lungo l'asse y. Quest'ultima approssimazione fa in modo che la soluzione non risenta di tutto il dominio di integrazione, ovverosia ciò che succede a valle non perturba la soluzione a monte, semplificando notevolmente il calcolo. Il sistema delle equazioni di conservazione della massa e della quantità di moto diviene perciò:

${\begin{cases}u_{x}+v_{y}=0\\uu_{x}+vu_{y}=-{\frac {1}{\rho }}p_{x}+{\mathcal {D}}_{p}u_{yy}\\p_{y}=0\end{cases}}$

Dette equazioni di Prandtl, anche se più propriamente solo la seconda equazione del sistema è l'equazione di Prandtl.

Accanto allo strato limite fluidodinamico, nel caso siano presenti fenomeni di trasmissione del calore nel fluido per via convettiva, può anche essere definito uno strato limite termico, come lo strato di fluido all'interno del quale la temperatura varia dal valore di parete fino al 99% (o altre percentuali a seconda della convenzione, valgono infatti le stesse considerazioni presentate per lo strato limite meccanico) del valore di temperatura esterno. L'ordine di grandezza del rapporto fra lo spessore dello strato limite fluidodinamico e quello dello strato limite termico è fornito dal numero di Prandtl.

Lo strato limite, di rilevante importanza nei problemi di aerodinamica ed in particolare nello studio delle superfici alari, è estendibile anche alla meteorologia ad allo studio delle correnti atmosferiche. È dovuto all'effetto della viscosità del fluido ed è collegato all'effetto Leidenfrost e, in particolare, al valore del numero di Reynolds (Re).

All'esterno dello strato limite il valore del numero di Reynolds è sufficientemente alto; in queste condizioni la regione esterna del flusso può essere considerata come ideale (viscosità trascurabile), che soddisfa cioè le equazioni di Eulero.

Nella regione vicino ai contorni solidi, invece questa approssimazione non può essere fatta perché bisogna rispettare le condizioni al contorno, che sono la condizione di aderenza (la velocità delle particelle di fluido a contatto con la parete deve essere nulla, in inglese è indicata come no slip condition) per lo strato limite della velocità, e la continuità del campo di temperatura sulla parete per lo strato limite termico.

Si considera quindi una lastra piana di apertura b e lunghezza $\ell$. Partendo dalle equazioni di conservazione della massa e della quantità di moto, e considerando che nello strato limite $p_{x}=0$, si ha

$\left\{{\begin{array}{l}u_{x}+v_{y}=0\\uu_{x}+vu_{y}={\mathcal {D}}_{p}u_{yy}\end{array}}\right.$

si introduce la variabile $\eta =y{\sqrt {\frac {V_{\infty }}{{\mathcal {D}}_{p}x}}}$ e la funzione di corrente $\Psi ={\sqrt {V_{\infty }x{\mathcal {D}}_{p}}}f(\eta )$, in modo che la soluzione trovata sarà simile per ogni x.

Si avrà

${\begin{array}{l}{\displaystyle u={\frac {\partial \Psi }{\partial y}}={\frac {\partial \Psi }{\partial \eta }}{\frac {\partial \eta }{\partial y}}={\sqrt {V_{\infty }x{\mathcal {D}}_{p}}}f'{\sqrt {\frac {V_{\infty }}{{\mathcal {D}}_{p}x}}}=V_{\infty }f'}\\{\displaystyle v=-{\frac {\partial \Psi }{\partial x}}=-\left[{\frac {1}{2}}{\sqrt {\frac {V_{\infty }{\mathcal {D}}_{p}}{x}}}f(\eta )-{\frac {1}{2}}{\sqrt {\frac {V_{\infty }{\mathcal {D}}_{p}}{x}}}f'\left(y{\sqrt {\frac {V_{\infty }}{{\mathcal {D}}_{p}x}}}\right)\right]={\frac {1}{2}}{\sqrt {\frac {V_{\infty }{\mathcal {D}}_{p}}{x}}}\left[\eta f'-f\right]}\end{array}}$

Derivando e sostituendo nelle equazioni di conservazione si ottiene

$f\cdot f''+2f'''=0$

Le condizioni al contorno sono

${\begin{array}{crclr}{\mbox{in }}x=0:&u=V_{\infty }\quad \forall y&\rightarrow &\eta =\infty :&f'=1\\{\mbox{in }}y=0:&u=v=0\quad \forall x&\rightarrow &\eta =0:&f=f'=0\end{array}}$

Per un fluido newtoniano gli sforzi viscosi di taglio in prossimità della parete risultano uguali a $\tau =\left.\mu {\frac {\partial u}{\partial y}}\right|_{y=0}.$

Si ha che

$u_{y}=V_{\infty }{\frac {\partial f'}{\partial y}}=V_{\infty }{\frac {\partial f'}{\partial \eta }}{\frac {\partial \eta }{\partial y}}=V_{\infty }{\sqrt {\frac {V_{\infty }}{{\mathcal {D}}_{p}x}}}f''$

quindi

$\tau =\mu u_{y}(0)=\mu V_{\infty }{\sqrt {\frac {V_{\infty }}{{\mathcal {D}}_{p}x}}}f''(0)$

Integrando lungo la lastra (e considerando le due facce della lastra bagnata dal fluido) è possibile trovare il valore della resistenza d'attrito

$D=2b\int _{0}^{\ell }\tau {\mbox{d}}x=2b\mu V_{\infty }{\sqrt {\frac {V_{\infty }}{{\mathcal {D}}_{p}}}}f''(0)\int _{0}^{\ell }{\frac {{\mbox{d}}x}{\sqrt {x}}}=2b\mu V_{\infty }{\sqrt {\frac {V_{\infty }}{{\mathcal {D}}_{p}}}}2{\sqrt {\ell }}f''(0)$

o in termini adimensionali

$C_{D}={\frac {D}{{\frac {1}{2}}\rho V_{\infty }^{2}b\cdot \ell }}={\frac {2\cdot 4\cdot f''(0)}{\sqrt {Re_{\ell }}}}$

Dove si è introdotto il numero di Reynolds locale:

$Re_{\ell }={\frac {V_{\infty }\ell }{{\mathcal {D}}_{p}}}$

inoltre il valore di $f''$ è noto e pari a $f''(0)={\frac {1{,}328}{4}}.$

Il coefficiente di forza medio sulla superficie della lastra risulta uguale a

${\bar {C}}_{f}={\frac {C_{D}}{2}}={\frac {1{,}328}{\sqrt {Re_{\ell }}}}.$

Distribuzione delle velocità secondo Pohlhausen

Come accennato in precedenza l'approssimazione di Pohlhausen è data da un polinomio di quarto grado del tipo

$\sum _{i=0}^{4}a_{i}\eta _{i}=u_{1}{.}$

Per potere quindi ricavare la distribuzione delle velocità all'interno dello strato limite sarà necessario disporre di cinque condizioni al contorno. Per lo strato adiacente alla parete, per $\eta =0$ le condizioni sono date da

${\begin{aligned}u_{1}&=0\\{\frac {\partial ^{2}u_{1}}{\partial \eta ^{2}}}&=0\end{aligned}}$

mentre le condizioni al bordo esterno per η = 1 sono

${\begin{aligned}u_{1}&=1\\{\frac {\partial u_{1}}{\partial \eta }}&=0\\{\frac {\partial ^{2}u_{1}}{\partial \eta ^{2}}}&=0{.}\end{aligned}}$

La distribuzione di velocità che si ottiene nella seguente maniera è infine

$u_{1}=2\eta -2\eta ^{3}+\eta ^{4}{.}$

Inserendo a questo punto la seguente espressione nell'equazione di von Karman è quindi possibile svolgere l'integrale dal quale si ottiene

${\frac {\tau }{\rho U_{0}}}=-{\frac {\partial }{\partial x}}\delta \int _{0}^{1}(2\eta -2\eta ^{3}+\eta ^{4})(1-2\eta +2\eta ^{3}-\eta ^{4})d\mathbf {\eta } ={\frac {36}{315}}{\frac {\partial \delta }{\partial x}}$

${\frac {\tau }{\rho U_{0}}}={\frac {36}{315}}{\frac {\partial \delta }{\partial x}}{.}$

Tuttavia l'espressione precedente contiene ancora due incognite $\delta$ e $\tau$ e necessita pertanto ancora di una seconda equazione per essere risolta. Questa è appunto fornita dal fatto che il fluido sia newtoniano. Si applica quindi l'equazione di Newton

$\tau =\mu {\frac {\partial u}{\partial x}}=\mu {\frac {U_{0}}{\delta }}{\frac {\partial u_{1}}{\partial \eta }}=\mu {\frac {U_{0}}{\delta }}2{.}$

Inserendo quindi il risultato così ottenuto nella relazione

${\frac {\tau }{\rho U_{0}}}={\frac {36}{315}}\,{\frac {\partial \delta }{\partial x}}$

ed integrando si ricava

${\frac {1}{2}}\delta ^{2}=17{,}03\,{\frac {\mu }{\rho U_{0}}}x+cost$

dove la costante di integrazione può essere posta 0 in virtù del fatto che la velocità sul bordo inferiore adiacente alla parete è nullo. Rimanipolando l'espressione appena ricavata dividendo per x e facendo la radice quadra

${\frac {\delta }{x}}=5{,}836{\sqrt {\frac {\mu }{\rho U_{0}x}}}$

si può quindi riscrivere l'espressione per lo sforzo tangenziale come

$\tau =\mu {\frac {\partial u}{\partial x}}=2\mu U_{0}{\frac {1}{5{,}836{\sqrt {\frac {\mu }{\rho U_{0}x}}}}}=0{,}343\,\mu U_{0}{\frac {1}{\sqrt {\frac {\mu }{\rho U_{0}x}}}}$

ossia

${\frac {\tau }{{\frac {1}{2}}\rho U_{0}^{2}}}=0{,}685\,{\frac {1}{\sqrt {\frac {\mu }{\rho U_{0}x}}}}$

dove il termine

$C_{r}={\frac {\tau }{{\frac {1}{2}}\rho U_{0}^{2}}}$

è anche noto come coefficiente di resistenza locale. Si noti che nel caso si fosse applicata la trattazione di Blasius il valore ricavato per il coefficiente di resistenza locale espresso come nel caso precedente sarebbe stato

$C_{r}=0{,}664\,{\frac {1}{\sqrt {\frac {\mu }{\rho U_{0}x}}}}.$

A questo punto noto il coefficiente di resistenza locale e nota la geometria della parete è possibile calcolare la resistenza di una piastra immersa in un fluido.

Strato limite turbolento

Per valori del numero di Reynolds superiori a 3 · 10⁵ lo strato limite inizia a divenire turbolento e la trattazione vista in precedenza non è più approssimata. Nel caso in cui il numero di Reynolds sia costante, il problema può comunque essere studiato nell'ipotesi di moto statisticamente stazionario, applicando la seguente scomposizione della velocità

$u_{i}\left({\vec {x}},t\right)={\overline {u_{i}}}\left({\vec {x}},t\right)+u_{i}'\left({\vec {x}},t\right)$

ovvero scomponendo la velocità in un termine velocità media che si identifica con l'espressione ${\overline {u_{i}}}\left({\vec {x}},t\right)$ e con un secondo termine che invece rappresenta la variazione media che si esprime come $u_{i}'\left({\vec {x}},t\right)$, nelle equazioni di Navier-Stokes, si ottengono le equazioni di Navier-Stokes mediate. Omettendo la sovralineatura sui termini velocità media e supponendo sempre che il moto sia di tipo bidimensionale, la prima equazione si scriverà dunque

$\rho \left(u_{x}{\frac {\partial u_{x}}{\partial x}}+u_{y}{\frac {\partial u_{x}}{\partial y}}\right)=-{\frac {\partial p}{\partial x}}+{\frac {\partial }{\partial x}}\left(\mu {\frac {\partial u_{x}}{\partial x}}-\rho u_{x}'u_{x}'\right)+{\frac {\partial }{\partial y}}\left(\mu {\frac {\partial u_{x}}{\partial y}}-\rho u_{x}'u_{y}'\right).$

Si possono adesso fare le seguenti ipotesi

$u_{x}\gg u_{y}\$

${\frac {\partial u_{x}}{\partial y}}\gg {\frac {\partial u_{x}}{\partial x}}$

${\frac {\partial p}{\partial x}}=0$

e semplificare l'equazione, dove resta solo il seguente termine

$0={\frac {\partial }{\partial y}}\left(\mu {\frac {\partial u_{x}}{\partial y}}-\rho u_{x}'u_{y}'\right).$

Nel primo termine all'interno della parentesi tonda è immediato riconoscere lo sforzo viscoso mentre il secondo termine può essere interpretato come lo sforzo turbolento. Sempre dalla stessa equazione si evince che lo sforzo totale è la somma dei due sforzi e che questo si mantiene costante lungo l'asse verticale y. Mentre questa approssimazione si verifica sperimantalmente in prossimità della parete essa appare inaccettabile in prossimità del limite esterno dello strato limite, dove gli sforzi tendono asintoticamente a 0.

Lo strato limite viene quindi suddiviso in tre regioni che vengono studiate separatamente:

1. regione interna
2. buffer
3. regione esterna.

Regione interna

La regione interna risulta essere la più importante sia per i bilanci delle forze che per i bilanci energetici. Inoltre in questa regione prevalgono gli sforzi viscosi, mentre quelli turbolenti possono essere trascurati. Lo sforzo agente risulta quindi essere

$\tau =\mu {\frac {\partial u_{x}}{\partial y}}.$

Dividendo entrambi i membri per ρ e ricordando che τ diviso ρ ha le dimensioni di una velocità al quadrato è possibile porre

${\frac {\tau }{\rho }}=u_{0}^{2}$

dove u₀ è una velocità fittizia. La relazione precedente si riscriverà quindi come

$u_{0}^{2}={\mathcal {D}}_{v}{\frac {\partial u_{x}}{\partial y}}$

che si può esprimere come

$\left({\frac {u_{0}}{{\mathcal {D}}_{v}}}\right)={\frac {\partial }{\partial y}}\left({\frac {u}{u_{0}}}\right)$

dove ${\mathcal {D}}_{v}$ è la viscosità cinematica. Infine, rendendo adimensionale la precedente equazione

$d\left({\frac {u}{u_{0}}}\right)=d\left({\frac {u_{0}y}{{\mathcal {D}}_{v}}}\right)$

ed esprimendo nelle nuove variabili adimensionali anche dette variabili interne,

${\frac {u}{u_{0}}}=u^{+}$

${\frac {u_{0}y}{{\mathcal {D}}_{v}}}=y^{+}$

la legge per lo strato interno si riassume come

$y^{+}=u^{+}.\$

Da osservazioni sperimentali si è osservato che tale legge è approssimata per valori di y⁺ < 5.

Buffer

Per valori di y⁺ compresi tra 5 e 30 l'effetto degli sforzi turbolenti non è più trascurabile e gli sforzi viscosi sono ugualmente importanti. Questa regione è chiamata buffer o semplicemente cuscinetto. Come nel caso della regione interna anche in questa regione ci sono teorie matematiche che permettono di ricavare la distribuzione delle velocità, tuttavia questa regione è caratterizzata da vortici intermittenti e da moti non stazionari, fatto che rende la trattazione relativamente complessa.

Pertanto nel corso degli anni sono state proposte diverse distribuzioni di velocità, molte delle quali dedotte sperimentalmente. Una delle distribuzioni utilizzate più di frequente è quella proposta da Spalding nel 1961 che normalmente viene scritta in forma implicita

$y^{+}=u^{+}+e^{-kB}\left[e^{(ku^{+})}-1-ku^{+}-{\frac {(ku^{+})^{2}}{2}}-{\frac {(ku^{+})^{3}}{6}}\right]$

dove il valore della costante B è fissato con 5,5, mentre per la costante di von Karman si utilizza il valore 0,4.

Regione esterna

Al contrario di quanto visto per la regione interna, nella regione esterna prevalgono gli sforzi di tipo turbolento. Essa si estende per valori della y⁺ compresi tra 35 e 50 ed è pertanto la regione più grande dello strato limite. Gli sforzi in questa regione possono essere indicati tramite la seguente espressione

${\frac {\tau }{\rho }}=-u_{x}'u_{y}'.$

Osservato quindi che in virtù del fatto che il moto è unidirezionale si osserva che il tensore degli sforzi è composto da un singolo termine. Pertanto appare lecito applicare all'equazione il modello diffusivo di Boussinesq riscrivendo il termine a destra dell'equazione precedente come segue

$u_{x}'u_{y}'=-\varepsilon \left({\frac {\partial u_{i}}{\partial x_{j}}}+{\frac {\partial u_{j}}{\partial x_{i}}}\right)$

che sempre in virtù del fatto che il moto sia unidemensionale si riscrive semplicemente come

$u_{x}'u_{y}'=-\varepsilon \left({\frac {\partial u_{x}}{\partial y}}\right).$

Inserendo quindi l'espressione ricavata nella prima equazione di partenza nella quale si erano trascurati il termine dello sforzo viscoso si ottiene

${\frac {\tau }{\rho }}=-\varepsilon \left({\frac {\partial u}{\partial y}}\right).$

Si osservi inoltre che in assenza di componenti della velocità in direzione y si è trascurato il pedice della velocità x. A questo punto per integrare l'equazione è necessario introdurre un opportuno modello di chiusura. Si può quindi ad esempio utilizzare il modello di Prandtl secondo il quale si ha

$\varepsilon =l^{2}\left({\frac {\partial u}{\partial y}}\right)$

con

$l^{2}=(ky)^{2}.\$

l'equazione di partenza diviene

$(ky)^{2}\left({\frac {\partial u}{\partial y}}\right)^{2}=u_{0}^{2}$

dove k indica la costante di von Karman il cui valore è compreso tra 0,4 e 0,41. Appurato che tutte le grandezze prese in esame sono positive e dopo aver estratto la radice quadrata positiva si ottiene

$(ky)\left({\frac {\partial u}{\partial y}}\right)=u_{0}.$

Rimanipolando l'equazione

$d\left({\frac {u}{u_{0}}}\right)={\frac {1}{k}}{\frac {dy}{y}}$

e ricordando la definizione data in precedenza per y^+ e u^+ si ottiene

$d\left({\frac {u}{u_{0}}}\right)={\frac {1}{k}}{\frac {d\left({\frac {u_{0}y}{{\mathcal {D}}_{v}}}\right)}{\left({\frac {u_{0}y}{{\mathcal {D}}_{v}}}\right)}}.$

Integrando infine l'espressione precedente si giunge all'equazione che descrive l'andamento delle velocità nello strato esterno

$u^{+}={\frac {1}{k}}\ln(y^{+})+B.$

Per calcolare la costante di integrazione B si dovrebbe conoscere il valore della velocità in una delle due estremità della regione esterna. Tuttavia nella maggior parte delle applicazioni non si conosce né la velocità del Buffer né quella del moto all'esterno dello strato limite e non è quindi possibile calcolare il valore della costante. Dal punto di vista applicativo però si potrebbe immaginare di trascurare il buffer ed immaginare di collegare la regione interna con quella esterna, supponendo che esista una transizione diretta dal sottostrato di turbolenza di parete al sottostrato viscoso. Le risultanze sperimentali di Nikuradse indicano che questo passaggio avviene per valori di y⁺ = 11,6, per cui questo valore è stato fissato con 5,5. Studi più recenti condotti nel 1986 effettuati da Coles e Hirst hanno invece corretto il valore calcolato da Nikuradse fissandolo a 5.

Approssimazione dello strato limite turbolento con un'unica legge di potenza

Per molte applicazioni, come il calcolo della resistenza, l'utilizzo delle tre distribuzioni di velocità può risultare poco agevole. In passato sono quindi state proposte diverse leggi, come approssimare con un'unica legge la distribuzione di velocità. Una tra le più note è la legge proposta da Schlichting nel 1968, che propone la seguente espressione

$u^{+}=8{,}74\,\left(y^{+}\right)^{\frac {1}{7}}.$

Si noti che anche questa legge è scritta in variabili interne. Inoltre la costante dell'esponente dovrebbe dipendere dal numero di Reynolds. Pertanto per la formula proposta da Schlichting sono stati consigliati dei correttivi. Hinze ad esempio propose per il valore dell'esponente di utilizzare il valore 8,3, mentre Clauser propose un valore compreso tra 0,3 e 0,1. Va inoltre precisato che la legge può essere anche estesa alla regione esterna.

Voci correlate

• Aerodinamica
• Fluidodinamica
• Strato limite
• Controllo dello strato limite

Altri progetti

Altri progetti

• Wikimedia Commons

• Wikimedia Commons contiene immagini o altri file su strato limite di quantità di moto

Collegamenti esterni

• Strato limite: principali caratteristiche ed effetti (PDF), su diiar.polimi.it. URL consultato il 22 giugno 2009 (archiviato dall'url originale il 29 dicembre 2009).

Controllo di autorità	J9U (EN, HE) 987007283990105171 · NDL (EN, JA) 00567209