Introduction :

Dans la pratique industrielle, les réacteurs se comportent rarement comme des réacteurs idéaux (piston ou
parfaitement agité) : une géométrie complexe (présence d’équipement de contrôle et régulation…etc.), une
mauvaise distribution du fluide …etc peuvent perturber l’écoulement ; ce qui en résulte que les performances d’un
réacteur réel s’éloignent de ses performances prédites par les modèles idéaux.

Dans un réacteur réel, les molécules séjournent dans le volume réactionnel pendant ses temps « ts» qui dépendent
notamment du profil hydrodynamique, et de la géométrie du réacteur.Ces temps peuvent s’écarter notablement du
temps de séjour théorique, donc il existe une distribution des temps de séjour (DTS) et cette dispersion exerce
une influence sur les performances chimiques du réacteur .

Deux approches ont été développées pour décrire les écoulements dans les réacteurs .D’une part, une approche
fondée sur les principes de la physique théorique avec leséquations de Navier Stokes par exemple, qui vise à
décrire la position, la vitesse et lapression d’un élément fluide à tout instant et en tout point de l’espace. D’autre
part, uneapproche systémique qui cherche à décrire cet écoulement par des informations plusglobales mais
indispensables à la description des comportements transitoires des solutésréactifs que le fluide en écoulement
transporte. La première approche n’est possible quepour des systèmes de géométrie simple et connue. De plus son
développement demande des moyens de calculs importants. Quant à l’approche systémique, elle nécessite de
définir a priori le degré de complexité désiré de l’écoulement. Elle est largement utiliséeen génie chimique.

But De La Manipulation:
L’objectif principal de ce TP est de Modéliser l’écoulement dans une cascade de deux réacteurs agités.

Généralités :
Dans la réalité, les hypothèses ne sont pas respectées et ceci d'autant plus si le réacteur est grand. Ainsi dans
un réacteur de plusieurs milliers de litres, la température n'est pas uniforme, ce qui engendre des profils de
concentrations, car la vitesse de réaction, dépendante de la température, n'est pas uniforme dans la masse
réactionnelle.

Afin de mesurer l'idéalité d'un réacteur, il existe une technique, la distribution de temps de séjour (DTS). Cette
technique permet, via la concentration d'un traceur à différents endroits du réacteur, de comparer un réacteur
avec les modèles ci-dessus et le cas échéant de corriger le modèle pour tenir compte des déviations.

L’expression de distribution de temps de séjour s'utilise en génie des procédés. Elle permet de rendre compte du
type d'écoulement d'un fluide dans un réacteur (réservoir, conduite…).

La distribution de temps de séjour est un modèle qui permet de caractériser l'hydrodynamique d'un réacteur
chimique et de déterminer quel modèle de réacteur définit le mieux l'installation étudiée (réacteur continu ou
réacteur tubulaire). Cette caractéristique est importante pour pouvoir calculer la performance d'une réaction
avec une cinétique connue.

Les modèles des réacteurs idéaux sont construits sur un certain nombres d'hypothèses: le mélange d'un flux
entrant d'un réacteur continu est considéré comme complet et instantané avec le milieu réactionnel et dans un
réacteur tubulaire, l'écoulement est défini comme piston (pas de rétro-mélange). Or dans la réalité, il est
impossible d'obtenir de telles conditions, notamment pour des réacteurs industriels qui ont en général une taille
comprise entre 1 et plusieurs dizaines de m3. La distribution des temps de séjour permet donc de déterminer
dans quelle mesure une installation dévie du modèle idéal et d'apporter les corrections nécessaires dans le
fonctionnement du procédé afin de compenser cette non-idéalité.
Pour caractériser l'écoulement, on utilise dès lors la distribution de temps de séjour qui est une approche statistique. En
effet, on considère un élément du fluide à son entrée dans le réacteur et on mesure le temps que ce dernier met pour
atteindre la sortie. Si on répète l'expérience ou on considère plusieurs éléments en même temps, on constatera que les
résultats ne sont pas identiques. Les causes principales de cette différence est l'existence de zones stagnantes ou de
court-circuits dans l'installation étudiée.

On peut dès lors établir une distribution des temps de séjour, le plus souvent représentée par une distribution de
fréquences appelée habituellement E. Pour ce faire, trois hypothèses sont posées:

. Le réacteur est à l'état stationnaire . Le fluide est incompressible à l'entrée et à la sortie du réacteur,
. Le transport a lieu uniquement par convection .

Définitions générales :

Définition du réacteur :

Par définition, un
réacteur est un appareil
dans lequel on effectue
une transformation sur Réacteur fermé :
la matière, en réglant
les conditions pour C’est un appareil dans lequel on introduits les réactifs et on laisse la
obtenir un produit réaction se dérouler au cours du temps.

