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 :
Équation de mouvement dans la couche limite (Théorie de Prandtl)
Analysons l'écoulement bidimensionnel stationnaire dans le plan près d'une plaque plane placée en , pour un écoulement extérieur potentiel qu'on suppose parallèle à la paroi.
La dimension caractéristique dans la direction parallèle à l'écoulement est une longueur arbitraire que l'on suppose très grande devant la dimension caractéristique perpendiculaire à l'écoulement (épaisseur de la couche limite) :
Le raisonnement suivant est basé sur l'existence de ces 2 échelles.
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 :
Démonstration
La définition même de la couche limite réside dans le fait qu'elle représente la région de l'écoulement où les effets visqueux sont aussi importants que les effets inertiels. Ainsi nous allons comparer séparément les termes diffusifs et les termes convectifs. Commençons par utiliser l'équation de continuité :
Poursuivons avec les termes de diffusion, on remarque qu'on peut donc aisément négliger le phénomène dans la direction longitudinale :
Par contre concernant la convection, les termes sont comparables :
Enfin, comme convection et diffusion sont du même ordre de grandeur, on en déduit :
Soit en faisant intervenir le nombre de Reynolds global on a la relation : . On peut ajouter qu'avec l'équation de continuité on a la relation entre les ordres de grandeurs des vitesses tangentielles et normales :
Avec l'étude des ordres de grandeurs on ne peut considérer uniquement l'équation pour la vitesse tangentielle (sans diffusion tangentielle). En effet, les termes en sont tous très faibles réduisant la deuxième projection à . Cela signifie que la pression est constante dans l'épaisseur de la couche limite (à abscisse fixée) et elle peut donc être déterminée via l'équation de Bernoulli en dehors de la couche limite. Il faut noter que dans le cas de la couche limite de Blasius, l'écoulement extérieur est supposé uniforme de vitesse . Ainsi on montre que le gradient de pression est nul :
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 à il est nécessaire de comparer à . Il en est de même pour les variables d'espace. Les nouvelles variables sont donc :
Avec ces nouvelles variables, le système devient simplement :
Clairement on recherche le champ de vitesse tangentielle de la forme . Or il demeure une imprécision : la longueur caractéristique . 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 et afin qu'elle soit indépendante de cette longueur . Posons simplement :
Cette variable est indépendante de la longueur .
Démonstration
On en déduit que la solution cherchée est de la forme :
Avec cette nouvelle variable, il suffit d'exprimer chaque terme de l'équation de mouvement pour aboutir à l'équation de Blasius :
Démonstration
Dans un premier temps, il faut remarquer que la variable . Ainsi nous aurons besoin par la suite des dérivées partielles :
Il s'agit désormais d'exprimer chacun des termes de l'équation de mouvement. Commençons par le terme de diffusion :
Le premier terme convectif est simple à exprimer :
Contrairement au second qui nécessite le passage par l'équation de continuité :
Pour trouver la vitesse on intègre cette équation selon à fixé:
Puis à l'aide d'une intégration par parties :
On rappelle que lorsque c'est lorsque soit au niveau de la paroi. Or la vitesse normale est nulle à ce niveau (le fluide ne pénètre pas dans la paroi) ce qui donne finalement :
Le second terme de convection s'écrit alors :
On regroupe tous les termes ainsi exprimés :
Après simplifications, nous obtenons l'équation de Blasius :
On vérifie que est bien l'unique variable du problème.
Rappelons qu'à fixé, varie de la même manière que , 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 :
Nous savons déjà que , cela permet d'affirmer, à l'aide de l'équation de Blasius, que . On peut montrer qu'il en est de même pour et que de plus
Démonstration
Dérivons une première fois l'équation de Blasius :
On obtient donc directement que . Et si l'on dérive une nouvelle fois :
De la même manière on obtient :
Profil des vitesses proche de la paroi
Pour des faibles valeurs de nous avons donc l'approximation suivante :
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 )
Des grandes valeurs de signifient que l'on s'éloigne de plus en plus de la paroi ainsi le champ de vitesse dans la couche limite . Autrement dit, comme on a donc . L'équation limite de Blasius s'écrit alors :
On intègre aisément :
Cette expression montre que converge exponentiellement vers impliquant que atteint de la même manière sa valeur asymptotique . 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 . Ainsi le domaine de la fonction est donnée : pour
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.