La prévention des risques industriels nécessite de simuler la dispersion turbulente de polluants. Cependant, les outils majoritairement utilisés à ce jour ne permettent pas de traiter les champs proches dans le cas de géométries complexes, et il est nécessaire d'utiliser les outils de CFD (“ Computational Fluid Dynamics ”) plus adaptés, mais plus coûteux. Afin de simuler les écoulements atmosphériques avec dispersion de polluants, les modèles CFD doivent modéliser correctement d'une part, les effets de flottabilité, et d'autre part les effets de la turbulence. Plusieurs approches existent, notamment dans la prise en compte des effets de flottabilité et la modélisation de la turbulence, et nécessitent des méthodes numériques adaptées aux spécificités mathématiques de chacune d'entre elles, ainsi que des schémas numériques précis pour ne pas polluer la modélisation. Une formulation d'ordre élevé en volumes finis, sur maillages non structurés, parallélisée, est proposée pour simuler les écoulements atmosphériques avec dispersion de polluants. L'utilisation de schémas d'ordre élevé doit permettre d'une part de réduire le nombre de cellules et diminuer les temps de simulation pour atteindre une précision donnée, et d'autre part de mieux contrôler la viscosité numérique des schémas en vue de simulations LES (Large Eddy Simulation), pour lesquelles la viscosité numérique des schémas peut masquer les effets de la modélisation. Deux schémas d'ordre élevé ont été étudiés et implémentés dans un solveur 3D Navier Stokes incompressible sur des maillages volumes finis non structurés. Nous avons développé un premier schéma d'ordre élevé, correspondant à un schéma Padé volumes finis, et nous avons étendu le schéma de reconstruction polynomiale de Carpentier (2000) aux écoulements incompressibles. Les propriétés numériques des différents schémas implémentés dans le même code de calcul sont étudiées sur différents cas tests bi-dimensionnels (calcul de flux convectifs et diffusifs sur une solution a-priori, convection d'une tâche gaussienne, décroissance d'un vortex de Taylor et cavité entraînée) et tri-dimensionnel (écoulement autour d'un obstacle cubique). Une attention particulière a été portée à l'étude de la précision et du traitement des conditions limites. L'implémentation proposée du schéma polynomial permet d'approcher, pour un maillage identique, les temps de simulation obtenus avec un schéma décentré classique d'ordre 2, mais avec une précision supérieure. Le schéma compact donne la meilleure précision. En utilisant une méthode de Jacobi sans calcul implicite de la matrice pour calculer le gradient, le temps de simulation devient intéressant uniquement lorsque la précision requise est importante. Une alternative est la résolution du système linéaire par une méthode multigrille algébrique. Cette méthode diminue considérablement le temps de calcul du gradient et le schéma Padé devient performant même pour des maillages grossiers. Enfin, pour réduire les temps de simulation, la parallélisation des schémas d'ordre élevé est réalisée par une décomposition en sous domaines. L'assemblage des flux s'effectue naturellement et différents solveurs proposés par les librairies PETSC et HYPRE (solveur multigrille algébrique et méthode de Krylov préconditionnée) permettent de résoudre les systèmes linéaires issus de notre problème. / The prevention of industrial risks requires simulating turbulent dispersion of pollutants. However, the tools mostly used so far do not allow near fields treated in the case of complex geometries, and it is necessary to utilize the tools of CFD (Computational Fluid Dynamics ") more suitable but more expensive. To simulate atmospheric flows with dispersion of pollutants, the CFD models must correctly model the one hand, the effects of buoyancy, and secondly the effects of turbulence. Several approaches exist, including taking into account the effects of buoyancy and turbulence modeling, and require numerical methods adapted to the specific mathematics of each, and accurate numerical schemes to avoid pollution modeling. A formulation of high order finite volume on unstructured meshes, parallelized, is proposed to simulate the atmospheric flows with dispersion of pollutants. The use of high order schemes allow one hand to reduce the number of cells and decrease the simulation time to achieve a given accuracy, and secondly to better control the viscosity numerical schemes for simulation LES (Large Eddy Simulation), for which the numerical viscosity patterns may mask the effects of modeling. Two high-order schemes have been studied and implemented in a 3D Navier Stokes solver on unstructured mesh finite volume. We developed the first high-order scheme, corresponding to a Padé finite volume scheme, and we have extended the scheme of reconstruction polynomial Carpentier (2000) for incompressible flows. The numerical properties of the various schemes implemented in the same computer code are studied different two-dimensional test cases (calculation of diffusive and convective flow on a solution a priori, a task Gaussian convection, decay of a vortex of Taylor and driven cavity) and tri-dimensional (flow past an obstacle cubic). Particular attention has been paid to the study of the accuracy and treatment of boundary conditions. The implementation of the polynomial allows to obtain quasi identical simulation time compared to a classical upwind scheme of order 2, but with higher accuracy. The compact layout gives the best accuracy. Using a Jacobi method without calculation implied matrix to calculate the gradient, the simulation time becomes interesting only when the required accuracy is important. An alternative is the resolution of linear system by an algebraic multigrid method. This method significantly reduces the computation time of the gradient and the Padé scheme is effective even for coarse meshes. Finally, to reduce simulation time, the parallelization schemes of high order is achieved by a decomposition into subdomains. The assembly flow occurs naturally and different solvers provided by PETSc libraries and HYORE (algebraic multigrid solver and preconditioned Krylov method) used to solve linear systems from our problem. The work was to identify and determine the parameters that lead to lowest time resolution simulation. Various tests of speed-up and scale-up were used to determine the most effective and optimal parameters for solving linear systems in parallel from our problem. The results of this work have been the subject of a communication in an international conference "Parallel CFD 2008" and an article submitted to "International Journal for Numerical Methods in Fluids" (Analysis of high-order finite volume schemes for the incompressible Navier Stokes equations)

