In [23]:
from pint import UnitRegistry
u=UnitRegistry()
Método de Coe e Clevenger¶
Qual o diâmetro de um decantador do tipo Dorr para tratar 50 $m^3/h$ de uma suspensão com 236 g/L de sólidos. Pretende-se uma lama com a concentração de 550 g/L.
Os testes em descontínuo deram os seguintes resultados:
In [24]:
QA=50*u.m^3/u.hr
CA=0.236*u.t/u.m^3
CE=0.550*u.t/u.m^3
C=[265,285,325,415,465,550]*u.g/u.L
U=[20,5,1.5,1.1,1,0.9]*u.cm/u.hr
print 'C=',C,'\n','u=',U
In [25]:
A=[]
c=[]
for i in range(len(C)):
A.append((QA*CA*(C[i]^-1-CE^-1)/U[i]).to(u.m^2).magnitude)
c.append(C[i].magnitude)
In [26]:
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=('$c\ (g/L)$','$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))
