13.07.2015 Views

Analyse expérimentale et modélisation du transfert de matière et du ...

Analyse expérimentale et modélisation du transfert de matière et du ...

Analyse expérimentale et modélisation du transfert de matière et du ...

SHOW MORE
SHOW LESS
  • No tags were found...

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Chapitre 2: Modélisation à <strong>de</strong>ux flui<strong>de</strong>s <strong>de</strong> l’hydrodynamique <strong>et</strong> <strong>du</strong> <strong>transfert</strong> <strong>de</strong> masse dans les écoulements à bullespseudo-turbulente basé sur la solution potentielle qui génère <strong>de</strong> l'énergie cinétique in<strong>du</strong>itepar les bulles par les seuls eff<strong>et</strong>s <strong>de</strong> la force <strong>de</strong> masse ajoutée est incompl<strong>et</strong> pour lesécoulement à taux <strong>de</strong> vi<strong>de</strong> relativement élevés.Chahed <strong>et</strong> al. (2004) ont ainsi proposé <strong>de</strong> prendre en compte ces eff<strong>et</strong>s <strong>de</strong> sillage enconsidérant que l’énergie cinétique totale in<strong>du</strong>ite par les bulles k BIest la somme <strong>de</strong>l'énergie pro<strong>du</strong>ite dans le sillage k WT<strong>et</strong> <strong>de</strong> la pseudo-turbulence engendrée par les eff<strong>et</strong>s<strong>de</strong> masse ajoutée kS:k = k + k(2.115)BISWToù kSest l’énergie cinétique pro<strong>du</strong>ite par la masse ajoutée dans le mouvement relatifmoyen. En écoulement homogène kSest donnée par la solution <strong>de</strong> l’écoulement potentiel<strong>de</strong> Biesheuvel & Wijngaar<strong>de</strong>n (1984) :1 2kS= αU R(2.116)4En extrapolant la solution auto-similaire <strong>de</strong> l’énergie cinétique turbulente dans le sillagelointain d’une bulle isolée à celui d’une bulle évoluant en milieu confiné, Chahed <strong>et</strong> al.(2004) expriment ainsi l’énergie cinétique turbulente intégrée dans l’espace inter-bullespar :kWT=1−3 1ηα3( + 1) −14 3.7523 2 4mαπ χ mU3R0RETη (1 − α)(2.117)13⎡πRappelons que6 ⎥ ⎤η = ⎢, que γ est le coefficient d'aplatissement défini par le rapport <strong>de</strong>⎣ ⎦la longueur <strong>du</strong> p<strong>et</strong>it axe sur celle <strong>du</strong> grand axe, que R ETest le nombre <strong>de</strong> Reynoldsturbulent, UR0<strong>et</strong> URsont respectivement la vitesse limite d’une bulle en ascension libre<strong>et</strong> la vitesse relative dans l’écoulement diphasique à fort taux <strong>de</strong> vi<strong>de</strong> donnée par laformulation (2.61).C<strong>et</strong>te formulation a été confrontée aux données expérimentales <strong>de</strong> Garnier <strong>et</strong> al. (2002),la comparaison indique une excellente concordance jusqu’à <strong>de</strong>s taux <strong>de</strong> vi<strong>de</strong> <strong>de</strong> l’ordre <strong>de</strong>15%. Rappelons qu’au <strong>de</strong>là <strong>de</strong> c<strong>et</strong>te valeur, la distance inter-bulle <strong>de</strong>vient trop faible pourque la représentation même <strong>du</strong> sillage confiné ait un sens.Il est intéressant <strong>de</strong> remarquer, que l’énergie cinétique est normée par0 RU . C’est unrésultat qui est en accord avec les expériences <strong>de</strong> Larue <strong>de</strong> Tournemine (2001) quiindiquent que c’est en eff<strong>et</strong> la vitesse <strong>du</strong> champ lointain qui normalise bien l’énergiecinétique turbulente in<strong>du</strong>ite par les bulles.47

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!