Paper 3

Etude du comportement aérodynamique d'une pale flexible
par interaction fluide structure

FAKHFEKH Marwaa, BEN MIRA Waêl a,b, ABID Malek c, MAALEJ Arefa
a
Laboratoire des systèmes éléctromécaniques LASEM, Ecole nationale d’ingénieurs de Sfax, Route
soukra km 4 Sfax, Tunisie, fakhfekh.marwa@enis.tn
b
Département de génie mécanique, Institut supérieur des systèmes industriels de Gabès -
Université de Gabès, Gabes, 6029, Tunisie, benamira.wael@gmail.com
c
Institut de recherche sur les phénomènes hors équilibre IRPHE, Université Aix-Marseille, Marseille,
France, E-mail abid@irphe.univ-mrs.fr
Résumé – Dans ce papier, nous présentons une investigation numérique basée sur l’interaction fluide
structure FSI pour l’étude de l’effet de l’élasticité sur les performances aérodynamiques d’un rotor
à pales flexibles. Les simulations sont développées sous le logiciel Ansys Workbench avec un
couplage fort entre la dynamique des fluides et la structure transitoire. Comme résultats, nous
discutons l’effet de l’élasticité sur la variation de la force de trainée en fonction de la vitesse de
l’écoulement et la vitesse de rotation du rotor. Pour une pale rigide, la force de traînée varie en
fonction du carré de la vitesse mais pour une pale flexible cet exposant est modéré par un coefficient
υ appelé nombre de Vogel.
Mots clés : Force de trainée/ Pale flexible/ ANSYS Fluent / FSI
Abstract – In this paper, we present a numerical investigation based on the fluid-structure interaction
FSI to study the effect of elasticity on the aerodynamic performance of a flexible blade rotor. The
simulations are carried out under Ansys Workbench software with a 2 way strong coupling between
fluid dynamics and transient structure. As results, we discuss the impact of the elasticity on the
variation of the drag force as a function of the flow velocity and the rotor rotation frequency. In
literature, drag force is a function of the body’s geometry, velocity of the fluid, and the fluid
properties. For rigid blade, drag force varies according to the square of the speed but for flexible
body this exponent is moderated by a coefficient υ called Vogel number
Key words: Drag force / Flexible blade/ ANSYS Fluent / FSI
Monastir 17-19 mars 2023, Tunisie

6e Congrès Tunisien de Mécanique COTUME’2023
1 Introduction
  2sε  tr (ε).I (4)
Au cours des dernières décennies et pour répondre à
Avec:
la demande croissante d'énergie, les recherches sont
focalisées sur le développement de grandes
éoliennes fiables et efficaces. E E.
 
La conception des éoliennes de très grande tailles 2(1   ) (1   )(1  2 )
E est le module de young,  est le champs de
pousse les chercheurs à concevoir des matériaux de
contrainte,  est la déformation du solide, I est le

plus en plus performants. Pour cela, la plupart des
pales d'éoliennes sont fabriquées à partir des
matériaux composites qui combinent la légèreté et la tenseur d’identité et  est le coefficient de poisson.
résistance. [1] Les équations de couplage sont modélisées par une
Le développement de taille des rotors, et l’utilisation condition cinématique et une condition dynamique :
des matériaux légers font que les pales sont soumises u s
à des déformations qui sont parfois importantes. Ces u f ( x, t )  (5)
déformations perturbent l’écoulement autour de ces t
rotors et affectent leurs performances. Pour cette n. int
solid
erface  n. intfluid
erface (6)
raison, il est intéressant d’étudier les problèmes
d’aéroélasticité. Dans ce cadre, différentes études n est le vecteur unitaire
sont présentées afin de caractériser l’effet de
l’élasticité sur le comportement aérodynamique des Pour tenir compte de la turbulence de l’écoulement,
corps en écoulement. On peut citer les travaux de les équations de Navier Stokes sont complétés par
Kang et al [2] et Dai et al [3] qui ont étudié l'effet de l’un des modèles de turbulence (le modèle K dans
la flexibilité structurelle sur les performances des notre cas)
ailes battantes et Gosselin et al. [4,5] qui ont examiné
l'effet de la flexibilité et de la reconfiguration sur la 3 Simulations numériques :
réduction de la traînée. Toutes fois, ces études
s’intéressent à l’effet de la vitesse d’écoulement sur La simulation numérique est réalisée sous le logiciel
les performances aérodynamiques. Nous nous Ansys Workbench en faisant appel à trois modules :
intéressons dans ce travail à l’étude de l’effet de la un pour la modélisation numérique de la dynamique
vitesse de rotation sur la force de trainée d’un rotor des fluides, un deuxième pour la résolution du la
à pales flexibles. structure transitoire, et un troisième pour assurer le
Pour cette étude, nous développons une approche couplage entre ces deux derniers.
numérique basée sur l’interaction fluide-structure
pour caractériser l'effet de l'élasticité sur les Le domaine fluide est constitué de deux sous
propriétés aérodynamiques d’un rotor à pales domaines : un domaine rectangulaire stationnaire qui
flexible dans l’eau. représente le canal d'écoulement, et un domaine
cylindrique rotatif pour caractériser la rotation du
2 Modélisation rotor. Les dimensions du domaine de calcul (figure
1) sont choisies pour des exigences numériques et
Les problèmes d'interaction fluide-structure pour pour la convergence de la solution.
une structure transitoire dans un écoulement
turbulent sont modélisés par les équations
incompressibles de Navier Stokes, et les équations
linéaires de l'élasticité et qui sont modélisées par les
systèmes suivants :
u f (1)
f   f (u f .)u f  P   f .(u f )  Fv
t
u f  0 (2)
 us
