Bilan de matière et thermique dans un réacteur

Introduction
La réaction dans le réacteur est : DMA + FM → DMF + MOL
Nous allons calculer le débit molaire du DMF.

Données des masses molaires

• Masse molaire de DMA : MDMA = 45, 08 g/mol = 0, 04508 kg/mol

• Masse molaire de FM : MFM = 60, 05 g/mol = 0, 06005 kg/mol

• Masse molaire de DMF : MDMF = 73, 09 g/mol = 0, 07309 kg/mol

• Masse molaire de MOL : MMOL = 32, 04 g/mol = 0, 03204 kg/mol

Calcul du débit molaire du DMF

Nous avons les données suivantes :

• Débit massique du DMF : ṁDMF = 1250 kg/h

Le débit molaire (ṅDMF ) peut être calculé à l’aide de la formule :

ṁDMF
ṅDMF =
MDMF
A.N. :
1250
ṅDMF = = 17171, 2 mol/h
0, 07309
Le débit molaire du DMF est donc de 17 171,2 mol/h.

Tableau d’avancement pour une réaction totale

Quantité DMA FM DMF MOL

Initial (mol) ṅ0 ṅ0 0 0
Intermédiaire (mol) ṅ0 − ẋ ṅ0 − ẋ ẋ ẋ
Final (mol) 0 0 ṅ0 ṅ0
D’après le tableau, les débits molaires des réactifs DMA et FM sont égaux, et ceux
des produits DMF et MOL le sont également :

ṅ(DMA) = ṅ(FM) = ṅ(DMF) = ṅ(MOL) = 17171, 2mol/h

1
Tableau récapitulatif des débits massiques

Élément Entrée (Débit massique) Sortie (Débit massique)

DMA 774, 63 kg/h 0
FM 1034, 99 kg/h 0
DMF 0 1250 kg/h
MOL 0 549, 36 kg/h
Total 1809, 62 kg/h 1799, 36 kg/h

Calculs des débits massiques réels

Pour calculer les débits massiques réels, nous appliquons les formules suivantes :
Formule pour les réactifs :
ṁ
ṁréel =
ηchimique × ηphysique

Formule pour les produits :

ṁréel = ṁ × ηphysique

Valeurs des coefficients η :

- DMA : ηchimique = 0, 976, ηphysique = 0, 995
- FM : ηchimique = 0, 958, ηphysique = 0, 995
- DMF : ηphysique = 0, 995
- MOL : ηphysique = 0, 984

Tableau des résultats des débits massiques réels

Élément Entrée (Débit massique réel) Sortie (Débit massique réel)

DMA 799, 99 kg/h 0
FM 1080, 65 kg/h 0
DMF 0 1250 kg/h
MOL 0 541, 23 kg/h
Total 1880, 64 kg/h 1791, 23 kg/h

Calcul du titre x
Nous avons la formule générale de l’enthalpie :

∆H = n0 ·x·∆rH(25◦ C)+[n0 · x · (Cp (DMF) + Cp (Mol)) + n0 · (1 − x) · (Cp (DMA) + Cp (FM))]·∆T

Étant donné qu’il n’y a pas d’échange de chaleur, cela signifie que ∆H = 0. La formule
devient donc :

0 = n0 ·x·∆rH(25◦ C)+[n0 · x · (Cp (DMF) + Cp (Mol)) + n0 · (1 − x) · (Cp (DMA) + Cp (FM))]·∆T

En isolant x, on obtient l’expression suivante :

2
 − (Cp (DMA) + Cp (FM)) · ∆T
x=
∆rH(25◦ C) + (Cp (DMF) + Cp (Mol) − Cp (DMA) − Cp (FM)) · ∆T
Nous utilisons maintenant les valeurs suivantes :

• Cp (DMF) = 148.16 J/mol · K

• Cp (Mol) = 83.3 J/mol · K

• Cp (DMA) = 136.77 J/mol · K

