Utilisateur:


Écoulement turbulent à travers des tubes

Storyboard

Si le nombre de Reynolds dépasse 2000, l'écoulement dans un tube devient toujours instable et finit par devenir complètement turbulent. Par conséquent, il n'est plus possible d'utiliser l'approximation d'écoulement visqueux laminaire qui donne lieu à la loi de Hagen-Poiseuille, et un modèle alternatif est nécessaire.

Le modèle qui décrit un écoulement où la viscosité est négligeable est celui qui donne naissance à l'équation de Bernoulli. Cependant, ce modèle suppose que la densité d'énergie soit conservée. Une alternative consiste à supposer que les turbulences conduisent à un mélange de telle manière que la densité d'énergie ne soit pas conservée mais reste constante. Dans ce cas, l'écoulement peut être modélisé à l'aide d'une équation similaire à celle de Bernoulli, mais avec une correction pour prendre en compte l'homogénéisation due aux effets de mélange.

>Modèle

ID:(1970, 0)



Modélisation des turbulences

Définition

L'écoulement laminaire est décrit par des "feuilles" qui se déplacent à différentes vitesses de manière coordonnée. En revanche, dans l'écoulement turbulent, ces "feuilles" n'existent pas. Les éléments du fluide sont déviés, peuvent changer de direction et participer à des mouvements circulaires, souvent de manière chaotique.

La première conséquence de cela est que les paramètres du fluide ont tendance à se mélanger, provoquant la disparition des différences de vitesse et l'émergence d'un type de vitesse moyenne. En moyennant le mouvement, des motifs structurés émergent, mais ils ne proviennent plus des éléments individuels du fluide, mais plutôt d'un mélange temporel. En conséquence, un profil relativement constant similaire aux profils définis dans l'écoulement laminaire émerge à nouveau, mais dans ce cas, ce sont des valeurs moyennes et présentant moins de grands gradients, les rendant plus uniformes.

ID:(14525, 0)



Diagramme de Moody

Image

En 1944, Lewis Ferry Moody a mesuré le facteur de frottement de Darcy-Weisbach en fonction du nombre de Reynolds et de la rugosité relative de la paroi, ce qui a conduit à la création du diagramme suivant :

La rugosité relative peut être estimée en considérant la taille des irrégularités de surface (hauteur des saillies ou profondeur des rainures) par rapport au diamètre hydraulique.

Deux comportements distincts sont observés :

• Pour des nombres de Reynolds inférieurs à 2000, le facteur de frottement de Darcy-Weisbach dépend uniquement du nombre de Reynolds, suivant une relation de $64/Re$. Cela correspond au régime d'écoulement laminaire.

• Pour des nombres de Reynolds supérieurs à 2000, un comportement est observé qui dépend à la fois du nombre de Reynolds et de la rugosité relative de la surface du tube.

ID:(14528, 0)



Écoulement à la limite laminaire

Noter

Si nous remplaçons le facteur de frottement de Darcy-Weisbach dans la limite laminaire, donné par

$ f_D = \displaystyle\frac{64}{ Re }$



dans l'équation de Darcy-Weisbach, exprimée comme

$ \Delta p = \displaystyle\frac{ \Delta L }{ D_H } f_D \displaystyle\frac{1}{2} \rho_w j_s ^2 $



et utilisons la définition du nombre de Reynolds $Re$, nous pouvons démontrer que l'écoulement est régi par

$ J_V =-\displaystyle\frac{ \pi R ^4}{8 \eta }\displaystyle\frac{ \Delta p }{ \Delta L }$

ce qui correspond à l'équation de Hagen-Poiseuille.

ID:(14530, 0)



Profil de vitesse d'écoulement turbulent dans le tube

Citation

Le profil de vitesse de l'écoulement turbulent à travers un tube présente deux zones distinctes en fonction de la distance à la surface ($z$), où $\delta$ représente l'épaisseur de la couche limite. Dans la région proche de la surface ($z < \delta$), l'écoulement est essentiellement laminaire, tandis que dans la région plus éloignée de la surface ($z > \delta$), l'écoulement devient turbulent.



