Génie de la Réaction Chimique

extrait : Distribution des Temps de Séjour

CGP215
Génie de la réaction chimique & Évaluation économique des procédés

Marie DEBACQ
marie.debacq@cnam.fr

version du cours de GRC de février 2009 – dernière révision février 2016

attribution, pas d’utilisation commerciale, partage dans les mêmes conditions
 écoulement dans les réacteurs réels :
Distribution des Temps de Séjour

Ce chapitre est consacré à une méthode très efficace pour l'étude des écoulements dans les réacteurs réels : la
Distribution des Temps de Séjour (DTS). Il s'agit en réalité d'une méthode d'étude applicable à bien d'autres cas
que les réacteurs : des appareillages de l'industrie chimique, aux réseaux urbains d'évacuation des eaux usées, en
passant par les zones portuaires. Il s'agit aussi de l'exemple le plus connu d'une méthodologie d'étude très efficace
en génie des procédés (et ailleurs) : la méthode systémique, qui consiste à obtenir des informations macroscopiques
utiles à partir de mesures expérimentales simples.
Une séance de Travaux Pratiques illustrera ce chapitre.

On se propose d'étudier l'écoulement de liquides, de gaz ou de solides divisés en circulation dans des installations
industrielles (réacteurs, séparateurs) ou dans des milieux naturels (milieux poreux, bassins fluviaux, zones
côtières…). Dans le cas des installations industrielles, il s'agit de disposer des informations nécessaires sur
l'écoulement pour prévoir leur capacité de production, et détecter d'éventuels défauts. Dans le cas des milieux
naturels il s'agit de prévoir les risques de pollution, de comprendre les ensablements, etc. Dans tous les cas, on
souhaite obtenir une information GLOBALE sur l'écoulement dans le système considéré.
On dispose aujourd'hui de codes de calcul permettant de déterminer l'écoulement dans toutes sortes de systèmes.
Mais dans de nombreuses situations rencontrées dans la pratique (systèmes polyphasiques, fluides aux propriétés
mal connues, géométries très complexes), les résultats de ces simulations sont quelque peu hasardeux et nécessitent
des moyens et des temps de calcul sans commune mesure avec l'objectif poursuivi. En outre, même lorsqu'il est
faisable, ce calcul conduit à une description inutilement détaillée (champ de vitesse en tout point du système) au
regard des besoins simples. Enfin toute approche théorique devant être confrontée à des données expérimentales,
la mesure de Distribution des Temps de Séjour reste très pertinente, malgré les progrès constants et rapides en
mécanique des fluides numérique.

Il est rare que tous les éléments fluides ou solides qui pénètrent au même instant dans un système (réacteur, port,
etc.) sortent ensemble de ce volume. Ces éléments n'ont pas tous le même "temps de séjour" individuel : comme
l'illustre la Figure 1, certains  passent directement de l'entrée à la sortie, d'autres  y séjournent plus ou moins
longtemps ; il peut exister des zones de recirculation , ou encore des zones immobiles  (zones mortes). La
distribution des temps de séjour dans un système est une caractéristique importante de celui-ci, car elle influence
de façon notable les performances d'un réacteur ou l'efficacité d'une roche comme barrière contre la pollution.


 

Figure 1 : Distribution des temps de séjour.

 fonction de distribution des temps de séjour E

C'est le mérite de DANCKWERTS d'avoir imposé dès 1953 une approche systémique à ce problème en définissant une
fonction de distribution aisément accessible à l'expérience. La fonction de distribution des temps de séjour E est
telle que "la fraction de fluide qui séjourne dans l'installation pendant un temps compris entre ts et ts+dts égale le
produit E(ts).dts", comme illustré sur la Figure 2.
 E

ts
dts
Figure 2 : Fonction de distribution des temps de séjour.

Le but de ce chapitre est d'étudier les distributions des temps de séjour (DTS) dans les systèmes en écoulement, de
voir quelles informations on peut en tirer quant à leur fonctionnement hydrodynamique et de rendre compte de
ces distributions à l'aide de modèles représentatifs. Dans tout ce qui suit la notion de "réacteur" est prise au sens
large. Les méthodes décrites s'appliquent en effet à toute portion d'espace traversée par des éléments mobiles qui y
séjournent pendant des durées variables : aussi bien les appareils industriels, que les bassins fluviaux, les organes
d'un être vivant, ou une autoroute la veille d'un week-end prolongé.

