Une méthode de frontières immergées est développée pour la simulation d'écoulements compressibles et validée au travers de cas-tests spécifiques (réflexion d'ondes acoustiques et quantification de la conservation de la masse dans des canaux inclinés). La simulation aux grandes échelles (LES) d'une cavité transsonique est ensuite présentée. Le bouclage aéro-acoustique, très sensible aux conditions aux limites, est reproduit avec précision par la LES dans le cas où les parois sont immergées dans un maillage structurée. La comparaison des stratégies de modélisation de sous-maille pour cet écoulement transsonique et l'adaptation des filtres en présence de frontières immergées sont également discutées. Le rôle, souvent sous-estimé, du schéma de viscosité artificiel, est quantifié.Dans la dernière partie du manuscrit, des études sont réalisées pour aider au dimensionnement d'un nouveau concept de chambre de combustion où la flamme est stabilisée par la recirculation de gaz brûlés dans une cavité (chambre TVC pour Trapped Vortex Combustor). La modélisation de la combustion turbulente est basée sur une chimie tabulée, couplée à une fonction densité de probabilité présumée (PCM-FPI). L'étude de la dynamique de la flamme est réalisée pour diverses conditions de fonctionnement (débit de l'écoulement principal et présence ou non d'un swirl). Les spécificités de mise en œuvre de la simulation d'un écoulement de ce type sont discutées et un soin particulier est apporté au traitement de la condition de sortie, qui constitue un point sensible de la chaîne de modélisation. Les phénomènes d'instabilités et de retour de la flamme sont mis en évidence ainsi que les modifications à apporter au dispositif afin de minimiser ces effets. L'existence d'un cycle limite acoustique est souligné et une formule permettant d'anticiper le niveau des fluctuations de pression est proposée et validée. Une correction au modèle PCM-FPI est présentée afin de préserver la vitesse de flamme et d'assurer une reproduction plus précise de la dynamique de flamme.

Les travaux menés, depuis plusieurs années, au CORIA ont abouti à la construction d'un outil numérique (ARCHER) permettant la simulation numérique directe d'écoulements diphasiques et notamment l'atomisation d'un jet liquide à haute vitesse. Ce type de simulation permet de capturer les phénomènes d'atomisation au voisinage de l'injecteur difficilement caractérisables par les outils expérimentaux actuels. Ces simulations requièrent des conditions d'injection délicates à évaluer a priori car elles dépendent des caractéristiques de l'écoulement au sein de l'injecteur. Or, certains jets présentent une grande sensibilité à ces conditions d'injection. Dès lors, il est nécessaire de simuler l'écoulement au sein de l'injecteur afin d'appréhender la nature de cette sensibilité. L'utilisation d'un maillage cartésien par le code ARCHER conjuguée à la volonté de simuler le système d'atomisation dans son ensemble ont orienté ces travaux vers l'utilisation d'une méthode de frontière immergée. Ces travaux ont ainsi permis de reproduire des écoulements au sein d'injecteurs de forme quelconque tout en conservant le maillage cartésien d'origine, précieux tant pour l'efficacité du solveur que pour sa précision. Dans un premier temps, l'implantation dans le code ARCHER d'une méthode de frontière immergée a été réalisée et testée sur des configurations de canal et de conduite et de l'écoulement autour d'un cylindre. L'application de cette méthode a porté sur la simulation de l'écoulement au sein d'un injecteur triple disque mono-trou et a notamment permis de caractériser l'origine de l'écoulement secondaire formé dans l'orifice de décharge. Afin d'évoluer vers la construction d'un outil numérique capable de simuler le système d'atomisation dans son ensemble, un couplage entre la méthode de frontière immergée et la méthode Ghost fluid a été nécessaire. La version bi-dimensionnelle développée a été testée sur la relaxation d'une goutte posée sur une paroi. Cette version a permis de simuler des écoulements au sein de canaux à différents rapports de longueur sur diamètre et l'écoulement au sein d'une buse convergente. La simulation simultanée de l'écoulement interne et externe a permis de lier les fluctuations de vitesses des écoulements internes à la création de surface engendrée sur les écoulements externes.

