THESE_EL HAMMAMI.pdf - Toubkal

3.2 Discrétisation 60∆r +∆r −= α (3.11)Les pas de calcul dans la direction axiale et dans la direction radiale (dans les deuxphases) sont donnés par des séries géométrique comme suit :∆z i = β∆z i−1 i = {1, ........ni} (3.12)∆r mj = α m ∆r mj−1 j = {1, ........nj} (3.13)∆r lj = α l ∆r lj−1 j = {nj + 1, ........nk − 1} (3.14)Avec nk et nj représentent respectivement le nombre de noeuds dans le sens radial etle nombre de noeuds dans la phase gazeuse. Le nombre de noeuds dans la phase liquideet dans la direction axiale sont données respectivement par nl = nk − nj et niChaque noeuds du maillage est repéré par rapport à l’axe des z et des r dans la phaseliquide et la phase gazeuse respectivement par les indices i et j. Les coordonnées z et rd’un noeud (i, j) sont données par :z i+1 − z i = β(z i − z i−1 ) = β 2 (z i−1 − z i−2 ) = ... = β i−1 (z 2 − z 1 )= β i−1 ∆z 1 (3.15)z i − z i−1 = β i−2 ∆z 1z i−1 − z i−2 = β i−3 ∆z 1z 2 − z 1 = β∆z 1Donc on obtient :z i = ∆z 1β i−1 − 1β − 1De la même façon on obtient :dans la direction axiale i = {1, ........ni} (3.16)r j = ∆r m1α j−1m − 1α g − 1dans la phase gazeuse j = {1, ........nj} (3.17)