1.1. MÉTHODES DE MESURE

Pour caractériser l'hydrodynamique d'un système on peut imaginer de suivre l'histoire d'un traceur introduit dans
l'écoulement.

1.1.1. principe

La méthode des traceurs consiste à "marquer" des molécules entrant dans le système, et à les dénombrer dans le
courant de sortie en fonction du temps. Ainsi on applique, à l'aide d'un traceur, un certain signal à l'entrée du
système et on examine sa réponse en sortie, pour en déduire des renseignements sur le comportement interne du
fluide sous la forme d'une fonction de transfert : c'est donc une méthode systémique.
L'introduction du traceur, qui ne doit pas perturber l'écoulement dans le milieu étudié, peut se faire suivant une
injection de type impulsion (Figure 3-a) ou échelon (Figure 3-b).
C/C0 C/C0

ts ts
(a) (b)
Figure 3 : Injections de type impulsion (a) ou échelon (b).

L'injection impulsion consiste à introduire N moles de traceur pendant un temps très court (au moins 100 fois
inférieur au temps de passage ). Le volume de traceur introduit doit être très faible par rapport au volume V du
système afin de ne pas perturber son écoulement. On note C0 = N / V, la concentration qu'aurait le traceur
uniformément réparti dans tout le volume. La réponse C(ts) en sortie du système, une fois normalisée par son
intégrale, est la fonction de distribution E :
C (t s )
E (t s ) 

 Ct   dt
0
L'injection échelon consiste à faire passer brusquement la concentration en traceur dans le flux d'alimentation du
système d'une concentration nulle à une concentration C0, sans perturber l'écoulement dans le système. La réponse
normalisée est appelée "fonction F".

1.1.2. conditions d'application de la méthode

L'application de cette méthode est soumise à 3 hypothèses principales : l'écoulement doit être en régime
permanent, le fluide doit être incompressible et l'écoulement doit être déterministe, c'est à dire qu'il ne fait pas
intervenir de processus aléatoires macroscopiques, du moins à grande échelle de temps et d'espace (un processus
aléatoire stationnaire tel que la turbulence est quasi-déterministe à l'échelle de l'expérimentateur).
Toutefois, les résultats présentés ici peuvent être généralisés aux cas d'écoulements en régime transitoire, aux
fluides compressibles ou dilatables, ou au cas d'alimentations multiples, mais ces aspects ne seront pas détaillés ici.
Des conditions supplémentaires sont imposées en entrée et en sortie du système : l'écoulement dans les conduites
d'entrée et de sortie doit être SANS diffusion (type piston), et le fluide doit être homogène dans la zone de mesure.

1.1.3. choix du traceur

Le traceur choisi doit respecter certaines caractéristiques : avoir les mêmes propriétés hydrodynamiques que le
fluide, être décelable par une propriété physique, être stable et parfaitement neutre pour le système (mesures dans
les milieux naturels, ou dans les installations de production alimentaire…). Il faut enfin être sûr de pouvoir se
débarrasser de tout le traceur.
Plusieurs types de traceurs et de méthodes de détection peuvent être envisagés : les plus courantes sont listées
dans le Tableau 1.

Tableau 1 : Exemples de couples traceur / moyen de détection.

traceur détection
sel conductimétrie
colorant colorimétrie
isotope radioactif gammamétrie

1.2. PROPRIÉTÉS DE LA FONCTION DE DISTRIBUTION

Considérons le signal mesuré en sortie du système (Figure 4).

C

ts
ts
Figure 4 : Réponse à une injection.
1.2.1. réponse à une impulsion : la fonction de distribution E

La fonction de distribution des temps de séjour E possède les propriétés usuelles des distributions. En particulier
on peut définir les moments de cette distribution.

moments de la distribution
Le moment d'ordre n de la fonction de distribution des temps de séjour est :

 n   t s n  E t s   dt s
0
C (t s )
Le moment d'ordre 0, est le facteur de normation. Il est unitaire puisque E (t s )  :

 Ct   dt
0
0  1
Le moment d'ordre 1 représente la moyenne de la distribution. Il est relié au temps de séjour moyen par :
1
ts   1
0
Le moment d'ordre 2 représente l'étalement autour du temps de séjour moyen. Il est lié à la variance par :
2
  
  2   1  2  12

0  0 
Le moment d'ordre 3 représente l'asymétrie de la distribution et le moment d'ordre 4 son aplatissement.

