dt=60.0; #ID of tube in cm
dp=300; #Size of particles of bed in micrometers
umf=3; #Velocity at minimum fluidization condition in cm/s
ephsilonmf=0.5; #Void fraction of bed at minimum fluidization condition
db=5; #Diameter of bubble in cm
g=980; #Acceleration due to gravity in cm/s**2
#CALCULATION
#Computation of rise velocity of bubble
if (db/dt)<0.125:
ubr=(0.711*((g*db)**0.5));#Rise velocity by Eqn.(3)
elif (db/dt)<0.6:
ubr=(0.711*((g*db)**0.5))*1.2*exp(-1.49*(db/dt));#Rise velocity by Eqn.(4)
#Computation of cloud thickness
Rb=db/2.0; #Radius of bubble
uf=umf/ephsilonmf; #Velocity of emulsion gas
Rc=Rb*((ubr+(2*uf))/(ubr-uf))**(1/3.0);#Radius of cloud by Eqn.(6)
#OUTPUT
print 'The rise velocity of the bubble=%.1f cm/s'%ubr
print 'The cloud thickness=%.2f cm'%(Rc-Rb)
print 'From Fig.8(page 124)comparing fw vs dp, for dp = %.0f micrometer, wake fraction = 0.24'%dp
import math
#Variable declaration
uo=15; #Superificial gas velocity in cm/s
umf=1; #Velocity at minimum fluidization condition in cm/s
lor=2.0; #Pitch of perforated plate in cm
g=980; #Acceleration due to gravity in cm/s**2
#CALCULATION
#Case(a) For porous plate
dbo1=(2.78/g)*(uo-umf)**2;#Initial bubble size using Eqn.(19)
#Case(b) For Perforated plate
Nor=(2/math.sqrt(3))*(1/lor)**2;#Number of orifices in cm**-2
dbo2=(1.30/(g**0.2))*((uo-umf)/Nor)**0.4;#Initial bubble size using Eqn.(15) assuming inital bubble size is smaller than hole spacing
#OUTPUT
print 'Case(a) For porous plate'
print '\tInitial bubble size=%.2fcm'%dbo1
print 'Case(b) For Perforated plate'
print '\tInitial bubble size=%.2fcm'%dbo2
print '\tSince %f<%f, the equation used is correct.'%(dbo2,lor)