Perda de carga em tubulação: como calcular com Darcy-Weisbach e Colebrook, passo a passo

A resposta curta: as três equações

A perda de carga distribuída em uma tubulação cheia se calcula pela equação de Darcy-Weisbach, com o fator de atrito dado pela equação de Colebrook-White, e se converte em pressão pela hidrostática:

h_f = f · (L/D) · (v²/2g)
Δp  = ρ · g · h
1/√f = −2 · log₁₀( ε/(3,7·D) + 2,51/(Re·√f) )

Em ASCII puro, do jeito que costuma ser digitado:

h_f = f * (L/D) * v^2/(2*g)
dp = rho * g * h
1/sqrt(f) = -2*log10( (e/D)/3.7 + 2.51/(Re*sqrt(f)) )

Onde h_f é a perda de carga em metros de coluna do fluido, f o fator de atrito de Darcy (adimensional), L o comprimento reto, D o diâmetro interno, v a velocidade média, g = 9,80665 m/s², ρ a massa específica, ε a rugosidade absoluta e Re = v·D/ν o número de Reynolds. As perdas localizadas entram à parte, como h_loc = ΣK·v²/(2g).

O resto deste artigo é o passo a passo, com um caso resolvido em que todos os números podem ser refeitos.

Passo 1 — Velocidade com o diâmetro interno real

A velocidade é v = Q/A, com A = π·D²/4. O erro clássico é usar o diâmetro nominal: um tubo de 4” Sch 40 (ASME B36.10M) tem diâmetro externo de 114,3 mm e interno de 102,26 mm. Como a perda varia com aproximadamente 1/D⁵, errar o diâmetro em 10% erra a perda em cerca de 60%.

Passo 2 — Reynolds com a viscosidade na temperatura de operação

Re = v·D/ν, com ν = μ/ρ. A viscosidade da água cai à metade entre 10 °C e 40 °C (de 1,306 para 0,653 mPa·s), então “água a 1 cSt” só vale perto de 20 °C. Valores a 1 atm, do NIST Chemistry WebBook, que implementa a formulação IAPWS-95 (massa específica) e a IAPWS 2008 (viscosidade):

Re < 2000 é laminar; Re ≥ 4000 é turbulento; entre os dois fica a zona de transição.

Passo 3 — Rugosidade relativa ε/D

A rugosidade absoluta ε depende do material e do estado do tubo. Aço carbono novo sem costura tem ε típico de 0,045 mm; o mesmo aço com mais de 5 anos de serviço pode chegar a 1 mm ou mais. A tabela de rugosidade da LOGOS traz mínimo, típico e máximo para 31 materiais, novos e usados — é o mesmo conjunto de dados das calculadoras.

Passo 4 — Fator de atrito por regime

• Laminar (Re < 2000): f = 64/Re (Hagen-Poiseuille). Não depende da rugosidade.
• Transição (2000 ≤ Re < 4000): f é fisicamente mal definido. A LOGOS usa a interpolação cúbica de Dunlop (1991), a mesma descrita no manual do EPANET 2.2, que parte de 64/Re em Re = 2000 e chega ao valor de Colebrook em Re = 4000. Trate qualquer resultado nessa faixa como incerto.
• Turbulento (Re ≥ 4000): Colebrook-White (1939), a base do diagrama de Moody (1944). Como f aparece dos dois lados, a equação é implícita. A LOGOS a resolve com um estimador de Serghides seguido de dois passos de Newton, o que chega à precisão de máquina — no caso abaixo, o resíduo da equação é de 10⁻¹⁵.

Passo 5 — Perdas localizadas: K ou comprimento equivalente

Cada válvula ou conexão dissipa h = K·v²/(2g). Somam-se os K de todas as peças do trecho (ΣK) e aplica-se a altura de velocidade do trecho. Os valores de K estão na tabela de coeficientes K da LOGOS.

O método do comprimento equivalente diz o mesmo de outra forma: L_eq = K·D/f é o comprimento de tubo reto que perde o mesmo que a peça. A diferença prática é que L_eq tabelado foi calculado com um f fixo; o método do K usa o f do seu caso. Para fluidos viscosos (acima de cerca de 20 cSt), o K de tabela subestima a perda, e a LOGOS aplica o método 3-K de Darby, que aumenta o coeficiente quando o Reynolds cai.

Passo 6 — Desnível não é perda

O desnível Δz entra no balanço de energia (Bernoulli), mas não é perda de carga: é energia potencial que o fluido ganha ao subir e devolve ao descer. A distinção importa porque a perda varia com aproximadamente Q², enquanto o desnível não depende da vazão. Somar os dois num número só esconde a forma da curva do sistema.

Exemplo resolvido: 80 m³/h de água a 40 °C em 4” Sch 40

Entradas: água a 40 °C (ρ = 992,22 kg/m³; μ = 0,6527 mPa·s); Q = 80 m³/h; aço carbono novo, 4” Sch 40, D = 102,26 mm; ε = 0,045 mm; L = 250 m; desnível +12 m (descarga em tanque elevado). Acessórios:

Cálculo:

1. A = π·0,10226²/4 = 0,008213 m²; Q = 0,02222 m³/s → v = 2,706 m/s.
2. ν = 0,6527·10⁻³ / 992,22 = 0,6578·10⁻⁶ m²/s → Re = 2,706·0,10226/0,6578·10⁻⁶ ≈ 420.600 (turbulento).
3. ε/D = 0,045/102,26 = 4,40·10⁻⁴.
4. Colebrook → f = 0,01745.
5. v²/(2g) = 0,3733 m → h_f = 0,01745·(250/0,10226)·0,3733 = 15,93 m.
6. h_loc = 10,1·0,3733 = 3,77 m (equivale a 59,2 m de tubo reto de 4”).

