20
8.3.20
12Dialog
{FEDBAT}
{Fermentation with batch start up}
{Flow rate is initially zero and is turned on at time=tfeed.}
{ Constants}
UM=0.3; 1/h
KS=0.1; kg/m3
K1=0.03;gP/kgX h
K2=0.08;gP/kgX
Y=0.8; kg X/kg S
X0=0.01; Initial biomass inoculum, kg/m3
S0=10; Initial substrate conc., kg/m3
P0=0; Initial product conc.,g/m3
SF=10; Feed conc. ,kg/m3
F1=1.5; Feed flow rate, m3/h
tfeed=22.5; Start time for the feed
{Initial Conditions}
init V=1
init VX=V*X0
init VS=V*S0
init VP=V*P0
{Mass balances, kg/h}
d/dt(V)=F
d/dt(VX)=RX*V
d/dt(VS)=F*SF+RS*V
d/dt(VP)=RP*V {g/h}
{Calculation of concentrations}
X=VX/V
S=VS/V
P=VP/V
{Kinetics}
RX=U*X
U=UM*S/(KS+S)
RS=-RX/Y
RP=(K1+K2*U)*X
D=F/V {nominal dilution rate, 1/h}
{Turning the feed on at time = tfeed}
F=if time=tfeed then F1 else 0 {batch start up}
12Dialog
rk4
STOPTIME 50.000000
F1
1.500000
0.000000
3.000000
0.100000
0
0
1.0
1
0
1
Graph 0
Graph
1
0
0
1
false
0
1
0
0
1
1
0
0
0
0.000000 0.000000
TIME
X, P
S
0 1
1 1
2 1
0.796875 0.545313
0.333333 0.500000
TIME
V
0
0
VX
0
1
VS
0
0
VP
0
1
X
1
0
S
1
1
P
1
0
RX
0
1