Une capsule est une goutte de liquide enveloppée par une membrane fine et déformable. Les propriétés mécaniques de la membrane sont essentielles pour le mouvement de la capsule. L'analyse de l'écoulement d'une suspension de capsules dans un canal microfluidique au moyen d'un modèle mécanique est une technique permettant de déterminer les propriétés élastiques de la membrane. Un modèle numérique tridimensionnel a été développé pour résoudre ce problème d'interaction fluide-structure en écoulement confiné. Il couple une méthode des intégrales de frontières pour les écoulements des fluides et une méthode éléments finis pour la déformation de la membrane. Le modèle est utilisé pour étudier l'écoulement d'une capsule initialement sphérique dans des canaux de différentes sections. Dans un canal cylindrique, on montre que l'effet de confinement du canal conduit à la compression de la capsule. Cela engendre la formation de plis sur la membrane autour de l'axe de l'écoulement, phénomène également observé expérimentalement. Dans un canal de section carrée, les effets de la loi constitutive de la membrane, du rapport de taille et du débit d'écoulement sur la déformation de la capsule sont systématiquement étudiés. La comparaison entre les résultats expérimentaux et numériques nous permet de déduire les propriétés mécaniques de la membrane d'une population de capsules artificielles. Ce travail démontre la faisabilité de la mesure de propriétés mécaniques d'une membrane en utilisant une technique microfluidique en canal carré. Il pourrait être étendu par l'étude d'écoulements instationnaires dans un canal de section variable ou avec bifurcations.

Le bassin de santé de la Haute Casamance est à cheval sur quatre Etats que sont le Sénégal, la Guinée Bissau, la Guinée et la Gambie au destin lié par un passé de pratique territoriale commune. Toutefois, la construction des Etats au lendemain de la conquête coloniale a instauré une coupure de fait en balkanisant cet ensemble territorial par l'établissement de frontières artificielles en déphasage avec la manière dont les populations vivent l'espace à travers une différenciation propre à la typologie naturelle entre zones sahélienne, soudanienne et forestières. Le vécu spatial séculaire des populations toujours valable favorise un eternel recommencement de couture de cet ensemble territorial en faisant fi des frontières modernes. Ainsi des polarités transfrontalières se construisent en arborant un statut de centre gravité autour d'un espace commun entre plusieurs Etats. Le secteur de la santé n'y échappe pas. Des localités comme Kolda au coeur de cet espace transfrontalier attirent des populations d'une origine très diversifiée et non prévues au sein des infrastructures sanitaires. Des dynamiques de recours aux soins se sont mises en place et fortement caractérisées par une mobilité de soins qui façonne un territoire sanitaire original. Celui-ci entretient des flux de population à la recherche de meilleurs soins en relation avec l'accessibilité physique, la qualité des soins proposée en certains endroits, les activités économiques. Ce pendant, l'offre de soins ne suit pas cette logique de continuité du fait de l'absence de connexion des systèmes de soins des différents pays. C'est la conséquence d'une absence de coopération transfrontalière sanitaire qui fait de cet espace de partage un cadre sanitaire frappé par des disparités notables entre les Etats aux problèmes commun en matière d'endémicité.

