import math
#initialisation of variables
w= 78.9 #lbf.ft**3 weight
d= 0.01 #in bore
u= 8.67*10**-9 #lbf/ hr ft**2 viscosity
h= 18. #ft
l= 10. #ft length
#CALCULATIONS
Q= math.pi*w*(d/12)**4*(h+l)/(l*128*u)
#RESULTS
print 'Flow rate = %.2e ft**3/hr'%(Q)
#initialisation of variables
x= 0.1 #ft
w= 62.4 #lbf/ft**3 fluid density
v1= 10. #ft/sec velocity
u= 2.4*10**-5 #lbf/ft viscosity
g= 32.2 #ft/sec**2
k= 4.91
#CALCULATIONS
s= k*x*(w*v1*x/(u*g))**-0.5
Tw= 0.332*w*v1**2*(u*g/(w*x*v1))**0.5/g
R= 0.332*6*Tw
#RESULTS
print 'Thickness = %.2e*ft'%(s)
print ' Shear stress = %.3f lbf/ft**2'%(Tw)
print ' Shear stress = %.3f lbf/ft'%(R)
#initialisation of variables
r = 1
r1 = 1
#CALCULATIONS
e1= r+r1
e2= r-r1
#RESULTS
print ' vorticity in a forced vortex = %.f*k'%(e1)
print ' vorticity in a free vortex = %.f'%(e2)