Modélisation de la turbulence

Un article de Wikipédia, l'encyclopédie libre.
représentation théorique et très simplifiée des turbulences générées en aval d'un arbre par un vent régulier.

La modélisation de la turbulence est une branche de la mécanique des fluides utilisée pour prédire le comportement d'un écoulement dans lequel tout ou partie du fluide est turbulent.

Introduction[modifier | modifier le code]

La présence d'une vorticité dans un écoulement ne fait pas nécessairement de celui-ci un écoulement turbulent. Le terme est réservé aux situations où de nombreuses échelles tourbillonnaires sont présentes et interagissent dans la cascade turbulente. Celle-ci est bornée aux petites échelles par la dimension de Kolmogorov, en-dessous de laquelle les tourbillons sont dissipés par viscosité.

Un tel écoulement est décrit par les équations de Navier-Stokes mais la faible taille de la dimension de Kolmogorov interdit en pratique une simulation numérique directe (en anglais DNS pour Direct Numerical Simulation), sauf pour des expériences numériques destinées à comprendre les mécanismes mis en jeu.

Outre la simulation directe les méthodes mises en œuvre pour résoudre ce problème reposent sur la physique statistique : la turbulence est considérée comme un processus statistique dont on suppose qu'il peut être décrit par la seule distribution temporelle en chaque point. L'approche repose sur un certain nombre d'étapes :

  • l'écriture d'équations décrivant valeurs moyennes et fluctuations,
  • modélisation des termes liés aux fluctuations,
  • s'il y a lieu, raccorder ces termes aux descriptions standard des lois décrivant l'écoulement au voisinage de la paroi.

Il est également possible d'utiliser des méthodes hybrides dites de simulation des grandes échelles (en anglais LES pour Large Eddy Simulation) dans lesquelles on filtre le spectre de turbulence : les grandes échelles sont capturées par le calcul, les petites modélisées comme ci-dessus.

Les équations de Navier-Stokes moyennées[modifier | modifier le code]

On s'intéresse à un fluide incompressible décrit par les équations de Navier-Stokes correspondantes

On note p la pression, ρ la masse volumique, μ la viscosité dynamique et

 le tenseur des déformations
  le tenseur des contraintes visqueuses

En tenant compte de l'équation d'incompressibilité on remarque que

On définit l'opérateur Υ (ui) pour l'équation de conservation de la quantité de mouvement (le changement d'indice servira par la suite)

Le milieu est décrit par une distribution statistique des vitesses et l'on suppose que l'on peut caractériser ce milieu par la moyenne temporelle et la fluctuation de la vitesse en un point r

L'énergie cinétique de la turbulence est

En introduisant cette expression de la vitesse dans les équations de Navier-Stokes on obtient les équations moyennées introduites par Osborne Reynolds en 1895[1] :

On a défini le tenseur des contraintes de Reynolds :

Comme tout tenseur de contraintes ce tenseur est symétrique[2]. Le problème de la turbulence consiste à exprimer les 6 quantités indépendantes qu'il contient.

Équation de transport des contraintes[modifier | modifier le code]

Julius C. Rotta[3] a introduit en 1951 une équation de transport sur les contraintes de Reynolds. Pour y parvenir on utilise l'opérateur défini plus haut en écrivant[4],[5]

soit

avec

Expression Signification physique
 Production : transfert de l'énergie de l'écoulement moyen vers la turbulence
 Transport de la turbulence (contient une corrélation triple)
 Redistribution de l'énergie turbulente (retour vers l'état isotrope)
 Diffusion de la contrainte
 Dissipation visqueuse

où δij est le symbole de Kronecker.

Ces 6 équations contiennent 22 nouvelles inconnues. Il faut donc simplifier (modéliser) en remplaçant ces termes par des expressions des variables déjà présentes comme les composantes de τij. L'approche classique a été introduite par Kemal Handjalić et Brian Launder (1972)[6],[7].

Modèles à N équations de transport[modifier | modifier le code]

Ces modèles sont appelés en anglais Reynolds Averaged Navier-Stokes ou en abrégé RANS.

L'hypothèse de Boussinesq[modifier | modifier le code]

En 1877 Joseph Boussinesq a proposé d'écrire ce tenseur comme le tenseur des contraintes dans le cas d'un fluide newtonien en faisant intervenir une viscosité de turbulence μt[8]

Le problème est réduit à la connaissance de k et μt, cette dernière valeur n'étant pas une propriété du fluide.

Les modèles à deux équations[modifier | modifier le code]

En prenant la trace de l'équation des contraintes de Reynolds ci-dessus on obtient une équation de transport pour k

où ε est la dissipation

Celle-ci peut être obtenue en écrivant l'équation

soit

En fait cette expression comporte au second membre des termes très difficiles à modéliser et on se contente d'écrire un second membre analogue à l'équation sur l'énergie cinétique turbulente[4].

La viscosité turbulente est déduite de l'analyse dimensionnelle

où Cμ est une constante de modélisation.

Les modèles le plus connu utilisé dans ce domaine est le modèle k - ε de William P. Jones et Brian Launder[9] publié en 1972 et reformulé ultérieurement[10].

Il est également possible de travailler sur le taux de dissipation

Ce type de modèle dit modèle k - ω a été introduit[4] par Andreï Kolmogorov en 1942 à une époque où il n'était pas possible de le résoudre[11]. Sa forme actuelle est due à David C. Wilcox[12].

Modèle à une équation de transport[modifier | modifier le code]

Ce type de modèle a été introduit dans les années 60[4]. On part de la viscosité turbulente ci-dessus avec Cμ = 1 et l'on dérive

