La turbulence est l'un des phénomènes les plus complexes observés dans la nature, ce qui rend difficile de définir avec précision. Dans un écoulement turbulent, le fluide suit des trajectoires courbes irrégulières appelées des tourbillons. Le flux est généralement entrelacé et crée des structures de nombreuses tailles différentes. Ils se déplacent et pivotent instantanément, interagissent entre eux et avec le champ d'écoulement principal et changent rapidement de forme et de taille. Le mélange est important et affecte la transmission de quantité de mouvement et, par conséquent, les forces aérodynamiques dans un obstacle fluide et environnant, comme les bâtiments. Si vous souhaitez étudier ce phénomène complexe et regarder sous le bonnet, nous vous recommandons cette Introduction à la turbulence. [1]
Les structures turbulentes provoquent dans le fluide le tourbillon, qui est souvent utilisé pour décrire la turbulence plutôt que la vitesse.
Le tourbillon est principalement généré aux limites des solides. Dans les couches de contour formées le long des limites de solide, la vitesse varie de zéro au niveau de la limite (condition de non-glissement) à une valeur qui n'est pas affectée par la limite, déterminée par le flux. Une turbulence se produit lorsque des instabilités, telles que la rugosité de la surface limite, provoquent une instabilité du tourbillon, soutenue par un nombre de Reynolds suffisamment élevé. La couche limite se décompose de la limite, le tourbillon et la turbulence sont ainsi balayés dans les zones de fluide éloignées des limites. Les grandes turbulences sont généralement anisotropes (par exemple, l'écoulement au-delà d'un cylindre provoque le détachement des tourbillons). Les turbulences provoquent des instabilités qui provoquent l'élargissement, la compression et la suppression des tourbillons. Les structures à écoulement cohérent se désintègrent rapidement en une masse de turbulence qui se développe à grande échelle. Les grandes turbulences deviennent plus petites jusqu'à ce qu'elles atteindre une taille où la dissipation de leur énergie cinétique due à la viscosité est importante. La perte d'énergie cinétique fait apparaître ces tourbillons. [3]
Dans le cas d'un fluide incompressible, le tourbillon respecte l'équation de transport.
Numerical Modélisation de Turbulence
Afin de saisir complètement la turbulence par modélisation numérique, les équations du mouvement pour le flux du fluide à toutes les échelles spatiales et temporelles doivent être résolues. Il n'existe pas de méthode universelle appropriée.
La méthode exacte de calcul du flux à l'aide des équations mentionnées ci-dessus pour toutes les échelles, appelée « Simulation numérique directe » (DNS), n'est pas applicable pour la CFD pratique en raison des coûts de calcul. Les ressources de calcul requises par le DNS dépassent de loin la capacité des supercalculateurs les plus puissants actuellement disponibles.
La « Simulation des grandes structures de la turbulence » (LES) utilise plutôt des schémas numériques précis tels que DNS pour les grandes échelles, tandis que pour les petites échelles, la modélisation de la turbulence est utilisée (dite de modélisation à l'échelle de sous-grille). Elle a de graves limitations dans les zones proches des parois, car l'effort de calcul requis pour la couche limite, où l'échelle de longueur de turbulence devient très faible, augmente rapidement. Cependant, pour les écoulements de cisaillement libre, où les grandes turbulences sont du même ordre que la couche de cisaillement et sont fortement anisotropes, les résultats peuvent être extrêmement fiables. Il est utile de résoudre des problèmes tels que le flux provoque des vibrations, etc.
Pour la plupart des problèmes de CFD pratiques, le coût de calcul du DNS et, dans une mesure inférieure, sont trop importants. La méthode d'équations « Reynolds-Aveated Navier-Stokes » (RANS) est beaucoup plus économique (voir le sous-chapitre {%}#RansModelfForTurbulence Modèles RANS pour la turbulence]]).
Pour des problèmes plus complexes, où des avantages des méthodes mentionnées ci-dessus sont requis, mais les coûts de calcul doivent rester raisonnables, les « méthodes d' hybride globales » peuvent être utilisées (voir le sous-chapitre {%}#GlobalHybridModelsForTurbulence Modèles hybrides globaux pour la turbulence ]]). Les méthodes hybrides globales sont basées sur la combinaison de méthodes LES et RANS qui les modifient en fonction du niveau de résolution. RANS est appliqué pour une partie de la couche limite et les grandes turbulences sont résolues en dehors de ces régions par les. Les modèles les plus populaires sont « Detached Eddy Simulation » (DES) ou « Delayed Detached Eddy Simulation » (DDES).
RANS Modèles pour Turbulence
Pour les flux stationnaires, RWIND 3 utilise le modèle de turbulence RANS. Les RANS sont basés sur la décomposition de Reynolds selon laquelle une variable d'écoulement est décomposée en composantes moyennes et fluctuantes. Lorsque la décomposition est appliquée aux équations de Navier-Stokes, un terme supplémentaire appelé « Tenseur des contraintes de Reynolds » apparaît et un système d'équations doit être « fermé ». Les niveaux des modèles de turbulence RANS sont liés au nombre d'équations différentielles ajoutées aux équations RANS afin de les « fermer ». [2]
Les modèles à deux équations k-ε a k-ω les plus populaires sont également disponibles dans RWIND 3. Le modèle de turbulence à une équation « Spalart-Allmaras » (SA) a été développé spécifiquement pour les flux aérodynamiques et est également souvent utilisé dans les méthodes d' hybride globales. Il est utilisé dans RWIND 3 Pro pour la modélisation de la turbulence dans les flux transitoires (voir le sous-chapitre [#GlobalHybridModelsForTurbulence Modèles hybrides globaux pour la turbulence])
k-ε Variation Model
Le modèle k-ε a été le premier modèle de turbulence à être largement utilisé pour une variété de flux dans la CFD. Elle est basée sur une analogie du mouvement aléatoire des turbulences dans un écoulement fluide turbulent avec des particules à l'échelle structurale, suggérée par Boursinesq. Il a introduit le concept de la viscosité de turbulence, qui n'est pas une propriété du fluide mais qui est proportionnelle à une échelle de vitesse et de longueur caractéristique de la turbulence. Les modèles sont requis pour représenter chacune de ces échelles. Le spectre de vitesse est représenté par l'énergie cinétique turbulente k, décrite par une équation de transport. L'équation k inclut un terme pour son taux de dissipation ε; une équation de transport pour ε fournit un modèle pour ce terme – qui représente également l'échelle de longueur de la turbulence. [2], [3]
Le modèle k-ε est robuste et calcule peu. Celle-ci n'est valide que pour les écoulements turbulents. Par conséquent, il convient aux itérations initiales et aux études paramétriques. Elle fonctionne mal pour les écoulements complexes présentant un gradient de pression extrême ou défavorable, des séparations et de fortes courbures de ligne. Il se comporte également difficile aux limites.
k-ω Variation Model
Le modèle k-ω « ferme » le système RANS par deux équations différentielles partielles pour k et ω, avec la première variable étant à nouveau l'énergie cinétique de turbulence et la seconde étant le taux de dissipation spécifique (de l'énergie cinétique turbulente k dans la température interne énergie). Son meilleur terme de dissipation donne au modèle k-ω un avantage par rapport au modèle k- ε dans la zone de la paroi. Il se caractérise également par de bonnes performances pour le cisaillement libre et les flux à faible nombre de Reynolds. Il convient mieux aux flux et à la séparation des couches limites complexes en aérodynamique externe (mais la séparation des flux est généralement prévue trop importante et précoce, ce qui nécessite donc une résolution de maillage élevée près de la paroi). Elle peut également être utilisée pour les écoulements de transition.
Les modèles à deux équations contiennent de nombreuses hypothèses et sont calibrés pour fonctionner uniquement selon les fonctionnalités bien connues des applications qu'ils sont conçus pour résoudre. Néanmoins, leur résistance a prouvé son utilisation et sont largement utilisés dans le secteur des calculs CFD.
Global Modèles hybrides pour Turbulence
L'objectif de ces modèles hybrides globaux est de bénéficier des avantages des modèles RANS et LES disponibles. La méthode RANS est appliquée pour une partie de la couche limite, où RANS aurait des coûts de calcul élevés, et le reste de l'écoulement avec de grandes turbulences est résolu par Les, où RANS ne peut pas bien modéliser des structures turbulentes anisotropes. En d'autres termes, les régions où l'échelle de longueur turbulente est inférieure à la dimension maximale de grille se voient attribuer le mode de solution RANS. Comme l'échelle de longueur turbulente dépasse la dimension de la grille, les régions sont résolues à l'aide du mode LES, réduisant ainsi considérablement les coûts de calcul tout en offrant certains des avantages de la méthode LES dans des régions séparées.
Spalart-Allmaras DDES Model
Pour l'analyse des flux transitoires (uniquement dans RWIND 3 Pro), un modèle hybride global « Spalart-Allmaras Delayed Detached Eddy Simulation » est utilisé, voir turbulence-des-spalart-allmaras-ddes.html Openfoam®.
L'amélioration principale de la « simulation Delayed Detached Eddy Simulation » (DDES) est d'inclure les informations de viscosité turbulente dans le mécanisme de basculement afin de retarder ce basculement dans les couches limites. Le système RANS est « fermé » par une équation de transport de la viscosité selon le « modèle Spalart-Allmaras » avec une échelle de longueur du modèle intégrée à la distance au mur.
du modèle turbulent à équation unique de Spalart-Allmaras qui résout l'équation de transport modélisée pour la viscosité turbulente du tourbillonνT. L'équation résout une vibration de type Spalart-Allmaras, ṽ. En termes simples, la variable ṽ est plus facile à calculer que νT directement, de sorte que la variable ṽ est d'abord calculée numériquement. La viscosité turbulente turbulente νT est ensuite calculée (corrigée) à l'aide de ṽ et enfinνT est ajoutée aux équations de quantité de mouvement pour fermer le système d'équations et peut être résolue. Une description détaillée est disponible https://cfdisrael.blog/2017/04/26/understand-the-spalart-allmaras-model/ici.
La turbulence est l'un des phénomènes les plus complexes observés dans la nature, ce qui rend difficile de définir avec précision. La littérature donne de nombreuses définitions, par exemple celle incluse dans [1] : « Un mouvement de fluide est décrit comme turbulent s’il est 3D, de rotation, d’continu, de flux très entier et de dissipation ». Si vous souhaitez étudier ce phénomène complexe et regarder sous le bonnet, nous vous recommandons cette Introduction à la turbulence.
Afin de capturer complètement la turbulence par modélisation numérique, il est nécessaire de résoudre les équations du mouvement pour le flux de fluide à toutes les échelles spatiales et temporelles. Cette approche est appelée « Simulation numérique directe » (DNS). Pour les applications industrielles, les ressources de calcul requises par le DNS dépassent de loin la capacité des supercalculateurs les plus puissants actuellement disponibles.
À la place, RWIND 3 utilise une technique différente, telle que la vitesse ou la pression décomposée en composantes moyennes (moyennes) et en composantes fluctuantes. En d'autres termes, les équations déterminantes du mouvement des fluides sont moyennées afin de supprimer les petites échelles, ce qui permet d'obtenir un ensemble d'équations modifiés qui est moins difficile à résoudre en calcul. Ces équations sont appelées « équations de Navier-Stokes à moyenne de Reynolds » (RANS).
Afin de résoudre les RANS dans RWIND 3, le modèle de turbulence k–ε [2] est utilisé, qui introduit deux équations de transport pour les propriétés de turbulence : La première est l'équation de transport de l'énergie cinétique turbulente k et la seconde équation détermine le transport de la vitesse de dissipation ε de k. Cette méthode est le modèle le plus utilisé et le plus testé pour les calculs CFD. Sa robustesse, son caractère économique et sa précision raisonnable pour un large éventail d'applications des flux turbulents expliquent sa popularité dans les simulations de flux industriels. De plus, RWIND 3 propose le modèle de turbulence k–ω comme alternative (voir cet article de Wikipédia).
La « Simulation des grandes structures de la turbulence » (LES) permet de résoudre les structures turbulentes à relativement grande échelle comme dans (DNS). Les structures à petite échelle, appelées échelles de sous-grille, sont modélisées.
Dans l'« Analyse des flux transitoires », une modification d'une équation « RANS », le modèle « Spalart-Allmaras Delayed Detached Eddy Simulation » est utilisé, voir Openfoam®. Ce modèle tente de traiter les zones proches des parois à l'aide de la méthode RANS et de traiter le reste de l'écoulement de la manière LES. En d'autres termes, les régions où l'échelle de longueur turbulente est inférieure à la dimension maximale de grille se voient attribuer le mode de solution RANS. Comme l'échelle de longueur turbulente dépasse la dimension de la grille, les régions sont résolues à l'aide du mode LES.