REACTEURS NON IDÉAUX

Les modèles d'écoulement des réacteurs idéaux représentent les deux cas extrêmes de
mélange total et instantané (récipient uniforme) et absence totale de mélange (gradient à la
longueur du tube produite uniquement par la réaction chimique). Les réacteurs idéaux sont simples
de résoudre, et elles montrent souvent un contact entre les réactifs et un comportement
optimaux. Par exemple, le réservoir continu idéal est un système parfaitement agité avec
composición uniforme en el reactor y en el efluente [mezcla completa], cuyo balance de
la matière à l'état stationnaire est : =cox/v.
Dans de nombreux cas, une agitation efficace rapproche suffisamment ces
réacteurs à leur modèle idéal. Mais il y en a d'autres où, pour diverses raisons, l'écart
peut être significatif : grands réservoirs avec agitation insuffisante, réactions rapides
face au temps de mélange, canalisations, recirculation ou stagnation du fluide. Les
les problèmes d'écoulement non idéal sont étroitement liés aux changements d'échelle ; à
menudo, le facteur non contrôlé lors du passage d'un réacteur de laboratoire à un autre pilote ou
industriel est la magnitude de non-idéalité du flux.

TYPE DE FLUX DANS UN RÉACTEUR

Pour caractériser le type d'écoulement dans un réacteur réel, on étudie la réponse transitoire à un
changement brusque de concentration à l'entrée, en l'absence de réaction chimique. Pour effectuer le
On utilise des substances compatibles avec le fluide de mélange et facilement analysables.
(traceurs), et selon le type de perturbation, il existe deux courbes typiques de réponse :

