Como o regime permanente é aplicado em sistemas de tubulações?

Passo 1

O que o enunciado pergunta: “Como o regime permanente é aplicado em sistemas de tubulações?”
O pedido é conceitual e prático: explicar o que significa regime permanente (steady‑state) num sistema de tubulação e mostrar, passo a passo, como se usa esse conceito para calcular vazões, perdas de carga, e resolver redes de tubulação. Não foram fornecidos números no enunciado, então vou explicar as equações, as hipóteses, os procedimentos e ainda mostrar um exemplo numérico didático para deixar tudo concreto.
Informações dadas (implícitas no enunciado): nenhum dado numérico específico; pressupõe‑se escoamento incompressível típico (água), tubulações com perdas por atrito e acessórios, possibilidade de bombas/motores e redes em série/paralelo.
O que vou entregar: definição do regime permanente; equações essenciais (continuidade, energia com perdas); como calcular perdas por atrito (Darcy‑Weisbach), fator de atrito (Reynolds, Colebrook/Swamee‑Jain), perdas localizadas; procedimento para solução de um tronco simples e para redes (Hardy Cross e método nodal); e um exemplo numérico passo a passo com iterações até convergência.

Passo 2

Definição e hipóteses do regime permanente (steady‑state) em tubulações.
  • Em regime permanente, todas as grandezas que interessem ao problema (vazão \,Q\,, velocidade \,V\,, pressão/altura piezométrica \,H) são constantes no tempo em cada ponto do sistema. Ou seja, não há variação temporal: \partial/\partial t = 0.
  • Hipóteses comuns para água: fluido incompressível ( \rho constante), propriedades constantes (viscosidade \mu ou cinemática \nu constantes), escoamento unidimensional médio nas seções.
  • Consequência prática: podemos aplicar conservação de massa em forma estacionária e equação de energia/ Bernoulli com termos de perda de carga dependentes apenas da vazão.

Passo 3

Conservação de massa (continuidade) em regime permanente.
  • Para um único tubo de seção uniforme, a vazão volumétrica é constante ao longo do tubo:
Q = A \, V
onde \,Q é a vazão volumétrica ( \text{m}^3/\text{s}), \,A é a área da seção transversal ( \text{m}^2) e \,V é a velocidade média ( \text{m/s}).
  • Se houver junções, a continuidade implica soma algébrica das vazões: por exemplo, numa junção 1 onde entram Q_1,Q_2 e sai Q_3, vale \,Q_1+Q_2-Q_3=0.

Passo 4

Equação da energia em regime permanente (forma de cabeça/piezométrica).
  • Entre dois pontos 1 e 2 no escoamento, a versão prática da equação de Bernoulli incluindo perdas e bombas é:
H_1 = H_2 + h_f + h_m - h_{bomba}
onde:
  • \,H é a cabeça piezométrica ( \text{m}), \,H = z + \dfrac{p}{\rho g} (altura geométrica z mais pressão convertida em altura);
  • \,h_f é a perda por atrito distribuída (major losses) ( \text{m});
  • \,h_m são as perdas localizadas (válvulas, cotovelos, bocais) ( \text{m});
  • \,h_{bomba} é a elevação fornecida por bomba (se houver) ( \text{m}). Sinal positivo quando a bomba adiciona energia.
  • Em problemas práticos, isolamos \Delta H = H_1-H_2 e equacionamos as perdas em função de Q.

Passo 5

Perdas por atrito (Darcy‑Weisbach) — fórmula e significado das variáveis.
  • A perda por atrito distribuída num comprimento \,L de tubo é dada por:
h_f = f \,\frac{L}{D}\,\frac{V^2}{2g}
onde:
  • \,f é o fator de atrito de Darcy (adimensional);
  • \,L é o comprimento do tubo ( \text{m});
  • \,D é o diâmetro interno do tubo ( \text{m});
  • \,V é a velocidade média ( \text{m/s});
  • \,g é a aceleração da gravidade ( \approx 9{,}81\,\text{m/s}^2).
  • Substituindo \,V = Q/A e \,A = \pi D^2/4 obtemos perda expressa diretamente em função de \,Q:
h_f = f \,\frac{L}{D}\,\frac{1}{2g}\left(\frac{Q}{A}\right)^2.

Passo 6

Fator de atrito \,f: dependência de Reynolds e rugosidade.
  • O fator \,f depende do regime hidráulico, caracterizado pelo número de Reynolds:
Re = \frac{V D}{\nu}
onde \,\nu é a viscosidade cinemática ( \text{m}^2/\text{s}).
  • Para escoamento laminar ( Re < 2000) tem‑se:
f = \frac{64}{Re}.
  • Para escoamento turbulento (comuns em tubulações práticas) o \,f depende de Re e da rugosidade relativa \,\varepsilon/D. A equação implícita de Colebrook é:
