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.

L'utilisation croissante des sources d'énergie renouvelable avec une grande volatilité de production, comme l'énergie éolienne et solaire, conduit à des fluctuations dans le réseau électrique qui doivent être compensées. Pour cette raison les machines hydrauliques, turbines et pompes, sont plus souvent opérées dans les régimes de fonctionnement hors fonctionnement nominal et la fréquence des phases de démarrage et arrêt augmente. Ce type de fonctionnement peut avoir des conséquences importantes sur le cycle de vie des machines. Il est donc essentiel de prendre en compte l'écoulement dans les phases transitoires lors de la conception de la machine et la simulation numérique des écoulements est un outil adapté pour cela. La présente étude a pour objectif de développer une méthode de couplage flexible qui combine la méthode à maillage volumes finis (VF) et la méthode sans maillage Smoothed Particle Hydrodynamics - Arbitrary Lagrange Euler (SPH-ALE). Cette méthode couplée peut être utilisée comme outil pour l'investigation des phénomènes transitoires dans les machines hydrauliques. SPH-ALE est particulièrement bien adapté aux simulations des écoulements fortement dynamiques avec des géométries mobiles mais elle a des difficultés pour calculer des forts gradients de pression et vitesse. Un raffinement de particules est difficile à implémenter, surtout si les particules doivent être raffinées de manière anisotrope. Les méthodes volumes finis (VF) sont établies pour les simulations numériques d'écoulements grâce à leur stabilité et précision. Par contre, elles peuvent être lourdes pour les simulations avec des géométries mobiles et demandent souvent une interface entre des parties mobiles et statiques du maillage ce qui génère des erreurs supplémentaires. Pour combiner les deux approches complémentaires, une méthode de couplage a été développée qui décompose le domaine de calcul en zones où la vitesse et la pression sont calculées par la méthode VF, en zones où elles sont obtenues par SPH-ALE et en zones de recouvrement où les informations sont transférées de la zone VF à la zone SPH et inversement. Dans les zones de recouvrement les points de calcul VF sont utilisés comme voisins pour l'intégration en espace des particules SPH. Aux limites du maillage VF la vitesse et la pression sont interpolées des particules SPH, similairement aux méthodes Chimére des maillages recouvrants. Un logiciel SPH-ALE existant du groupe ANDRITZ est utilisé pour cette étude. Un solveur VF faiblement compressible est implémenté dans ce logiciel. Le solveur discrétise la même forme des équations de Navier-Stokes que le solveur SPH-ALE. Des solveurs de Riemann avec des états reconstruits par la méthode MUSCL sont employés. En outre, le solveur SPH-ALE est amélioré et adapté aux écoulements internes. Pour cette raison des conditions à l'entrée et à la sortie du type subsonique sont implémentées. Du plus, une méthode de correction du gradient de la fonction kernel est présentée qui améliore la précision du champ de pression, notamment si les particules ne sont pas distribuées régulièrement. La méthode couplée est validée à l'aide des cas test académiques en unidimensionnel et en bidimensionnel, comme le cas de tube à choc, les tourbillons de Taylor-Green et l'écoulement autour d'une aube symétrique du type NACA avec des particules en description eulérienne. En outre, le couplage offre la possibilité d'imposer des conditions à la sortie aux particules lagrangiennes. La méthode est appliquée aux simulations d'écoulement transitoire en 2D avec des particules qui se déplacent en suivant les géométries mobiles. / The increased use of intermittent forms of renewable energy like wind and solar energy produces fluctuations in the electric grid that have to be compensated. For this reason, hydraulic machines like turbines and pumps are more often operated under non-conventional operating conditions and are submitted to frequent starts and stops. This type of operating conditions has important consequences on the life cycle of the machines. It is thus of paramount importance that transient flows at off-design conditions are properly taken into account in the design phase and numerical simulation is an appropriate way to do so. The present study aims at developing a flexible coupling method of the meshbased Finite Volume Method (FVM) and the meshless Smoothed Particle Hydrodynamics - Arbitrary Lagrange Euler (SPH-ALE) method, which can be used as a tool for the investigation of transient phenomena in hydraulic machines. SPH-ALE is very well adapted for the simulation of highly dynamic flows with moving geometries but has difficulties to correctly represent rapidly changing gradients of the field variables. Particle refinement is difficult to implement, especially if particles are refined in an anisotropic way. FV methods are well established in CFD because of their accuracy and stability. However, they can be tedious for simulations with moving geometries and often necessitate an interface between moving and static parts of the mesh which introduces additional errors. To overcome the shortcomings of both methods, a coupling method is developed that uses a decomposition of the computational domain into regions where the physical field variables are computed by the FV method, regions where they are computed by SPH-ALE and overlapping regions where the information is transferred from the FV domain to the SPH domain and vice versa. In the overlapping regions FV calculation points are used as neighbors for the SPH integration in space. At the boundaries of the FV mesh, velocity and pressure are interpolated from the SPH particles by means of scattered data interpolation techniques, similarly to Chimera methods for overlapping grids. For this study, an existing SPH-ALE software of the ANDRITZ Group is used. A weakly compressible FV solver is implemented into this software that discretizes the same form of the Navier-Stokes equations than the SPH-ALE solver. Similar to the present SPH-ALE method, Riemann solvers with reconstructed states, obtained by a MUSCL scheme, are employed. Moreover, adaptations and improvements of the SPH-ALE solver itself are made, which are important for the coupling and for the simulation of internal flows in general. Thus, subsonic inlet and outlet conditions are implemented. Furthermore, a correction method of the kernel gradient is presented that ensures zeroth order consistency of the SPH-ALE approximation of the divergence of the convective fluxes. The correction improves greatly the SPH pressure field on non-uniform particle distributions. The implemented coupled method is successfully validated by means of inviscid academic one-dimensional and two-dimensional testcases like a shock tube case, Taylor-Green vortices and the flow around a symmetric NACA airfoil with particles in Eulerian description. Furthermore, the coupling provides a possibility to implement outlet boundary conditions to Lagrangian moving SPH particles. It is then applied to the simulation of transient flows in rotor stator systems in 2D with moving particles.