Ce travail de thèse concerne l'étude des écoulements de fluides gazeux compressibles laminaires subsonique-supersonique dans une tuyère de type convergent-divergent. Les écoulements étudiés sont à nombres de Reynolds modérés et s'affranchissent de l'hypothèse de condition adiabatique de paroi couramment utilisée afin de mieux prendre en compte les phénomènes de transfert de chaleur par convection et rayonnement avec le milieu extérieur. Cette étude des phénomènes de transferts conjugués a permis de déterminer le comportement dynamique simultané du fluide et de la paroi de la tuyère. Enfin, compte tenu des niveaux élevés de températures mis en jeu dans ces systèmes, une analyse concernant le comportement thermomécanique de l'ensemble de la structure de paroi avec des matériaux monocouches et multicouches a été réalisé. De nombreuses configurations géométriques, propriétés physiques et conditions aux limites sur le fluide et la paroi ont été analysées. Les résultats présentés montrent, la structure des écoulements à travers les iso-contours de vitesses, des nombres de Mach, des pressions dans le fluide, des températures dans le fluide et dans la paroi ainsi que les déformations et les contraintes de la paroi qui résultent des couplages thermomécaniques. Une analyse des performances de la tuyère, en termes de force de poussée et de coefficient de débit spécifique, est largement discutée dans ce travail.

Les ondes gravito-inertielles jouent un rôle majeur dans les échanges d'énergie globaux sur la planète. Si la génération des ondes est bien connue dans l'atmosphère et l'océan, le devenir de ces ondes au cours de leur propagation n'est pas complètement défini aujourd'hui. Ces ondes peuvent interagir de façon non-linéaire avec elles-mêmes et créer des structures de plus petite échelle qui vont se dissiper plus facilement. Ainsi, le phénomène d'instabilité paramétrique sous-harmonique (PSI), a été étudié de façon expérimentale. Nous avons effectué la première mise en évidence expérimentale de l'interaction de trois ondes planes inertielles bi-dimensionnelles, sous la forme d'une triade résonnante. Cette étude améliore en outre la compréhension de la turbulence en rotation. Les ondes internes peuvent aussi créer, ou interagir avec des écoulements lents de grande échellequi peuvent modifier la biodiversité au fond des océans. Nous avons mis en évidence une situation expérimentale à l'origine d'un tel écoulement moyen induit par les ondes et, à l'aide d'un modèle théorique simplifié, nous avons expliqué la formation de ces écoulements. Enfin, on étudie également des tourbillons en fluide stratifié pour permettre de futures études sur l'interaction d'ondes gravito-inertielles avec des tourbillons.

La dynamique des courants de gravité turbulents sur une pente uniforme en rotation est analysée par des expériences. Ces écoulements denses jouent un rôle important en océanographie dans le renouvellement des eaux profondes qui font parties du cycle de convection thermohaline global. La grande cuve tournante Coriolis (LEGI) permet d'étudier à l'échelle du laboratoire un courant de gravité en similitude avec les échelles océaniques.<br /> La propagation de ces courants est fortement influencée par la rotation et le mélange turbulent avec le milieu ambiant, et sa dynamique est instable. Les interactions entre le courant et le milieu ambiant liées aux phénomènes d'entraînement et d'extraction de fluide, contrôlent la position de stabilisation du courant le long des côtes, sa densité et sa vitesse. Des lois d'échelle basées sur la conservation du flux de flottabilité sont proposées, décrivant l'évolution des principales caractéristiques de l'écoulement le long de la pente. <br />Elles montrent le caractère auto-similaire du courant de gravité. L'analyse des résultats expérimentaux permet de vérifier ces lois et d'étudier la structure turbulente de l'écoulement.<br />Des instabilités de plus grande échelle sont observées, générant des tourbillons cycloniques qui sont également visibles dans les écoulements océaniques. Leur étude permet de suggérer le mécanisme de leur formation. Les résultats expérimentaux sont comparés par similitude aux mesures in-situ du courant dense s'écoulant par le détroit du Danemark.

