O foco principal do presente trabalho é estender uma metodologia numérica embasada no uso de uma técnica de refinamento adaptativo de malha (AMR - Adaptive Mesh Refinement) e no uso de esquemas temporais multipasso implícitos-explícitos (IMEX) a aplicações envolvendo escoamentos reativos com baixo número de Mach. Originalmente desenvolvida para escoamentos incompressíveis, a formulação euleriana daquela metodologia emprega as equações de Navier-Stokes como modelo matemático para descrever a dinâmica do escoamento e o Método da Projeção, baseado no divergente nulo da velocidade do escoamento, para tratar o acoplamento pressão-velocidade presente na formulação com variáveis primitivas. Tal formulação euleriana original é estendida para acomodar novas equações agregadas ao modelo matemático da fase contínua: conservação de massa, fração de mistura (para representar as concentrações de combustível e oxidante), e energia. Além disso, uma equação termodinâmica de estado é integrada ao modelo matemático estendido e é empregada juntamente com a equação de conservação de massa para produzir uma nova restrição (não nula desta vez) ao divergente do campo de velocidade. Assume-se que o escoamento ocorre a baixo número de Mach (hipótese principal). O Método de Diferença Finita é empregado na discretização espacial das variáveis eulerianas de estado, empregando-se uma malha AMR. As vantagens e dificuldades desta extensão são cuidadosamente investigadas e reportadas. Pela importância, do ponto de vista de aplicações práticas, alguns estudos numéricos preliminares envolvendo escoamentos incompressíveis turbulentos com sprays são realizados (as gotículas compõem a fase dispersa). Num primeiro momento, apenas sprays com gotículas inertes são considerados. Embora ainda apenas iniciais, tais estudos já se mostram importantes pois identificam com clareza, em primeira instância, algumas das dificuldades inerentes a serem enfrentadas ao se tratar dentro desta nova metodologia um conjunto relativamente grande de gotículas lagrangianas. No caso de escoamentos incompressíveis turbulentos com sprays, a integração temporal se dá com métodos IMEX para a fase contínua e com o Método de Euler Modificado para a fase dispersa. A turbulência, em todos os casos que a envolvem, é tratada pelo modelo de Simulação das Grandes Escalas (LES - Large Eddy Simulation). As simulações computacionais se dão em um domínio tridimensional, um parelelepípedo, e empregam uma extensão (resultante do presente trabalho) do código AMR3D, um programa de computador sequencial implementado em Fortran90, oriundo de uma colaboração de longa data entre o IME-USP e o MFLab/FEMEC-UFU (Laboratório de Dinâmica de Fluidos da Universidade Federal de Uberlândia). O processamento foi efetuado no LabMAP (Laboratório da Matemática Aplicada do IME-USP). / It is the main goal of the present work to extend a numerical methodology based on both the use of an adaptive mesh refinement technique (AMR) and the use of a multistep, implicit-explicit time-step strategy (IMEX) to applications involving low Mach number reactive flows. Originally developed for incompressible flows, the Eulerian formulation of that methodology employs the Navier-Stokes equations to model the flow dynamics and the Projection Method, based on the vanishing divergence of the velocity field, to tackle the pressure-velocity coupling present when using primitive variables. That Eulerian formulation is extended by adding a new set of equations to the original mathematical model, describing the various properties of the continuous phase: mass conservation, mixture fraction (to represent concentrations of fuel and oxidizer) and energy. Also, a thermodynamic equation of state is included into the extended mathematical model which is employed, along with the equation for the conservation of mass, to derive a new restriction (this time, different from zero) to the divergence of the velocity field. It is assumed that one is dealing with a low Mach number flow (the main hipothesis). The discretization in space employs the Finite Difference Method for the Eulerian variables on a AMR mesh. Advantages and difficulties of such an extension of the previous methodology are carefully investigated and reported. For its importance in the real-world applications, few preliminary numerical studies involving incompressible turbulent flows with sprays are performed (the droplets form what it is called the dispersed phase). Only sprays formed by inert droplets are considered. Even though initial yet, such studies are most important because they clearly identify, first hand, certain difficulties in handling relatively large sets of Lagrangian droplets in the context of this new AMR methodology. In the context of turbulent incompressible flows with sprays, the overall time-step scheme is given by IMEX methods for the continuous phase and by the Improved Euler Method for the dispersed phase. In all the cases in which it is considered, turbulence is modeled by the Large Eddy Simulation (LES) model. The computational simulations are held in a tridimensional domain given by a paralellepiped and all of them employ the extention (resulting of the present work) of the AMR3D code, a sequencial computer program implemented in Fortran90, whose origin is the collaborative work between IMEUSP and MFLab/FEMEC-UFU (Fluid Dynamics Laboratory, Federal University of Uberlândia). Computations were performed at LabMAP (Applied Mathematics Laboratory at IME-USP).