Dans la plage laminaire, la vitesse est proportionnelle à la distance normalisée :

$u^+ = y^+$



Cette relation est similaire au profil de vitesse de Hagen-Poiseuille près de la paroi et représente une approximation linéaire.

Dans la plage turbulente, le profil de vitesse normalisée prend une forme logarithmique :

$u^+ = \displaystyle\frac{1}{\kappa} \ln\left(\displaystyle\frac{y^+}{y_0}\right)$

Ici, $\kappa$ est la constante de Karman (environ $0,41$), et $y_0\sim 1/8$ est la distance normalisée à laquelle la vitesse serait nulle. Cependant, il est important de noter que la largeur de la couche laminée, comme indiqué sur le graphique, est d'environ 7,072 fois plus grande que $y_0$.

ID:(14536, 0)



Écoulement turbulent à travers des tubes

Description

Si le nombre de Reynolds dépasse 2000, l'écoulement dans un tube devient toujours instable et finit par devenir complètement turbulent. Par conséquent, il n'est plus possible d'utiliser l'approximation d'écoulement visqueux laminaire qui donne lieu à la loi de Hagen-Poiseuille, et un modèle alternatif est nécessaire. Le modèle qui décrit un écoulement où la viscosité est négligeable est celui qui donne naissance à l'équation de Bernoulli. Cependant, ce modèle suppose que la densité d'énergie soit conservée. Une alternative consiste à supposer que les turbulences conduisent à un mélange de telle manière que la densité d'énergie ne soit pas conservée mais reste constante. Dans ce cas, l'écoulement peut être modélisé à l'aide d'une équation similaire à celle de Bernoulli, mais avec une correction pour prendre en compte l'homogénéisation due aux effets de mélange.

Variables

Symbole
Texte
Variable
Valeur
Unités
Calculer
Valor MKS
Unités MKS
$S$
S
Coupe des pores
m^2
$\rho$
rho
Densité
kg/m^3
$j_s$
j_s
Densité de flux
m/s
$\rho_w$
rho_w
Densité du liquide
kg/m^3
$D_H$
D_H
Diamètre hydrodynamique
m
$\Delta p$
Dp
Différence de pression
Pa
$R$
R
Dimension typique du système
m
$f_D$
f_D
Facteur de frottement Darcy-Weisbach
-
$\Delta h$
Dh
Hauteur de la colonne de liquide
m
$\Delta h_d$
Dh_d
Hauteur pour compenser la perte d'énergie
m
$\epsilon$
e
Inégalité
m
$Re$
Re
Le numéro de Reynold
-
$\Delta L$
DL
Longueur du tube
m
$P$
P
Périmètre
m
$P_H$
P_H
Périmètre hydrodynamique
m
$h$
h
Profondeur dans un tube non rempli
m
$R$
R
Rayon du tube
m
$R_H$
R_H
Rayon hydraulique
m
$\eta$
eta
Viscosité
Pa s
$J_V$
J_V
Volumique flux
m^3/s

Calculs


D'abord, sélectionnez l'équation:   à ,  puis, sélectionnez la variable:   à 

Symbole
Équation
Résolu
Traduit

Calculs

Symbole
Équation
Résolu
Traduit

 Variable   Donnée   Calculer   Cible :   Équation   À utiliser



Équations

S'il existe a différence de pression ($\Delta p$) entre deux points, comme le d termine l' quation :

$ dp = p - p_0 $



nous pouvons utiliser a pression de la colonne d'eau ($p$), qui est d finie comme suit :

$ p_t = p_0 + \rho_w g h $



Cela donne :

$\Delta p=p_2-p_1=p_0+\rho_wh_2g-p_0-\rho_wh_1g=\rho_w(h_2-h_1)g$



Comme a différence de hauteur ($\Delta h$) est d finie comme suit :

$ \Delta h = h_2 - h_1 $



a différence de pression ($\Delta p$) peut tre exprim e comme suit :

