The far-field underwater explosion (UNDEX) problem is a complicated problem dominated by two phenomena: the shock wave traveling through the fluid and the cavitation in the fluid. Both of these phenomena have a significant effect on the loading of ship structures subjected to UNDEX. An approach to numerically modeling these effects in the fluid and coupling to a structural model is using cavitating acoustic finite elements (CAFE) and more recently cavitating acoustic spectral elements (CASE). The use of spectral elements in CASE has shown to offer the greater accuracy and reduced computational expense when compared to traditional finite elements. However, spectral elements also increase spurious oscillations in both the fluid and structural response. This dissertation investigates the application of CAFE, CASE, and a possible improvement to CAFE in the form of a finite element flux-corrected transport algorithm, to the far-field UNDEX problem by solving a set of simplified UNDEX problems. Specifically we examine the effect of increased oscillations on structural response and the effect of errors in cavitation capture on the structural response which have not been thoroughly explored in previous work. The main contributions of this work are a demonstration of the problem dependency of increased oscillations in the structural response when applying the CASE methodology, the demonstration of how the sensitivity of errors in the structural response changes with changes in the structural model, a detailed explanation of how error in cavitation capture influences the structural response, and a demonstration of the need to accurately capture the shape and magnitude of cavitation regions in the fluid in order to obtain accurate structural response results. / Ph. D.

Ultrasonic humidifiers use sound vibrations at frequencies higher than can be heard by humans (> 20,000 Hz) to generate aerosolized water also have potential for inducing sonochemical reactions for chemicals present in water. This research focuses on examining oxidants formed within ultrasonic humidifiers, as well as the oxidants effects of contaminants in water used in the systems. Hydroxyl radicals were found using DMPO as a spin trap. Caffeine and 17β-estradiol, as pharmaceutical contaminants of drinking water, were both emitted from the humidifier when present in the water reservoir and would enter breathing air. Emitted 17β-estradiol was found at 60% of the initial concentration filled in the ultrasonic humidifier after 480 minutes. Caffeine exhibited less degradation than 17β-estradiol. Degradation of both pharmaceuticals was attributed to ultrasonic processes, most likely oxidation with hydroxyl radicals produced. Bromide as a contaminant of the fill water was found to remain constant over time. / MS

