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

Chapitre 4 : Simulation de l’hydrodynamique et du transfert de masse dans les écoulements de types zone de mélange diphasique4.5.5 Sensibilité aux divers termes du transfert interfacial dequantité de mouvementOn se propose dans cette partie de tester les effets des différents termes de la formulationdu terme d’échange interfacial de quantité de mouvement (éq. 4.6) sur la distribution desphases. Nous analysons dans cette étude numérique l’importance des termes non linéairesissus de la masse ajoutée et du lift dans le phénomène de distribution des phases. Dansune première série de simulations de la zone de mélange à bulles à faible taux de vide,nous analysons l’effet des termes non linéaires issus de la masse ajoutée, ensuite dans uneseconde série de simulations, nous étudions l’effet du coefficient de portance.4.5.5.1 Sensibilité de la distribution des phases aux termes nonlinéaires de masse ajoutéeNous testons dans cette section la sensibilité du modèle d’échange interfacial aux termesnon-linéaires issus de la masse ajoutée. Nous avons remarqué, sur la figure 3.19, que lerapport des variances dans les deux phases n’est pas constant dans la zone des gradients etvarie entre les valeurs 1 à 2 et 4 à 5. À l’heure actuelle, on n’a pas de modèle précis pourcalculer l’agitation des bulles mais nous allons nous appuyer sur ces mesuresexpérimentales pour tester l’effet du terme non linéaire de la masse ajoutée avecdifférentes valeurs de C ;C 11 22allant de 2 à 9 parce que la valeur 9 correspond à la valeurmaximale obtenue par la théorie de la dispersion de Tchen.Nous avons réalisé une série de simulations décrites ci après.Dans la simulation ZDM1100, on ne tient compte, dans le terme d’échange interfacial dequantité de mouvement, que des termes moyens de traînée. Dans les simulations(ZDM1211 ; ZDM1311 et ZDM1911), nous supposons que les tenseurs des contraintesturbulentes dans le gaz sont proportionnelles à celles du liquide en utilisant les relations(4.13) avec respectivement les valeurs 2, 3 et 9 pour les coefficients des termes nonlinéaires ( C ;C 11 22). On prend la valeur C = 121 pour le coefficient des composantes nondiagonales. Nous introduisons aussi dans ces trois simulations un effet de la forcemoyenne de portance et un effet de la force moyenne de masse ajoutée. La force deportance est exprimée avec un coefficient CL= 0. 25 conformément aux expériences deLance & Naciri (1992) pour une bulle dans un écoulement en rotation. Et la force demasse ajoutée est exprimée avec le coefficient CM= 0. 5 . Le tableau 4.1 récapitule lesconditions de calcul dans ces simulations. Les profils de taux de vide sont reportés sur lesfigures (4.24) à (4.27).123