Développement d'une méthode de couplage partitionné fort en vue d'une application aux turbomachines

Development of a partitioned strong coupling procedure with the aim of turbomachinery application

• Auteur / Autrice : Guillaume Bénéfice
• Direction : Pascal Ferrand, Fabrice Thouverez
• Type : Thèse de doctorat
• Discipline(s) : Mécanique
• Date : Soutenance le 11/12/2015
• Etablissement(s) : Ecully, Ecole centrale de Lyon
• Ecole(s) doctorale(s) : Ecole doctorale Mécanique, Energétique, Génie Civil, Acoustique (Villeurbanne ; 2011-....)
• Partenaire(s) de recherche : Laboratoire : Laboratoire de mécanique des fluides et acoustique (Rhône) - Laboratoire de tribologie et dynamique des systèmes (Écully, Rhône ; 1970-)
• Jury : Président / Présidente : Lazhar Houas
• Jury : Examinateurs / Examinatrices : Yannick Rozenberg, Eric Lippinois, Pierre Argoul
• Jury : Rapporteurs / Rapporteuses : Lazhar Houas, Dominique Pelletier

Mots clés

Mots clés contrôlés

Turbomachines -- Compresseurs

Interaction fluide-structure

Mots clés libres

Interaction fluide-structure

Couplage partitionné fort

Schéma temporel implicite

Compresseur

Ondes de choc

Modélisation de la turbulence

Fluid-structure interaction

Partitionned strong coupling

Implicit temporal scheme

Compressor

Shockwaves

Turbulence model

Résumé

Pour améliorer la conception des turbomachines, les industriels doivent appréhender des phénomènes aéroélastiques complexes présents dans les compresseurs comme les cycles limites d’interaction fluide-structure des fans. La compréhension et la modélisation de ces phénomènes impliquent de développer des modèles numériques complexes intégrant des phénomènes multi-physique et de valider ces modèles à l’aide de bancs d’essais. Le banc d’essai du compresseur CREATE est instrumenté pour étudier des instabilités aérodynamiques couplées à des vibrations, notamment sur le rotor du premier étage, et permet de valider des modèles numériques. La modélisation de l’écoulement en amont du premier étage du compresseur à l’aide du logiciel Turb’Flow, développé pour l’étude des écoulements dans les compresseurs aéronautiques, a permis de mettre en évidence l’importance des conditions limites d’entrée pour l’obtention de résultats précis. En particulier, il a été possible de modéliser correctement l’ingestion d’une alimentation non-homogène en entrée de la roue directrice d’entrée. Ce phénomène peut se produire en amont des fans et interagir avec un mode de la structure. Une stratégie de couplage partitionné fort explicite dans le domaine temporel a été introduite dans le logiciel Turb’Flow. Comme cette méthode présente un risque de décalage temporel à l’interface fluide-structure, une attention particulière a été portée à la modélisation de la conservation de l’énergie à cette interface. La conservation de l’énergie à l’interface est cruciale quand les déplacements sont importants et quand un comportement non-linéaire fort apparaît entre le fluide et la structure (onde de choc et amortissement structurel nonlinéaire). Parallèlement au développement du module aéroélastique, le schéma implicite de Runge- Kutta d’ordre 3 en temps (RKI-3) a été développé et évalué sur un cas de dynamique (vibration d’une aube de turbine transsonique) et sur un cas de propagation d’onde de choc. L’utilisation du schéma RKI-3 permet d’augmenter, à iso-précision, d’un ordre le pas de temps par rapport aux schémas de Gear et de Newmark. S’il apporte un gain en temps CPU pour l’étude de la dynamique des structures, il est pénalisant dans le cadre de simulation URANS. Cependant, le schéma RKI-3 est utilisable dans le cadre de simulations couplées fluide-structure.

To increase turbomachinery design, manufacturers have to comprehend complex aeroelastic phenomena involving compressors like fluid-structure interaction limit cycles of fans. The understanding and the modeling of these phenomena involve developing complex solvers coupling techniques and validating these techniques with bench tests. The bench test of the CREATE compressor is instrumented to study the coupling between aerodynamic instabilities and structure vibration, in particular on the first stage rotor, and allows to validate numerical techniques. The flow modeling upstream to the first stage with the Turb’Flow flow solver (targeting turbomachinery applications) shows that, to have accurate results, inlet limit conditions must take into account. The ingestion of non-homogeneous flow upstream to the inlet guide vane is accurately modeled. This phenomenon can appear upstream to fans and interact with structure Eigen-modes. Explicit partitioned strong coupling considered in time domain was implemented in a Turb’Flow flow solver. As there is a risk of time shift at the fluid-structure interface, careful attention should be paid to energy conservation at the interface. This conservation is crucial when displacements are large and when strong non-linear behaviors occur in both fluid and structure domains, namely shock waves, flow separations and non-linear structural damping. In parallel with coupling technique development, the three-order implicit Runge-Kutta scheme (RKI-3) was implemented and validated on a structure dynamic case (transonic turbine blade vibration) and on a case of shock waves propagation. The RKI-3 scheme allows increasing the time step of one order of magnitude with the same accuracy. There is a CPU time gain for structure dynamics simulations, but no for URANS simulations. However, the RKI-3 scheme can be to use for fluid-structure coupling simulations. The coupling technique was validated on a test case involving tube in which the shock wave impinges on a cross flow flexible panel, initially at rest. This case allows modeling an interaction between sonic flow and a panel movement with a tip clearance. Some numerical simulations were carried out with different temporal schemes. The RKI-3 scheme has no influence on results (compared with Gear and/or Newmark scheme) on the energy conservation at the fluid-structure interface. Compared to experimental results, pressure is in fairly good ix Liste des publications agreement. The analysis of numerical results highlighted that a vertical shock tube with up and down waves creates pressure fluctuation. Frequency is under predicted and amplitude is not in fairly good agreement. The panel root modeling might be questionable.