La cavitation peut avoir lieu dans les turbines hydrauliques. Ce phénomène se produit lorsque les bulles de vapeur s’effondrent à proximité de la surface de la machine. Ceci entraîne des conséquences négatives, telles que l’érosion, affectant ainsi les performances de la machine. L’effondrement d’une bulle de gaz non-condensable dans l’eau est simulé en utilisant la méthode sans maillage SPH-ALE, qui intègre un modèle pour simuler les écoulements compressibles et multiphases. Le modèle résout les équations de conservation de masse, de quantité de mouvement et d’énergie du système d’Euler, en utilisant l’équation d’état de Stiffened Gas pour l’eau et l’équation d’état de gaz parfait pour le gaz non-condensable à l´ıntérieur de la bulle. Les deux phases sont modélisées comme compressibles et le changement de phase n’est pas considéré. La caractéristique sans maillage de la méthode SPH-ALE permet le calcul des écoulements diphasiques où l’interface est nettement définie. Pour les applications de cavitation, où le nombre de Mach atteint des valeurs de 0.5, la distribution de particules doit être corrigée. Cela est réalisé grâce à la fonctionnalité ALE. Le modèle compressible a été validé à l’aide de configurations monodimensionnelles, comme le cas du tube à choc pour des écoulements monophase et multiphases. L’effondrement de la bulle près d´une paroi a été abordé comme le mécanisme fondamental qui produit des dégâts. Son comportement général se caractérise par la formation d’un micro jet d’eau et par l’effondrement de la bulle sur elle-même. Le phénomène est analysé en tenant compte des principaux paramètres qui le régissent, comme la distance initiale entre le centre de la bulle et la paroi (H0), la taille de la bulle (R0) et le taux de pression qui entraîne l’effondrement (pw/pb). Il est démontré que l’intensité de l’effondrement dépend principalement du rapport de pression entre le liquide et la bulle (pw/pb). De plus, quatre indicateurs, comme la pression en paroi, l’impulsion, la pression du coup de bélier et la vitesse du micro jet d’eau, servent à déterminer le chargement. Cette analyse indique qu’une bulle initialement située à une distance inférieure à H0/R0 = 2 présente un haut potentiel d’endommagement. Afin de prédire cet endommagement, la mécanique du solide est analysée à l’aide de simulations d’interaction fluide-structure. On obtient que le matériau réagit aux charges hydrauliques en ayant des zones de compression et de traction. Ceci suggère qu’un mécanisme de fatigue entraîne le phénomène d’endommagement. En plus, on constate que les contraintes les plus importantes sont situées sous la surface du matériau, indiquant que cette zone peut être sujette à une déformation plastique. / Hydraulic turbines can experience cavitation, which is a phenomenon occurring when vapor bubbles collapse in the vicinity of the machine’s surface. This phenomenon can lead to negative consequences, such as erosion, that affect the machine’s performance. The compression of a non-condensable gas bubble in water is simulated with the Smoothed Particle Hydrodynamics method following the Arbitrary Lagrange Euler approach (SPHALE), where a compressible and multiphase model has been developed. The model solves the mass, momentum and energy conservation equations of the Euler system using the Stiffened Gas EOS for water and the ideal gas EOS for the non-condensable gas inside the bubble. Both phases are modeled as compressible and the phase change is not considered. The meshless feature of the SPH-ALE method allows the calculation of multiphase flows where the interface is sharply defined. For cavitation applications, where the Mach number reaches values of 0.5, the distribution of particles must be corrected, which is achieved by the ALE feature. The compressible model was validated through monodimensional configurations, such as shock tube test cases for monophase and multiphase flows. The bubble compression close to the wall has been addressed as the fundamental mechanism producing damage. Its general behavior is characterized by the formation of a water jet and by the collapse of the bubble by itself. The phenomenon is analyzed by considering the major parameters that govern the bubble collapse dynamics, such as the initial distance between the bubble center and the wall (H0), the bubble size (R0), and the collapse driven pressure ratio (pw/pb). It is shown that the intensity of the collapse depends mainly on the pressure ratio between the liquid and the bubble (pw/pb). As well, four indicators, such as the pressure at the wall, the impulse, the water-hammer pressure and the water jet velocity, are used to determine the loading. This analysis gives that the bubble initially located at a distance lower than H0/R0 = 2 presents high potential to cause damage. In order to predict the damage due to the bubble collapse, the solid mechanics is analyzed through fluid-structure interaction simulations. It is obtained that the material reacts to the hydraulic loads by having compression and traction zones, suggesting that a fatigue mechanism drives the damage phenomenon. Additionally, it is found that the highest stresses are located below the material surface, indicating that this zone may reach plastic deformation.

Chaque année, le nombre de demandeur d’organe augmente en France comme dans le reste du monde. Pour combattre ce fléau, il existe aujourd’hui des technologies permettant d’imprimer du vivant, telle que la Bioimpression Assistée par Laser (LAB). Robuste et précise, cette méthode s’appuie sur les propriétés d’interaction laser-matière pour éjecter une bio-encre constituée de cellules vivantes. Pour éviter l’utilisation d’une couche absorbante sacrificielle, généralement utilisée, on focalise directement un faisceau laser dans la bioencre afin de générer un plasma puis une bulle de cavitation. La position de cette bulle est essentiellement maitrisée pas la longueur d’onde, et sa taille est gérée par l’énergie et la durée d’impulsion du laser. Ce sont les facteurs clés pour maîtriser l’éjection de matière biologique. Cependant, l’inhomogénéité locale apportée par les cellules perturbe l’impact du laser et donc la reproductibilité des jets, mais une fois imprimées, ces cellules sont viables et permettent de reconstruire des tissus vivants. / Every year, the transplant waiting list gets bigger in France as in all over the world. To fight this curse, Bioprinting makes organ printing possible, especially with Laser Assisted Bioprinting (LAB). Robust and precise, this method use laser-matter interaction to eject a bioink made of living cells. To avoid the use of absorbing sacrificial layer, we directly focalize a laser beam into a living cells bioink, to create a plasma then a cavitation bubble. Its position, which is mainly driven by laser wavelength, and its size, managed by the energy and pulse duration, are the most important keys to control liquid jet ejection. However, the laser energy deposition and jet ejection is disturbed because of cells local concentration disparity, but when cells are printed, they are still viable and able to reconstruct living tissues.