• Cp (FM) = 120.57 J/mol · K

• ∆T = 35 °C

• ∆f H ◦ (DMF) = −239.0 kJ/mol

• ∆f H ◦ (Mol) = −238.4 kJ/mol

• ∆f H ◦ (DMA) = −49.8 kJ/mol

• ∆f H ◦ (FM) = −378.0 kJ/mol

Calcul de ∆rH(25◦ C) :

∆rH(25◦ C) = ∆f H ◦ (DMF) + ∆f H ◦ (Mol) − ∆f H ◦ (DMA) − ∆f H ◦ (FM)

∆rH(25◦ C) = (−239.0)+(−238.4)−(−49.8)−(−378.0) = −477.4+49.8+378.0 = −49.6 kJ/mol

Ensuite, nous avons calculé x :

−(136.77 + 120.57) · 35
x=
(−49.6 × 1000) + (148.16 + 83.3 − 136.77 − 120.57) · 35
Ce qui donne :

x ≈ 0.178 ou 17.8%
Ainsi, le titre x est d’environ 17.8 %.

Calcul de Q
Le calcul de Q repose sur les deux formules suivantes :

Q = ∆H = n0 · (1 − x) · ∆rH(60◦ C)
∆rH(60◦ C) = ∆rH(25◦ C) + ∆Cp · ∆T
Où :

∆rH(25◦ C) = −49.6 kJ/mol, ∆Cp = −25.88 J/mol·K , ∆T = 35 K

3
 Le calcul de ∆rH(60◦ C) donne :
∆rH(60◦ C) = −50.505 kJ/mol
Ensuite, avec les données suivantes :
n0 = 17101.267 mol/h, x = 17.8, ∆rH(60◦ C) = −50.505 kJ/mol
Le résultat obtenu pour Q est :
Q ≈ 197.21 kW
article amsmath graphicx float

Tableau des données détaillées

Nous avons calculé la masse et la fraction massique de chaque élément du mélange en
fonction de x, et ensuite, nous avons calculé la masse volumique du mélange en utilisant
la formule suivante :
1
ρmélange = P  
ωi
ρi

Où :
• ωi est la fraction massique de l’élément i,
• ρi est la masse volumique de l’élément i.

Table 1: Données : Conversion, Masse, Fraction massique et Masse volumique

X (%) DMA FM DMF MOL ρ (kg/m3 )
m (Kg) ω m (Kg) ω m (Kg) ω m (Kg) ω
17.8 636.3 0.352 847.6 0.469 223.4 0.123 97.9 0.054 828.25
54.97 348.56 0.1931 464.32 0.257 689.9 0.382 302.43 0.167 855.78
89.6 80.5 0.0446 107.24 0.059 1124.52 0.623 492.95 0.273 883.35
100 0 0 0 0 1250 0.694 549.35 0.305 892.643

Calcul de la Viscosité du Mélange

La viscosité du mélange ηmélange est calculée par la relation suivante :
X
ln(ηmélange ) = yi ln(ηi )
Où :
• yi est la fraction molaire de l’élément i,
• ηi est la viscosité de l’élément i,
• ηmélange est la viscosité du mélange.

4
 Table 2: Tableau des données de DMA et viscosité
X (%) DMA et FM DMF et MOL Viscosité
n ξ n ξ µ (x 10−4 )
17.8 14114.7 0.411 3056.47 0.089 2.246
54.97 7732.19 0.225 9439.01 0.275 3.075
89.6 1785.8 0.052 15325.39 0.448 4.198
100 0 0 17171.2 0.5 4.497

1 Algorithme de calcule:
Nous avons calculé la viscosité dynamique, la conductivité thermique et la capacité calori-
fique à l’entrée et à la sortie du mélange. Les valeurs moyennes ont été déterminées en
prenant la moyenne de ces propriétés entre l’entrée et la sortie. Nous avons également
calculé le nombre de Reynolds et le nombre de Prandtl à partir des propriétés du mélange.
Ces paramètres sont essentiels pour l’analyse des écoulements et des transferts de chaleur
dans le système.