Le plus connu de ces modèles est sans doute le modèle Spalart-Allmaras (1992) de Philippe R. Spalart et Steven R. Allmaras pour les problèmes de couche limite en écoulement compressible[13].

Modèle à longueur de mélange[modifier | modifier le code]

Le modèle utilisant une longueur de mélange, également nommé à zéro équation de transport, a été introduit par Ludwig Prandtl en 1925[4],[14]. Par analogie avec la théorie cinétique des gaz il a supposé que l'on pouvait construire une viscosité cinématique à partir du produit d'une vitesse caractéristique u par une longueur de mélange lm et que le temps caractéristique formé à partir de ces deux quantités était du même ordre de grandeur que celui associé au cisaillement moyen

d'où la composante correspondante du tenseur de Reynolds

Cette expression peut être généralisée par :

L'expression de lm est spécifique d'un problème donné[6],[5],[4].

Modèles de simulation des grandes échelles[modifier | modifier le code]

La méthode SGS[15] ou en anglais LES consiste à séparer les échelles de turbulence en

  • grandes échelles calculées directement,
  • petites échelles, modélisées.

La première étape du processus consiste à définir un filtre passe-bas par l'intermédiaire du produit de convolution

Le filtre est normalisé :

Ce n'est pas un projecteur : . De plus cet opérateur ne commute pas avec la dérivée.

L'exemple le plus simple est le filtre « chapeau » (en anglais top hat) basé sur la taille de maille Δ

On écrit la solution sous la forme de la somme de la valeur filtrée et d'une perturbation de petite échelle, laquelle n'a la signification d'une fluctuation temporelle.

on peut alors écrire les équations de Navier-Stokes filtrées :

où tij est le tenseur introduit par Anthony Leonard[16] :

On remarquera que si G était l'opérateur moyenne de Reynolds les quatre premiers termes s'annuleraient. Par ailleurs si tij respecte l'invariance galiléenne, ce n'est pas vrai pour chacun des termes qui le composent.

Pour fermer le problème il faut définir une approximation dans la maille, par exemple du type longueur de mélange (voir ci-dessus) comme l'a fait Joseph Smagorinsky (1963)[17]

où Cs ~ 0.1 est une constante de modélisation liée à la constante de Kolmogorov.

Références[modifier | modifier le code]

  1. (en) John W. S. Rayleigh, « On the Dynamical Theory of Incompressible Viscous Fluids, and the Determination of the Criterion », Philosophical Transactions of the Royal Society A, vol. clxxxiv,‎ (lire en ligne)
  2. (en) Rutherford Aris, Vectors, Tensors, and the Basic Equations of Fluid Mechanics., Dover Publications, , 286 p. (ISBN 0-486-66110-5, lire en ligne)
  3. (de) J. C. Rotta, « Statistische Theorie nichthomogener Turbulenz », Zeitschrift fur Physik, vol. 129,‎ , p. 547-572
  4. a b c d e et f (en) David C. Wilcox, Turbulence Modeling for CFD : CD-ROM, DCW Industries, , 522 p. (ISBN 1-928729-08-8, lire en ligne)
  5. a et b (en) Christophe Bailly et Geneviève Comte-Bellot, Turbulence, Springer, (ISBN 978-3-319-16159-4)
  6. a et b (en) S. B. Pope, Turbulent Flows, Cambridge University Press,
  7. (en) K. Hanjalić et B. E. Launder, « A Reynolds stress model of turbulence and its application to thin shear flows », Journal of Fluid Mechanics, vol. 52, no 4,‎ , p. 609-638
  8. J. Boussinesq, « Essai sur la théorie des eaux courantes », Comptes rendus de l'Académie des Sciences, vol. 23,‎ , p. 1-680 (lire en ligne)
  9. (en) W. P. Jones et B. E. Launder, « The prediction of laminarization with a two-equation model of turbulence », International Journal of Heat and Mass Transfer, vol. 15, no 2,‎ , p. 301-314
  10. (en) B. E. Launder et D. B. Spalding, « The Numerical Computation of Turbulent Flows », Computer Methods in Applied Mechanics and Engineering, vol. 3, no 2,‎ , p. 269-289
  11. (ru) A. Kolmogorov, « Equation of Turbulent Motion of an Incompressible Fluid », Doklady Akademii Nauk,‎
  12. (en) D. C. Wilcox, « Reassessment of the Scale-Determining Equation for Advanced Turbulence Model », AIAA Journal, vol. 26, no 11,‎ , p. 1299-1310
  13. (en) P. R. Spalart et S. R. Allmaras, « A One-Equation Turbulence Model for Aerodynamic Flows », AIAA Paper, nos 92-0439,‎ (lire en ligne)
  14. (de) L. Prandtl, « Bericht uber Untersuchungen zur ausgebildeten Turbulenz », Zeitschrift für Angewandte Mathematik und Mechanik,‎ , p. 136-139
  15. (en) P. Sagaut, Large Eddy Simulation for Incompressible Flows : An Introduction, Springer-Verlag, , 556 p. (ISBN 978-3-540-26344-9, lire en ligne)
  16. (en) A. Leonard, « Energy Cascade in Large-Eddy Simulations of Turbulent Fluid Flows », Advances in Geophysics, vol. A 18,‎ , p. 237–248
  17. (en) J. S. Smagorinsky, « General Circulation Experiments with the Primitive Equations: I. The basic Experiment », Monthly Weather Revue, vol. 91, no 3,‎ , p. 99-164 (lire en ligne)

Liens externes[modifier | modifier le code]