Dans cette thèse, nous présentons une nouvelle approche pour la simulation numérique des phénomènes électroconvectifs et électro-therrno-convectifs. La principale difficulté réside dans la détermination du champ électrique et de la distribution de densité volurnique de charges électriques. Dans cette approche, des schémas de type TVD (Total Variation Dirninishing) et IDC (ImprovedDeferred Correction) sont utilisées dans la discrétisation des flux convectifs et diffusifs par la méthode des volumes finis. La première partie de cette thèse présente certains aspects numériques liés à l'implémentation de ces schémas. Une approche unifiée pour les schémas convectifs TVD de type borné à haute résolution est présentée et diverses fonctions limiteur sont comparées. Dans une deuxième partie, l'électro-convection entre deux plaques parallèles est simulée. La méthodologie a étéévaluée et validée par la détermination des Critères de stabilité linéaire et non linéaire. Les différents scenarii d'évolution du développement de cette instabilité électroconvective vers l'état chaotique ont été définis. L'effet du mécanisme de diffusion la densité volumique de charge sur la boucle d'hystérésis et sur la structure de l'écoulement est étudié. L'influence du rapport d'aspect de la cavité est analysé. Enfin dans une dernière partie, nous étudions l'électro-thermo-convection lorsque le fluide est soumis simultanément à une injection unipolaire et à un gradient thermique. L'augmentation des transferts de chaleur a été caractérisée. / In this thesis, a numerical approach is presented to simulate the electro- and electro-thermo convection in dielectric liquids. The total variation diminishing (TVD) scheme and improved deferred correction (IDC) scheme are used to compute the convective and diffusive respectively. The aim of TVD scheme is to avoid non-physical oscillations and to capture high gradient of charge density. Some fundarnental aspects related to TVD and LDC schemes are investigated firstly. A unified approach for TVD schemes is explained and various limiter functions are compared. The connection among three methods for diffusive flux computation has been revealed. The original IDC scheme is improved by the application of 2nd order gradient evaluation method.The electro-convection between two parallel plates is then simulated. The methodology was assessed by the determination of the linear and nonlinear stability criterion. By continuously increasing the driving parameter, the successive instabilities and route to chaotic state has been defined. The effects of the diffusion mechanism for the charge density and vertical walls on the hysteresis 100p and the structure are also investigated. The last part is to simuiate electro-thermo-convection when injection and thermal gradient are simultaneously applied. Our solver was verified with a stationary and an overstable stability problem.The case that both heating and injection are from a bottom electrode has been analyzed in details. The neutral stabiliïy curve was reproduced. The existence of nonlinear phenornenon and the structure are highlighted.

