Study of a high Reynolds number flow around a two dimensional airfoil at stall; an approach coupling a RANS framework and bifurcation theory

Denis Busquet

2020-06-24

Abstract:Le phenomene de decrochage est souvent decrit comme une chute soudaine de portance lorsque l'angle d'incidence augmente. Ce phenomene est prejudiciable aux avions et aux helicopteres et limite leur enveloppe de vol. Plusieurs etudes numeriques et experimentales, particulierement centrees sur le decrochage statique (i.e. pour des ailes fixes), ont revele des phenomenes apparaissant proche de l'angle de decrochage : des oscillations basses frequences et une hysteresis des coefficients aerodynamiques. Le premier phenomene se traduit par une oscillation de la portance entre une valeur maximale et une valeur minimale obtenues quand l'ecoulement est respectivement attache ou detache. Le nombre de Strouhal associe (St ~ 0.02) est habituellement un ordre de grandeur plus faible que le nombre de Strouhal (St ~ 0.2) du lâcher tourbillonnaire qui apparait pour de plus grandes incidences. Le second phenomene est caracterise par l'existence de solutions moyennees en temps autour de l'angle de decrochage qui different selon que l'angle d'attaque est augmente ou diminue.L'objectif de cette these est d'avoir une meilleure comprehension de l'origine du decrochage et de ces deux phenomenes grâce a des simulations numeriques d'ecoulements turbulents modelises par une approche RANS (Reynolds-Averaged Navier-Stokes). Une combinaison de diverses approches numeriques et theoriques (simulations instationnaires, continuation de solutions stationnaires, stabilite lineaire et analyse de bifurcation) est developpee et appliquee dans le cas du decrochage d'un profil 2D de pale d'helicoptere, le OA209, a bas nombre de Mach (M~0.2) et haut nombre de Reynolds (Re~1.8x10^6).Des solutions stationnaires sont calculees pour differents angles d'attaques en considerant le modele de turbulence de Spalart-Allmaras et en utilisant des methodes de continuation (continuation naive et methode du pseudo-arclength). Les resultats mettent en evidence une branche superieure (a haute portance), une branche inferieure (a basse portance) et, entre les deux, une branche du milieu. Pour un meme angle d'attaque, des solutions coexistent proche de l’angle de decrochage sur chacune des branches, ce qui est caracteristique d’un phenomene d'hysteresis. Des analyses de stabilite lineaire realisees autour de ces etats d'equilibres revelent l'existence d'un mode instable basse frequence associe au decrochage. L'evolution des valeurs propres associees a ce mode le long des branches stationnaires nous permet d'etablir une premiere version du diagramme de bifurcation. Afin de le completer, des calculs RANS instationnaires sont realises et des cycles limites basse frequence sont identifies sur une plage reduite d'angles d'attaque proches du decrochage. Ces solutions periodiques sont caracterisees par des valeurs de portance maximales et minimales plus grandes et plus petites que celles des solutions stationnaires a haute et basse portance associees, respectivement. Pour clarifier la formation et la disparition de ce cycle limite basse frequence et permettre une meilleure comprehension du scenario de bifurcation, un modele a une equation reproduisant les caracteristiques lineaires du phenomene est propose. Ce modele non-lineaire du decrochage statique est calibre sur les etats stationnaires et leur comportement lineaire obtenus par calculs RANS. Une etude du comportement non-lineaire de ce modele revele un scenario possible qui pourrait conduire a l'apparition et a la disparition du cycle limite basse frequence. Finalement, le cas d'un NACA0012 a nombre de Reynolds Re~1.0x10^6 est considere pour valider la robustesse du scenario identifie.

Physics