Mecánica de Fluidos (25022) - 1’2023
Control 1
Prof. Leonardo Gordillo
Fecha: 30 de marzo de 2023
Inicio: 11:25 - Fin: 12:55
Apellidos, Nombres:
1. Sumidero de reservorio
El campo de velocidades 2D durante la descarga de un reservorio a través de un
sumidero (por ejemplo al retirarse un tapón en un lavaplatos) puede ser aproximado
por la siguiente expresión:
Q Γ
u=− r̂ + θ̂,
2πr 2πr
donde Q y Γ son constantes que representan la tasa y la rotación espontánea de la des-
carga. La velocidad está expresada con los vectores unitarios de las coordenadas polares.
1. Calcule la divergencia y vorticidad de este campo de velocidades. . . . . . . . . . [1 pto.]
2. Determine las trayectorias que siguen las particulas dada una condición inicial (r0 ,
θ0 ) en t = 0. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . [2 ptos.]
3. Obtenga las líneas de corriente del campo de velocidades e interprete los resultados
gráfica y conceptualmente. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . [3 ptos.]
1
Solución
1.1. Divergencia y vorticidad
En coordenadas cilíndricas, el operador gradiente se escribe ∇ = r̂∂r + 1r θ̂∂θ + ẑ∂z .
Considerando que ∂θ r̂ = θ̂ y que ∂θ θ̂ = −r̂, evaluamos la divergencia
Q Γ 1 Q Γ
∇ · u = r̂ · ∂r − r̂ + θ̂ + θ̂ · ∂θ − r̂ + θ̂ ,
2πr 2πr r 2πr 2πr
Q Q
= − = 0.
2πr2 2πr2
De manera similar, la vorticidad es
Q Γ 1 Q Γ
∇ × u = r̂ × ∂r − r̂ + θ̂ + θ̂ × ∂θ − r̂ + θ̂ ,
2πr 2πr r 2πr 2πr
Γ Γ
=− 2
ẑ + ẑ = 0,
2πr 2πr2
donde hemos usado que el sistema es una base ortogonal tal que r̂ × θ̂ = ẑ.
Es posible demostrar que la divergencia y vorticidad son cero usando también coorde-
nadas cartesianas.
1.2. Trayectorias
Para determinar las trayectorias, conviene usar coordenadas polares e integrar la ecua-
ción de movimiento (descripción Lagrangiana). Esto es resolver,
Q Γ
∂t r = − , r∂t θ = .
2πr 2πr
Se debe primero resolver la ecuación para r (t).
Q Q
2r∂t r = − −→ ∂t r2 = − .
π π
Usando las condiciones iniciales, se puede despejar
r
Q
r (t, r0 , θ0 ) = r02 − t.
π
Para la coordenada polar, debemos reemplazar el valor de r (t) antes de integrar la
ecuación diferencial
Γ Γ 2 Q
∂t θ = −→ ∂t θ = ∂t − ln r0 − t .
2π r02 − Q t 2Q π
π
Finalmente, se reemplaza la condición inicial, para obtener
Γ Q
θ (t, r0 , θ0 ) = θ0 − ln 1 − 2 t
2Q πr0
2
1.3. Líneas de corriente
Para calcularlas, basta con resolver la ecuación diferencial
∂r ∂t r Q
= =− r
∂θ ∂t θ Γ
cuya solución es
Q
r = r0 exp − [θ − θ0 ] .
Γ
Este resultado también se puede obtener resolviendo t en la expresión de r (t) o θ (t) y
reemplazando dicho valor en θ (t) o r (t) respectivamente.
Gráficamente, las trayectorias se ven así
Finalmente, se asignará puntaje por alguna de las siguientes interpretaciones:
Las partículas son atraídas hacia el origen, decayendo en espiral.
La tasa de decaimiento está gobernada por la cantidad Q/Γ. A mayor tasa o menor
rotación de descarga, más rápido es el decaimiento.
p
Las partículas llegan al origen r = 0 en tiempo finito r0 π/Q, pero el ángulo con
el que se aproximan diverge.
Las líneas de corriente corresponden a espirales logarítmicas.
Dada una línea de corriente, se puede construir el resto por rotaciones en el plano.
Dada una línea de corriente, se puede construir el resto por escalamientos en el
plano.
