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

Cuesta no suponer que una pelota lisa corta el aire con más limpieza que una rugosa a la misma velocidad. La pelota de golf contradice de frente esa intuición. Lleva cientos de hoyuelos (dimples, las pequeñas hendiduras de la pelota de golf) y aun así vuela casi el doble de lejos que una esfera lisa. Este artículo persigue esa paradoja hasta su origen: la crisis de arrastre (drag crisis). Veremos por qué $C_D$ (el coeficiente de arrastre) cae en picado a cierto número de Reynolds, la separación de capa límite que se esconde detrás, y por qué los hoyuelos adelantan ese precipicio, hasta llegar a un cálculo en Python al final.

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

Una pelota siente resistencia desde dos fuentes. Una es el arrastre de fricción (skin friction) que genera la viscosidad al rozar la superficie. La otra es el arrastre de forma (form drag) que genera la diferencia de presión entre el frente y la parte trasera.

En un cuerpo romo (bluff body) como una pelota, el arrastre de forma domina por completo. La cara frontal se vuelve de alta presión al chocar el aire. La cara trasera queda en baja presión, donde el flujo se desprende. Esa brecha de presión tira de la pelota hacia atrás.

Todo el arrastre se agrupa en un único coeficiente adimensional.

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

Aquí $\rho$ es la densidad del fluido, $U$ la velocidad relativa, $A$ el área frontal proyectada y $C_D$ el coeficiente de arrastre. El secreto de la pelota de golf está en reducir $C_D$. Y $C_D$ varía de forma drástica con el número de Reynolds.

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

$D$ es el diámetro y $\mu$ la viscosidad. $Re$ es la razón entre fuerzas inerciales y viscosas.

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

Lo grande que crece el arrastre de forma lo fija cuándo el flujo deja la pared. Esa es la separación de capa límite (separation).

Sobre la mitad delantera de la pelota, el aire acelera y la presión baja. Viento a favor. Sobre la mitad trasera ocurre lo contrario: el aire desacelera y la presión vuelve a subir. Es un gradiente de presión adverso (adverse pressure gradient, la presión crece en la dirección del flujo).

El aire cercano a la pared dentro de la capa límite ya ha perdido mucho impulso por la viscosidad. No tiene fuerza para resistir ese viento de cara hasta el final. En algún punto se detiene y lo empujan hacia atrás. Allí la capa límite deja la pared, y detrás queda una ancha estela (wake) de baja presión.

Cuanto antes ocurre la separación, más ancha es la estela. Cuanto más ancha la estela, mayor el arrastre de forma.

El reparto entre capa laminar y capa turbulenta#

Aquí viene el giro. El punto de separación depende de si la capa límite es laminar o turbulenta.

Una capa límite laminar es ordenada. Sus capas no se mezclan. La zona junto a la pared es lenta y pobre en impulso, así que cede pronto al encontrar el gradiente adverso, separándose hacia los 82° desde el punto de estancamiento frontal.

Una capa límite turbulenta se revuelve. Arrastra el aire rápido del exterior hacia la pared, de modo que la zona junto a la pared es más firme. Por eso aguanta más tiempo el mismo viento de cara, empujando la separación hasta unos 120°. La estela se estrecha.

Pruébalo en la simulación de abajo.

Cambia la capa límite a turbulenta y el punto de separación retrocede, mientras la región sombreada de la estela se encoge a ojos vistas. Una estela más estrecha reduce el arrastre de forma a menos de la mitad. Contra la intuición, la capa límite más desordenada reduce el arrastre.

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

Si dibujas el $C_D$ de una esfera lisa frente a $Re$, ves una larga meseta. Desde $Re \sim 10^3$ hasta $3{\times}10^5$, $C_D$ ronda 0,4–0,5.

Entonces, hacia $Re \approx 3{\times}10^5$, se desploma de golpe por debajo de 0,1. La dirección del flujo no ha cambiado. La capa límite ha transitado (transition) de laminar a turbulenta, empujando la separación hacia atrás. La estela estrechada arrastra el coeficiente por un precipicio. Esa caída es la crisis de arrastre.

Pruébalo en el gráfico de abajo.

Sube el deslizador de roughness y el precipicio se desplaza a la izquierda. Una superficie más rugosa dispara la transición a turbulencia a un $Re$ menor. Los hoyuelos son superficie rugosa 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 la banda de $C_D$ alto. Los hoyuelos adelantan la crisis para que la pelota ya use un $C_D$ bajo a la velocidad de vuelo.

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

Calculemos directamente el arrastre de una pelota lisa y una con hoyuelos. Hallamos $Re$ a la velocidad de vuelo y evaluamos el coeficiente de arrastre según la rugosidad.

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

A una velocidad de vuelo de 70 m/s, el $Re$ de la pelota de golf es de unos $2{\times}10^5$. Mismo $Re$, y sin embargo ambas se sitúan en lados opuestos de la curva de arrastre. La pelota lisa sigue antes de la crisis ($C_D\approx0{,}49$), así que su estela es ancha. La pelota con hoyuelos, gracias a su rugosidad, ya ha cruzado al otro lado de la crisis ($C_D\approx0{,}12$), así que su estela es estrecha. El arrastre difiere en más de cuatro veces. Esa separación es justo la razón por la que una pelota de golf real supera en vuelo a una lisa.

Hoyuelos, y otros escenarios#

El mismo principio aparece por todas partes. Los diseñadores de coches suavizan la cola para retrasar la separación y encoger la estela. Los lanzadores de críquet pulen solo un lado de la pelota para hacer asimétricos los dos puntos de separación y curvar su vuelo. Enrollar un cable de disparo (trip wire, un alambre fino que vuelve turbulenta la capa límite a propósito) en un cilindro adelanta la crisis de arrastre a voluntad.

Comparten un mismo hilo. Lo grande que crece el arrastre lo fija menos la fricción que la viscosidad genera directamente, y más cuándo el flujo deja la pared.

Tres líneas para recordar#

• El arrastre de un cuerpo romo es en su mayoría arrastre de forma, y su tamaño lo fija dónde se separa la capa límite.
• Una capa límite turbulenta retrasa la separación y estrecha la estela, así que $C_D$ se desploma cerca de $Re \approx 3{\times}10^5$: la crisis de arrastre.
• Los hoyuelos, la rugosidad y los cables de disparo adelantan la transición y permiten usar un $C_D$ bajo a menor velocidad. La pelota más rugosa vuela más lejos.