- Courbe F : réponse à un échelon (changement soudain et soutenu de la concentration de

trazador à l'entrée), qui a tendance à se généraliser par F=c/coy =t/(V/Vo).

Pour un réacteur parfaitement agité, le bilan au traceur est : Vo(co-c)=V· c/ t, e intégrant

le résultat de l'exponentielle asymptotique : F = 1 - e- .

- Courbe C : réponse à une impulsion (changement brusque et instantané à l'entrée du

réacteur), qui se généralise dans ce cas avec co=m/V, siendomla cantidad de
trazador inyectado.

Mathématiquement, la fonction d'impulsion instantanée est la dérivée de la fonction échelon.

antérieur, donc si le mélange est linéaire, les fonctions de réponse seront également :
C= F/ = e- , qui correspond à une exponentielle décroissante vers l'axe des abscisses.

Pour un réacteur tubulaire idéal, les courbes seraient la reproduction exacte de la perturbation, avec le
retard du temps de résidence.

Le dispositif expérimental peut être une vanne à trois voies, qui se change de position à temps.
zéro, et un prélèvement périodique ou un détecteur en ligne pour déterminer la concentration
saliente avec le temps ; dans le deuxième cas, la valve est remplacée par une seringue avec le traceur.

solvants traceur F =1
Vo Vo,co perturbation
perturbation respuesta

RÉACTEUR
V 0
C
( perturbation
Vo réponse

réponse c(t)

0
PARTIE EXPÉRIMENTALE

Expérience en escalier.

On choisit un volume de réacteur (V), on introduit un traceur, on remplit de

eau et se homogénéise en sélectionnant une vitesse d'agitation. On alimente
un flux stationnaire d'eau (Vot=0) et la concentration du traceur est mesurée dans le
effluents, répéter l'expérience pour plusieurs débits calibrés. Pour obtenir F(t)
normalisée nous changeons l'échelle pour que la courbe monte de 0 à 1, ayant
en tenant compte que l'escalier a été réalisé à l'envers ; et en calculant les pentes
la courbe C(t) peut être obtenue :
F(t)= 1-(At-A ) / (A0-A )= (A0-At/(A0-A ) C(t)=dF(t)/dt

Expérience avec pulsation.

Des flux d'eau stationnaires sont fournis au réacteur (Vo), injectant à

t=0 un pulse de solution concentrée de traceur. Comme auparavant se
représente la concentration en fonction du temps ou sous forme de la propriété mesurée;
l'aire sous la courbe se calcule par la règle des trapèzes (Qt= t ½(At-1+At),
pour obtenir les fonctions normalisées :
C(t)=At/Q ; F(t)=Qt/Q

En utilisant des températures et des temps spatiaux similaires, on peut étudier l'effet de
changement d'échelle avec des réacteurs plus grands.
Distribution des temps de résidence.

Le temps de résidence est le temps qu'un élément de fluide reste dans un réacteur continu;
La distribution des temps de résidence (RTD) doit être connue pour la conception correcte du réacteur.
La RTD est donnée par la fonction E(t) normalisée, qui peut être obtenue par des expériences.
estímulo-réponse de traceurs, que ce soit par des signaux de type impulsion (courbe C) ou en échelon
(courbe F) : la relation entre elles est E=C=dF/dt. La courbe E fournit des résultats plus réalistes,
puisque la F, en intégrant les effets, peut cacher le comportement réel.
Para la normalización ( Edt=1) on divise la réponse instantanée par l'aire sous la courbe
ou la concentration maximale stationnaire dans votre cas : C=c/ c·dt F=c/cmax.

Les courbes de distribution des temps de résidence des molécules qui sortent ou sont à l'intérieur
del réacteur, sont respectivement données par les fonctions extérieure (E=C) ou intérieure (I=1-F); pour un réacteur
de mélange complet, celles-ci sont identiques et correspondent à la fonction exponentielle e- .

Le temps moyen de résidence peut être calculé à partir des deux fonctions : tm(obs)= tEdt = tdF ;
Les essais doivent être réalisés avec des temps similaires à ceux de l'opération du réacteur.

Pour les réacteurs de mélange idéaux, la RTD théorique s'avère être comme nous l'avons vu : E=e-t/tmy F=1-e-t/tm ;
pour les systèmes à densité constante en régime stationnaire, le temps de résidence moyen coïncide avec le
temps spatial : tm= =V/VoUn temps observé inférieur peut indiquer des canalisations (et des zones stagnantes);
s'il est supérieur, cela ne peut pas s'expliquer par le bilan matière (cela peut être attribué à des erreurs dans la mesure du flux)
y, avec une probabilité moindre, du volume du réacteur ou du caractère non inerte du traceur).
En ce qui concerne les courbes du traceur, une déviation au-dessus de F(t) par rapport à l'idéal peut
peut également être un indice de stagnation. Un écart au-dessus de E(t) reflète une sortie plus rapide
(plus d'éléments de fluide que le flux idéal pour un tmdeterminado) y por tanto una menor conversión; un
Le retard dans la partie initiale ascendante de la courbe indique l'écoulement du piston en série avec le mélange complet.
Un increment du débit d'alimentation a pour effet l'approche des profils au flux
de piston, et parfois l'approche des courbes réelles et idéales. Pour sa part, les ondulations
peuvent provenir des fluctuations du flux.

Mélange de fluides.

Il y a deux aspects liés au mélange des fluides qui concernent également l'écoulement non
idéal : le degré de ségrégation (entre macro-microfluide), et le temps de mélange (important avec
deux courants de réactifs, où il peut y avoir des zones riches en l'un ou l'autre). Le premier concerne
le mélange à l'échelle moléculaire, et le second au comportement global du flux à travers le réservoir.
Lorsque un macrofluide entre dans un RCTA, la concentration dans les agrégats individuels
ne change pas instantanément -comme cela se produit avec un microfluide- mais de manière similaire à un système
discontinu
1-x= cm/co= [c/co]disque E dt [c/co]disque=(1+cokt)-1

Voici le fondement de la méthode d'attribution : à partir de paires de valeurs {[c/co]disque;F}(t) se

représente la fonction [c/co]disquevs F; ainsi peut-on obtenir xasigà partir de la F(t) réelle ou idéale (sa différence en
plus ou moins dépend de la position des courbes respectives au-dessus ou en dessous.

Cette méthode permet de relier la cinétique (qui a à voir avec les réactions chimiques)
et les expériences avec des traceurs (qui concernent l'écoulement des fluides), mais cela ne sert que
pour les macrofluides ; lorsque cela n'explique pas les résultats expérimentaux (xobs) il ne reste que les
"modèles empiriques de flux" comme alternative.

Ceci repose sur l'hypothèse d'un mécanisme pour le processus de mélange qui conduise à des RTD
suffisamment compatibles avec les courbes expérimentales ; alors on peut simuler le comportement
du réacteur, introduisant l'effet de la réaction (cinétique chimique) dans les équations du modèle d'écoulement.
Réacteur tubulaire avec dispersion axiale

Ce modèle assimile le réacteur réel à un flux de piston, mais avec un mélange axial contrôlé par le
diffusion moléculaire n/ t= -D·A· c/ l ; les bilans de matière pour le processus de mélange du traceur dans un
élément de volume de réacteur tubulaire, permettent d'obtenir les courbes C et F de ce modèle en fonction de
module D/uL, qui peut également être estimé à partir du nombre de Re et Sc, et du facteur géométrique (dt/L).

Modèle discret
Ce modèle est basé sur l'équivalence du réacteur tubulaire idéal à une batterie d'infinités éléments
parfaitement agités du même volume : L(C) 1
e
 S·
L(Co) lim n
n 1 S· /n)

similitude qui se démontre avec la dynamique des deux systèmes, dont la fonction de transfert est identique (quotient de
transforms de Laplace de la réponse et perturbation), obtenue à partir de l'équation différentielle en régime transitoire,
et correspond à un retard (combinaison de systèmes de 1erordre de constante de temps V/Votemps de résidence).

Dans les figures, on peut observer comment, en augmentant le nombre de réacteurs en série, les courbes de réponse
se rapprochent de l'impulsion ou de l'escalier (réacteur tubulaire idéal), et de tout ce qui précède, on déduit que cela peut
emploi de la fonction de transfert d'un processus de mélange comme l'expression mathématique du type d'écoulement
dans un réacteur, de sorte que si deux systèmes ont la même fonction, leur comportement sera identique.

En observant les réponses C du premier modèle de réacteur tubulaire, nous voyons qu'à mesure que
diminue la dispersion axiale (Bo=uL/D 100) les courbes coïncident pratiquement avec la distribu-
tion de Gauss (si étant la variance 2 2/Bo=1/n); cela nous permet de traiter un réacteur réel comme une batterie
finita de réacteurs agités : n=Bo/2. Les deux modèles emploient un seul paramètre.
Modèles pour réacteurs agités
On utilise généralement des modèles à deux paramètres, comme celui de Cholette-Cloutier, qui se
basa en admettre des zones stagnantes (volume mort, Vm) et courts-circuits (derivation, V o2), ou celui de
Van-de-Vusse, basé sur la circulation du fluide avec agitateur axial dans un récipient cylindrique
(équivalent à un réacteur tubulaire avec recirculation a=R/Voy mélange axialBo 0); dans la figure se
représentent les courbes de réponse de ce dernier, pour un recyclage défini (a=5) et une série
de nombres de réservoirs (n). Connus les paramètres de ces modèles, il est assez simple
simuler le comportement du réacteur réel correspondant.

REACTEUR AGITÉ (Modèle de Cholette et Cloutier)

Vo2
Voco c V= VM+Vm
Vo=Vo1+Vo2 Vo, co Vo1 c´ c
VM
V c= (Vo1·c´+Vo2·c)/Vo

c´=fM(Vo1,VM)
Vm

REACTEUR AGITÉ (Modèle de Van de Vusse)

R
Voco c
5 Vo, co c´ V (n,Bo) c
1 4
V
G(s)=L(c)/L(co)= G´(s)/[(1+a) -a·G´(s)]
2 3
G´(s)=L(c)/L(c´)

a=R/Vo= 5
EJERCICIOS DE FLUJO NO IDEAL

On utilise un traceur afin de caractériser le comportement d'un réservoir de mélange.

continuo, par un flux d'eau Q= 30 l/min. Le réservoir dispose de plusieurs baffles où se
ont produit des incrustations au cours des années où l'équipe a été opérationnelle, de sorte que
il ne connaît pas exactement le volume effectif de celui-ci (V, dm 3); dans l'expérience, on injecte
un pouls mo(g) du traceur à l'entrée (t=0), suivie de la réponse dans le courant de sortie-
da (C, g/l) :
t (min) 1 2 3 4
C (µg/l) 0,223 0,050 0,011 0,0025

Vérifiez graphiquement que le réservoir se comporte comme un mélangeur parfait et déterminez son
volume effectif.
C(t) = mo/V·e-Q·t/V

C (µg/l)

0,25
0,223
0,2

0,15

0,1

0,05 0,05

0,011 0,0025
0
1 2 3 4
t (min)

Ln C = ln (mo/V) - (Q/V)·t

-1

-2

-3

-4

-5

-6

-7
1 2 3 4 t (min)
ln C -1,5006 -2.9957 -4,5099 -5.9915

Q/V = 1,497 min-1 V = (30 l/min)/(1,497 min-1) = 20,0 l

Le tableau montre la réponse du traceur à une impulsion d'entrée pour caractériser le flux d'un
réacteur où se déroule une réaction de 1er ordonnance dont la constante cinétique est 0,307 min-1.
L'objectif est d'estimer la conversion en tenant compte des différentes situations et modèles de
flux non idéal.

Tiempo (min) 0 5 10 15 20 25 30 35
Conc. (g/dm3) 0 3 5 5 4 2 1 0

Pour obtenir la courbe E normalisée, ainsi que le temps de résidence moyen et l'écart type, on
integra el área bajo la curva:
c· t = (3+5+5+4+2+1) = 100 E = c / c· t

tm t·E· t = 15 min 2 1/tm2· t2·E· t - 1 = 0,211

Flux de piston : = Co x/v xp= 1 - e-k = 0,99

Mélange complet : = Cox /v xm= 0,82

Réacteur réel : xun 1 - e-kt·E· t = 0,953

Flux dispersé dans le piston : D/uL = f( ) 20,12 graf (k =4,6): x = 0,96

Réservoirs en série : N = 1/ 2 4,76 x = 1-C/Co= (1+k /N)N 0,96

Conception de mélangeurs pour maximiser le rendement

Une réaction compétitive-consécutive de deuxième ordre se déroule dans un réacteur continu

type de réservoir agité industriel
A b R b S

oA=
Données : C 0,1 mol/dm3k1= 35 dm3/mol·s viscosité du mélange, = 1,08 csk (1,16·10-5ft2/s)
Longueur du récipient : Lf= 1,5 ft (diámetro del agitador: D= 0,5 ft)

Utiliser les corrélations de Paul-Treybal pour estimer la vitesse d'agitation nécessaire avec
objet de minimiser les effets de mélange et d'atteindre le rendement maximal attendu (Y).

temps de micromélange : = 10-5/CoAk1= 2,86·10-6s

velocidad turbulenta: u = [ 0,882·( ·Lf)3/4/ ]4/7= 12,5 pi/s

N = u / 0,45 D =17,7 rps

Estimer le rendement par rapport à l'attendu (Yexp), si la concentration initiale est CoA= 0,2 mol/l
et la vitesse de l'agitateurs N = 200 tr/min.

u = 0,45 ·D·N = 2,36 pi/s

= 0,882·( ·Lf)3/4/u7/4= 5,3·10-5s

k1· .CoA= 3,7·10-4 (A/Aexp 0,92)