TP – Initiation sur Matlab
Le procédé de conversion de soufre en acide sulfurique se réalise selon 3 étapes chimiques :
La combustion de soufre liquide dans un four de combustion :
S+O 2 → SO 2
L’oxydation de dioxyde de sulfure :
1
SO 2+( )O2 → SO 3
2
Et la réaction entre l’eau et le SO3 (absorption chimique) :
H 2O+ SO 3 → H 2 SO 4
Dans le procédé par contact, la conversion de SO 2 en SO3 se réalise dans un réacteur à quatre lits
catalytiques, utilisant comme catalyseur l’oxyde de vanadium. La cinétique de conversion est
donnée par le modèle de Foment et Bischoff :
( )
P SO
−( K 1 PO P SO 1− 3
0.5
)
2 2
K p P SO PO
r SO =
2 2
2
2
(22.414 (1+ K 2 P SO + K 3 PSO ) )
2 3
Avec:
K p =exp (
11300
T
−10.68); K 1=2.1125 105 exp
T (
−3599.78
)
; K 2=14.641 et K 3=6.5775.
La vitesse de réaction est exprimée en (kmol/kg de cat./sec.).
La température d’activation du catalyseur est de 400°C (Température minimale d’entrée du
premier réacteur). La composition du mélange gazeux à la sortie du four de combustion est
donnée dans le tableau suivant:
Constituant Fraction molaire (%)
O2 9
SO2 12
N2 79
Débit molaire total 2.15 kmol/sec
Pression des Gaz 1.5 bars
Les spécification et propriétés thermodynamiques des différents constituants mis en jeu sont
données comme suit :
% Données de Physico-chimiques
Molar_H2SO4=98; % Molar mass (kg/kmol)
Molar_SO2=64; % Molar mass (kg/kmol)
Molar_SO3=80; % Molar mass (g/mol)
Molar_O2=32; % Molar mass (g/mol)
Molar_N2=28; % Molar mass (g/mol)
Molar_H2O=18; % Molar mass (g/mol)
Molar_S=32; % Molar mass (g/mol)
Df_H2SO4=-811.3; % Enthalpie de formation standard (kJ/mol)
Df_SO3=-395.2; % Enthalpie de formation standard (kJ/mol)
Df_SO2=-296.9; % Enthalpie de formation standard (kJ/mol)
Df_H2O=-285.8; % Enthalpie de formation standard (kJ/mol)
%capacités thermiques massiques des constituants
Cp_H2SO4=1.34; % kJ/K/kg
Cp_H2O=4.18; % kJ/K/kg
Cp_SO2=0.64; % kJ/K/kg
Cp_SO3=0.63; % kJ/K/kg
Cp_S_liq=0.72; % kJ/K/kg
Cp_S_solid=0.71; % kJ/K/kg
Cp_N2=1.04; % kJ/K/kg
Cp_O2=0.919; % kJ/K/kg
Cp_air=0.79*Cp_N2+0.21*Cp_O2;
%capacités thermiques molaires des constituants
Cpm_H2SO4=Cp_H2SO4*Molar_H2SO4; % kJ/K/kmol
Cpm_H2O=Cp_H2O*Molar_H2O; % kJ/K/kmol
Cpm_SO2=Cp_SO2*Molar_SO2; % kJ/K/kmol
Cpm_SO3=Cp_SO3*Molar_SO3; % kJ/K/kmol
Cpm_S_liq=Cp_S_liq*Molar_S; % kJ/K/kmol
Cpm_S_solid=Cp_S_solid*Molar_S; % kJ/K/kmol
Cpm_N2=Cp_N2*Molar_N2; % kJ/K/kmol
Cpm_O2=Cp_O2*Molar_O2; % kJ/K/kmol
1) Créer sur Matlab les deux fonctions de bilan de matière et thermique caractérisant la
conversion de SO2 en SO3. (Pour un volume élémentaire dV sous forme de disque).
2) Donner la procédure d’intégration des deux équations de bilan.
3) En intégrant les équations de bilan, déterminer la quantité de catalyseur pour atteindre la
conversion maximale de SO2 au niveau de chaque lit catalytique. On suppose que T1=400°C,
T2=405°C, T3=410°C et T4=410°C.
4) Représenter dans un graphique l’effet de la température des gaz à l’entrée de chaque lit
catalytique sur le taux de conversion et aussi sur la puissance thermique récupérée. On suppose
que la quantité du catalyseur est suffisante pour atteindre l’équilibre au niveau de chaque lit
catalytique.
Algorithme de calcul