L'instabilité de Rayleigh-Taylor (IRT) est notamment rencontrée lors des expériences de Fusion par Confinement Inertiel, et son développement est un obstacle à la réussite de ces expériences. L'objet de cette thèse est d'étudier la croissance de l'IRT pour différents régimes de compressibilité, au moyen de simulations numériques directes réalisées à l'aide d'un code pseudo-spectral multidomaine de type Chebyshev-Fourier-Fourier.La méthode du développement asymptotique permet d'établir des modèles à bas nombre de Mach pour lesquels la contribution acoustique est éliminée. L'implantation dans le code de simulation du modèle anélastique, qui met en jeu des fluides stratifiés et capture les effets thermiques, est améliorée. Le modèle de Boussinesq est ajouté au code. La précision de la méthode numérique est étudiée pour différents découpages en sous-domaines. Plusieurs éléments de validation sont présentés, dont la comparaison avec une expérience.La première simulation présentée, réalisée avec le modèle de Boussinesq, s'intéresse à la croissance auto-semblable de l'IRT. Les lois d'échelle de la vorticité et de la dissipation sont dégagées. La structure de la turbulence et du mélange entre les deux fluides est discutée. Certaines propriétés de la turbulence homogène et isotrope sont retrouvées, mais on note la persistance d'anisotropie aux petites échelles. Les premières simulations 3D de l'IRT avec le modèle anélastique sont présentées. L'influence des effets de compressibilité sur les premières phases de la croissance est étudiée. En outre, une couche de mélange anélastique en faible stratification est analysée et présente des effets de compressibilité non négligeables. / The Rayleigh-Taylor instability (RTI) is especially observed in inertial confinement fusion experiments, and its development prevents the success of these experiments. The purpose of this work is to study the growth of the RTI for different compressibility regimes by using a multidomain pseudospectral Chebyshev-Fourier-Fourier simulation code. The asymptotic expansion method allows to establish several low Mach number models which do not contains acoustics. The implantation of the anelastic model, which deals with stratified fluids and captures thermal effects, has been improved. Moreover, the Boussinesq model is added to the simulation code. The accuracy of the entire numerical method is studied, as a function of the subdomain separation, and several validation elements are shown, including a comparison with an experimental study. The first simulation to be analyzed is achieved with the Boussinesq model. We focus on the self-similarity of the RTI growth. The temporal scalings of vorticity and dissipation are displayed, and the structures of turbulence and mixing are discussed. Some properties of isotropic and homogeneous turbulence are observed, however some anisotropy remains at small scales. The first three-dimensional anelastic simulations are presented. The influence of compressibility effects on the first stages of the growth is studied. Finally, a developed anelastic mixing layer involving weakly stratified fluids is described and was found to display non-negligible compressibility effects.

This work represents the initial steps in a wider project that aims to map out the sensitive areas in fuel injectors and combustion chambers. Direct numerical simulation (DNS) using a Low-Mach-number formulation of the Navier–Stokes equations is used to calculate direct-linear and adjoint global modes for axisymmetric low-density jets and lifted jet diffusion flames. The adjoint global modes provide a map of the most sensitive locations to open-loop external forcing and heating. For the jet flows considered here, the most sensitive region is at the inlet of the domain. The sensitivity of the global-mode eigenvalues to force feedback and to heat and drag from a hot-wire is found using a general structural sensitivity framework. Force feedback can occur from a sensor-actuator in the flow or as a mechanism that drives global instability. For the lifted flames, the most sensitive areas lie between the inlet and flame base. In this region the jet is absolutely unstable, but the close proximity of the flame suppresses the global instability seen in the non-reacting case. The lifted flame is therefore particularly sensitive to outside disturbances in the non-reacting zone. The DNS results are compared to a local analysis. The most absolutely unstable region for all the flows considered is at the inlet, with the wavemaker slightly downstream of the inlet. For lifted flames, the region of largest sensitivity to force feedback is near to the location of the wavemaker, but for the non-reacting jet this region is downstream of the wavemaker and outside of the pocket of absolute instability near the inlet. Analysing the sensitivity of reacting and non-reacting variable-density shear flows using the low-Mach-number approximation has up until now not been done. By including reaction, a large forward step has been taken in applying these techniques to real fuel injectors.

