Lire la première partie

Bourguet

Description

Niveau: Secondaire, Lycée, Première

• mémoire

Lire la première partie de la thèse

• simulation numérique


• turbulents instationnaires


• étape de modélisation de la turbulence


• équations de navier-stokes en moyenne


• contexte industriel


• modélisation d'écoulements pariétaux


• problème pratique de la fermeture des équations de navier-stokes en moyenne


• modèle physique

Sujets

• Première
• Mémoire
• Simulation
• Viscosity
• Simulation informatique

Informations

• Publié par: profil-mazor-2012
• Nombre de lectures: 27
• Langue: Français
• Poids de l'ouvrage: 4 Mo

Extrait

Lire lapremièrepartie de la thèse

Chapitre 4

Mod´elisationd’´ecoulements parie´tauxturbulents instationnaires

Aperçu 1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45 1.1 Un premier aperçu des grandes classes de modélisation . . . . . . . . . . . . . . 46 1.2 Contribution de la présente étude . . . . . . . . . . . . . . . . . . . . . . . . . . 48 2 Différentes approches dans la littérature . . . . . . . . . . . . . . . . . . . . 48 2.1 Simulation numérique directe des équations de Navier-Stokes . . . . . . . . . . 49 2.2 Simulations aux grandes échelles . . . . . . . . . . . . . . . . . . . . . . . . . . 49 2.3 Approches statistiques fondées sur la moyenne de Reynolds . . . . . . . . . . . 51 2.4 Approches statistiques avancées - Organised Eddy Simulation . . . . . . . . . . 52 2.5 Approches hybrides - Detached Eddy Simulation . . . . . . . . . . . . . . . . . 55 3 Modèles de fermeture dans le cadre d’une modélisation statistique de la turbulence . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56 3.1 Modélisation au second ordre . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56 3.2 Modèles au premier ordre linéaires fondés sur une viscosité de turbulence . . . 59 3.3 Modèles non-linéaires et algébriques explicites . . . . . . . . . . . . . . . . . . . 65 4 Modèle OES à viscosité de turbulence tensorielle . . . . . . . . . . . . . . . 71 4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72 4.2 Physical analysis of turbulent stress-strain relation in the near-region . . . . . . 73 4.3 Anisotropic OES modelling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75 4.4 Numerical simulation of strongly detached flows around bodies . . . . . . . . . 77 4.5 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81 4.6 Appendix A. Transport equations for the stress-strain projection coefficients . . 81 4.7 Appendix B. Summary of the OES anisotropic first-order model . . . . . . . . 82 4.8 Compléments à l’article . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83 5 Prédiction d’écoulements turbulents compressibles . . . . . . . . . . . . . . 87 5.1 Quelle moyenne statistique dans le cas d’écoulements compressibles ? . . . . . . 87 5.2 Equations de Navier-Stokes compressibles en moyenne de Favre . . . . . . . . . 88 5.3 Modèles de fermeture . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88 6 Conclusion - modélisation “haute-fidélité” des écoulements turbulents . . 89

1 Introduction Malgré de nombreuses avancées théoriques dans l’analyse physique des mécanismes fondamentaux de la turbulence et le développement d’outils de calcul numérique de plus en plus performants, la prédiction d’écoulements pariétaux instationnaires à grands nombres de Reynolds reste un problème difficile. En effet, la résolution numérique directe des équations de Navier-Stokes tridimensionnelles et instationnaires

45

Chapitre 4. Modélisation d’écoulements pariétaux turbulents instationnaires

n’est généralement pas envisageable pour la prédiction d’écoulements réalistes et la mise en œuvre de stratégies de modélisation de la turbulence est inévitable. De multiples méthodologies, notamment des méthodes de macrosimulation, ont été développées dans ce sens. Plusieurs approches théoriques différentes peuvent être considérées selon la configuration physique d’intérêt et la nature des résultats attendus, es- sentiellement la prédiction de quantités physiques moyennes ou filtrées, stochastiques ou déterministes. Dans le cadre de la présente étude, l’étape de modélisation de la turbulence par une méthode de macro- simulation peut être considérée comme une première phase de “réduction de la dimension” du problème. En effet, dans le contexte du développement d’une modélisation hiérarchique d’écoulements réalistes, les simulations mettant en jeu une modélisation de l’effet de la turbulence non résolue constituent une approximation du modèle physique “haute-fidélité” que représentent les équations de Navier-Stokes ins- tationnaires tridimensionnelles. A l’inverse, les approches de simulation décrites dans ce chapitre peuvent également être considérées comme des modèles physiques détaillés dont la complexité numérique peut encore être réduite grâce à une méthode de modélisation d’ordre faible fondée sur la décomposition or- thogonale aux valeurs propres. Ce point sera abordé dans les chapitres suivants. Dans une optique de modélisation hiérarchique, l’étape de modélisation de la turbulence est donc cruciale pour le développe- ment de modèles d’ordre réduit pertinents, représentatifs de l’écoulement réel.

