Por qué una pelota más rugosa vuela más lejos — La crisis de arrastre y la separación de capa límite

Una rugosidad de 0.85 deja el Cd de vuelo en 0.116, pero subirla a 1.0 devuelve ese mismo punto de 70 m/s a 0.124.

A la misma velocidad, una pelota lisa debería cortar el aire con menos arrastre. La pelota de golf rompe esa idea. Lleva cientos de hoyuelos (dimples, las hendiduras pequeñas) y aun así vuela casi el doble de lejos. La brecha es la crisis de arrastre (drag crisis): un tramo de Reynolds donde el coeficiente $C_D$ cae en picado.

Al pasar el botón de capa límite de laminar a turbulenta, el punto de separación se corre hacia atrás y la estela sombreada se estrecha.

Toggle the boundary layer. The turbulent one carries momentum closer to the wall, separates later, and leaves a thinner wake — the shaded region behind the cylinder shrinks and the pressure drag drops.

Conviene mirar cuánto se encoge la estela en el lado turbulento. Cuando ese ancho baja, el arrastre de forma cae a menos de la mitad. La capa más desordenada es la que recorta el arrastre.

Arrastre de presión y arrastre de fricción#

Una pelota siente dos fuerzas. El arrastre de fricción (skin friction), que deja la viscosidad al rozar la superficie. El arrastre de forma (form drag), que deja la diferencia de presión entre frente y cola.

En un cuerpo romo (bluff body) el término de presión es casi toda la carga. El morro queda a alta presión donde choca el aire. La base queda a baja presión donde el flujo ya se fue. Esa brecha tira de la pelota hacia atrás.

Un solo coeficiente carga el conjunto.

$F_D = \tfrac{1}{2}\, \rho\, U^2\, A\, C_D$

$\rho$ es la densidad, $U$ la velocidad relativa, $A$ el área frontal, $C_D$ el coeficiente de arrastre. Lo que cambia la pelota de golf es $C_D$ . El valor se queda plano con el Reynolds y después se desploma.

$Re = \frac{\rho\, U\, D}{\mu}$

$D$ es el diámetro, $\mu$ la viscosidad. $Re$ es la razón entre inercia y fuerza viscosa.

El instante en que la capa límite deja la pared#

El arrastre de forma lo fija el momento en que el flujo deja la pared. Esa salida es la separación de capa límite (separation).

En la mitad delantera el aire acelera y la presión baja. En la mitad trasera ocurre lo contrario. El aire frena y la presión vuelve a subir. Un gradiente adverso de presión (adverse pressure gradient, la presión crece en la dirección del flujo).

El aire junto a la pared ya gastó mucho impulso en la viscosidad. No sostiene ese viento de cara hasta el final. Se detiene y lo empujan hacia atrás. Allí se va la capa, y detrás queda una estela (wake) ancha y de baja presión.

Separación más temprana, estela más ancha. Estela más ancha, más arrastre de forma.

El reparto entre capa laminar y capa turbulenta#

El ángulo de separación depende de si la capa es laminar o turbulenta.

Una capa laminar no mezcla sus capas. El aire junto a la pared es lento y pobre en impulso, así que el gradiente adverso gana pronto. Unos 82° desde el punto de estancamiento frontal.

Una capa turbulenta arrastra hacia la pared el aire rápido de afuera. La zona junto a la pared es más firme y aguanta más el mismo viento de cara. La separación espera hasta unos 120° y la estela se estrecha. El cambio de ancho del botón de arriba es este cambio de ángulo. Por qué la zona de pared aguanta es la misma estructura en capas que un flujo de tubería al desordenarse.

La crisis de arrastre — $C_D$ cae en picado#

Una esfera lisa mantiene $C_D$ cerca de 0.4–0.5 desde $Re \sim 10^3$ hasta $3{\times}10^5$ .

Cerca de $Re \approx 3{\times}10^5$ baja de 0.1. El flujo no ha dado la vuelta. La capa ha transicionado (transition) de laminar a turbulenta y la separación se ha ido hacia atrás. La estela estrecha arrastra el coeficiente por el precipicio. Esa caída es la crisis de arrastre.

Una superficie más rugosa dispara la transición a un $Re$ menor. Los hoyuelos son rugosidad puesta a propósito. El $Re$ de vuelo de una pelota de golf (unos $2{\times}10^5$ ) dejaría a una pelota lisa todavía antes de la crisis, en el lado de $C_D$ alto. Los hoyuelos adelantan el precipicio para que la velocidad de vuelo ya use el valor bajo. La curva que se arrastra a mano queda justo después de los números impresos.

Python — el arrastre de una pelota de golf a ambos lados de la crisis#

La misma velocidad para una pelota lisa y una con hoyuelos. Diámetro 0.0427 m, velocidad 70 m/s.

