Hydraulic

Multirramal en circuito cerrado de unificación: colectores de impulsión/retorno y equilibrado

En un circuito cerrado de unificación, el caudal de cada ramal no es libre: emerge del acoplamiento de todos los ramales a través de colectores comunes de impulsión y retorno. El método resuelve simultáneamente la conservación de masa y el balance de presión por Newton-Raphson y equilibra los ramales con placas de orificio.

Cuándo usar

Úsalo cuando varios consumidores (intercambiadores, camisas de reactor, cojinetes, puntos de uso) se alimentan desde un colector de distribución (impulsión) y retornan a un colector de recogida (retorno), formando un anillo cerrado de recirculación. Es la topología típica de agua helada, agua de enfriamiento, aceite térmico, sellado y lubricación. A diferencia de un ramal único, aquí la presión en el punto de derivación es compartida: abrir o estrangular un ramal redistribuye el caudal en todos los demás. El dimensionamiento sirve para garantizar que cada consumidor reciba su caudal de proyecto, calcular el caudal total que la bomba del circuito debe impulsar y definir las placas de orificio que equilibran los ramales cortos frente a los largos.

Qué es un circuito cerrado de unificación

Buena parte de los servicios industriales — agua helada, agua de enfriamiento, aceite térmico, sistemas de sellado y lubricación — no alimenta un único punto, sino varios consumidores en paralelo desde un tronco común. La bomba impulsa hacia un colector de impulsión (distribuidor), del que derivan los ramales que atienden cada intercambiador, camisa o cojinete; los retornos confluyen en un colector de retorno que reúne, o unifica, todos los caudales de vuelta al tanque o a la aspiración de la bomba. El resultado es un anillo cerrado de recirculación.

La característica que vuelve difícil este problema — y distinto de dimensionar un ramal aislado — es el acoplamiento: la presión en el punto de derivación de cada ramal no es un dato fijo, depende del caudal que circula por todos los ramales a la vez. Estrangular un consumidor reduce la pérdida en el colector y empuja caudal hacia los demás. Por eso no existe solución ramal a ramal: la red entera debe resolverse de una sola vez.

Por qué el caudal se distribuye de forma desigual

Considera dos consumidores idénticos, uno cerca de la bomba y otro al final del colector. Ambos ven aproximadamente la misma diferencia de presión entre impulsión y retorno, pero el ramal lejano recorre un tramo de colector más largo — con más pérdida de carga acumulada. Como el caudal se ajusta para que cada camino consuma la altura disponible, el ramal cercano, de menor pérdida, capta más caudal, y el lejano queda subatendido. Es habitual que el ramal corto reciba de 1,5 a 3 veces el caudal del ramal largo antes de cualquier corrección.

Ese desequilibrio es puramente hidráulico e inevitable en sistemas ramificados. La solución es el equilibrado: introducir pérdida adicional en los ramales privilegiados hasta que todos consuman la misma diferencia de carga.

El modelo: leyes de los nodos y de las mallas

La red se describe como un grafo, con un caudal incógnita por tramo. Dos familias de ecuaciones cierran el sistema:

  • Conservación de masa (ley de los nodos): en cada derivación, confluencia, entrada y salida de equipo, Σ Q_entra = Σ Q_sale. Como cada tramo tiene caudal independiente, hasta los nodos en serie generan una ecuación.
  • Balance de presión (ley de las mallas): todo nodo alcanzado por más de un camino — una confluencia en el colector de retorno, el cierre del anillo — debe tener la misma altura H calculada por cualquier camino. Equivale a imponer que la suma de las pérdidas en cada malla cerrada sea nula.

La altura en cada nodo se propaga desde la descarga de la bomba, H = z_nivel − h_aspiración + H_bomba(Q_total), descontando la pérdida h_f de cada tramo aguas abajo. La pérdida sigue Darcy-Weisbach, h_f = [f·(L/D) + ΣK]·v²/(2g), con el factor de fricción según el método de Colebrook-White: laminar (f = 64/Re) por debajo de Re 2000, interpolación cúbica de Dunlop en la transición y la estimador de Serghides + Newton sobre Colebrook-White en el régimen turbulento (Re ≥ 4000).

Cómo resuelve la red el Newton-Raphson