L'objectif de cette thèse est l'accélération des solveurs de Gradient Conjugué avec déflation utilisés pour la résolution de l'équation de Poisson pour la pression, dans le cas de la simulation d'écoulements à faible nombre de Mach sur des maillages non structurés. Une méthode de redémarrage basée sur une estimation de l'effet des erreurs numériques a été mise en œuvre et validée. Par la suite, une méthode à trois niveaux de maillage a été créée, et deux techniques ont dû être développées pour réduire le nombre d'itérations sur les niveaux grossiers : l'une permet la création de solutions initiales grâce à une méthode de projection adaptée, l'autre consiste en une adaptation du critère de convergence sur les niveaux grossiers. Les résultats numériques sur des simulations massivement parallèles montrent entre autres une réduction considérable du temps de calcul global. D'autres pistes de recherche sont introduites, notamment concernant l'équilibrage dynamiques de charge de calcul. / The main objective of this thesis is to accelerate deflated Conjugate Gradient solvers used for solving the pressure Poisson equation, for the simulation of low-Mach number flows on unstructured meshes. A restart method based on an estimation of the effect of numerical errors has been implemented and validated. Then, a three-level deflation method has been created, and two techniques are developed in order to reduce the number of iterations on the coarse levels : one of them is the creation of initial guesses thanks to a well-suited projection method, the other one consists in adapting the convergence criterion on the coarse grids. Numerical results on massively parallel simulations show, among others, a drastic reduction of the computational times of the solver. Other lines of research are introduced, especially regarding dynamic load balancing.

L'objectif de cette thèse est de développer un nouveau schéma numérique du type volumes finis pour la dynamique des gaz. Dans deux articles, F.Berthelin, T.Goudon et S.Minjeaud proposent de résoudre le système des équations d'Euler barotrope en dimension 1 d'espace, avec un schéma d'ordre 1 fonctionnant sur grilles décalées et dont la conception des flux est inspirée des schémas cinétiques. Nous proposons d'enrichir ce schéma afin qu'il puisse résoudre le système des équations d'Euler barotrope ou complet, en dimension 2 d'espace sur maillage cartésien ou non structuré, possiblement à l'ordre 2 et le cas échéant à bas nombres de Mach. Nous commencerons par développer une version 2D du schéma sur grilles cartésiennes (ou MAC) à l’ordre 2 via une méthode de type MUSCL, d'abord pour les équations barotropes puis pour les équations complètes. Ces dernières demandent de traiter une équation d’énergie supplémentaire et l’un des problèmes -résolu- est de trouver une définition discrète convenable de l’énergie totale telle qu'elle satisfasse une équation conservative locale. Dans un troisième chapitre nous étudierons le passage à la limite du compressible vers l'incompressible et nous verrons comment utiliser les atouts de notre schéma afin de le modifier et d'en faire un schéma Asymptotic Preserving pour des écoulements à bas nombres de Mach. Dans un quatrième temps nous proposerons une adaptation du schéma sur des maillages non structurés. Notre approche sera fortement inspirée des méthodes DDFV et pourra présenter des avantages dans les régimes à faibles nombres de Mach. Cette thèse se termine par un cinquième chapitre issu d’une collaboration lors du CEMRACS 2017, où le point de vue considéré n’est plus macroscopique mais microscopique. Nous commencerons par étudier un modèle micro/macro idéalisé auquel un processus stochastique a été ajouté puis nous tenterons d'en déduire un modèle à grande échelle pour un système fortement couplé, qui soit consistant avec la description micro/macro sous-jacente du problème physique. / The objective of this thesis is to develop a new numerical scheme of finite volume type for gas dynamics. In two articles, F.Berthelin, T.Goudon and S.Minjeaud propose to solve the barotropic Euler system in dimension 1 of space, with a first order scheme that works on staggered grids and of which fluxes are inspired by kinetic schemes. We propose to enhance this scheme so that it can solve the barotropic or complete Euler systems, in dimension 2 of space on Cartesian or unstructured grids, possibly at order 2 and at Low Mach numbers where appropriate. We begin with the development of a 2D version of the scheme on Cartesian (or MAC) grids, at order 2 via a MUSCL type method, for the barotropic equations at first and then for the complete equations. The latter require to handle with an additional energy equation and one of the -solved- problems is to find a suitable discrete definition of the total energy such that it satisfies a local conservative equation. In a third chapter we study the transition from the compressible case to the incompressible limit and we shall see how to use the advantages of our initial scheme in order to make it an Asymptotic Preserving scheme at low Mach numbers. In a fourth chapter we propose an adaptation of the scheme on unstructured meshes. Our approach is strongly inspired by the DDFV methods and may have advantages in low-Mach regimes.This thesis ends with a fifth chapter issued from a collaboration during CEMRACS 2017, where the considered point of view is no longer macroscopic but microscopic. We begin by studying a simplified micro/macro model with an added stochastic process and then we attempt to deduce a large-scale model for a strongly coupled system which has to be consistent with the underlying micro / macro description of the physical problem.