$ \Delta p = \rho_w g \Delta h $

(ID 4345)

Le flux est d fini comme le volume le élément de volume ($\Delta V$) divis par le temps le temps écoulé ($\Delta t$), ce qui est exprim dans l' quation suivante :

$ J_V =\displaystyle\frac{ \Delta V }{ \Delta t }$



et le volume est gal la section transversale a section de tube ($S$) multipli e par la distance parcourue le élément tubulaire ($\Delta s$) :

$ \Delta V = S \Delta s $



tant donn que la distance parcourue le élément tubulaire ($\Delta s$) par unit de temps le temps écoulé ($\Delta t$) correspond la vitesse, elle est repr sent e par :

$ j_s =\displaystyle\frac{ \Delta s }{ \Delta t }$



Ainsi, le flux est une densité de flux ($j_s$), qui est calcul l'aide de :

$ j_s = \displaystyle\frac{ J_V }{ S }$

(ID 4349)

La solution originale de S.E. Haaland est la suivante:

$\displaystyle\frac{1}{\sqrt{f_D}}=-1.8\log\left(\left(\displaystyle\frac{\epsilon/R_H}{12}\right)^{1.11} + \displaystyle\frac{6.9}{Re}\right)$



Elle peut tre r organis e pour obtenir l'expression du facteur de friction de Darcy-Weisbach $f_D$ comme suit :

$ f_D = \displaystyle\frac{25}{81 \log\left(\left(\displaystyle\frac{ \epsilon }{12 R_H }\right)^{1.11} + \displaystyle\frac{6.9}{ Re }\right)^2}$

(ID 14534)

La solution originale de S.E. Haaland est la suivante:

$\displaystyle\frac{1}{\sqrt{f_D}}=-1.8\log\left(\left(\displaystyle\frac{\epsilon/D_H}{3.7}\right)^{1.11} + \displaystyle\frac{6.9}{Re}\right)$



Elle peut tre r organis e pour obtenir l'expression du facteur de friction de Darcy-Weisbach $f_D$ comme suit :

$ f_D = \displaystyle\frac{25}{81 \log\left(\left(\displaystyle\frac{ \epsilon }{3.7 D_H }\right)^{1.11} + \displaystyle\frac{6.9}{ Re }\right)^2}$

(ID 14535)

A coupe ou surface ($S$) du tube contenant le liquide peut tre exprim en fonction de a profondeur dans un tube non rempli ($h$) en int grant sur le rayon jusqu' Le rayon du tube ($R$) :

$S = 2\displaystyle\int_0^h dz \sqrt{2Rz - z^2}=\displaystyle\frac{1}{2}(R-h)\sqrt{2Rh - h^2}+\displaystyle\frac{1}{2}R^2\arcsin\left(\displaystyle\frac{1-h}{R}\right)$



Si nous d veloppons cette expression en fonction du facteur $h/R$ dans la limite $h\ll R$, nous obtenons, au premier ordre :

$S = \sqrt{\displaystyle\frac{2^3}{3^2} R h ^3}$



Si nous r solvons pour la profondeur, nous obtenons finalement:

$ h = \left(\displaystyle\frac{3^2 S^2 }{2^3 R }\right)^{1/3}$

(ID 14541)

tant donn que l'angle entre le rayon au bord de la surface du liquide et la verticale peut tre calcul avec a profondeur dans un tube non rempli ($h$) et le rayon du tube ($R$) via :

$\phi =\arccos\left(1-\displaystyle\frac{h}{R}\right)$



L'arc correspondant est $R\phi$, donc l'arc total est :

$2 R \arccos\left(1-\displaystyle\frac{h}{R}\right)$



De mani re analogue, la moiti de la surface peut tre d termin e en utilisant le th or me de Pythagore, ce qui donne :

$\sqrt{2Rh - h^2}$



Ainsi, le périmètre hydrodynamique ($P_H$) s'exprime comme suit :

$P_H = 2\sqrt{2Rh - h^2}+2R\arcsin\left(1-\displaystyle\frac{h}{R}\right)$