This thesis presents the author's contribution <p>to the design and implementation of COOLFluiD,<p>an object oriented software platform for <p>the high performance simulation of multi-physics phenomena on unstructured grids. In this context, the final goal has been to provide a reliable tool for handling high speed aerothermodynamic <p>applications. To this end, we introduce a number of design techniques that have been developed in order to provide the framework with flexibility<p>and reusability, allowing developers to easily integrate new functionalities such as arbitrary mesh-based data structures, numerical algorithms (space discretizations, time stepping schemes, linear system solvers, ),and physical models. <p>Furthermore, we describe the parallel algorithms <p>that we have implemented in order to efficiently <p>read/write generic computational meshes involving <p>millions of degrees of freedom and partition them <p>in a scalable way: benchmarks on HPC clusters with <p>up to 512 processors show their effective suitability for large scale computing. <p>Several systems of partial differential equations, <p>characterizing flows in conditions of thermal and <p>chemical equilibrium (with fixed and variable elemental fractions)and, particularly, nonequilibrium (multi-temperature models) <p>have been integrated in the framework. <p>In order to simulate such flows, we have developed <p>two state-of-the-art flow solvers: <p>1- a parallel implicit 2D/3D steady and unsteady cell-centered Finite Volume (FV) solver for arbitrary systems of PDE's on hybrid unstructured meshes; <p>2- a parallel implicit 2D/3D steady vertex-centered Residual Distribution (RD) solver for arbitrary systems of PDE's on meshes with simplex elements (triangles and tetrahedra). <p>The FV~code has been extended to handle all <p>the available physical models, in regimes ranging from incompressible to hypersonic. <p>As far as the RD code is concerned, the strictly conservative variant of the RD method, denominated CRD, has been applied for the first time in literature to solve high speed viscous flows in thermochemical nonequilibrium, yielding some preliminary outstanding results on a challenging double cone flow simulation. <p>All the developments have been validated on real-life testcases of current interest in the aerospace community. A quantitative comparison with experimental measurements and/or literature has been performed whenever possible. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

This thesis describes analysis and optimization of existing heating system aboveground fire sprinkler water tanks. The purpose is to design a heating system so as to achieve optimal operating and economic conditions for the operation of heating in winter. I will use available computer software and adhere to the applicable regulations and standards. Using experi-mental methods and computer simulations I will assess the current state of the heating system which in turn will provide an optimizing its use and operation.

Doctor of Philosophy / Department of Mechanical and Nuclear Engineering / Daniel V. Swenson / The objective of the thesis is to provide a fully coupled thermo-hydro-mechanical (THM) tool, T2STR, which enables quantitative understanding and prediction of thermal as well as mechanical effects on flow in the porous media under multiphase conditions. This is achieved by incorporating a finite element based hydro-thermo-mechanical stress capability into the well-established IFDM (Integrated Finite Difference Method) based flow simulation code TOUGH2. TOUGH2 is a program for calculation of multi-phase, multi-component, non-isothermal flow in porous media. It implements several equation of state modules to represent different fluid mixtures. The dual mesh technique is natural for combining both discretization methods and is used innovatively and effectively. A generalized approach is developed to accommodate the switching of variables implemented in TOUGH2 to adapt the phase changes. The forward coupling is achieved by using the thermal, hydrostatic, and poroelastic effects in the stress calculations. The backward coupling includes the effect of strain on the fluid flow. T2STR also allows the user to study the variation in porosity, permeability and capillary pressure as function of mean effective stress in the porous media. Multiple materials can be used to model the reservoir in T2STR, parallel to the implementation in TOUGH2. T2STR is implemented to carry out as a fully coupled, one way coupled (only deformation as function of hydro-thermal effects), or original TOUGH2 implementation. It provides the ability to switch on and off the thermal and/or poroelastic effects. T2STR is developed to model the fractured porous media using discrete fractures. The modeling of fractured porous media is limited to a staggered coupling approach. The fluid parameters like permeability, porosity are modified based on the stresses and/or aperture changes due to deformation. A set of verification problems, used to validate the code and display the capabilities of the code, are discussed. A graphical user interface is designed to pre-process the necessary data. Macros are developed for excel and Tecplot to post-process the results for easy visualization.