Como h_f ∝ Q², el sistema es no lineal y requiere un método iterativo:

  1. Estimación inicial de caudal para cada tramo (mayor en los colectores, menor en los ramales).
  2. Cálculo de los residuos — cuánto se desvía de cero cada ecuación de masa y de presión.
  3. Jacobiano numérico — las derivadas de los residuos respecto a cada caudal se estiman por diferencias finitas centradas (perturbación de 0,001 m³/h), armando la matriz que describe cómo reacciona la red ante pequeños cambios de caudal.
  4. Paso de Newton — se resuelve el sistema lineal J·ΔQ = −residuo por eliminación de Gauss con pivoteo, obteniendo la corrección de todos los caudales a la vez.
  5. Búsqueda en línea (line search) — se aplica una fracción α del paso (reducida a la mitad hasta que el residuo realmente baje), garantizando estabilidad incluso lejos de la solución.
  6. Convergencia cuando la norma del residuo cae por debajo de 10⁻⁷ (típicamente en menos de 15 iteraciones). Hay un límite de 200 iteraciones como salvaguarda.

La ventaja frente al clásico Hardy-Cross — que corrige una malla por vez — es la convergencia global y casi cuadrática: todos los caudales se actualizan en cada paso, lo que importa en circuitos con muchos ramales y equipos en serie.

Equilibrado automático con placas de orificio

Resuelta la red, aparece el desequilibrio: algunos ramales por encima y otros por debajo del caudal de proyecto. El autoequilibrado actúa sobre las hojas (los consumidores reales, nunca sobre los colectores), insertando una placa de orificio en el ramal con presión de sobra. La placa añade pérdida ΔP ∝ Q², dimensionada para disipar exactamente la diferencia de carga entre ese ramal y el ramal crítico (el de mayor pérdida):

K_placa = 2g·Δh_objetivo / v_ramal²

Impuesto el caudal de proyecto en cada hoja, se calcula el exceso de carga de cada ramal y se convierte en K. El resultado es una distribución en la que todos los ramales consumen la misma diferencia de presión entre colectores — y por tanto la misma proporción de caudal solicitado. Los colectores (ramales-padre) no reciben placa automática: si el proyectista quiere aliviar las placas de las hojas (escalonamiento para evitar cavitación), inserta una placa manual en el colector, que entra de forma natural en la propagación de presiones.

Consideraciones prácticas de proyecto

  • Trata la impulsión y el retorno como dos colectores: la pérdida del colector de retorno penaliza a los ramales más lejanos — ignorarla subestima el desequilibrio.
  • No exageres con la placa: un ΔP excesivo genera ruido, vibración y riesgo de cavitación justo aguas abajo; cuando el exceso es grande, repártelo entre la placa del ramal y una placa manual en el colector.
  • Verifica los ramales inalcanzables: una cota alta o una placa agresiva pueden dejar a un consumidor sin altura; el solver lo señala, y la corrección es diámetro, placa, cota o curva de la bomba.
  • Mantén las velocidades dentro de los rangos: 1,5–2,5 m/s en los colectores y 1,0–3,0 m/s en los ramales, con los ramales largos en el límite inferior para no acumular pérdida.

El encadenamiento — armar el grafo, imponer las leyes de los nodos y de las mallas, resolver por Newton-Raphson y equilibrar con placas — entrega una distribución de caudal que respeta el proyecto de cada consumidor y se resuelve con rigor numérico por la conservación de masa y energía.

Fórmulas y fundamentos

Conservación de masa en el nodo (Ley de los Nodos) Σ Q_entra − Σ Q_sale = 0

En cada nodo interno (derivación en el colector de impulsión, confluencia en el colector de retorno, entrada/salida de equipo) la suma de los caudales que llegan iguala a la suma de los que salen. Q en m³/h. Genera una ecuación por nodo de unión y es la base del balance de la red.

Balance de presión en el anillo (Ley de las Mallas) Σ_anillo h_f(Q) = 0 ⇔ H_A(camino 1) = H_A(camino 2)

Para un nodo alcanzado por más de un camino (confluencia de ramales en el colector de retorno, o malla), la altura H calculada por todos los caminos debe coincidir. Equivale a decir que la suma algebraica de las pérdidas a lo largo de cualquier malla cerrada es nula. h_f en mca; H en mca.

Pérdida de carga del tramo (Darcy-Weisbach) h_f = [ f·(L/D) + ΣK ] · v²/(2g) , v = 4Q/(π·D²)

Pérdida total de un tramo sumando la pérdida de carga distribuida y la localizada (accesorios). f factor de fricción [adim.], L longitud [m], D diámetro interno [m], ΣK suma de los coeficientes de los accesorios [adim.], v velocidad media [m/s], g = 9,81 m/s². Define la relación no lineal h_f ∝ Q² que acopla los ramales.

Factor de fricción (Serghides explícito / Colebrook-White) 1/√f = −2·log₁₀( ε/(3,7·D) + 2,51/(Re·√f) ) , Re = ρ·v·D/μ

Régimen turbulento (Re ≥ 4000): estimador de Serghides + 2 pasos de Newton sobre la ecuación de Colebrook-White (precisión de máquina). Para Re < 2000 se usa f = 64/Re (laminar) y entre 2000 y 4000 la interpolación cúbica de Dunlop. ε rugosidad absoluta [m], μ viscosidad dinámica [Pa·s], ρ densidad [kg/m³]. Es el método de factor de fricción de Colebrook-White (Serghides).

