Cette thèse traite plusieurs aspects scientifiques inhérents à la simulation numérique de problèmes d'interaction fluide-structure impliquant de fines membranes déformables. Deux cas spécifiques relatifs à la biomécanique cardiovasculaire sont considérés : l'interaction de l'écoulement sanguin avec la valve aortique (qui se produit à l'échelle macroscopique), et l'interaction de la membrane des globules rouges avec ses fluides interne et externe (qui se produit à l'échelle microscopique). Dans les deux cas, le couplage fluide-structure est géré par l'intermédiaire d'un formalisme de frontières immergées, en représentant la membrane par un maillage Lagrangien se mouvant au travers d'un maillage fluide Eulérien. Lorsque l'on traite la dynamique des globules rouges, la membrane est considérée comme étant une structure sans masse et infiniment fine. La première question à laquelle on s'intéresse dans cette thèse est la manière de modéliser la microstructure complexe de la membrane des globules rouges. Un moyen possible pour caractériser un modèle de membrane adapté est de simuler l'expérience des pinces optiques, qui consiste en une configuration expérimentale bien contrôlée qui permet d'étudier la mécanique individuelle d'un globule rouge isolé dans une large gamme de déformations. Plusieurs modèles pertinents sont identifiés, mais les caractéristiques de déformation mesurées durant l'expérience des pinces optiques se révèlent n'être pas assez sélectives pour être utilisées dans un contexte de validation. Des mesures de déformation additionnelles sont proposées, qui pourraient permettre une meilleure caractérisation de la mécanique de la membrane des globules rouges. En ce qui concerne les configurations macroscopiques, une méthode numérique innovante est proposée afin de gérer des simulations numériques de membranes 3D continues, en conservant le formalisme de frontières immergées. Dans cette méthode, appelée méthode des frontières immergées épaisses, la membrane a une épaisseur finie. La précision et la robustesse de la méthode sont démontrées par l'intermédiaire d'une variété de cas tests bien choisis. La méthode proposée est ensuite appliquée à un problème d'interaction fluide-structure réaliste, à savoir l'interaction d'un écoulement (sanguin) pulsé avec une valve aortique biomimétique. Une étude combinée expérimentale et numérique est menée, montrant que la méthode est capable de capturer la dynamique globale de la valve, ainsi que les principales caractéristiques de l'écoulement en aval de la valve. Tous les développements ont été effectués dans le solveur YALES2BIO (http://www.math.univ-montp2.fr/~yales2bio/) développé à l'IMAG, qui est donc disponible pour toutes autres améliorations, validations et études applicatives. / This thesis deals with several scientific aspects inherent to the numerical simulation of fluid-structure interaction problems involving thin deformable membranes. Two specific cases relevant to cardiovascular biomechanics are considered: the interaction of the blood flow with the aortic valve (which occurs at the macroscopic scale), and the interaction of the red blood cells membrane with its inner and outer fluids (which occurs at the microscopic scale). In both cases, the fluid-structure interaction coupling is handled using an immersed boundary formalism, representing the membrane by a Lagrangian mesh moving through an Eulerian fluid mesh.When dealing with red blood cells dynamics, the membrane is considered to be an infinitely thin and massless structure. The first question which is addressed in the present thesis work is how to model the complex microstructure of the red blood cells membrane. A possible way to characterize a suitable membrane model is to simulate the optical tweezers experiment, which is a well-controlled experimental configuration enabling to study the individual mechanics of an isolated red blood cell in a large range of deformation. Some relevant membrane models are identified, but the deformation characteristics measured during the optical tweezers experiment reveal to be not selective enough to be used in a validation context. Additional deformation measurements are proposed, which could allow a better characterization of the red blood cell membrane mechanics.Regarding the macroscopic configurations, an innovative numerical method is proposed to handle numerical simulations of 3D continuum membranes, still within the immersed boundary formalism. In this method, called immersed thick boundary method, the membrane has a finite thickness. The accuracy and robustness of the method are demonstrated through a variety of well-chosen test cases. Then, the proposed method is applied to a realistic fluid-structure interaction problem, namely the interaction of a pulsatile (blood) flow with a biomimetic aortic valve. A combined experimental and numerical study is led, showing that the method is able to capture the global dynamics of the valve, as well as the main features of the flow downstream of the valve.All the developments were performed within the YALES2BIO solver (http://www.math.univ-montp2.fr/~yales2bio/) developed at IMAG, which is thus available for further improvements, validations and applicative studies.

Cette thèse porte sur l’analyse mathématique de l’interaction d’un fluide non-visqueux avec des structures immergées. Plus précisément, elle est structurée autour de deux axes principaux. L’un d’eux est l’analyse asymptotique du mouvement d’une particule infinitésimale en milieu liquide. L’autre concerne l’interaction entre des vagues et une structure immergée. La première partie de la thèse repose sur l’analyse mathématique d’un système d’équations différentielles ordinaires non-linéaires d’ordre 2 modélisant le mouvement d’un solide infiniment petit dans un fluide incompressible en 2D. Les inconnues du modèle décrivent la position du solide, c’est-à-dire la position du centre de masse et son angle de rotation. Les équations proviennent de la deuxième loi de Newton avec un prototype de force de type Kutta-Joukowski. Plus précisément, nous étudions la dynamique de ce système lorsque l’inertie du solide tend vers 0. Les principaux outils utilisés sont des développements asymptotiques multiéchelles en temps. Pour la dynamique de la position du centre de masse, l’étude met en évidence des analogies avec le mouvement d’une particule chargée dans un champ électromagnétique et la théorie du centre-guide. En l’occurrence, le mouvement du centreguide est donné par une équation de point-vortex. La dynamique de l’angle est quant à elle donnée par une équation de pendule non-linéaire lentement modulée. Des régimes très différents se distinguent selon les données initiales. Pour de petites vitesses angulaires initiales la méthode de Poincaré-Lindstedt fait apparaitre une modulation des oscillations rapides, alors que pour de grandes vitesses angulaires initiales, un movement giratoire bien plus irrégulier est observé. C’est une conséquence particulière et assez spectaculaire de l’enchevêtrement des trajectoires homocliniques. La deuxième partie de la thèse porte sur le problème des vagues dans le cas où le domaine occupé par le fluide est à surface libre et avec un fond plat sur lequel un objet solide se translate horizontalement sous l’effet des forces de pression du fluide. Nous avons étudié deux systèmes asymptotiques qui décrivent le cas d’un fluide parfait incompressible en faible profondeur. Ceux-ci correspondent respectivement aux équations de Saint-Venant et de Boussinesq. Grâce à leur caractère bien-posé en temps long, les modèles traités permettent de prendre en compte certains effets de la mécanique du solide, comme les forces de friction, ainsi que les effets non-hydrostatiques. Notre analyse théorique a été complétée par des études numériques. Nous avons développé un schéma de différences finies d’ordre élevé et nous l’avons adapté à ce problème couplé afin de mettre en évidence les effets d’un solide (dont le mouvement est limité à des translations sur le fond) sur les vagues qui passent au dessus de lui. A la suite de ces travaux, nous avons souligné l’influence des forces de friction sur ce genre de systèmes couplés ainsi que sur le déferlement des vagues. Quant à l’amortissement dû aux effets hydrodynamiques, une vague ressemblance avec le phénomène de l’eau morte est mise en évidence. / This PhD thesis concerns the mathematical analysis of the interaction of an inviscid fluid with immersed structures. More precisely it revolves around two main problems: one of them is the asymptotic analysis of an infinitesimal immersed particle, the other one being the interaction of water waves with a submerged solid object. Concerning the first problem, we studied a system of second order non-linear ODEs, serving as a toy model for the motion of a rigid body immersed in a two-dimensional perfect fluid. The unknowns of the model describe the position of the object, that is the position of its center of mass and the angle of rotation; the equations arise from Newton’s second law with the consideration of a Kutta-Joukowski type lift force. It concerns the detailed analysis of the dynamic of this system when the solid inertia tends to 0. For the evolution of the position of the solid’s center of mass, the study highlights similarities with the motion of a charged particle in an electromagnetic field and the wellknown “guiding center approximation”; it turns out that the motion of the corresponding guiding center is given by a point-vortex equation. As for the angular equation, its evolution is given by a slowly-in-time modulated non-linear pendulum equation. Based on the initial values of the system one can distinguish qualitatively different regimes: for small angular velocities, by the Poincaré-Lindstedt method one observes a modulation in the fast time-scale oscillatory terms, for larger angular velocities however erratic rotational motion is observed, a consequence of Melnikov’s observations on the presence of a homoclinic tangle. About the other problem, the Cauchy problem for the water waves equations is considered in a fluid domain which has a free surface on the upper vertical limit and a flat bottom on which a solid object moves horizontally, its motion determined by the pressure forces exerted by the fluid. Two shallow water asymptotic regimes are detailed, well-posedness results are obtained for both the Saint-Venant and the Boussinesq system coupled with Newton’s equation characterizing the solid motion. Using the particular structure of the coupling terms one is able to go beyond the standard scale for the existence time of solutions to the Boussinesq system with a moving bottom. An extended numerical study has also been carried out for the latter system. A high order finite difference scheme is developed, extending the convergence ratio of previous, staggered grid based models. The discretized solid mechanics are adapted to represent important features of the original model, such as the dissipation due to the friction term. We observed qualitative differences for the transformation of a passing wave over a moving solid object as compared to an immobile one. The movement of the solid not only influences wave attenuation but it affects the shoaling process as well as the wave breaking. The importance of the coefficient of friction is also highlighted, influencing qualitative and quantitative properties of the coupled system. Furthermore, we showed the hydrodynamic damping effects of the waves on the solid motion, reminiscent of the so-called dead water phenomenon.

La réduction des émissions polluantes et l'amélioration des performances des turboréacteurs nécessitent une connaissance détaillée des phénomènes physiques mis en jeu dans une chambre de combustion. L'atomisation du carburant résulte du cisaillement engendré par un fort écoulement d'air généré dans l'injecteur. La simulation numérique directe d'écoulements avec interface permet de simuler l'ensemble du processus d'atomisation. L'utilisation de maillages Cartésiens permet la réalisation de calculs HPC efficaces et précis. Mais, une des complexités de l'atomisation vient d'une interaction forte entre le comportement de la nappe liquide et l'écoulement gazeux dans les conduites de l'injecteur, rendant impératif la simulation de l'injecteur complet. Ceci étant impossible avec des maillages Cartésiens structurés, l'objectif de cette thèse est de développer une méthode d'interface immergée permettant l'inclusion d'objets solides dans un domaine de calcul, indépendamment du maillage, afin de réaliser des DNS du système d'injection complet. Les équations de Navier-Stokes incompressibles diphasiques sont résolues à l'aide d'un algorithme de projection, l'interface liquide-gaz étant transportée avec une méthode CLSVOF conservative en masse et quantité de mouvement. La présence du solide est prise en compte grâce à la méthode d'interface immergée. Cette méthode a été appliquée à la simulation numérique de nappes liquides cisaillées pour une configuration d'injecteur utilisée en essais à l'ONERA et a permis une meilleure prédiction de la fréquence de battement de la nappe. / The reduction of polluting emissions and improvement of aeronautical engines efficiency depends on the detailed knowledge of the physical phenomena encountered in a combustion chamber. Fuel atomization results from the shearing effect induced by the high velocity airflow generated inside the injector. The Direct Numerical Simulation of interfacial flows allows the simulation of the whole atomization process, while Cartesian structured meshes allows efficient and accurate HPC computations. However, the complexity of atomization comes from a strong interaction between the jet behavior and the injector internal flow, which makes essential to simulate the whole injector system. As that is impossible on Cartesian structured grids, the main objective of this thesis is to develop an Immersed Interface Method (IIM) allowing the inclusion of solid objects in the computational domain, independently of the mesh. The incompressible two-phases Navier-Stokes equations are solved using a projection algorithm with the CLSVOF method, conservative in mass and momentum. The solid presence is taken into account using the IIM. The proposed IIM has been applied to the numerical simulation of sheared liquid sheets corresponding to an ONERA experimental configuration and allows a better prediction of the flapping frequencies of the liquid sheet.

L’objectif de cette thèse est de répondre au besoin d’accélérer la restitution des résultats de calcul d’un code CFD pour la simulation d’écoulements hydrodynamiques quasi-incompressibles. Ce code présente l’originalité de résoudre explicitement les équations de Navier-Stokes sous l’hypothèse de faible compressibilité avec des schémas numériques d’ordre élevé. Les développements effectués visent à réduire les temps de calcul à précision équivalente.Une première partie est consacrée à l’implémentation d’une formulation purement incompressible avec une résolution implicite de la pression par un schéma de projection. La formulation incompressible autorise des pas de temps plus grand en s’affranchissant de la vitesse du son, mais au prix d’une algorithmique plus complexe et de la nécessité de résoudre un système linéaire. La comparaison des deux formulations,faiblement-compressible et incompressible, tend à montrer la pertinence du schéma de projection pour les écoulements laminaires instationnaires.Un deuxième axe de développement a consisté en la proposition d’une amélioration de la méthode de frontière immergée initialement présente dans le code.Si les résultats obtenus ne sont pas encore pleinement satisfaisants, ils montrent que la montée en ordre d’une méthode de frontière immergée peut être moins contraignante en formulation incompressible.Enfin la dernière partie traite de l’immersion rapide et robuste de géométries complexes telles qu’elles peuvent être rencontrées dans l’industrie. La localisation géométrique par arbre octal permet d’évaluer rapidement une fonction de distance signée indispensable pour la méthode de frontière immergée. / An original strategy to address hydrodynamic flow was recently proposed through a high-order weakly-compressible Cartesian grid approach. The method is based on a fully-explicit temporal scheme for solving the Navier-Stokes equations. The present thesis aims to reduce the computational time required to obtain the results without deteriorating the accuracy.A first part is dedicated to the implementation of a truly incompressible formulation with an implicit solution for the pressure field through a projection scheme. The incompressible solver allows larger time step size for time integration since the speed of sound tends to infinity. In return the algorithms are no longer straight forward and a linear system has to be solved through the Pressure Poisson Equation. Comparisons carried out between both formulations show that the projection scheme can be better adapted to efficiently simulate unsteady viscous flows. Then an improvement of the immersed boundary method has been proposed. Results are not fully satisfactory for now. However, it seems easier to develop a numerical scheme for the incompressible approach rather than the weakly-compressible one.Finally, the last part addresses the setup up of complex triangulations in immersed boundary simulations. A fast and robust procedure is developed for distance computation with an octree data structure.

Le sujet de la thèse est le développement d'un outil numérique qui permet de modéliser l'écoulement autour des pales d'éoliennes. Nous nous sommes intéressés à la résolution des équations de Navier-Stokes en incompressible sur des maillages de type octree où les échelles plus petites en proche parois ont été modélisées par la méthode dite des wall functions. Un procédé d'adaptation automatique du maillage (AMR) a été développé pour affiner le maillage dans les zones où la vorticité est plus importante. Le modèle de structure d'une pale d'éolienne a été également implémenté et couplé avec le modèle fluide car une application de l'outil numérique est l'étude des effets des rafales de vent sur les pales d'éolienne. Un travail expérimental a été mené sur une éolienne avec une mesure de vent en amont. Ces données permettent ainsi de calibrer et valider les modèles numériques développés dans la thèse. / The subject of the thesis is the development of a numerical tool that allows to model the flow around wind blades. We are interested in the solving of incompressible Navier-Stokes equations on octree grids, where the smallest scales close to the wall have been modelled by the use of the so-called Wall Functions. An automatic Adaptive Mesh Refinement (AMR) process has been developed in order to refine the mesh in the areas where the vorticity is higher. The structural model of a real wind blade has also been implemented and coupled with the fluid model. Indeed, an application of the numerical tool is the study of the effects of wind gusts on blades. An experimental work has been conducted with an in-service wind turbine with the measurement of wind speed upstream. This data will allow to calibrate and validate the numerical models developed in the thesis.

La simulation d'écoulements dans des géométries complexes nécessite la création de maillages parfois difficile à réaliser. La méthode de pénalisation proposée dans ce travail permet de simplifier cette étape. En effet, la résolution des équations qui gouvernent l'écoulement se fait sur un maillage plus simple mais non-adapté à la géométrie du problème. Les conditions aux limites sur les parties du domaine physique immergées dans le maillage sont prises en compte à travers l'ajout d'un terme de pénalisation dans les équations. Nous nous sommes intéressés à l'approximation du terme de pénalisation pour une discrétisation par volumes finis sur maillages décalés et colocatifs. Les cas tests de vérification réalisés attestent d'un ordre de convergence spatial égal à 2 pour la méthode de pénalisation appliquée à la résolution d'une équation de type Poisson ou des équations de Navier-Stokes. Enfin, on présente les résultats obtenus pour la simulation d'écoulements turbulents autour d'un cylindre à Re=3900 et à l'intérieur d'une partie d'un assemblage combustible à Re=9500. / The simulations of fluid flows in complex geometries require the generation of body-fitted meshes which are difficult to create.The penalty method developed in this work is useful to simplify the mesh generation task.The governing equations of fluid flow are discretized using a finite volume method on an unfitted mesh.The immersed boundary conditions are taken into account through a penalty term added to the governing equations.We are interested in the approximation of the penalty term using a finite volume discretization with collocated and staggered grid.The penalty method is second-order spatial accurate for Poisson and Navier-Stokes equations.Finally, simulations of turbulent flows around a cylinder at Re=3900 and turbulent motions in a rod bundle at Re=9500 are performed.

Cette thèse s’inscrit dans le contexte de la simulation et du contrôle des vibrations de structures montées sur ressort qui peuvent apparaître sous l’effet de l’interaction avec l’écoulement de sillage instationnaire. Le contrôle de ce phénomène, appelé vibrations induites par vortex (VIV), est un enjeu critique dans l’optimisation de nombreux systèmes. Une méthode de frontières immergées (IBM) a été intégrée dans l’algorithme PISO du code OpenFOAM, dédié à la simulation d’écoulements fluides incompressibles. La méthode IBM permet une représentation précise de corps fixes ou en mouvement, tout en conservant des maillages structurés conduisant à des algorithmes plus précis et efficaces en termes de performances numériques. Pour calculer la divergence de l’équation de quantité de mouvement dans une boucle PISO et l’interpolation des flux, un calcul hybride orignal a été proposé avec une résolution analytique utilisant l’équation de la fonction noyau des quantités impliquant le terme force de l’IBM (quantités singulières). La méthode a été étendu au formalisme d’écoulements en régimes turbulents. Une loi de paroi a été intégrée permettant de modéliser la couche limite à grand nombre de Reynolds. Le travail de validation a été réalisé au regard des données expérimentales et numériques disponibles dans la littérature pour l’étude d’écoulements autour de cylindres et de sphères, sur une large gamme de nombres de Reynolds. Avec l’objectif de développer des lois de contrôle optimal pour le VIV, basées sur les mécanismes d’instabilité linéaire du système couplé dans le cadre de la théorie du contrôle, un solveur adjoint a été développé et validé. / This thesis is related to the simulation and the control of the vortex induced vibrations phenomenon (VIV), which can result from the fluid structure interactions between an unsteady wake and the body, when the shedding frequency in the wake is close to the natural frequency of the body. The control of VIV is a critical issue when optimizing many systems. An Immersed Boundaries Method (IBM) was implemented into the PISO algorithm as a new library of OpenFOAM, in order to perform reliable simulations of incompressible flows around bluff bodies.To compute the divergence of the momentum equation and the interpolation of the fluxes, an hybrid calculation with an analytical resolution of the quantities involving the force term (singular quantities) has been proposed. The mesh convergence of several errors was shown by means of a manufactured solution, allowing to analyze both the errors irelated to the discretization and to the IBM. The new algorithm was subsequently extended to the RANS and DDES formalism proposed in OpenFOAM for the simulation of turbulent flows. A wall law was integrated into theIBM method to model the boundary layers that develop around the bodies at large Reynolds numbers. Various 2D and 3D well-documented test cases of academic flows around fixed or moving solid bodies (cylinderand sphere) have been simulated and carefully validated against existing data from the literature in a large range of Reynolds numbers. With the objective of developing optimal control laws for VIV, based on the linear instability mechanisms of the coupled system within the framework of the control theory, a new adjoint solver was also developed and validated in OpenFOAM.

We study the turbulent flow of a fluid over a canopy, that we model as a deformable porous medium. This porous medium is more precisely a carpet of fibres that bend under the hydrodynamic load, hence initiating a fluid-structure coupling at the scale of a fibre's height (honami). The objective of the thesis is to develop a macroscopic model of this fluid-structure interaction in order to perform numerical simulations of this process. The volume averaging method is implemented to describe the large scales of the flow and their interaction with the deformable porous medium. An hybrid approach is followed due to the non-local nature of the solid phase; While the large scales of the flow are described within an Eulerian frame by applying the method of volume averaging, a Lagrangian approach is proposed to describe the ensemble of fibres. The interface between the free-flow and the porous medium is handle with a One-Domain- Approach, which we justify with the theoretical development of a mass- and momentum- balance at the fluid/porous interface. This hybrid model is then implemented in a parallel code written in C$++$, based on a fluid- solver available from the \openfoam CFD toolbox. Some preliminary results show the ability of this approach to simulate a honami within a reasonable computational cost. Prior to implementing a macroscopic model, insight into the small-scale is required. Two specific aspects of the small-scale are therefore studied in details; The first development deals with the inertial deviation from Darcy's law. A geometrical parameter is proposed to describe the effect of inertia on Darcy's law, depending on the shape of the microstructure of the porous medium. This topological parameter is shown to efficiently characterize inertia effects on a diversity of tested microstructures. An asymptotic filtration law is then derived from the closure problem arising from the volume averaging method, proposing a new framework to understand the relationship between the effect of inertia on the macroscopic fluid-solid force and the topology of the microstructure of the porous medium. A second research axis is then investigated. As we deal with a deformable porous medium, we study the effect of the pore-scale fluid-structure interaction on the filtration law as the flow within the pores is unsteady, inducing time-dependent fluidstresses on the solid- phase. For that purpose, we implement pore-scale numerical simulations of unsteady flows within deformable pores, focusing for this preliminary study on a model porous medium. Owing to the large displacements of the solid phase, an immersed boundary approach is implemented. Two different numerical methods are compared to apply the no-slip condition at the fluid-solid interface: a diffuse interface approach and a sharp interface approach. The objective is to find the proper method to afford acceptable computational time and a good reliability of the results. The comparison allows a cross-validation of the numerical results, as the two methods compare well for our cases. This numerical campaign shows that the pore-scale deformation has a significant impact on the pressure drop at the macroscopic scale. Some fundamental issues are then discussed, such as the size of a representative computational domain or the form of macroscopic equations to describe the momentum transport within a soft deformable porous medium.

Les énergies fossiles sont largement utilisées depuis les années 1900 pour satisfaire l’augmentation mondiale de la demande d’énergie. Cependant, la combustion est un procédé qui libère des polluants comme le CO2 et les NOx. Un des principaux challenges du 21ème siècle est de réduire ces émissions et les constructeurs automobiles sont impliqués dans cette course. Pour augmenter le rendement des moteurs à pistons, des solutions techniques, tels que le "downsizing", sont développées. Cette technique consiste à réduire la cylindrée des moteurs tout en maintenant leurs performances grâce à un turbocompresseur qui permet d’augmenter la masse enfermée dans la chambre de combustion. Malheureusement, l’augmentation de la pression dans les cylindres induite par le turbocompresseur est à l’origine de combustions anormales : des variations cycles à cycles importantes apparaissent, les gaz frais peuvent s’auto-allumer (allumage avant le claquage de la bougie) entrainant des phénomènes de cliquetis ou de rumble. La Simulation aux Grandes Echelles (SGE) a déjà prouvé qu’elle était un outil fiable pour prédire ces combustions anormales. Cependant ces calculs reposent sur des modèles pour prédire la propagation de la flamme dans la chambre de combustion. Ces modèles sont généralement issus de corrélations réalisées dans des cas où la turbulence est supposée homogène et isotrope. Définir théoriquement ou numériquement une telle turbulence est relativement simple mais expérimentalement la tâche est plus délicate. Cette thèse s’intéresse à un dispositif classiquement utilisé: une enceinte fermée dans laquelle la turbulence est générée par des ventilateurs. L’objectif de ce travail est donc double: 1. caractériser la turbulence générée dans ce type d’enceinte pour vérifier si elle est homogène et isotrope. 2. caractériser finement la combustion, laminaire et turbulente, afin d’enrichir les connaissances dans ce domaine et ainsi améliorer les modèles utilisés. Une première étude sur la propagation des flammes laminaires a été menée. Elle présente les effets de l’étirement et du confinement sur la vitesse de flamme laminaire. La principale difficulté pour la simulation de l’enceinte complète consiste à trouver une méthode numérique permettant de reproduire précisément l’écoulement généré par un ventilateur mais aussi d’en gérer plusieurs simultanément. Deux méthodes ont alors été testées. Premièrement, une méthode type Frontières Immergées a été implémentée dans le code de calcul AVBP. Malgré les bons résultats obtenus sur des cas tests simples, cette méthode ne s’est pas montrée adaptée pour reproduire précisément l’écoulement généré par un seul ventilateur. Une autre approche, provenant du monde du calcul des turbomachines, et basée sur le couplage de codes (appelée MISCOG), a quant à elle démontré ses capacités à le faire et est donc utilisée pour calculer l’écoulement généré par les six ventilateurs à l’intérieur de l’enceinte. L’écoulement non réactif est d’abord analysé: les résultats montrent qu’il existe une zone d’environ 6 cm de diamètre au centre de l’enceinte dans laquelle la vitesse moyenne de l’écoulement est proche de zéro et dans laquelle la turbulence est quasiment homogène et isotrope. Enfin, le pré-mélange de gaz frais est allumé en déposant un noyau de gaz chauds au centre de l’enceinte et la phase de propagation turbulente est analysée. En particulier, il est montré que la température des gaz brulés déposés au moment de l’allumage est un paramètre critique. / Fossil energy is widely used since the 1900s to satisfy the global increasing energy demand. However, combustion is a process releasing pollutants such as CO2 and NOx. One of the major challenges of the 21th century is to reduce these emissions and car manufacturers are involved in this race. To increase fuel efficiency of piston engines, some technical solutions are developed such as ‘downsizing’. It consists in reducing the engine size while maintaining its performances using a turbocharger to increase the trapped mass in the combustion chamber. Unfortunately, downsizing can lead to abnormal combustions: intense cycle to cycle variations can appear, the fresh mixture can auto-ignite (ignition before spark-plug ignition) leading to knock or rumble. Large Eddy Simulation has proven to be a reliable tool to predict these abnormal combustions in real engines. However, such computations are performed using models to predict the flame propagation in the combustion chamber. Theses models are generally based on correlations derived in cases where turbulence is assumed to be homogeneous and isotropic. Defining theoretically or numerically such a turbulence is a simple task but experimentally it is more challenging. This thesis focuses on a apparatus used in most experimental systems: fans stirred vessel. The objective of this work is twofold: 1. characterize the turbulence generated inside the vessel to check wether it is homogeneous and isotropic or not, 2. finely characterize laminar and turbulent combustion in this setup in order to increase the knowledge in this field, and thereby improve models used. First, a laminar flame propagation study has been conducted to address both confinement and curvature effects on the laminar flame speed in a spherical configuration. The main difficulty to perform the simulation of the whole configuration consists in finding a numerical method able to compute accurately the flow generated by one fan and able to handle six fans simultaneously too. Two numerical methodologies have been tested. First an Immersed Boundaries method was implemented. Despite good results obtained on academic test cases, this method was shown to be unadapted to compute accurately the flow generated by one fan. On the other hand, a numerical approach, coming from turbomachinery calculations and based on code coupling (called MISCOG), demonstrates its ability to do it and it is used to compute the flow generated by the six fans inside the closed vessel. Non-reacting flow is first analyzed and reveals a zone at the vessel center of around 6 cm of diameter where mean velocity is near zero and turbulence is almost homogeneous and isotropic. After that, the premixed fresh mixture is ignited depositing a hot gases kernel at the vessel center and the turbulent propagation phase is analyzed. In particular, it is shown that the amount of energy deposited at ignition is a critical parameter.

Cette thèse porte sur l'étude des ballotements dans les réservoirs de satellites à poste, lors des phases de manoeuvre à faible accélération. En effet la bulle de gaz d'hélium servant à pressuriser le réservoir se met en mouvement, générant ainsi des perturbations sur la stabilité globale du satellite. Afin de mener à bien cette étude, des méthodes numériques ont été développées, avec une méthode de frontières immergées pour prendre en compte les parois du réservoir.Le code est utilise la méthode Level-Set pour capturer l'interface, et gère les sauts à l'aide de la méthode Ghost-Fluid. Un solveur BlackBox Multigrid est également développé pour améliorer lesperformances de calcul. Une étude est présentée dans le dernier chapitre pour définir quelques lois de comportements en fonction des vitesses et accélérations générées lors des manoeuvres. / The core study of this PhD thesis is the sloshing in satellite tanks, during low acceleration maneuvers. Indeed the helium bubble used to pressurize the tank moves, thus generating perturbations on the global stability of the satellite. In order to understand this problem, numerical schemes have been developed, such as an immersed boundary method to model the tank wall. The numerical tool uses a Level-Set function coupled to a Ghost Fluid Method to track the interface and to account for the jump conditions.A BlackBox Multigrid Solver have been developed to improve computational cost. Finally a study is presented in the last chapter to predict the behaviour of the fluids with a varying rotational speed generated during some classical maneuvers.