Équation de Blasius

En physique, et plus particulièrement en mécanique des fluides, l'équation de Blasius décrit l'écoulement stationnaire et incompressible en 2 dimensions dans la couche limite se formant sur une plaque plane semi-infinie parallèle à l'écoulement. Plus précisément, le champ de vitesse tangentielle adimensionné est solution de cette équation :

$f^{″}(\theta )+\frac{1}{2}{f}^{\prime }(\theta )\underset{0}{\overset{\theta }{\int }}f(t)\mathrm{d}t=0$

Équation de mouvement dans la couche limite (Théorie de Prandtl)

Généralités

Phénomène de couche limite

Analysons l'écoulement bidimensionnel stationnaire dans le plan $(xOy)$ près d'une plaque plane placée en $y=0$, pour un écoulement extérieur potentiel $U(x)$ qu'on suppose parallèle à la paroi.

La dimension caractéristique dans la direction parallèle à l'écoulement est une longueur arbitraire $L$ que l'on suppose très grande devant la dimension caractéristique perpendiculaire à l'écoulement $\delta$ (épaisseur de la couche limite) :

$\delta \ll L$

Le raisonnement suivant est basé sur l'existence de ces 2 échelles.

Comparaison des termes

Nous considérons qu'un fluide de masse volumique $\rho$ et de viscosité dynamique $\eta$ (viscosité cinématique $\nu$) s'écoule le long de notre plaque plane. Partons de l'équation de Navier-Stokes en régime stationnaire avec la condition d'incompressibilité :

$\rho (\overrightarrow{v}.\overrightarrow{\mathrm{\nabla }})\overrightarrow{v}=-\overrightarrow{\mathrm{\nabla }}p+\eta \mathrm{\Delta }\overrightarrow{v}$

$\overrightarrow{\mathrm{\nabla }}.\overrightarrow{v}=0$

Avec le vecteur $\overrightarrow{v}=\left(\begin{array}{c}u\\ v\end{array}\right)$

Écrivons ce système sous sa forme projetée :

$\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0$

$u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=-\frac{1}{\rho }\frac{\partial p}{\partial x}+\nu (\frac{{\partial }^{2}u}{\partial x^2}+\frac{{\partial }^{2}u}{\partial y^2})$

$u\frac{\partial v}{\partial x}+v\frac{\partial v}{\partial y}=-\frac{1}{\rho }\frac{\partial p}{\partial y}+\nu (\frac{{\partial }^{2}v}{\partial x^2}+\frac{{\partial }^{2}v}{\partial y^2})$

Par une étude d'ordres de grandeurs, on peut montrer que les équations à résoudre peuvent se simplifier (en négligeant certains termes devant d'autres). Explicitons ce nouveau système :

$\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0$

$u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=\nu \frac{{\partial }^{2}u}{\partial y^2}$

Modèle:Démonstration

Équations adimensionnées

Il s'agit désormais d'adimensionner ces équations. Pour cela, nous allons réduire les variables judicieusement.

Pour plus de détails concernant les choix de réduction, ils s'inspirent de l'étude des ordres de grandeurs. Par exemple, lorsque l'on compare $u$ à $U$ il est nécessaire de comparer $v$ à $\frac{U}{\sqrt{{\mathrm{R}\mathrm{e}}_L}}$. Il en est de même pour les variables d'espace. Les nouvelles variables sont donc :

$u^{\prime }=\frac{u}{U}{v}^{\prime }=\frac{v\sqrt{{\mathrm{R}\mathrm{e}}_L}}{U}$

$x^{\prime }=\frac{x}{L}{y}^{\prime }=\frac{y\sqrt{{\mathrm{R}\mathrm{e}}_L}}{L}$

Avec ces nouvelles variables, le système devient simplement :

$\frac{\partial u^{\prime }}{\partial x^{\prime }}+\frac{\partial v^{\prime }}{\partial y^{\prime }}=0$

$u^{\prime }\frac{\partial u^{\prime }}{\partial x^{\prime }}+v^{\prime }\frac{\partial u^{\prime }}{\partial y^{\prime }}=\frac{{\partial }^2{u}^{\prime }}{\partial {y^{\prime }}^2}$

Obtention de l'équation de Blasius

Nouvelle variable