Dans le cas d'une petite hauteur, o $h\ll R$, cette expression peut tre d velopp e, aboutissant :

$ P_H = \sqrt{2^5 R h }$

(ID 14542)


Exemples

Lorsqu'on mod lise l' coulement dans un tube en supposant que la densit d' nergie est conserv e, on obtient l' quation de Bernoulli, qui d crit l' coulement l'aide de a différence de pression ($\Delta p$) en fonction de a densité du liquide ($\rho_w$), a vitesse moyenne du fluide ($v$) et a différence de vitesse entre les surfaces ($\Delta v$) :

$ \Delta p = - \rho \bar{v} \Delta v $



Dans le cas dun coulement turbulent, le processus de m lange agit comme une friction qui r duit le gradient de vitesse pr sent dans l coulement laminaire entre le centre et les parois du tube. Si lon suppose que cet effet de m lange peut tre mod lis par un simple facteur de correction, on arrive empiriquement l quation de Darcy-Weisbach, qui implique le facteur de frottement Darcy-Weisbach ($f_D$), le longueur du tube ($\Delta L$) et a diamètre hydrodynamique ($D_H$) :

$ \Delta p = \displaystyle\frac{ \Delta L }{ D_H } f_D \displaystyle\frac{1}{2} \rho_w j_s ^2 $



Le facteur le facteur de frottement Darcy-Weisbach ($f_D$) a t d termin empiriquement pour diff rentes conditions d' coulement et est exprim en fonction de le le numéro de Reynold ($Re$).

(ID 15893)

En 1944, Lewis Ferry Moody a mesur le facteur de frottement de Darcy-Weisbach en fonction du nombre de Reynolds et de la rugosit relative de la paroi, ce qui a conduit la cr ation du diagramme suivant :

La rugosit relative peut tre estim e en consid rant la taille des irr gularit s de surface (hauteur des saillies ou profondeur des rainures) par rapport au diam tre hydraulique.

Deux comportements distincts sont observ s :

• Pour des nombres de Reynolds inf rieurs 2000, le facteur de frottement de Darcy-Weisbach d pend uniquement du nombre de Reynolds, suivant une relation de $64/Re$. Cela correspond au r gime d' coulement laminaire.

• Pour des nombres de Reynolds sup rieurs 2000, un comportement est observ qui d pend la fois du nombre de Reynolds et de la rugosit relative de la surface du tube.

(ID 14528)

Dans le contexte de l' quation de Darcy-Weisbach, on utilise une diamètre hydrodynamique ($D_H$), qui correspond une g n ralisation du diam tre traditionnel d'un cercle. Cela permet de consid rer une section non circulaire et de calculer un diam tre quivalent bas sur la surface de a section de tube ($S$) et son le périmètre ($P$) l'aide de la formule suivante :

$ D_H = \displaystyle\frac{ 4 S }{ P }$



Pour une section circulaire, nous obtenons le diam tre traditionnel d'un cercle de la mani re suivante :

$D_H = \displaystyle\frac{4 S}{P} = \displaystyle\frac{4 \pi R^2}{2 \pi R} = 2R$

(ID 15894)

Dans le contexte du facteur de friction de Darcy-Weisbach, on utilise un rayon hydraulique ($R_H$), qui est une g n ralisation du rayon traditionnel d'un cercle. Ainsi, il est possible de calculer un diam tre bas sur la surface de a section de tube ($S$) et son p rim tre en contact avec le périmètre hydrodynamique ($P_H$) l'aide de la formule suivante :

$ R_H = \displaystyle\frac{ S }{ P_H }$



Pour une section circulaire, nous pouvons obtenir le rayon hydraulique traditionnel d'un cercle de la mani re suivante :

$R_H = \displaystyle\frac{S}{P} = \displaystyle\frac{\pi R^2}{2 \pi R} = \displaystyle\frac{1}{2} R$

(ID 15895)

Dans un tube cylindrique, la profondeur est li e au d bit de la mani re suivante :



