Simulation de l'érosion de cavitation par une approche CFD-FEM couplée

This research is devoted to understanding the physical mechanism of cavitation erosion in compressible liquid flows on the fundamental scale of cavitation bubble collapse. As a consequence of collapsing bubbles near solid wall, high pressure impact loads are generated. These pressure loads are believed to be responsible for the erosive damages on solid surface observed in most applications. Our numerical approach begins with the development of a compressible solver capable of resolving the cavitation bubbles in the finite-volume solver YALES2 employing a simplified homogenous mixture model. The solver is extended to Arbitrary Lagrangian-Eulerian formulation to perform fluid structure interaction simulation with moving mesh capabilities. The material response is resolved with the finite element solver Cast3M, which allowed us to perform one-way and two-way coupled simulations between the fluid and solid domains. In the end, we draw comparisons between 2D and 3D vapor bubble collapse dynamics and compare them with experimental observations. The estimated pressure loads on the solid wall and different responses of materials for attached and detached bubble collapses are discussed. Finally, the damping of pressure loads by different materials is identified with two-way coupled fluid-structure interaction.Ce travail de recherche est dédié à la compréhension des mécanismes physiques de l’érosion de cavitation dans un fluide compressible à l’échelle fondamentale de l’implosion d’une bulle de cavitation. Suite à l’implosion d’une bulle de vapeur à proximité d’une surface solide, des très hautes pressions sont générées. Ces pressions sont considérées responsables de l’endommagement (érosion) des surfaces solides observé dans la plupart des applications. Notre approche numérique démarre avec le développement d’un solveur compressible capable de résoudre les bulles de cavitation au sein du code volumes finis YALES2 en utilisant un simple modèle de mélange homogène des phases fluides. Le solveur est étendu à une approche ALE (Arbitraire Lagrangien Eulérien) dans le but de mener des simulations d’interaction fluide-structure sur un maillage mobile. La réponse du matériau solide est calculée avec le code de calcul éléments finis Cast3M, et nous a permis de mener des simulation avec un couplage d’abord monodirectionnel, ensuite bidirectionnel, entre le fluide et le solide. On compare des résultats obtenus à deux dimensions, puis à trois, avec des observations expérimentales. On discute les chargements de pression estimés, et les réponses de différents matériaux pour des implosions de bulle à des différentes distances de la surface. Enfin, à travers l’utilisation de simulations avec couplage bidirectionnel entre fluide et solide, on identifie l’amortissement des chargements de pression pour les différents matériaux