L’émergence de bruits auparavant inaudibles dans les réservoirs à carburants automobiles requiert des constructeurs une meilleure compréhension des phénomènes physiques intervenants au sein de leurs produits. Dans cette thèse, différents travaux ont été conduits autour de l’étude du ballottement de fluide dans une cuve rigide rectangulaire partiellement remplie de fluide et soumise à une excitation extérieure. La première partie présente un état de l’art sur le sloshing suivant trois approches complémentaires - approche analytique, approche numérique et approche expérimentale - permettant d’orienter les travaux. Dans une deuxième partie, une étude préliminaire sur le sloshing dans une cuve rectangulaire soumise à une excitation harmonique forcée est réalisée. La confrontation des résultats numériques entre une approche linéaire - basée sur la théorie d’écoulement potentiel tenant compte de la viscosité du fluide [Schotté et Ohayon, 2013] - et une approche non linéaire commerciale – basée sur la résolution des équations de Navier-Stokes - permet de définir un paramètre de linéarité. Ce dernier permet de déterminer les cas de sloshing qui nécessitent une résolution non linéaire et ceux pour lesquels la théorie linéaire suffit pour prédire le phénomène. La troisième partie de ce document présente une étude expérimentale du ballottement de fluide dans une cuve rectangulaire rigide soumise à un freinage automobile. Deux niveaux de remplissage créant deux types d’impacts contre les parois (avec et sans enfermement de poche d’air) ont été analysés. Les essais menés ont permis de mesurer les forces engendrées par le mouvement du fluide, les pressions d’impact en paroi ainsi que le champ de vitesse par méthode Particle Image Velocimetry (PIV). Ce chapitre constitue une importante base de données expérimentales ayant permis d’étudier précisément le phénomène physique. L’étude est complétée par une confrontation des résultats expérimentaux avec des résultats Computational Fluid Dynamics (CFD). Enfin, pour conclure ce mémoire, une étude du sloshing dans un réservoir en tenant compte de la Fluid-Structure Interaction (FSI) est présentée. Le choix du couplage a été porté sur un schéma partitionné itératif faible avec, dans un premier temps, une approche potentielle instationnaire, puis avec une approche Volume Of Fluid (VOF) pour la physique fluide. Les limites d’un tel couplage dans le cas d’étude d’un réservoir partiellement rempli de fluide et attaché de manière flexible en fonction du rapport de masse fluide-réservoir ont été mises en évidence. La correction du schéma de couplage par l’effet de masse ajoutée présentée dans [Song et al., 2013] permet la résolution d’un système couplé quel que soit le rapport de masse en jeu et améliore de manière significative la convergence en réduisant également fortement le temps de calcul. / The present thesis focuses on an investigation of the sloshing phenomenon in a partially filled fuel tank submitted to a harmonic excitation motion. In the first part, the confrontation of numerical results between a linear approach - taking into account viscosity - and a nonlinear approach based on a commercial code leads to define a parameter of linearity. This parameter allows determining cases of sloshing who require non-linear resolution and those who need a linear theory to predict the phenomenon. An experimental study of fluid sloshing in a rectangular tank submitted to an automotive braking is conducted. Tests leaded allow measuring global forces engendered by the motion of the fluid, pressure of fluid impact and velocity field by PIV. This chapter provides an important data base and helps to investigate on the physical phenomenon. This study is completed by CFD results. To conclude, a numerical model for fluid-structure interactions is presented. Limits of this segregated partitioned coupling in case of sloshing in tank flexibly attached are highlighted, depending mostly on the mass ratio between fluid and tank structure. An added-mass term is integrated to the corrected staggered scheme ensuring systematically the convergence of the coupled solution and reducing significantly the iterations required.