1.1 Un premier aperçu des grandes classes de modélisation Dans le cadre de la prédiction d’écoulements instationnaires autour d’obstacles à grands nombres de Reynolds, la simulation numérique directe ouDirect Numerical Simulation(DNS) nécessite des discréti- sations spatiale et temporelle extrêmement fines pour capturer toutes les échelles des quantités physiques aléatoires mises en jeu. Les approches les plus couramment considérées pour contourner cette difficulté sont brièvement présentées. La méthodologie la plus proche de la simulation directe sur le plan théo- rique est la simulation aux grandes échelles ouLarge Eddy Simulation(LES). Cette approche, initiée par Smagorinsky(1963) dans le domaine météorologique, consiste à appliquer un filtre spatial aux différentes quantités physiques et à ne simuler que les grandes échelles, les échelles non résolues étant prises en compte par l’introduction d’un modèle de fermeture dans les équations de Navier-Stokes filtrées. Cette approche est, par nature, tridimensionnelle et comme cela sera explicité dans ce chapitre, elle donne accès, comme la DNS, à des réalisations des variables aléatoires que constituent les quantités physiques mises en jeu. La démarche proposée par la LES s’applique a priori à tous les types d’écoulements, qu’ils présentent ou non une cohérence spatio-temporelle marquée par la présence de structures organisées. Une des difficultés liées à cette approche est le traitement des résultats de simulation. En effet, une réalisation de l’écoulement ne peut être considérée comme représentative de la moyenne d’ensemble que sous des hypothèses d’ergodicité et d’homogénéité délicates à démontrer dans la pratique. Dans le cas général, plu- sieurs résolutions successives doivent être envisagées avant d’accéder aux grandeurs statistiques. De plus, la modélisation des échelles non résolues ou échelles de sous-maille passe en général par l’introduction d’un concept de viscosité de turbulence qui est essentiellement adapté à la modélisation d’une partie du spectre d’énergie cinétique turbulente dans le cas d’une turbulence homogène isotrope. Etant donné que la taille du filtre spatial est dans la pratique liée à la finesse de la discrétisation spatiale, ce point implique que l’approche LES doive tendre vers une simulation numérique directe dans les régions de proche-paroi. Cela rend la LES difficilement applicable pour la modélisation d’écoulements pariétaux à grands nombres de Reynolds dans un contexte industriel (Davidsonet al.,2003) et justifie le développement de méthodes de macrosimulation hybrides, couplant la simulation aux grandes échelles à des méthodes statistiques plus adaptées à la simulation des écoulements en proche-paroi. Les méthodes statistiques constituent une alternative très largement répandue aux simulations di- recte et aux grandes échelles. D’une manière générale, elles conduisent à décomposer l’ensemble des variables physiques en termes moyen et fluctuant. L’approche la plus classique qui consiste à considérer une moyenne d’ensemble 1 conduit aux équations de Navier-Stokes en moyenne de Reynolds ouRey- nolds Averaged Navier-Stokes(RANS). Les variables considérées sont alors déterministes et non plus stochastiques comme dans le cas de la DNS ou de la LES, ce qui simplifie l’exploitation des résultats de simulation. Les méthodes dérivées de cette approche ne sont en théorie valables que pour des écoulements en équilibre statistique autrement dit des écoulements dont les quantités physiques caractéristiques sont 1 La moyenne d’ensemble peut être remplacée par une moyenne temporelle dans le cas de processus ergodiques.

46

1. Introduction

