Atravesando un río a nado

Un nadador intenta cruzar un río moviéndose con velocidad constante c sobre el agua. El río tiene una anchura a, la velocidad de la corriente en su centro es máxima v0 y nula en las orillas. Supondremos que el perfil de la velocidad es

v(x)= v 0 ( 14 x 2 a 2 )

que se corresponde con perfil de velocidades de un líquido viscoso que circula en régimen laminar en un tubo de radio a/2 en función de la distancia radial x a su eje

El nadador parte de la posición (-a/2,0) y se mueve sobre el agua hacia la derecha con velocidad constante c, haciendo un ángulo θ=0, con la horizontal

Cuando se encuentra en la posición (x,y), la velocidad del nadador (respecto de tierra) será la suma vectorial de la velocidad del nadador sobre el agua, c =c·cosθ i ^ +c·sinθ j ^ y la velocidad de la corriente en dicha posición v =v(x) j ^

c + v =c·cosθ i ^ +( c·sinθ+v(x) ) j ^

Vamos a determinar la trayectoria que deberá seguir el nadador para emplear un tiempo mínimo T en cruzar el río.

T= a/2 a/2 dx c·cosθ

cosθ es la velocidad horizontal (a lo largo del eje X) del nadador. Dado que la velocidad del nadador c + v , es tangente a la trayectoria. Calculamos la pendiente dy/dx

y ˙ = dy dx = c·sinθ+v(x) c·cosθ = sinθ+ϑ(x) cosθ ,ϑ(x)= v 0 c ( 14 x 2 a 2 )

Despejamos cosθ. Tenemos en cuenta que para que el nadador alcance la orilla derecha, -π/2≤θ≤π/2, cosθ tiene que ser positivo

sinθ= y ˙ cosθϑ(x) ( y ˙ 2 +1 ) cos 2 θ2 y ˙ ·ϑ(x)cosθ+ ϑ 2 (x)1=0 cosθ= y ˙ ·ϑ(x)+ y ˙ 2 +1 ϑ 2 (x) y ˙ 2 +1

El tiempo que tarda el nadador en cruzar el río es

T[ y(x) ]= 1 c a/2 a/2 y ˙ 2 +1 y ˙ ·ϑ(x)+ y ˙ 2 +1 ϑ 2 (x) dx

Ecuación de la trayectoria

Se tratará de encontrar la ecuación de la trayectoria y=y(x) que hace que el tiempo que tarda el nadador T en cruzar el río sea mínimo.

Aplicamos la ecuación de Euler-Lagrange,

f y d dx ( f y ˙ )=0

a un funcional que no depende de la variable y

f(x,y, y ˙ )= y ˙ 2 +1 y ˙ ·ϑ(x)+ y ˙ 2 +1 ϑ 2 (x)

Por este motivo, la obtención de la ecuación de la trayectoria es más sencilla. Derivamos el funcional f respecto a y ˙

f y ˙ = C 1 1 1 ϑ 2 (x) { y ˙ y ˙ 2 +1 ϑ 2 (x) ϑ(x) }= C 1

La constante C1 de determina sabiendo que en la posición inicial x=-a/2, y=0, dy/dx=tanθ=0. Dado que ϑ(-a/2)=0, se concluye que C1=0

Despejamos dy/dx e integramos para obtener la ecuación de la trayectoria y=y(x)

y ˙ =ϑ(x) y ˙ 2 +1 ϑ 2 (x) y ˙ =ϑ(x) dy dx = v 0 c ( 14 x 2 a 2 ) y= v 0 c ( x 4 3 x 3 a 2 )+ C 2

Determinamos la constante C2 sabiendo que para x=-a/2, y=0. La ecuación de la trayectoria es

y= v 0 c ( x 4 3 x 3 a 2 + a 3 )

Cuando el nadador llega a la orilla opuesta x=a/2, su desviación es

y= v 0 c ( 2 3 a )

Dibujamos las trayectorias que describe un nadador que parte de la orilla izquierda de un río de anchura a=2, para tres valores del cociente v0/c: 0.5, 1 y 1.5.

a=2; %anchura del río
hold on
for k=[0.5,1,1.5] %cociente v0/c
    fplot(@(x) k*(x-4*x.^3/(3*a^2)+a/3),[-a/2,a/2], 'displayName',num2str(k))
end
hold off
grid on  
legend('-DynamicLegend','location','northwest')
xlabel('x')
ylabel('y')
title('nadador que atraviesa un río')

La trayectoria es simétrica respecto de la recta horizontal que pasa por el punto ( 0, v 0 c a 3 ) , este es el punto de inflexión de la trayectoria.

Velocidad de la corriente constante

La velocidad de la corriente en el río es máxima en el centro y nula en las orillas. La velocidad media de la corriente de agua en el río es

<v>= a/2 a/2 v 0 ( 14 x 2 a 2 )dx= 2 3 v 0

Supongamos un río cuya corriente lleve una velocidad uniforme e igual a la velocidad media 2v0/3. La velocidad c del nadador en el agua es constante y se dirige hacia la orilla opuesta partiendo de x=-a/2.

El nadador sigue la trayectoria rectilínea de la figura, empleando un tiempo t=a/c, en alcanzar la orilla opuesta y desviándose

y= 2 3 v 0 ( a c )= v 0 c ( 2 3 a )

El mismo resultado que hemos obtenido previamente, suponiendo que la corriente tiene un perfil de velocidades parabólico, la velocidad es máxima en el centro y nula en las orillas

Referencias

Vladimir Ivchenko. Using Poiseuille flow: a refined solution of the riverboat problem. Eur. J. Phys. 40 (2019) 015002