Mecânica dos fluidos/Exercícios resolvidos/E4

• Teoria
• Voltar para a lista de exercícios resolvidos

Enunciado

Deduzir o perfil de velocidades para escoamento laminar de um líquido Newtoniano em um tubo cilíndrico horizontal de diâmetro D.

Solução

Aqui, não podemos modelar as paredes do tubo como duas placas paralelas, e será preciso empregar as equações em coordenadas cilíndricas.

$\frac{1}{r}\frac{\partial }{\partial r}(r{v}_r)+\frac{1}{r}\frac{\partial }{\partial \theta }(v_{\theta })+\frac{\partial }{\partial z}(v_z)=0$

${\rho }_0(\frac{\partial v_r}{\partial t}+v_r\frac{\partial v_r}{\partial r}+\frac{v_{\theta }}{r}\frac{\partial v_r}{\partial \theta }-\frac{v_{\theta }^2}{r}+v_r\frac{\partial v_r}{\partial z})=$

$={\rho }_0{g}_r-\frac{\partial p}{\partial r}+{\mu }_0[\frac{\partial }{\partial r}(\frac{1}{r}\frac{\partial }{\partial r}(r{v}_r))+\frac{1}{r^2}\frac{{\partial }^2{v}_r}{\partial {\theta }^2}-\frac{2}{r^2}\frac{\partial v_{\theta }}{\partial \theta }+\frac{{\partial }^2{v}_r}{\partial z^2}]$

${\rho }_0(\frac{\partial v_{\theta }}{\partial t}+v_r\frac{\partial v_{\theta }}{\partial r}+\frac{v_{\theta }}{r}\frac{\partial v_{\theta }}{\partial \theta }+\frac{v_r{v}_{\theta }}{r}+v_r\frac{\partial v_{\theta }}{\partial z})=$

$={\rho }_0{g}_{\theta }-\frac{1}{r}\frac{\partial p}{\partial \theta }+{\mu }_0[\frac{\partial }{\partial r}(\frac{1}{r}\frac{\partial }{\partial r}(r{v}_{\theta }))+\frac{1}{r^2}\frac{{\partial }^2{v}_{\theta }}{\partial {\theta }^2}+\frac{2}{r^2}\frac{\partial v_{\theta }}{\partial \theta }+\frac{{\partial }^2{v}_{\theta }}{\partial z^2}]$

${\rho }_0(\frac{\partial v_z}{\partial t}+v_r\frac{\partial v_z}{\partial r}+\frac{v_{\theta }}{r}\frac{\partial v_z}{\partial \theta }+v_z\frac{\partial v_z}{\partial z})=$

$={\rho }_0{g}_z-\frac{\partial p}{\partial z}+{\mu }_0[\frac{1}{r}\frac{\partial }{\partial r}\left(r\frac{\partial v_z}{\partial r}\right)+\frac{1}{r^2}\frac{{\partial }^2{v}_z}{\partial {\theta }^2}+\frac{{\partial }^2{v}_z}{\partial z^2}]$

Por simetria, para qualquer propriedade η, $\frac{\partial \eta }{\partial \theta }=0$. O movimento é unidimensional, portanto $v_r=v_{\theta }=0$. Em regime estacionário, $\frac{\partial \eta }{\partial t}=0$. Assim, as equações se tornam

$0+0+\frac{\partial }{\partial z}(v_z)=0\Rightarrow \frac{\partial v_z}{\partial z}=0$

${\rho }_0(0+0-0+0)={\rho }_0{g}_r-\frac{\partial p}{\partial r}+{\mu }_0[0+0-0+0]\Rightarrow \frac{\partial p}{\partial r}={\rho }_0{g}_r$

${\rho }_0(0+0+0+0+0)={\rho }_0{g}_{\theta }-\frac{1}{r}\frac{\partial p}{\partial \theta }+{\mu }_0[0+0+0+0]\Rightarrow \frac{1}{r}\frac{\partial p}{\partial \theta }={\rho }_0{g}_{\theta }$

${\rho }_0(0+0+0+v_z\frac{\partial v_z}{\partial z})=0-\frac{\partial p}{\partial z}+{\mu }_0[\frac{1}{r}\frac{\partial }{\partial r}\left(r\frac{\partial v_z}{\partial r}\right)+0+\frac{{\partial }^2{v}_z}{\partial z^2}]$

$\frac{\partial v_z}{\partial z}=0\Rightarrow \frac{{\mu }_0}{r}\frac{\partial }{\partial r}\left(r\frac{\partial v_z}{\partial r}\right)=\frac{\partial p}{\partial z}$

As três primeiras equações nos trazem informação trivial: o fluxo deve ser constante ao longo do eixo Z (horizontal) e a distribuição de pressões ao longo de cada seção circular é hidrostática. Integrando a última equação, teremos

$\frac{\partial }{\partial r}\left(r\frac{\partial v_z}{\partial r}\right)=\frac{1}{{\mu }_0}\frac{\mathrm{\Delta }p}{L}r\Rightarrow r\frac{\partial v_z}{\partial r}=k_1{r}^2+k_2(k_1=\frac{1}{2{\mu }_0}\frac{\mathrm{\Delta }p}{L})$

$\frac{\partial v_z}{\partial r}=k_{1}r+k_2\frac{1}{r}$

mas k₂ deve ser nulo, caso contrário a velocidade seria infinita em r = 0 (centro do tubo). Assim,

$\frac{\partial v_z}{\partial r}=k_{1}r\Rightarrow v_z=\frac{k_1}{2}{r}^2+k_3$

Como v_z(R) = 0, temos

$v_z(R)=0\Rightarrow \frac{k_1}{2}{R}^2+k_3=0\Rightarrow k_3=-\frac{k_1}{2}{R}^2$

Assim

$v_z=\frac{k_1}{2}(r^2-R^2)=\frac{1}{4{\mu }_0}\frac{\mathrm{\Delta }p}{L}(r^2-R^2)$

Finalmente, é bom lembrar que, para saber se o escoamento é laminar, é preciso examinar o número de Reynolds, que aqui tem a forma:

$N_{Re}=\frac{\rho \overline{v}D}{{\mu }_0}$

Por isso, é importante determinar a velocidade média. Para isso, calcula-se primeiro a vazão na seção circular A

$\mathrm{\Phi }={\int }_A{v}_{z}dA={\int }_0^R\frac{1}{4{\mu }_0}\frac{\mathrm{\Delta }p}{L}(r^2-R^2)(2\pi rdr)=\frac{\pi }{2{\mu }_0}\frac{\mathrm{\Delta }p}{L}{\int }_0^R(r^3-R^{2}r)dr$

$\mathrm{\Phi }=\frac{\pi }{2{\mu }_0}\frac{\mathrm{\Delta }p}{L}{(\frac{r^4}{4}-\frac{R^2{r}^2}{2}))|}_0^R$

$\mathrm{\Phi }=-\frac{\pi R^4}{8{\mu }_0}\frac{\mathrm{\Delta }p}{L}$

$\overline{v}=\frac{Phi}{A}=-\frac{\pi R^4}{8{\mu }_0}\frac{\mathrm{\Delta }p}{L}\frac{1}{\pi R^2}=-\frac{R^2}{8{\mu }_0}\frac{\mathrm{\Delta }p}{L}$

Os valores negativos de vazão e velocidade média devem-se ao fato de o fluxo ocorrer do ponto de maior pressão para o ponto de menor pressão (Δp < 0). Predefinição:AutoCat