Développement spatial d'une couche de mélange plane 3D incompressible :  influences de la modélisation de sous-maille et des conditons d'entrée pour la simulation des grandes échelles

_____________________

S. Pellerin, A. Dulieu, L. Ta Phuoc, C. Tenaud

Figure

Objet

L'étude du développement des couches de mélange permet de mieux appréhender les phénomènes physiques se produisant par exemple en aval d'aéronefs ou dans les chambres de combustion. Les simulations des grandes échelles effectuées au LIMSI [1] ont permis de décrire la complexité d'une couche de mélange issue d'un bord de fuite d'une plaque plane séparant deux écoulements incompressibles, turbulents, à vitesses différentes. Ce type de simulation permet de prédire spatialement l'apparition et le développement des tourbillons dans la couche de mélange, en bon accord quantitatif avec l'expérience de référence, réalisée au CEAT-Poitiers. Au regard des premiers résultats de validations [1], deux aspects ont plus particulièrement été étudiés au cours de ce travail : (i) l'influence des conditions amont et (ii) l'importance de la modélisation de la contribution de sous-maille sur le développement et l'établissement de la couche de mélange.

Description

La simulation des grandes échelles (L.E.S.) du développement spatial de la couche de mélange est obtenue par résolution des équations de Navier-Stokes écrites sous forme vitesse-tourbillon et filtrées spatialement. Pour tenir compte de la contribution des échelles de sous-maille, le vecteur de sous-maille est modélisé grâce à la théorie du transfert de vorticité (Taylor, 1932), en utilisant une viscosité de sous-maille. Deux modèles de viscosité de sous-maille ont été utilisés pour relier cette viscosité au champ des quantités résolues : un modèle de vorticité et un modèle d'échelles mixtes (LIMSI). Pour tenir compte de l'anisotropie du maillage dans la définition de la longueur caractéristique du filtre, une fonction de Scotti a été utilisée. Du fait du caractère convectif des couches de mélange, les résultats de simulations dépendent des conditions à l'entrée du domaine de calcul. Une étude paramétrique de sensibilité aux perturbations amont a conjointement été menée.

Résultats et perspectives

La Figure 1 montre l'évolution de l'épaisseur de vorticité $\delta_{\omega}$ en fonction de la distance au bord de fuite de la plaque pour différentes perturbations. Un bruit blanc de 5 à 7,5 % de la vitesse supérieure donne un très bon accord avec l'expérience. Même si l'amplitude des profils est assez bien prédite, on observe que le modèle de vorticité est très dissipatif puisque les profils de tension longitudinale (Fig. 2) ne se superposent pas dans le plan transformé. Ce modèle ne permet pas de retrouver une zone de similitude. Un très bon accord avec les expériences est obtenu avec le modèle d'échelles mixtes sur les profils des tensions de Reynolds (Fig. 3). Ces profils tracés pour différentes abscisses dans une région s'étendant sur 5 $\delta_{\omega_0}$ ( $\delta_{\omega_0} = 30$ mm) montrent la qualité de la similitude retrouvée. De plus, une organisation spatiale et temporelle à grande échelle est bien marquée puisque cette L.E.S. prédit correctement des fréquences réduites (St = 0,3) et des longueurs d'onde caractéristiques (Figs. 4 et 5) en accord à la fois avec l'expérience et la théorie de Barnal et Roshko (1986). Une analyse modale (P.O.D.) a confirmé le bon accord expérience - calculs, en ce qui concerne les directions privilégiées du tenseur de Reynolds. Nous envisageons d'étudier l'influence d'un cisaillement transversal sur le développement d'une couche cisaillée libre turbulente, représentative de l'intersection entre deux jets coaxiaux en rotation. Références

[1] Lardat R., Dulieu A., Ta Phuoc L., and Tenaud C. : ``L.E.S. of a spatially developing 3-D incompressible mixing layer with velocity vorticity formulation.'', Lectures Notes in Physics, vol. 515, C.-H. Bruneau (ed.), Springer, pp. 183-188, 1998.
[2] S. Pellerin, A. Dulieu, L. Ta Phuoc and C. Tenaud : ``Incompressible 3-D mixing layer using LES: Influences of subgrid scale models and of upstream perturbations.'', ISCFD'99, Bremen, Germany, September 1999.
Article à paraître dans CFD Journal : Vol. 9, 1, April 2000.

Gpe Dynamique des Fluides

Dpt Mécanique

+ Sommaire

Présentation