Calculateur de distance y+

Lors de la création d'un maillage CFD, le dimensionnement correct des cellules proches de la paroi est essentiel pour une résolution précise de la couche limite. La distance adimensionnelle à la paroi y+ détermine si la première cellule se situe dans la sous-couche visqueuse (y+ ≈ 1, requis pour les modèles de turbulence à bas nombre de Reynolds et la LES résolue à la paroi) ou dans la région de la loi logarithmique (30 < y+ < 300, utilisée avec les lois de paroi). Cette page fournit les formules de l'estimation de la distance y+ et un calculateur qui renvoie la hauteur de première cellule requise pour toute combinaison de conditions d'écoulement et de y+ souhaité.

Formules pour estimer la hauteur de première cellule à partir de y+

Le nombre de Reynolds pour une vitesse en écoulement libre \(U_\infty\) et une longueur de référence \(L\) est :

$$Re = \frac{\rho \cdot U_{\infty} \cdot L}{\mu}$$

La corrélation de Schlichting pour le coefficient de frottement est utilisée pour estimer le coefficient de frottement \(C_f\). Elle est valide pour \(Re_x < 10^9\) :

$$C_f = \left[ 2 \cdot \log_{10}(Re_x) - 0.65 \right]^{-2.3} \quad \text{pour} \quad Re_x < 10^9$$

La contrainte de cisaillement pariétale découle du coefficient de frottement :

$$\tau_w = C_f \cdot \tfrac{1}{2} \, \rho \, U_{\infty}^2$$

La vitesse de frottement \(u^*\) est dérivée de la contrainte de cisaillement pariétale et de la masse volumique du fluide :

$$u^* = \sqrt{\frac{\tau_w}{\rho}}$$

Enfin, la distance à la paroi (hauteur de première cellule) pour un y+ souhaité est :

$$y = \frac{y^+ \cdot \mu}{\rho \cdot u^*}$$

Le calculateur ci-dessous évalue toutes ces équations et renvoie les résultats intermédiaires (Re, Cf, τw, u*) ainsi que la distance à la paroi en millimètres.

Calculateur de distance y+

Les valeurs par défaut de masse volumique et de viscosité dynamique sont pour l'air à 1 atm et 20 °C. Pour l'eau dans les mêmes conditions, utilisez ρ = 998,21 kg/m3 et μ = 1,002×10−3 Pa·s.

Entrée
Sortie

Avertissement : Les outils, calculateurs et formules fournis sur ce site sont destinés uniquement à des fins éducatives et informatives. Bien que nous nous efforçions de garantir l'exactitude, nous ne pouvons pas garantir que les résultats seront applicables à vos circonstances spécifiques. Les utilisateurs sont encouragés à vérifier les résultats de manière indépendante et à consulter un professionnel qualifié si nécessaire. En utilisant ces outils, vous reconnaissez que l'utilisation de toute information obtenue à partir de ce site se fait à vos propres risques.

Pour les géométries complexes où une estimation de plaque plane n'est pas suffisante, notre équipe CFD peut réaliser des simulations préliminaires pour vérifier la résolution du maillage proche paroi et ajuster les couches prismatiques de la couche limite en conséquence. Pour en savoir plus sur les fondamentaux, consultez notre cours Introduction au Computational Fluid Dynamics.

Questions fréquentes

Questions courantes sur y⁺ et le maillage proche paroi.

Cela dépend du modèle de turbulence. Les modèles à bas nombre de Reynolds (par ex. SST k-ω sans lois de paroi) et la LES résolue à la paroi nécessitent y+ ≈ 1, ce qui signifie que la première cellule doit se situer dans la sous-couche visqueuse. Les approches avec lois de paroi standard (par ex. k-ε standard) fonctionnent mieux avec 30 < y+ < 300, plaçant la première cellule dans la région de la loi logarithmique. Les valeurs dans la zone tampon (5 < y+ < 30) doivent être évitées, car ni le modèle de sous-couche visqueuse ni le modèle de loi de paroi n'est valide dans cette zone.

La longueur de référence est la longueur caractéristique de la surface sur laquelle la couche limite se développe. Pour l'aérodynamique externe, il s'agit typiquement de la corde (profil aérodynamique) ou de la longueur du corps (véhicule). Pour les écoulements internes, c'est généralement le diamètre hydraulique du conduit. Le choix de la longueur de référence affecte le nombre de Reynolds et donc le coefficient de frottement estimé — utilisez la longueur qui représente le mieux le développement de la couche limite dans votre application.

La formule de Schlichting fournit une estimation empirique rapide du coefficient de frottement pour une couche limite turbulente sur plaque plane, valide jusqu'à Re ≈ 109. Elle est suffisamment précise pour le dimensionnement du maillage, où l'objectif est une estimation de l'ordre de grandeur de la hauteur de première cellule. Pour les géométries avec de forts gradients de pression, du décollement ou de la courbure, la contrainte de cisaillement pariétale réelle différera de l'estimation sur plaque plane, donc la valeur y calculée doit être considérée comme un point de départ et vérifiée après l'exécution de la simulation.