Cette thèse vise à une meilleure compréhension de la dynamique du voilier et à la validation des outils numériques de prédiction de performances et d'optimisation par l'étude expérimentale in situ du problème aéro-élastique d'un gréement. Une instrumentation est développée sur un voilier de 8m de type J80 pour la mesure dynamique des efforts dans le gréement, de la forme des voiles en navigation, du vent et des attitudes du bateau. Un effort particulier est apporté à la mesure des caractéristiques géométriques et mécaniques des éléments du gréement, la calibration des capteurs et au système d'acquisition des données. Les principaux résultats montrent que le voilier instrumenté est un outil adapté pour les mesures instationnaires et soulignent l'amplitude de variation d'effort rencontrée en mer (20 à 50% de l'effort moyen dans une houle modérée). En outre, les variations du signal d'effort sont déphasées avec l'angle d'assiette, créant un phénomène d'hystérésis. Le comportement dynamique d'un voilier en mouvement diffère ainsi de l'approche quasi-statique. Les simulations numériques proviennent du code ARAVANTI, couplage implicite d'un code structure éléments finis ARA et d'un code fluide parfait, limitant son domaine de validité aux allures de près Les résultats de simulation sont très proches des cas stationnaires et concordent bien avec les mesures en instationnaire dans une houle de face. L'expérimentation numérique d'un gréement soumis à des oscillations harmoniques en tangage souligne l'importance de l'approche Interaction Fluide Structure (IFS) et montre que l'énergie échangée par le système avec la houle est reliée à la fréquence réduite et l'amplitude du mouvement. Certaines informations n'étant pas disponibles sur le voilier instrumenté, une expérience contrôlée en laboratoire est développée. Elle consiste en un carré de tissu tenu par deux lattes en oscillation forcée. Les mesures sur cette " voile oscillante " permettent d'étudier les phénomènes IFS avec décollement et sont utilisées pour la validation du couplage ARA-ISIS entre un code fluide Navier-Stokes (RANS) et le même code structure.

La simulation numérique des écoulements turbulents cavitants revêt de nombreuses difficultés tant dans la modélisation des phénomènes physiques que dans le développement de méthodes numériques robustes. En effet de tels écoulements sont caractérisés par un changement de phase associé à des gradients de la masse volumique, des variations du nombre de Mach causées par une chute de la vitesse du son, des zones de turbulence diphasique et la présence d'instationnarités.Les travaux de la présente thèse s'inscrivent dans la continuité des études expérimentales et numériques menées au sein du Laboratoire des Ecoulements Géophysiques et Industriels (LEGI),qui visent à améliorer la compréhension et la modélisation d'écoulements cavitants. Les simulations s'appuient sur un code compressible associé à une technique de pré-contionnement bas-Mach qui permet de traiter les zones incompressibles. Les écoulements diphasiques sont reproduits à l'aide d'un modèle de mélange homogène 1-fluide avec discrétisation implicite en pas de temps dual. Enfin la résolution adopte l'approche moyennée RANS qui couple le système des équations de conservation avec des modèles de turbulence du premier ordre basés sur la notion de viscosité turbulente.Dans les zones diphasiques, le calcul des variables thermodynamiques nécessite l'introduction d'équations d'état. La pression au sein du mélange est ainsi reliée aux grandeurs conservatives soit à partir d'une équation d'état de mélange des gaz raides, soit par une relation sinusoïdale incorporant la fraction volumique de vapeur (le taux de vide). La valeur ajoutée de ces travaux de thèse repose sur l'introduction d'une équation de transport pour le calcul du taux de vide. Celle-ci incorpore un terme source dont le transfert de masse entre les phases est fermé grâce à une hypothèse de proportionnalité à la divergence du champ de vitesse. Outre l'amélioration des phénomènes de convection, de dilatation et de collapse, cette équation supplémentaire permet de relaxer l'équilibre thermodynamique local et d'introduire un état métastable pour la phase vapeur.Les simulations 2D et 3D sont réalisées sur des géométries de type Venturi caractérisées par le développement de poches de cavitation partielle instables. L'objectif consiste à reproduire les instationnarités inhérentes à chaque profil telles que la formation d'un jet rentrant liquide à proximité de la paroi ou la production de nuages de vapeur convectés par l'écoulement principal.Les résultats numériques mettent en avant une variation de la fréquence des instationnarités en fonction du calcul de la vitesse du son en zone de mélange. D'autre part, la prise en compte de déséquilibre de la phase vapeur amplifie les phénomènes de propagation d'ondes de pression générées par le collapse des structures cavitantes et participe à la déstabilisation de la poche. Enfin, l'influence de l'équation de transport de taux de vide est analysée en confrontant les résultats des simulations à ceux obtenus ultérieurement à partir d'un modèle à seulement trois équations de conservation. / The computation of turbulent cavitating flows involves many difficulties both in modeling the physical phenomena and in the development of robust numerical methods. Indeed such flows are characterized by phase transitions and large density gradients, Mach number variation due to speed of sound decrease, two-phase turbulent areas and unsteadiness.This thesis follows experimental and numerical studies led at the Laboratoire des Ecoulements Géophysiques et Industriels which aim to improve the understanding and modeling of cavitating flows. Simulations are based on a compressible code coupled with a pre-conditionning technique which handles low-Mach number areas. The two-phase flows are reproduced using a one-fluid homogeneous model and temporal discretisation is performed using an implicit dual-time stepping method . The resolution is based on the RANS approach that couples conservation equations with firts-order closure models to compute eddy viscosity.In two-phase flows areas, the computation of thermodynamic quantities requires to close the system with equations of state (EOS). Thus, two formulations are investigated to determine the pressure in the mixture. The stiffened gas EOS is written with conservative quantities while a sinusoidal law deduces the pressure from the volume fraction of vapor (the void fraction). The present study improves the homogeneous equilibrium models by including a transport equation for the void ratio. The mass transfer between phases is assumed to be proportional to the divergence of the velocity. In addition to a better modeling of convection, expansion and collapse phenomenon, this added transport equation allows to relax the local thermodynamic equilibrium and to introduce a mestastable state to the vapor phase.2D and 3D simulations are performed on Venturi type geometries characterized by the development of unstable partial cavitation pockets. The goal is to reproduce unsteadiness linked to each profile such as the formation of a re-entrant jet or the quasi-periodic vapor clouds shedding. Numerical results highlight frequency variations of unsteadiness depending on the speed of sound computation. Moreover, the simulation conducted with a relaxed vapor density increase the pressure wave propagation magnitude generated by the collapse of cavitating structures. It contributes to the destabilization of the pocket. Finally, the role of the void ratio equation is analyzed by comparing the simulation results to those obtained subsequently from a model involving only three conservation equations.

Pour être utilisé en réfrigération secondaire comme matériau à changement de phase dans un fluide frigoporteur, un hydrate doit posséder une enthalpie de dissociation élevée, à une température adaptée à l’application (positive pour de la climatisation) et à une pression acceptable pour une installation industrielle. De plus, le mélange diphasique doit présenter de bonnes conditions d’écoulement et d’échanges thermiques. Un dispositif dynamique constitué d'un réacteur agité et d'une boucle de circulation, placé dans une chambre climatique, a été développé pour l’étude de coulis d’hydrates de dioxyde de carbone. Pour ce dispositif, un protocole a été mis en œuvre pour le refroidissement d’une solution chargée en CO2, générant des coulis d'hydrate à taux d’hydrate allant jusqu'à 25°% véhiculant une chaleur latente importante (374 kJ.kg-1). Une étude cinétique sur la formation des hydrates de CO2 dans le réacteur agité à été réalisé et montre que la cinétique de formation de ces hydrates est de 2 %.h-1 dans le réacteur pour une pression, une température, et une vitesse d'agitation de départ de 3 MPa, 12 °C et 1042 tr/min et pour une température de consigne du groupe froid de -1 °C. L’étude rhéologique du coulis a permis de mettre en place une corrélation empirique de type Herschel-Bulkley en fonction du taux d’hydrates qui permet une caractérisation satisfaisante du comportement rhéologique du fluide. Les résultats de l’étude thermique ont permis d’établir des corrélations entre les nombres de Nusselt et de Reynolds qui caractérisent les régimes de fonctionnement thermique du système. / To be used as phase change material in a secondary refrigerant, a hydrate must enclose a high dissociation enthalpy, at a temperature adapted to the application (positive for air-conditioning) and a pressure suitable for an industrial facility. Moreover, the two-phase mixture must present good conditions of flowing and heat exchange. A dynamical set-up with a stirred tank reactor with a loop of circulation, placed in a thermostated room, was developed for the study of carbon-dioxide-hydrate slurries. In this set-up, one protocol was carried out: by cooling of a water-CO2 solution, generating slurries hydrate rates 25 %, conveying a high latent heat (374 kJ.kg-1). A study of CO2 hydrate formation kinetics was lead in the tank reactor and the results showed the CO2 hydrate formation kinetics is 2 %.h-1 for a pressure, temperature and a stirring rate at the begin of 3 MPa, 12°C and 1042°rpm and a temperature of the chiller about -1°C. A study rheological allows the development of an empirical Herschel-Bulkley-type correlation as a function of the hydrate rate was achieved for rheological characterization of the fluid behavior. The thermal results obtained allow establishing correlation as a function of Nusselt and Reynolds number which characterize operational regimes of the system thermal.

Dans la conception des turbomachines modernes, la prédiction des phénomènes aéroélastiques est devenue un point clé. La tendance à réduire la masse et à augmenter la charge des composants aérodynamiques accroit le risque de rupture. Dans un tel contexte, la compréhension et la bonne prédiction des diverses instabilités constituent un enjeu industriel et scientifique majeur. Le présent travail de recherche a pour objectif d’améliorer la prédiction des phénomènes instationnaires intervenant dans les problèmes d’aéroélasticité en turbomachines. Cette thèse est plus particulièrement axée sur la simulation de l’interaction onde de choc/couche limite. Le support d’étude est une tuyère transsonique présentant un écoulement avec des zones décollées. L’oscillation forcée de l’onde de choc est simulée grâce à une méthode de petites perturbations instationnaires couplée avec une hypothèse de turbulence variable. Cette approche est validée par comparaison avec des mesures. Elle permet une prédiction tout à fait satisfaisante du premier harmonique de pression sur la paroi de la tuyère. Ce travail a montré la nécessité de linéariser le modèle de turbulence. Le besoin de dériver le modèle de turbulence nous a amené à investiguer la modélisation faite pour prédire l’interaction onde de choc/couche limite. Un modèle de turbulence à deux équations complété par une équation de « retard » est implémenté afin de capter un déséquilibre de la turbulence. Les résultats obtenus en tuyère sont cohérents avec la théorie mais une surproduction d’énergie turbulente en présence de bord d’attaque rend le modèle inefficace pour des configurations de turbomachines. Au final, l’introduction d’un limiteur de viscosité turbulente dans un modèle de turbulence à deux équations s’avère donner de bons résultats. La méthode de dérivation du modèle est alors présentée sur le modèle de Wilcox proposé en 2008. Enfin, la technique de linéarisation est étendue à la problématique aéroélastique. Une approche de couplage fluide-structure faible est adoptée. L’oscillation structurelle des aubages suivant les modes propres est considérée mais en laissant la fréquence évoluer au cours du couplage. La nouvelle méthode utilisée s’appuie sur la construction d’un méta-modèle du comportement dynamique du fluide afin de résoudre directement le système fluide-structure couplé. Cette technique est validée sur une configuration de grille annulaire de turbine en haut subsonique et présente l’avantage d’un temps de calcul réduit. / In modern turbomachinery design, predicting aerolastic phenomena has become a key point. The development of highly loaded components, while reducing their weight, increases the risk of failure. In this context, good understanding and prediction of various instabilities are a major industrial and scientific challenge. This research work aims to improve the prediction of unsteady phenomena involved in turbomachinery aeroelasticity. This study focuses especially on the simulation of shock wave/boundary layer interaction. To begin with, a transonic nozzle separated flow is investigated. Forced oscillation of the shock wave system is simulated through a small unsteady perturbation method combined with the assumption of variable turbulence. This approach is validated against exprimental measurements. The first harmonic of pressure on the wall of the nozzle is predicted quite satisfactorily. The need to linearize the turbulence model was shown of high importance. Deriving the turbulence model, leads us to investigate the turbulence modeling performed to predict the shockwave/boundary layer interaction. A two equations turbulence model supplemented by a "time-lagged" equation is implemented to capture non-equilibrium effects of turbulence. All achieved results for a nozzle are consistent with theory, but overproduction of turbulent kinetic energy at leading edge makes the model useless for turbomachinery configurations. However, the introduction of an eddy viscosity stress limiter inside a two-equation turbulence model proves to give good results. The derivation method is thus presented on such a model, precisely on Wilcox model proposed in 2008. Finally, the linearization technique is extended to aeroelastic problems. A loose fluid-structure coupling strategy is adopted. The structural oscillation of the blades is considered for eigen-modes but frequency is free to change during coupling resolution. The new approach is based on the building of a meta-model to describe the fluid dynamic behavior in order to solve directly the coupled fluid-structure system. This technique is validated on a standard high subsonic turbine configuration and takes advantage of a reduced computation time.

Dans cette thèse, nous proposons l'analyse numérique et le développement d'algorithmes partitionnés pour coupler l'écoulement du sang dans différents comparti- ments cardiovasculaires (3D-3D, 3D-0D) Dans une première partie, un problème couplé fluide-fluide est introduit. Sur l'interface qui sépare les domaines, des conditions aux limites de type Robin-Robin dérivées de la formulation d'interface de Nitsche sont considérées. Nous proposons différents schémas explicites dont la stabilité est analysée dans la norme de l'énergie. Des simulations numé- riques illustrent le potentiel des méthodes présentées. La deuxième partie propose des applications cardiovasculaires plus réalistes. Tout d'abord, un modèle d'ordre réduit pour les valves cardiaques est décrit. Sans traiter l'inter- action fluide-structure avec le sang, les valves sont remplacées par des surfaces agissant comme des résistances immergées dans le fluide. Des simulations numériques montrent l'efficacité et la robustesse de ce modèle. Pour finir, une formulation ALE est utilisée pour la résolution d'un modèle fluide sur un domaine mobile. Nous montrons qu'en ajoutant un terme consistent, une inégalité d'éner- gie stable peut être obtenue sans considérer aucune hypothèse de Loi de Conservation Géométrique. Le travail se termine avec des simulations numériques sur la dynamique du sang dans le ventricule gauche, couplé avec l'écoulement du sang dans l'aorte. / In this thesis we present the numerical analysis and the development of parti- tioned algorithms in order to couple the blood dynamics in different cardiovascular compart- ments (3D-3D, 3D-0D). In the first part a fluid-fluid coupled problem is introduced. On the interface between the domains Robin-Robin boundary conditions, derived from the interface Nitsche’s formulation, are considered. We suggest different staggered explicit schemes whose stability is analyzed in the energy norm. Extensive numerical experiments illustrate the accuracy of the methods presented. The second part deals with more realistic cardiovascular applications. First a reduced order model for the heart valves is described. Without dealing with fluid-structure interaction with the blood flow, the valves are replaced by immersed surfaces acting as resistances on the fluid. Numerical simulations show the efficiency and the robustness of this model in the framework of a fluid-fluid interaction scheme. In the end, an ALE formulation is used to solve a fluid model in a moving domain. We show that adding a suitable consistent term, a stable energy inequality can be obtained without considering any Geometric Conservation Laws. The work ends with numerical sim- ulations on blood dynamics in the left ventricle coupled with the blood flowing in the aorta.

Solides élastiques au repos, les fluides à seuil s’écoulent comme un liquide au-delà d’une certaine contrainte. Plusieurs applications industrielles concernent l’écoulement de ces fluides dans des milieux poreux. On peut citer par exemple les émulsions dans le processus de récupération du pétrole, les opérations de cimentation dans le sol, ou le nettoyage d’un sol contaminé par une boue. Pour ces applications, il est nécessaire de connaitre la pression nécessaire pour un débit voulu à la sortie du milieu poreux. Dans de tels cas, l’écoulement est perturbé par la complexité de la géométrie. Les modèles développés pour décrire la loi de Darcy supposent une loi rhéologique appliquée localement, mais ces modèles décrivent mal ce type d’écoulement. De plus, des effets complexes peuvent s’ajouter comme le glissement à la paroi ou la thixotropie. Dans cette thèse, nous étudions l’écoulement de carbopol (ETD 2050) à travers différentes géométries. Tout d’abord au rhéomètre, nous montrons que le fluide, sous certaines conditions, correspond bien à un fluide à seuil modèle. Nous démontrons que le protocole expérimental utilisé est très important et qu’un comportement thixotropique peut apparaitre s’il n’est pas respecté. Ce comportement apparait notamment lorsque le fluide reste sous le seuil, l’impact augmentant avec le temps d’attente. Ensuite, nous comparons la loi d’écoulement obtenue au rhéomètre à l’écoulement dans un canal droit obtenu par microfabrication. Nous montrons alors l’importance du glissement proche du seuil et ses conséquences sur la loi d’écoulement. Enfin nous étudions l’écoulement du carbopol dans un milieu poreux. Le milieu poreux de 5x5cm est obtenu par microfabrication. La largeur moyenne des canaux est égale à celle du canal droit. Nous avons développé une nouvelle méthode de mesure des champs de vitesse. Nous montrons l’apparition d’une chenalisation de l’écoulement à travers quelques canaux du milieu poreux. Nous comparons ensuite la loi d’écoulement du milieu poreux à celle obtenue dans le canal droit. On remarque que la vitesse d’écoulement est plus faible dans le milieu poreux que dans le canal droit. / Elastic solids at rest, yield stress fluids flow like a liquid beyond a certain stress. Many industrial applications required the flow of these fluids in porous media, for example: the emulsion flow in oil recovery processes, the cementing operations in the ground, or the cleaning of sludge in a contaminated soil. For many applications, it could be interesting to know the pressure required for a desired flow rate. In such cases, the flow behavior of the fluid is complicated by the complexity of the geometry. The models developed to describe Darcy's law assume a rheological law applied locally, but these models poorly describe this type of flow. Furthermore, complex effects can be added like the wall slip or the thixotropy. In this thesis, we study the flow of carbopol (ETD 2050) through different geometries. First we show that the fluid, for some conditions, corresponds to model yield stress fluids. The experimental protocol used is very important and a thixotropic behavior can appear if it is not respected. This behavior appears especially when the fluid remains below the yield stress, the impact increases with the waiting time. We then compare the flow law obtained by rheometer in a straight channel obtained by microfabrication. We show the importance of the wall slip near the yield stress and the impact on the flow law. Finally, using a new method to measure the velocity fields developed during this thesis, we study the flow of carbopol in a porous medium. This porous medium of 5x5cm is obtained by microfabrication. The mean width of the channels is equivalent to the one of the straight channel. We show the emergence of a channeling flow through some channels of the porous medium. We then compare the flow law of the porous medium to the one obtained in the straight channel. It can be observed that the flow rate is lower in the porous medium than in the straight channel.

Les tubes des générateurs de vapeur des centrales nucléaires vibrent sous l'effet d'écoulement eau/vapeur. Pour appréhender ce phénomène et le comprendre, des expériences à échelles réduites sont réalisées. La simulation numérique a montré son habilité à reproduire l'interaction fluide-structure sur ce type de géométrie pour des écoulements monophasiques. L'objectif est désormais de faire de même en écoulement diphasique et de caractériser les propriétés physiques du mélange liquide/gaz influant sur la vibration.Pour se faire, un code CFD avec une approche bi-fluide est utilisé. Une méthode dite de "Discrete forcing" est implémentée pour permettre le mouvement imposé de corps solides au sein d'un écoulement à plusieurs phases. Celle-ci est alros validée sur des cas simples et intégraux avec une comparaison systématique à des résultats expérimentaux ou théoriques.En se basant sur un algorithme implicite existant dans la littérature, un couplage fluide-structure utilisant cette méthode de suivi d'interface est implémenté. Validé sur des cas monophasiques et diphasiques, ce couplage offre désormais la possibilité de déplacer un solide en fonction des forces fluides diphasiques qui lui sont appliquées.Les différentes méthodes numériques présentes dans NEPTUNE_CFD sont ensuite évaluées pour un écoulement fréon/fréon au travers d'un faisceau de tubes inclinés. La nécessité d'utiliser des modèles dit "multi-régime" est mis en avant.Afin de déterminer l'influence sur l'écoulement des différentes propriétés physiques d'un mélange diphasique, plusieurs cas simples sont réalisés.Finalement, l'application industrielle cible, un écoulement eau/fréon dans un faisceau de tubes à pas carré, est simulée et comparée à un écoulement en conditions réelles (eau/vapeur à 70 bar). Les vibrations induites par écoulement monophasique puis diphasique sont correctement reproduites sur des cas dit de "faisabilité". / In nuclear power plants, steam generator tubes vibrate because of steam/water cross-flows. In order to understant this phenomenon, reduced-scale experiments are performed. Numerical simulations have shown their ability to accurately reproduce the vibration induced by a single phase flow in a tube bundle. The aim of the present work is to do the same with two-phase flow and to characterize the effect of the mixture physical properties on vibration.To do so, a CFD code based on a two-fluid approach is used. A "discrete forcing" method is implemented in order to allow solid body motion in a two-phase flow. The validation is performed with simple and industrial cases using experimental and theoretical results.Using an existing implicit algorithm, a fluid-structure coupling based on the developed interface tracking method is implemented. Validated for single and two-phase flows, it is now possible to have solid motion induced by fluid forces.The different numerical models dedicated to two-phase flows are then evaluated on a freon/freon flow across an inclined tube bundle. The use of a multi-regime model is required. In order to investigate the role of the different physical properties on the vibration, three simple studies are performed.Finally, the industrial application, a freon/water flow across a square pitch tube bundle, is performed. First, it is compared to a steam/water flow in order to characterize the discrepancies when we are using a modeling mixture. Then, the vibration induced by single- and two-phase flows is reproduced by the developed method on feasibility test cases.