جامعة وهران للعلوم و التكنولوجيا محمد بوضياف
République Algérienne Démocratique et Populaire
Ministère de l’Enseignement Supérieur et de la Recherche
Scientifique
Université des Sciences et de la Technologie d’Oran « Mohamed BOUDIAF »
Faculté de Chimie Département de Génie Chimique
Domaine : Sciences et Technologie
Spécialité : Génie chimique
Thème
Teste 2 MPD
Présenté par :
RAMDANI SIDAHMED
HARRAG Kada
Groupe : 01
Exercice TD :
Une colonne garnie est utilisée pour la distillation sous vide d'un mélange d'isomères
de masse moléculaire 155kg/kmole
La T° moyenne régnant dans la colonne est de 373 °K et les pressions en bas et en haut de la
colonne sont respectivement de 1,3KN/m² et 0,13 kN/m².
En admettant que l'écoulement a contre-courant n'affecte pas la chute de pression du gaz et
que cette dernière (ΔP) peut être calculée par l'équation de Carman.
Montrer que dans ce domaine de pression de fonctionnement la perte de
charge est directement proportionnelle au débit massique de Vapeur.
Calculer approximativement le débit massique
de Vapeur.
Données :
Z= 9m, Diamètre = 1,2m, ε= 0.71
ag= 656m2/m3, Viscosité =0,018 10-3Ns/m2
Volume molaire de vapeur 22,4 m3/Kmole.
Solution de l’EXO TD par Matlab :
1-on calcule Reynolds :
R’ep= Qm
η ag (1- ε)
R’ep= 292,03Qm > 0,2Régime Turbulant
Equation de Carman :
f’p= 5Rep-1+0,4Rep-0,1
f’p= 5(292,03Qm) -1+0,4(292,03Qm) -0,1
f’p=0,017Qm-1 +0,227Qm-0,1
Calculer ρ :
ρ=PM = P 155 .
RT 0,082 373 101325
ρ=5 10-5 p
f’p= ΔP . ε3 . = ΔP . ε3 ρ . = ρ Um= Qm
Z ag(1- ε) ρ Um2 Z ag(1- ε) Qm2
f’p= ΔP P 5 10-5(0.71)3
Qm2 9 656 (1-0.71)
0,017Qm-1+0,227Qm-0,1=P ΔP . 1,045 10-8
Qm2
PΔP =1,627 .106 Qm+21,72. 106 Qm
On fait l’intégrale PΔP
On trouve :
P12 – P22=1,627 .106 Qm+21,72. 106 Qm2
2
0,83655=1,627 .106 Qm+21,72. 106 Qm2
On résoudre cette équation par logiciel Matlab comme suivants :
Le code utilisée :
dp = 0.83655 \% Dd́éfinir la fonction à résoudre
equation = @(qm) 1.627e6 * qm + 21.72e6 * qm.^1.9 - dp;
qm_initial_guess = 1;
qm_solution = fsolve(equation, qm_initial_guess); disp([’La valeur de
qm est : ’, num2str(qm_solution)]);
La valeur de Qm est : 5.1415e-0
Exercice similaire :
Dans le cadre de l'exploitation d'une colonne de distillation sous vide utilisée pour séparer un
mélange de composés organiques de masse moléculaire égale à 180 kg/kmol, on maintient une
température moyenne de 350K. Les pressions mesurées à la base et au sommet de la colonne
sont, respectivement, de 50 kN/m² et de 0,15 kN/m².
On néglige l'impact de l'écoulement en contre-courant du liquide sur la chute de pression du
gaz, et on suppose que cette dernière peut être déterminée à l'aide de l'équation de carman.
Établir la relation entre la chute de pression dans la colonne et le flux massique des
vapeurs, en justifiant pourquoi elle est supposée linéaire dans la gamme de pressions
opérationnelles données.
Dériver une expression permettant de calculer le débit massique des vapeurs en se
basant sur les données opérationnelles et vérifier cette expression avec les valeurs
fournies.
Données : Hauteur de la colonne Z = 10m Diamètre interne Φ = 0.15m Porosité ε = 0.80 Aire
spécifique de la garniture ag = 700 m²/m³ Viscosité du gaz μGaz = 0,012 x 10^-3 Ns/m²
Volume molaire du gaz à l'état standard 22.4 m³/kmol.
Solution de l’exercice similaire :
On résoudre cette équation par logiciel Matlab comme suivants :
dp = 1249.98
% D´efinir la fonction `a r´esoudre
equation = @(qm) 3.587e6 * qm + 94.61883e6 * qm.^1.9 - dp;
% Supposons une valeur initiale de qm qm_initial_guess = 1;
% Utiliser fsolve pour trouver la valeur de qm qm_solution = fsolve(equation,
qm_initial_guess);
% Afficher la valeur de qm
disp([’La valeur de qm est : ’, num2str(qm_solution)]);
La valeur de Qm est : 0.00034164