Coeficiente de la placa de equilibrado K_placa = 2g·Δh_objetivo / v_ramal² ⇒ ΔP_placa ∝ Q_ramal²

Pérdida adicional impuesta a un ramal con exceso de presión para igualar la altura disponible a la pérdida del ramal crítico. Δh_objetivo es el exceso de presión a disipar [mca], v_ramal la velocidad en el ramal con el caudal de proyecto [m/s]. El autoequilibrado calcula K_placa para que todos los ramales consuman la misma diferencia de carga entre colectores.

Normas y métodos

  • ABNT NBR 12214 — Proyecto de sistemas de bombeo de agua
  • Colebrook-White — estimador de Serghides (turbulento), transición de Dunlop, laminar 64/Re
  • Hydraulic Institute (HI) — análisis de sistemas de tuberías y curva del sistema
  • Crane TP-410 — Flow of Fluids (coeficientes K de accesorios y placas de orificio)
  • ISO 5167 — Medición de caudal con placas de orificio (base del K de placa)

Valores típicos de referencia

Magnitud Rango típico Observación
Velocidad en el colector (impulsión/retorno) 1,5 a 2,5 m/s El colector es un tubo común: una velocidad algo más baja reduce la pérdida acumulada en los ramales lejanos.
Velocidad en los ramales 1,0 a 3,0 m/s Ramales largos en el límite inferior para no penalizar la pérdida; los cortos toleran mayor velocidad.
Desequilibrio antes de equilibrar el ramal cercano recibe 1,5 a 3× el caudal del ramal lejano El ramal más cercano a la bomba tiene más presión disponible y 'roba' caudal; de ahí la necesidad de placas.
Número de ramales admitidos hasta 20 ramales (hojas) Incluye colectores intermedios (subramales-padre) que agrupan consumidores.
Tolerancia de convergencia ‖residuo‖ < 1×10⁻⁷ Newton-Raphson con jacobiano numérico; típicamente converge en menos de 15 iteraciones.
Margen de la placa de equilibrado disipar del 10 % al 40 % de la carga del colector Una placa demasiado agresiva eleva el ruido y el riesgo de cavitación aguas abajo; prefiera repartir entre ramal y colector.

Ejemplo resuelto

Anillo de agua de enfriamiento con 4 consumidores

Datos de entrada

Caudal total de proyecto
Q_total = 80 m³/h
Nº de ramales (hojas)
n = 4 ramales
Caudal objetivo por ramal
Q_ramal = 20 m³/h
Diámetro del colector (impulsión/retorno)
D_h = 150 mm
Diámetro de los ramales
D_r = 80 mm
Longitud ramal corto / largo
L = 12 / 60 m

Resultados

Velocidad en el colector
v_h ≈ 1,26 m/s
Velocidad en el ramal (20 m³/h)
v_r ≈ 1,11 m/s
Caudal del ramal corto sin equilibrar
≈ 27 m³/h
Caudal del ramal largo sin equilibrar
≈ 13 m³/h
K de la placa en el ramal corto
K ≈ 35 adim.

Con 80 m³/h en el colector de 150 mm, v = 4Q/(πD²) ≈ 1,26 m/s, y cada ramal de 80 mm a 20 m³/h da ≈ 1,11 m/s — dentro de los rangos típicos. Al resolver la red sin placas, sin embargo, el Newton-Raphson revela el desequilibrio: el ramal corto (12 m) tiene menos pérdida y capta ≈ 27 m³/h, mientras que el largo (60 m) queda en ≈ 13 m³/h. La diferencia de pérdida entre ambos caminos es el 'exceso' de carga del ramal corto. El autoequilibrado calcula K ≈ 35 para una placa en el ramal corto (y valores menores en los intermedios), disipando ese exceso con ΔP ∝ Q² hasta que todos los ramales consuman la misma diferencia de carga entre los colectores — restaurando 20 m³/h en cada uno. El caudal total que la bomba del circuito impulsa permanece en 80 m³/h, ahora repartido conforme al proyecto.

Haz tu propio cálculo Crea una cuenta gratis y obtén esto con tus propios datos en segundos.