L'application d'ultrasons sur une tumeur, où des liposomes se sont accumulés, permet potentiellement de libérer le médicament mais aussi d'en favoriser l'absorption dans les cellules. La cavitation inertielle ultrasonore est le phénomène pressenti pour la libération sous ultrasons de médicament encapsulé dans de petits liposomes solides. Elle est très dépendante des conditions expérimentales et peut-être intense et imprévisible. L'objectif principal du travail réalisé dans le cadre de cette thèse est de contrôler et quantifier la cavitation inertielle, pour induire le largage de médicaments encapsulés dans des liposomes. Dans cette optique, une dose de cavitation inertielle (DC), basée sur le filtrage du bruit large bande émis lors de ce régime de cavitation, est mise au point in vitro pour suivre le largage de médicament encapsulé. Sous divers régimes d'ultrasons pulsés, la DC a été validée en dosant chimiquement les radicaux hydroxyles générés lors de l'implosion des bulles. Les tests menés sur diverses formulations de liposomes contenant de la doxorubicine (dox) ont montrés une haute corrélation entre le taux de largage de dox et la DC permettant de conclure que la cavitation inertielle est impliquée dans ce largage. Le rôle de la température sur la production de radicaux hydroxyles et la libération de dox a également été exploré. Les expériences réalisées ont permis de sélectionner les formulations les plus sensibles aux ultrasons pour les tester sur des rats implantés avec des tumeurs prostatiques. Après plusieurs expériences in vivo menées avec différents dispositifs ultrasonores et formulations de liposomes, le bénéfice du traitement combiné a pu être démontré.

Pour assurer le transport de la sève des racines vers les feuilles, les plantes vasculaires génèrent de très fortes dépressions dans le liquide, pouvant atteindre -200 bar, au niveau des feuilles. Cette dépression « tire » sur la colonne d'eau contenue dans l'appareil vasculaire de l'arbre. La cohésion de l'eau maintient la sève sous forme liquide. Cet état métastable peut se rompre : des bulles de cavitation apparaissent. Elles créent un « bouchon » d'air dans le réseau hydraulique de la plante et gênent la circulation de la sève. C'est ce que l'on appele l'embolie. Si ce phénomène se généralise, il peut provoquer la mort de la plante.Ce travail de thèse est consacré à 'invasion d'air dans des réseaux hydrauliques naturels ou artificiels initialement à pression négative. Nous avons d'abord étudié l'embolie dans les feuilles. Nous avons développé une technique novatrice permettant de relever la propagation spatiale de l'embolie dans le réseau hydraulique des feuilles. Nous montrons que l'embolie, quelque soit l'espèce, se propage par à-coups des plus grosses nervures aux plus petites.Afin de comprendre les lois physiques sous-jacentes, nous utilisons deux systèmes modèles. Nous réalisons d'abord des réseaux artificiels dans un hydrogel reproduisant les caractéristiques de la circulation de la sève ascendante. Après la relaxation de la tension dans le réseau par l'apparition de la bulle, nous observons des oscillations de surface et une croissance lente de la bulle, liée à l'évacuation de l'eau à travers l'hydrogel. Cette croissance peut atteindre un régime quasi-stationnaire. Ce systèmes ne nous permettant pas de reproduire toutes les caractéristiques géométriques du xylème, nous présentons une modélisation informatique reposant sur l'analogie entre réseaux hydrauliques et électrocinétique. Nous reproduisons les caractéristiques du xylème dans lequel circule la sève : les éléments conducteurs sont reliées par les ponctuations, des valves protégeant la plante de l'embolie. Nous retrouvons les à-coups caractéristiques de la propagation de l'embolie dans les feuilles.Enfin, nous discutons l'application des résultats précèdents dans le cas du bois et nous présentons quelques résultats obtenus sur du pin sylvestre. / To assure the transport from the roots to the leaves, vascular plants create strong depressions in the sap, next to -200 bars. This depression pulls the water column contained by the tree vascular system. The water cohesion keeps the sap under liquid state. This metastable state can breaks: cavitation bubbles appear. They create an air plug inside the plant hydraulic network and impede sap flow. This phenomena called embolism could lead to the plant death by preventing the sap transport.This thesis is dedicated to the air invasion into hydraulics networks under negative pressure. First, we study the leaf embolism. We developed a new technique which allows us to record the spatial propagation of embolism in leaves hydraulic network. We show that the embolism propagates by steps from biggest veins to smallest veins.Next, in order to understand the underlying physical laws, we use two model systems. We build artificial networks in a hydrogel which mimics the sap flow characteristics. After the relaxation of the negative pressure in the network by the nucleation of a bubble, we observe surface oscillations and the slow growth of the bubble. This growth is linked to the water transport through the hydrogel and can reach a stationary regime.As we are not able to reproduce all the characteristics of the leaf network with the hydrogel, we create a computer modeling based on the Ohm analogy between hydraulics networks and electrical circuits. We reproduce the specific features of the xylem which transport the sap: the conduits are linked by pits, small valves which limit the progression of the embolism. We were able to recover the distinctiveness steps in embolism.Finally, we discuss the application of the preceding results to wood and we present some results on Pinus sylvestris.