En int grant a coupe ou surface ($S$), nous pouvons calculer comment la surface varie en fonction de a profondeur dans un tube non rempli ($h$) travers l'int grale sur le rayon jusqu' Le rayon du tube ($R$), ce qui donne :

$S = 2\displaystyle\int_0^h dz \sqrt{2Rz - z^2}=\displaystyle\frac{1}{2}(R-h)\sqrt{2Rh - h^2}+\displaystyle\frac{1}{2}R^2\arcsin\left(\displaystyle\frac{1-h}{R}\right)$



Pour les petits d bits, o la profondeur est significativement inf rieure au rayon, la relation entre la section transversale et la profondeur se simplifie consid rablement. En r solvant l' quation de profondeur, nous obtenons :

$ h = \left(\displaystyle\frac{3^2 S^2 }{2^3 R }\right)^{1/3}$

(ID 15896)

Le périmètre hydrodynamique ($P_H$) dans un tube partiellement rempli correspond aux bords de la section en contact avec le liquide, cest- -dire larc qui touche la fois la paroi du tube et la surface :



Ainsi, on peut lexprimer en fonction de le rayon du tube ($R$) et a profondeur dans un tube non rempli ($h$) de la mani re suivante :

$P_H = 2 R \arccos\left(1-\displaystyle\frac{h}{R}\right) + 2\sqrt{2Rh-h^2}$



Pour les faibles d bits, lorsque la profondeur est nettement inf rieure au rayon, cela simplifie la relation entre la section transversale et la profondeur :

$ P_H = \sqrt{2^5 R h }$

(ID 15897)

Si nous rempla ons le facteur de frottement de Darcy-Weisbach dans la limite laminaire, donn par

$ f_D = \displaystyle\frac{64}{ Re }$



dans l' quation de Darcy-Weisbach, exprim e comme

$ \Delta p = \displaystyle\frac{ \Delta L }{ D_H } f_D \displaystyle\frac{1}{2} \rho_w j_s ^2 $



et utilisons la d finition du nombre de Reynolds $Re$, nous pouvons d montrer que l' coulement est r gi par

$ J_V =-\displaystyle\frac{ \pi R ^4}{8 \eta }\displaystyle\frac{ \Delta p }{ \Delta L }$

ce qui correspond l' quation de Hagen-Poiseuille.

(ID 14530)

L' quation de Darcy-Weisbach permet de calculer a différence de pression ($\Delta p$) en fonction de a densité du liquide ($\rho_w$), a diamètre hydrodynamique ($D_H$), le facteur de frottement Darcy-Weisbach ($f_D$), le longueur du tube ($\Delta L$) et a vitesse moyenne du fluide ($v$) travers :

$ \Delta p = \displaystyle\frac{ \Delta L }{ D_H } f_D \displaystyle\frac{1}{2} \rho_w j_s ^2 $

(ID 14526)

Une densité de flux ($j_s$) peut tre exprim en termes de le volumique flux ($J_V$) l'aide de a coupe ou surface ($S$) par la formule suivante :

$ j_s = \displaystyle\frac{ J_V }{ S }$

(ID 4349)

A diamètre hydrodynamique ($D_H$) peut tre calcul partir de a section de tube ($S$) et le périmètre ($P$) via :

$ D_H = \displaystyle\frac{ 4 S }{ P }$

(ID 14527)

Le rayon hydraulique ($R_H$) peut tre calcul l'aide de a section de tube ($S$) et le périmètre hydrodynamique ($P_H$) via :

$ R_H = \displaystyle\frac{ S }{ P_H }$

(ID 14531)

A profondeur dans un tube non rempli ($h$) peut tre estim lorsque le niveau du liquide est bas, c'est- -dire lorsque ce niveau est bien inf rieur au rayon du tube, en fonction de a coupe ou surface ($S$) et le rayon du tube ($R$) :

$ h = \left(\displaystyle\frac{3^2 S^2 }{2^3 R }\right)^{1/3}$

(ID 14541)