Dans ce travail, nous étudions un schéma élément fini que nous appliquons à un modèle décrivant l'intrusion saline dans les aquifères côtiers confinés et libres. Le modèle est basé sur l'approche hydraulique qui consiste à moyenner verticalement le problème initial 3D, cette approximation repose sur une hypothèse d'écoulement quasi-hydrostatique qui, loin des épontes et des sources, est vérifiée. Pour modéliser les interfaces entre l'eau douce et l'eau salée (respectivement entre la zone saturée et la zone sèche), nous combinons l'approche 'interface nette' à l'approche avec 'interface diffuse' ; cette approche est déduite de la théorie de champ de phase, introduite par Allen-Cahn, pour décrire les phénomènes de transition entre deux zones. Compte tenu de ces approximations, le problème consiste en un système fortement couplé d'edps quasi-linéaires de type parabolique dans le cas des aquifères libres décrivant l'évolution des profondeurs des 2 surfaces libres et de type elliptique-prabolique dans le cas des aquifères confinés, les inconnues étant alors la profondeur de l'interface eau salée/eau douce et la charge hydraulique de l'eau douce. Dans la première partie de la thèse, nous donnons dans le cas d'un aquifère confiné, des résultats d'estimation d'erreur d'un schéma semi-implicite en temps combiné à une discrétisation en espace de type élément fini Pk Lagrange. Ce résultat utilise entre autre un résultat de régularité du gradient de la solution exacte dans l'espace Lr(ΩT), r > 2, ce qui permet de traiter la non-linéarité et d'établir l'estimation d'erreur sous des hypothèses de régularité raisonnables de la solution exacte. Dans la seconde partie de la thèse, nous généralisons l'étude précédente au cas de l'aquifère libre. La difficulté principale est liée à la complexité du système d'edps paraboliques mais à nouveau, grâce au résultat de régularité Lr(ΩT), r > 2 établi pour les gradients des surfaces libres, nous montrons que le schéma est d'ordre 1 en temps et k en espace pour des solutions suffisamment régulières. Nous concluons ce travail par des simulations numériques dans différents contextes (impact de la porosité et de la conductivité hydraulique sur l'évolution de l'interface, pompage et injection d'eau douce, effet des marées) validant ainsi le modèle et le schéma. Puis nous comparons les résultats à ceux obtenus avec un schéma volume fini construit à partir d'un maillage structuré. / In this work, we study a finite element scheme we apply to a model describing saltwater intrusion into coastal aquifers confined and free. The model is based on the hydraulic approach of vertically averaging the 3D original problem, this approximation is based on a quasi-hydrostatic flow hypothesis which, instead of the walls and springs, is checked. To model the interface between freshwater and salt water (respectively between the saturated zone and dry zone), we combine the approach net interface (approach with the diffuse interface) ; This approach is derived from the phase field theory introduced by Allen-Cahn, to describe the phenomena of transition between two zones. Given these approximations, the problem consists of a strongly couple to edps parabolic quasi-linear system in the case of unconfined aquifers describing the evolution of the depths of two free surfaces and elliptical-parabolic type in the case confined aquifer, the unknowns being then the depth of salt water / fresh water and the hydraulic load of fresh water. In the first part of the thesis, we give in the case of a confined aquifer, error estimation results of a semi-implicit scheme in a combined time discretization space finite element type Pk Lagrange. This result among other uses a regularity result of the gradient of the exact solution in the space Lr(ΩT), r > 2, which can handle the non-linearity and to establish the error estimate under assumptions reasonable regularity of the exact solution. In the second part of the thesis, we generalize the previous study to the case of the free aquifer. The main difficulty is related to the complexity of the system of parabolic edps but again, thanks to regularity result Lr(ΩT), r > 2 gradients established for the free surfaces, we show that the scheme is of order 1 time and space k for sufficiently regular solutions. We conclude this work by numerical simulations in different contexts (impact of porosity and hydraulic conductivity of the evolution of the interface, and pumping fresh water injection, tidal effects) thus validating the model and diagram. The we compare the results with those obtained using a finite volume scheme constructed from a structured mesh.