\frac{1}{\sqrt{f}} = -2\log_{10}\!\left(\frac{\varepsilon}{3{,}7D} + \frac{2{,}51}{Re\sqrt{f}}\right).
  • Uma fórmula explícita prática (Swamee‑Jain) conveniente para cálculo direto é:
f = \frac{0{,}25}{\left[\log_{10}\!\left(\dfrac{\varepsilon}{3{,}7D} + \dfrac{5{,}74}{Re^{0{,}9}}\right)\right]^2}.
Expliquei cada termo: \varepsilon é a rugosidade absoluta do material do tubo ( \text{m}), D o diâmetro e Re o número de Reynolds. O uso da Swamee‑Jain evita solução iterativa da Colebrook, mas pode-se usar Colebrook junto com iterações.

Passo 7

Perdas localizadas (válvulas, curvas, entradas/saídas).
  • Perdas localizadas têm a forma:
h_m = K \,\frac{V^2}{2g},
onde K é coeficiente adimensional do acessório (tabelas ou fabricantes).
  • No cálculo total de perdas para um trecho, soma‑se todas as perdas localizadas mais a perda por atrito distribuída.

Passo 8

Procedimento prático para um tubo isolado entre dois reservatórios (fluxo determinado por diferença de cabeça). Passos gerais:
  1. Conhecer \Delta H = H_{upstream} - H_{downstream} (m), comprimento L, diâmetro D, rugosidade \varepsilon, e quaisquer perdas locais (coeficientes K).
  2. Escrever equação de energia: \Delta H = h_f + h_m, com h_f em função de Q como em Passo 5 e h_m em função de V (ou Q).
  3. Como h_f depende de f que depende de Re que depende de Q, a equação é geralmente implicita em Q. Resolver por iteração:
    • Chute inicial para f (ex.: f=0{,}02) ou chute inicial para Q.
    • Calcular Q a partir da equação isolada (ver fórmula abaixo).
    • Calcular V = Q/A e Re; atualizar f (Swamee‑Jain ou Colebrook).
    • Repetir até convergência do Q (ou f).
  4. Fórmula útil para atualizar Q dado f:
Q = A \sqrt{\frac{2g\,\Delta H}{\,f\,(L/D)\,}} = A\sqrt{\frac{2g\,\Delta H D}{f L}}.

Passo 9

Exemplo numérico passo a passo (tronco simples) — resolvendo por iteração. Dados escolhidos (números simples para cálculo manual razoável):
  • Diferença de cabeça: \Delta H = 10\\ \text{m}.
  • Comprimento: L = 200\\ \text{m}.
  • Diâmetro: D = 0{,}10\\ \text{m}.
  • Rugosidade: \varepsilon = 0{,}00015\\ \text{m}.
  • Viscosidade cinemática da água: \nu = 1{,}0\times 10^{-6}\\ \text{m}^2/\text{s}.
  • Área: \displaystyle A = \frac{\pi D^2}{4} = \frac{\pi (0{,}10)^2}{4} \approx 0{,}00785398\\ \text{m}^2.
Iteração 1 (chute inicial de f):
  • Chute \,f_0 = 0{,}02.
  • Calcular \,Q_1 usando
Q_1 = A\sqrt{\frac{2g\,\Delta H D}{f_0 L}}.
Substituindo numericamente:
2g\Delta H D = 2\cdot 9{,}81\cdot 10\cdot 0{,}1 = 19{,}62,
f_0 L = 0{,}02\cdot 200 = 4.
\sqrt{\frac{19{,}62}{4}} = \sqrt{4{,}905} \approx 2{,}214.
Q_1 = 0{,}00785398\cdot 2{,}214 \approx 0{,}01738\\ \text{m}^3/\text{s}.
  • Velocidade: \,V_1 = Q_1/A \approx 2{,}214\\ \text{m/s}.
  • Reynolds: \displaystyle Re_1 = \frac{V_1 D}{\nu} \approx \frac{2{,}214\cdot 0{,}1}{1\times 10^{-6}} \approx 2{,}21\times 10^{5} (turbulento).
Atualizar f com Swamee‑Jain:
f = \frac{0{,}25}{\left[\log_{10}\!\left(\dfrac{\varepsilon}{3{,}7D} + \dfrac{5{,}74}{Re^{0{,}9}}\right)\right]^2}.
Calcule cada termo (aproximações numéricas suficientes):
\frac{\varepsilon}{3{,}7D} = \frac{0{,}00015}{3{,}7\cdot 0{,}1} \approx 4{,}05\times 10^{-4}.
Re_1^{0{,}9}\approx (2{,}21\times 10^5)^{0{,}9}\approx 6{,}45\times 10^{4}\quad\Rightarrow\quad \frac{5{,}74}{Re_1^{0{,}9}}\approx 8{,}9\times 10^{-5}.
Soma interna \approx 4{,}944\times 10^{-4}. Logo
f_1 \approx \frac{0{,}25}{(\log_{10}(4{,}944\times 10^{-4}))^2} \approx \frac{0{,}25}{(-3{,}306)^2} \approx 0{,}0229.
Iteração 2:
  • Com f_1 = 0{,}0229 recalcule Q_2:
