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}\)