Resumen I3 - Mecánica de Fluidos

Análisis Puntual:
1. Conservación de la Masa: Usando el Teorema de la Divergencia (Gauss-Ostrogradsky)
sobre la conservación de la masa en análisis global, y descartando la integral debido a un
volumen cualquiera tenemos:

∂ρ ~ ~ ∂ρ ~ ~ ~ ·V
~ =0
+ ∇(ρV ) = + V · ∇ρ + ρ∇
∂t ∂t
~ = ∂ρ/∂t = 0) tenemos entonces:
Si el fluido es incompresible (∇ρ
~ ·V
∇ ~ =0

2. Conservación de la cantidad de movimiento:

a) Aceleración: Recordemos que desde el punto de vista de Euler:

~
∂V
~a = ~ · ∇)
+ (V ~ V~
∂t
Donde el primer término recibe el nombre de aceleración local y el segundo de aceleración
convectiva o advectiva. Notemos que si: V ~ = u(x, y, z, t)î + v(x, y, z, t)ĵ + w(x, y, z, t)k̂
entonces:
~
∂V ∂u ∂v ∂w        
= î + ĵ + k̂ ∧ ~ ·∇
V ~ V ~ = V~ · ∇u
~ î + V~ · ∇v
~ ĵ + V~ · ∇w
~ k̂
∂t ∂t ∂t ∂t
b) Rotación y Vorticidad: Definimos la Rotación (tendencia del elemento de fluido a
girar sobre si mismo) como
1~
ζ~ = ∇ ×V ~
2
y definimos la Vorticidad como el doble de la rotación:
~ =∇
Ω ~ ×V
~

Ası́:
~ = ωx î + ωy ĵ + ωz k̂
Ω
     
∂w ∂v ∂u ∂w ∂v ∂u
= − î + − ĵ + − k̂
∂y ∂z ∂z ∂x ∂x ∂y
Y en cilı́ndricas:
~ = ωr r̂ + ωθ θ̂ + ωz k̂
Ω
     
1 ∂uz ∂uθ ∂ur ∂uz 1 ∂
1 ∂ur
= − r̂ + − θ̂ + r uθ − k̂
r ∂θ ∂z ∂z ∂r r ∂r r ∂θ
Además como sabemos:
~ · (∇
∇ ~ ×V
~)=0→ ∇
~ ·Ω
~ =0

por lo que podemos definir lı́neas de vorticidad (la vorticidad es conservativa).

1 Rodrigo Henrı́quez Auba

 c) Circulación: Definimos la circulación como (y usando el teorema de Stokes):
I Z  
Γ= ~ · d~l =
V ~ ×V
∇ ~ · dS ~
∂S S

Ası́ la Conservación de momentum en fluidos ideales (µ = 0), usando el Teorema de