Cette thèse s’applique à rendre compte de la manière dont des tensions s’incarnent et s’expriment dans les activités quotidiennes d’une organisation humanitaire, Médecins sans frontières. Pour ce faire, nous faisons appel à une approche dite « constitutive » de la communication, approche qui vise à montrer que les tensions organisationnelles que vivent les intervenants humanitaires dans leur travail peuvent s’analyser à partir des figures qui animent leurs conversations et leurs activités de tous les jours. Par tension, nous entendons l’expérience d’une incompatibilité, antinomie ou opposition vécue ou perçue entre deux ou plusieurs logiques d’action, tandis que la notion de figure renvoie à ce qui anime ces mêmes logiques d’action, soit toute préoccupation, valeur ou intérêt au nom desquels un acteur s’exprime ou agit. Dans cette thèse, nous avons ainsi identifié puis analysé, grâce à une approche d’inspiration ethnographique, cinq figures, trois tensions et trois « scènes communicationnelles » qui, selon nous, illustrent d’une manière fidèle les activités humanitaires de MSF. La présente thèse se veut donc une démonstration empirique, via les pratiques communicationnelles des acteurs, de ces figures et tensions à gérer au quotidien, comblant du coup certaines lacunes dans la littérature portant sur les tensions organisationnelles et l’aide humanitaire. Ces lacunes s’articulent principalement autour d’un manque d’intérêt pour, d’une part, la réalité des pratiques des acteurs en tant que telles, et, d’autre part, l’application d’une vision de la communication proprement constitutive de l’organisation, c’est-à-dire une vision qui reflèterait et exprimerait ce qui anime ses principaux acteurs au quotidien. Comme nous le montrons, ce qui semble animer ou préoccuper MSF (et ses acteurs), ce sont, d’une part, certains principes ou valeurs tels l’indépendance de pratique, la neutralité ou encore l’expérience et le scrupule. Ces idéaux, que les volontaires de MSF cultivent dans leurs conversations et leurs activités, sont donc constamment pris en considération lors de la mise sur pied et la réalisation des missions de cette célèbre organisation. D’autre part, une autre source d’animation se retrouve aussi au niveau des lieux et environnements complexes dans lesquels les soins médicaux sont prodigués aux populations en danger. Au final, et grâce à l’approche communicationnelle mise à profit dans cette thèse, il nous a été possible de dresser un portrait réaliste de ces mêmes pratiques - les « sources » de ce qui les anime -, reflétant à la fois le mode d’être et d’agir de MSF, dans toutes ses intensités et toutes ses tensions. / This thesis illustrates the ways by which some tensions enbody and expresses themselves in the daily activities of a humanitarian organization, Doctors without borders (MSF). To do so, we mobilize a "constitutive" approach to communication, an approach that allows us to show that organizational tensions faced by (humanitarian) actors can be analyzed from the figures that animate their conversations and activities on a daily basis. By tension, we mean every practices of inconsistency, contradiction or opposition experienced or perceived between two or more logics of action, while the notion of figure refers to what drives these logics of action, that is, any concern, value or interest on behalf of which an actor comes to speak or act. In this thesis, we thus identified and analyzed, through an ethnographic approach, five figures, three tensions and three "communicative scenes" that we believe illustrate quite faithfully the humanitarian activities of MSF. This thesis is therefore an empirical demonstration, via the communication practices of actors, of these figures and tensions that have to be managed on a daily basis. Consequently, we were able to fill some “gaps” in the literature on organizational tensions and humanitarian aid. These gaps correspond with a lack of interest for the actual practices of organizational actors as such, and secondly, as well as for an approach to communication conceived as being constitutive of an organization, that is to say, an approach that reflects and expresses what drives or guides organizational actors on a daily basis. As we show, what seem to animate or preoccupy MSF (and its actors) are, on one side, some principles and values such as independence of practice, neutrality or experience and scruples. Thus, these ideals that actors cultivate in their conversations and activities are constantly taken in account when missions are planned and carried out for this famous organization. On the other side, another source of animation lies in the complex environments in which medical care is provided to populations in distress. Through our constitutive approach to communication, we were thus able to draw a realistic picture of these practices - the "sources" of what animates humanitarian action - reflecting both the mode of being and acting of MSF, in all its intensities and tensions.

Les vésicules sont utilisées d'une manière extensive comme modèle pour comprendre les dynamiques et les déformations des globules rouges au niveau individuel, mais aussi concernant les phénomènes collectives et la rhéologie. La membrane de la vésicule résiste à la flexion mais pas au cisaillement, contrairement aux globules rouges, néanmoins elles partagent plusieurs propriétés dynamiques avec les globules rouges, comme le tank-treading (mouvement en chenille de char) et le tumbling (mouvement de bascule) sous écoulement de cisaillement, ou les formes parachutes et slippers (pantoufles) sous un écoulement de Poiseuille. Les globules rouges sont connus pour former des trains de cellules (clusters) dans la microcirculation attribués à la nature attractive des interactions hydrodynamiques. Nous avons étudié numériquement plusieurs types de problème comme:(i) les dynamiques de cellules isolées, (ii) le couplage hydrodynamique entre globules rouges (en utilisant les vésicules comme modèle) soumis à un écoulement de Poiseuille sous différent confinements; (iii) l'agrégation des globules rouges et la formation de rouleaux; et (iv) le rôle des macromolécules dans la formation de clusters sous écoulement. les résultats obtenus apportent un nouveau regard à la physique des objets déformables et sont transposables au cas de l'écoulement des globules rouges dans la microcirculation. / Vesicles are extensively used as a model for understanding dynamicsand deformation of red blood cells at the individual level but also regarding collective phenomena and rheology. Vesicles' membranes withstand to bending butdo not have a shear resistance, unlike red blood cells, but they still share several dynamical properties with red blood cells, like tank-treading and tumbling under linear shear flow, or parachute and slipper shapes under Poiseuille flow. The red blood cells are known to form train of cells in the microcirculation attributed to attractive hydrodynamic interactions. We investigate numerically several kind of problems such as: (i) the dynamics of isolated cells; (ii) the hydrodynamic coupling between the red blood cells (by using vesicles as a model) subject to a Poiseuille flow under different confinements; (iii) the aggregation of red blood cells and formation of rouleaux; and (iv) the contribution of macromolecules in the formation of clusters under flow condition. The obtained results give a new insight into thephysics of deformable objects under confinement that are transposable to the flow of red blood cells in the microcirculation.

