Modélisation de l'évaporation de gouttes multi-composants

CHAPITRE 4 MODELISATION NUMERIQUEEtude du profil de compositionLes conditions sont toujours les mêmes que précédemment avec Δt=100µs, NT=20 et lesvaleurs utilisées pour NC sont 60, 100, 150 et 300. Les résultats obtenus sont présentés sur laFigure 4.11.10080250200Dg (µ m)Ydecane,l s60402000.50.40.30.20.10diff. limitée NC=60diff. limitée NC=100diff. limitée NC=150diff. limitée NC=3000 0.01 0.02 0.03 0.04 0.05t (s)0 0.01 0.02 0.03 0.04 0.05t (s)Tl s (°C)dm/dt_total (µ g/s)150100141210864205000 0.01 0.02 0.03 0.04 0.05t (s)0 0.01 0.02 0.03 0.04 0.05t (s)dm/dt_decane (µ g/s)8765432100 0.01 0.02 0.03 0.04 0.05t (s)dm/dt_hexadecane (µ g/s)98765432100 0.01 0.02 0.03 0.04 0.05t (s)Figure 4.11 : Comparaison des calculs pour différentes valeurs de NCLe diamètre, la température à la surface de la goutte et le débit total de vapeur ne sont pasaffectés par la valeur de NC, ce qui n'est pas le cas des fractions massiques des espècesliquides à la surface de la goutte et donc les débits massiques de vapeur pour chaquecomposant.Il apparaît clairement que la faible diffusion des espèces liquides contraint à choisir unevaleur de NC plus importante que celle de NT. D'ailleurs, le calcul ne converge pas avecNC=50. D'après ces résultats, 60 nœuds de discrétisation semblent suffisants pour obtenir unebonne précision, mais le calcul est à la limite de la stabilité. Des calculs effectués pour137