Tourbillon de Burgers

Tourbillon de Burgers. Lignes de courant (rouge) et toubillon (vert).

En dynamique des fluides, le tourbillon de Burgers ou toubillon de Burgers–Rott est une solution exacte des équations de Navier-Stokes dans un écoulement incompressible. Il est ainsi nommé d'après Johannes Martinus Burgers et Nicholas Rott. Il décrit un écoulement stationnaire auto-similaire en point d'arrêt en symétrie cylindrique où l'écoulement radial tend à concentrer le toubillon dans une colonne étroite autour de l'axe de symétrie tandis qu'un étirement axial provoque l'augmentation du diamètre. En même temps la dissipation visqueuse tend à étendre le tourbillon. Le tourbillon de Burgers stationnaire apparaît lorsque les trois effets sont en équilibre.

Le tourbillon de Burgers, en plus de servir d'illustration du mécanisme d'étirement tourbillonnaire peut décrire des écoulements tels que des tornades, où la vorticité est alimentée par la convection.

Champ d'écoulement

Vitesse azimutale et composante $z$ de vorticité dans un tourbillon de Burgers.

L'écoulement est décrit dans des coordonnées cylindriques $(r,\theta ,z)$. En supposant une symétrie axiale (pas de dépendance à $\theta$), le champ d'écoulement associé à l'écoulement axisymétrique de point d'arrêt s'écrit :

$v_r=-\alpha r$

$v_z=2\alpha z$

$v_{\theta }=\frac{\mathrm{\Gamma }}{2\pi r}g(r)$

où le taux de déformation $\alpha >0$ et la circulation de la vitesse $\mathrm{\Gamma }>0$ sont donnés (paramètres du problème). L'écoulement satisfait l'équation de continuité par les deux premières équations ci-dessus. L'équation de conservation de quantité de mouvement azimutale se réduit alors à :

$r\frac{d^{2}g}{d{r}^2}+(\frac{\alpha r^2}{\nu }-1)\frac{dg}{dr}=0$

où $\nu$ est la viscosité cinématique du fluide. L'équation est intégrée avec la condition $g(\mathrm{\infty })=1$ de sorte qu'à l'infini la solution se comporte comme un écoulement potentiel. Le choix $g(0)=0$ garantit $v_{\theta }=0$ sur l'axe. La solution est :

$g=1-\exp (-\frac{\alpha r^2}{2\nu })$

L'équation de vorticité ne donne qu'une composante non triviale dans la direction $z$, donnée par :

${\omega }_z=\frac{\alpha \mathrm{\Gamma }}{2\pi \nu }\exp (-\frac{\alpha r^2}{2\nu })$

Intuitivement, l'écoulement peut être compris en regardant les trois termes de l'équation de vorticité pour ${\omega }_z$ :

$-\alpha r\frac{d{\omega }_z}{dr}=2\alpha {\omega }_z+\frac{\nu }{r}\frac{d}{dr}\left(r\frac{d{\omega }_z}{dr}\right)$

Le premier terme du côté droit de l'équation ci-dessus correspond à l'étirement du tourbillon qui intensifie la vorticité du noyau du tourbillon en raison de la composante de vitesse axiale $v_z=2\alpha z$. Cette vorticité tente de se diffuser radialement vers l'extérieur en raison du deuxième terme du côté droit, mais est empêchée par la convection de vorticité radiale due à $v_r=-\alpha r$ présent du côté gauche de l'équation ci-dessus. L'équilibre des trois mécanismes conduit à une solution stable.

L'une des propriétés importantes du tourbillon de Burgers est que le taux de dissipation visqueuse total $\mathrm{\Phi }$ par unité de longueur axiale est indépendant de la viscosité, indiquant que la dissipation est non nulle même dans la limite $\nu \to 0$. Pour cette raison, il sert de candidat approprié pour la modélisation et la compréhension des tubes tourbillonnaires étirés observés dans les écoulements turbulents. Le taux de dissipation total par unité de longueur axiale est égal à l'enstrophie totale par unité de longueur, qui est donnée par :

$\mathrm{\Phi }=\nu \underset{0}{\overset{2\pi }{\int }}\underset{0}{\overset{\mathrm{\infty }}{\int }}{\omega }_z^{2}rdrd\theta =\frac{\alpha {\mathrm{\Gamma }}^2}{4\pi }$

Évolution instationnaire vers le tourbillon de Burgers

Une solution exacte des équations de Navier-Stokes dépendantes du temps pour une fonction arbitraire $\alpha =\alpha (t)$ peut être obtenue. En particulier, lorsque $\alpha$ est constant, le champ de vorticité $\omega (r,t)$ avec une distribution initiale arbitraire $\mathrm{\Omega }(r)=\omega (r,0)$ est donné par :

$\omega (r,t)=\frac{\alpha }{2\pi \nu (1-e^{-2\alpha t})}\iint \mathrm{\Omega }\left(\sqrt{{\xi }_1^2+{\eta }_1^2}\right)\exp \left[\frac{-(x-{\xi }_1{e}^{-\alpha t}{)}^2-(y-{\eta }_1{e}^{-\alpha t}{)}^2}{(2\nu /\alpha )(1-e^{-2\alpha t})}\right]d{\xi }_{1}d{\eta }_1$

Lorsque $t\to \mathrm{\infty }$, le comportement asymptotique est donné par :