Ce travail de thèse porte sur l'extension et l'analyse d'un solveur volumes finis eulérien, co-localisé avec capture d'interfaces pour la simulation des écoulements multi-matériaux non miscibles. Les extensions proposées s'inscrivent dans la volonté d'élaborer un outil de simulation multi-physiques. Dans le cadre de ce mémoire, le caractère multi-physiques recouvre les champs que nous allons détailler. Nous traitons le cas des écoulements radiatifs modélisés par un système à deux températures qui couple les phénomènes purement hydrodynamiques aux phénomènes radiatifs. Nous proposons un solveur permettant la prise en compte des effets de tension superficielle à l'interface entre deux fluides. Nous développons un solveur implicite permettant la simulation précise d'écoulements faisant intervenir de faibles nombres de Mach par le biais d'une méthode de renormalisation de la diffusion numérique. Enfin, les effets tri-dimensionnels sont considérés ainsi que la possibilité d'étendre le schéma de base aux écoulements à un nombre quelconque de matériaux. A chaque étape, les solveurs développés sont validés sur des cas-tests. / This work is devoted to the extension of a eulerian cell-centered finite volume scheme with interfaces capturing for the simulation of multimaterial fluid flows. Our purpose is to develop a simulation tool which could be able to handle multi-physics problems in the following sense. We address the case of radiating flows, modeled by a two temperature system of equations where the hydrodynamics are coupled to radiation transport. We address a numerical scheme for taking surface tension forces into account. An implicit scheme is proposed to handle low Mach number fluid flows by means of a renormalization of the numerical diffusion. Eventually, the scheme is extended to three-dimensional flows and to multimaterial flows, that is with an arbitrary number of materials. At each step, numerical simulations validate our schemes.

Le modèle numérique que nous avons développé au cours de cette thèse présente deux caractéristiques principales : un modèle dilatable pour l'eau et la prise en compte de domaines ouverts. Les difficultés associées au premier aspect concernent l'adaptation de la loi d'état de l’eau au modèle dilatable sous l’approximation à faibles nombres de Mach, tandis que celles associées au second sont relatives à la mise en œuvre de conditions aux limites numériques de sortie compatibles avec l'algorithme de projection utilisé. Les résultats de simulations d'écoulement de convection mixte en canal horizontal chauffé par le bas ont été confrontés à celles utilisant l'approximation de Boussinesq et aux expériences. / The 3D numerical model which we developed in this thesis presents two main features: a Low-Mach-Number approximation for water along with an open boundary condition formulation. Indeed, the difficulties related to the former point stand in a computationally efficient adaptation of the water equation of state in the framework of Low Mach number approximation, whereas the difficulties related to the latter concern the introduction of Open Boundary Conditions in the projection algorithm used. We have computed a mixed convection flow in a horizontal channel uniformly heated from below and compared the results obtained with both the Boussinesq approximation and experimental results.

