In [4]:
from pint import UnitRegistry
u=UnitRegistry()
Método de Coe e Clevenger¶
Qual o diâmetro de um decantador do tipo Dorr para tratar 4800 $lb_{\text{sólido}}/h$ de uma suspensão com 5 $lb_{\text{água}}/lb_{\text{sólido}}$ de sólidos. Pretende-se uma lama com a concentração de 1,5 $lb_{\text{água}}/lb_{\text{sólido}}$.
$\rho_{\text{água}}=\ 62,3\ lb/ft^3$
Os testes em descontínuo deram os seguintes resultados:
In [5]:
Q=4800*u.lb/u.hr
CA=5*u.lb/u.lb
CE=1.5*u.lb/u.lb
C=[5.0,4.2,3.7,3.1,2.5]*u.lb/u.lb
U=[2.0,1.2,0.94,0.7,0.5]*u.ft/u.hr
rho=62.3*u.lb/u.ft^3
print 'C=',C,'\n','u=',U
In [6]:
A=[]
c=[]
for i in range(len(C)):
A.append((Q*(C[i]/rho-CE/rho)/U[i]).to(u.m^2).magnitude)
c.append(C[i].magnitude)
In [10]:
dados=zip(c,A)
G=scatter_plot(dados,figsize=4)
S=spline(dados)
G+=plot(S,min(c),max(c),color='gray',linestyle='--')
x=[S[i][1] for i in range(len(S))]
x=[S((c[-1]-c[0])/100*i+c[0]) for i in range(100)]
ymax=max(x)
xmax=x.index(ymax)/100*(c[-1]-c[0])+c[0]
G+=arrow((xmax*1.1,ymax),(min(c),ymax),arrowsize=2,linestyle='--',
color='red')
G.show(gridlines="Major",
title=u'Área do decantador é de %.0f $m^2$'%ymax,
axes_labels=(u'$c\ (lb_{água}/lb_{sólido})$','$A\ (m^2)$'),frame=True,
ymax=1.1*ymax)
print(u'Factor de segurança 2\nDiâmetro do decantador é de %.1f m'%((2*ymax*4/pi())^0.5))