L’adaptation de maillage est un outil puissant pour l’obtention de simulations aérodynamiques précises à coût limité. Dans le cas particulier des simulations visant au calcul de fonctions aérodynamiques (efforts, moments, rendements...), plusieurs méthodes dites de raffinement ciblé (ou, en anglais, « goal-oriented ») basées sur le vecteur adjoint de la fonction d’intérêt ont été proposées. L’objectif de la thèse est l’extension d’une méthode de ce type basée sur la dérivée totale dJ/dX de la grandeur aérodynamique d’intérêt, J, par rapport aux coordonnées du maillage volumique, X. Les trois méthodes usuelles de calcul de gradient discret – la méthode de différentiation directe, la méthode adjointe-"paramètres" et la méthode adjointe-"maillage" évaluant dJ/dX – ont tout d’abord été étudiées et codées dans le logiciel elsA de l’ONERA pour des maillages non-structurés, pour des écoulements compressibles de fluide parfait et des écoulements laminaires. La seconde étape du travail a consisté à créer un senseur local θ basé sur dJ/dX qui identifie les zones du maillage volumique où la position des nœuds a une forte incidence sur l’évaluation de la fonction J. Ce senseur sert d’indicateur pour l’adaptation de différents maillages, pour différents régimes d’écoulement (subsonique, transsonique, supersonique), pour des configurations d’aérodynamique interne (aube et tuyère) et externe (profil d’aile). La méthode proposée est comparée à une méthode de raffinement ciblée très populaire (Venditti et Darmofal, 2001) et à une méthode de raffinement basée sur les caractéristiques de l’écoulement (ou, en anglais, « feature-based ») ; elle conduit à des résultats très satisfaisants. / Mesh adaptation is a powerful tool to obtain accurate aerodynamic simulations with limited cost. In the specific case of computation of aerodynamic functions (forces, moments, efficiency ...), goal-oriented methods based on the adjoint vector have been proposed. The aim of the thesis is the extension of a method of this type based on the total derivative dJ/dX of the aerodynamic output of interest, J, with respect to the volume mesh coordinates, X. The three common methods for calculating discrete gradient – the direct differentiation method, the parameter-adjoint method and mesh-adjoint method evaluating dJ/dX – have been studied first and coded in the elsA ONERA software for unstructured grids, for compressible inviscid and laminar flows. The second part of this work was has been to define a local sensor θ based on dJ/dX in order to identify zones where the volume mesh nodes position has a strong impact on the evaluation of the function J. This sensor is the selected indicator for different mesh adaptations for different flow regimes (subsonic, transonic, supersonic) for internal (blade and nozzle) and external (wing profile) aerodynamic configurations. The proposed method is compared to a well-known goal-oriented method (Darmofal and Venditti, 2001) and to a feature-based method ; it leads to very consistent results. very consistent results.