des processus stochastiques stationnaires 2 en équilibre dans le sens de la théorie deet Kolmogorov(1941) tels que la zone lointaine d’un sillage, d’un jet ou d’une zone de mélange par exemple. Le fait de considérer les équations de Navier-Stokes en moyenne fait apparaître un terme supplémentaire traduisant la “perte d’information” par rapport aux équations originelles définies dans un cadre stochastique : les corrélations doubles des fluctuations de vitesse en un point, appelées tensions de Reynolds lorsque cette moyenne est mise en œuvre et plus généralement contraintes turbulentes. Ainsi les équations de Navier-Stokes en moyenne constituent un système ouvert et leur résolution nécessite une représentation (simulation ou modélisation) des contraintes turbulentes. Deux types d’approches sont généralement distinguées : d’une part l’approche dite “au second ordre” qui assure la prédiction des contraintes turbulentes au moyen d’équations de transport supplémentaires associées à ces quantités où les termes inconnus tels que les corrélations triples sont modélisés, d’autre part, les modèles “au premier ordre” qui relient algébrique- ment les tensions de Reynolds aux grandeurs physiques moyennes. Les modèles au premier ordre les plus répandus sont les modèles linéaires qui utilisent une hypothèse de fermeture fondée sur une analogie avec la loi constitutive des fluides newtoniens : la loi deBoussinesq(1877). Le tenseur des contraintes est alors relié linéairement au taux de déformation moyen grâce à un concept de viscosité de turbulence. Comme cela sera détaillé dans ce chapitre, une relation constitutive linéaire peut conduire à une mo- délisation erronée de certaines propriétés structurales de la turbulence entraînant d’importantes erreurs d’estimation des quantités globales d’intérêt telles que les coefficients aérodynamiques par exemple. Afin de pallier certaines limitations du modèle linéaire, des lois constitutives non-linéaires ont été développées. Elles consistent à inclure, dans la relation constitutive, des tenseurs d’ordre supérieurs issus de combi- naisons des tenseurs des taux de déformation et rotation moyens (Pope,1975, par exemple). Plusieurs approches ont été envisagées par analogie aux lois constitutives considérées en rhéologie pour les fluides visco-élastiques ou à effet de mémoire donnant lieu auxNon Linear Eddy Viscosity Models(NLEVM) (Shihet al.,1993, par exemple) ou encore en dégénérant les équations aux dérivées partielles de modèles au second ordre pour en dériver des expressions algébriques des tensions turbulentes,Explicit Algebraic Stress Modelling(EASM) (Pope,1975;Gatski & Speziale,1993;Wallin & Johansson,2000). Ce type d’approche a notamment pour objectif une meilleure capture de l’anisotropie des contraintes turbulentes normales en proche-paroi par rapport au modèle linéaire fondé sur l’hypothèse deBoussinesq(1877) qui prédit un comportement isotrope de ces tensions. Dans le cas où les processus physiques étudiés ne sont pas statistiquement stationnaires, par exemple lorsque des structures organisées apparaissent dans l’écoulement, il semble que ces phénomènes puissent être considérés non plus comme des fluctuations aléatoires mais comme des évolutions déterministes des propriétés statistiques de ces processus. Cela conduit à l’approche RANS instationnaire ouUnsteady RANS(URANS) qui est la méthodologie de modélisation des écoulements turbulents la plus largement mise en œuvre. Les équations simulant l’évolution des quantités physiques moyennes sont les mêmes que les équations RANS, au terme temporel près. Par ailleurs, les mêmes méthodes de fermeture que dans le cas stationnaire sont généralement utilisées dans ce contexte qui ne correspond pas, a priori, à leur domaine d’application. En particulier, les échelles caractéristiques de la turbulence évaluées par ces approches ne sont plus nécessairement adaptées et peuvent conduire à des prédictions erronées. Une reconsidération des approches statistiques classiques dans le cas de processus instationnaires est donc cruciale. Afin d’étendre de manière rigoureuse les approches statistiques précédentes aux cas d’écoulements ins- tationnaires caractérisés par la présence de structures organisées, une méthodologie fondée sur la décom- position triple distinguant pour chaque processus, une composante moyenne, une composante fluctuante organisée et une composante aléatoire a été imaginée (Reynolds & Hussain,1972). Cette approche s’est développée, à la fin des années 1970, alors que de nombreux travaux expérimentaux tentaient de quan- tifier la partie cohérente du mouvement turbulent par des mesures conditionnelles telles que la moyenne de phase (Cantwell,1981;Boissonet al.,1983par exemple). Dans la pratique, une distinction est faite, entre la partie cohérente du processus (composantes moyenne et fluctuante organisée) qui est considérée comme une variable déterministe et sa partie aléatoire selon la décomposition proposée parCantwell & Coles(1983le cas où l’écoulement d’intérêt présente un fort caractère périodique, l’utilisation de). Dans l’opérateur de moyenne de phase conduit à une dissociation efficace des quantités cohérentes et aléatoires 2 La stationnarité d’un processus stochastique traduit l’indépendance de ses différents moments statistiques (moyenne, corrélations doubles...) par rapport à l’origine des temps.