L'objectif de cette thèse est d'améliorer la compréhension des phénomènes jouant un rôle dans la mesure effectuée dans un analyseur sanguin, en particulier le comptage et la mesure de volumétrie d'une population de globules rouges reposant sur l'effet Coulter. Des simulations numériques sont effectuées dans le but de prédire la dynamique des globules rouges dans les zones de mesure et pour reproduire la mesure électrique associée, servant au comptage et à la volumétrie des cellules. Ces simulations sont effectuées à l'intérieur de configurations industrielles d'analyseur sanguin, en utilisant un outil numérique développé à l'IMAG, le solveur YALES2BIO. En utilisant la méthode des frontières immergées avec suivi de front, un modèle de particule déformable est introduit, celui-ci prend en compte le contraste de viscosité ainsi que les effets mécaniques de la courbure et de l'élasticité sur la membrane. Le solveur est validé grâce à de nombreux cas tests parcourant différents régimes et effets physiques. L'écoulement fluide dans cette géométrie d'analyseur sanguin est caractérisée par un fort gradient de vitesse axial dans la direction de l'écoulement, impliquant la présence d'un écoulement extensionnel au niveau du micro-orifice, là où a lieu la mesure. La dynamique des globules rouges est étudiée par des simulations numériques pour différentes conditions initiales, telles que sa position ou son orientation. Il est observé que les globules rouges vont se réorienter selon l'axe principal de l'analyseur sanguin dans tous les cas. Pour comprendre le phénomène, des modèles analytiques sont adaptés au cas des écoulements extensionnels et reproduisent correctement les tendances de réorientation.Cette thèse présente également la reproduction de la mesure électrique utilisée pour le comptage et la mesure de la distribution des volumes de globules rouges. De nombreuses simulations de la dynamique des globules rouges sont effectuées et utilisées pour générer l'impulsion électrique correspondant au passage du globule rouge dans le micro-orifice. Les amplitudes d'impulsions électriques résultantes permettent la caractérisation de la réponse électrique en fonction des paramètres initiaux de la simulation par une approche statistique. Un algorithme de Monte-Carlo est utilisé pour la quantification des erreurs de mesure liées à l'orientation et la position des globules rouges dans le micro-orifice. Ceci permet la génération d'une distribution de volume mesurée pour une population de globules rouges bien définie et la caractérisation des erreurs de mesure associées. / The aim of this thesis is to improve the understanding of the phenomena involved in the measurement performed in a blood analyzer, namely the counting and sizing of red blood cells based on the Coulter effect. Numerical simulations are performed to predict the dynamics of red blood cells in the measurement regions, and to reproduce the associated electrical measurement used to count and size the cells. These numerical simulations are performed in industrial configurations using a numerical tool developed at IMAG, the YALES2BIO solver. Using the Front-Tracking Immersed Boundary Method, a deformable particle model for the red blood cell is introduced which takes the viscosity contrast as well as the mechanical effects of the curvature and elasticity on the membrane into account. The solver is validated against several test cases spreading over a large range of regimes and physical effects.The velocity field in the blood analyzer geometry is found to consist of an intense axial velocity gradient in the direction of the flow, resulting in a extensional flow at the micro-orifice, where the measurement is performed. The dynamics of the red blood cells is studied with numerical simulations with different initial conditions, such as its position or orientation. They are found to reorient along the main axis of the blood analyzer in all cases. In order to understand the phenomenon, analytical models are adapted to the case of extensional flows and are found to reproduce the observed trends.This thesis also presents the reproduction of the electrical measurement used to count red blood cells and measure their volume distribution. Numerous dynamics simulations are performed and used to generate the electrical pulse corresponding to the passage of a red blood cell inside the micro-orifice. The resulting electrical pulse amplitudes are used to characterize the electrical response depending on the initial parameters of the simulation by means of a statistical approach. A Monte-Carlo algorithm helps quantifying the errors on the measurement of cell depending on its orientation and position inside the micro-orifice. This allows the generation of a measured volume distribution of a well defined red blood cell population and the characterization of the associated measurement errors.

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.

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.