Cuánto debe fluir el agua antes de volverse una parábola — Longitud de entrada y Hagen–Poiseuille

La región de entrada en tubería, la parábola de Hagen–Poiseuille y el factor de fricción 64/Re

Al abrir del todo un grifo, el chorro se mantiene liso durante un tramo. Dentro de una tubería ocurre algo parecido. En la entrada, el agua casi no siente la pared. Su perfil es casi plano. Pero tras fluir un poco, la velocidad de la sección se acomoda en una parábola suave. El tramo intermedio es el protagonista de hoy.

La distancia que tarda la pared en difundir la noticia#

En la pared el fluido no se desliza (no-slip: velocidad nula en la superficie). Esa condición no alcanza toda la sección de inmediato en la entrada.

Al principio solo se forma una capa límite delgada (la capa cercana a la pared donde la velocidad cambia bruscamente). A medida que el flujo avanza, esa capa crece hacia adentro. En el instante en que las capas límite superior e inferior se encuentran en el centro, el perfil deja de cambiar. Desde ese punto el flujo se llama completamente desarrollado (fully developed).

La distancia desde la entrada hasta el desarrollo completo es la longitud de entrada hidrodinámica $L_h$ . En flujo laminar la fijan el diámetro $D$ y el número de Reynolds (razón entre fuerzas inerciales y viscosas).

\frac{L_h}{D} \approx 0.05\,\mathrm{Re}_D, \qquad \mathrm{Re}_D = \frac{\rho\,u_m D}{\mu}

Aquí $u_m$ es la velocidad media de la sección, $\rho$ la densidad y $\mu$ la viscosidad. Con $\mathrm{Re}_D=2000$ la longitud de entrada llega a 100 diámetros. Si en una línea de laboratorio delgada se elige mal el punto de medición, se mide un flujo que aún se está desarrollando.

Conviene manipular el número de Reynolds en la simulación de abajo.

Reynolds Re_D800

Al subir Re de 100 hasta 2000, la línea naranja "fully developed" se desplaza a la derecha. El perfil plano de entrada necesita un tramo más largo para asentarse en su parábola.

La parábola de la región desarrollada — Hagen–Poiseuille#

Una vez completamente desarrollado, la velocidad axial deja de cambiar. Se cumple $\partial u/\partial x = 0$ y la velocidad radial $v=0$ . La ecuación de cantidad de movimiento en coordenadas cilíndricas se simplifica.

\frac{1}{r}\frac{d}{dr}\!\left(r\frac{du}{dr}\right) = \frac{1}{\mu}\frac{dp}{dx}

El gradiente de presión $dp/dx$ es constante. Se imponen $u(R)=0$ en la pared y $du/dr=0$ en el centro, y se integra dos veces.

u(r) = -\frac{1}{4\mu}\frac{dp}{dx}\left(R^2 - r^2\right) = 2\,u_m\!\left(1 - \frac{r^2}{R^2}\right)

Aquí $R=D/2$ es el radio de la tubería y $r$ la distancia al centro. El resultado es limpio. La velocidad en el eje es exactamente el doble de la media. Esta parábola es el flujo de Hagen–Poiseuille.

Por qué la presión cae a un ritmo constante — el factor de fricción#

Al promediar la parábola sobre la sección, la velocidad media queda ligada directamente al gradiente de presión.

\frac{dp}{dx} = -\frac{32\,\mu\,u_m}{D^2}

Un gradiente constante significa que la presión baja de forma lineal a lo largo de la tubería. Al empaquetar esa pérdida en un número adimensional se obtiene el factor de fricción de Darcy $f$ .

f = \frac{-\,(dp/dx)\,D}{\tfrac{1}{2}\rho u_m^2} = \frac{64}{\mathrm{Re}_D}

El factor de fricción laminar depende solo del número de Reynolds. La rugosidad de la tubería no interviene. Abajo, deslizar el número de Reynolds mueve el punto a lo largo de la curva del factor de fricción.

Re_D1,000

regime: laminarf = 0.0640L_h/D ≈ 50.0

Para $\mathrm{Re}_D<2300$ el punto sigue con exactitud la recta $64/\mathrm{Re}$ . Al cruzar el umbral, la curva se dobla y cambia a otra ley.

El umbral de 2300 — de laminar a turbulento#

Cuando $\mathrm{Re}_D$ supera unos 2300, el flujo antes liso se rompe. Se mezclan remolinos y transiciona a turbulencia.

En flujo turbulento la longitud de entrada se acorta. $L_h/D$ se queda entre 10 y 60, frente a los cien del laminar. A cambio, la fricción sube. Una correlación empírica para tubería lisa es la ley de Blasius.

f = 0.316\,\mathrm{Re}_D^{-1/4} \quad (4000 \lesssim \mathrm{Re}_D \lesssim 10^5)

El exponente se suaviza de $-1$ a $-1/4$ . Aunque suba el número de Reynolds, la fricción ya no cae tan rápido como antes. Hay que empujar el mismo caudal con más presión.

Reproducir en Python la longitud de entrada y la fricción#

Conviene seguir una tubería pequeña con números.

import numpy as np
 
def entrance_length(Re, D):
    # longitud de entrada hidrodinamica laminar [m]
    return 0.05 * Re * D
 
def poiseuille_profile(r, R, u_mean):
    # perfil de velocidad laminar desarrollado (Hagen-Poiseuille)
    return 2.0 * u_mean * (1.0 - (r / R) ** 2)
 
def darcy_friction(Re):
    # factor de friccion laminar (Re < 2300)
    return 64.0 / Re
 
D = 0.02          # diametro de tuberia [m]
rho, mu = 1000.0, 1e-3
u_mean = 0.08     # velocidad media [m/s]
 
Re = rho * u_mean * D / mu
print(f"Re = {Re:.0f}")                                 # Re = 1600
print(f"longitud de entrada = {entrance_length(Re, D):.2f} m")   # 1.60 m
print(f"velocidad en el eje = {poiseuille_profile(0, D/2, u_mean):.3f} m/s")  # 0.160
dpdx = -32 * mu * u_mean / D**2
print(f"gradiente de presion = {dpdx:.2f} Pa/m")        # -6.40 Pa/m
print(f"factor de friccion f = {darcy_friction(Re):.4f}")  # 0.0400

Una tubería de 2 cm y, aun así, la longitud de entrada es de 1.6 m: 80 diámetros. La mayoría de los tubos cortos de mesa manejan flujo que todavía se está desarrollando.

La próxima vez que aparezca un flujo en tubería#

La esencia del flujo en tubería es "la información de la pared filtrándose hacia el centro". Conviene retener tres cosas.

En la región de entrada el perfil cambia sin parar. La longitud de entrada se estira hasta $0.05\,\mathrm{Re}\,D$ .
El flujo laminar desarrollado es una parábola. El eje va al doble de la media; el factor de fricción es $64/\mathrm{Re}$ .
La transición ocurre cerca de $\mathrm{Re}\approx2300$ . La turbulencia se desarrolla rápido pero cuesta más fricción.