5
 Table 3: Données expérimentales pour différentes conditions
Variable TRONOANA0 TRONCON1 TRONCON2 TRONCON3
Volumétrie 0.1974 × 10−3 0.5185 × 10−3 0.5185 × 10−3 0.5185 × 10−3
V /V t 11.26% 40.84% 70.42% 100%
x (Entrée) 0% 17.8% 55.47% 86.23%
x (Sortie) 17.8% 55.47% 86.23% 100%
Q (KW) 90.74 27.91 3.84
µentrée 2.246 × 10−4 3.075 × 10−4 4.418 × 10−4
µsortie 3.075 × 10−4 4.418 × 10−4 4.497 × 10−4
µmoyenne 2.66 × 10−4 3.59 × 10−4 4.307 × 10−4
λentrée 0.1587 0.1708 0.2354
λsortie 0.1708 0.2354 0.1853
λmoyenne 0.1647 0.2031 0.2103
Cp (Entrée) 127.41 127.89 128.95
Cp (Sortie) 127.89 128.95 128.22
Cp (Moyenne) 127.65 128.42 128.58
Reynolds 1557.22 1153.81 961.74
Prandt 0.206 0.2269 0.2633

1.1 methodologie de calcule:

Dans cette section, les calculs ont été effectués pour déterminer les principales car-
actéristiques thermiques du système. Nous avons utilisé les formules suivantes pour
obtenir les résultats.
La surface d’échange totale S est calculée à partir de la chaleur transférée Q et du
coefficient d’échange global U ainsi que de la différence de température ∆Te :
Q
S=
U · ∆Te
Calcul du nombre de tubes Nt
La surface d’échange par tube Se est calculée comme suit :

Se = πDext L
Ensuite, le nombre de tubes Nt est déterminé en fonction de la surface d’échange totale
:

6
 S
Nt = ⌈ ⌉
Se
Calcul du pas tubulaire Pt
Le pas tubulaire est donné par :

Pt = 1.25 · Dext
Calcul du diamètre du faisceau Dc
Le diamètre du faisceau est calculé par :
 1
 2.207
Nt
Dc = Dext + 0.068
0.249
Calcul du nombre de Reynolds Re
Le nombre de Reynolds pour le fluide est donné par :

4 · Qm /3600
Re =
π · µ · Dint
Calcul du nombre de Prandtl P r
Le nombre de Prandtl pour le fluide est calculé par :
Cp · µ
Pr =
λ
Calcul du nombre de Nusselt N u
Le nombre de Nusselt est obtenu à partir des nombres de Reynolds et de Prandtl :

N u = 0.023 · (Re0.8 ) · (P r0.4 )

Calcul du coefficient d’échange convectif interne hm
Le coefficient d’échange convectif interne est donné par :
Nu · λ
hm =
Dint
Calcul du coefficient d’échange convectif externe he
Pour l’eau, le nombre de Reynolds Reeau et le nombre de Prandtl P reau sont calculés
de la manière suivante :
G · Deq
Reeau =
µeau
Cpeau · µeau
P reau =
λeau
Ensuite, le coefficient d’échange convectif externe he est calculé par :
0.8 0.4
λeau · 0.023 · (Reeau ) · (P reau )
he =
Deq
Calcul du coefficient d’échange global U
Le coefficient d’échange global U est donné par :

7
  
Dext
1 1 Rtube Rcal 1 ln Dint
= + + + +
U he · Aext Aext Aint hm · Aint 2πλcycle · L
Où :
- Aext = π · Dext · L - Aint = π · Dint · L
En inversant cette relation, on obtient U .

Tableau des Résultats

Les résultats des calculs sont présentés dans le tableau ci-dessous :

