Écoulement de Hiemenz

Écoulement en point d'arrêt.

L'écoulement de Hiemenz est un écoulement potentiel de point d'arrêt en symétrie plane dont la solution est analytique au sens où elle se ramène à la résolution d'une simple équation différentielle. Il a été décrit par Karl Hiemenz dans sa thèse à l'université de Göttingen en 1911 et a été étendu au cas de la symétrie de révolution par Fritz Homann en 1936.

Solution de Hiemenz

Écoulement de point d'arrêt plan. Fonction $\varphi (\eta )$

Écoulement de point d'arrêt plan : profil de vitesse.

Le problème posé est celui d'un écoulement irrotationnel impactant un cylindre perpendiculairement à celui-ci. Le potentiel est $\psi =axz$, $z$ étant l'axe portant l'écoulement amont. Les composantes de la vitesse dans le milieu amont sont :

$U=ax,V=-az$

Si $p_0$ est la pression au point d'arrêt, la pression en tout point est donnée par la conservation de quantité de mouvement :

$p_0-p=\frac{1}{2}\rho (U^2+V^2)=\frac{1}{2}\rho a^2(x^2+z^2)$

Au voisinage de la paroi il apparaît une couche limite et la solution générale est recherchée sous la forme suivante :

$u=x{f}^{\prime }(z),v=-f(z),p_0-p=\frac{1}{2}\rho a^2(x^2+F(z{)}^2)$

Ces équations satisfont par construction à l'équation de continuité. La conservation de quantité de mouvement conduit à :

$f{\text{'}}^2-f{f}^{\prime }=a^2+\nu f^{‴}$

$f{f}^{\prime }=\frac{1}{2}{a}^2{F}^{\prime }-\nu f^{″}$

où υ est la viscosité cinématique.

Les conditions aux limites sont :

$f(0)=0,f^{\prime }(0)=0,F(0)=0,\lim {\mathrm{\backslash limits}}_{z\to \mathrm{\infty }}{f}^{\prime }=a$

La première équation est indépendante et peut être transformée en posant :

$\eta =\alpha z,f(z)=A\varphi (\eta )$

Elle devient :

${\alpha }^2{A}^2(\varphi {\text{'}}^2-\varphi {\varphi }^{″})=a^2+\nu A{\alpha }^3{\varphi }^{‴}$

En posant :

${\alpha }^2{A}^2=a^2,\nu A{\alpha }^3=a^2\Rightarrow \eta =\sqrt{\frac{a}{\nu }}z,f(z)=\sqrt{a\nu }\varphi (\eta )$

Elle devient l'équation adimensionée :

${\varphi }^{‴}+\varphi {\varphi }^{″}-{\varphi }^2+1=0$

avec les conditions aux limites :

$\varphi (0)=0,{\varphi }^{\prime }(0)=0,\lim {\mathrm{\backslash limits}}_{\eta \to \mathrm{\infty }}{\varphi }^{\prime }=1$

La vitesse relative parallèle à la paroi (dans la « couche limite ») est indépendante de x :

$\frac{u}{U}=\frac{1}{a}{f}^{\prime }(z)={\varphi }^{\prime }(\eta )$

L'accélération pariétale est ${\varphi }^{″}(0)=1,2326$.

Si l'on prend pour définition de l'épaisseur de couche limite $\frac{u}{U}=0.99$, l'épaisseur de celle-ci (voir courbe) est donnée par $\delta \approx 2.4\sqrt{\frac{\nu }{a}}$.

Généralisations

La méthode a été généralisée à un point d'arrêt axisymétrique, à un jet en incidence et à des parois en mouvement, par exemple un cylindre en rotation.