Le périmètre hydrodynamique ($P_H$) peut tre estim partir de a profondeur dans un tube non rempli ($h$) et le rayon du tube ($R$) en utilisant :

$ P_H = \sqrt{2^5 R h }$

(ID 14542)

Le crit re cl pour d terminer si un milieu est lamin ou turbulent est le num ro de Reynolds, qui compare l' nergie associ e l'inertie celle associ e la viscosit . La premi re d pend de a densité ($\rho$), a vitesse moyenne du fluide ($v$) et a dimension typique du système ($R$), tandis que la seconde d pend de a viscosité ($\eta$), le d finissant ainsi :

$ Re =\displaystyle\frac{ \rho R v }{ \eta }$

(ID 3177)

Dans la limite de faibles le le numéro de Reynold ($Re$), le diagramme de Moody indique que le facteur de frottement Darcy-Weisbach ($f_D$) est gal :

$ f_D = \displaystyle\frac{64}{ Re }$

o $Re$ est le nombre de Reynolds. Ceci s'applique pour les nombres de Reynolds jusqu' 2000. Au-del de cette valeur, la rugosit de la paroi commence influencer la d stabilisation de l' coulement et la formation de turbulences.

(ID 14529)

L' quation de Colebrook-White pour le cas d'un tube ferm :

$\displaystyle\frac{1}{\sqrt{ f_D }}=-2\log\left(\displaystyle\frac{\epsilon}{3.7 D_H}+\displaystyle\frac{2.51}{ Re \sqrt{ f_D }}\right)$



est une quation implicite utilis e pour d terminer le coefficient de friction de Darcy-Weisbach ($f_D$). Pour r soudre cette quation, plusieurs approximations ont t d velopp es, variant en termes de complexit et de pr cision. L'une des approximations les plus efficaces couvrant une large gamme de nombres de Reynolds $Re$ est celle propos e par S.E. Haaland :

$ f_D = \displaystyle\frac{25}{81 \log\left(\left(\displaystyle\frac{ \epsilon }{3.7 D_H }\right)^{1.11} + \displaystyle\frac{6.9}{ Re }\right)^2}$

(ID 14535)

L' quation de Colebrook-White pour le cas d'un tube ferm :

$\displaystyle\frac{1}{\sqrt{ f_D }}=-2\log\left(\displaystyle\frac{\epsilon}{12 R_H}+\displaystyle\frac{2.51}{ Re \sqrt{ f_D }}\right)$



est une quation implicite utilis e pour d terminer le coefficient de friction de Darcy-Weisbach ($f_D$). Pour r soudre cette quation, plusieurs approximations ont t d velopp es, variant en termes de complexit et de pr cision. L'une des approximations les plus efficaces couvrant une large gamme de nombres de Reynolds $Re$ est celle propos e par S.E. Haaland :

$ f_D = \displaystyle\frac{25}{81 \log\left(\left(\displaystyle\frac{ \epsilon }{12 R_H }\right)^{1.11} + \displaystyle\frac{6.9}{ Re }\right)^2}$

(ID 14534)

La diff rence de hauteur, repr sent e par a différence de hauteur ($\Delta h$), implique que la pression dans les deux colonnes est diff rente. En particulier, a différence de pression ($\Delta p$) est une fonction de a densité du liquide ($\rho_w$), a accélération gravitationnelle ($g$), et a différence de hauteur ($\Delta h$), comme suit :

$ \Delta p = \rho_w g \Delta h $

(ID 4345)

L' quation de Darcy-Weisbach permet de calculer a différence de pression ($\Delta p$) en fonction de a densité du liquide ($\rho_w$), a diamètre hydrodynamique ($D_H$), le facteur de frottement Darcy-Weisbach ($f_D$), le longueur du tube ($\Delta L$) et a vitesse moyenne du fluide ($v$) travers :

$ \Delta h_d = \displaystyle\frac{ \Delta L }{ 2 g D_H } f_D j_s ^2 $

(ID 15958)


ID:(1970, 0)