la Divergencia en la conservación de momentum en análisis global queda determinada la
ecuación de Euler: " #
∂V~
ρ ~ · ∇)
+ (V ~ V ~ + ρf~
~ = −∇P
∂t

De esto tenemos 3 ecuaciones y 4 incógnitas. Ahora consideremos un fluido incompresible,

y pongamos la ecuación de continuidad (conservación de la masa) ∇ ~ · V = 0 y usemos la
siguiente identidad vectorial:

~ · ∇)
(V ~ V~ = 1 ∇(
~ V~ ·V
~)−V
~ × (∇
~ ×V
~)
2 | {z }
Ω

~ /∂t = 0) e irrotacional (Ω
Por lo que si ahora asumimos un régimen permanente (∂ V ~ = 0)
~ ~
y un campo de fuerzas terrestre f = −g k̂ = −g ∇z nos queda (son 5 condiciones, en
negritas) la ecuación de Bernoulli:

V2 V2
 
~ P P
∇ z+ + =0→z+ + = cte
γ 2g γ 2g

3. Ecuación de Euler en coordenadas naturales: En coordenadas naturales (2D) ŝ y n̂

podemos obtener que en una linea de corriente (dirección ŝ) la ecuación de Bernoulli se
cumple, incluso si el fluido tiene vorticidad. Ası́ tenemos en el eje ŝ en régimen permanente:
∂ V2 P V2 P
 
+ +z =0→ + + z = cte
∂s 2g γ 2g γ
Y en eje n̂:
V2 V2
 
∂ P
= κ= z+
Rg g ∂n γ
Donde κ es la curvatura y R = κ1 es el radio de curvatura. En las rectas la curvatura es 0 y
por ende el radio de curvatura infinito. Ası́ en las corrientes rectas, la cota piezométrica se
conserva.

En el caso que el régimen sea impermanente, tenemos en el eje ŝ:

∂ V2 P
 
∂V
− =g + +z
∂t ∂s 2g γ
Ası́ podemos simplemente integrar desde una posición (1) a otra posición (2) si la tuberı́a
tiene largo L en ese intervalo:
P2 − P1 V22 − V12
 
dV
− L = g (z2 − z1 ) + +
dt γ 2g
Lo que resta es resolver la EDO, y usar intervalos apropiados para cancelar cosas.

2 Rodrigo Henrı́quez Auba

Teorı́a del Flujo Potencial (2D): Para esta teorı́a asumiremos 4 cosas:

~ ·V
1. Fluido Incompresible: ρ = cte → ∇ ~ =0

2. Fluido ideal: µ = 0
~ =∇
3. Flujo Irrotacional: Ω ~ ×V
~ =0

4. Régimen Permanente (no es necesario, solo para efectos de este curso): ∂(·)/∂t = 0

Definiremos 2 funciones de distintas formas:

~ ×V
1. Función Potencial Φ(x, y) : Dado que ∇ ~ = 0 entonces podemos escribir: V
~ = ∇Φ.
~ Ası́:

∂Φ ∂Φ
u= ∧ v=
∂x ∂t
~ ·V
Y si luego usamos la compresibilidad (∇ ~ = 0) tenemos:

~ · ∇Φ
∇ ~ = ∇2 Φ = 0

La ecuación de Laplace, que es lineal.

~ ·V
2. Función de Corriente Ψ(x, y) : Si ahora usamos que ∇ ~ = 0 entonces tenemos que existe
Ψ(x, y) que cumple:
∂Ψ ∂Ψ
u= ∧ v=−
∂y ∂t
que en efecto cumple que:
2 2
~ = ∂ Ψ − ∂ Ψ =0
~ ·V
∇
∂x∂y ∂y∂x
~ ×V
Ası́ que si ahora usamos que el fluido es irrotacional ∇ ~ = 0 tendremos:

∂v ∂u ∂2Ψ ∂2Ψ
− =0→ + = ∇2 Ψ = 0
∂x ∂y ∂x2 ∂y 2
Nuevamente la ecuación de Laplace.

Dado que las lı́neas de corriente están determinadas usando:

~ × d~s = 0 → dx = dy → −vdx + udy = 0 → dΨ = 0

V
u v
Ası́ por lo tanto Ψ no varı́a en las lı́neas de corriente y por lo tanto Ψ(x, y) = C define una
~ ⊥ ∇Φ
lı́nea de corriente. Ası́ ∇Ψ ~ y en 2D: Ψ ⊥ Φ.
Y como en fluido ideal las superficies solidas las puedo modelar como lı́neas de corriente,
tenemos que lı́neas de corriente son posibles sólidos.

Funciones analı́ticas complejas para solucionar la Ec. de Laplace: Dado que las funciones
complejas analı́ticas F (z) cumplen las condiciones de Cauchy-Riemann; luego como z = x+iy = reiθ
con i2 = −1 se tendrá:
F (z) = Φ(x, y) + iΨ(x, y)

3 Rodrigo Henrı́quez Auba

que es armónica, es decir:
∂2F ∂2F
+ =0
∂x2 ∂y 2
válidas tanto para la parte real e imaginaria.

Ası́ funciones como F (z) = U z, U cte corresponde a un flujo uniforme horizontal y F (z) = iU z
corresponde a un flujo uniforme vertical. Funciones como F (z) = az corresponden al doblete,
F (z) = iA ln z corresponde a un fluido que rota y F (z) = A ln z corresponde a un flujo de lı́neas
que salen del origen (ignorar posibles polos en el infinito o en el origen que hacen que la función no
sea analı́tica).

La superposición de funciones (aprovechando la linealidad de la ecuación de Laplace) me per-

mite generar distintos flujos que existan.

Flujo Viscoso: Si definimos T como la matriz de roces viscosos en todas las combinaciones (3×3)
tendremos las ecuaciones de Navier-Stokes para fluidos incompresibles (3D) dada por:
!
∂V~
ρ ~ · ∇)
+ (V ~ V~ = −∇P ~ + ρf~ + ∇T~
∂t

Si el fluido es Newtoniano:
!
~
∂V
ρ ~ · ∇)
+ (V ~ V~ ~ + ρf~ + µ∇
= −∇P ~ 2V
~
∂t

donde
~ 2V
∇ ~ = (∇2 u, ∇2 v, ∇2 w)
No existe solución analı́tica para esto ni obviamente para fluidos compresibles. ¡Resuelvalos, y
podrı́a ganar 1 millón de dolares!
1. Fluido que se mueve en una dirección (en 2D) V ~ = uî, si dividimos por ρ las NSE,
tenemos (ν = µ/ρ):
 2
∂ u ∂2u

∂u ∂u 1 ∂P
+u =− − g k̂ + ν + 2
∂t ∂x ρ ∂x ∂x2 ∂y
~
Usando la ecuación de continuidad ∇·V = 0 → ∂u/∂x = 0 y si asumimos régimen permanente
(∂u/∂t = 0) tenemos en el eje x:

−1 ∂P ∂2u ∂ 2 u ∂P
+ν 2 =0→ µ 2 − =0
ρ ∂x ∂y ∂x ∂x
En la dirección y tenemos 0 = −∂P/∂y entonces la presión no depende de y y asumimos que
la velocidad no depende de x (por continuidad), solo de y y entonces ∂P/∂x no depende de
x por lo que:
∂P
= cte = λ
∂x
ası́ podemos resolver la EDO integrando 2 veces:
λ 2
u(y) = y + yC1 + C2
2µ

4 Rodrigo Henrı́quez Auba

 si usamos las condiciones de borde u(y = 0) = 0 ∧ u(y = e) = V resolvemos:
 
1 ∂P 2 V e ∂P
u(y) = y + − y
2µ ∂x e 2µ ∂x

La ecuación anterior explica que para mover un fluido de izquierda a derecha es moviendo la
placa de arriba hacia la derecha. Otra forma es tener un gradiente de presión, por ejemplo
si V = 0 podemos también generar un movimiento hacia la derecha con un gradiente de
presiones negativo ∂P/∂x < 0, esto debido que el fluido se mueve de zonas de mayor presión
a menor presión. Combinaciones de V y ∂P/∂x me permite generar distintos tipos de u, pero
siempre serán parabólicos.

2. Flujo Laminar en una tuberı́a circular: Si usamos en las NSE coordenadas cilı́ndricas
~ V
para un fluido laminar Newtoniano que se mueve solo en el eje z (∇· ~ = ∂w/∂z = 0), además
∂P/∂r = ∂P/∂θ = 0 → P ≡ P (z) y asumimos régimen permanente (∂w/∂t = 0) se tiene en
la coordenada del eje z la ecuación:

∂P ∂w ∂2w
0 = −r +µ + µr 2
∂z ∂r ∂r
donde nuevamente ∂P/∂z no depende de r ni θ (mirar las otras componentes de la ecuación).
∂w
Ası́ para τ = µ tenemos:
∂r
dP 1 ∂τ 1 ∂
= τ+ = (τ r)
dz r ∂r r ∂r
Que integrando y despejando nos deja:
dP r C1
τ= +
dz 2 r
Dado que el esfuerzo de corte no es infinito en 0, |τ | < ∞ no queda otra que C1 = 0, luego
reemplazando τ :
dw dP r
µ =
dr dz 2
Que si resolvemos:
1 dP 2
w(r) = r + C2
4µ dz
Si usamos la condición de borde que w(r = R) = 0 nos quedará:

1 dP 2
r − R2


w(r) =
4µ dz

5 Rodrigo Henrı́quez Auba