déterminé, mais le
terme de réacteur peut
fort bien s’étendre à
des appareils qui ne sont
le siége d’aucune
Notions d’âge, d’espérances de vie, et de temps de séjour :
Avant de définir précisément la distribution des temps de séjour, il est intéressant dedonner quelques définitions
qui caractérisent chaque molécule présente dans un réacteur desection d’entrée «E» et de section de sortie «S».

Considérons une molécule «M» à un instant donné «t» son temps de présence dans leréacteur à partir du moment
où elle a franchie l’entrée «E» s’appelle l’age interne α de cette molécule, le temps qui s’écoulera entre le
moment α et la sortie du réacteur de la molécule
E(t), distribution des temps de séjour, est représente son espérance de vie λ. Son temps de séjour
telle que E(t).dt est lafraction du débit de réel, ts qui représente la duréetotale de présence de la
sortie qui est restée dans le réacteur un molécule au sein du réacteur peut s’écrire de la manière
tempscompris entre t et t+dt.
suivante

Propriétés de E(t) :
Définitions de E(t) et F(t):
-E(t) est toujours positive∞

-E(t) est normée :

fraction de débit restée entre t1 et t2
dans le reacteur F(t), distribution cumulée des temps de séjour, est la
- E(t) [T-1] fraction dedébit qui est restée moins de t dans le
réacteur (plus jeune que t).



Réacteur parfaitement agité :

C’est un appareil dans lequel a lieu

simultanément uneintroduction et une
soustraction de matière. L’agitation dumélange
réactionnel est parfaite au point ou les
conditionsrégnant à l’intérieur du réacteur
correspondent à celles desortie. Hypothèses fondamentales :
La théorie de DTS est essentiellement définie et utilisée
pour des milieux continus etmonophasiques, elle repose sur
 les hypothèses suivantes :
1- Ecoulement en régime permanent.
2- Ecoulement déterministe : une même cause produira les mêmes effets, difficultéd’application par exemple dans
les régimes intermédiaires entre laminaire et turbulent stabilisé.
3- Le fluide est incompressible.
4- L’écoulement à travers les sections d’entrée et de sortie se fait uniquement par convection forcée, à
l’exclusion de toute diffusion et de tout mélange en retour.
5- Les conduits d’entrée/sortie sont de petit diamètre devant les dimensions du réacteur.

Détermination expérimentale des

temps de séjour (traceur) :
Pour mesurer une distribution des temps de séjour
on utilse la methode de traceur, La méthode des
traceurs consiste à marquer les molécules entrant
dans le réacteur et à suivre leur histoire jusqu’à la
sortie, en particulier en les dénombrant dans le courant de sortie en fonction du temps. Pour cela, on utilise un
traceur qui peut être un traceur radioactif, un colorant, un sel ou toute substance de mêmes propriétés
hydrodynamiques que le fluide, mais décelable par une propriété physique caractéristique telle que conductivité
électrique, thermique, etc.

On applique, à l’aide du traceur, un certain signal à l’entrée du réacteur, de manière à ne pas perturber le
régime d’écoulement, et on examine sa réponse à la sortie. On en déduit des renseignements sur le
comportement interne du fluide. Dans les systèmes ou l’entrée et la sortie n’ont pas une section négligeable, le
marquage aussi bien que la détection doivent prendre en compte la pondération en flux :

injection homocinétique dans toute la section droite avec la vitesse locale du fluide, détection représentative de
la moyenne sur les débits (concentration de mélange). Lorsque cela est possible, il est préférable d’injecter et de
détecter le traceur dans des zones rétrécies où la valeur
du nombre de Reynolds est élevée et où règne donc une
2. Injection impulsion (ou Dirac) :
turbulence spontanée assurant la condition du mélange. Soit un réacteur alimenté en continue par un
Des mesures qui ne respecteraient pas ces précautions débit Q (régime permanent).
seraient dépourvues de signification. Introduction d’un nombre de moles n0 pendant
un temps très court.Mesure de la concentration
Il existe plusieurs méthodes de déterminer la DTS dont cI (t) en sortie de réacteur.
on donne seulement deuxprincipales et simples utilisés
généralement dans les études expérimentales qui sont :
-injection échelon
-injection impulsion (ou Dirac)
On note que la 2éme méthode qui est la plus utilisée
généralement par les chercheurs et quisera utilisé dans
notre recherche.

Traceur :
Il marque le fluide, sans
le perturber.
-Mêmes propriétés que
le fluide
-Inerte
-Mesurable en continu
-Injection et détection
Bilans en traceur :

Distribution des temps de séjour dans les réacteurs idéaux :

