Escoamento em tubos 6 Defina, diferencie e explique as características de escoamentos laminares e turbulentos considerando o escoamento de um fluido numa tubulação de área transversal circular. Discuta a respeito das tensões viscosas em escoamentos laminares e turbulentos. Em qual destes escoamentos tem-se maiores tensões? Por quê? Explique como elas atuam para determinar os perfis de velocidade de um fluido escoando em uma tubulação de seção reta circular em regime permanente.
Passo 1
O que o enunciado pede e as informações a usar:
Pedem: (i) definir e diferenciar escoamentos laminar e turbulento em uma tubulação circular; (ii) discutir as tensões (tensões viscosas) em cada caso e dizer em qual há maiores tensões, com justificativa; (iii) explicar como essas tensões atuam para determinar os perfis de velocidade em regime permanente.
Informação física básica que vamos usar: fluxo plenamente desenvolvido, regime permanente, seção transversal circular de raio R (diâmetro D=2R).
Objetivo da solução: dar definições claras, mostrar as fórmulas relevantes (Reynolds, leis de balanço de momento), derivar o perfil de velocidade no caso laminar (Hagen–Poiseuille) e explicar, com a mesma equação de balanço, como surgem as tensões em escoamento turbulento (introduzindo a tensão de Reynolds e a lei do muro).
Passo 2
Definição e critério quantitativo (Reynolds): o que é a quantidade que controla a transição.
A grandeza que costuma indicar o tipo de escoamento é o número de Reynolds
\mathrm{Re}=\dfrac{\rho U D}{\mu}=\dfrac{U D}{\nu},
onde \rho é a massa específica do fluido, U é a velocidade média do escoamento na seção, D é o diâmetro da tubulação, \mu é a viscosidade dinâmica e \nu=\mu/\rho a viscosidade cinemática.
Regiões típicas (valores orientativos para tubo circular):
Caracterização física do escoamento laminar (descrição conceitual).
Em escoamento laminar as partículas do fluido movem-se em camadas quase-paralelas, sem mistura macroscópica entre elas. O transporte de quantidade de movimento entre camadas é dominado pela viscosidade molecular.
Consequência prática: o perfil de velocidade é suave e determinístico (para tubo circular, simetricamente parabólico no caso plenamente desenvolvido e newtoniano). A dissipação de energia é menor que em fluxo turbulento, para a mesma velocidade média, até a estabilidade do caso laminar.
Passo 4
Caracterização física do escoamento turbulento (descrição conceitual).
Em escoamento turbulento há flutuações caóticas sobrepostas à velocidade média; ocorrendo mistura turbulenta que transporta quantidade de movimento de forma muito mais eficiente.
Consequência prática: o perfil de velocidade é “mais plano” no núcleo e apresenta um gradiente muito acentuado próximo à parede; a perda de carga (e o atrito) são muito maiores do que no caso laminar para a mesma velocidade média.
Passo 5
O que é tensão viscosa num fluido newtoniano e como ela aparece no escoamento tubular.
Para um fluido newtoniano, a tensão de cisalhamento viscosa local relaciona-se com o gradiente de velocidade. Num escoamento unidimensional (componente axial x e coordenada radial r), a componente relevante é \tau_{r x} e, localmente,
onde u(r) é a velocidade média axial como função do raio r (sinal e orientação devem ser tratados conforme convenção). Essa é a contribuição viscosa molecular para o transporte de momento entre camadas.
Passo 6
Balanço de momento em um tubo e expressão geral da tensão total em fluxo plenamente desenvolvido. (Derivação curta e didática.)
Considere um cilindro de fluido de raio r e comprimento L. Para regime permanente plenamente desenvolvido, o balanço de forças na direção x (axial) entre pressão e forças de cisalhamento leva a
\text{força de pressão} + \text{força de cisalhamento} =0.
A força de pressão sobre a seção transversal de raio r é -\dfrac{\mathrm{d}p}{\mathrm{d}x}\,\pi r^{2}\,L (onde \mathrm{d}p/\mathrm{d}x é constante no trecho). A força de cisalhamento na superfície lateral vale \tau_{r x}(r)\,2\pi r\,L. Igualando (em módulo e cuidando dos sinais) obtém-se
Observações: \mathrm{d}p/\mathrm{d}x<0 para escoamento no sentido positivo, então a tensão varia linearmente com r e atinge seu máximo em magnitude na parede r=R: \tau_{w}=\tau_{r x}(R)= -\dfrac{R}{2}\dfrac{\mathrm{d}p}{\mathrm{d}x}.
Passo 7
Perfil de velocidade no caso laminar (Hagen–Poiseuille): derivação completa e interpretação.
No laminar plenamente desenvolvido o transporte de momento é puramente viscoso: \tau_{r x}^{(\text{visc})}=\mu\dfrac{\mathrm{d}u}{\mathrm{d}r}. Igualando à expressão do balanço de momento do passo anterior:
onde C é constante de integração. Pela condição de não deslizamento na parede u(R)=0, obtemos
0 = -\dfrac{1}{4\mu}\dfrac{\mathrm{d}p}{\mathrm{d}x}\,R^{2} + C \quad\Rightarrow\quad C = \dfrac{1}{4\mu}\left(-\dfrac{\mathrm{d}p}{\mathrm{d}x}\right)R^{2}.
Observação prática: no laminar a tensão viscosa local é responsável por todo o transporte de momento; ela cresce linearmente com r e o perfil de velocidade sai parabolicamente.
Passo 8
Como o turbulento aparece nas equações médias: decomposição de Reynolds e tensão de Reynolds.
Em fluxo turbulento introduz-se a decomposição de Reynolds u = U + u', v = V + v' (onde U é a média axial, u' a flutuação, e v componente radial com média V=0 em plenamente desenvolvido). Ao fazer a média das equações de Navier–Stokes aparece um termo extra de transporte de momento por flutuações, a chamada tensão de Reynolds:
Interpretação: a distribuição total de esforço que equilibra a pressão é a mesma (linha em r), porém em regime turbulento a contribuição viscosa \mu \mathrm{d}U/\mathrm{d}r é pequena exceto muito próxima à parede; a maior parte do transporte de momento é feita por -\rho\overline{u'v'} (as “tensões turbulentas” ou tensões de Reynolds).
Passo 9
Perfil de velocidade em turbulência: conceitos úteis e a “lei do muro”.
Por causa da forte contribuição de -\rho\overline{u'v'}, o perfil médio U(r) em turbulência é muito mais plano no núcleo e com um gradiente grande numa camada muito fina próxima à parede. Essa estrutura pode ser descrita por escalas de parede introduzidas a seguir.
Defina velocidade de fricção
u_{*}=\sqrt{\dfrac{\tau_{w}}{\rho}},
onde \tau_{w}=\tau_{r x}(R) é a tensão na parede. As grandezas adimensionais clássicas são
já na região de sobreposição (camada log) a lei empírica é a lei logarítmica:
U^{+}=\dfrac{1}{\kappa}\ln y^{+} + B,
onde \kappa\approx 0{,}41 (constante de von Kármán) e B\approx 5{,}0 (constante empírica). Estas leis mostram explicitamente: perto da parede a viscosidade fixa o gradiente; fora dela a mistura turbulenta domina e o perfil é mais plano.
Passo 10
Comparação direta: em qual escoamento as tensões são maiores e por quê? (resumo físico e fórmula prática de atrito)
Para uma mesma velocidade média U, o escoamento turbulento exige uma diferença de pressão muito maior (portanto maior \tau_{w}) do que o escoamento laminar. Em outras palavras, as tensões efetivas (especialmente a tensão na parede) são maiores em turbulência porque as flutuações transportam momento de modo muito eficaz — isso aumenta as taxas de dissipação e a perda de carga.
Fórmula prática de atrito:
Laminar: fator de atrito de Darcy (para escoamento interno, superfície lisa) é
f=\dfrac{64}{\mathrm{Re}}.
Turbulento (empírico, regime hidráulico liso): por exemplo a fórmula de Blasius (válida aproximadamente para 4\times 10^{3}\lesssim\mathrm{Re}\lesssim10^{5}):
f\approx 0{,}3164\,\mathrm{Re}^{-0{,}25}.Para a mesma \mathrm{Re} (acima da transição), f turbulento é tipicamente maior que o valor que o escoamento laminar teria (lembrando que, de fato, o laminar só existe estável até \mathrm{Re}\approx2{,}300). Portanto a tensão na parede \tau_{w}=\dfrac{f}{8}\rho U^{2} é maior em turbulência.
Passo 11
Como as tensões (viscosa + turbulenta) determinam o perfil de velocidade — resposta final integrada e prática para fazer na prova.
Procedimento prático para “achar” o perfil médio:
Escrever o balanço médio de momento: a tensão total deve satisfazer \tau_{r x}^{(\text{tot})}(r) = -\dfrac{r}{2}\dfrac{\mathrm{d}p}{\mathrm{d}x}. Isso é um resultado geométrico/geométrico do equilíbrio de forças.
No regime laminar: substituir \tau_{r x}^{(\text{tot})}=\mu\,\mathrm{d}u/\mathrm{d}r e integrar para obter a parabólica de Hagen–Poiseuille. (Feito no Passo 7.)
No regime turbulento: \tau_{r x}^{(\text{tot})}=\mu\,\mathrm{d}U/\mathrm{d}r -\rho\,\overline{u'v'}. Como -\rho\overline{u'v'} não é conhecido a priori, usa-se modelos (viscosidade turbulenta/hipótese de Boussinesq) ou leis empíricas (lei do muro) para relacioná-lo ao gradiente médio e integrar de forma aproximada. O resultado qualitativo é um perfil plano no núcleo e uma camada de grande gradiente junto à parede.
Dica prática para prova: derive sempre a expressão de \tau_{r x}(r) do balanço de momento (Passo 6); a partir daí, se for laminar, use \tau=\mu\,\mathrm{d}u/\mathrm{d}r e integre. Se for turbulento, explique a presença da tensão de Reynolds e use as escalas de parede e a lei logarítmica para descrever o comportamento próximo à parede e no núcleo.
Passo 12
Resumo final (motivacional e prático) — o que recordar para a prova:
Número de Reynolds para decidir o regime; balanço de momento para obter distribuição total de tensão; laminar → integração direta → perfil parabólico (Poiseuille); turbulento → decomposição de Reynolds, tensão de Reynolds domina fora da subcamada viscosa, uso da lei do muro/log-law para descrever o perfil médio; tensões (em particular \tau_{w}) maiores em turbulência, explicação física pela mistura e transporte de momento pelas flutuações.
Resposta
Escoamento laminar: fluxo em camadas ordenadas, sem mistura macroscópica; ocorre tipicamente para \mathrm{Re}\lesssim2{,}300; a tensão viscosa local é \tau_{r x}=\mu\,\dfrac{\mathrm{d}u}{\mathrm{d}r} e, pelo balanço de momento, \tau_{r x}(r)=-\dfrac{r}{2}\dfrac{\mathrm{d}p}{\mathrm{d}x}; integrando obtém-se o perfil parabólico
com vazão Q=\dfrac{\pi R^{4}}{8\mu}\left(-\dfrac{\mathrm{d}p}{\mathrm{d}x}\right).
Escoamento turbulento: existe flutuação e mistura turbulenta; típico para \mathrm{Re}\gtrsim4{,}000; o termo extra é a tensão de Reynolds \tau_{r x}^{(\text{Re})}=-\rho\,\overline{u'v'} e a tensão total média satisfaz
Em turbulência a contribuição de -\rho\overline{u'v'} domina fora da subcamada viscosa, produzindo um perfil médio mais “plano” no núcleo e um gradiente forte próximo à parede; nas proximidades da parede usa-se escalas de parede U^{+}=U/u_{*}, y^{+}=y u_{*}/\nu, com U^{+}\approx y^{+} (subcamada viscosa) e U^{+}=(1/\kappa)\ln y^{+}+B (camada log).
Tensões maiores: o escoamento turbulento apresenta maiores tensões efetivas (particularmente maior \tau_{w} para a mesma velocidade média) porque as flutuações turbulentas promovem transporte de momento muito mais eficaz do que a viscosidade molecular, aumentando a perda de carga; empiricamente, o fator de atrito evidencia isso (f=64/\mathrm{Re} para laminar versus f\sim0{,}3164\,\mathrm{Re}^{-0{,}25} de Blasius para turbulência em tubo liso em faixa adequada).