La simulation des grandes structures de la turbulence (SGS ou en anglais LES pour Large Eddy Simulation) est une méthode utilisée en modélisation de la turbulence. Elle consiste à filtrer les petites échelles qui sont modélisées et en calculant directement les grandes échelles de la cascade turbulente.
Cette méthode a été introduite par Joseph Smagorinsky en 1963 et utilisée pour la première fois par James W. Deardoff en 1970. Elle permet de calculer un écoulement turbulent en capturant les grandes échelles pour un coût raisonnable. Le milieu aux petites échelles étant raisonnablement supposé isotrope peut être décrit par une méthode simple. Les problèmes de couche limite, de flamme de diffusion ou problèmes diphasiques requièrent des modèles adaptés.
La première étape du processus consiste à définir un filtre passe-bas G pour la vitesse u par l'intermédiaire du produit de convolution. On traite ici d'un filtre spatial mais il est possible d'utiliser des filtres spatio-temporels.
Le filtre est normalisé
Ce n'est pas un projecteur : .
Cet opérateur commute avec la dérivée en temps mais ne commute avec la dérivée en espace que G est le même en tout point.
On traite ci-dessous les équations de Navier-Stokes dans un fluide incompressible. On écrit la vitesse 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 comme dans la moyenne de Reynolds
Dans les équations de Navier-Stokes dans un fluide incompressible on peut alors écrire les équations de Navier-Stokes filtrées
où ρ est la masse volumique, p la pression et
le tenseur des déformations | |
le tenseur des contraintes visqueuses |
μ est la viscosité dynamique et tij est le tenseur introduit par Anthony Leonard
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.
L'énergie cinétique filtrée est
où kf est l'énergie cinétique correspondant à la vitesse filtrée. On obtient une équation de transport en multipliant l'équation de quantité de mouvement filtrée par
avec
taux de dissipation visqueuse de la vitesse filtrée (petit si la coupure est faite dans la région inertielle du spectre turbulent) | |
« taux de dissipation » de la vitesse dans la maille (en fait un terme inertiel pouvant être négatif) |
Π représente l'énergie transférée des grandes vers les petites échelles.
Diverses méthodes sont utilisables. L'exemple le plus simple est le filtre « chapeau » (en anglais top hat) basé sur une longueur Δ de l'ordre de la taille de maille
Cela correspond à un filtre coupant les longueurs d'onde λ telles que
Un autre fltre couramment utilisé est le filtre gaussien
Pour fermer le problème il faut définir une approximation dans la maille, par exemple du type longueur de mélange
où Cs ~ 0.1 est une constante de modélisation parfois nommée constante de Smagorinsky. Elle est liée à la constante de Kolmogorov.
D'où la viscosité turbulente