Formulations ALE et SPH en Interaction Fluide-Structure : impact d’ondes de choc

ALE and SPH formulations for Fluid Structure Interaction : shock waves impact

• Auteur / Autrice : Ramzi Messahel
• Direction : Mhamed Souli
• Type : Thèse de doctorat
• Discipline(s) : Mécanique, Énergétique et Sciences des matériaux
• Date : Soutenance le 30/03/2016
• Etablissement(s) : Lille 1
• Ecole(s) doctorale(s) : École doctorale Sciences pour l'ingénieur (Lille)
• Partenaire(s) de recherche : Laboratoire : Laboratoire de mécanique de Lille (LML)

Mots clés

Mots clés contrôlés

Coup de bélier

Écoulement polyphasique

Ondes de choc

Interaction fluide-structure

Mots clés libres

Arbitraire Lagrangien-Eulérien (ALE)

Hydrodynamique des particules lisses (SPH)

Résumé

Ce travail de thèse porte sur l’étude numérique de la propagation d’ondes de choc dans les écoulements compressibles multiphasiques et en interaction (fluide-structure). Deux approches sont étudiées pour la résolution numérique de la partie fluide : L’approche ALE (Arbitrary Lagrangian Eulerian) et l’approche lagrangienne SPH (Smoothed Particle Hydrodynamics) ; la partie structure, quant à elle, est résolue par une approche classique EF (Éléments finis). L’étude des méthodes ALE et SPH constituent les deux principaux axes de recherche. La problématique des coups de bélier dans l’ingénierie nucléaire est abordée dans cette thèse. Lors d’un coup de bélier, les nombreuses réflexions d’ondes de choc dans les tuyauteries nucléaires peuvent faire baisser la pression de l’eau en dessous de sa pression de saturation et générer localement de la cavitation. Le modèle HEM (Homogeneous Equilibrium Model) de changement de phase proposé par Saurel et al. (1999) à trois équations est étudié et appliqué aux coups de bélier. Les résultats obtenus sont comparés aux données expérimentales. Malgré l’utilisation des techniques de renormalisation en SPH, des instabilités (oscillations numériques) se développent à l’interface entre les particules de matériaux différents. Ces instabilités restreignent l’utilisation des schémas SPH classiques pour des problèmes à faible ratio de densité. Afin de résoudre les problèmes de choc, le schéma proposé par Hu et Adams (2006) est adapté au régime fortement compressible en considérant le couplage entre la densité et la longueur de lissage. Les différents schémas SPH sont comparés entre eux pour les problèmes de chocs multiphasiques en 1-D et 2-D. Les résultats SPH sont validés avec la solution exacte pour les problèmes 1-D et la solution ALE pour les problèmes 2-D.

This thesis focuses on the numerical study of the propagation of shock waves in compressible multiphase flows and fluid structure interaction. Two approaches are being studied for the numerical solution of the fluid part: the ALE approach (Arbitrary Lagrangian Eulerian) and the Lagrangian SPH (Smoothed Particle Hydrodynamics) approach; while the structure part is solved by a conventional FE (Finite Element). The numerical investigation of the ALE and SPH methods are the two main areas of research.Water Hammers phenomena occuring in nuclear industries are investigated in this thesis. During a Water Hammer, the shock waves reflections in nuclear piping may drop locally the water pressure below its saturation pressure and generate cavitation. The three equations HEM (Homogeneous Equilibrium Model) phase change model proposed by Saurel et al. (1999) is studied and applied to solve water hammers. The obtained results are compared with experimental data. Despite the use of renormalization techniques in SPH, instabilities (numerical oscillations) are developed at the interface between particles from different materials. These instabilities restrict the use of traditional SPH schemes to problems with low density ratio. In order to solve the shock problems in the compressible regime, the scheme originally proposed by Hu and Adams (2006) is adapted to fully compressible regime (FC-SPH) by considering the coupling between the density and the smoothing length. The different SPH schemes are compared for 1-D and 2-D multiphase shock problems. Validation is performed in comparison with exact solutions for 1-D problems and ALE solution for 2-D problems.