A conversão usa Δp = ρ·g·h com ρ da água a 40 °C (1 mca = 9,730 kPa). A perda de carga propriamente dita é 19,70 m (191,6 kPa); o restante é o desnível. Dentro da perda, o atrito responde por 81% e os acessórios por 19% — e só os seis cotovelos pesam mais que a válvula de retenção.

Uma leitura de projeto: 2,7 m/s é velocidade alta para uma linha de recalque desse comprimento. É o que a próxima seção mostra.

E se o diâmetro fosse outro?

Mesma vazão (80 m³/h), mesma água a 40 °C, aço carbono Sch 40, perda por atrito a cada 100 m:

A regra “atrito ~ 1/D⁵” confere: de 3” para 4”, (102,26/77,92)⁵ = 3,89 e a perda cai 4,03 vezes; de 4” para 6”, (154,08/102,26)⁵ = 7,77 e a perda cai 7,95 vezes. O expoente efetivo fica em 5,1, ligeiramente acima de 5 porque o tubo menor também tem ε/D maior. Na prática: um degrau de diâmetro muda a perda por um fator de 4 a 8, e nenhuma margem de 10% no cálculo compensa escolher o diâmetro errado.

Colebrook × Swamee-Jain × Hazen-Williams no mesmo caso

Hazen-Williams (forma SI: h_f = 10,67·L·Q^1,852 / (C^1,852·D^4,87)) não tem viscosidade nem temperatura. O C que reproduziria Colebrook neste caso é 141,9 a 40 °C, mas seria 137,9 a 10 °C e 144,1 a 80 °C: o “coeficiente do material” está, na verdade, absorvendo o Reynolds. Por isso Hazen-Williams é razoável para água fria em redes de abastecimento, quando o C foi calibrado para aquela rede, e inadequada para água quente, condensado, óleo, glicol, soluções ou qualquer fluido em regime laminar ou de transição. Swamee-Jain é uma boa aproximação explícita de Colebrook (erro da ordem de 1% na faixa para a qual foi ajustada); com Colebrook resolvida exatamente, ela deixa de ser necessária.

Limites e cuidados

• Monofásico, incompressível, tubo cheio. Darcy-Weisbach com Colebrook pressupõe líquido sem vaporização e seção plena. Escoamento bifásico, canal parcialmente cheio ou flashing exigem outros métodos.
• Gás só com queda pequena. Para gás, o tratamento incompressível só é aceitável quando a queda de pressão é pequena frente à pressão absoluta (a literatura, como o Crane TP-410, usa ordens de 10% da pressão de entrada); acima disso, use formulação compressível.
• Envelhecimento manda mais que a temperatura. No caso resolvido, ε = 0,2 mm (aço com 1–2 anos) leva h_f de 15,93 m para 21,65 m (+36%); ε = 1,0 mm (mais de 5 anos) leva a 34,46 m (+116%). Já a temperatura, de 10 °C a 80 °C, muda h_f só de 16,80 m para 15,47 m. Projete com a rugosidade do fim de vida, não a do catálogo.
• A transição 2000–4000 é incerta. Na faixa de Dunlop, o f deste tubo vai de 0,0286 (Re 2500) a 0,0411 (Re 3999), e há um degrau de cerca de 1,8% na junção com Colebrook em Re = 4000. Nenhum método é “exato” ali.
• K de tabela é valor típico. Fabricantes publicam K ou Cv próprios para suas válvulas; quando houver, use o dado do fabricante.

Do Δp ao dimensionamento da bomba

A perda de carga é uma das parcelas da altura manométrica total: AMT = desnível + perdas na sucção + perdas no recalque (+ diferença de pressão entre os reservatórios, quando houver). Como a perda cresce com Q², o conjunto forma a curva do sistema, e o ponto de operação real é a interseção dessa curva com a curva da bomba.

A calculadora de dimensionamento de bombas da LOGOS faz essa sequência na área logada: você informa ρ, μ e ε do fluido, os trechos de sucção e recalque com diâmetro, comprimento e acessórios, e ela calcula a perda de cada trecho pelo mesmo método deste artigo, encontra o ponto de operação, verifica o NPSH disponível e gera o memorial de cálculo em Word (.docx). A física roda em Edge Functions versionadas no servidor; durante o beta aberto, o acesso completo é gratuito. Os números deste artigo foram calculados com o mesmo código de fator de atrito dessas calculadoras.

Para uma primeira verificação antes do dimensionamento completo — a velocidade numa linha dada, ou o menor tubo que mantém uma vazão abaixo de uma velocidade-alvo —, os cálculos rápidos de hidráulica da LOGOS rodam direto na página pública, grátis e sem cadastro.

O resultado de qualquer ferramenta, inclusive a da LOGOS, não substitui a responsabilidade técnica do engenheiro que assina o projeto.

Referências

• Colebrook, C. F. (1939). Turbulent flow in pipes, with particular reference to the transition region between the smooth and rough pipe laws. Journal of the Institution of Civil Engineers, 11(4).
• Moody, L. F. (1944). Friction factors for pipe flow. Transactions of the ASME, 66.
• Crane Co. Flow of Fluids Through Valves, Fittings, and Pipe — Technical Paper No. 410.
• Swamee, P. K.; Jain, A. K. (1976). Explicit equations for pipe-flow problems. Journal of the Hydraulics Division, ASCE, 102(5).
• Rossman, L. A. et al. EPANET 2.2 User Manual, US EPA — fórmula de transição atribuída a Dunlop (1991).
• ASME B36.10M — Welded and Seamless Wrought Steel Pipe.
• IAPWS-95 (Wagner e Pruß, 2002) e IAPWS 2008 para viscosidade (Huber et al., 2009), consultadas via NIST Chemistry WebBook, SRD 69.