La présente étude a été réalisée dans le cadre d'une collaboration entre le groupe de recherche GREM3 du laboratoire LAPLACE de Toulouse et la société de recherche contractuelle SINAPTEC située à Lille. SINAPTEC conçoit et développe des systèmes ultrasonores piézoélectriques résonants pour un grand nombre d'applications industrielles et a pour volonté d'améliorer sa compréhension des phénomènes mis en jeu ainsi que les performances et les fonctionnalités de sa technologie actuelle. L'objectif de cette étude est d'exploiter la sensibilité des transducteurs piézoélectriques résonants afin de leur conférer une fonction de capteur en vue de mieux contrôler leur action sur le milieu à traiter par le suivi de l'évolution de l'impédance électromécanique. Actuellement, pour le pilotage et le contrôle des transducteurs US de puissance dédiés aux applications en milieu liquide - qui font l'objet de la présente étude - plusieurs difficultés se posent. Tout d'abord la première difficulté, inhérente à la nature même des transducteurs piézoélectriques de puissance, est la dérive de leurs caractéristiques avec l'échauffement. L'échauffement dans le transducteur provoque une dérive de sa fréquence de résonance et de ses propriétés mécaniques qui est partiellement solutionnée par un asservissement en fréquence (et/ou en puissance) mais il n'en reste pas moins qu'une dégradation des performances de l'actionneur est observée. La deuxième difficulté réside dans le fait que le transducteur plongé dans une charge liquide est difficile à contrôler car l'application des ondes ultrasonores de puissance dans ce milieu provoque l'apparition de microbulles de vapeur, c'est le phénomène de cavitation qui est extrêmement complexe et fortement non linéaire. L'objectif est donc de trouver une méthode simple, peu consommatrice en ressources pour être implémentée dans la nouvelle génération de générateurs électroniques de SINAPTEC et qui permette d'extraire des informations pertinentes sur le milieu liquide - comme l'intensité de la cavitation par exemple - mais aussi sur l'état du transducteur et son évolution pendant qu'il fonctionne à forte puissance dans des conditions normales d'utilisation. Le principe de la méthode proposée repose sur le calcul de l'impédance de charge du transducteur en fonctionnement et son suivi en « temps réel »par mesure des signaux électriques et à partir de ses caractéristiques à vide. Pour montrer la faisabilité d'une telle méthode, une modélisation analytique par le modèle des matrices de transfert ainsi qu'une modélisation numérique, à l'aide du logiciel ANSYS, ont été réalisées afin de qualifier l'effet de la charge liquide en régime stable et cavitant sur l'impédance électrique du transducteur. Deux transducteurs ont été modélisés : un système ultrasonore standard de SINAPTEC, classiquement utilisé pour les applications en milieu liquide, et un transducteur spécialement conçu pour être plus sensible à la charge et plus facilement pilotable. Cette modélisation a été complétée par des expérimentations afin de confirmer la faisabilité de la méthode en milieu cavitant. La modélisation ainsi que les expérimentations ont permis de confirmer que les propriétés du milieu liquide en contact avec le transducteur sont corrélées à l'impédance électrique de ce dernier (et à sa fréquence de résonance). De ce fait, il devient alors possible, en suivant/observant l'impédance électrique du transducteur, de déterminer certaines propriétés du milieu et ainsi permettre un meilleur suivi de l'action du transducteur sur sa charge. / This study is the result of a collaboration between the research group GREM3 of the LAPLACE Laboratory in Toulouse and SINAPTEC, a research company located in Lille. SINAPTEC designs and develops resonant piezoelectric ultrasonic systems for many industrial applications and wants to improve the performance and functionality of its current technology. The objective of this study is to exploit the sensitivity of resonant piezoelectric transducers to use them as sensors in order to better control their action on the medium to be treated by monitoring the evolution of their electromechanical impedance. ABSTRACT : Currently, several difficulties arise for the control of high power ultrasonic transducers dedicated to applications in liquid media - which are the subject of this study. The first difficulty - inherent in the nature of piezoelectric transducers- is the shift of their characteristics with heating. Heating in the transducer causes a drift of the resonance frequency and its mechanical properties which is partially solved by a frequency control loop, but a performance degradation of the actuator is still observed. The second difficulty lies in the fact that the transducer immersed in a liquid load is difficult to control because the application of high power ultrasonic waves in the medium generates micro-bubbles in the liquid: it is the phenomenon of cavitation, extremely complex and highly non-linear. The objective is to find a simple method, to be implemented in the next generation of electronic generators of SINAPTEC, that allows to extract relevant information on the liquid medium (such as the intensity of cavitation for instance) but also on the state of the transducer and its evolution during operation at high power under normal conditions of use. The principle of the proposed method is based on the calculation of the transducer load impedance while functioning by measuring electrical signals and using its unloaded characteristics. To show the feasibility of such an approach, the analytical model of the transfer matrices as well as a finite element model ( using ANSYS software) were implemented to characterize the effect of stable and cavitating liquid loads on the electrical impedance of the transducer. Two transducers were modeled: a standard ultrasound system of SINAPTEC, typically used for applications in liquid media, and a transducer specifically designed to be more sensitive to the load and more easily controllable. These model were complemented by experiments to confirm the feasibility of the method in a cavitating environment. Modeling and experiments have confirmed that the properties of the liquid in contact with the transducer are correlated with the electrical impedance of the latter (and its resonance frequency). Thus, it becomes possible to determine certain properties of the medium by monitoring the electrical impedance of the transducer and it should allow a better control of the action of the transducer on the load.