f_1 L = 0{,}0229\cdot 200 \approx 4{,}58,
\sqrt{\frac{19{,}62}{4{,}58}} \approx \sqrt{4{,}29} \approx 2{,}07,
Q_2 = A\cdot 2{,}07 \approx 0{,}01629\\ \text{m}^3/\text{s}.
  • Velocidade V_2 \approx 2{,}07\\ \text{m/s}, Re_2 \approx 2{,}07\times 10^5.
  • Reaplicando Swamee‑Jain obtém‑se f_2 \approx 0{,}02295, muito próximo de f_1. Convergência alcançada.
Resultado final do exemplo:
Q \approx 0{,}0163\\ \text{m}^3/\text{s}\quad(\approx 16{,}3\\ \text{L/s}),\qquad f\approx 0{,}023.
Comentário: este é o fluxo em regime permanente que equilibra a diferença de cabeça com as perdas por atrito e locais. O procedimento prático é sempre esse: escrever balanço de energia, expressar perdas em função de Q, iterar até convergência do fator de atrito (ou da vazão).

Passo 10

Como tratar redes (várias tubulações em série/paralelo, anéis e nós). Resumo de métodos práticos:
Método direto para ligações em série:
  • Em série, a vazão é a mesma em todos os trechos; somam‑se as perdas: \Delta H = \sum h_{f,i} + \sum h_{m,i}. Procedimento idêntico ao do tronco, com h_f total em função de Q.
Método direto para ligações em paralelo:
  • Em paralelo, a cabeça entre entradas e saídas é a mesma, cada ramo tem sua perda em função da sua própria vazão Q_i. Impõe‑se:
\Delta H = h_{sistema}(Q_i) \quad\text{(igual para todos os ramos)} \quad\text{e}\quad \sum Q_i = Q_{total}.
Resolve‑se tipicamente por iteração até satisfazer conservação de massa no nó.
Método de Hardy‑Cross (método iterativo clássico para redes com malhas):
  1. Chutar vazões iniciais que respeitem continuidade nos nós (sinais arbitrários e orientações de malha).
  2. Para cada laço, calcular o somatório das quedas de cabeça como função das vazões atuais: \sum \text{sign}\cdot h(Q).
  3. Calcular correção de vazão \Delta Q para o laço usando fórmula baseada em derivada das perdas em relação a Q. Ajustar vazões nos ramos do laço.
  4. Repetir para todos os laços até convergir.
    Hardy‑Cross é muito usado em provas e em cálculos manuais porque transforma o problema global em várias correções locais.
Método nodal / matriz (métodos computacionais):
  • Transforma o sistema em equações não lineares nas vazões ou nas pressões nodais. Lineariza (por exemplo, pelas derivadas das perdas) e resolve iterativamente (Newton‑Raphson ou métodos de relaxamento). Método usado por softwares comerciais (EPANET, etc.).

Passo 11

Como integrar bombas e curvas sistema/pump: ponto de operação.
  • A bomba tem uma curva característica h_b = H_{bomba}(Q) fornecida pelo fabricante. O sistema tem sua curva de resistência H_{sistema}(Q) = H_{estático} + h_{perdas}(Q)á.
  • O ponto de operação em regime permanente é a interseção:
H_{bomba}(Q) = H_{sistema}(Q).
  • Procedimento prático: para cada Q calcula‑se H_{sistema}(Q) (somando perdas), e busca‑se a Q que iguale H_{bomba}. Novamente, isto geralmente é feito por iteração.

Passo 12

Dicas práticas para prova e aplicação rápida:
  • Sempre escrever a equação de energia com sinal claro e identificar \Delta H.
  • Calcule área \,A=\pi D^2/4 antes de tudo e mantenha unidades consistentes.
  • Para estimativa inicial, usar f\approx 0{,}02 a 0{,}03 para tubos comuns; depois refine com Swamee‑Jain/Colebrook.
  • Use Reynolds para decidir laminar ( f=64/Re) vs turbulent.
  • Em redes, divida em malhas e aplique Hardy‑Cross se for manual. Para provas, explique o método e mostre 1–2 iterações.
  • Ao apresentar resultados, mostre convergência (ex.: mudança de Q < 1% entre iterações).

Resposta

Em regime permanente aplica‑se conservação de massa e energia estacionárias: escrever \Delta H = h_f(Q)+h_m(Q)-h_{bomba}(Q), usar \,h_f = f(L/D)\,V^2/(2g) com V=Q/A, determinar f por Re (laminar: f=64/Re; turbulento: Colebrook ou Swamee‑Jain) e resolver a equação implícita em Q por iteração; em redes usa‑se Hardy‑Cross ou métodos nodais, e para bombas encontra‑se o ponto de operação pela interseção da curva da bomba com a curva de resistência do sistema.

Exercícios relacionados