Réacteur piston :
Dans un réacteur piston, le fluide avance en bloc sans se mélanger. La bouffée detraceur ressort donc au bout
d’un temps t = τ (le temps de passage dans le réacteurpiston sous la forme d’une impulsion étroite.
Réacteur parfaitement agité :
Dans un réacteur parfaitement agité le traceur se répartit instantanément dans toute lamasse, puis il est
graduellement élué suivant une décroissance exponentielle.
La réponse à une injection-impulsion est :

Dans un réacteur de ce type, le mélange est uniforme et instantané. A chaque t lescaractéristiques du fluide à
une injection impulsion à l’intérieur et à la sortie du réacteursont donc identiques. Pour établir l’expression de E
(t ) , cherchons la réponse à une injectionimpulsion, introduisons une bouffée de traceur dans le réacteur la
concentration monte instantanément à C0puis évolue suivant l’équation :

1- Réacteur parfaitement mélangé :

1. Injection échelon :
Dans ce type d’injection la concentration en
traceur passe brusquement de 0 à C0 dans 2- Réacteur piston :
l’alimentation à l’instant t = 0,la concentration en
sortie est notée C (t).
Le suivi type de cette concentration est présentée
sur la figure .
C0 : Concentration initiale supposée répartie dans
tout le volume.

formules :
 PARTIE EXPERIMENTALE

MODE OPERATOIRE :

Courbe D étalonnage  :
On a preparé des solution en pesant respectivement– 3 – 4 – 10 Grammes de NaCl .
Et les intoduire dans des fioles puis ona completé avec de l ’eau distillée jusqu’au 100 ml.
Aprés on a mesuré la conductivité de chaque solution preparée par un conductimètre .
Les valeurs de la conductivité en fonction de la concentration pour l’étalonnage sont portées dans le tableau ci-
dessous :

NaCl 3,0% 4,0% 10%

λ mS 70.3 113.2 194.9

Concentratio 0.51 0.68 1.71

n

C=m/M.V
Construction de la courbe d’étalonnage :
250

200

150
Axis Title

Conductiv(ms)
100
Linear (Conductiv(ms))

50

0
0 1 2 3 4 5 6 7 8
Axis Title

Calcul des volumes des deux réacteurs  :

-Eventuellement ona vindagé le réacteur 1 .en ouvrant le robinet de purge du réacteur et en établirant un
débit d’eau qu’on a maintiendra constant . ce débit est noté Q1.

- Puis on a fermé le robinet de purge et on a déclenché le chronomètre pour mésurer le temps de passaget 1.

Puis ona effectué les même étapes pour le second réacteur .

Expérience :
-On a arreté l’agitation du 1er réacteur et on a fermé son robinet de vidange . Puis on a ajusté le niveau
de l’eau dans le 2ème réacteur de manière à remplir presque à moitié .

- On a dissout 10 gr de NaCl dans un petit volume d’eau et l’introduire dans le 2 ème réacteur tout en
réglant le débit constant .

- On a effectué ensuite des prélèvement de NaCl du second réacteur .Telque le prélèvement se fera dans
un premier temps toutes les 30 s puis l’on espacera à partir de 3ème de 1 min .

- Enfin on a mésuré la conductivité de chaque prélèvement .

Remarque  :
A l’aide de la courbe d étalonnage « Conductivité = F( Concentration ) » , nous déterminons les
concentrations des solutions

.
 Partie calcul :
 Calcul des volumes des deux réacteurs :
V eau 400
Q1 = = =30 .77 ml/ s
 Le débit de l’eau :
t 1 13
 Temps de passage = τ1 = 78 sec.

 Le volume du premier réacteur = 1V =Q ⋅τ =30 . 77⋅78=2. 4 l

1 1
Etant donner que le deuxième réacteur a les mêmes dimensions que le premier, donc sont volume est de 2.4 L
Les valeur experimentales de C sont deduites a partir de la courbe d’etalonnage : la courbe est une droite tq :
f(c)=a.c= a= ( 70 – 0 )/(0.5-0)= 140 C= λ/140
Les valeurs theorique de C sont obtenues apartir de l’equation suivante : C/Co=Exp(t/T )
 Les résultats des calculs sont portés dans le tableau ci-
dessous :

t (sec) 0 30 60 120 180 240 300

Λ 16.7 13.5 9.4 6.8 5.6 4.6

Cexp 0.119 0.096 0.067 0.048 0.04 0.033

C/C0exp 1 1.67 1.35 0.94 0.67 0.56 0.46

C/C0th 1 0.68 0.46 0.215 0.1 0.046 0.021

Cthe 0.048 0.032 0.015 0.007 0.0033 0.0015

Eabs 0 0.99 0.89 0.725 0.57 0.514 0.44

Erela 0 1.455 1.93 3.37 5.7 11.17 20.95

 Co= m/M.V = 10/58.8 . 2.4 = 0.071

 L’erreur absolue : Eabs =| (C/C0exp )-(C/C0th) Eabs(moy)
=0.688

 L’erreur relative : Erel = Eabs / C/C0th
Erela(moy)=7.43
 Calcule le temps de sejour moyen et la variance :

1-Le temps de sejour moyen :
 emps de sejour moyen : Ts=∑
 T Ci.ti.∆ti/∑Ci.∆ti ∆ti :est le pas
 Ts(th) =
 Ts(exp)=
 La différence entre les deux temps est due à l’imperfections
des appareils
 2− La variance :б²/Ts²=1/N N: nombre de réacteur
2
 б²= Ts²/N =
  La figure – 2 représente le graphe C/C0) exp= f (t)