La simulation des écoulements cavitants est confrontée à des difficultés de modélisation et de résolution numérique provenant des caractéristiques particulières de ces écoulements : changement de phase, gradient de masse volumique important, variation du nombre de Mach, turbulence diphasique, instationnarités. Dans cette thèse, nous nous sommes appliqués à dériver proprement le modèle de mélange homogène 1-fluide couplé à une modélisation RANS de la turbulence. A partir des termes contenus dans ces équations et de la nature des écoulements cavitants étudiés, plusieurs modèles de turbulence basés sur la notion de viscosité turbulente ont été testés : modèles faiblement non-linéaires (corrections SST et de réalisabilité), ajout des termes de turbulence compressible, application de la correction de Reboud, modèles hybrides RANS/LES (DES, SAS). Ces modèles ont été incorporés dans un code compressible qui fait appel à une résolution implicite en pas de temps dual des équations de conservation avec une technique de pré-conditionnement bas-Mach pour traiter les zones incompressibles. Les simulations 2D et 3D ont porté sur deux géométries de type Venturi caractérisées par la présence d’une poche de cavitation instationnaire due à l’existence d’un jet rentrant liquide/vapeur le long de la paroi. Elles montrent que l’ensemble des modèles sont capables de capturer le jet rentrant. En revanche, la dynamique de la poche varie entre les modèles et le manque de données expérimentales ne permet pas de discriminer les modèles entre eux. Il apparaît à la vue des résultats que les approches avec la correction de Reboud ou les modèles SAS améliorent la simulation des écoulements. / The computation of cavitating flows is a challenging issue due to the characteristics of these flows : phase transition, large density gradient, Mach number variation, interaction between phases and turbulent flow, unsteadiness. In the present study, we performed a derivation of the one-fluid compressible homogenous model coupled with a RANS approach for the turbulent flow. From these equations and the nature of the cavitating flows, several models based on the eddy viscosity assumption have been tested : weakly non-linear models (SST and realisability corrections), compressible turbulence models, hybrid RANS/LES turbulence models (SAS, DES) and the Reboud correction. All the models are implemented in a compressible code, which solves the equations using an implicit dual-time stepping method coupling with a pre-conditionning technique for the incompressible area. 2D and 3D computations are performed on two Venturi geometries characterized by an unsteady cavitation sheet with a liquid/vapor re-entrant jet. All the models are able to capture the re-entrant jet. Nevertheless, the dynamic behaviour differs from one model to another and the lack of experimental data prevents to discriminate the models between them. From the results, the computations with the SAS model and the Reboud correction improve the prediction of the flow.

