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

Chapitre 2: Modélisation à deux fluides de l’hydrodynamique et du transfert de masse dans les écoulements à bullesOn remarque que ce modèle ne dépend pas du nombre de Reynolds mais dépendseulement du diamètre de la bulle. D’autres modèles sont proposés dans la littérature entenant compte de la contamination des interfaces (Tomiyama, 2004).Force de masse ajoutéeL’expression de la force de masse ajoutée prend des formes différentes selon la définitiondes opérateurs de dérivée matérielle qui sont utilisés pour exprimer les effets inertielsengendrés par une bulle en accélération (Voinov, 1973 ; Auton et al., 1988). On peutconsidérer que cette question, qui a fait l’objet d’analyses et de discussions, (Gatignol,1983 ; Maxey & Riley, 1983 ; Minier, 1992), est dépassée grâce à la simulationnumérique. Rappelons, en effet, que les expériences numériques de Rivero et al. (1991)indiquent non seulement que le coefficient de masse ajoutée conserve sa valeur 0,5 pourdes bulles sphériques dans une gamme importante de nombres de Reynolds (de 0.1 à300), mais encore que les dérivées matérielles des deux phases interviennent dansl’expression de cette force. L’expression générale du terme moyen d’échange interfacialde quantité de mouvement dû à cette force s’écrit sous la forme suivante :< χ fGavec1MAi>= CM⎡⎛⎢⎜∂UαρL⎢⎣⎝∂tGi+ UGj∂U∂xGij⎞ ⎛⎟ ⎜∂U−⎠ ⎝∂tLi+ ULj∂U∂xCMle coefficient de masse ajoutée égal à 0.5 pour le cas d’une sphère.Lij⎞⎤⎟⎥⎠⎥⎦(2.49)Force de portanceLa force de portance est, par définition, la force orthogonale à la vitesse de glissement.Auton et al. (1988) ont calculé les forces qu’exerce un fluide parfait en écoulementfaiblement cisaillé rotationnel sur une sphère ; ils expriment la force de portance sous laforme :G1Li< χ f >= ε C αρ UijkLLRj ( rot( U L) k(2.50)avec εijkle symbole d’antisymétrie, CLreprésente le coefficient de portance égal à 0.5(Auton et al., 1988). Legendre (1996) et Legendre & Magnaudet (1998) ont calculé laforce de portance sur une bulle sphérique en écoulement visqueux linéairement cisaillépour une large gamme de nombres de Reynolds ( Re ∈ [10,500]). Ces calculs quis’appuient sur la simulation numérique, montrent que le coefficient de portance CLn’estpas constant mais dépend fortement du nombre de Reynolds, avec notamment pour lesgrands nombres de Reynolds:1 6.5CL = −(2.51)2 ReLance & Naciri (1992) ont mesuré la force de portance dans le cas d'une bulle placée dansun écoulement en rotation. Ces expériences montrent que la valeur de CLest de l’ordrede 0.25. En s’appuyant sur des analyses expérimentales et numériques, Tomiyama et al.20