47

Chapitre 4. Modélisation d’écoulements pariétaux turbulents instationnaires

dans ce sens. L’effet de cette distinction sur le plan spectral est détaillé dans ce chapitre au§2.4. Substi- tuer cette moyenne conditionnelle à la moyenne d’ensemble conduit aux mêmes équations que le système URANS. Cette approche a été utilisée pour la mesure et la prédiction d’écoulements autour de profils d’aile oscillants à grand nombre de Reynolds (McCroskeyet al.,1982, notamment). Sur le plan théorique, une telle méthodologie implique le calcul des quantités organisées correspondant à la turbulence résolue et la modélisation des effets des processus aléatoires. Par rapport à l’approche RANS, la modélisation du tenseur des contraintes turbulentes apparaissant dans les équations en moyenne doit être reconsidérée, y compris dans le cas où une loi constitutive linéaire est mise en œuvre. Pour cela, l’approcheOrganised Eddy Simulation(OES) a été proposée parDervieuxet al.(1998),Braza(2000),Abalakin & Dervieux (2000),Hoarau(2002) etBrazaet al.(2006). La modification des fermetures turbulentes dans le contexte OES est décrite au§2.4. Le principal écueil rencontré quant à l’utilisation de la LES pour la simulation d’écoulements autour de corps est la finesse de discrétisation de la région de proche-paroi. Afin de limiter l’usage de la simu- lation aux grandes échelles à son domaine optimal d’application, des méthodes hybrides couplant LES et approches statistiques peuvent être considérées. En particulier, la méthodeDetached Eddy Simulation (DES) (Spalartet al.,1997;Travinet al.,2000) combine les approches LES et RANS par une sélection locale de l’échelle de longueur de la turbulence. Les approches de modélisation des écoulements turbulents présentées dans cette section introductive correspondent à un premier aperçu des méthodes détaillées dans ce chapitre. Il est important de noter que le tour d’horizon ici proposé n’est en aucun cas exhaustif.

1.2 Contribution de la présente étude Les développements effectués dans le cadre de cette thèse s’inscrivent dans le contexte de l’approche de modélisationOrganised Eddy Simulationde la présente étude concerne l’améliora-. La contribution tion de la capture des propriétés structurales de la turbulence grâce à l’introduction d’une nouvelle loi constitutive des tensions turbulentes incluant une viscosité de turbulence tensorielle. Plus précisément, les non-linéarités apparaissant dans un écoulement instationnaire fortement détaché entre le tenseur d’ani- sotropie des contraintes turbulentes et le tenseur des taux de déformations moyens sont examinées sur la base de résultats expérimentaux détaillés. Cette analyse illustre les limitations des approches linéaires fondées sur l’hypothèse deBoussinesq(1877). Afin de prendre en compte la non-colinéarité des tenseurs d’anisotropie des contraintes turbulentes et de déformation ainsi mise en évidence, une alternative aux modèles non-linéaires les plus répandus (NLEVM et EASM) est suggérée sous la forme d’une loi consti- tutive modifiée. De plus, un modèle de fermeture permettant la mise en œuvre pratique de cette loi de comportement des contraintes turbulentes est proposé. Ces développements sont détaillés au§4de ce chapitre et ont été rapportés dans deux articles publiés dans le cadre de cette thèse (Bourguetet al., 2007b,2008). Les différentes approches de modélisation précédemment évoquées sont plus précisément décrites au§2, notamment du point de vue de l’introduction des opérateurs de filtrage ou de moyenne dans les équations de Navier-Stokes. Après avoir mis en évidence le problème pratique de la fermeture des équations de Navier-Stokes en moyenne, le choix de modèles de turbulence adaptés dans le contexte des approches statistiques est détaillé au§3particulier, les modifications induites par la reconsidération des échelles. En caractéristiques de la turbulence dans le cadre de l’approche OES sont présentées. La contribution de la présente étude à l’amélioration des capacités prédictives des modèles OES est décrite au§4. Pour plus de clarté, l’ensemble des approches de modélisation ainsi que les développements menés dans cette thèse sur ce thème sont en premier lieu présentés dans le contexte des écoulements incompressibles. La transposition au cas d’écoulements compressibles est proposée au§5.