Ce manuscrit de thèse porte sur l'analyse numérique de schémas volumes finis pour la discrétisation de deux systèmes particuliers d'équations. Dans un premier temps nous étudions l'équation de Cahn-Hilliard associée à des conditions aux limites dynamiques dont l'une des principales difficultés est que cette condition aux limites est une équation parabolique, non linéaire, posée sur le bord et couplée avec l'intérieur du domaine. Nous proposons une discrétisation de type volumes finis en espace qui permet de coupler naturellement l'équation dans le domaine et celle sur sa frontière par un terme de flux et qui s'adapte facilement à la géométrie courbe du domaine. Nous montrons l'existence et la convergence des solutions discrètes vers une solution faible du système. Dans un second temps nous étudions la stabilité Inf-Sup du problème de Stokes pour un schéma volumes finis de type dualité discrète (DDFV). Nous donnons une analyse complète de la stabilité Inf-Sup inconditionnelle dans certains cas et de la stabilité de codimension 1 dans le cas de maillages cartésiens. Nous mettons également en place une méthode numérique permettant de calculer la constante Inf-Sup associée à ce schéma pour un maillage donné. On peut ainsi observer le comportement stable ou instable selon les cas en fonction de la géométrie des maillages. Dans une dernière partie nous proposons un schéma DDFV pour un modèle couplé Cahn-Hilliard/Stokes ce qui nécessite l'introduction de nouveaux opérateurs discrets. Nous démontrons la décroissance de l'énergie au niveau discret ainsi que l'existence d'une solution au problème discret. L'ensemble de ces travaux est validé par de nombreux résultats numériques. / This manuscript is devoted to the numerical analysis of finite-volume schemes for the discretization of two particular equations. First, we study the Cahn-Hilliard equation with dynamic boundary conditions whose one of the main difficulties is that this boundary condition is a non-linear parabolic equation on the boundary coupled with the interior of the domain. We propose a spatial finite-volume discretization which is well adapted to the coupling of the dynamics in the domain and those on the boundary by the flux term. Moreover this kind of scheme accounts naturally for the non-flat geometry of the boundary. We prove the existence and the convergence of the discrete solutions towards a weak solution of the system. Second, we study the Inf-Sup stability of the discrete duality finite volume (DDFV) scheme for the Stokes problem. We give a complete analysis of the unconditional Inf-Sup stability in some cases and of codimension 1 Inf-Sup stability for Cartesian meshes. We also implement a numerical method which allows us to compute the Inf-Sup constant associated with this scheme for a given mesh. Thus, we can observe the stable or unstable behaviour that can occur depending on the geometry of the meshes. In a last part we propose a DDFV scheme for a Cahn-Hilliard/Stokes phase field model that required the introduction of new discrete operators. We prove the dissipation of the energy in the discrete case and the existence of a solution to the discrete problem. All these research results are validated by extensive numerical results.

Ce travail est une contribution à la simulation numérique d'écoulements turbulents dans un but d'application industrielle. Nous nous intéressons dans un premier temps à une nouvelle combinaison "VMS-LES/procédure dynamique" pour la simulation d'écoulements autour de cylindres circulaire et carré. L'approche VMS-LES mise en œuvre a pour originalité d'utiliser une procédure de moyennage sur des volumes finis agglomérés dans le but de séparer les échelles, l'approche dynamique étant celle introduite par Germano en LES. Une approche hybride RANS/VMS-LES est ensuite évaluée sur le problème du cylindre circulaire à des nombres de Reynolds élevés. Cette approche introduit un paramètre d'hybridation qui selon la résolution locale de grille privilège le modèle RANS ou celui VMS-LES. Enfin, les performances d'un algorithme de Schwarz deux-niveau, qui utilise les méthodes de déflation et de balancing, sont examinées d'un point de vue efficacité et scalabilité dans le cas de simulations VMS-LES. / This work is a contribution to the numerical simulation of turbulent flows with the aim of industrial application. At first, we focus on a new combination "VMS-LES/dynamic procedure" for the simulation of flows around circular and square cylinders. The VMS-LES approach adopted in this work is original in using an averaging procedure over agglomerated finite volumes in order to separate the scales, the dynamic approach being the one introduced by Germano in LES. A RANS/VMS-LES hybrid approach is then evaluated on the circular cylinder test case at high Reynolds numbers. This approach introduces a hybridization parameter which privileges the RANS model or the VMS-LES model according to the grid resolution. Finally, the performance of a two-level Schwarz algorithm, which uses the deflation and balancing methods, are examined in terms of efficiency and scalability in the context of VMS-LES simulations.