4.2 Development of a Mathematical Model
73
PROGRAMME 1 ANALYSIS FOR HYBRID BIOREACTOR
WRITE(*,*)'INPUT THE VALUES OF S0,BK,AK,THETA,P,BS,BT,BD, D,DF,AL,XF,A,Y'
READ(*,*)S0,BK,AK,THETA,P,BS,BT,BD, D,DF,AL,XF,A,Y
THETA=THETA/24
WRITE(*,*)'IF PURELY ATTACHED GROWTH, THEN IATGR=1'
WRITE(*,*)'IF BOTH SUSPENDED AND ATTACHED GROWTH, THEN IHYBRD=1'
WRITE(*,*)'INPUT IATGR,IHYBRD'
READ(*,*)IATGR,IHYBRD
EPSI=10.**(-3)
SMIN=(AK*BT)/(Y*BK-BT)
SW=SMIN+0.0001
10 SW=SW
WRITE(*,*)'SW=',SW
KOUNT=100000
DO I=1,KOUNT
CONTINUE
THETAC=P*THETA
TEMP3=BK*SW*THETA*Y*A*(BS/BT)
TEMP4=(AK+SW)*((1./THETAC)+BD-((Y*BK*SW)/(AK+SW)))
A1=TEMP3/TEMP4
IF (IATGR.EQ.1) THEN
A1=0
ELSE
END IF
S1=SMIN+0.2*(SW-SMIN)
TEMP1=2.*BK*XF*DF
TEMP21=AK*ALOG((AK+SMIN)/(AK+S1))
TEMP31=(TEMP1*(S1-SMIN+TEMP21))
AJ1=SQRT(TEMP31)
S2=SMIN+0.4*(SW-SMIN)
TEMP22=AK*ALOG((AK+S1)/(AK+S2))
TEMP32=(TEMP1*(S2-S1+TEMP22))+((AJ1)*(AJ1))
AJ2=SQRT(TEMP32)
S3=SMIN+0.6*(SW-SMIN)
TEMP23=AK*ALOG((AK+S2)/(AK+S3))
TEMP33=(TEMP1*(S3-S2+TEMP23))+((AJ2)*(AJ2))
AJ3=SQRT(TEMP33)
S4=SMIN+0.8*(SW-SMIN)
TEMP24=AK*ALOG((AK+S3)/(AK+S4))
TEMP34=(TEMP1*(S4-S3+TEMP24))+((AJ3)*(AJ3))
AJ4=SQRT(TEMP34)
S5=SW
TEMP25=AK*ALOG((AK+S4)/(AK+S5))
TEMP35=(TEMP1*(S5-S4+TEMP25))+((AJ4)*(AJ4))
AJ5=SQRT(TEMP35)
AJ0=0
AJAVG=(AJ0+AJ1+AJ2+AJ3+AJ4+AJ5)/6
SWCAL=S0-AJAVG*(A1+A*THETA)
WRITE(*,*)'SWCAL=',SWCAL,'AJAVG=',AJAVG
IF(ABS(SWCAL-SW).LT.EPSI)THEN
WRITE(*,*)'ITERATION CONVERGES'
WRITE(*,*)'SW=',SW,'AJAVG=',AJAVG
STOP 'NORMAL END'
ELSE
SW=SW+0.00001
GO TO 10
END IF
KOUNTC=MOD(I,100)
IF (KOUNTC.EQ.1)THEN
WRITE(*,*)'ITERATION DOES NOT CONVERGE'
WRITE(*,*)'I=',I,'SW=',SW
END IF
IF(I.GE.KOUNT)THEN
EXIT
END IF
END DO
STOP'ABNORMAL END'
END
Program 2 For finding out L e