Dans cette thèse on s'intéresse aux jets en écoulement transverse dans une configuration générique de celle du refroidissement par effusion de chambres de combustion aéronautiques. L'amélioration des modèles de paroi avec transfert de masse passe par une meilleure connaissance de l'interaction entre les jets et l’écoulement principal. Nous avons donc réalisé la simulation numérique directe d'un jet issu d'un perçage incliné avec ou sans giration, pour des écoulements isothermes, turbulents et à bas nombre de Mach, dans un contexte compressible. Pour cela nous avons travaillé avec la bibliothèque AeroSol d'éléments finis continus et discontinus sur maillage hybride. En particulier nous nous sommes intéressés à la stabilité des flux numériques pour le compressible instationnaire associés à la méthode de Galerkin discontinue lorsque le nombre de Mach tend vers zéro. Nous avons pu mettre en évidence des comportements instables lors de l'utilisation de discrétisation temporelle explicite que nous avons corrigés en proposant un nouveau flux. Dans un deuxième temps, nous avons effectué les développements nécessaires à la réalisation des calculs. Nous nous sommes en particulier intéressés à la génération d'un champ de vitesse turbulent synthétique par la méthode SEM (Synthetic Eddy Method) que nous avons implantée dans AeroSol et validée. Grâce aux outils de post-traitement développés, nous avons conduit l'analyse de nos résultats. Dans le cas sans giration, les comparaisons avec les résultats expérimentaux et les résultats de simulations RANS que nous avons obtenus en parallèle sur la configuration du banc d'essai MAVERIC sont encourageants. La structure moyenne d'ensemble du jet est notamment correctement reproduite. En ce qui concerne la cas avec giration, le comportement attendu de déflexion successive du jet dans les deux plans caractéristiques (plan d'injection et plan de l'écoulement transverse) est bien reproduit et illustre tout le potentiel prévisionnel de la librairie de calcul que nous avons contribué à développer. / In this work we are interested in jet in crossflow in a generic configuration to the one used in effusion cooling for combustion chambers. Improved wall models with mass transfer requires a better knowledge of the interaction between the jets and the main flow. We therefore carried out the direct numerical simulation of a jet issuing from an inclined hole with or without gyration, for isothermal turbulent flow at low Mach number, in a compressible context. To achieved this, we worked with the continuous and discontinuous finite element library : AeroSol on hybrid grid. In particular we studied the stability of numerical flux for the unsteady compressible flow associated with discontinuous Galerkin method when the Mach number tends to zero. We were able to demonstrate unstable behavior when using explicit time discretization and we corrected them by providing a new flux. In a second time, we have performed the necessary development to achieve the calculations. We have been especially interested in the generation of a synthetic turbulent velocity field using the SEM method (Synthetic Eddy Method) that we have implemented in aerosol and validate. Thanks to the developed post-processing tools, we have conducted an analysis of our results. In the case without gyration, comparisons with experimental results and the results of RANS simulations we obtained on the Maveric test-bench configuration are encouraging. The mean flow of the jet is correctly reproduced. In the case with gyration, the expected behavior of successive deflection of the jet in both planes (injection plane and transverse plane of the flow) is reproduced and shows all the potential of the AeroSol library we helped to develop.