C/Co
1.8

1.6

1.4

1.2

1 C/Co

0.8

0.6

0.4

0.2

0
0 50 100 150 200 250 300 350

C
 La figure – 3 représente le graphe
( )
C0 the
= f (t).

C/Co
1.2

0.8
C/Co
0.6

0.4

0.2

0
0 1 2 3 4 5 6 7 8 9

CAS 2 :

T(s) 0 30 60 120 180 240 300

Conductivite 0 7.1 7.3 7.7 7.5 6.9 6.3
(ms)
(C/Co)exp 1 0.05 0.052 0.055 0.0535 0.049 0.045

(C/Co)th 1 0.26 0.356 0.33 0.229 0.141 0.082

Calcule de Cexp ; comme le cas 1

Calcule de Cth ; C/Co = t/T exp(-t/T) T : le temps de passage 78s

Co = 1mol

2-Le temps de sejour moyen ts ;

a- Ts(th)=
b- Ts(exp)=

3- La variance б²=

La variation de (C/Co)exp :

Valeur des Y
0.35

0.3

0.25

0.2 Valeur des Y

0.15

0.1

0.05

0
0 1 2 3 4 5 6 7 8

La variation de (C/Co)th :

Valeur des Y
0.4

0.35

0.3

0.25
Valeur des Y
0.2

0.15

0.1

0.05

0
0 1 2 3 4 5 6 7 8 9

4-La modélisation de l’écoulement :

μ = Ts/г tq :г=г1 + г2 =78*2 = 156

Interpretation :
D’après Les résultats que nous avons obtenus nous remarquons que les graphes obtenus ne sont pas parfaits
Elles n’ont pas la forme voulue hyperbole sous la forme d’exponentielle.
Nous remarquons que la conductivité diminue avec la dilution de notre solution en fonction du temps de passage du
1er réacteur au 2eme.
La valeur du concentration diminue en fonction du temps et cela à cause de la dilution provoqué par le premier
réacteur qui est liée avec le second par cascade et cela nous permet de dire que pour améliorer la conversion il
faut ajuster les réacteurs en série.

Conclusion :
La présente étude a pour but la modélisation numérique des écoulements dans une cascade de deux réacteurs
agitésutilisant l’analyse des courbes de distribution des temps de séjour (DTS) et selon les deux cas : Purge selon
le cas 1 et 2 respectivement tout en faisant une comparaison entre les deux cas afin d’examiner si les modelé
est convenable ou pas et afin de choisir le bon modèle pour améliorer les performances du réacteur .
la modélisation numérique a permis d’obtenir des résultats qualitatives par rapport à l’approche physique
tels que : la possibilité de l’obtention de la courbe d’étalonnage ainsi que les différents courbe de l’évolution des
concentration au cours du temps afin de déterminer la bonne déroulement du phénomène, la possibilité de
visualiser des zones de recirculations ,ainsi l’obtention des résultats en temps court et optimal.

Les résultats obtenus restes acceptable malgré les erreurs commissent pendant la manipulation et cela est
visualisé d’après les graphes obtenus ce qui implique le mauvaise déroulement du manipulation quel que soit du part
des étudiant ou bien des matériels utilisés tel que le chronométré Les résultats obtenus sont en bonne
concordance avec d’autres résultats numériques.

Les conclusions tirées de ce travail sont les suivantes :

1-le rôle important des dispositifs de contrôles dans l’augmentation du temps de séjour moyen.
2- la méthode des traceurs est une bonne méthode pour étudier les phénomènes qui se déroulent au niveau des
réacteurs.
3- la DTS nous a permis de déterminer la configuration optimale et qui est le répartiteur avec déflecteur.
4-le rôle important que joue la DTS dans la compréhension du phénomène d’écoulement du fluide au sein du
répartiteur
5- le temps de séjour influe sur les performances chimiques des réacteurs ainsi que l’agitation des molécules
entre elles.

Enfin en peut dire que le modèle utilisée pendant ce TP reste un peu convenable comme un réacteur.