$\omega (r,t)=\frac{\alpha }{2\pi \nu }\exp (-\frac{\alpha r^2}{2\nu })[\mathrm{\Gamma }+e^{-2\alpha t}\underset{0}{\overset{\mathrm{\infty }}{\int }}\mathrm{\Omega }(s)(\frac{\alpha r^2}{2\nu }-1)(\frac{\alpha s^2}{2\nu }-1)2\pi sds+O(e^{-4\alpha t})],\mathrm{\Gamma }=\underset{0}{\overset{\mathrm{\infty }}{\int }}\mathrm{\Omega }(s)2\pi sds$

Ainsi, à condition que $\mathrm{\Gamma }\ne 0$, une distribution de vorticité arbitraire tend vers le tourbillon de Burgers. Si $\mathrm{\Gamma }=0$, par exemple dans le cas où la condition initiale est composée de deux tourbillons égaux et opposés, alors le premier terme est nul et le second terme implique que la vorticité décroît jusqu'à zéro lorsque $t\to \mathrm{\infty }$.

Couche tourbillonnaire de Burgers

Vitesse de la composante $y$ et vorticité de la composante $z$ dans une couche toubillonnaire de Burgers.

La couche tourbillonnaire ou nappe de tourbillon de Burgers est une couche de cisaillement contrainte qui est un analogue bidimensionnel plan du tourbillon de Burgers. C'est également une solution exacte des équations de Navier–Stokes, décrites pour la première fois par Albert Townsend en 1951. Le champ de vitesse $(v_x,v_y,v_z)$ exprimé dans les coordonnées cartésiennes est :

$v_x=-\alpha x$

$v_z=\alpha z$

$v_y=U\mathrm{erf}\left(\frac{\sqrt{\alpha }x}{\sqrt{2\nu }}\right)$

où $\alpha >0$ est le taux de déformation, $v_y(+\mathrm{\infty })=U$ et $v_y(-\mathrm{\infty })=-U$. La valeur $2U$ est interprétée comme la force de la nappe tourbillonnaire. L'équation de vorticité ne donne qu'une composante non triviale dans la direction $z$ :

${\omega }_z=2U\sqrt{\frac{\alpha }{2\pi \nu }}\exp (-\frac{\alpha x^2}{2\nu })$

K. N. Beronov et S. Kida ont montré que la nappe tourbillonnaire de Burgers est instable subissant ainsi d'abord une instabilité de Kelvin-Helmholtz suivie d'autres instabilitésModèle:, et peut être en transition vers des tourbillon de Kerr–Dold à des nombres de Reynolds modérément élevés, mais devenant pleinement turbulente pour des nombres de Reynolds élevés.

Tourbillons de Burgers non axisymétriques

Les tourbillons de Burgers non axisymétriques apparaissent dans les écoulements contraints non axisymétriques. Leur existence pour les petits nombres de Reynolds $Re=\mathrm{\Gamma }/(2\pi \nu )$ a été établie par A. C. Robinson et Philip Saffman en 1984 tandis que Keith Moffatt, S. Kida et K. Ohkitani ont développé la théorie pour $Re\gg 1$ en 1994. La structure des tourbillons de Burgers non axisymétriques pour des valeurs arbitraires du nombre de Reynolds tourbillonnaire peut être discutée par le biais d'intégrations numériques. Le champ de vitesse prend la forme :

$v_x=-\alpha x+u(x,y)$

$v_y=-\beta y+v(x,y)$

$v_z=\gamma z$

avec $\gamma =\alpha +\beta$. Sans perte de généralité, on suppose $\alpha >0$ et $\gamma >0$. La section transversale du tourbillon se situe dans le plan $xy$, ce qui donne une composante de vorticité non nulle dans la direction $z$ :

${\omega }_z=\frac{\partial u}{\partial y}-\frac{\partial v}{\partial x}$

Le tourbillon de Burgers axisymétrique est récupéré lorsque $\alpha =\beta =\gamma /2$ tandis que la couche de tourbillon de Burgers est récupérée lorsque $\alpha =\gamma$ et $\beta =0$.

Tourbillon de Burgers avec source centrale

La solution explicite en géométrie cylindrique avec source sur l'axe z de rayon $r_s\ge 0$ a été résolue par Prabakaran Rajamanickam et Adam D. Weiss. La solution est exprimée dans le système de coordonnées cylindriques s'écrit :

$v_r=-\alpha (r-\frac{r_s^2}{r})$

$v_z=2\alpha z$

$v_{\theta }=\frac{\mathrm{\Gamma }}{2\pi r}P(1+\frac{\alpha r_s^2}{2\nu },\frac{\alpha r^2}{2\nu })$

où $\alpha >0$ est le taux de déformation, $\mathrm{\Gamma }>0$ est la circulation de la vitesse et $P$ est la fonction gamma incomplète. Cette solution n'est autre que le tourbillon de Burgers en présence d'un source linéaire avec une intensité de source $Q=2\pi \alpha r_s^2$. L'équation de vorticité ne donne qu'une composante non triviale dans la direction $z$, donnée par :

${\omega }_z=\frac{\alpha \mathrm{\Gamma }}{2\pi \nu \tilde{\mathrm{\Gamma }}(1+\alpha r_s^2/2\nu )}{\left(\frac{\alpha r^2}{2\nu }\right)}^{\alpha r_s^2/2\nu }\exp (-\frac{\alpha r^2}{2\nu })$

où $\tilde{\mathrm{\Gamma }}$ dans l'expression ci-dessus est la fonction gamma. Lorsque $r_s\to 0$, la solution se réduit à la solution du toubillon de Burgers et lorsque $r_s\to \mathrm{\infty }$ la solution devient la solution de la couche tourbillonnaire de Burgers.