2 Différentes approches dans la littérature Dans cette section les principales approches de modélisation d’écoulements turbulents à grands nombres de Reynolds sont présentées. Pour chacune d’elles, les équations de Navier-Stokes filtrées ou en moyenne

48

2. Différentes approches dans la littérature

sont détaillées et, le cas échéant, la nécessité d’introduire des modèles de fermeture est mise en évidence dans l’optique d’une résolution numérique de ce système.

2.1 Simulation numérique directe des équations de Navier-Stokes La simulation directe consiste à résoudre numériquement les équations de Navier-Stokes tridimen- sionnelles instationnaires qui correspondent, dans le cas d’écoulements incompressibles, à l’équation de continuité et aux trois équations de conservation de la quantité de mouvement, pouri= 1,2,3: ( u u α,α t,i (+ = u 0 i u α ) ,α =− ρ 1 p ,i +νu i,αα , (4.1) oùu i désigne la i ème composante de la vitesse,ρla masse volumique du fluide supposée constante dans le cas incompressible,pla pression etνla viscosité cinématique également supposée constante dans ce contexte. ,t et ,i représentent respectivement les dérivées temporelle et spatiale dans la directioni. Pour plus de clarté, les indices grecs sont utilisés pour indiquer les sommations implicites. La résolution directe de ces équations, munies de conditions initiales et aux limites, conduit à la simulation de toutes les échelles de structure dans l’écoulement et ne nécessite donc pas de modélisation supplémentaire. En ce sens, ce type d’approche peut être considéré comme une véritable “expérience numérique” (Chassaing,2000) utile pour la compréhension des propriétés de la turbulence, notamment dans l’optique du développement d’approches de modélisation. Le principal obstacle à la mise en œuvre pratique d’une telle méthode pour des écoulements réalistes à grands nombres de Reynolds est la finesse de discrétisation nécessaire à la capture de l’ensemble du spectre comme cela a déjà été évoqué. Par ailleurs, l’accès aux propriétés statistiques d’écoulements simulés par DNS ne peut être envisagé qu’a posteriori, chaque simulation correspondant à une réalisation d’un processus aléatoire. Dans la pratique, la simulation directe reste encore un outil d’investigation scientifique dans le contexte de la prédiction des écoulements pariétaux à grands nombres de Reynolds ; des méthodes de modélisation de la turbulence doivent nécessairement être mises en œuvre pour la simulation d’écoulements turbulents instationnaires autour de géométries complexes.

2.2 Simulations aux grandes échelles En considérant que la principale limitation de la simulation directe est la finesse de discrétisation nécessaire à la capture de l’ensemble des structures de l’écoulement, la simulation aux grandes échelles constitue une alternative où seuls les processus de grande taille sont simulés alors que l’effet des petites structures est modélisé par une loi de “sous-maille”. L’idée directrice de cette méthode est que la partie du spectre négligée correspond aux petites structures obéissant aux hypothèses d’équilibre de la turbulence homogène isotrope dont l’effet peut, par conséquent, être modélisé simplement par l’introduction d’un terme assurant la dissipation de l’énergie provenant des structures résolues de plus grande taille. Un schéma comparatif de la simulation aux grandes échelles et de l’approche OES sera présenté par la suite dans le plan spectral (figure4.1issues du filtrage spatial et un exemple de fermeture). Les équations classiquement appliquée en LES sont décrits dans cette section. La LES n’a pas été utilisée dans la présente étude ; la présentation proposée est par conséquent succincte, l’objectif étant de situer l’approche mise en œuvre dans le cadre de cette thèse par rapport aux autres méthodes existantes. Pour plus de détails concernant la simulation aux grandes échelles, le lecteur pourra se référer aux articles et ouvrages de référence deRogallo & Moin(1984),Lesieur & Métais(1996) etSagaut(2002).

Equations filtrées spatialement La mise en œuvre de la simulation aux grandes échelles implique la définition d’un filtre spatial permettant de distinguer les structures résolues de celles dont les effets devront être modélisés. Pour une quantité stochastique de dépendance spatio-temporellev, la moyenne spatiale filtréevest définie comme suit : v(x, t) = Z Ω G(x,x 0 ,Δ)v(x 0 , t)dx 0 ,(4.2) oùG(x,x 0 ,Δ)désigne l’opérateur de filtrage spatial au pointx,Δreprésente la plus petite échelle de structure résolue ou plus généralement un paramètre caractéristique de la coupure spectrale de ce filtre.

49

Chapitre 4. Modélisation d’écoulements pariétaux turbulents instationnaires

L’opérateur de filtrage est supposé normé : Z G(x,x 0 ,Δ)dx 0 = 1.(4.3) Ω L’ensemble des variables de l’écoulement peuvent alors être décomposées comme suit : v(x, t) =v(x, t) +v 0 (x, t),(4.4) oùv 0 désigne la fluctuation devpar rapport à la moyenne spatiale filtrée.vétant un processus stochas- tique, les termesvetv 0 sont également des variables aléatoires. Comme le souligneChassaing(2000), l’accès aux propriétés statistiques devpose donc un problème méthodologique conséquent. Par ailleurs, afin de conduire à une forme simple des équations de Navier-Stokes filtrées, l’opérateurGdoit commuter avec les opérations de dérivations spatiale et temporelle. Ainsi, un filtre gaussien normé dont l’écart-type est proportionnel à la taille de la maille locale de discrétisation peut par exemple être considéré. Il est important de noter les propriétés suivantes, liées à la définition de la moyenne spatiale filtrée : v6=vet par conséquentv 0 =v−v6= 0 (4.5)dans le cas général. L’application du filtre spatial aux équations de Navier-Stokes (4.1) conduit au système filtré suivant : ( u u αα, it, = (+u 0 i u α ) ,α =− ρ 1 p ,i +νu i,αα . (4.6) Afin d’exprimer la moyenne filtrée du produit des vitesses en fonction du produit des vitesses filtrées, la décomposition suivante est généralement utilisée : u i u j =u i u j +L ij +R ij ,(4.7) avec L ij =u i u j −u i u j etR ij =u i u 0 j +u 0 i u j +u i 0 u 0 j .(4.8) Le termeL ij correspond aux tensions de Leonard qui peuvent être évaluées directement et le termeR ij représente la contribution des tensions de sous-maille. Le système (4.6) ainsi obtenu est ouvert du point de vue des variables filtrées et une représentation de l’effet des échelles non résolues doit être introduite via une modélisation du tenseurR ij . Un exemple de fermeture - le modèle de Smagorinsky Le premier modèle de fermeture ou modèle de sous-maille, proposé parSmagorinsky(1963), relie le tenseur des contraintes non résoluesR ij aux grandeurs filtrées simulées comme suit : R ij − 1 3 R αα δ ij =−2ν t S ij oùν t = (C s Δ) 2 q 2S αβ S αβ .(4.9) oùδ ij est le symbole de Kronecker,C s est la constante de Smagorinsky. L’échelle de longueur caractéris- tique est celle du filtreΔet la vitesse caractéristique est estimée grâce au tenseur des taux de déformation moyens : S ij ( 2 1 =u i,j +u j,i ).(4.10) Ce modèle est le plus largement utilisé mais présente le défaut de dissiper trop d’énergie sans assurer sa redistribution dans le cas où les structures non résolues après la troncature spectrale ne sont pas strictement dissipatives 3 . Dans la littérature, de très nombreux modèles de sous-maille ont été développés pour pallier les limitations du modèle original de Smagorinsky. Pour ne citer que quelques exemples, le modèle mixte deBardinaet al.(1983) autorise une redistribution de l’énergie vers les échelles résolues, le modèle deSchumann(1975) inclut une équation de transport de l’énergie cinétique de sous-maille utilisée pour déterminer la viscosité de turbulenceν t . Plus récemment,Germanoet al.(1991) etGermano (1992) ont introduit le concept de modèle dynamique qui autorise une variation du coefficientC s . Ce 3 discrétisations spatiales théoriquement nécessaires ne peuvent pas êtreEn particulier dans les régions pariétales, où les mises en œuvre pour des raisons de coût numérique de résolution.

50