Paramètre Valeur
Surface d’échange totale (S) 99.2804 m²
Surface d’échange par tube (Se) 0.3654 m²
Nombre de tubes (Nt) 272
Pas tubulaire (Pt) 0.0200 m
Diamètre du calandre (Dc) 0.4489 m
Nombre de Reynolds (Re) 646.9797
Nombre de Prandtl (Pr) 0.2062
Nombre de Nusselt (Nu) 2.1683
Espace entre chicanes (Ec) 0.0898 m
Surface d’écoulement croisé entre chicanes (A) 0.0081 m²
Diamètre hydraulique du calandre (Deq) 0.0114 m
Nombre de Reynolds (Re) de l’eau 1522.5247
Nombre de Prandtl (Pr) de l’eau 5.6949
Coefficient d’échange externe (he) 842.2040 W/m².K
Vitesse massique côté calandre (G) 107.2125 m²
Débit par tube (Qm) 6.6176 kg/s
Coefficient d’échange interne (hm) 26.2590 W/m².K
Surface extérieure (Aext) 0.3654 m²
Surface intérieure (Aint) 0.3106 m²
Coefficient d’échange global (U) 45.8134 W/m².K
Table 4: Tableau récapitulatif des résultats calculés pour le premier trancon

Tableau de deusieme tronçon :

on a calculer les resultat de deusime trancons :

8
 Paramètre Valeur
Surface d’échange totale (S) 27.0532 m²
Surface d’échange chaque tube (Se) 0.3654 m²
Nombre de tubes (Nt) 75
Pas tubulaire (Pt) 0.0200 m
Diamètre du calandre 0.2805 m
Nombre de Reynolds (Re) 1738.5433
Nombre de Prandtl (Pr) de fluide 0.2270
Nombre de Nusselt (Nu) 4.9692
Espace entre chicanes 0.0561 m
Surface d’écoulement croisé entre chicanes (A) 0.0031 m²
Diamètre hydraulique calandre (Deq) 0.0114 m
Nombre de Reynolds (Re) de l’eau 1199.6411
Nombre de Prandtl (Pr) de l’eau 5.6949
Coefficient d’échange externe (he) 695.9963 W/m².K
Vitesse massique côté calandre (G) 84.4758 m²
Débit par tube (Qm) 24.0000 kg/s
Coefficient d’échange interne (hm) 74.2086 W/m².K
Surface extérieure (Aext) 0.3654 m²
Surface intérieure (Aint) 0.3106 m²
Coefficient d’échange global (U) 51.8564 W/m².K
Table 5: Tableau récapitulatif des résultats calculés

Tableau de troisième tronçon:

et voici les resultats de troisime tronçon:

9
Paramètre Valeur
Surface d’échange totale (S) 3.7221 m²
Surface d’échange chaque tube (Se) 0.3654 m²
Nombre de tubes (Nt) 11
Pas tubulaire (Pt) 0.0200 m
Diamètre du calandre 0.1570 m
Nombre de Reynolds (Re) 9880.3802
Nombre de Prandtl (Pr) de fluide 0.2633
Nombre de Nusselt (Nu) 21.1714
Espace entre chicanes 0.0314 m
Surface d’écoulement croisé entre chicanes (A) 0.0010 m²
Diamètre hydraulique calandre (Deq) 0.0114 m
Nombre de Reynolds (Re) de l’eau 526.4876
Nombre de Prandtl (Pr) de l’eau 5.6949
Coefficient d’échange externe (he) 360.1441 W/m².K
Vitesse massique côté calandre (G) 37.0740 m²
Débit par tube (Qm) 163.6364 kg/s
Coefficient d’échange interne (hm) 327.3789 W/m².K
Surface extérieure (Aext) 0.3654 m²
Surface intérieure (Aint) 0.3106 m²
Coefficient d’échange global (U) 51.4376 W/m².K
Table 6: Tableau récapitulatif des résultats calculés

10