DEPARTAMENTO DE FÍSICA MATEMÁTICA Y DE FLUIDOS

Examen de Fı́sica de Fluidos

4º curso de Grado en Fı́sica
(códigos: carrera 6104, asignatura 6104 4052)

curso 2024-2025 Convocatoria de Febrero (primera semana)

TEST (0,5 puntos cada una)

1. En un flujo a temperatura constante
a) No hay variaciones de densidad.
b) Las variaciones de densidad serán despreciables si las velocidades son suficientemente altas.
c) Las variaciones de densidad nunca son despreciables.
d) Las variaciones de densidad serán despreciables si el tiempo caracterı́stico en el que varı́a la velocidad
es suficientemente pequeño.
e) Ninguna de las anteriores.
2. Los flujos a altos números de Reynolds se caracterizan porque en ellos
a) las aceleraciones debidas a variaciones espaciales de la velocidad son importantes.
b) la vorticidad no puede conservarse a lo largo de las trayectorias de las partı́culas fluidas.
c) las fuerzas debidas a variaciones de presión son del orden de las fuerzas viscosas.
d) las fuerzas viscosas son siempre poco importantes.
3. Dado un fluido con densidad ρ y velocidad v (con componentes cartesianas vi ) explique brevemente
(máximo 2 lı́neas) el significado fı́sico de esta expresión:
ρvi vi .

Solución:
Es una componente del tensor de flujo de momento lineal: componente i de momento lineal por unidad
de volumen transportada por el movimiento del fluido en la dirección i.
4. En la superficie libre de un lı́quido inicialmente en reposo se establece un gradiente de tensión superficial
en la dirección positiva del eje x, y como consecuencia
a) si no se ejercen fuerzas externas o interviene un gradiente de presiones el fluido permanece en reposo.
b) se induce una corriente superficial en la dirección del gradiente de tensión superficial.
c) se induce una corriente superficial en la dirección opuesta al gradiente de tensión superficial.
d) se induce una corriente superficial en la dirección perpendicular al gradiente de tensión superficial.
e) ninguna de las anteriores.

CUESTIONES (2 puntos cada una)

1. Indique de manera razonada cuál es el orden de magnitud de la energı́a cinética turbulenta por unidad
de masa en función de la escala (λ) y la disipación de energı́a (ε) a lo largo del rango inercial, según la
teorı́a de Kolmogorov para la turbulencia homogénea isótropa.
Solución:
La energı́a cinética por unidad de masa en la escala λ es kλ = (1/2)vλ2 ∼ vλ2 . El tiempo caracterı́stico de
esta escala es τλ ∼ λ/vλ . La disipación de energı́a en función de la escala está dada por ε = kλ /τλ = vλ3 /λ.
Como la disipación de energı́a cinética turbulenta es constante a lo largo del rango inercial deducimos
vλ = (ελ)1/3 , de donde
kλ = (ελ)2/3
 2. El potencial de velocidades de un flujo viene dado por
Q
ϕ= ln(r) + U∞ r cos θ
2π
∂ϕ 1 ∂ϕ
donde Q y U∞ pueden considerarse constantes. Sabiendo que vr = ∂r y vθ = r ∂θ :

a) (0.5 puntos) Determine el campo de velocidades en coordenadas cilı́ndricas y en cartesianas.

b) (0.5 puntos) Calcule las lı́neas de corriente.
c) (0.5 puntos) Calcule la trayectoria seguida por la partı́cula fluida que pasa por el punto r = 1,
θ = 0.
d ) (0.5 puntos) compruebe si el flujo es incompresible y/o irrotacional.
Solución:

a) El campo de velocidades se calcula a partir del potencial de velocidades por medio de v = ∇ϕ.
Para obtener las velocidades en coordenadas polares basta hacer
∂ϕ Q 1 ∂ϕ
vr = = + U∞ cos θ y vθ = = −U∞ sin θ
∂r 2πr r ∂θ
El campo de velocidades se puede expresar en cualquier sistema de coordenadas,
v = vr u r + vθ u θ = vx i + vy j
por simple trigonometrı́a se puede observar que
ur = cos θi + sin θj, uθ = − sin θi + cos θj
y por tanto
vx = vr cos θ − vθ sin θ, vy = vr sin θ + vθ cos θ,
resultando en:
Q x Q y
vx = + U ∞ y vy = .
2π (x2 + y 2 ) 2π (x2 + y 2 )
b) (0.5 puntos) Para obtener las lı́neas de corriente basta con imponer
dr rdθ dr rdθ
= que se puede escribir como Q
=
vr vθ 2πr + U∞ cos θ −U∞ sin θ

Con un poco de manipulación algebraica se llega a

d Q
(U∞ r(θ) sin θ) + =0
dθ 2π
y por lo tanto
Q
U∞ (r(θ) sin θ − r(θ0 ) sin θ0 ) + (θ − θ0 ) = 0
2π
que expresa las lı́neas de corriente r = r(θ) en coordenadas polares.
c) Como el flujo es estacionario las trayectorias seguidas por las partı́culas fluidas coinciden con las
lı́neas de corriente. Para calcular la trayectoria indicada en el enunciado sustituimos θ0 = 0, r(θ0 ) = 1
en la expresión general de las lı́neas de corriente. El resultado es:
 
Qθ
r(θ) = − + sin θ0 / sin θ
2πU∞
d ) Como el campo de velocidades está dado por el gradiente de un potencial es evidente que el flujo
es irrotacional. Para ver si es incompresible calculamos su divergencia. A partir del campo de
velocidades en cartesianas tenemos v = U∞ i + (Q/2π)r/r2 . Está claro que el sumando constante
U∞ i tiene divergencia nula y que el factor constante Q/2π no influye. Calculamos entonces la
divergencia de r/r2
2π ∂vx 1 x2 2π ∂vy 1 y2
= 2 −2 4, = 2 −2 4
Q ∂x r r Q ∂y r r
de donde (2π/Q)∇ · v = 2/r2 − 2/r2 = 0 y por tanto el flujo es incompresible.

Examen de Fı́sica de Fluidos

curso 2024-2025 Convocatoria de Febrero (primera semana) 2/5
PROBLEMA (4 puntos)
El sistema que se muestra en la Figura está concebido para mover un fluido, con densidad ρ y viscosidad
µ constantes, entre dos depósitos, separados una distancia L, que se encuentran a presiones constantes P0 y
P0 + ∆P , siendo ∆P > 0.
Dicho sistema esta formado por dos cintas, A y B, separadas una distancia h ≪ L, que se mueven a velocidades
constantes VA y VB , arrastrando el fluido entre los depósitos. Suponiendo que el problema es estacionario, se
pide:
a. (1 punto) Proporcionar la condición necesaria para que el flujo esté dominado por las fuerzas viscosas.

b. (1 punto) Determinar la velocidad del flujo entre las cintas, asumiendo v = (vx (x, y), 0, 0), ası́ como la
distribución de presiones.
c. (1 punto) Calcular el caudal Q arrastrado por las cintas ası́ como la fuerza ejercida sobre las mismas.
d. (0.5 puntos) Si asumimos VA = VB = 0, ¿cuál serı́a la velocidad caracterı́stica del problema?

e. (0.5 puntos) Calcule qué relación que debe cumplirse entre VA , VB y la velocidad caracterı́stica del
apartado anterior para que el sistema funcione, es decir, para que Q > 0.

A VA A
P0 P0+∆P
h
B VB B

Figura 1: Configuración del flujo

AYUDA: Considere la velocidad caracterı́stica del flujo arrastrado por las cintas como V = VA + VB .
Solución:

a. Para que la flujo esté dominado por la viscosidad, una vez asumido que el flujo es estacionario, los
esfuerzos convectivos deben ser despreciables frente a los viscosos. Por la geometrı́a del problema (en
particular la condición h ≪ L) lejos de las regiones de entrada y salida el flujo será aproximadamente
monodireccional (v ≃ (u, 0, 0)), en cuyo caso el término de inercia en NS serı́a nulo, ya que el flujo
es incompresible. Como en este caso el flujo solo es aproximadamente monodireccional, para obtener la
condición que debe cumplirse hacemos uso de las ecuaciones de Navier-Stokes para flujo incompresible
estacionario:
ρ (v · ∇) v = −∇p + µ∇2 v.
Analizando la componente horizontal de velocidad tenemos
 2 
∂u ∂p ∂ u
ρu =− +µ
∂x ∂x ∂y 2

Examen de Fı́sica de Fluidos

curso 2024-2025 Convocatoria de Febrero (primera semana) 3/5
 Podemos obtener el ordenes de magnitud del términos convectivo escalando con el término de esfuerzos
viscosos y usando la velocidad caracterı́stica del flujo arrastrado por las cintas, V :
ρV 2 µV ρV h2
≪ 2 que resulta en ≪ 1.
L h µL
Si definimos el número de Reynolds por medio de la anchura h del canal (que es la longitud caracterı́stica
en la que la velocidad tiene variaciones del orden de la velocidad) la condición anterior queda como
h
Re ≪1
L
y como puede verse la condición h ≪ L, impuesta por la geometrı́a del problema, es determinante para
que el flujo esté dominado por viscosidad. De manera similar, una vez despreciado el término convectivo
µV
se puede deducir que ∆P L ∼ h2 y por lo tanto la condición que debe cumplirse, expresada en términos
de la diferencia de presiones, es
ρ∆P h4
≪ 1.
µ2 L2
b. Al despreciar los términos convectivos podemos simplificar la componente según x de NS
∂p ∂2u
− +µ 2 =0
∂x ∂y
Asumiendo que el flujo es aproximadamente monodireccional, por continuidad u no depende de x, lo
cual obliga a que p varı́e de forma lineal con x, de manera que aplicando las condiciones de contorno del
problema encontramos
p = P0 + ∆P x/L
que expresa el campo de presiones de esta configuración de flujo, ya que al ser h ≪ L las variaciones de
presión con y son despreciables. Sustituyendo en la componente x de NS encontramos
∆P ∂2u
− +µ 2 =0
L ∂y
que se debe integrar con las condiciones de contorno u(y = 0) = VB y u(y = h) = VA . El resultado es
y ∆P
u = VB + (VA − VB ) − y (h − y)
h 2µL

c. El flujo volumétrico Q y la distribución de presiones están relacionadas por:

Z h
1 ∆P h3
Q= u dy = (VA + VB ) h −
0 2 12µL
Esta relación define la velocidad promedio del flujo, dada por
Q 1 ∆P h2
V = = (VA + VB ) −
h 2 12µL
Por otro lado, si despreciamos los efectos verticales de la presión, los esfuerzos cortantes en cada cinta
′
pueden escribirse como τxy = µ∂u/∂y resultando en

′ ∂u VA − VB ∆P h
τxyA =µ =µ + para la cinta A, y
∂y y=h h 2L
′ ∂u VA − VB ∆P h
τxyB =µ =µ − para la cinta B.
∂y y=0 h 2L
RL ′
Sabiendo que la fuerza ejercida sobre las cintas está dada por F = 0
τxy dx tenemos:
L ∆P h
FA = µ (VA − VB ) + para la cinta A y
h 2
L ∆P h
FB = µ (VA − VB ) − para la cinta B.
h 2

Examen de Fı́sica de Fluidos

curso 2024-2025 Convocatoria de Febrero (primera semana) 4/5
 d. De todo lo anterior es obvio deducir que cuando la cinta está parada el flujo va de derecha a izquierda,
dominado por la viscosidad y el gradiente de presiones, por lo tanto en ese caso la velocidad caracterı́stica
es
∆P h2
V0 =
12µL
e. La condición para que el caudal sea positivo es

1 ∆P h2
(VA + VB ) >
2 12µL

Examen de Fı́sica de Fluidos

curso 2024-2025 Convocatoria de Febrero (primera semana) 5/5