Cette thèse est consacrée à l'analyse mathématique de quelques modèles hétérogènes intervenants en mécanique des fluides. En particulier, elle est consacré a l'étude théorique des systèmes d'équations aux derivées partielles décrivants les trois modèles principaux que nous voulons présenter dans la suite. Le premier modèle étudié dans cette thèse est consacré à l'étude de l'écoulement d'un fluide visqueux newtonien et incompressible dans un bassin avec bathymétrie qui dégenère proche du bord. Le modèle mathématique étudié provient alors des équations de Navier-Stokes incompressible. On cherche à montrer que le problème de Cauchy correspondant est bien posé, au sens qu'on peut garantir l'existence globale et l'unicité de solutions faibles. Nous discuterons aussi la régularité de la solution faible. Finalement,nous établissons la convergence de la solution du modèle visqueux vers celle du modèle non visqueux quand le coefficient de viscosité tend vers zéro.La deuxième partie est dédiée à l'étude d'un modèle issu du système de Navier-Stokes dispersif ( il contient une correction dispersive) qui est lui meme obtenu à partir de la théorie cinétique des gaz. Notre modèle mathématique est dérivé a partir de ce dernier en supposant que le nombre de Mach est très faible. Le modèle obtenu est nommé effet fantôme (ou ghost effect an anglais), puisqu'il ne peut pas être obtenu à partir du modèle de Navier-Stokes compressible classique. L'objectif dans ce cadre sera d'étendre un résultat concernant l'existence locale d'une solution forte vers l'existence globale d'une solution faible. L'ingrédient principal dans notre analyse est une nouvelle inégalité fonctionnelle de type Log-Sobolev.La troisième partie de ce document est une contribution à une thématique de recherche se proposant d'analyser la compréhension des phénomènes rencontrés en géophysique qui font intervenir des milieux granulaires. Le modèle mathématique choisi est de type Bingham incompressible dont on fait dépendre le seuil de plasticité et le coe fficient de viscosité de la pression. On montre un résultat d'existence globale d'une solution faible du problème de Cauchy associé. / This thesis is devoted to the mathematical analysis of heterogeneous models raised by fluid mechanics. In particular, it is devoted to the theoretical study of partial differential equations used to describe the three main models that we present below.Initially, we are interested to study the motion of a compressible newtonienfluids in a basin with degenerate topography. The mathematical model studied derives from incompressible Navier-Stokes equations. We are interested to prove that the Cauchy problem associated is well posed. Well-posedness means that there exists a solution, that it is unique. In the meantime, we prove that the solution of the viscous model converges to the one of the inviscid limit model when the viscosity coe cient tends to zero.The second part in my thesis is devoted to study a model that arises from dispersive Navier-Stokes equations (that includes dispersive corrections to the classical compressible Navier-Stokes equations). Our model is derived from the last model assuming that the Mach number is very low. The obtained system is a Ghost eect system, which is so named because it can be derived from Kinetic theory. The main goal of this part is to extend a result concerning the local existence of strong solution to a global-in time existence of weak solutions. The main ingredient in this work is a new functional inequality of Log-Sobolev type.The last part of my thesis is a part of a research theme intends to analyze the understanding of phenomena encountered in geophysics which involves granular media. The mathematical model is of Bingham incompressible type with viscosity and placticity depending on the pressure. We provide a global existence of weak solutions of the Cauchy problem associated.

Dans les centrales solaires à tour, le flux solaire est concentré vers un récepteur solaire où son énergie est transférée à un fluide caloporteur. L'écoulement au sein du récepteur solaire est turbulent, fortement anisotherme et à bas nombre de Mach. L'optimisation du récepteur solaire exige une meilleure compréhension et modélisation de l'interaction entre la température et la turbulence. Cette thèse cherche à y contribuer selon deux approches. Tout d'abord, on étudie les échanges énergétiques entre les différentes parties de l'énergie totale. On propose pour cela une nouvelle représentation des échanges énergétiques, fondée sur la moyenne de Reynolds. Cette représentation permet la caractérisation, à partir de simulations numériques directes d'un canal plan bipériodique anisotherme, de l'effet du gradient de température sur les échanges énergétiques associées à l'énergie cinétique turbulente dans les domaines spatial et spectral. Ensuite, on étudie la simulation des grandes échelles des équations de bas nombre de Mach. En utilisant les résultats de simulations numériques directes, on identifie les termes sous-mailles spécifiques à modéliser lorsque l'on utilise le filtre classique, non pondéré, et lorsque l'on utilise le filtre de Favre, pondéré par la masse volumique. Dans les deux cas, on évalue a priori la performance de différents modèles sous-mailles. La pertinence des modèles est vérifiée a posteriori par la réalisation de simulation des grandes échelles. / In solar power towers, the solar flux is concentrated towards a solar receiver, wherethrough its energy is transferred to a heat transfer fluid. The flow in the solar receiver is turbulent, strongly anisothermal and at low Mach number. The optimisation of the solar receiver requires a better understanding and modelling of the interaction between temperature and turbulence. In this thesis, this is investigated following two approaches. First, we study the energy exchanges between the different parts of total energy. To this end, a new representation of the energy exchanges, based on the Reynolds averaging, is established. The representation allows the characterisation, from direct numerical simulations of a strongly anisothermal channel flow, of the effect of the temperature gradient on the energy exchanges associated with turbulence kinetic energy in the spatial and spectral domains. Second, we study the large-eddy simulation of the low Mach number equations. Using the results of direct numerical simulations, we identify the specific subgrid terms to model when the unweighted classical filter is used and when the density-weighted Favre filter is used. In both cases, the performance of different subgrid-scale models is assessed a priori. The relevance of the subgrid-scale models is then verified a posteriori by carrying out large-eddy simulations.