Dimensionamento de rede de bombeamento multi-ramal com solução por Newton-Raphson

O problema da rede ramificada

Quando uma única bomba alimenta vários pontos de consumo ao mesmo tempo, surge um problema que o dimensionamento de um trecho isolado não captura: as vazões dos ramais são interdependentes. Abrir ou fechar um ramal muda a pressão no nó de derivação, e essa mudança redistribui a vazão entre todos os outros. Não existe “a vazão do ramal 3” calculável isoladamente — existe um estado de equilíbrio da rede inteira que precisa ser encontrado de uma vez.

A intuição comum — dividir a vazão da bomba igualmente, ou em proporção à área dos tubos — está errada. A vazão se reparte segundo a resistência hidráulica de cada caminho: ramais curtos, largos e em cota baixa escoam mais; ramais longos, estreitos e altos escoam menos. O que é igual entre eles não é a vazão, é a pressão no ponto onde se derivam.

Os dois princípios físicos

A solução repousa sobre dois pilares clássicos, os mesmos da hidráulica de redes de abastecimento:

• Conservação de massa nos nós. Em todo nó interno, a soma das vazões que entram iguala a soma das que saem (Σ Q_entra = Σ Q_sai). É a versão hidráulica da 1ª lei de Kirchhoff. Fluido incompressível não se acumula.
• Conservação de energia nos trechos. Ao longo de cada tubo, a altura piezométrica cai pela perda de carga de Darcy-Weisbach: h_f = f·(L/D)·(v²/2g) + ΣK·(v²/2g). O fator de atrito f segue o método de Colebrook-White (Serghides no regime turbulento, laminar abaixo de Re 2000, transição cúbica entre eles).

O nó complica o problema: como h_f ∝ Q², as equações são não-lineares. Não há fórmula fechada para a repartição das vazões. É por isso que precisamos de um método iterativo.

Como o Newton-Raphson resolve a rede

O solver monta um vetor de incógnitas com a vazão de cada trecho e um vetor de resíduos F(Q) que reúne todas as equações que devem valer no equilíbrio:

1. Resíduos de continuidade — um por nó interno: Σ Q_entra − Σ Q_sai. Quando vale zero, a massa se conserva ali.
2. Resíduos de pressão — um por ramal-folha: H_nó(Q) − (z + P_min). Quando vale zero, o consumidor recebe exatamente a pressão exigida.

A cada iteração, o algoritmo:

1. Propaga a altura a partir da bomba: H_saída = z_nível + H_bomba(Q_tot) − h_sucção(Q_tot), e desce a rede subtraindo a perda de cada trecho até cada consumidor. A curva da bomba H_bomba(Q) vem de um ajuste de Lagrange em 3 pontos do catálogo.
2. Avalia os resíduos F(Q) com as vazões atuais.
3. Monta o Jacobiano J = ∂F/∂Q por diferença finita central, perturbando cada vazão para cima e para baixo.
4. Resolve o passo de Newton J·ΔQ = −F por eliminação de Gauss com pivotamento parcial e atualiza as vazões.
5. Amortece o passo (line-search): se o passo cheio não reduz a norma dos resíduos, reduz-se o fator até que reduza — garante robustez quando o chute inicial está longe.

A iteração para quando ‖F‖ < 1e-7, o que costuma acontecer em 5 a 15 passos. Como Q_total é uma das incógnitas, a bomba escorrega na própria curva durante a convergência: o ponto de operação real não é imposto, ele emerge da solução conjunta de bomba e rede.

Ramais inalcançáveis e balanceamento

Nem todo ramal recebe vazão. Se um consumidor está alto ou distante demais para a altura disponível, a solução tenderia a uma vazão negativa — fisicamente impossível. O solver detecta esses ramais (vazão que cai a quase zero), marca-os como inalcançáveis e refaz o equilíbrio com os ramais restantes, em vez de produzir um resultado sem sentido. Isso sinaliza, na prática, que a bomba é insuficiente ou que aquele ponto precisa de booster próprio.

No outro extremo está o ramal privilegiado, que “rouba” vazão dos demais. A correção é o balanceamento hidráulico: inserir uma perda localizada deliberada — uma placa de orifício ou válvula de balanceamento — no ramal favorecido, dimensionada (via seu K) para dissipar o excesso de pressão e empurrar a vazão de volta aos ramais desfavorecidos. O resultado é uma rede onde cada ponto recebe a vazão de projeto, não a que o acaso da geometria entregaria.

Considerações práticas de projeto

• O gargalo dita a AMT. A altura manométrica da bomba é governada pelo ramal mais desfavorável (alto, longo, estreito), não pela média. Dimensione a bomba para garantir pressão mínima nele.
• Velocidades de referência: mantenha a linha principal entre 1,5 e 3,0 m/s e os ramais entre 1,0 e 2,5 m/s — controla perda, ruído e golpe de aríete.
• A perda da principal é compartilhada: ela depende da soma das vazões dos ramais. Por isso não se pode dimensionar ramal a ramal de forma independente.
• Cota importa tanto quanto comprimento: dois ramais com a mesma tubulação, mas em cotas diferentes, recebem vazões diferentes — a energia geométrica entra direto no balanço.
• Reavalie ao mudar a configuração: ligar/desligar ramais muda Q_total, desloca o ponto de operação da bomba e redistribui tudo. O cálculo precisa ser refeito, não escalado.

Em resumo: dimensionar uma rede multi-ramal é encontrar o estado de equilíbrio em que a massa se conserva em cada nó e a energia se conserva em cada tubo, simultaneamente. O Newton-Raphson entrega esse estado — vazão de cada ramal, pressão de cada nó e ponto de operação real da bomba — e revela onde o balanceamento é necessário para que o projeto faça o que promete.