import math
 
def cd_subcritical(Re):
    # Correlacion de Clift-Gauvin (valida para Re < 3e5)
    return (24.0 / Re) * (1 + 0.15 * Re**0.687) + 0.42 / (1 + 42500 * Re**-1.16)
 
def smoothstep(x, a, b):
    t = max(0.0, min(1.0, (x - a) / (b - a)))
    return t * t * (3 - 2 * t)
 
def drag_coefficient(Re, roughness):
    # roughness 0 (lisa) a 1 (hoyuelos): desplaza el Re critico a la izquierda
    sub = cd_subcritical(Re)
    log_re = math.log10(Re)
    log_crit = 5.55 - 1.35 * roughness
    t = smoothstep(log_re, log_crit - 0.18, log_crit + 0.18)
    sup = 0.08 + 0.13 * smoothstep(log_re, log_crit, 7.0)  # rama tras la crisis
    return sub * (1 - t) + sup * t
 
def drag_force(U, D, roughness, rho=1.2, mu=1.8e-5):
    Re = rho * U * D / mu
    Cd = drag_coefficient(Re, roughness)
    A = math.pi * (D / 2) ** 2          # area frontal proyectada
    Fd = 0.5 * rho * U**2 * A * Cd
    return Re, Cd, Fd
 
D = 0.0427   # diametro de la pelota de golf [m]
U = 70.0     # velocidad justo tras el driver [m/s]
 
for label, rough in [("smooth", 0.0), ("dimpled", 0.85)]:
    Re, Cd, Fd = drag_force(U, D, rough)
    print(f"{label}: Re={Re:.2e}  C_D={Cd:.3f}  F_D={Fd:.3f} N")

Salida:

smooth: Re=1.99e+05  C_D=0.487  F_D=2.049 N
dimpled: Re=1.99e+05  C_D=0.116  F_D=0.488 N

roughness = 0.00Re marker = 1e5.0

cyan: sphere C_D · green: Stokes 24/Re. Raise roughness and the drag-crisis cliff slides toward lower Re — a dimpled ball trips the boundary layer early.

Al subir el deslizador de roughness el precipicio se corre a la izquierda. Con el marcador de logRe cerca de 5.3, la curva lisa sigue en la meseta y la curva rugosa ya está en el fondo. Esos dos asientos son el 0.487 y el 0.116 impresos.

Los 2.049 N y los 0.488 N comparten una presión dinámica. $0.5 \times 1.2 \times 70^{2} = 2940$ Pa por un área frontal cercana a $1.432 \times 10^{-3}$ m² da 4.210 N. Por 0.487 son 2.049 N. Por 0.116 son 0.488 N. Cociente 4.20.

El logaritmo crítico es $5.55 - 1.35 \times$ roughness: 5.55 en liso ( $Re \approx 3.55 \times 10^{5}$ ) y 4.40 con roughness 0.85 ( $Re \approx 2.53 \times 10^{4}$ ). El logaritmo de vuelo es 5.299. La ventana lisa se abre en 5.37, así que 5.299 queda fuera y la rama subcrítica devuelve 0.487. Los hoyuelos ya están en la rama supercrítica, de 0.08 hacia 0.21 en el logaritmo 7. Con roughness 0.85 el smoothstep vale 0.275, y $0.08 + 0.13 \times 0.275 = 0.116$ .

Roughness 1 en ese mismo 5.299 da un smoothstep de 0.341, así que $C_D$ pasa a 0.124 y la fuerza sube a 0.52 N. El centro del precipicio, $Re \approx 2.53 \times 10^{4}$ , son 8.9 m/s. Si la pelota baja a esa velocidad, la superficie rugosa regresa a un $C_D$ alto.

Hoyuelos, y otros escenarios#

El mismo ángulo aparece en otro equipo. Un auto recorta la cola para retrasar la separación y encoger la estela. Una pelota de críquet se pule de un solo lado para desfasar los dos puntos de separación. Un cable de disparo (trip wire, un alambre fino que vuelve turbulenta la capa límite) enrollado en un cilindro adelanta la crisis.

El arrastre lo fija menos la fricción que la viscosidad escribe directo que el momento en que el flujo deja la pared.

Donde más rugosidad devuelve el arrastre#

70 m/s ya está pasado la crisis. Subir allí la roughness por encima de 0.85 devuelve $C_D$ de 0.116 a 0.124. El factor 4.20 impreso no es el extremo del deslizador.

Poner hoyuelos en una superficie que además tiene que pasar calor sube también la fricción de pared. Cómo esa fricción se convierte en transferencia de calor es la pregunta de leer el calor a partir del arrastre. ¿La superficie que adelantó la transición para bajar el arrastre de forma recibe después una factura térmica aparte?