Here it is an idea for having a very first approx of the solution that you can use for more precise calculations:
Hypotesis:
1)water like system (cp,k,d water like)
2)simple Michaelis Menten cinetic enzyme-substrate expression
3)CSTR reactor (well mixed)
4)fixed U derived from tables: 200-500 w/(m2°C)
known quantities:
Feed to the reactor
Exchange surface A
T_cin (cooling temperature in)
T_r (reaction temperature)
Volume V (wet surface)
DeltaH (heat of reaction)
HEAT TO REMOVE:
E=DeltaH*r=W*cp*(T_out-T_in) (1)
E=DeltaH*r=A*U*DeltaT_log (2)
DeltaT_log=
=[(T_r-T_out)-(T_r-T_in)]/ln[(T_r-T_out)/(T_r-T_in)] (3)
cp: water specific heat
W: cooling water capacity
r (mol/time): rate of reaction in Michaelis Menten form:
r=(v_max*cs)/(Km*cs) (4)
where cs is the glucose concentration, Km is the M-M costant and v_max is the maximum rate of reaction
r is evaluated at stationary state: r(cs_ss)
cs_ss is found numerically letting the material balance go to zero:
dcs/dt=Q/V*(cs_in-cs_ss)-(v_max*cs_ss)/(Km*cs_ss)=0 (5)
Km and v_max can be found here for the alcohol fermentation:
Q is the volumetric feed
(they're at 20°C, but i couldn't find anything better)
once you have cs_ss you can evaluate r(cs_ss) and evaluate DeltaH*r(cs_ss) (heat generated) (6)
now we can find an approx solution for W and T_out:
- fix W and evaluate T_out in eq. (1)
- evaluate DeltaT_log in eq. (3)
- evaluate E=A*U*DeltaT_log
(with U fixed,for example at 350)
- check if the value |(E-DeltaH*r(cs_ss))/E| is less then a value you chose (for example 0.1)
- if not, choose another W and try again and so on..
i hope it helped somewhat