1.2.2. réponse à un échelon : la fonction F

La fonction F est simplement l'intégrale de la fonction de distribution des temps de séjour E, elle apporte donc la
même information. La méthode de mesure de la DTS par injection échelon est donc dite "méthode intégrale".
ts
F t s    Et   dt
0
Ceci se retrouve simplement en se rappelant que l'impulsion est la fonction dérivée de l'échelon.

1.2.3. méthode à 2 points de mesure

Il est souvent difficile de réaliser pratiquement l'injection "proprement", en particulier d'avoir un temps d'injection
très faible. Il peut également arriver que l'entrée du système que l'on souhaite étudier ne soit pas accessible, pour
des raisons de sécurité par exemple.
La méthode dite des "deux points de mesure" permet, à partir d'une injection de forme quelconque, de déduire la
DTS à l'aide de 2 mesures réalisées en entrée et en sortie du système (Figure 5).

points de mesure
injection
Figure 5 : Méthode à 2 points de mesure.

Le dépouillement s'effectue alors à partir des concentrations Ce(ts) et Cs(ts) mesurées en entrée et en sortie, à l'aide
du produit de convolution :
ts
C s t s    C e t   E t s  t   dt
0
 1.3. EXPLOITATION QUALITATIVE : DIAGNOSTIC DES ÉCOULEMENTS

La simple observation de la forme de la courbe E obtenue permet d'effectuer un diagnostic de l'écoulement.

1.3.1. les réacteurs idéaux

Il existe deux situations limites d'écoulement :

 L'écoulement piston, caractérisé par le fait que tous les éléments situés dans une section droite
progressent à la même vitesse, comme poussés par un piston. Le temps de séjour est identique
pour tous les éléments et la composition évolue le long du réacteur : on obtient de la sorte un
traitement très uniforme de la matière.
 Le mélange parfait, qui se traduit par le fait que des propriétés physico-chimiques (concentration,
température, etc.) sont les mêmes en tout point du réacteur et identiques à celles de la sortie. Dans
le cas d'un fonctionnement continu, le mélange parfait donne lieu à une importante distribution
des temps de séjour des éléments : certains quittent le système dès l'introduction (court-circuit),
d'autres séjournent très longtemps (par suite de recirculations), et les temps de séjour s'échelonnent
de zéro à l'infini.
L'intérêt de ces deux concepts simples et extrêmes est que leur fonctionnement est bien connu. Par exemple on sait
très bien calculer le taux de conversion chimique dans un réacteur pison ou un réacteur parfaitement agité (RPA).
E réacteur piston
1

quelconque

RPA

0 ts
ts
Figure 6 : Réacteur piston - réacteur quelconque - réacteur parfaitement agité.

1.3.1.1. réacteur piston

Un réacteur piston se comporte comme un retard pur. Toutes les molécules ont le même temps de séjour, c'est le
temps de séjour moyen, et aussi le temps de passage. La réponse à une impulsion est donc simplement une
impulsion décalée dans le temps de t s   :
Ct s   C 0   t s    *

 t s   
On voit que  C t s  dt s  C 0  . Ainsi, Et s  

0
Cette fonction distribution des temps de séjour dans un réacteur piston est représentée sur la Figure 6.

1.3.1.2. RPA
Dans le cas d'une injection impulsion de N moles de traceur, le bilan en traceur sur le réacteur parfaitement

1 pour x  0
*  est la fonction Dirac :  ( x )  
0 pour x  0
agité fonctionnant en régime permanent s'écrit :
dC traceur
0  Q  C traceur  V 
dt
Avec comme condition initiale la concentration C0 = N / V, la solution de cette équation différentielle du premier
ordre sans second membre est :
 t 
C  C 0  exp   s 
  
La concentration en traceur croît donc instantanément de 0 à C0 à l'instant d'injection du traceur, puis décroît
exponentiellement. Elle ne devient inférieure à 1% de C0 qu'au bout de 5  t s .
Notons que le temps de séjour moyen et le temps de passage sont ici aussi égaux : t s    V / Q *.
 t  1
Ainsi,  exp   s  Et s  
   
Cette fonction distribution des temps de séjour dans un RPA est représentée sur la Figure 6.

1.3.2. défauts d'écoulement typiques

Comme on le voit sur la Figure 6, un système réel présente un cas intermédiaire entre, d'une part la situation du
mélangeur parfait, et d'autre part la situation de l'écoulement piston.
Le temps de séjour moyen t s et le temps de passage   V /Q calculé à partir du débit d'alimentation Q du