2
s .  .  Fs (3) Figure 1. Dimensions du domaine de calcul
t 2
Avec:uf est le champs du fluide incompressible, us Le rotor étudié est de rayon R=88mm. Il est constitué
est le champs de déplacement de la pale, t est le d'une seule pale, d’un moyeu et d’un arbre. L’arbre
temps, μf est la viscosité dynamique, Fv est la force est un axe de diamètre 15 mm en nylon. Le moyeu
du corps sur le fluide, et Fs est la force par unité de est conçu en fibre de carbone. La pale est une lame
  rectangulaire avec une corde constante c=20mm et
volume sur la lame. f et s sont la densité du
une épaisseur uniforme de 0,7 mm. Elle est flexible
fluide et la densité de la structure respectivement.
en polyéthylène basse densité (LDPE) ayant un
Monastir 17-19 mars 2023, Tunisie 1

6e Congrès Tunisien de Mécanique COTUME’2023
module de Young de 3,5 GPa, une densité de 1070 rotor pour des pales rigides et flexibles et pour deux
kg/m3 et un module de cisaillement de 2.1 Gpa. Les configurations de vitesse d’entrée U0 = 18 cm/s et
dimensions du rotor sont résumées dans la Figure 2. U0=0.36m/s.
Cette évolution montre bien l’effet la flexibilité sur
les performances aérodynamiques. En effet, la
flexibilité diminue d’avantage la force de trainée et
l’évolution de cette dernière n’est plus monotone
dans le cas flexible. Ce changement de loi est
strictement lié au comportement élastique de la pale.
Figure 2. Dimensions du rotor
Le couplage entre le domaine fluide et le domaine

structure transitoire est effectué par des simulations
couplées avec un transfert de données bidirectionnel.
Le couplage est assuré par le maillage dynamique au
niveau de l’interface fluide- structure.
Dans cette étude, le rotor est immergé dans l'eau où

les pressions du fluide sont plus importantes que
dans l'air et peuvent générer des déformations Figure 4. Evolution de la force de trainée en
beaucoup plus importantes donc un effet plus visible. fonction de la fréquence de rotation
4 Résultats et discussion : 5 Conclusion :
Notre objectif est de caractériser l’effet de l’élasticité Nous avons étudié le comportement aérodynamique
sur les performances aérodynamiques. Donc comme d'un rotor à pales flexibles. Le couplage entre la
résultats, on représente dans un premier lieu structure transitoire et la dynamique des fluides,
l’évolution de la force de trainée en fonction de la nous a permis d'étudier l'effet de l’élasticité sur les
vitesse linéaire d’entrée de l’écoulement pour deux performances aérodynamiques des pales.
cas de pale : flexible et rigide. La figure 3 montre Les paramètres d’écoulement génèrent des
bien les lois de puissance pour l'évolution de la force déformations dans les pales, pouvant être
de traînée en fonction de la vitesse d'écoulement. importantes et affectent les performances des rotors.
Cette loi ayant un exposant 2 pour le cas rigide et un
exposant modéré par un coefficient υ appelé nombre 6 Bibliographie
de Vogel (υ -0.3 dans ce cas) [4].
[1] J Manwell, J. F., McGowan, J. G. & Rogers, A.
1 (2002), Wind Energy Explained: Theory,
F  CD AU 2 Design and Application,
2 (7) [2] C.-K. Kang, H.Aono, C. E. S. Cesnik and W.
Shyy, Effects of flexibility on the aerodynamic
2 performance of flapping wings. J. Fluid Mech.
F=1,09U2
Rigide
689 (2011), 32–74.
1,5 Flexible
F=0,88U1,7
[3] H. Dai, H. Luo, J. F. Doyle, Dynamic pitching
F(N)
1 of an elastic rectangular wing in hovering

motion, Journal of Fluid Mechanic 693
0,5 (2012)473{499.
[4] F.Gosselin, E.Delangre, B.A.Machado-
0
0 0,5
Almeida, Drag reduction of flexible plates by
U (m/s) 1 1,5
reconfiguration, J.Fluid Mech.650, (2010) 319.
[5] E. De Langre, A. Gutierrez and J. Cossé, On the
Figure 3. Evolution de la force de trainée en
fonction de la vitesse d’entrée scaling of drag reduction by reconfiguration in
plants. C. R. Mecanique 340 (2012) 35–40.
La figure 4 représente l’évolution de la force de
trainée en fonction de la fréquence de rotation du
Monastir 17-19 mars 2023, Tunisie 2

Paper 3

Transféré par

Droits d'auteur :

Formats disponibles

Paper 3

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Paper 3

Transféré par

Droits d'auteur :

Formats disponibles

Etude du comportement aérodynamique d'une pale flexible

par interaction fluide structure

Mots clés : Force de trainée/ Pale flexible/ ANSYS Fluent / FSI

Key words: Drag force / Flexible blade/ ANSYS Fluent / FSI

Monastir 17-19 mars 2023, Tunisie

contrainte,  est la déformation du solide, I est le

Monastir 17-19 mars 2023, Tunisie 1

Figure 2. Dimensions du rotor

Le couplage entre le domaine fluide et le domaine

Dans cette étude, le rotor est immergé dans l'eau où

4 Résultats et discussion : 5 Conclusion :

1 of an elastic rectangular wing in hovering

Monastir 17-19 mars 2023, Tunisie 2

Vous aimerez peut-être aussi