THESE_EL HAMMAMI.pdf - Toubkal

3.3 Résolution numérique 683. À l’entrée, l’épaisseur du condensat est considérée nulle, on fixe dp ddzl’équation de quantité de mouvement pour calculer u.puis on résoud4. Intégrer numériquement l’équation de continuité pour trouver vv = 1 ρr∂∂z∫ r0ρurdr (3.44)5. Corriger la pression en utilisant la méthode de Newton [45], puis corriger les vitesses6. Résoudre les équations de l’énergie et de la diffusion7. Évaluer les nouvelles valeurs des variables dépendantes des propriétés thermophysiquedes fluides8. Justifier la satisfaction de la conservation de la masse dans l’écoulement gazeux etdans le condensatSi les critères suivants :∣ ∫ zπ(R − δ 0 z)ρ m v I dz − ∫ R−δ z∣πrρ0 m u m dr − m˙0 < ǫ (3.45)m˙0∣ Ṁ + ∫ z0 π(R − δ z)ρ m v I dz − ∫ RR−δ zπrρ l u l dr ∣ ∣Ṁ< ǫ (3.46)sont vérifiés, passer à la convergence des champs de la vitesse, de la température etde la fraction massique.Si l’erreur relative entre deux itérations successives pour u, T , et W satisfont auxcritèresmax ∣ φki,j − φ k+1 ∣i,jmax ∣ ∣ φki,j∣ ∣< ǫ φ : {u, T , et W } (3.47)La solution à chaque point de calcul est complète.Si non, répéter les étapes (2)-(7) jusqu’à ce que la condition de l’équation 3.47 soitjustifiée.