système, peuvent être différents ; en outre on peut avoir  C t s  dt s  C 0  t s ; en particulier si tout le volume
0
interne n'est pas réellement accessible au fluide.

volume mort court-circuit recirculations

ts <  ts >  t s < 

  C t s  dt s  C 0  t s 
0

Figure 7 : Défauts d'écoulement dans un réacteur agité.

 
 ts 
 
C (t s )
 t s  C t s  dts  t s  C 0  exp   
  dts

* En effet t s   1   t s  E t s   dts   t s   dts  0  0
  
0 0  t
 C (t )  dt  C (t )  dt  C 0  exp      dt
0 0 0
 
 t   t 
Or  C 0  exp      dt   C 0   exp      C 0   0  1  C 0  .
0 0
1 1
D'autre part  x  exp  k  x   dx  k
 x  exp  k  x  
k2
 exp  k  x  [ceci s'obtient par intégration par parties].


  t   t 
C 0     t s  exp   s    2  exp   s  
       0   t   t 
donc t s    t s  exp   s     exp   s   0  0    0  1  
C 0         0
 volume mort court-circuit recirculations
t s <  t s >  t s < 

  C t s  dt s  C 0  t s 
0

Figure 8 : Défauts d'écoulement dans un système quelconque.

C (t s )
 En l'absence de ces défauts d'écoulement, on a E (t s ) 
C 0 
.

1.4. EXPLOITATION QUANTITATIVE : MODÉLISATION

Ayant obtenu par la méthode des traceurs une courbe expérimentale de DTS, on se propose de rendre compte de
l'allure observée au moyen d'un modèle hydrodynamique simple, utilisable commodément par la suite pour
prédire la conversion chimique d'un réacteur ou pour servir de base à son extrapolation.

1.4.1. modèle à dispersion axiale

Il s'agit de représenter les faibles écarts à l'écoulement piston. Le modèle superpose un écoulement piston
purement convectif de vitesse u et une dispersion aléatoire obéissant formellement à une loi de Fick, dans le
système de coordonnées x se déplaçant à la vitesse moyenne u du fluide :
C  2C
D 
t 2
z
D est ici un coefficient de dispersion phénoménologique qui rend compte à la fois du fluide, du régime
d'écoulement et de la géométrie du système.
~
Pour N moles de traceur par unité de surface déposées en x = 0, la solution de l'équation de dispersion est :
~ 
N x 2 
C ( z)   exp  
4  D  t s  4  Dts 
 
L'abscisse axiale z (= x + u.t) est fixe par rapport au système ; ainsi z = 0 en entrée et z = L  u  t s en sortie.
L'observateur posté à la sortie voit la concentration en fonction du temps en écrivant :
x  u  ts  ts  
~
N C t s 
Comme C 0  et E t s   , on a :
L C0  ts

Et s  
1  Pe
  exp   s 
1 /2
 Pe  t  t 2 
 s   
2    t s  t s
 
4  ts  ts 

 
Cette expression fait intervenir le critère de Péclet Pe (produit d'une vitesse caractéristique de l'écoulement par une
uL
distance caractéristique, divisé par un coefficient de dispersion) : Pe  .
D

L'écoulement piston parfait est obtenu pour D = 0, soit Pe infini (dans la pratique pour Pe > 200). Pour les faibles
critères de Péclet, la DTS s'affaisse et devient dissymétrique (cf. Figure 9). On retrouve la DTS d'un RPA pour
Pe  0.
 E
Pe = 20

Pe = 5
0,5

Pe = 1

Pe = 0,1
0
0 1 2 3
t /
Figure 9 : Courbes de DTS en fonction du critère de Péclet.

Les moments de la DTS sont reliés au critère de Péclet, pour les grandes valeurs de celui-ci, selon :
2
t 
Pe  2   s 

 

1.4.2. modèle de cascade de RPA

Lorsque J réacteurs parfaitement agités sont associés en cascade (c'est-à-dire en série), la DTS est la suivante :
J 1  t 
J ts  exp   J s 
 J   t s 
Et s      
t   J  1 !
 s
E
10 RPA

5 RPA
1

2 RPA

0,5

1 RPA

0
0 1 2 3
t /
Figure 10 : Courbes de DTS en fonction du nombre de RPA en cascade.

On retrouve le cas du mélangeur parfait pour J = 1, et le piston parfait pour J infini (dans la pratique pour J > 100).
1.4.3. lien entre les deux modèles

Ces deux modèles sont de natures différentes, en particulier la matière peut remonter le courant (diffusion) dans le
cas "dispersion axiale" mais pas dans le cas "cascade de RPA". Toutefois, à Pe et J élevés, ces modèles coïncident,
comme on le remarque sur la Figure 11. On a alors Pe  2  J .
3

Pe = 100
E

50 RPA

0
0 1 2
t /
Figure 11 : Lien entre le modèle à dispersion axiale et la cascade de RPA.

Ces deux modèles permettent une interpolation continue entre les deux cas idéaux que sont le piston et le
mélangeur parfaits.

1.4.4. représentation d'un réacteur réel

Pour rendre compte du comportement de réacteurs réels où l'écoulement est complexe, on est amené à compliquer
les modèles élémentaires (figurés sur la Figure 12) en les associant (selon les possibilités de la Figure 13).
piston
avec ou sans
dispersion axiale
RPA
Figure 12 : Motifs élémentaires.

série parallèle

échange

court-circuit recyclage

Figure 13 : Modes de liaison.

On dit qu'une association de réacteurs idéaux (par exemple celle de la Figure 14) est représentative d'un réacteur
réel, et on l'appelle modèle, si les courbes de répartition des temps de séjour des molécules qui les traversent sont
identiques. La courbe relative à l'association est calculée, et la courbe relative au réacteur est obtenue
expérimentalement en utilisant la technique des traceurs. Ce modèle n'est pas toujours unique, notamment pour les
écoulements complexes.

Figure 14 : Exemple de représentation d'un réacteur réel comme l'association de réacteurs idéaux.

représentation des défauts d'écoulement

La Figure 15 présente le modèle de Cholette et Cloutier. Le paramètre (1-)Q
 permet de rendre compte des courts-circuits, et le paramètre  des
zones stagnantes (recirculations). V
La fonction DTS associée est : Q
2   t 
Et   1      t    exp   
  ts  t 
 s  (1-)V
La Figure 16 représente la DTS pour des valeurs de  et  données.
Figure 15 : Modèle de Cholette et Cloutier.
E

0,5

1 RPA

0
0 1 2 3 4 5
t /
Figure 16 : Modèle de Cholette et Cloutier ( = 0,8 et  = 0,5).

1.5. CONCLUSIONS

Les mesures de distribution des temps de séjour permettent de détecter les circulations défectueuses. Dans les
installations industrielles, ces défauts peuvent apparaître lors de réutilisations d'équipements existants, ou bien
lorsque la capacité de production varie, ou encore quand la conception a été imprécise. Lorsque l'on a détecté un
défaut d'écoulement sur un réacteur industriel, on n'effectue pas forcément d'amélioration, car le gain de
productivité escompté est parfois trop faible pour justifier une modification coûteuse de l'installation : mais il est
important de connaître ces défauts et d’en tenir compte dans la marche de l’appareillage.
L'association adéquate de réacteurs idéaux avec les modes de liaison appropriés permet de modéliser l'écoulement
réel dans un système. Ce modèle pouvant ensuite être utilisé par exemple pour calculer la conversion chimique
dans un réacteur.

Lorsque l'on est confronté à un système trop complexe pour être immédiatement "calculé" à l'aide d'une approche
analytique classique, on doit définir une approche quantitative plus globale. L'approche systémique s'appuie sur la
relation entrée/sortie et sur la réponse du système à des stimuli extérieurs. Elle est complémentaire de l'approche
analytique détaillée et fournit un autre type d'information : globale et souvent très efficace. La Distribution des
Temps de Séjour est l'exemple le plus célèbre d'approche systémique, car elle fut extrêmement féconde.

 limites
Cette approche, souvent très efficace, a cependant des limites. En particulier du fait des conditions d'application
(régime permanent, fluide incompressible), toutefois la généralisation est souvent possible.

La courbe de répartition des temps de séjour d'une association de réacteurs est indépendante de l'ordre dans lequel
on les place, comme symbolisé sur la Figure 17. En chimie, cette indépendance est limitée au cas des réactions du
premier ordre (les plus fréquentes toutefois). Pour les réactions d'un ordre différent, il convient d'étudier en outre
les macro- et micromélanges.

Figure 17 : Deux configurations donnant la même DTS.