Dimensionnement du système de refroidissement d'un bioréacteur
Technologies disponibles
Différentes technologies existent pour réguler la température d'un bioréacteur :
(a) : la double enveloppe
(b) : le serpentin externe
(c) : le serpentin hélicoïdal interne
(d) : la chicane interne de refroidissement
(f) : l'échangeur de chaleur externe
Source image : Bioprocess Engineering Principles, Pauline M.Doran, Academic Press, 2002
Notions générales sur le transfert thermique dans une cuve agitée
Bioréacteur équipé d'une double enveloppe, régulé thermiquement par un fluide caloporteur.
Cuve : vert, double enveloppe : orange
Tco : température de sortie du fluide caloporteur (°C)
Tci : température d'entrée du fluide caloporteur(°C)
\(\dot{M}_C\): flux de fluide caloporteur (kg.s-1)
Cpc: capacité calorifique du fluide caloporteur (J.kg-1.°C-1)
af : surface d'échange (m2)
U : coefficient global d'échange thermique(W.m-2.°C-1)
h1 : coefficient d'échange thermique à l'intérieur de la cuve (W.m-2.°C-1)
h2 : coefficient d'échange thermique au sein du fluide caloporteur (W.m-2.°C-1)
e: épaisseur de la paroi (m)
λ: conductivité thermique de la paroi (W.m-1.°C-1)
Cpb: capacité calorifique du milieu (J.kg-1.°C-1)
λb: conductivité thermique du milieu (W.m-1.°C-1)
Tb : température du milieu (°C)
Mb : quantité de milieu (kg)
μb : viscosité du milieu (Pa.s)
N : vitesse de rotation (s-1)
Ø : Flux de chaleur transféré (W)
Équation générale du bilan thermique dans le cas du refroidissement d'un bioréacteur
Pour avoir un état stationnaire au niveau de la température à l'intérieur de la cuve, il faudra que la chaleur produite par le bioréacteur ØT soit égale au flux transféré et au flux capté par le fluide caloporteur.
ØT=Flux transféré=Flux capté
\(Ø_T=U.\frac{T_{co}-T_{ci}}{ln\frac{T_b-T_{ci}}{T_b-T_{co}}}.a_f\)\(= \dot{M}_C.C_{pc}.(T_{co}-T_{ci})\)
Règle de calcul du coefficient global d'échange thermique U
\(\frac{1}{U}=\frac{1}{h_1}+\frac{e}{λ}+\frac{1}{h_2}\)
En milieu visqueux, comme c'est souvent le cas dans les bioréacteurs, h1<<h2, et par ailleurs h1<<e/λ, il sera donc possible d'approximer : \(\frac{1}{U}=\frac{1}{h_1}\)
Calcul de h1 :
Pour la calcul du coefficient d'échange thermique à l'intérieur de la cuve, on recherchera une corrélation empirique du Nusselt correspondant à la géométrie de la cuve, du type : \(Nu_1=k.{Re_A}^α.{Pr}^β\)
Avec :
- Nu1 : Nombre de Nusselt caractérisant l'échange thermique à l'intérieur de la cuve, \(Nu_1=\frac{h_1.D}{λ_b}\)
- ReA : Reynolds d'agitation, \(Re_A=\frac{ρ.N.d^2}{µ_b}\)
- Pr : Nombre de Prandlt, \(Pr=\frac{C_{pb}.µ_b}{λ_b}\)