Clairement on recherche le champ de vitesse tangentielle de la forme $u^{\prime }=f(x^{\prime },y^{\prime })$. Or il demeure une imprécision : la longueur caractéristique $L$. En effet, celle-ci étant arbitraire, il est indispensable que la solution n'en dépende pas. Dans cette démarche, il s’avèrerait que la solution ne dépende que d'une seule variable, combinant $x^{\prime }$ et $y^{\prime }$ afin qu'elle soit indépendante de cette longueur $L$. Posons simplement :

$\theta =\frac{y^{\prime }}{\sqrt{x^{\prime }}}=\frac{y}{\sqrt{x}}\sqrt{\frac{U}{\nu }}$

Cette variable est indépendante de la longueur $L$.

On en déduit que la solution cherchée est de la forme :

$\frac{u}{U}=f(\theta )\text{soit}u=Uf(\theta )$

Équation de Blasius

Avec cette nouvelle variable, il suffit d'exprimer chaque terme de l'équation de mouvement pour aboutir à l'équation de Blasius :

$f^{″}(\theta )+\frac{1}{2}{f}^{\prime }(\theta )\underset{0}{\overset{\theta }{\int }}f(t)\mathrm{d}t=0$

On vérifie que $\theta$ est bien l'unique variable du problème.

Résolution approchée

Cas des faibles valeurs de $\theta$

Rappelons qu'à $x$ fixé, $\theta$ varie de la même manière que $y$, ainsi la variable est représentative de la distance vis-à-vis de la paroi au niveau de la couche limite. Appliquons un développement limité de la fonction pour $\theta \simeq 0$ :

$f(\theta )=f(0)+\theta f^{\prime }(0)+\frac{{\theta }^2}{2}{f}^{″}(0)+\frac{{\theta }^3}{6}{f}^{(3)}(0)+\frac{{\theta }^4}{24}{f}^{(4)}(0)+o({\theta }^4)$

Nous savons déjà que $f(0)=0$, cela permet d'affirmer, à l'aide de l'équation de Blasius, que $f^{″}(0)=0$. On peut montrer qu'il en est de même pour $f^{(3)}(0)$ et que de plus $f^{(4)}(0)=-\frac{1}{2}{f^{\prime }}^2(0)$

Modèle:Démonstration

Profil des vitesses proche de la paroi

Pour des faibles valeurs de $\theta$ nous avons donc l'approximation suivante :

$f(\theta )=(\theta -\frac{{\theta }^4}{48}{f}^{\prime }(0))f^{\prime }(0)$

On en déduit donc que proche de la paroi, le profil de la vitesse tangentielle varie linéairement puis suit un profil plus concave (terme en ${\theta }^4$)

Cas des grandes valeurs de $\theta$

Des grandes valeurs de $\theta$ signifient que l'on s'éloigne de plus en plus de la paroi ainsi le champ de vitesse dans la couche limite $u\to U$. Autrement dit, comme $u=Uf(\theta )$ on a donc $f(\theta )\to 1$. L'équation limite de Blasius s'écrit alors :

$f^{″}(\theta )\approx -\frac{1}{2}\theta f^{\prime }(\theta )$

On intègre aisément :

$f^{\prime }(\theta )=f^{\prime }(0){\mathrm{e}}^{-\frac{{\theta }^2}{4}}$

Cette expression montre que $f^{\prime }(\theta )$ converge exponentiellement vers $0$ impliquant que $f$ atteint de la même manière sa valeur asymptotique $1$. De plus, l'exponentielle atteint 98 % de sa valeur finale lorsque son paramètre est de l'ordre de 6, on en déduit que l'on sort de la couche limite pour $\theta \approx 5$. Ainsi le domaine de la fonction $f$ est donnée : $f(\theta )\in [0,1]$ pour $\theta \in [0,5]$

Ces résultats sont en accord avec le principe de couche limite : dès qu'on s'éloigne de la paroi on retrouve l'écoulement uniforme et en combinant les deux comportements aux bords on s'imagine que le profil des vitesses va passer de manière abrupte de son comportement linéaire à son comportement asymptotique.

Voir aussi

Articles connexes

• Couche limite
• Nombre de Reynolds
• Écoulement laminaire

Bibliographie

• Hydrodynamique physique, Étienne Guyon, Jean Pierre Hulin, Luc Petit, (EDP Sciences, 2001) Modèle:ISBN