Le développement de structures portantes flexibles dans le domaine naval, telles que les hélices ou les safrans, pose de nouveaux problèmes de dimensionnement. Cette thèse a pour but de développer une méthode de dimensionnement validée par des essais pour des structures portantes déformables soumises à des écoulements, éventuellement diphasiques de type cavitant. Les essais sont réalisés sur un hydrofoil de type NACA66-312(mod.), fabriqué en polyacetate, au sein du tunnel hydrodynamique de l'Institut de Recherche de l'Ecole Navale. Lors des essais, des mesures de déformations du profil portant ainsi que de niveaux vibratoires sont réalisées. Une méthode numérique couplant un code structure éléments finis (ANSYS Mechanical) avec un code fluide volumes finis (ANSYS CFX) par une méthode partitionnée, itérative, synchrone et séquentielle, laquelle est validée en terme de prédiction du déplacement et des contraintes pour des écoulements subcavitants dans un premier temps, puis pour des écoulements cavitants stables et instables. / The design of flexible lifting bodies in the naval industry, such as propelleror rudders, create some new design problems. This thesis proposes a numerical method validated by experimental comparison for solving the case of lifting bodies loaded by flow with or whitout cavitation. The tests are carried out in the hydrodynamic tunnel of the French Naval Academy Research Institute, on a polyacetate flexible hydrofoil NACA66-312 (mod.). During tests, strains and vibrations are measured for comparisons with numerical results. The numerical method uses a sequential synchrone iterative partitionned coupling betweena structural finite-element code (ANSYS Mechanical) and a finite-volume code (ANSYS CFX). Good agreement between numerical and experimental results for displacements, and stresses of the structure is highlighted. For the cavitating flow, a good agreement for the cavitation dynamic is observed and the stresses are evaluated with satisfying accuracy.

La cavitation est l'un des phénomènes physiques les plus contraignants influençant les performances des machines hydrauliques. Il est donc primordial de savoir prédire son apparition et son développement, et de quantifier les pertes de performances qui lui sont associées. L'objectif de ce travail est de développer un algorithme 3D instationnaire pour la simulation numérique de la cavitation dans le code industriel « Code_Saturne ». Il est basé sur la méthode à pas fractionnaires et préserve le principe du minimum/maximum sur le taux de vide. Un solveur implicite, basé sur l'équation de transport du taux de vide couplée avec les équations Navier-Stokes est proposé. Un traitement numérique spécifique des termes sources de cavitation permet d'obtenir des valeurs physiques du taux de vide (entre 0 et 1) sans aucune limitation artificielle. L'influence des modèles de turbulence RANS sur la simulation de la cavitation est étudiée sur deux types de géométries 2D (Venturi et Hydrofoil). Cela confirme que la modification de Reboud et al. (1998) appliquée aux modèles à viscosité turbulente à deux équations, k-epsilon et k-omega-SST, permet de reproduire les principales caractéristiques du comportement instationnaire de la poche de cavitation. Le modèle du second ordre RSM-SSG, basé sur le transport des contraintes de Reynolds, se révèle capable de reproduire le comportement instationnaire de l'écoulement sans aucune modification arbitraire. Les effets tridimensionnels intervenant dans les mécanismes d'instabilité de la poche sont également analysés. Ce travail nous permet d'aboutir à un outil numérique, validé sur des configurations d'écoulements cavitants complexes, afin d'améliorer la compréhension des mécanismes physiques qui contrôlent les effets instationnaires tridimensionnels intervenants dans les mécanismes d'instabilité. / Cavitation is one of the most demanding physical phenomena influencing the performance of hydraulic machines. It is therefore important to predict correctly its inception and development, in order to quantify the performance drop it induces, and also to characterize the resulting flow instabilities. The aim of this work is to develop an unsteady 3D algorithm for the numerical simulation of cavitation in an industrial CFD solver « Code_saturne ». It is based on a fractional step method and preserves the minimum/maximum principle of the void fraction. An implicit solver, based on a transport equation of the void fraction coupled with the Navier-Stokes equations is proposed. A specific numerical treatment of the cavitation source terms provides physical values of the void fraction (between 0 and 1) without including any artificial numerical limitation. The influence of RANS turbulence models on the simulation of cavitation on 2D geometries (Venturi and Hydrofoil) is then studied. It confirms the capability of the two-equation eddy viscosity models, k-epsilon and k-omega-SST, with the modification proposed by Reboud et al. (1998) to reproduce the main features of the unsteady sheet cavity behavior. The second order model RSM-SSG, based on the Reynolds stress transport, appears able to reproduce the highly unsteady flow behavior without including any arbitrary modification. The three-dimensional effects involved in the instability mechanisms are also analyzed. This work allows us to achieve a numerical tool, validated on complex configurations of cavitating flows, to improve the understanding of the physical mechanisms that control the three-dimensional unsteady effects involved in the mechanisms of instability.