Errores comunes

  • Dimensionar cada ramal aislado, como si la presión en el punto de derivación fuera fija — en circuito cerrado la presión es compartida y el caudal se redistribuye entre todos los ramales.
  • Olvidar la pérdida del colector de retorno: el retorno también acumula pérdida de carga, y penaliza precisamente a los ramales que descargan más lejos en el colector.
  • No equilibrar: sin placas, los ramales cercanos a la bomba 'roban' caudal y los lejanos quedan por debajo del proyecto, incluso con la bomba correcta.
  • Colocar una única placa en el colector en lugar de placas por ramal, creyendo que equilibra — la placa en el colector estrangula todos los ramales (hojas) a la vez y no corrige el desequilibrio entre ellos.
  • Sobredimensionar la placa de equilibrado y generar un ΔP excesivo, causando ruido, vibración y riesgo de cavitación justo aguas abajo del orificio.
  • Ignorar la cota de los puntos de retorno: si el retorno asciende a la parte alta del tanque de origen, la componente geométrica entra en el balance y altera el caudal de equilibrio.

Preguntas frecuentes

¿Por qué no puedo dimensionar cada ramal por separado?

Porque en circuito cerrado todos los ramales comparten los mismos colectores de impulsión y de retorno. La presión en el punto de derivación no es un dato fijo: depende del caudal que pasa por todos los demás ramales. Abrir, cerrar o estrangular un ramal cambia la pérdida en el colector y redistribuye el caudal entre los demás. Por eso la red debe resolverse de forma acoplada, con todos los caudales como incógnitas simultáneas.

¿Qué resuelve aquí el método Newton-Raphson?

Resuelve el sistema no lineal formado por dos familias de ecuaciones: conservación de masa en cada nodo (el caudal que entra iguala al que sale) y balance de presión (la altura calculada por caminos distintos hasta un mismo nodo debe coincidir, es decir, la suma de las pérdidas en cada malla es cero). Como la pérdida h_f crece con Q², el sistema es no lineal; el Newton-Raphson lo lineariza con un jacobiano e itera hasta que el residuo cae por debajo de 10⁻⁷.

¿Cuál es la diferencia con el método de Hardy-Cross?

Hardy-Cross corrige una malla cada vez, iterando sobre los caudales de circulación — converge, pero lentamente en redes grandes o muy acopladas. El Newton-Raphson global arma el jacobiano de toda la red y actualiza todos los caudales en cada paso, con convergencia cuadrática cerca de la solución. Para circuitos de hasta 20 ramales, colectores de impulsión/retorno y equipos en serie, el enfoque global es más robusto y rápido.

¿Para qué sirven las placas de orificio en el equilibrado?

La placa de orificio añade una pérdida de carga controlada (ΔP ∝ Q²) al ramal que tiene presión de sobra — típicamente los ramales cortos y cercanos a la bomba. Sin ella, esos ramales captarían más caudal que el de proyecto. La placa disipa exactamente el exceso de carga para que la diferencia de presión consumida entre los colectores sea igual en todos los ramales, igualando los caudales. Es el equivalente fijo de una válvula de equilibrado.

¿El colector de retorno también cuenta en el cálculo?

Sí, e ignorarlo es un error frecuente. El retorno es un colector: el caudal se acumula a lo largo de él, y la pérdida de carga del colector de retorno penaliza a los ramales que descargan más lejos. El cálculo trata la impulsión y el retorno como dos colectores separados, cada uno con sus segmentos, derivaciones y accesorios, cerrando el anillo en el nodo de unificación.

¿Qué ocurre si un ramal queda inalcanzable?

Si la bomba no tiene altura suficiente para vencer la pérdida hasta un consumidor (cota alta, ramal muy largo, placa exagerada), ese ramal no recibe el caudal de proyecto. El solver detecta el caso, marca el ramal como inalcanzable y lo señala en el resultado — indicando que es necesario revisar el diámetro, la placa, la curva de la bomba o la cota del punto.

Glosario

Colector (header/distribuidor)
Tubería principal de mayor diámetro de la que derivan (impulsión) o a la que confluyen (retorno) los ramales. En circuito cerrado hay un colector de impulsión y uno de retorno.
Circuito cerrado de unificación
Topología en anillo donde la bomba alimenta el colector de impulsión, los ramales consumen el caudal y el colector de retorno reúne (unifica) los retornos de vuelta a la bomba o al tanque.
Equilibrado hidráulico
Acción de igualar el caudal real de cada ramal al de proyecto, normalmente insertando pérdidas (placas o válvulas) en los ramales con presión excedente.
Jacobiano numérico
Matriz de las derivadas parciales de los residuos respecto a cada caudal, estimada por diferencias finitas centradas. Es lo que el Newton-Raphson invierte en cada iteración para actualizar los caudales.
Ley de los nodos y de las mallas
Análogas hidráulicas de las leyes de Kirchhoff: conservación de masa en los nodos y suma nula de pérdidas en cada malla cerrada. Forman el sistema que la red resuelve.
Placa de orificio (equilibrado)
Disco con orificio calibrado que impone una pérdida de carga conocida (ΔP ∝ Q²) para estrangular un ramal y equilibrar la distribución de caudal.