Apuntes Mecánica de Fluidos Avanzada
Apuntes Mecánica de Fluidos Avanzada
Apuntes Mecánica de Fluidos Avanzada
Mecánica de Fluidos II
3
4 ÍNDICE GENERAL
3.2.2. Ecuaciones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
3.2.3. Estimación del espesor de la capa limite . . . . . . . . . . . . . . . . . . . 78
3.2.4. Estimación de los esfuerzos a la pared . . . . . . . . . . . . . . . . . . . . . 80
3.2.5. Estimación del arrastre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
3.2.6. Ejemplo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
3.3. Parámetros caracterı́sticos de la capa limite dinámica . . . . . . . . . . . . . . . . 81
3.3.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
3.3.2. Espesor convencional de la capa limite dinámica . . . . . . . . . . . . . . . 81
3.3.3. Espesor de desplazamiento . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
3.3.4. Espesor de perdida cantidad de movimiento . . . . . . . . . . . . . . . . . 83
3.3.5. Espesor de energı́a cinética . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
3.3.6. Relación entre espesor de perdida de cantidad de movimiento y el esfuerzo
de corte . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
3.3.7. Factor de forma H . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
3.4. Despegamiento de la capa limite . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
3.4.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
3.4.2. Estudio del despegamiento . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
3.4.3. Localización del punto de separación s . . . . . . . . . . . . . . . . . . . . 87
3.4.4. Ejemplos de despegamiento . . . . . . . . . . . . . . . . . . . . . . . . . . 88
3.5. Transición hacia la turbulencia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
3.6. Ecuaciones de la capa limite laminar: modelo de Prandtl . . . . . . . . . . . . . . 90
3.6.1. Introducción y hipótesis de Prandtl . . . . . . . . . . . . . . . . . . . . . . 90
3.6.2. Ecuaciones para describir el flujo al interior de la capa limite . . . . . . . . 91
3.6.3. Ecuaciones para describir el flujo al exterior de la capa limite . . . . . . . . 93
3.6.4. Ecuaciones de la capa limite: modelo de Prandtl . . . . . . . . . . . . . . . 94
3.7. Solución exacta de las ecuaciones de Prandtl: solución de Blasius . . . . . . . . . . 94
3.7.1. Hipótesis de afinidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94
3.7.2. Configuración del caso estudiado . . . . . . . . . . . . . . . . . . . . . . . 95
3.7.3. Al exterior de la capa limite . . . . . . . . . . . . . . . . . . . . . . . . . . 95
3.7.4. Al interior de la capa limite . . . . . . . . . . . . . . . . . . . . . . . . . . 96
3.7.5. Solución de Blasius . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
3.7.6. Resumen del método de resolución de Blasisus . . . . . . . . . . . . . . . . 99
3.8. Principales resultados obtenidos a partir de la solución de Blasius . . . . . . . . . 100
3.8.1. Calculo de f y de su derivadas . . . . . . . . . . . . . . . . . . . . . . . . . 100
3.8.2. Espesor de desplazamiento . . . . . . . . . . . . . . . . . . . . . . . . . . . 102
3.8.3. Espesor de perdida de cantidad de movimiento . . . . . . . . . . . . . . . . 103
3.8.4. El factor de forma . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104
3.8.5. Evaluación del esfuerzo parietal y del coeficiente de fricción . . . . . . . . . 104
3.8.6. Evaluación del arrastre y del coeficiente de arrastre de la placa . . . . . . . 105
3.8.7. Algunos resultados experimentales . . . . . . . . . . . . . . . . . . . . . . 105
3.8.8. Resumen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
3.9. Ecuaciones integrales para la capa limite . . . . . . . . . . . . . . . . . . . . . . . 107
3.9.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
3.9.2. Velocidad transversal a la frontera de capa limite . . . . . . . . . . . . . . 107
3.9.3. Expresión de la velocidad transversal a la frontera de capa limite en función
de δ ∗ y δ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
3.9.4. Balance global de masa: velocidad aspiración VE (x) . . . . . . . . . . . . . 109
3.10. Ecuación de von Kármán . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 110
3.10.1. Ecuación intégrale para la cantidad de movimiento: Ecuación de Von Kármán110
6 ÍNDICE GENERAL
6.5.7. Expresión de la variación de entropı́a a través una onda de choque derecha 232
6.5.8. Expresión a partir de las condiciones de estancamiento . . . . . . . . . . . 233
6.5.9. Condición de existencia de una onda de choque normal . . . . . . . . . . . 234
6.5.10. Valores para un choque normal (k 6= 1, 4) . . . . . . . . . . . . . . . . . . . 234
6.6. Resumen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 236
6.6.1. Ecuaciones fundamentales para las ondas derechas . . . . . . . . . . . . . . 236
6.6.2. Ecuaciones para las propiedades del flujo a través una onda de choque de-
recha (caso gas ideal) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 236
6.6.3. Relación de Hugoniot (Rankine-Hugoniot) . . . . . . . . . . . . . . . . . . 236
6.6.4. Variación de entropı́a a través una onda de choque derecha . . . . . . . . . 237
6.7. Ondas de choque oblicuas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 237
6.7.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 237
6.7.2. Ecuaciones fundamentales para los choques oblicuos . . . . . . . . . . . . . 237
6.7.3. Ecuaciones fundamentales de las ondas de choque oblicuas para los gases
ideales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 240
6.7.4. Ángulo de choque en función de la desviación y del numero de Mach aguas
arriba . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 241
6.7.5. Relación entre los números de Mach M1 y M2 . . . . . . . . . . . . . . . . 244
6.7.6. Relación de Prandtl del choque oblicuo . . . . . . . . . . . . . . . . . . . . 246
6.7.7. Utilización de tablas de choque derechos para el calculo de choque oblicuos 247
6.7.8. Ondas de choque débiles con desviación muy pequeña . . . . . . . . . . . . 249
6.7.9. Resumen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 250
7. Ayudantias 251
7.1. Capitulo 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 251
7.1.1. Problema 1: Flujo de Poiseuille . . . . . . . . . . . . . . . . . . . . . . . . 251
7.1.2. Problema 2: Flujo laminar sobre una placa plana . . . . . . . . . . . . . . 252
7.1.3. Problema 3: Transportar un fluido con cintas transportadoras . . . . . . . 253
7.1.4. Problema 4: Estudio de un flujo de Couette-Poiseuille modificado (Certamen
2014) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 254
7.1.5. Problema 5: Modelo de un tornado . . . . . . . . . . . . . . . . . . . . . . 255
7.1.6. Problema 6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 255
7.1.7. Problema 7: Ecuación de Laplace para el vórtice de Rankine (Certamen 2014)256
7.1.8. Problema 8: Estimación del valor de la viscosidad dinámica de la astenosfera
(Certamen 2015) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 256
7.2. Capitulo 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 259
7.2.1. Problema 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 259
7.2.2. Problema 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 259
7.2.3. Problema 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 261
7.2.4. Problema 4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 262
7.2.5. Problema 5: Descarga de un jet de agua turbulento en un lago . . . . . . . 263
7.2.6. Problema 6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 265
7.2.7. Problema 7 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 266
7.3. Capitulo 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 266
7.3.1. Problema 1: Hidróptero . . . . . . . . . . . . . . . . . . . . . . . . . . . . 266
7.3.2. Problema 2: Arrastre de una placa . . . . . . . . . . . . . . . . . . . . . . 266
7.3.3. Problema 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 266
7.3.4. Problema 4: Capa lı́mite sobre una pared . . . . . . . . . . . . . . . . . . . 267
7.3.5. Problema 5: Velocidad de sedimentación de una esfera (certamen 2014) . . 267
10 ÍNDICE GENERAL
Appendices 293
Figura 1.1: Volumen arbitrario V . Superficie cerrada S . Elemento de superficie dS con su normal
exterior ~n
Tenemos, entre los instantes, t y t + dt, una cantidad de fluido que entra y sale del volumen.
La variación de la masa m contenida al interior del volumen V es igual (y opuesto) al flujo de
masa que quita la superficie 1 .
Dm ∂ y x
= 0 =⇒ ρf dV + ρf ~v .~n dS = 0 (1.1)
Dt ∂t
V S
1
Se aplica el Teorema de Transporte de Reynolds
13
14 CAPÍTULO 1. COMPLEMENTO SOBRE LAS ECUACIONES DE NAVIER-STOKES
div(~v ) = 0 (1.12)
Existe una condición para considerar un fluido incompresible3 .
y y x
D
ρV~ dV = ρf~ dV + ~σ~n dS (1.16)
Dt
V V S
donde dV representa el volumen de un elemento material de fluido y dS es un elemento de la
superficie cerrada S que limita el volumen V . El tensor σ toma en cuenta las fuerzas de superficie
de presión y de viscosidad que se ejercen sobre el elemento dS . La fuerza de volumen por unidad
de masa f~ aplicada al fluido es por ejemplo: la gravedad o la fuerza electrostática cuando tenemos
un fluido cargado.
D
Comentarios La derivada Lagrangiana Dt está calculada en un referencial que sigue el
movimiento del fluido. En este referencial,
3
Una condición para que un fluido pueda ser considerado como incompresible es:
U c (1.13)
donde U representa la velocidad caracterı́stica del flujo y c es la velocidad de las ondas de presión en el fluido
(velocidad o más precisamente celeridad del sonido). La Ec. 1.13 puede ser escrita de manera adimensional usando
el número de Mach, Ma:
U
Ma = 1 (1.14)
c
Claramente, la condición de fluido incompresible no se verifica cuando estudiamos la dinámica de gas a alta velocidad
(aplicaciones aeronáuticas, ondas de choques).
4
d(m~v ) X
= Fi (1.15)
dt i
16 CAPÍTULO 1. COMPLEMENTO SOBRE LAS ECUACIONES DE NAVIER-STOKES
D y ~ y DV~
ρV dV = ρ dV (1.18)
Dt Dt
V V
Esta relación representa el flujo del vector de componentes (σix ,σiy ,σiz ) a través la superficie S .
Usando el teorema de Green-Ostrogradsky se puede transformar la integral doble en integral
triple:
x x
~σ~n dS = σij nj dS
S S
y ∂σ
ij
= dV
∂xj
V
y
= div~σ dV (1.20)
V
La ecuación
D y ~ y x
ρV dV = ρf~ dV + ~σ~n dS (1.21)
Dt
V V S
La Ec. 1.22 es válida para todo volumen V y si V → 0 tenemos una ecuación local:
DV~
= ρf~ + div~σ
ρ (1.23)
Dt
Es posible descomponer el tensor de los esfuerzos ~σ :
donde δij es la función delta de Kronecker5 . Entonces se puede escribir la matriz σij
vis vis vis
σ11 σ12 σ13 −p 0 0 τ11 τ12 τ13
vis vis vis
σ21 σ22 σ23 = 0 −p 0 + τ21 τ22 τ23 (1.26)
vis vis vis
σ31 σ32 σ33 0 0 −p τ31 τ32 τ33
Como σij y δij son tensores simétricos entonces el tensor viscoso τijvis es simétrico también,
es decir:
DV~
ρ = ρf~ + div~σ (1.31)
dt
puede ser reescrita trabajando sobre la divergencia de ~σ :
∂(pδij ) ∂p
div~σ = div~τ vis − = div~τ vis −
∂xj ∂xi
−−→
= div~τ vis − gradp (1.32)
DV~ −−→
= ρf~ − gradp + div~τ vis
ρ (1.33)
Dt
Por otra parte, también sabemos que la derivada sustancial (o material) se puede escribir para
un campo de velocidad V~ como:
DV~ ∂ V~ −−→
= + (V~ .grad)V~ (1.34)
Dt ∂t
5
La función delta de Kronecker es una función de dos variables, que vale 1 si son iguales, y 0 si son diferentes.
1 si i = j
δij = (1.25)
0 si i 6= j
6
Ver clase de Magister Mecánica de los Medios Continuos.
18 CAPÍTULO 1. COMPLEMENTO SOBRE LAS ECUACIONES DE NAVIER-STOKES
Ası́ tenemos:
∂ V~ −−→ −−→
ρ + (V~ .grad)V~ = ρf~ − gradp + div~τ vis (1.35)
∂t
Al primer término del lado izquierdo de la ecuación anterior se le conoce como aceleración
local y que representa la variación de la velocidad para cada uno de los puntos del espacio
(aceleración en un campo homogéneo no-estacionario V~ (~r, t)).
El segundo término del lado izquierdo de la ecuación anterior se conoce como aceleración
convectiva y representa el cambio de velocidad a medida que se cambia de posición para
un instante dado. Esta aceleración estará presente de igual modo incluso si el campo de
velocidad es estacionario: V~ (~r).
El término ρf~ representa el conjunto de todas las fuerzas de volumen aplicadas al fluido.
−−→
−gradp representa las fuerzas de presión correspondiente a los esfuerzos normales. Esta
fuerza existe en ausencia de movimiento (presión hidrostática).
div~τ vis representa las fuerzas de viscosidad generadas por la deformación del fluido. Este
término contiene los esfuerzos tangenciales y normales que pueden intervenir durante el
movimiento de un fluido compresible o visco-elástico.
Para flujos a altos números de Reynolds, la ecuación de Euler puede ser aplicada lejos de las
paredes si y solamente si el flujo no es turbulento. La ecuación de Euler no se aplica en el caso de
estudio de una capa limite cerca de una placa, o atrás un objeto.
Condición sobre la velocidad tangencial: esta condición depende de las propiedades (vis-
cosidad) del fluido.
Si el fluido es ideal (ninguna viscosidad: ν → 0, el flujo de fluido se describe mediante la
ecuación Euler) no hay otra condición cinemática porque no pueden existir gradientes de
velocidad entre el lı́quido y el sólido, entonces la velocidad tangencial, ~v .~t, a la pared
puede tomar cualquier valor.
Si el fluido es real (viscosidad: ν 6= 0, el flujo de fluido es descrito por la ecuación de Navier-
Stokes completa), se observa experimentalmente que las moléculas del fluido cerca de la pared
están, en promedio, inmóviles sobre todo si la superficie es un poco rugosa. La velocidad
tangencial de una partı́cula de fluido es nula. Tenemos una condición de adherencia
sobre la superficie (pared) sólida:
vk = ~v .~t = 0 (1.41)
20 CAPÍTULO 1. COMPLEMENTO SOBRE LAS ECUACIONES DE NAVIER-STOKES
Si la pared se desplaza, las condiciones en la frontera son fı́sicamente las mismas. Sin
embargo, se necesita escribir la nulidad de las velocidades normales y tangenciales en la pared en
el referencial de la pared, es decir:
Consideramos dos fluidos 1 y 2 presentados en la Figura 1.2. Se debe escribir también una
condición sobre las velocidades al nivel de la interfaz ζ(x, t). En particular, las velocidades
normales deben ser iguales ¡¡ si no, vamos a tener un espacio entre los dos fluidos !!
v1 ⊥ = v2 ⊥ (1.44)
Figura 1.2:
1 1
Pf luido 1 = Pf luido 2 + γ + 0 (1.48)
R R
0
donde R y R son los radios de curvatura de la superficie y γ la tension de superficie. La ley de
Laplace muestra que la presión es más alta del costado cóncavo de la superficie.
1.3.3.2.1. Caso de una superficie sólida Se puede mostrar que en el tensor de los gradientes
de velocidad, varios términos son nulos sobre una pared sólida.
Si la pared es perpendicular a la dirección Oz,
~
Cerca de un punto x tenemos V~ (x + dx) = V~ (x) + ∂V
∂x
dx.
Como la velocidad a la pared es nula en x y tambien en x + dx se deduce que:
∂ V~
=0 (1.49)
∂x
Se demuestra de igual manera que:
∂ V~
=0 (1.50)
∂y
como divV~ = 0 tenemos ∂vz
∂z
= 0 al nivel de la pared.
Ası́ tenemos solamente:
∂Vx ∂Vy
6= 0 y 6= 0 (1.51)
∂z ∂z
El tensor de los esfuerzos tiene las componentes siguientes:
∂Vx ∂Vy
σxz = µ y σyz = µ y σzz = −p (1.52)
∂z ∂z
22 CAPÍTULO 1. COMPLEMENTO SOBRE LAS ECUACIONES DE NAVIER-STOKES
1.3.3.2.2. Caso de una interfase entre dos fluidos sin efecto de tensión de superficie
En presencia de dos fluidos no miscible (Ver Figura 1.3), se puede escribir en la interfase:
En el caso particular donde hay solamente una componente de velocidad según x, la igualdad de
los esfuerzos tangenciales da:
∂v1x ∂v2x
µ1 = µ2 (1.54)
∂z ∂z
Figura 1.3:
Comentarios: Se llama superficie libre a la superficie en la cual los esfuerzos tangenciales son
nulos:
∂Vx
=0 (1.55)
∂z
Esta condición de superficie libre esta justificada, sobre todo, entre el aire y el agua por que
la viscosidad dinámica del aire es muy pequeña frente la viscosidad del agua.
La condición de esfuerzo tangencial nulo impone que la pendiente del perfil de velocidad sea
nulo.
1.3.3.2.3. Caso de una interfase entre dos fluidos con efecto de tension de superficie
Si existe una tensión de superficie, el caso se trata de manera igual al caso precedente y se
adjunta un esfuerzo normal proporcional a la curvatura de la interfase (como en el caso de dos
fluidos ideales) dado por la ley de Laplace. Eventualmente un esfuerzo es tangencial si tenemos
un gradiente de tensión de superficie.
1.3.4. Resumen
Se puede resumir los diferentes casos en la tabla abajo:
Interfase Fluido ideal (Euler) Fluido viscoso (NS)
Interfase
con curvatura
y tensión de • (v⊥ )1 = (v⊥ )2 • ~v1 = ~v2
superficie • pf luido 1 = pf luido 2 + γ R1 + 1
• (σij )1 −(σij )2 = −γ R1 + 1 ~ i
ni +(∇γ)
R0 R0
∂ V~
∀ (r, θ, z, t) =0 (1.58)
∂θ
es un flujo axisimétrico o flujo de revolución invariante por rotación alrededor del eje
polar.
24 CAPÍTULO 1. COMPLEMENTO SOBRE LAS ECUACIONES DE NAVIER-STOKES
∀ (x, y, z, t) W =0 (1.59)
1.4.2. ~
El potencial vector de velocidad ψ
1.4.2.1. Definición
~ y, z, t) tal que:
Llamamos potencial vectorial de velocidad el vector ψ(x,
−→~
V~ = rotψ (1.60)
Observamos también por definición que el potencial vector puede ser definido solo o con el (suma
de un) gradiente de una función arbitraria:
→ ~ −−→
−
V~ = rot ψ + gradf (1.62)
Teorema: El caudal que pasa a través la superficie S limitada por la curva cerrada C es igual
a la circulación del potencial vector de velocidad a lo largo de la curva C
Comentario: El valor del caudal se queda igual sea cual sea la superficie S considerada, si es
limitada por la misma curva C
1.4. FUNCIÓN DE CORRIENTE 25
Figura 1.4:
1.4.3.3. Comentarios
Para todo los problemas de flujo bidimensional, si el fluido es viscoso o no viscoso (ideal), se
puede siempre introducir la función de corriente Ψ(x, y). En la parte flujo potencial se va a definir
la función potencial de velocidad Φ, esta función es solamente definida para un fluido ideal (sin
viscosidad).
En el caso de problema donde las coordenadas polares son necesarias (r, θ) en el plano del
flujo, la relación entre las componentes de velocidad y la función de corriente Ψ se escribe:
1 ∂Ψ(r, θ) ∂Ψ(r, θ)
vr = vθ = − (1.71)
r ∂θ ∂r
Llamamos lı́nea de corriente, al lugar geométrico de los puntos tangentes al vector velocidad
de las partı́culas de fluido en un instante t determinado.
Si llamamos dM ~ el elemento de arco de una lı́nea de corriente, la ecuación vectorial diferencial
de esta linea se escribe:
∀t fijo, dM ~ × V~ (M, t) = ~0 (1.72)
1.4.4.2. Teorema
Para un flujo bidimensional plano, la función corriente es constante a lo largo de una lı́nea de
corriente.
~ × V~ = ~0
dOM (1.73)
dx dy
= (1.74)
U (x, y, t) V (x, y, t)
con
∂Ψ(x, y, t) ∂Ψ(x, y, t)
U (x, y, t) = y V (x, y, t) = − (1.75)
∂y ∂x
La ecuación 1.74 permite escribir:
∂Ψ(x, y, t) ∂Ψ(x, y, t)
dΨ = dx + dy = 0 (1.76)
∂x ∂y
Demostración Para establecer esta propiedad, evaluamos el caudal Q por una espesor unidad
en la dirección Oz (normal ~n = (0; 0; 1)). Tenemos:
Z M2 Z M2
Qv = ~ =
V~ .dS ~ × ~n)
V~ .(dl (1.79)
M1 M1
Figura 1.5:
Z M2 " Z M2 #
∂Ψ ∂Ψ
Qv = dy + (− (−dx)) = dΨ = Ψ2 − Ψ1 (1.81)
M1 ∂y ∂x M1
Miramos el perfil de velocidad y las lı́neas de corriente para varios valores del ratio: ab .
a
1.4.6.2. El caso: b
=0
Es un flujo generado por un cizallamiento (flujo de Couette simple). Esto flujo se obtiene cuando
un fluido viscoso es posicionado entre dos placas (distancia d) y que las placas son en movimiento
(velocidad U ). El gradiente de velocidad es Ud . Las lı́neas de corriente corresponden a y =Cte, son
paralelas a los planos (ver Figura 1.6).
a
1.4.6.3. El caso: b
= −1
La función de corriente es la ecuación de una hipérbola equilátera, sus ası́ntotas son: y = ±x. Si
calculamos el rotacional del campo de velocidad tenemos:
!
−→ ∂vy ∂vx
rotV~ = − ~z = 0 (1.85)
∂x ∂y
a
1.4.6.4. El caso: b
=1
Si ab = 1 tenemos Ψ = a(x2 +y 2 ) = Cte y las componentes del vector velocidad son: vx = 2ay
y vy = −2ax
Las lı́neas de corriente son cı́rculos concéntricos (Ver Figura 1.9):
Ψ = a x2 + y 2 = Cte (1.86)
−
→
Si calculamos el rotacional del campo de velocidad tenemos: rotV~ = −4~z. Este tipo de flujo
corresponde a una rotación pura según ~z. La velocidad angular
1− → ~
Ω = rotV = 2a
(1.87)
2
donde:
Ψciz corresponde a un cizallamiento puro (los coeficientes de x2 y y 2 son opuestos)
Ψrot corresponde a una rotación pura (los coeficientes de x2 y y 2 son iguales). El modulo Ω
de la velocidad de rotación angular es: b + a.
La relación entre la función de corriente Ψ y el campo de velocidad es lineal (derivación).
Entonces, el campo de velocidad es la suma de los campos de velocidad correspondiendo a una
rotación pura y a un esfuerzo cortante:
vx = (b + a)y + (b − a)y
vy = −(b + a)x + (b − a)x (1.89)
El análisis del caso general permite encontrar por supuesto, los casos presentados anteriormente:
flujo esfuerzo cortante: b = a
−1 ∂Ψ 1 ∂Ψ
vr = vz = (1.92)
r ∂z r ∂r
Se usa las coordenadas esféricas (r, ϕ, θ). La ecuación de conservación de masa se escribe:
En el caso de un flujo axisimétrico, invariante por rotación alrededor del eje Oz, las componentes
de la velocidad son independientes de ϕ, la Ec. 1.93 se escribe:
1 ∂Ψ 1 ∂Ψ
vr = vθ = − (1.95)
r2 sin θ ∂θ r sin θ ∂r
Varios ejemplos de utilización de la función de corriente Ψ(x, y) serán dadas en la parte: función
potencial.
Figura 1.10:
V~ = Vx (r)~ex (1.99)
∂p
0 = − − ρg sin θ
∂r
1 ∂p
0 = − − ρg cos θ
r ∂θ " #
∂p ∂ 2 Vx 1 ∂Vx
0 = − +µ + (1.100)
∂x ∂r2 r ∂r
Si suponemos que la presión se queda constante en toda sección del conducto, la presión es
entonces función de x solamente:
p = p(x) (1.102)
1.5. ALGUNAS SOLUCIONES ANALÍTICAS PARA FLUJOS LAMINARES 33
La condición de simetrı́a en r = 0:
dVx
Vx (r = 0) = Vmax sea =0 (1.103)
dr r=0
Vx (r = R) = 0 (1.104)
1.5.1.4. Integración
Integramos dos veces con respecto r, la ecuación 1.101:
!
∂ ∂Vx 1 ∂p
r = r (1.105)
∂r ∂r µ ∂x
sea:
∂Vx 1 dp r2
r = + C1
∂r µ dx 2
1 dp r2
Vx (r) = + C1 ln r + C2 (1.106)
µ dx 4
1 dP 2
Vx (r) = r + C1 ln r + C2 (1.107)
4µ dx
1 dP 2
Vx (r) = r + C1 ln r + C2 (1.108)
4µ dx
Las constantes C1 y C2 son las constantes de integración que se determinan usando las condiciones
de fronteras 1.103 y 1.104.
1 dP 2
0= R + C2 (1.109)
4µ dx
1 dP 2
Vx (r) = − R − r2 (1.110)
4µ dx
Llamamos solución de Hagen-Poiseuille, la velocidad Vx (r) obtenida a partir de Ec. 1.101.
34 CAPÍTULO 1. COMPLEMENTO SOBRE LAS ECUACIONES DE NAVIER-STOKES
1.5.1.7. El caudal Q
El caudal Q que atraviesa la sección es definido por:
Z R
Q = Vx (2πrdr)
0
Z R !
1 dP 2
= − R − r2 (2πrdr)
0 4µ dx
! !
π dP 4 π dP
= − R = − D4 (1.113)
8µ dx 128µ dx
donde D = 2R es el diámetro del conducto.
Pentrada −Psalida
Como el gradiente de presión axial es − dP
dx
= L
, podemos escribir:
πR4 Pentrada − Psalida
Q= (1.114)
8µ L
El caudal es proporcional al gradiente de presión y al radio con exponente 4.
1.5.2.2. Hipótesis
Tratamos el caso de un flujo laminar y hace varias hipótesis:
Es decir, un flujo donde el campo de velocidad y de presión son independientes de x y θ. El
cual corresponde a un flujo observado con pequeñas velocidades.
Supongamos, la ausencia de gradiente de presión axial (tenemos vx = 0). Entonces, tenemos
una simetrı́a del problema con respecto a los planos perpendiculares a Ox.
Además, hacemos una última simplificación, vθ es independiente de θ por simetrı́a.
36 CAPÍTULO 1. COMPLEMENTO SOBRE LAS ECUACIONES DE NAVIER-STOKES
Figura 1.11:
1 ∂
(rvr ) = 0 (1.126)
r ∂r
porqué tenemos: vθ = Cte y vz = 0.
vr (r = R1 ) = vr (r = R2 ) = 0 (1.128)
v2
La primera ecuación expresa el equilibrio entre la fuerza de inercia rθ provocado por la curva-
tura de las trayectoria de las partı́culas y el gradiente de presión radial.
La segunda ecuación corresponde al equilibrio entre las diferentes componentes de la fuerza de
viscosidad en coordenadas cilı́ndricas.
La última ecuación, Ec. 1.129, se escribe:
!
∂ 2 vθ 1 ∂vθ vθ ∂ ∂vθ vθ
2
+ − 2 =0⇒ + =0 (1.130)
∂r r ∂r r ∂r ∂r r
∂vθ vθ ∂(rvθ )
+ = C1 ⇒ = rC1 (1.131)
∂r r ∂r
Una segunda integración:
C1 2
rvθ = r + C2 (1.132)
2
Ası́:
1 C2
vθ (r) = rC1 + (1.133)
2 r
Las constantes de integración C1 y C2 son determinadas por las condiciones de frontera según θ:
entonces:
(Ω1 − Ω2 )R12 R22 Ω2 R22 − Ω1 R12
C2 = C1 = 2 (1.135)
R22 − R12 R22 − R12
Finalmente, tenemos:
(Ω1 − Ω2 )R12 R22 1 Ω2 R22 − Ω1 R12
vθ (r) = + r (1.136)
R22 − R12 r R22 − R12
1.5.2.5. Comentarios
Tres casos particulares pueden ser considerados:
Figura 1.12:
1.6.1.1. Definición 1
Llamamos la vorticidad ω
~ el campo vectorial definido por el rotacional del campo de veloci-
dades:
−
→
~ = rot~v
ω (1.137)
1.6.1.2. Definición 2
Un flujo es llamado irrotacional si:
−→
rotV~ = 0 (1.138)
es decir: la vorticidad es nula (~ω = 0), los elementos del fluido no están puestos a rotación.
Un ejemplo de un flujo rotacional y irrotacional es presentado sobre la Figura 1.13.
gradV~ = S + R (1.139)
!
1 ∂Ui ∂Uj
Rij = − (1.142)
2 ∂xj ∂xi
~ llamado vector rotación:
Usando esta analogı́a podemos definir un vector Ω,
con ijk = el pseudo-tensor alternada, 123 = 231 = 312 = 1; 321 = 213 = 132 = −1, ijk = 0 si a
los menos 2 indices son iguales.
Ası́, podemos notar:
~ = 1−
Ω
→ 1~
rot~v = ∇ × ~v (1.144)
2 2
Finalmente:
Se llama vorticidad al doble del tensor rotación o al rotacional del vector velocidad:
−→ ~
~ = rot~v = 2Ω
ω (1.145)
1.6.2. Comentarios
1.6.2.1. Comentario 1
De manera intuitiva, un flujo con vorticidad puede ser definido por dos parámetros: su tamaño
(diámetro D = 2r), su velocidad periférica (tipo rotación sólido).
Un perfil de velocidad vθ tı́pico en un vórtice es presentado en Figura 1.14.
La vorticidad ω tiene solamente una componente paralela al eje del vórtice. Se observa que la
vorticidad es concentrada en el corazón del vórtice.
1.6.2.2. Comentario 2
El cálculo de la vorticidad ω ~ es un parámetro importante porque localiza los vórtices y pro-
porciona sus intensidades.
Un mecanismo importante de amplificación de la vorticidad (aumento de la intensidad del
vórtice) es la elongación del vórtice. Si tiremos en el sentido de la longitud, el vórtice girara mas
40 CAPÍTULO 1. COMPLEMENTO SOBRE LAS ECUACIONES DE NAVIER-STOKES
1.6.2.3. Comentario 3
La vorticidad permite cuantificar como rota en función del tiempo el entorno local de un punto.
Para describir el campo vectorial ω ~ se usa el mismo enfoque empleado para el campo de
velocidad: linea de vorticidad, tubo de vorticidad, etc.
!
1 ∂Vθ Vr 1 ∂Vθ 1 ∂Vz
Sθθ = + Szθ = Sθz = + (1.148)
r ∂θ r 2 ∂z r ∂θ
!
∂Vz 1 ∂Vr ∂Vz
Szz = Szr = Srz = + (1.149)
∂z 2 ∂z ∂r
Llamamos lı́nea de vorticidad (o lı́nea de rotación) toda curva en el espacio donde la tangente
en cada punto es paralela a cada instante t fijo al vector vorticidad ω
~ (Ver Figura 1.15).
~ = ~0
d~x × ω (1.150)
sea:
dx1 dx2 dx3
= = (1.151)
ω1 ω2 ω3
donde d~x = (dx1 , dx2 , dx3 ) es tangente a la curva en coordenadas curvilı́neas.
Los tubos de vorticidad son definidos a partir de las lı́neas de vorticidad de la misma manera
que son definidos los tubos de corriente a partir de las lineas de corriente (Ver Figura 1.16).
42 CAPÍTULO 1. COMPLEMENTO SOBRE LAS ECUACIONES DE NAVIER-STOKES
Consecuencia importante de la propiedad Como, los flujos de vorticidad a través dos sec-
ciones (cualquiera) de un tubo son igual, todo tubo de vorticidad debe:
o extenderse al infinito,
1.6.6. Circulación
1.6.6.1. Definición
Consideramos un contorno cerrado C en un flujo. La circulación a lo largo de C es la integral
a lo largo de C de la velocidad tangente a C .
Si notamos d~l un vector tangente a C en un punto M de C , de norma d~l, la circulación vale:
I
Γ= ~v .d~l [m2 .s−1 ] (1.157)
C
1.6.6.2. Propiedad
−→
~ = rot~v . Esta relacionada con la circulación por el teorema
La vorticidad es definida como: ω
del rotacional: x −
I
→ − → x − →
Γ = ~v .d~l = rot~v .dS = ω
~ .dS [m2 .s−1 ] (1.158)
C
S(C ) S(C )
44 CAPÍTULO 1. COMPLEMENTO SOBRE LAS ECUACIONES DE NAVIER-STOKES
1.6.6.3. Comentarios
Se puede usar la relación 1.158 para cualquier contorno sobre las paredes del tubo de vorticidad
que delimite una superficie S. El flujo de vorticidad a través S es igual a la circulación a lo largo
de dC .
Además, como sabemos que el flujo de vorticidad no depende de S:
del vector vorticidad a través una superficie S que seguimos en su movimiento es contante.
d I ~ d x dΓ
~v .dl = ~ .~n dS =
ω =0 (1.159)
dt C dt S dt
Las condiciones sobre el fluido son las mismas que usamos para establecer el teorema de Bernoulli
Restricción Este teorema no es verdad en el caso de los fluidos disipativos, es decir viscosos.
La vorticidad ω
~ crece por un movimiento de elongación paralelamente a ω
~
1.7.6. Ejemplos
1.7.6.1. Primero ejemplo: las lı́neas de vorticidad son cı́rculos concéntricos
Consideramos por ejemplo un flujo donde las lı́neas de vorticidad son cı́rculos concéntricos. La
velocidad es dada por ~v = ω e~z × re~r
Calculamos la circulación a lo largo de un circulo C0 de radio r0 :
I I
Γ= ~v .d~l = ω r0 (e~z × ~er ).(r0~eθ dθ) = 2πωr02 (1.161)
C0 C0
1
1.7.6.2. Segundo ejemplo: un campo de velocidad en r
K
~v = ~eθ (1.162)
r
donde K es una constante. Calculamos la circulación a lo largo de un circulo C0 de radio r0 :
K K
I I
Γ= ~v .d~l = ~eθ . (r0~eθ dθ) = 2πr0 = 2πK (1.163)
C0 C0 r0 r0
1
Este resultado muestra que si el campo de velocidad es en r
la circulación Γ es constante.
Además,
−→
rot(∆~v ) = ∆~ω (1.169)
Ası́, a partir de Ec. 1.165 tenemos la ecuación de la vorticidad (Ec. 1.170).
1.8. ECUACIÓN DE VORTICIDAD (ECUACIÓN DE HELMHOLTZ) 47
2.1. Introducción
La mayorı́a de los flujos de interés tecnológicos son turbulentos: el jet que sale de un turbo-
rreactor, el penacho de humo de una chimenea, la capa lı́mite que se establece sobre el ala de un
avión, el agua que circula en el conducto de un embalse, el fluido que gira en un reactor quı́mico,
los flujos de convección natural y forzados, . . .
Ası́ mismo, la mayorı́a de los flujos naturales son turbulentos: flujos en canales, rı́os, flujos en
bordes costeros, flujos atmosféricos, . . .
Por lo que el estudio de la turbulencia es necesario para describir y comprender estos
flujos. Sin embargo, varias preguntas aún están abiertas, tales como: el origen de la turbulencia,
el papel que juegan las inestabilidades en el proceso de transición: laminar - turbulento, la presen-
cia de estructuras ordenadas en los flujos, los intercambios energéticos entre las diversas escalas
de fluctuaciones, el rol de las intermitencia . . . Todas estas preguntas no serán tratadas en este
capı́tulo.
En este capı́tulo, vamos a presentar las bases del tratamiento estadı́stico de la turbu-
lencia, dar nociones simples sobre la descomposición entre los campos promedios y fluctuantes,
entregar algunas propiedades de los flujos turbulentos (en particular su energı́a cinética turbulenta,
su difusividad, . . . ), introducir también la noción de longitud de mezcla. Finalmente, presentaremos
algunos modelos como RANS (Reynolds Averaged Navier Stokes) y LES (Large Eddy Simulation),
...
49
50 CAPÍTULO 2. INTRODUCCIÓN A LOS FLUJOS TURBULENTOS
La turbulencia siempre se extiende a través del fluido. Las transferencias de masa, cantidad
de movimiento y de energı́a crecen fuertemente. La turbulencia aumenta la fricción parietal, las
pérdidas de cargas y los flujos de calor. La turbulencia acelera los procesos de mezcla (como por
ejemplo, la dispersión de los contaminantes).
La turbulencia aparece en los flujos cuando el valor del número de Reynolds, Re, es elevado.
Generalmente, su origen se debe a una inestabilidad presente en el flujo laminar. Cuando el número
de Re alcanza un cierto valor crı́tico, las pequeñas perturbaciones del flujo son inestables
y, el crecimiento de estas, produce la transición hacia la turbulencia. Ası́, cuando tenemos una
capa lı́mite laminar que se desarrolla (Ver Figura 2.2), por ejemplo en el fondo de un rı́o (que
puede ser considerado liso o no liso), el flujo es turbulento cuando ReL > 5 × 105 . En un conducto
cilı́ndrico, el régimen aparece para un ReD > 2000
Figura 2.2: Transición del flujo laminar hacia el flujo turbulento sobre un placa plana.
2.1. INTRODUCCIÓN 51
2.1.1.8. Resumen
En resumen, los flujos turbulentos aparecen cuando el número de Reynolds es grande, son
fuertemente rotacionales, están caracterizados por su irregularidad, son muy disipativos, tienen
varias escalas, tienen estructuras organizadas y aceleran el proceso de mezcla.
52 CAPÍTULO 2. INTRODUCCIÓN A LOS FLUJOS TURBULENTOS
Figura 2.4: Flujo turbulento alrededor un esfera de diámetro D. La velocidad del flujo es U0 .
Una análisis dimensional muestra que la tasa de disipación de la energı́a cinética es:
Fx U0
ε= (2.2)
ρD3
Sabemos (ver Mecánica de Fluidos 1), que la fuerza arrastre es dependiente del número Reynolds
y se escribe: !
Fx U0 D
= Cd (Re) = Cd Re = (2.3)
ρU02 D2 ν
donde Cd es el coeficiente de arrastre.
Combinando las Ec. 2.2 y 2.3, tenemos:
U03
ε = Cd (2.4)
D
Es la cantidad de energı́a que se debe gastar para hacer avanzar el fluido alrededor de la esfera (o
de un perfil aerodinámico !!!). Se debe realizar experiencias para tener información sobre Cd (Re).
La Figura 2.5 presenta el coeficiente de arrastre Cd como función del número de Reynolds. Las
experiencias muestran que el coeficiente Cd es prácticamente constante cuando el número de Rey-
nolds es suficientemente grande. Este resultado no es predecible de forma teórica Alrededor
de Re ≈ 5 × 105 , tenemos un déficit de arrastre, la turbulencia contamina no solamente a la estela
sino que también las capas lı́mites de la esfera.
La potencia por unidad de masa absorbida es una fracción constante de la potencia por unidad
U3
de masa D0 que está en juego en la interacción entre la esfera de tamaño D y el flujo de velocidad
U0 .
Sobre la figura 2.6, se puede visualizar el comportamiento del flujo turbulento para varios
número de Reynolds. La estructura del flujo es más difuso cuando el número de Reynolds aumenta.
Este proceso permite adaptar la disipación en función de la potencia que impone el flujo libre.
Esto se traduce en un cambio del tamaño de las escalas de estos movimientos turbulentos: las
escalas son más pequeñas cuando la energı́a a disipar es grande.
Es a partir de esta observación que Richardson (1921) y Kolmogorov (1941) sugirieron y for-
malizaron el modelo que es capaz levantar la paradoja: la cascada de la energı́a.
La cascada de la energı́a es un proceso (disipativo) que transforma en calor una fracción
constante de energı́a transportada por el flujo (tasa de transferencia a ε =Cte), para cualquier
escala de la estructura.
Figura 2.5: Cd vs Re
54 CAPÍTULO 2. INTRODUCCIÓN A LOS FLUJOS TURBULENTOS
.
Figura 2.6: Para varios número de Reynolds, comportamiento del flujo turbulento alrededor una
esfera a) Re = 200 × 103 (arriba) b) Re = 400 × 103 (abajo).
Normalmente, el periodo de integración T debe ser infinito, pero nunca disponemos del histórico
completo de las evoluciones de una variable. En la práctica, se realiza el promedio sobre un periodo
2.2. TRATAMIENTO ESTADÍSTICO DE LOS FLUJOS TURBULENTOS 55
Por eso debemos tomar un valor de T suficientemente grande para hacer un promedio sobre las
fluctuaciones, pero suficientemente pequeño para el tiempo caracterı́stico de las variaciones lentas
de φ. La fluctuación φ0 (~x, t) se obtiene haciendo: φ(~x, t) − φ(~x).
Figura 2.7:
El promedio sobre un conjunto es muy útil sobre el plano conceptual, pero en la práctica es
generalmente inaccesible. Es difı́cil realizar N experiencias idénticas o de repetir N veces la misma
experiencia. En algunos casos, el flujo considerado es único y disponemos de una única realización
con las variables del sistema en estudio.
Estas ecuaciones promedios hacen aparecer un conjunto de momento de orden 2, los esfuerzos
turbulentos o esfuerzos de Reynolds, que representan las transferencias asociadas a las fluc-
tuaciones. La resolución de estas ecuaciones promedio necesitan de la modelación de los esfuerzos
turbulentos (Ver parágrafo 2.7.1).
1. Linealidad:
(φ1 (~x, t) + φ2 (~x, t)) = φ1 (~x) + φ2 (~x) (2.11)
∂φ(~x, t) ∂φ(~x)
= (2.14)
∂t ∂t
∂φ(~x, t) ∂φ(~x)
= i = 1, . . . , 3 (2.15)
∂xi ∂xi
2
Unicamente el promedio temporal será utilizado en todo este capı́tulo
2.2. TRATAMIENTO ESTADÍSTICO DE LOS FLUJOS TURBULENTOS 57
φ(~x, t)ψ(~x, t) = φ(~x) ψ(~x) + φ0 (~x, t) ψ(~x) + φ(~x) ψ 0 (~x, t) + φ0 (~x, t)ψ 0 (~x, t) (2.20)
Este resultado, clásico en estadı́stica, permite establecer las ecuaciones promedio de Navier-Stokes,
de la energı́a cinética . . . . Si las variables φ(~x, t) y ψ(~x, t) representan los componentes de velocidad
u(~x, t), v(~x, t) y w(~x, t), entonces se puede escribir3 :
u2 = (u)2 + u0 2 (2.22)
v 2 = (v)2 + v 0 2 (2.23)
w2 = (w)2 + w0 2 (2.24)
uv = u v + u0 v 0 (2.25)
uw = u w + u0 w0 (2.26)
vw = v w + v 0 w0 (2.27)
Todas estas expresiones se pueden escribir en una forma más compacta usando la siguiente nota-
ción:
0 0
ui vj = ui vj + ui vj (2.28)
3
Se usa una notación más liviana, es decir, se suprimió (~x, t). Por supuesto las fluctuaciones depende (~x, t) y el
valor promedio dependen de (~x).
58 CAPÍTULO 2. INTRODUCCIÓN A LOS FLUJOS TURBULENTOS
Para el campo de fluctuaciones de velocidad, es decir φ(~x, t) = u~0 (~x, t), tenemos:
T1 = T2 = T3 = T (2.39)
de la viscosidad del fluido, porque bajo un cierto tamaño, la viscosidad transforma las es-
tructuras en calor y
de la disipación ε.
La escala de Kolmogorov da una estimación del tamaño de las estructuras (vórtices), para la
cual, los efectos disipativos juegan un papel primordial. Eso se debe a que para escalas mayores,
el efecto viscoso es más débil que el efecto inercial. Estas estructuras tienen una escala (o tamaño)
lk que se conoce como escala de Kolmogorov y una velocidad uk . El tiempo caracterı́stico es:
lk
τk = (2.41)
uk
La energı́a de los vórtices de Kolmogorov es disipada por los esfuerzos viscosos en un tiempo
igual al tiempo caracterı́stico τk . Para que sea ası́, la escala lk debe ser del orden del espesor
de una capa viscosa que se desarrollará sobre un tiempo τk :
√
lk ∼ ντk (2.42)
Es decir que tenemos: s
√ lk
lk ∼ ντk = ν (2.43)
uk
que podemos escribir:
uk lk
∼1 (2.44)
ν
62 CAPÍTULO 2. INTRODUCCIÓN A LOS FLUJOS TURBULENTOS
Por lo tanto, se obtiene que el número de Reynolds de los vórtices de Kolmogorov es del orden
de la unidad. En términos de número de ondas, la escala de Kolmogorov le corresponde al valor
Kmax relacionado con la escala más pequeña lk :
2π
Kmax = (2.45)
lk
2 2u0i u0i
λ = 2 (2.47)
∂u0
∂r r=0
Por tanto no se identifica con ningún tamaño fı́sicamente relevante, aunque puede decirse que
corresponde a tamaños pequeños respecto a las escalas mayores del sistema a estudiar, sin ser tan
pequeños como la escala de Kolmogorov. Por tanto λ, se encuentra comprendida entre L y lk .
La ventaja de utilizar un número de Reynolds referido a este valor, radica en que nos permite
comparar la intensidad de la turbulencia en experimentos (numéricos o experimentales) diferentes,
pues es independiente de la geometrı́a de estos:
u0 λ
Reλ = (2.48)
ν
2. Los grandes vórtices (Macro-escala de Taylor) obtienen su energı́a del flujo promedio (velo-
cidad caracterı́stica Up y escala caracterı́stica Lp , tiempo caracterı́stico Tp ). El mecanismo
de transferencia de energı́a es eficiente cuando los tiempos caracterı́sticos del flujo promedio
y de las fluctuaciones turbulentas de gran tamaño (Macro-escala de Taylor) son similares:
Lp L
Tp = ∼ =T (2.49)
Up U
3. La energı́a del movimiento (Macro-escala de Taylor) es transferido hacia las pequeñas es-
tructuras desde las más grandes estructuras de escala L, de velocidad U y de tiempo de vida
T donde actúan las fuerzas exteriores. La transferencia se realiza con una tasa:
U2 U3
ε0 ∼ ∼ (2.50)
T L
2.5. EL MODELO DE LA CASCADA DE ENERGÍA: CASCADA DE KOLMOGOROV 63
4. La energı́a es disipada, a la misma tasa ε, al nivel de las estructuras que son suficientemente
pequeñas para poder ser eliminadas por el acción de la fricción viscosa. Por definición, las
escalas (lk , uk ) de estas estructuras verifican:
lk uk
∼1 (2.51)
ν
Es al nivel del régimen viscoso donde se realiza la transformación de energı́a cinética en
calor (efecto Joule).
u2l u3
εl ∼ ∼ l = ε0 = Cte (2.52)
τl l
es constante para cualquier l. Las propiedades dinámicas de las estructuras de la turbulencia
en el régimen inercial depende solamente de l y ε0
lk l Lp (2.53)
64 CAPÍTULO 2. INTRODUCCIÓN A LOS FLUJOS TURBULENTOS
Figura 2.10: Jet turbulento. a) bajo número de Reynolds. b) Alto número de Reynolds. Cambio
de textura entre los 2 casos.
La energı́a es transferida en cascada, es decir, de las estructuras más grandes hacia las
estructuras las más finas para ser disipada al nivel de los esfuerzos viscosos.
La cascada de Kolmogorov ofrece una descripción intuitiva de la repartición de las escalas en los
flujos turbulentos, es adecuada en el caso donde la turbulencia es cuasi-isotrópica. Los resultados
obtenidos sobre esta base son fácilmente usados para estimar las escalas tı́picas de flujos
turbulentos complejos.
Las diferentes escalas que fueron introducidas son resumidas en la tabla 2.11:
U Lp
Macro (grande vortice) Lp U T = Lp /u ReL = ν
donde εk es la tasa de disipación promedio. Una análisis dimensional muestra que la única posibi-
lidad es:
2/3
E(K) ∝ εk K −5/3 (2.68)
Esta ley es representada sobre la Figura 2.12.
Para K η, es decir para vórtices muy pequeños, el espectro decrece muy rápidamente, la
viscosidad elimina todos los movimientos de tamaño inferior.
Para K → 0, es decir para las grandes escalas, argumentos teóricos (más complejos que el
análisis dimensional) y las observaciones experimentales y numéricas indican que:
Este modelo describe relativamente bien la realidad, la experiencia (Ver Figura 2.13, en escala
logarı́tmica), muestra que el perfil de E(K) es único. Se observa un zona inercial E(K) ∝ K −5/3
cuando Re es suficientemente grande. El tamaño de la zona inercial depende de Re. La determina-
ción de la forma completa del espectro es un campo de investigación abierto. Numerosos modelos
fueron propuestos desde los años 1940.
Finalmente se puede encontrar la constante de proporcionalidad del espectro de energı́a E(k):
Figura 2.13: Perfil de E(K) obtenido experimentalmente para varios flujos: capa lı́mite, turbulencia
de rejilla, canales, . . .
1 1 Z +∞
k(t) = u0i u0i (t) = Rii (t) = E(K, t) dK (2.71)
2 2 0
El análisis de E(K, t) permite conocer la repartición de la energı́a cinética fluctuante entre las
diferentes escalas dinámicamente activas. Entonces, es posible saber cuales son los movimientos
turbulentos más intensos.
E(K)dK representa la energı́a de todos los vórtices con una longitud de onda entre 2π K
2π
y K+dK
2.7. MODELACIÓN DE LA TURBULENCIA 69
Z +∞ Z +∞
εk = D(K) dK = 2νK 2 E(K) dK (2.72)
0 0
2.6.3. Comentarios
La teorı́a de Kolmogorov está establecida sobre la consideración de una transferencia de ener-
gı́a homogénea en espacio. Solamente toma en cuenta un valor promedio de la disipación. Al
contrario, el mecanismo de estiramiento de la vorticidad indica que esta transferencia es principal-
mente realizada por objetos estirados, que llenan una pequeña parte del espacio. El mecanismo de
estiramiento parece poner en duda el enfoque de Kolmogorov. Además, recientemente fue demos-
trado que existe una transferencia de energı́a también desde las pequeñas escalas hacia las grandes
escalas.
2.7.1.1. Introducción
Consideramos los flujos turbulentos de fluidos incompresibles y newtonianos. Las ecuaciones
del movimiento en este caso se escriben como:
∂u ∂v ∂w
+ + =0 (2.73)
∂x ∂y ∂z
∂ ∂ 2 ∂ ∂ ∂p
ρu + ρu + ρuv + ρuw = − + µ∆u + ρgx (2.74)
∂t ∂x ∂y ∂z ∂x
∂ ∂ ∂ 2 ∂ ∂p
ρv + ρuv + ρv + ρvw = − + µ∆v + ρgy (2.75)
∂t ∂x ∂y ∂z ∂y
∂ ∂ ∂ ∂ ∂p
ρw + ρuw + ρvw + ρw2 = − + µ∆w + ρgz (2.76)
∂t ∂x ∂y ∂z ∂z
Remplazamos las variables u, v, w, p utilizando la descomposición de Reynolds para cada una
de las variables:
u(~x, t) = u(~x, t) + u0 (~x, t) (2.77)
Tomando el promedio de la Ec. 2.81 y usando las propiedades presentadas en la parte 2.2.5.1,
tenemos:
∂u(~x, t) ∂v(~x, t) ∂w(~x, t) ∂u0 (~x, t) ∂v 0 (~x, t) ∂w0 (~x, t)
+ + + + + =0 (2.82)
∂x ∂y ∂z ∂x ∂y ∂z
Como los promedios de las fluctuaciones son nulas:
∂u ∂v ∂w
+ + =0 (2.84)
∂x ∂y ∂z
Restando la Ec.2.84 de la Ec. 2.81
∂u0 ∂v 0 ∂w0
+ + =0 (2.85)
∂x ∂y ∂z
que constituye la ecuación de continuidad para las fluctuaciones
2.7. MODELACIÓN DE LA TURBULENCIA 71
DV~ −−→
ρ = ρf~ − gradp + div~τ vis (2.94)
Dt
De la misma manera se puede escribir la ecuación de cantidad de movimiento para el campo
promedio:
DV~ −−→
ρ = ρf~ − gradp + div~τ vis + div~τ turb (2.95)
Dt
Ası́ la ecuación de cantidad de movimiento promedio según el eje x se escribe:
∂ ∂ 2 ∂ ∂
ρ(u) + ρu + (ρu v) + (ρu w)
∂t ∂x ∂y ∂z
∂
= − p + ρgx
∂x
∂ vis ∂ vis ∂ vis
+ τxx + τxy + τxz
∂x ∂y ∂z
| {z }
esf uerzos viscosos
∂ turb ∂ turb ∂ turb
+ τxx + τxy + τxz (2.96)
∂x ∂y ∂z
| {z }
esf uerzos turbulentos
∂ ∂ ∂p ∂ vis ∂ turb
ρvi + ρvi vj = − + ρgi + τij + τ (2.100)
∂t ∂xj ∂xi ∂xj ∂xj ij
donde:
El tensor de los esfuerzos de Reynolds o de los esfuerzos turbulentos está representado por:
0 0
τijturb = −ρ (ui uj − ūi ūj ) = −ρvi vj (2.102)
Comentarios Los esfuerzos de Reynolds aparecen en las ecuaciones anteriores como parámetros
no conocidos suplementarios. El problema central del análisis estadı́stico de los flujos
turbulentos consiste en la modelación de estos esfuerzos de Reynolds. En efecto, para
resolver las ecuaciones que representan el flujo promedio se debe establecer un modelo que permita
el cálculo de los esfuerzos turbulentos a partir de las variables del flujo promedio.
No hay una representación única de los esfuerzos de Reynolds en función del flujo promedio
y las hipótesis de cierre (o los modelos de turbulencia) usados en la práctica son muy
variados. Los modelos de turbulencia se presentan en general de la forma de un conjunto
de ecuaciones en las cuales aparecen los esfuerzos turbulentos.
Este conjunto de ecuaciones completa las ecuaciones dinámicas promedio y debe permitir la reso-
lución simultánea de todas las ecuaciones del movimiento.
turb 0 0
τ12 = −ρv1 v2 (2.103)
t t 0 0
Los otros esfuerzos de Reynolds (τ11 = −ρv12 , τ22 = −ρv22 , . . . ) son despreciables o no par-
ticipan de la dinámica de la mezcla. La analogı́a con los mecanismos que producen los esfuerzos
viscosos conduce a escribir:
∂v1 turb 0 0
τ12 = −ρv1 v2 = µt (2.104)
∂x2
turb
A partir de lo anterior, el esfuerzo τ12 quedarı́a proporcional al gradiente transversal de
velocidad. El coeficiente de proporcionalidad µt , a diferencia de la viscosidad laminar, depende
directamente del flujo considerado. El problema se reduce entonces a determinar µt . Las
posibilidades para elegir µt son numerosas, pero todas tienen un campo de aplicación limitado.
Un método particularmente simple consiste en suponer que µt mantiene un valor constante en
cada sección de la capa de mezcla (Ver Figura 2.15). Una análisis dimensional muestra que:
µt ∼ ρU l (2.105)
Por ejemplo, se podrı́a tomar para U la diferencia de velocidad y para l el espesor local δ de la
capa de mezcla. Sea:
µt = Kρ(U1 − U2 )δ (2.106)
La constante K puede ser ajustada por comparación entre las soluciones teóricas y los resultados
experimentales.
con:
2.7. MODELACIÓN DE LA TURBULENCIA 75
∂ Ūx
ttxy = µt (2.113)
∂y
Ası́ entonces el primer y más simple modelo de turbulencia fue introducido por Boussinesq en
1877 y consistió en elegir una viscosidad turbulenta µt (o νt ) constante.
νt = Cte (2.114)
Este modelo, toma en cuenta las propiedades locales del flujo. Puede dar buenos resultados cuando
los campos turbulentos son homogéneos en espacio. Por ejemplo, en el centro de un conducto. No
es el caso de los flujos cerca de fronteras, por ejemplo el fondo de un rı́o o cerca de las paredes de
los tubos. En estos casos se necesita modelos de turbulencia más complejos.
la distancia a la pared.
∂v
2 1
µt = ρlm (2.115)
∂x2
Este método de modelación que utiliza el término de la longitud de mezcla da buenos resultados
para configuraciones simples.
lm = Ky (2.117)
Figura 2.16: (a) Perfil de una capa lı́mite turbulenta. (b) Longitud de mezcla lm según distancia
vertical x2 .
Muy muy cerca (esta distancia será definida posteriormente en el Capitulo 4) de la pared
(Ver Figura 2.16), se puede definir una longitud de mezcla proporcional a la distancia x2 de
la pared:
lm = κx2 donde κ ≈ 0, 41 (2.118)
Figura 2.17:
δ(x)
lm ≈ (2.119)
10
Después de los modelos básicos de Prandtl y Boussinesq se presentarán modelos más desa-
rrollados los modelos con ecuaciones de transporte. Estos modelos son más apropiados cuando el
numero de Reynolds es elevado y se cuenta con geometrı́as complejas.
Escribimos la Ec. 2.123 para 2 ı́ndices diferentes i (sea Ei ) y k (sea Ek ) y formamos las 2
combinaciones:
u0i Ek + u0k Ei
obtenemos después un cálculo no inmediato una ecuación para la correlación u0i u0j
a) términos de difusión.
b) términos de producción.
Términos de difusión:
• estos términos se caracterizan por una simple derivada de los parámetros en una direc-
ción del espacio. Si el flujo es adyacente a una zona no turbulenta, su contribución es
nula en el balance de la ecuación.
Términos de Producción:
• Estos términos son positivos o negativos o también conocidos como términos de fuente
o términos de sumidero, respectivamente.
• La producción es principalmente una contribución del flujo promedio, que interviene
por sus derivadas ∂u k
∂xj
, es decir por su cizalle, este representa la parte esencial de la
producción de la turbulencia.
2.8. MODELACIÓN DE LA TURBULENCIA ISOTRÓPICA 79
νt = νt (k) (2.125)
Ası́, es necesario escribir la ecuación de la energı́a cinética turbulenta k.
Para escribir una ecuación de k, tomamos i = k en la Ec. 2.96.
" #
∂k ∂k ∂ 1 0 0 0 1 ∂ui
+ uj =− (ui uj uk ) + (p0 u0j ) − (u0i u0j ) + ν∆k − ε (2.126)
∂t ∂xj ∂xj 2 ρ ∂xj
Se puede escribir:
∂k ∂k ∂ 2k
+ ūj = Diff t + P + ν −ε (2.127)
∂t ∂xj ∂xj ∂xj
80 CAPÍTULO 2. INTRODUCCIÓN A LOS FLUJOS TURBULENTOS
con
∂u0i ∂u0i
" #
∂ 1 0 0 0 1 0 0 ∂ ūi
Diff t = − uuu + up P = −u0i u0j ε=ν (2.128)
∂xj 2 i i j ρ j ∂xj ∂xj ∂xj
ε es la disipación viscosa (m2 .s−3 ) de la energı́a cinética turbulenta por las fricciones viscosas
(con producción de calor). Este término es siempre positivo.
∂k ∂kūj ∂ 2k
+ = P + Diff t + ν −ε (2.131)
∂t ∂xj ∂xj ∂xj
∂u0i ∂u0i
!
εij = ν (2.132)
∂xj ∂xj
νt = νt (k, ε) (2.133)
2.8. MODELACIÓN DE LA TURBULENCIA ISOTRÓPICA 81
Las ecuaciones que damos aquı́ son válidas solamente lejos de las paredes (cerca de las paredes se
debe usar leyes de pared o modelo de amortiguamiento de la viscosidad turbulenta a la pared):
" #
∂k ∂k ∂ ∂k
+ ūj =P −ε+ (ν + νt ) (2.134)
∂t ∂xj ∂xj ∂xj
" #
∂ε ∂ε ε ∂ νt ∂ε
+ ūj = (Cε1 P − Cε2 ε) + ν+ (2.135)
∂t ∂xj k ∂xj σε ∂xj
con
k2
νt = C µ (2.136)
ε
y
P = νt S 2 (2.137)
q
donde S = 2Sij Sij es la tasa de deformación escalar asociada al campo de velocidad promedio. El
tensor de deformación Sij es escribe Sij = 21 ∂x ∂ ūi
j
+ ∂∂xūji . Se observa que la ecuación de disipación
se inspira de la ecuación de la energı́a cinética turbulenta. Las constantes numéricas del modelo
Cµ , Cε1 , Cε2 , σε son
Figura 2.18:
Capı́tulo 3
Re → ∞ (3.1)
Esta situación permito ver paradojas que no puede ser explicado con la hipótesis de flujo
ideal. Los flujos que consideramos ahora son reales, es decir viscosos, incompresibles, a número
de Reynolds grande pero finito. De un punto práctico, estas condiciones permiten describir una
amplia variedad de situaciones aerodinámica, hidrodinámica, . . .
Para un fluido newtoniano, la respuesta, es a buscar en las regiones del flujo donde existe un gran
gradiente de velocidad.
Cuando el número global de Reynolds es importante, la variación espacial de velocidad entre el
obstáculo y el fluido al infinito no afecta de manera significativa la totalidad del campo del flujo,
pero se concentra en regiones particulares de fuertes gradientes.
Estas regiones son localizadas a proximidad inmediata de la pared y en la estela aguas arribas
del obstáculo. Corresponden a las zonas de influencia de la viscosidad.
83
84 CAPÍTULO 3. CAPA LIMITE LAMINAR DINÁMICA
Mas el numero de Reynolds global crece, mas el espesor relativo de la zona de influencia de la
viscosidad decrece. La región parietal, donde se manifiesta la influencia de la viscosidad del fluido
a gran numero de Reynolds, es llamada capa limite dinámica (Ver Figura 3.1). Los trabajos de
Prandtl (1904) permitieron una gran avance en el entendimiento del comportamiento de la capa
limite dinámica.
Además, vimos que el punto de separación de la capa limite varia según el régimen del flujo o
número de Reynolds (Ver Figura 3.2).
el otro cerca de la pared donde la cantidad de movimiento .es transportada.en raı́z del tiempo
por la viscosidad cinemática.
Además, una descripción espacial permite ver dos escalas de longitud muy distintas: una grande
en el sentido del flujo y una pequeña en el sentido transversal a la pared.
Todo flujo real (es decir viscoso) de velocidad ~v en contacto con un cuerpo impermeable móvil
(cualquier cuerpo) con una velocidad w ~ debe satisfacer a las siguientes condiciones:
~v .~n = w.~
~n (3.2)
Figura 3.2: Despegamiento de la capa limite alrededor de un cilindro en función del régimen del
flujo.
3.2.2. Ecuaciones
Para poner en evidencia el fenómeno de la capa limite, se considera las ecuaciones de conserva-
ción de la masa y de la cantidad de movimiento y se considera despreciable las fuerzas de volumen
(número de Froude grande). Sea de manera adimensional:
div~v = 0 (3.6)
La ecuación de Euler es del primer orden mientras que la ecuación de NS son del segundo
orden. Esto implica que una de las condiciones de frontera no sera satisfecha. La ecuación
modela los fluidos ideales, sin viscosidad, por lo que no se puede tener condición de adhe-
rencia.
Consideramos las fuerzas de inercia. Por unidad de volumen estas fuerzas de la forma:
∂u
ρu (3.10)
∂x
∂u
Para un placa plana de largo L, el gradiente de velocidad ∂x
puede ser estimado por:
∂u U∞
∼ (3.11)
∂x L
3.2. PRIMERAS ECUACIONES DE CAPA LIMITE DINÁMICA 87
Figura 3.3: Esquema de una capa lı́mite que se desarolla sobre un placa plana.
En esta condición, las fuerzas de inercia por unidad de volumen tiene por orden de magnitud:
∂u U2
ρu ∼ρ ∞ (3.12)
∂x L
Consideramos ahora, las fuerzas asociadas al esfuerzo cortante viscoso. Por unidad de volumen,
estas fuerzas se escriben:
∂τxy
(3.13)
∂y
En flujo laminar y por un fluido newtoniano y incompresible, este termino se expresa
∂ 2u
µ (3.14)
∂y 2
Si δ es el espesor de la capa limite, una estimación de las fuerzas viscosas por unidad de volumen
es:
∂τxy U∞
∼µ 2 (3.15)
∂y δ
En la capa limite, las fuerzas de inercia y de viscosidad se equilibran. Se puede escribir:
2
U∞ U∞
ρ ∼µ 2 (3.16)
L δ
Sea: s s
µL µ
δ∼ =L (3.17)
ρU∞ ρU∞ L
Ası́, el ratio Lδ es inversamente proporcional a la raı́z cuadrada del número de Reynolds basado
sobre el largo de la placa: s
δ µ 1
∼ =√ (3.18)
L ρU∞ L ReL
La ecuación:
δ 1
∼√ (3.19)
L ReL
implica que: el espesor relativa de la capa limite Lδ decrece como ReL −1/2 , El espesor δ varia como
la raı́z cuadrado de L y la frontera de la capa limite tiene la forma parabólica.
88 CAPÍTULO 3. CAPA LIMITE LAMINAR DINÁMICA
• Automóvil: δ ≈ 0, 9mm
• Avión: δ ≈ 1, 2mm
• Bote: δ ≈ 1, 3mm
y en estas condiciones: s s
U∞ ρU∞ 2 µ
τp ∼ µ = µU∞ = ρU∞ (3.23)
δ µL ρU∞ L
Ası́ se puede escribir que:
τp 1
2
∼√ (3.24)
ρU∞ ReL
También la resolución exacta de las ecuaciones de capa limite dinámica permite encontrar el
coeficiente de proporcionalidad. En el caso de una capa limite dinámica sobre una placa plana y
para un flujo laminar incompresible, tenemos:
τp 0, 332
2
= √ (3.25)
ρU∞ ReL
FD
CD = 1 2 A
(3.27)
2
ρU∞
Se obtiene fácilmente a partir de la Ec. 3.26 que:
1
CD ∼ √ (3.28)
ReL
La resolución exacta de las ecuaciones de capa limite dinámica permite encontrar el coeficiente de
proporcionalidad. En el caso de una capa limite dinámica sobre una placa plana y para un flujo
laminar incompresible, tenemos:
1, 328
CD = √ (3.29)
ReL
3.2.6. Ejemplo
Una placa plana de largo 1m se encuentra posicionada en un flujo de aire con una velocidad
U = 5 m.s−1 . La viscosidad dinámica del aire es 1, 76 × 10−5 Pa.s a la temperatura de 20◦ C y su
masa especifica ρ = 1, 2kg.m−3 .
Notamos ue (x) el valor local de la velocidad obtenida en fluido ideal en el punto de abscisa
x sobre el obstáculo.
u(x, δ(x))
= 0, 99 (3.30)
ue (x)
donde ue (x) es el valor local de la velocidad obtenida en fluido ideal en el punto de abscisa x
sobre el obstáculo también la velocidad al exterior de la capa limite en la abscisa x. Esta definición
(un poco arbitrario), permite determinar de manera precisa la frontera de la capa limite.
90 CAPÍTULO 3. CAPA LIMITE LAMINAR DINÁMICA
Figura 3.4: Campo de velocidad alrededor de un cilindro: a) para un fluido ideal b) para un fluido
viscoso.
sea: !
Z δ
∗ ρu
δ = 1− dy (3.34)
0 ρ e ue
El parámetro δ ∗ es llamado: espesor de desplazamiento. Este espesor corresponde a la distancia
de desplazamiento necesaria para conservar, en fluido ideal sobre el espesor δ − δ ∗ el mismo flujo
másico que en fluido viscoso sobre toda la sección de capa limite.
Figura 3.5: Ilustración del espesor de desplazamiento. a) para un fluido ideal b) para un fluido
viscoso c) desplazamiento δ ∗ .
3.3. PARÁMETROS CARACTERÍSTICOS DE LA CAPA LIMITE DINÁMICA 91
3.3.3.1. Demostración
Tomamos despreciable, en un fluido ideal, todo gradiente transversal a la escala δ. La igualdad
entre los flujos másicos se expresa:
Z δ Z δ Z δ
ρe ue dy = ρu dy ⇐⇒ ρe ue (δ − δ ∗ ) = ρu dy (3.35)
δ∗ 0 0
Comparamos el flujo definido por le Ec. 3.37 con un flujo que habrá el mismo flujo másico ṁv
pero transportado como un fluido ideal con una velocidad ue (o con un velocidad externa ue ).
El flujo másico en la capa limite se escribe:
Z δ
ṁv = ρu dy (3.38)
0
Ji − Jv = ρe u2e θ (3.41)
Figura 3.6: Capa limite sobre placa plana. Se indica el volumen de control elegido.
ρe u2e δ (3.44)
Este balance expresa simplemente que la variación de cantidad de movimiento del fluido es
producida por las fuerzas que se ejercen sobre la superficie del volumen de control ABCD.
En ausencia del gradiente de presión, la presión es constante en todo el fluido y tenemos:
Z δ
Fw = ρu (ue − u) dy (3.49)
0
Sea:
Fw = ρu2e θ (3.50)
Ası́, si conocemos (o medimos) θ se puede estimar directamente la fuerza parietal que se ejerce
por el fluido sobre la placa entre 0 y x. Aquı́ la fuerza parietal estaba calculada por unidad de
ancho de la placa.
3.3.7.1. Ejemplo
Consideramos un perfil de velocidad dado en la forma de una ley de potencia:
1/n
u y
= 0<y<δ (3.56)
ue δ
94 CAPÍTULO 3. CAPA LIMITE LAMINAR DINÁMICA
donde n es un constante.
Podemos mostrar que para un flujo incompresible (ρ = ρe ) con este tipo de perfil, el factor de
forma H se escribe:
n+2
H= (3.57)
n
Si tenemos un perfil lineal n = 1 entonces H = 3.
Si tenemos un perfil laminar real n1 = 0, 8 entonces H = 2, 59 (este resultado será demostrado)
El valor n ≈ 7 que corresponde a un perfil de velocidad muy cerca al perfil que se encuentra
en las capas limites turbulentas. Por este valor H = 1, 29
Vamos a ver más adelante que cuando la capa limite despega: el perfil de velocidad presenta un
punto de inflexión a la pared y el factor de forma tiene un valor H ≈ 4. En esta situación te-
nemos un criterio de despegamiento basado sobre la evaluación del factor de forma H. Finalmente,
en ausencia de capa limite H = 1
En una región cerca del borde de ataque el fluido exterior (fluido ideal) el movimiento es
acelerado ( du
dx
e
> 0) y el gradiente longitudinal de presión es negativo ( dP
dx
e
< 0) (Bernoulli).
dPe
Si dx < 0 el gradiente de presión es favorable (Tenemos una presión menor aguas abajo.).
Una partı́cula encuentra aguas abajo solamente una resistencia provocada por las fuerzas de
viscosidad.
Según la incidencia y/o la curvatura del perfile, un sentido opuesto del gradiente de presión
puede ser observado mas allá aguas abajo, En este caso, la deceleración del fluido por la
viscosidad se puede conjugar a una presión creciente dP dx
e
> 0 cerca del borde fuga. La
fuerza de inercia debe entonces equilibrar 2 fuerzas de resistencia al movimiento: presión y
viscosidad que actúan en el mismo sentido. Si dP
dx
e
> 0 el gradiente de presión es desfavorable.
3.4. DESPEGAMIENTO DE LA CAPA LIMITE 95
Figura 3.7: Flujo de un fluido alrededor de un cuerpo. Ilustración del despegamiento de la capa
limite
Figura 3.9: Ejemplos de despegamiento de capa limite sobre una superficie plana (arriba) y un
perfil alar (abajo).
!
∂u ∂u
= 0 y τp = µ =0 (3.58)
∂y y=0
∂y y=0
3.4.4.1. El cilindro
En el caso de un cilindro posicionado en un flujo uniforme (Ver Figura 3.10) se observa (para
un rango de valor del número de Reynolds), la generación de vórtices alternados (Calle de von
Karman). La caracterı́stica de este proceso depende del número de Reynolds. La frecuencia de
generación de los vórtices es generalmente determinada por una ley de Strouhal. Por ejemplo, en
el caso de cilindro circular posicionados en un flujo uniforme el numero de Strouhal:
fD
≈ 0, 2 (3.59)
ν
3.5. TRANSICIÓN HACIA LA TURBULENCIA 97
3.4.4.2. El difusor
El despegamiento se produce también en otros flujos, por ejemplo: en los conductos fuertemente
divergentes con en un difusor (Ver Figura 3.11). Antes, la contracción, la presión decrece, el flujo es
acelerado y la capa limite adhiere a la pared del convergente. Después el convergente, el gradiente
de presión es positivo, el flujo decelera y la capa limite se aleja de la pared. El despegamiento
es generalmente asimétrico con la formación de una zona de retorno cerca de la pared, el flujo
principal siguiendo la otra pared.
Para evitar, el despegamiento, se usa en práctica ángulos de divergentes relativamente pequeños
(< 15◦ ). Para limitar el despegamiento, se puede usar baffles. Esta herramienta divide el divergente
Figura 3.11: Despegamiento de la capa limite en un difusor: a) flujo simétrico b) flujo asimétrico.
√
δ Re
Figura 3.12: Espesor de la capa limite, x
, en función de número de Reynolds. Transición
laminar-turbulenta.
Cuando el flujo es turbulento, las fluctuaciones turbulentas en la capa limite intensifican los
intercambios de cantidad de movimiento. La pendiente del perfil de velocidad cerca de la pared
crece y genera un crecimiento de los esfuerzos parietales. La capa limite es turbulenta.
El espesor de la capa limite dinámica turbulenta sobre una placa plana puede ser estimada
usando la siguiente relación:
0,37x
δturbulento = δ = 5 × 105 ≤ Rex ≤ 107 (3.61)
Re1/5
x
Hacemos una estimación de los términos que aparecen en las Ec. 3.65, 3.66, 3.67. Para esta esti-
mación, supongamos que:
100 CAPÍTULO 3. CAPA LIMITE LAMINAR DINÁMICA
Los términos:
(1) y (2) son del mismo orden de magnitud
(4) es mas pequeño que el (5), en efecto:
!2
4 δ
∼ 1 (3.74)
5 L
Se puede tomar despreciable el termino (4) pero se debe guardar el (5) sino los efectos de la
viscosidad serán despreciables en el problema.
!
U2 Uδ U 1 ∂p U U
ν ν (3.75)
L
|{z} L δ ρ ∂x L2}
| {z δ2
|{z}
| {z } | {z }
1 2 3 4 5
Ası́, los dos primeros términos deben ser del mismo orden de magnitud que el termino (5):
U2 U
(1) ∼ (2) ∼ (5) ⇒ ∼ν 2 (3.76)
L δ
3.6. ECUACIONES DE LA CAPA LIMITE LAMINAR: MODELO DE PRANDTL 101
Sea:
L
δ∼√ (3.78)
Re
Ası́ para que Lδ 1 se necesita Re 1. Una consecuencia de la relación 3.78 no es válida cerca
del borde de ataque. En todo rigor, se debe usar las ecuaciones de Navier-Stokes completa en esta
región.
Consideramos la ecuación de cantidad de movimiento según el eje y: Ec.3.67
!
∂v ∂v 1 ∂p ∂ 2v ∂ 2v
u +v =− +ν + (3.79)
∂x ∂y ρ ∂y ∂x2 ∂y 2
1 ∂p U2
∼ (3.82)
ρ ∂x L
Ası́, vemos que el gradiente de presión transversal es inferior al gradiente de presión axial:
1 ∂p
ρ ∂y δ
1 ∂p
∼ 1 (3.83)
ρ ∂x
L
En la capa limite, la presión transversal varia muy poco. Ası́, en la capa limite, la presión
transversal varia muy poco y podemos escribir:
1 δp
=0 (3.84)
ρ δy
La presión p depende solo de x y entonces las derivadas parciales son reemplazadas por deri-
vadas:
due 1 dp
ue =− (3.86)
dx ρ dx
La integración de la Ec. 3.86 indica simplemente que la presión y la velocidad al exterior de la
capa limite están ligadas por la ecuación de Bernoulli:
1
p + ρu2e = Cte (3.87)
2
En el flujo interior:
∂u ∂v
+ =0 (3.88)
∂x ∂y
∂u ∂u 1 ∂p ∂ 2u
u +v =− +ν 2 (3.89)
∂x ∂y ρ ∂x ∂y
∂p
=0 (3.90)
∂y
En el flujo exterior:
due 1 dp 1
ue =− o p + ρue = Cte (3.91)
∂x ρ dx 2
Con las condiciones de frontera:
Figura 3.14: Capa limite sobre una placa plana: a) perfiles de velocidad en 2 secciones b) perfiles
de velocidad en forma adimensional Uu(x,y)
e (x)
y
en función de g(x) .
U (x, y) y
= Φ(η) con η= (3.93)
ue (x) g(x)
Varias verificaciones experimentales permitieron verificar y confirmar este resultado. Se toma esta
propiedad como hipótesis de resolución de las ecuaciones de Prandtl, llamada condición de afinidad.
En resumen: Existen situaciones donde es legitimo buscar una clase particular de funciones
solución de las ecuaciones de Prandtl que satisface la hipótesis de afinidad, es decir, tal que el perfil
adimensional de velocidad Uu(x,y)
e (x)
depende solamente de las variables x e y en el reagrupamiento
y
g(x)
. Usaremos entonces para la solución de Blasius:
U (x, y) y
= Φ(η) con η= (3.94)
ue (x) δ(x)
ue = Cte (3.96)
Para resolver las ecuaciones es cómodo introducir la función corriente Ψ(x, y):
∂Ψ ∂Ψ
u= v=− (3.100)
∂y ∂x
La ecuación de conservación de la masa
∂u ∂v
+ =0 (3.101)
∂x ∂y
es automáticamente satisfecha.
! !
∂ ∂Ψ ∂ ∂Ψ
+ − =0 (3.102)
∂x ∂y ∂y ∂x
La ecuación de la cantidad de movimiento
∂u ∂u ∂ 2u
u +v =ν 2 (3.103)
∂x ∂y ∂y
hace aparecer derivadas parciales de orden 3:
∂Ψ ∂ 2 Ψ ∂Ψ ∂ 2 Ψ ∂ 3Ψ
− = ν (3.104)
∂y ∂x∂y ∂x ∂y 2 ∂y 3
3.7. SOLUCIÓN EXACTA DE LAS ECUACIONES DE PRANDTL: SOLUCIÓN DE BLASIUS105
se escriben:
∂Ψ ∂Ψ ∂Ψ
(x, 0) = 0 (x, 0) = 0 (x, ∞) = ue (x) = Cte (3.106)
∂y ∂x ∂y
La condición:
∂Ψ
(x, 0) = 0 (3.107)
∂x
puede ser transformada notando:
Z
∂Ψ
(x, 0) dx = Ψ(x, 0) − Ψ(0, 0) = 0 (3.108)
∂x
Si elegimos:
Ψ(0, 0) = 0 ⇒ Ψ(x, 0) = 0 (3.109)
La elección de la condición Ψ(0, 0) = 0 no modifica el grado de generalidad
donde Φ es una función universal para todas las abscisas x y ue una constante.
Sabemos que s
νx νx 1/2
δ∼ = (3.113)
ue ue
Ası́ se puede usar como variable de similitud:
1/2
ue
η=y (3.114)
νx
Sabemos que la velocidad axial se escribe:
y ∂Ψ
u(x, y) = ue Φ u(x, y) = (3.115)
δ ∂y
entonces la función de corriente Ψ(x, y) debe calcula según:
y
Z
Ψ(x, y) = ue Φ dy (3.116)
δ
106 CAPÍTULO 3. CAPA LIMITE LAMINAR DINÁMICA
es decir:
Ψ(x, y) = ue δf (η) (3.117)
La función f que aparece en esta expresión debe ser tal que:
f 0 (η) = Φ (3.118)
Finalmente, como: s
νx
δ∼ y Ψ(x, y) = ue δf (η) (3.119)
ue
Blasius buscó perfiles autosimilares de la forma:
con:
ue 1/2
η=y (3.121)
νx
Tomamos la ecuación de cantidad de movimiento en forma de función de corriente:
∂Ψ ∂ 2 Ψ ∂Ψ ∂ 2 Ψ ∂ 3Ψ
− = ν (3.122)
∂y ∂x∂y ∂x ∂y 2 ∂y 3
Para resolver la Ec. 3.122, se debe determinar las derivadas sucesivas de Ψ en función de x y
y, con:
!
ue 1/2
1/2 1/2
Ψ(x, y) = (ue xν) f (η) = (ue xν) f y (3.123)
νx
Tenemos: 1/2
ue
η=y (3.124)
νx
Ası́: 1/2
∂η 1 ue 1η
=− y =− (3.125)
∂x 2 νx3 2x
1/2
∂η ue η
= = (3.126)
∂y νx y
Además:
∂Ψ 1 ue ν 1/2 ∂η
= f + (ue νx)1/2 f 0
∂x 2 x ∂y
1/2
1 ue ν
= (f − ηf 0 ) (3.127)
2 x
∂Ψ ∂η
= (ue νx)1/2 f 0
∂y ∂y
0
= ue f (3.128)
∂ 2Ψ ∂η
= ue f 00
∂x∂y ∂x
1η
= − ue f 00 (3.129)
2x
3.7. SOLUCIÓN EXACTA DE LAS ECUACIONES DE PRANDTL: SOLUCIÓN DE BLASIUS107
∂ 3Ψ
1/2
ue ∂η
3
= ue f 000
∂y νx ∂y
u e
= ue f 000 (3.131)
νx
Sustituimos todas las relaciones obtenidas anteriormente en la ecuación:
∂Ψ ∂ 2 Ψ ∂Ψ ∂ 2 Ψ ∂ 3Ψ
− = ν (3.132)
∂y ∂x∂y ∂x ∂y 2 ∂y 3
Tenemos:
u2e 000
1/2 1/2
1 η 1 ue ν ue
ue f 0 − ue f 00 − (f − ηf 0 ) ue f 00 =ν f (3.133)
2 x 2 x νx νx
Después simplificación se encuentra la ecuación de Blasius:
2f 000 + f f 00 = 0 (3.134)
∂Ψ ∂Ψ ∂Ψ
(x, 0) = 0 (x, 0) = 0 (x, ∞) = ue (x) (3.141)
∂y ∂x ∂y
Buscar soluciones auto-similares
Ψ(x, y) = (ue xν)1/2 f (η) (3.142)
ue 1/2 y q
η=y = Rex (3.143)
νx x
2f 000 + f f 00 = 0 (3.144)
f (0) = f 0 (0) = 0 f 0 (∞) = 1 (3.145)
El campo de velocidad es dado por
3.8.1.1. La función f 0
1/2
El comportamiento de la función f 0 en función del parámetro η = y νx
ue
es presentado en
0
la Figura 3.17. La función f permite determinar u(x, y) cuando se conocer η y ue :
u(x, y) = ue f 0 (η) (3.150)
Se puede notar que f 0 (η) varia casi linealmente cuando η → 0. Además f 0 (η) → 1 cuando η → ∞.
Se observa que tenemos f 0 (η) = 0, 99 cuando η = 5 es decir que en este punto u(x,y)
ue
= 0, 99.
El espesor de la capa limite δ definido como la distancia a la cual u(x, y) tiene 99 % de la
velocidad del flujo externo ue , entonces:
1/2
νx
δ=5 (3.151)
ue
es decir
δ 5
=√ (3.152)
x Rex
3.8. PRINCIPALES RESULTADOS OBTENIDOS A PARTIR DE LA SOLUCIÓN DE BLASIUS109
1/2
Figura 3.16: Comportamiento de las funciones f , f 0 y f 00 en función del parámetro η = y ue
νx
.
1/2
Figura 3.17: Comportamiento de la función f 0 en función del parámetro η = y ue
νx
.
3.8.1.2. La función f
1/2
El comportamiento de la función f 0 en función del parámetro η = y νx
ue
es presentado en la
Figura 3.18. El conocimiento de la función f permite el componente normal de velocidad, v(x, y),
puede ser deducido de f (η):
1/2
∂Ψ 1 ν
v(x, y) = − = ue (ηf 0 − f ) (3.153)
∂x 2 ue x
110 CAPÍTULO 3. CAPA LIMITE LAMINAR DINÁMICA
Si escribimos el producto:
1/2
v ue x vq 1
= Rex = (ηf 0 − f ) (3.154)
ue ν ue 2
Se observa cuando η → ∞ que la función 21 (ηf 0 − f ) es una constante. En efecto, 1
2
(ηf 0 − f ) →
0, 86 cuando η → ∞.
Ası́, sobre la frontera de la capa limite podemos escribir:
0, 86ue
v(x, δ) = ve = √ (3.155)
Rex
El aumento del espesor de la capa limite induce un desplazamiento del fluido en la dirección
transversal.
1/2
ue
Figura 3.18: Comportamiento de la función f en función del parámetro η = y νx
.
3.8.1.3. La función f 00
El comportamiento de la función f 00 es presenta en la Figura 3.16.
δ1∗ = 1, 72 (3.162)
Ası́, encontramos una expresión para determinar el espesor de desplazamiento de una capa limite
laminar sobre una placa plana:
1, 72x
δ∗ = √ (3.163)
Rex
Z η∞
f 02 dη = f (η∞ )f 0 (η∞ ) − f (0)f 0 (0) + 2f 00 (η∞ ) − 2f 00 (0)
0
= f (η∞ ) − 2f 00 (0) (3.171)
Finalmente tenemos: Z η∞
θ1 = f 0 (1 − f 0 ) dη = 2f 00 (0) (3.172)
0
δ∗
H= = 2, 59 (3.174)
θ
τw 2f 00 (0)
Cf = 1 2 = √ (3.179)
2
ρue Rex
Conocemos el valor de f 00 (0) = 0, 332. Ası́, tenemos un valor exacto del coeficiente de fricción:
0, 664
Cf = √ (3.180)
Rex
3.8. PRINCIPALES RESULTADOS OBTENIDOS A PARTIR DE LA SOLUCIÓN DE BLASIUS113
donde: 1/2
ue
τw = µue f 00 (0) (3.182)
νx
Ası́, debemos evaluar:
Z L
ue 1/2
00
FD = µue f (0) dx (3.183)
0 νx
La integración conduce la expresión de la fuerza de arrastre:
1/2
ue L
00
FD = 2µue f (0) (3.184)
ν
El coeficiente de arrastre está dado por:
FD 4f 00 (0)
CD = 1 = 1/2 (3.185)
ρLu2e
2 ue L
ν
Finalmente:
1, 328
CD = √ (3.186)
ReL
Esta expresión obtenida por la primera vez por Blasius es válido para números de Reynolds
ReL < 5 × 105 . Para Re > 5 × 105 , la capa limite es turbulenta y la fricción a la pared es
mucho mas importante.
3.8.8. Resumen
Resultados principales par una capa limite laminar sobre una placa plana.
Figura 3.19: Distribución de velocidad comparación entre los resultados experimentales y numéri-
cos.
δ ∗ (x) 1, 72
=√ (3.189)
x Rex
Espesor de cantidad de movimiento
θ(x) 0, 664
=√ (3.190)
x Rex
Esfuerzo parietal
τw 0, 332 τw 0, 664
2
=√ Cf = 1 =√ (3.191)
ρue Rex 2
ρu2e Rex
3.9. ECUACIONES INTEGRALES PARA LA CAPA LIMITE 115
Velocidad transversal
0, 86ue
v(x, δ) = √ (3.192)
Rex
Numero de Reynolds
xue
Re = (3.193)
ν
3.9.1.1. Uso
El método de las ecuaciones integrales pueden ser utilizado para estudiar otros problemas
de mecánica de fluido como el desarrollo de flujos libres del tipo jet, las estallas, la propagación
de llama en los conductos, . . .
Vamos a mirar las ecuaciones integrales para una capa limite en el caso de un flujo bidimensional
incompresible.
A la frontera de la capa limite, el vector velocidad tiene una componente transversal no nula.
116 CAPÍTULO 3. CAPA LIMITE LAMINAR DINÁMICA
d Zδ dδ
ve = − u dy + ue (3.201)
dx 0 dx
Además, vamos a ver después que no tenemos:
Ve dδ
≈ (3.202)
Ue dx
y
d Zδ dδ
ve = − u dy + ue (3.206)
dx 0 dx
Podemos escribir:
d ∗ dδ
ve = (δ ue − δue ) + ue (3.207)
dx dx
Sea:
d due
ve = (ue δ ∗ ) − δ (3.208)
dx dx
Esta expresión da la velocidad transversal sobre la frontera de la capa limite en función de ue , δ ∗
y δ.
Figura 3.22: Capa limite sobre placa plana. Esquema del volumen de control.
d Z δ(x) d
VE (x) = u(x, y) dy = ue (δ − δ ∗ ) (3.215)
dx 0 dx
es decir:
d
VE (x) = ue (δ − δ ∗ ) (3.216)
dx
Se puede simplificar con:
d due
ve = (ue δ ∗ ) − δ (3.217)
dx dx
Es decir:
dδ
VE (x) = ue − ve (3.218)
dx
La velocidad de aspiración VE (x) no es igual a la velocidad transversal ve .
dθ 1 due τw
+ (δ ∗ + 2θ) = 2 (3.219)
dx ue dx ρue
3.11. RESOLUCIÓN APROXIMADA LA ECUACIÓN INTEGRAL POR EL MÉTODO DE KÁRMÁN-PO
Esta ecuación integral, desarrollado por Von Kármán, relaciona los espesores: de desplazamiento y
de deficit de cantidad de movimiento, a la variación de la velocidad exterior y al esfuerzo parietal.
La ecuación 3.219 puede también se escribir:
due d 2 τw
δ ∗ ue + uθ = (3.220)
dx dx e ρ
Se encuentra 3.220 a desarrollando 3.219. Además, a partir de la ecuación:
dθ 1 due τw Cf
+ (δ ∗ + 2θ) = 2 ≡ (3.221)
dx ue dx ρue 2
se puede escribir
dθ θ due Cf
+ (H + 2) = (3.222)
dx ue dx 2
Esta ecuación integral relaciona: los espesores: de desplazamiento y de déficit de cantidad de mo-
vimiento, a la variación de la velocidad exterior y al coeficiente de fricción local.
Comentario
La ecuación de Von Karman es aplicable a las capas limites incompresibles, sometidas o no
a un gradiente de presión. El flujo de la capa limite puede ser laminar o turbulento. Notamos
que estas expresiones generaliza la relación presentada en la parte 1 de esta clase sobre la base de
un balance de las cantidades de movimiento por un volumen de control elemental.
dθ 1 due τw
+ (δ ∗ + 2θ) = 2 (3.225)
dx ue dx ρue
para obtener de manera simple pero aproximativa la mayorı́a de las propiedades importantes del
dθ τw
flujo en la capa limite: δ, θ, dx , ρu2 .
e
Para simplificar, consideramos el caso de una capa limite laminar sobre una placa plana en
ausencia de gradiente de presión:
due
=0 (3.226)
dx
con el sistema de referencia presentado en Figura 3.23.
Figura 3.23: Capa limite sobre placa plana con su sistema de referencia.
∂u ∂u ∂ 2u
u +v =ν 2 (3.229)
∂x ∂y ∂y
Para una capa limite sobre una placa plana, en ausencia de gradiente de presión, la curvatura
del perfil de velocidad es nulo a la pared.
La idea de Kármán y Polhausen consiste en representar el perfil de velocidad en una forma
aproximada satisfaciendo a lo menos parcialmente a las condiciones de frontera del problema.
se puede entonces definir, por ejemplo, un polinomio de grado 3 para definir el perfil de velocidad
según u:
u = a0 + a1 y + a2 y 2 + a3 y 3 (3.233)
∂u
= a1 + 2a2 y + 3a3 y 2 (3.234)
∂y
∂ 2u
= 2a2 + 6a3 y (3.235)
∂y 2
u e = a1 δ + a3 δ 3
0 = a1 + 3a3 δ 2 (3.236)
se deduce fácilmente:
3 ue ue
a1 = a3 = − (3.237)
2 δ 2δ 3
El perfil de velocidad toma entonces la forma:
3
u 3y 1 y
= − (3.238)
ue 2δ 2 δ
122 CAPÍTULO 3. CAPA LIMITE LAMINAR DINÁMICA
u 3 1
= η − η3 (3.242)
ue 2 2
Después integración:
3 9 1 3 1
θ=δ − − + − (3.244)
4 12 8 10 28
El espesor de cantidad de movimiento se expresa:
39
θ=δ (3.245)
280
Ası́ el perfil de velocidad dado por la Ec. 3.242 entrega un esfuerzo de corte:
3 µue
τw = (3.247)
2 δ
En esta condición tenemos:
τw 3 ν
2
= (3.248)
ρue 2 δue
Para determinar θ y τw debemos determinar ahora el espesor de la capa limite δ.
3.12. UN EJEMPLO DE UTILIZACIÓN DEL MÉTODO DE KÁRMÁN-POLHAUSEN 123
Figura 3.24:
Divergente (Ver Figura 3.25): La capa limite en esta entrada de motor de avión se desarrolla
en presencia de un gradiente de presión positivo entonces desfavorable.
Las capas limites sobre las alabes de una turbo-maquina son sometidas a gradientes de
presión. También, sobre el estrados de una ala, el gradiente de presión es positivo a partir de
una distancia del borde de ataque del perfile. Este gradiente puede generar el despegamiento
de la capa limite (Ver Figura 3.26).
3.13. CAPA LIMITE LAMINAR INCOMPRESIBLE CON GRADIENTE DE PRESIÓN EXTERNO. ECU
3.13.2. Objetivos
El interés principal del análisis teórico que vamos realizar es para permitir entender porque un
gradiente de presión positivo, correspondiendo a una deceleración del flujo, es desfavorable, y
puede conducir al despegamiento de la capa limite.
Porque un gradiente de presión negativo, correspondiendo a una aceleración del flujo es al
contrario favorable.
El estudio de los efectos de los gradientes de presión externos es basado sobre las ecuaciones
de Falkner-Skan.
Estas ecuaciones son obtenidas a partir de la búsqueda de soluciones auto-similares de las
ecuaciones de la capa limite.
La solución existe cuando el campo de velocidad externo es de la forma:
u e = c0 x m (3.260)
1 dp 1 dp due
= −c0 mx2m−1 porque = −ue (3.261)
ρ dx ρ dx dx
126 CAPÍTULO 3. CAPA LIMITE LAMINAR DINÁMICA
Figura 3.27: Flujo sobre una placa plana que desarolla una capa limite con una pared exterior.
∂u ∂v
+ =0 (3.262)
∂x ∂y
∂u ∂u 1 dP ∂ 2u
u +v =− +ν 2 (3.263)
∂x ∂y ρ dx ∂y
u(x, 0) = 0 v(x, 0) = 0 u(x, ∞) = ue (x) (3.264)
∂u ∂u due ∂ 2u
u +v = ue +ν 2 (3.267)
∂x ∂y ∂x ∂y
Como para el problema de Blasius es cómodo introducir la función de corriente Ψ(x, y) tal que:
∂Ψ ∂Ψ
u= v=− (3.268)
∂y ∂x
3.14. ESTABLECIMIENTO DE LA ECUACIÓN DE FALKNER-SKAN 127
∂Ψ ∂ 2 Ψ ∂Ψ ∂ 2 Ψ due ∂ 3Ψ
− = u e + ν (3.269)
∂y ∂x∂y ∂x ∂y 2 ∂x ∂y 3
con las 3 condiciones a la frontera:
∂Ψ ∂Ψ ∂Ψ
(x, 0) = 0 (x, 0) = 0 (x, ∞) = ue (x) (3.270)
∂y ∂x ∂y
Aquı́, tenemos:
due
ue 6= 0 (3.271)
∂x
Además, la condición:
∂Ψ
(x, 0) = 0 (3.272)
∂x
puede ser transformada notando:
Z
∂Ψ
(x, 0) dx = Ψ(x, 0) − Ψ(0, 0) = 0 (3.273)
∂x
Si elegimos:
Ψ(0, 0) = 0 ⇒ Ψ(x, 0) = 0 (3.274)
La elección de la condición Ψ(0, 0) = 0 no modifica el grado de generalidad pero sabemos que
Ψ es siempre definida como una constante.
Ası́ vamos a buscar una solución al problema:
∂Ψ ∂ 2 Ψ ∂Ψ ∂ 2 Ψ due ∂ 3Ψ
− = u e + ν (3.275)
∂y ∂x∂y ∂x ∂y 2 ∂x ∂y 3
con las 3 condiciones a la frontera:
∂Ψ ∂Ψ
(x, 0) = 0 Ψ(x, 0) = 0 (x, ∞) = ue (x) (3.276)
∂y ∂y
h0 ∂
! !
∂ ∂ ∂ζ ∂ ∂η ∂
= + = +η (3.279)
∂x ∂ζ ∂x ∂η ∂x ∂ζ h ∂η
! !
∂ ∂ ∂ζ ∂ ∂η ∂
= + =h (3.280)
∂y ∂ζ ∂y ∂η ∂y ∂η
Por ejemplo, tenemos:
h0 ∂Ψ
! " #
∂ ∂Ψ ∂ ∂Ψ
=h +η (3.281)
∂y ∂x ∂η ∂ζ h ∂η
h0
y como h
depende solamente de x, es decir de ζ se puede escribir:
!
∂ ∂Ψ ∂ 2Ψ ∂Ψ ∂ 2Ψ
=h + h0 + ηh0 2 (3.282)
∂y ∂x ∂η∂ζ ∂η ∂η
De la misma manera, podemos encontrar:
∂ 2Ψ 2
2∂ Ψ
= h (3.283)
∂y 2 ∂η 2
∂ 3Ψ 3
3∂ Ψ
= h (3.284)
∂y 3 ∂η 3
Ası́, la ecuación
∂Ψ ∂ 2 Ψ ∂Ψ ∂ 2 Ψ due ∂ 3Ψ
− = ue + ν (3.285)
∂y ∂x∂y ∂x ∂y 2 ∂x ∂y 3
puede ser escrito:
h0 ∂Ψ
" # " #" #
∂ ∂ 2Ψ 0 ∂Ψ
2
0∂ Ψ ∂Ψ 2
2∂ Ψ
h h +h + ηh − +η h
∂η ∂η∂ζ ∂η ∂η 2 ∂ζ h ∂η ∂η 2
due ∂ 3Ψ
= ue + νh3 3 (3.286)
∂x ∂η
Finalmente, la ecuación:
∂Ψ ∂ 2 Ψ ∂Ψ ∂ 2 Ψ due ∂ 3Ψ
− = ue + ν (3.287)
∂y ∂x∂y ∂x ∂y 2 ∂x ∂y 3
se escribe en variables η y ζ:
" #
∂ ∂ 2Ψ ∂Ψ ∂Ψ ∂ 2 Ψ due 3
3∂ Ψ
h h + h0 − h2 = ue + νh (3.288)
∂η ∂η∂ζ ∂η ∂ζ ∂η 2 ∂x ∂η 3
Sea:
2
(α + β) c21 c22 ζ 2α+2β−1 (f 0 ) − c21 c22 βζ 2α+2β−1 f f 00
= mc20 ζ 2m−1 + νc31 c2 ζ 3α+β f 000 (3.294)
La ecuación
2
(α + β) c21 c22 ζ 2α+2β−1 (f 0 ) − c21 c22 βζ 2α+2β−1 f f 00
= mc20 ζ 2m−1 + νc31 c2 ζ 3α+β f 000 (3.295)
debe estar homogénea en ζ (longitud):
2α + 2β − 1 = 2m − 1
2α + 2β − 1 = 3α + β (3.296)
Sea:
2α + 2β − 1 = 2m − 1
α−β−1 = 0 (3.297)
La relación 3.297 conduce a elegir los siguientes parámetros:
m−1 m+1
α= β= (3.298)
2 2
Si reemplazamos en la ecuación:
2
(α + β) c21 c22 ζ 2α+2β−1 (f 0 ) − c21 c22 βζ 2α+2β−1 f f 00
= mc20 ζ 2m−1 + νc31 c2 ζ 3α+β f 000 (3.299)
m−1 m+1
α= β= (3.300)
2 2
130 CAPÍTULO 3. CAPA LIMITE LAMINAR DINÁMICA
2 m − 1 2 2 00
mc21 c22 (f 0 ) − c1 c2 f f = mc20 + νc31 c2 f 000 (3.301)
2
2 m − 1 2 2 00
mc21 c22 (f 0 ) − c1 c2 f f = mc20 + νc31 c2 f 000 (3.302)
2
En general, c0 es dado (se conoce ue ) pero c1 y c2 son constantes que podemos especificar
arbitrariamente.
Es cómodo imponer:
c2 = νc1 (3.304)
Obtenemos:
c0 √
r
c1 = c2 = νc0 (3.305)
ν
Las ecuaciones donde las variables c1 y c2 son expresadas en función de c0 permiten simplificar
la ecuación:
2 m − 1 2 2 00
mc21 c22 (f 0 ) − c1 c2 f f = mc20 + νc31 c2 f 000 (3.306)
2
m + 1 00
f 000 + f f + m(1 − f 02 ) = 0 (3.307)
2
ue = c0 xm (3.310)
Ψ = c2 ζ β f (η) (3.312)
√ m+1
sea como c2 = νc0 y β = 2
:
√ √
Ψ= νc0 x(m+1)/2 f (η) = νc0 xm xf (η) (3.313)
Ası́:
√
Ψ= νue xf (η) (3.314)
Blasius (sin gradiente de presión) tuvimos:
Ψ(x, y) = (ue xν)1/2 f (η) (3.315)
f 0 (∞) = 1 (3.322)
3.14.3.4. Importante
Si m = 0 tenemos Blasius.
Esfuerzo parietal
τw f 00 (0)
Cf = 1 = 2 √ (3.328)
2
ρu2e Rex
Velocidad transversal en función de f
∂Ψ ∂ h i
v(x, y) = − =− (ue xν)1/2 f (η) (3.329)
∂x ∂x
es decir:
ue η 0 m+1
v(x, y) = √ f − f (3.330)
Rex 2 2
f 000 + f f 00 − f 02 + 1 = 0 (3.339)
134 CAPÍTULO 3. CAPA LIMITE LAMINAR DINÁMICA
τw f 00 (0)
Cf = 1 = 2 √ (3.340)
2
ρu2e Rex
√
Rex 00
τw = µue f (0) (3.341)
x
3.15.2.2. Caso m = 0
Presentamos otro perfile de velocidad reducida u
ue
= f 0 (η) sobre la Figura 3.33.
3.15.2.3. Caso m = 1
Presentamos otro perfile de velocidad reducida u
ue
= f 0 (η) sobre la Figura 3.34.
axial (m > 0), aparece como favorable y la pendiente de la velocidad a la pared aumenta cuando
m crece.
Figura 3.36: Comportamiento de una capa limite sobre una placa plana en una sección convergente.
Figura 3.37: Comportamiento de una capa limite sobre una placa plana en una sección divergente.
Capı́tulo 4
Figura 4.1:
Existe lo que se conoce como una capa externa (80 % del espesor total δ), donde el gradiente
∂U
∂y
es pequeño y donde el perfil de velocidad U (y) está poco influenciando por la pared.
Mientras que también se puede diferenciar una capa interna muy cerca de la pared (20 %
del espesor total δ) donde el flujo es controlado por las condiciones a la pared (rugosidad) y
por la naturaleza del fluido. En esta zona, la ley que describe el perfil de velocidad U (y) se
conoce como ley de pared.
En la capa interna podemos distinguir dos 2 sub-capas (ver gráfico en log de y distancia a la
pared):
139
140 CAPÍTULO 4. INTRODUCCIÓN A LA CAPA LIMITE TURBULENTA DINÁMICA
Figura 4.2:
la zona intermediaria 5 < y + < 70 es una zona de transición del flujo entre un régimen
dominado por la viscosidad molecular y un régimen dominado por los esfuerzos turbulentos.
• No se dispone de formulas analı́ticas de la ecuación de NS para describir el flujo en esta
zona.
la zona logaritmica (el flujo es turbulento) donde
1
u+ = ln y + + C cuando 70 > y + (4.4)
κ
donde κ = 0, 4 es la constante de von Kárman y 3, 7 < C < 5
4.1. PERFIL UNIVERSAL DE VELOCIDAD EN LA CAPA LIMITE 141
Figura 4.3:
• Winckelmans and Bricteux (2008) dan κ = 0,37 y C = 3,7 a partir de una DNS.
Comentario:
Esta ultima relación da una ecuación que relaciona el coeficiente de perdida de carga f , la
velocidad de fricción y la velocidad promedio:
1 Z 2e
U= ū(y) dy (4.6)
2e 0
s
f
uτ = U (4.7)
8
lm = κy si 0 ≤ y ≤ yl (4.8)
lm = κyl si y ≥ yl (4.9)
Sin embargo,
∂ p̄ −ρu2τ
= (4.11)
∂x h
Comentario:
Se puede extender a cualquier sección que no tenga una singularidad (sección elı́ptica por
ejemplo).
subcapa viscosa
0 < y+ < 8 : u+ = y + (4.14)
dū uτ
µt = ρκuτ = (4.16)
dy κy
Las constantes que aparecen en las formulas son determinadas experimentalmente o con
simulaciones numéricas directas de la turbulencia (DNS)
La zona intermediaria 8 < y + < 30 es una zona de transición del flujo entre un régimen
dominado por la viscosidad molecular y un régimen dominado por los esfuerzos turbulentos.
Figura 4.4:
εuτ
Reτ = (4.18)
ν
Régimen liso para Reτ < 3
1 y+
u+ = ln + 8, 5 (4.19)
κ ε
Régimen de transición 5 < Reτ < 70
144 CAPÍTULO 4. INTRODUCCIÓN A LA CAPA LIMITE TURBULENTA DINÁMICA
Capı́tulo 5
Flujo potencial
5.1. Generalidades
5.1.1. Introducción
Primero vamos a definir que es un flujo potencial. Sabemos que el movimiento de los fluidos
ideales (no viscosos) es descrito por la ecuación de Euler:
−→
rotV~ (~r, t) = ∇ × V~ (~r, t) = 0 (5.2)
5.1.1.1. Definición
Llamamos flujo potencial, un flujo que posee la siguiente propiedad:
−→
~ = rotV~ (~r, t) = 0 ∀~r, t
ω (5.3)
5.1.1.2. Definición
La condición de flujo irrotacional constituye una condición necesaria y suficiente
para tener una función potencial cuya deriva el campo de velocidad.
Se nota φ a esta función llamada potencial de velocidad tal que, por definición:
−−→
V~ (~r, t) = gradφ(~r, t) (5.4)
5.1.1.3. Comentarios
Evidentemente un fluido ideal no existe en la naturaleza. Sin embargo, existe condicio-
nes para que un flujo tenga (aproximadamente) esta propiedad. Por ejemplo: los flujos lejos de
obstáculos o frontera.
145
146 CAPÍTULO 5. FLUJO POTENCIAL
5.1.1.4. Expresión de V~
Sistema de coordenadas ~ =−
V
−→
gradφ
∂φ ∂φ ∂φ
Cartesianas vx = ∂x , vy = ∂y , vz = ∂z
∂φ 1 ∂φ ∂φ
Cilı́ndricas vr = ∂r , vθ = r ∂θ , vz = ∂z
∂φ 1 ∂φ 1 ∂φ
Esféricas vr = ∂r , vθ = r ∂θ , vϕ = r sin θ ∂ϕ
5.1.1.5. Propiedades
Supongamos que el flujo es incompresible, se puede escribir:
h−−→ i
divV~ = div gradφ(~r, t) = ∆φ(~r, t) = 0 (5.5)
Se obtiene una ecuación
∆φ(~r, t) = 0 (5.6)
llamada: Ecuación de Laplace
Si, tenemos una interfaz S entre dos fluidos ideales 1 y 2, la condición a la frontera es V~n1 = V~n2 .
Esta condición se escribe:
" # " #
∂φ1 ∂φ2
= (5.8)
∂n S
∂n S
5.1.1.7. En resumen
Para determinar la función φ potencial de velocidades se necesita:
Como el problema planteado es totalmente cinemático por que se debe determinar el potencial
Φ. Si el flujo es cinematicamente admisible este flujo es admisible dinámicamente. El problema
dinámico se traduce usando la ecuación de Bernoulli (fluido ideal!!!) para un flujo irrotacional:
∂φ 1 Z
dp
+ |∇φ|2 + Φ + = C(t) (5.9)
∂t 2 ρ(p)
∆φ3 = 0 (5.11)
5.2.1.1. Principio
Si φ1 y φ2 definidas sobre el mismo dominio, representan los potenciales de velocidad de 2 flujos
irrotacionales de un fluido incompresibles, entonces toda combinación lineal de estas funciones son
también solución de un otro flujo.
Figura 5.1: Vector velocidad, funciones potencial y corriente para dos flujos potenciales 1 y 2.
Estos 2 flujos de base pueden generan un flujo 3.
1 y −−→ 2
E = ρ gradφ3 dV (5.14)
2
V
1 y −−→ 2
E = ρ gradφ3 dV (5.17)
2
V
1 y −−→
E = ρ div φ3 gradφ3 dV (5.18)
2
V
1 x ∂φ3
E = ρ φ3 dS (5.19)
2 ∂n
S
∂φ3
Pero sabemos que ∂n
= 0 sobre la superficie S y por consecuencia:
5.2. PROPIEDADES GENERALES DE LA ECUACIÓN DE LAPLACE 149
E =0 (5.20)
−−→ 2 −−→ 2
Como E es la integral de gradφ3 que no es negativa. Tenemos entonces que gradφ3 no
puede anularse si no es igual a zero en todo V .
Finalmente, para que E = 0 se necesita que:
−−→
gradφ3 = 0 en todo S ⇒ V3 = 0 ⇒ V2 − V1 = 0 (5.21)
5.2.2.2. Comentarios
−−→
Sin producir cambio en el campo de velocidad V~ = gradφ, φ1 y φ2 pueden ser diferentes de
una constante (C1 o C2 ). Es decir, podemos tener:
φ1 + C1 φ2 + C2 (5.22)
por que:
−−→
V~ = gradφ (5.26)
y
divV~ = divV~ 0 = 0 (5.27)
150 CAPÍTULO 5. FLUJO POTENCIAL
Finalmente, usando el teorema de Ostrogradsky, la ultima integral puede escribirse bajo la forma
integral de superficie sobre las fronteras:
y h i x
ρ div φ(V~ 0 − V~ ) dV = ρ φ(V~ 0 − V~ ) dS = 0 (5.28)
V S
por que V~ y V~ 0 son las mismas componentes normales a S1 y S2 . Ası́, la diferencia entre las dos
energı́as se escribe:
1 y ~0
E0 −E = ρ (v − ~v )2 dV (5.29)
2
V
Esta relación es positiva como v~0 es diferente de ~v . Entonces E 0 es superior a E , y por consecuencia
la energı́a cinética del flujo potencial es la menor entre todas las otras soluciones admisibles por
el flujo al interior V .
Vamos después analizar como los campos de velocidad de estos flujos simples pueden ser
combinados para resolver problemas mas complejos.
5.3. FLUJOS ELEMENTALES BIDIMENSIONALES PLANOS IRROTACIONALES ESTACIONARIOS1
∂ 2Ψ ∂ 2Ψ
∆Ψ = + =0 (5.32)
∂x2 ∂y 2
Prueba
∂ 2Ψ ∂ 2Ψ ∂ ∂
∆Ψ = 2
+ 2
= (−v) + (u) = 0 (5.33)
∂x ∂y ∂x ∂y
−→
porqué el flujo es irrotacional: rotV~ (~r, t) = 0
∂Ψ
Cartesianas vx = ∂y
, vy = − ∂Ψ
∂x
1 ∂Ψ
Cilı́ndricas vr = r ∂θ
, vθ = − ∂Ψ
∂r
Invariante en ϕ y uϕ = 0
1 ∂Ψ 1 ∂Ψ
Coordenadas esféricas vr = r2 sin θ ∂θ
, vθ = − r sin θ ∂r
eso implica que existe una función φ(x, y), potencial de velocidad tal que:
∂φ ∂φ
u= v= (5.35)
∂x ∂y
Llamamos las lineas φ =Cte las equipotenciales del campo de velocidad (u, v).
∂φ
= ∂Ψ
(
u = ∂x ∂y
~v = ∂φ (5.37)
v = ∂y
= − ∂Ψ
∂x
Sabemos que:
−−→ −−→
~v = gradφ = gradΨ × e~z (5.39)
según las coordenadas (x, y, z)
∂φ ∂Ψ ∂φ ∂Ψ
= = vx = U y =− = vy = 0 (5.40)
∂x ∂y ∂y ∂x
entonces por integración tenemos:
φ = Ux Ψ = Uy (5.41)
Figura 5.4: Flujo uniforme, las lineas equipotenciales bidimensional son derechas perpendicular (o
normal) a las lineas de corriente
∂φ 1 ∂Ψ ∂φ 1 ∂Ψ
=− = vz = U y = = vr = 0 (5.45)
∂z r ∂r ∂r r ∂z
es decir después integración:
r2
φ = Uz Ψ = −U (5.46)
2
1
φ = U r cos θ Ψ = U r2 sin2 θ (5.50)
2
154 CAPÍTULO 5. FLUJO POTENCIAL
5.3.5.4. Resumen
Flujo Potencial φ de Función de corriente Ψ
velocidad
Flujo uniforme según Ox 2D φ = Ux Ψ = Uy
cartesianas
Flujo uniforme según Ox φ = U r cos θ Ψ = U r sin θ
cilı́ndricas
Flujo uniforme según Oz φ = Uz Ψ = − 12 U r2
cilı́ndricas
Flujo uniforme según Oz φ = U r cos θ Ψ = 12 U r2 sin2 θ
esféricas
k = 0 si x > 0
k = 1 si x < 0
φ = κθ y Ψ = −κ ln r (5.62)
Γ
con κ = 2π
.
a partir de un punto
• Q > 0 por la fuente
hacia un punto
• Q < 0 por el sumidero
En coordenadas cilı́ndricas (r, θ, z), como el flujo es plano, la única componente no nula del
campo de velocidad es:
Q
vr (r) = y vθ = 0 vz = 0 (5.63)
2πr
Se observa que la velocidad no es definida en r = 0.
Si calculamos el flujo de velocidad a través un circulo de radio r centrado sobre el origen, se
obtiene: I Z 2π
~
~v .dl = rvr dθ = Q (5.64)
C 0
Q es la intensidad de la fuente o del sumidero (m2 .s−1 ).
A partir del expresión del vector velocidad, se puede determinar las funciones φ y Ψ:
−−→ −−→
~v = gradφ = gradΨ × e~z (5.65)
según las coordenadas (r, θ, z).
∂φ 1 ∂Ψ Q 1 ∂φ ∂Ψ
= = vr = y =− = vθ = 0 (5.66)
∂r r ∂θ 2πr r ∂θ ∂r
Damos la expresión del potencial de velocidad y de la función de corriente:
Q r Qθ
φ= ln y Ψ= (5.67)
2π r0 2π
Como en el caso precedente, φ y Ψ son definidas con una constante arbitraria.
Q r Qθ
φ= ln y Ψ= (5.68)
2π r0 2π
las lineas de corriente son rectas radiales: θ = Cte
las curvas equipotenciales son circunferencias: r = Cte
Figura 5.6:
5.3.7.3. Comentarios
Notamos la correspondencia entre la solución de este problema y la solución de un vórtice.
Las funciones φ y Ψ son intercambiadas. Las lineas de corriente radiales, en el caso de la fuente,
son idénticas a las equipotenciales del flujo en un vórtice. Las lineas de corrientes del vórtice son
158 CAPÍTULO 5. FLUJO POTENCIAL
cı́rculos centrados alrededor del origen, y corresponde a las equipotenciales del flujo creado por la
fuente.
Esta correspondencia será analizado de manera mas directa cuando vamos a tratar la noción
de potencial compleja de la velocidad.
Calculemos el campo de velocidad del flujo dipolar en el caso plano (dipolo lineal). Para
simplificar tomamos el origen de coordenadas en S1 . Las dos funciones potencial de velocidad φ1
y φ2 se escriben:
~
Q |~r − d|
φ1 = ln (5.74)
2π r0
y !
Q |~r|
φ2 = − ln (5.75)
2π r0
5.3. FLUJOS ELEMENTALES BIDIMENSIONALES PLANOS IRROTACIONALES ESTACIONARIOS1
~
!
Q |~r − d| − ln
|~
r |
φ = φ1 + φ2 = ln
2π r0 r0
Q |~r − d| ~
= ln (5.76)
2π |~r|
Q~r.d~
φ=− (5.79)
2πr2
y en el limite d~ → 0 resulta:
~r.~p
φ=− (5.80)
2πr2
Usando coordenadas polares en el plano del flujo y tomando p~ como origen de θ tenemos:
p cos θ
φ=− (5.81)
2π r
Para obtener las componentes del vector velocidad en coordenadas polares, sabemos que:
!
−−→ −−→ p cos θ
~v = gradφ = grad − (5.82)
2π r
es decir
∂φ p cos θ
vr = = (5.83)
∂r 2πr2
y
1 ∂φ p sin θ
vθ = = (5.84)
r ∂θ 2πr2
A partir de las componentes de la velocidad, se obtiene la función de corriente Ψ después una
integración:
1 ∂Ψ p cos θ
vr = = (5.85)
r ∂θ 2πr2
∂Ψ p sin θ
vθ = − = (5.86)
∂r 2πr2
Sea:
p sin θ p~ × ~r
Ψ= = (5.87)
2π r 2πr2
160 CAPÍTULO 5. FLUJO POTENCIAL
Sobre la Figura 5.9, las lineas de corriente y equipotenciales son representadas. Las lineas de
corrientes son una familia de circunferencias tangentes a p~ y que pasan por el origen. Las equipo-
tenciales son una familia de circunferencias ortogonales a las primeras, que pasan por el origen y
son tangentes a la recta θ = ± π2 . El campo de velocidad es representado sobre la Figura 5.10.
5.3.11. Resumen
5.3.11.1. Flujo uniforme
Flujo Potencial φ de Función de co-
velocidad rriente Ψ
Flujo uniforme según Ox 2D φ = Ux Ψ = Uy
cartesianas
Flujo uniforme según Ox φ = U r cos θ Ψ = U r sin θ
cilı́ndricas
Flujo uniforme según Oz φ = Uz Ψ = − 12 U r2
cilı́ndricas
Flujo uniforme según Oz φ = U r cos θ Ψ = 12 U r2 sin2 θ
esféricas
5.3.11.3. Dipolo
Flujo Potencial φ de Función de corriente
velocidad Ψ
Dipolo p~ = Qd~ p
~.~
φ = − 2πrr
2
p sin θ
Ψ = 2π r
cilı́ndricas
p sin2 θ
Dipolo p~ = Qd~ p
~.~
φ = − 4πrr
3 Ψ= 4π r
esféricas
5.3.12. Ejercicios
5.3.12.1. Ejercicios 1
Se considera, en un sistema de referencia cartesiano ortonormado, un campo de vector V~ =
(U, V, W ) tal que:
U = kx V = ky W = −2kz (5.90)
3. Determinar el potencial de velocidad de este flujo. Podemos asociar a este flujo una función
de corriente? Por qué?
162 CAPÍTULO 5. FLUJO POTENCIAL
4. Reescribir el campo de velocidad en coordenadas cilı́ndricas de eje Oz. Mostrar que este flujo
es axisimetrico y describir con un esquema este flujo.
5.3.12.2. Ejercicio 2
Se considera el flujo generado por la superposición:
un flujo uniforme a la velocidad V0 alineado con el eje x;
el efecto Magnus.
p cos θ
φdipolo = − (5.92)
2πr
p cos θ
φ = φunif orme + φdipolo = U r cos θ −
2πr
p
= Ur − cos θ (5.93)
2πr
Figura 5.12: Esquema modelando un flujo alrededor de un cilindro circular sin circulación: el flujo
uniforme + dipolo .
164 CAPÍTULO 5. FLUJO POTENCIAL
∂φ p
vr = = U+ cos θ (5.94)
∂r 2πr2
1 ∂φ p
vθ = =− U− sin θ (5.95)
r ∂θ 2πr2
~
~v = U para r → ∞ ∀ θ (5.96)
vr (r = R) = 0 ∀ θ (5.97)
Comentario: La segunda condición impone que la componente de velocidad normal es nula sobre
el cilindro (fluido ideal).
La primera condición
~
~v = U para r → ∞ (5.98)
es verificada, porque el campo de velocidad creado por el dipolo:
p
vr = U + cos θ (5.99)
2πr2
p
vθ = − U − sin θ (5.100)
2πr2
1
decrece en r2
La segunda condición da
vr (r = R) = 0 ∀ θ (5.101)
Es decir por:
p p
vr = U + 2
cos θ o vθ = − U − sin θ (5.102)
2πr 2πr2
p
= −U ⇒ p = −U 2πR2 (5.103)
2πR2
Finalmente se determinó el valor de p y tenemos la función potencial de velocidad que se
escribe:
!
p R2
φ = Ur − cos θ = U r cos θ 1 + 2 (5.104)
2πr r
5.5. FLUJO ALREDEDOR DE UN CILINDRO CIRCULAR SIN CIRCULACIÓN: PARADOJA DE ALEM
∂Ψ
vθ = − (5.107)
∂r
después de la integración tenemos:
!
R2
Ψ = U sin θ r − (5.108)
r
En caso general, las funciones de corrientes no son curvas muy simples.
Podemos notar que para Ψ = 0 tenemos dos semirrectas que parten de dos puntos de estanca-
miento situados en r = R con θ = 0 o θ = π.
A1 = (−a = R, 0) A2 = (a = R, 0) (5.109)
flujo de fluido ideal, las únicas fuerzas de superficie son las fuerzas de presión. La resultante de la
fuerza sobre el cilindro se obtiene por integración sobre el contorno:
x
F~ = − P ~ndS (5.112)
S
A partir del campo de velocidad podemos calcular la presión en cada punto. Por eso se usa la
ecuación de Bernoulli, considerando que tenemos ninguna fuerza de volumen.
Escribimos el teorema de Bernoulli, entre un punto de estancamiento (PA1 , VA1 = 0) y un punto
corriente donde la presión y la velocidad son (P, V ), tenemos:
1 1
P + ρV~ 2 = PA1 + ρV~A1 = Cte = PA1 (5.114)
2 2
Sea:
1
P = PA1 − ρV~ 2 (5.115)
2
En r = R, las componentes de la velocidad se escriben:
Ası́, tenemos:
1
P = PA1 − ρV~ 2 = PA1 − 2ρU 2 sin2 θ (5.117)
2
La fuerza ejercida sobre el cilindro se puede escribir entonces:
Z 2π Z 2π
F~ = − PA1~nRdθ + 2ρU 2 sin2 θ~nRdθ (5.118)
0 0
Figura 5.14:
Figura 5.15: Esquema modelando un flujo alrededor de un cilindro circular con circulación: el flujo
uniforme + dipolo + vortice.
Entonces se deduce:
!
∂φ R2
vr = = U 1 − 2 cos θ (5.123)
∂r r
y
!
1 ∂φ R2 Γ
vθ = = −U 1 + 2 sin θ + (5.124)
r ∂θ r 2πr
!
R2 Γ r
Ψ = U sin θ r − − ln (5.126)
r 2π R
Figura 5.16: Circulación débil: 0 < |Γ| < 4πR|U |. Dos sentidos de circulación son presentados.
Figura 5.17: Lı́neas de corriente y las equipotenciales de velocidad para un cilindro en un flujo
uniforme, con una circulación Γ no nula: UΓR = −6
Figura 5.18: Isobaras correspondientes al flujo presentado en la figura 5.17. Se han representado
∆p
las curvas de igual p = ρU 2
2
Tenemos trayectorias cerradas que encierran el cilindro (fluido detenido) y trayectorias abiertas
a mayor distancia. Según 5.124 podemos escribir:
!
R2 |Γ|
−|U | 1 + 2 + =0 (5.130)
r 2πr
Las dos condiciones 5.129 y 5.130 permiten tener respectivamente: vr = 0 y vθ = 0.
172 CAPÍTULO 5. FLUJO POTENCIAL
Figura 5.20: Circulación fuerte: |Γ| > 4πR|U |. Dos sentidos de circulación son presentados.
Para determinar las componentes de esta fuerza, se calcula la resultante de las fuerzas de
presión sobre el cilindro a partir del campo de presión p(r = R, θ). La presión verifica la
relación de Bernoulli 2 Tomando como referencia un punto al infinito (presión p∞ y velocidad U ),
se escribe:
1 1
p(r = R, θ) + ρvθ2 (r = R, θ) = p∞ + ρU 2 (5.134)
2 2
donde la presión sobre el cilindro:
!2
1 Γ
p = p∞ + ρU 2 1 − −2 sin θ + (5.135)
2 2πRU
La sustentación por unidad de longitud del cilindro (dirección z) es dada por la componente Fs :
−
→ 2π Z 2π Z 2π
Z
Fs = − p~nRdθ = − p cos θRdθ ~ex + − p sin θRdθ ~ey (5.136)
0 0 0
Usando la relación 5.135 es fácil de ver que el único término no nulo es:
" # " !#2π
−
→ Z 2π ρU Γ 2 ρU Γ 1 sin(2θ)
Fs = − sin θdθ ~ey = − θ− ~ey
0 π π 2 2 0
ρU Γ
= − π~ey
π
= −ρU Γ~ey (5.137)
Ası́,
−
→
Fs = −ρU Γ~ey + 0e~x (5.138)
−
→
La fuerza Fs es dirigida hacia arriba. Esta expresión se llama sustentación (también llamado
fuerza de Magnus).
El Alcyone es un buque oceanográfico que utiliza el efecto Magnus (o rotor Flettner) operado
por la sociedad Cousteau (Ver Figura 5.21).
2
se puede aplicar la relación de Bernoulli a todo el espacio por que el flujo de fluido es potencial.
174 CAPÍTULO 5. FLUJO POTENCIAL
Figura 5.22: Esquema representando la fuerza de F~ ejercitada sobre el obstáculo, cuando el obs-
~.
táculo es sometido a un flujo de fluido con una velocidad uniforme U
5.7. TEOREMA DE JOUKOVSKI 175
Sin embargo, se puede también encontrar el signo de la fuerza, aplicando la relación de Bernoulli
entre dos puntos opuestos situados arriba y a abajo del obstáculo. Con el sentido de circulación
Figura 5.23:
5.7.1.2. Comentarios
La fuerza de Magnus es responsable de la sustentación:
de los aviones: la geometrı́a del perfil alar (curvatura) es lo que origina la aparición de la
circulación necesaria.
pero esta fuerza permite también la propulsión o de sustentación de objetos por un hélice
giratorio: barco, aviones, helicópteros
finalmente, esta fuerza permite explicar la trayectoria de una pelota de tenis o de tenis de
mesa cuando se le ha dado un efecto.
V1
1. Topspin: Mano hacia arriba al momento del impacto (Ω1 ≈ d
), alta velocidad inicial
pero la trayectoria sera hacia abajo.
2. No hay influencia sobre la curvatura de la trayectoria
3. Slice el contrario del topspin (mano abajo): una velocidad tangencial V3 que crea de la
sustentación, la trayectoria sera plana.
176 CAPÍTULO 5. FLUJO POTENCIAL
Figura 5.24:
En coordenadas cartesianas:
y
Ψ(x, y) = U y + m arctan (5.141)
x
En coordenadas cilı́ndricas:
Ψ(r, θ) = U r sin θ + mθ (5.142)
Q
donde m = 2πb
y a su vez b es el espesor según z
Figura 5.26: Lineas de corriente generadas por una fuente en un flujo uniforme.
Este flujo tiene una simetrı́a de revolución alrededor de la dirección del flujo no perturbado.
Las lineas de estancamiento permiten definir una forma imaginaria de semi-cuerpo, un cuerpo
semi-infinito solido. Podemos determinar algunos puntos particular:
θ r
π πa
2 2
3π πa
2 2
θ → 0 con r → ∞ r sin θ → ±πa
178 CAPÍTULO 5. FLUJO POTENCIAL
y y
Ψ = U∞ y + m kπ + arctan − arctan (5.146)
x+a x−a
Ası́: !
2ay
Ψ = U∞ y − m arctan 2 (5.148)
x + y 2 − a2
5.9. EL OVALO DE RANKINE 179
" #
∂Ψ a (x2 − y 2 − a2 )
U= = U∞ − 2m (5.149)
∂y (x2 + y 2 + a2 )2 − 4a2 x2
" #
∂Ψ 2axy
V =− = 2m (5.150)
∂x (x2 + y 2 + a2 )2 − 4a2 x2
V = 0 es obtenido para x = 0 o y = 0
Si x = 0 la condición U = 0 da:
Q a
0 = V0 + × 2 (5.151)
π y + a2
El segundo miembro es siempre positivo (no nulo). Si y = 0 la condición U = 0 da:
Q a(x2 − a2 ) Q a(x2 − a2 )
0 = U∞ − × 2 = U∞ − × 2 (5.152)
π (x + a2 )2 − 4a2 x2 π (x − a2 )2
Figura 5.28:
m
Figura 5.29: Forma del ovalo de Rankine en función del parámetro Ua
~
Además, si el flujo es incompresible, sabemos que existe una función potencial vector Ψ(M, t) tal
que en todo punto del campo y cada instante, el vector velocidad verifica:
−
→~
∀M, ∀t V~ (M, t) = rotΨ(M, t) (5.158)
~
En geometrı́a bidimensional plana, vimos que la función vectorial Ψ(M, t) se reduce a una sola
componente Ψ(M, t) llamada función corriente.
La expresión del campo de velocidad se deduce por:
−−→
~v = gradΨ × e~z (5.159)
donde e~z es el vector unitario de la dirección ortogonal al plano del flujo:
En componentes, las relaciones 5.157 y 5.159 se escriben:
∂φ
= ∂Ψ
(
U = ∂x ∂y
~v = ∂φ (5.160)
V = ∂y
= − ∂Ψ
∂x
5.10.4.2. Demostración
Evaluamos la derivada por un incremento dz real de la variable compleja. dz = dx (el afijo z
se desplaza al largo del eje real), encontramos:
∂ ∂ (φ + iΨ) ∂φ ∂Ψ
f (z) = = +i = U − iV (5.166)
∂z ∂x ∂x ∂x
Un desplazamiento a lo largo del eje imaginario (dz = idy) da el mismo resultado:
∂ ∂ (φ + iΨ) ∂φ ∂Ψ
f (z) = = −i + = −iV + U (5.167)
∂z i∂y ∂y ∂y
De manera general, para cualquiera dirección del incremento del afijo de z, siempre tenemos:
∂
f (z) = U − iV (5.168)
∂z
5.10.4.3. Definición
∂
Llamamos por w(z, t) a la función velocidad compleja la derivada ∂z
f (z) y escribimos:
∂
w(z, t) ≡ f (z) = U (x, y, t) − iV (x, y, t) (5.169)
∂z
El modulo de la velocidad en todo punto es:
√ √
|w(z)| = U 2 + V 2 = ww∗ (5.170)
= ~ +i
~v .dl ~v .~ndl
C C
= Γ + iQ (5.174)
~
donde ~n representa el vector unitario normal al elemento de contorno dl.
El segundo termino representa el caudal del fluido (medido por unidad de longitud perpen-
dicular al plano del flujo) a partir de las fuentes situadas al interior de dC .
f (z) = V∞ z (5.176)
184 CAPÍTULO 5. FLUJO POTENCIAL
Q
f (z) = log (z − z0 ) (5.177)
2π
donde Q es una constante real y z0 un parámetro complejo, representa el flujo:
o un sumidero Q < 0
posicionado al punto z = z0 .
Tenemos el siguiente corte en el eje Ox, en razón del carácter multiforme del logaritmo
complejo
5.11.3. Vórtice
Bajo la restricción arg(z) ∈ ]0, 2π[, la función potencial compleja3
Γ
f (z) = −i log (z − z0 ) (5.178)
2π
donde Γ es una constante real y z0 un parámetro complejo, representa el flujo descrito por un
vórtice (Ver Figura 5.33):
orientado en el sentido trigonométrico Γ > 0
5.11.4. Doblete
Se llama doblete o dipolo un conjunto fuente-sumidero del mismo caudal |Q| (Ver Figura ??):
3
Cuidado: el sı́mbolo log es el logaritmo del número complejo, sea, con z = reiθ , log(z) = ln r + iθ.
186 CAPÍTULO 5. FLUJO POTENCIAL
5.11.5. Ejercicio
5.11.5.1. Flujo cerca de un punto de estancamiento
El flujo cerca de un punto de estancamiento es dada por la función compleja:
a
f (z) = z 2 (5.181)
2
1 1
donde a es una constante. El valor de a es a
= 5
s (segundos).
x = 3 cm y y = 0, 2 cm
Para el valor limite n = 12 , tenemos el contorno de un medio plano. Un ultimo caso importante
se obtiene para n = 2. Tenemos un flujo normal a un plano, donde tenemos localmente un punto
de estancamiento (Ver Figura 5.35).
Figura 5.35: (a) Flujo alrededor de una placa. (b) Flujo impactando sobre una placa.
Figura 5.36:
5.13.1.1. Definición
La representación del dominio d en el dominio D es dicho conforme si y solamente si la función
H(z) es: definida y continua en d, biyectiva, holormorfa y con derivada no nula en todo d.
Las hipótesis: definida y continua en d, biyectiva, holormorfa y con derivada no nula en todo
d aseguran la existencia de una correspondencia inversa de D sobre d que podemos notar:
z = h(Z) (5.195)
Figura 5.37:
190 CAPÍTULO 5. FLUJO POTENCIAL
5.13.3.1. Definición
La transformación conforme de Joukovski es definida por la relación:
R2
Z =z+ (5.198)
z
donde R es un número real no nulo. Esta transformación transforma un circulo, centrado al origen
y de radio R en el plano (x, y), en una elipse en el plano (X, Y ).
5.13.3.2. Demostración
Sea z = ρeiφ , entonces tenemos:
R2 R2 −iφ
Z = X + iY = z + = ρeiφ + e
z ρ
! !
R2 R2
= ρ+ cos φ + i ρ − sin φ (5.199)
ρ ρ
Eliminando el ángulo φ entre las expresiones de X y Y , obtenemos:
X2 Y2
2 2
+ 2 =1 (5.200)
R R2
ρ+ ρ
ρ− ρ
5.13. TRANSFORMACIÓN CONFORME 191
La Ec. 5.200 es una elipse donde los dos focos P1 y P2 son posicionados sobre el eje Ox en u = ±2R.
Sobre la Figure 5.39 se presenta la terminologı́a usada en aeronáutica para describir un perfil
alar.
192 CAPÍTULO 5. FLUJO POTENCIAL
Figura 5.41:
Estos perfiles presentan: un borde de ataque redondo, un borde de fuga de ángulo nulo, al
contrario de los perfiles alar reales. El perfil de Joukovski presenta un espesor y una curvatura en
función:
5.14. TEORÍA DE LOS PERFILES A PUNTA 193
Figura 5.42:
cuerda c ≈ 4k
e
√
Espesor relativa c
≈ −0, 75 3 xkc
f
flecha relativa c
≈ 0, 5 ykc
Γ
sin (β − α) = (5.201)
4πaV0
5.14. TEORÍA DE LOS PERFILES A PUNTA 195
Figura 5.44:
un único valor Γ para el cual la velocidad permanece finita en el borde de fuga, que se puede
encontrar precisamente imponiendo esta condición.
5.14.4.1. Definición
La sustentación de un perfil alar de Joukovski posicionado en un flujo bidimensional plano
estacionario de fluido ideal incompresible con circulación adaptada es único. Su modulo, por unidad
de envergadura, es:
|F⊥ | = ρV∞ |Γ| = 4πρV∞2 a| sin (β − α) | (5.203)
Ası́, la teorı́a de Kutta-Joukovski de los perfiles a punta permite dar un valor y una sola de la
sustentación de un perfil para una incidencia dada.
196 CAPÍTULO 5. FLUJO POTENCIAL
Figura 5.45: Para varias circulaciones, comportamiento del flujo alrededor de un perfil alar.
Figura 5.46:
realidad no tenemos rotación del solido. Cuales son los mecanismos que pueden en fluido viscoso
justificar igualmente una circulación de las velocidades alrededor del perfile.
El borde de fuga de un perfil a punta obtenido por transformación de Joukovski proviene de
un punto de estancamiento aguas arriba del flujo alrededor del circulo generador.
Figura 5.47:
Para fijar las ideas, supongamos que el perfil que se mueve en el fluido al instante t0
En un fluido viscoso, el desplazamiento del solido genera el desplazamiento del fluido que adhiera
a la pared y después, por viscosidad del fluido mas lejos.
Para los primeros instantes muy cercanos de t0 :
Figura 5.48:
Sin embargo, inicialmente el fluido esta en reposo, un circulación de las velocidades al largo
del contorno C1 es nula Γ = 0.
c) configuración adaptada
198 CAPÍTULO 5. FLUJO POTENCIAL
Figura 5.49:
Figura 5.50: Para un perfile NACA 0012, comparación entre los resultados analı́ticos y los resul-
tados experimentales (en función del angulo de ataque α).
Figura 5.52: Para varios ángulos de ataque α, comparación entre la fuerza de arrastre (CD ) y la
fuerza de sustentación (CL ).
Capı́tulo 6
6.1. Introducción
6.1.1. Objetivos
En estas clases vamos estudiar los flujos unidimensional y suponer que no hay: intercambio
de calor, de fricción a la pared, de inyección de masa o de reacciones quı́micas en el seno del flujo
de gas.
Vamos a tratar sucesivamente de la dinámica de los flujos isentropicos unidimensional estacio-
nario, de los flujos compresibles a grande velocidad (ondas de choques derechas y oblicas) y de los
flujos en las toberas y túneles de viento.
Vamos a tratar primero la dinámica de los flujos isentropicos unidimensionales estacionarios.
Estos flujos juegan un papel importante en numerosos problemas tecnológicas: los flujos en con-
ductos y los flujos externos en los tubos de corriente elementales.
201
202 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
Figura 6.2:
Figura 6.3:
6.1.4. Comentarios
En estas clases vamos a considerar los flujos isentrópicos:
de gas reales
• en efecto se puede obtener varios resultados cualitativos sin precisar las formas de la
ecuación de estado y de la entalpı́a.
204 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
v2
h+ = Cte (6.4)
2
Generalmente, llamamos la constante:
v2 v2
hi = h + = Cp T + = Cp Ti (6.5)
2 2
la entalpı́a de remanso o de estancamiento.
la cantidad Cp Ti representa la entalpı́a de estancamiento isentrópica que corresponde al
2
valor que toma la entalpı́a total h + v2 en un punto donde la velocidad se anula de manera
adiabática y reversible.
Es importante notar que la entalpı́a h no puede superar la entalpı́a de estancamiento
hi .
Comentarios: A partir de la ecuación:
v2 v2
hi = h + = Cp T + = Cp Ti (6.6)
2 2
se puede determinar la velocidad local en función de la entalpı́a y de la entalpı́a de estancamiento:
q
v= 2 (hi − h) (6.7)
206 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
Para este elemento se puede expresar la fuerza parietal que se ejerce sobre δΣ según x:
dp + ρvdv = 0 (6.15)
6.2. FLUJO ISENTRÓPICO UNIDIMENSIONAL DE GAS REALES EN CONDUCTOS DE SECCIONES
dp + ρvdv = 0 (6.17)
Ecuación de energı́a:
v2
h+ = hi = Cte (6.18)
2
Las ecuaciones 6.16, 6.17, 6.18 deben ser completadas en el caso de flujos isentrópicos (adia-
báticos reversibles) para las ecuaciones que indica que la entropı́a es constante:
dS = 0 (6.19)
Además, se debe tener: Ecuaciones de estado:
h = h(s, p)
ρ = ρ(s, p) (6.20)
Las ecuaciones de estado del fluido pueden ser especificadas en forma de tablas, de diagramas o de
expresiones algébricas. Cuando las ecuaciones de estado tienen formas algebraicas simples (caso
de los gases ideales), se puede obtener soluciones completas.
v12
hi = h1 + (6.21)
2
208 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
Sabemos que esta ultima cantidad es constante si el fluido es adiabático. Además, en este caso
(por definición) las fuerzas de fricción y el carácter reversible o irreversible del flujo no interviene.
Consideramos ahora una expansión a partir del estado 1 (Ver Figura 6.9). Si la expansión es
isentrópica, el punto 2 se encuentra vertical del punto 1. El estado termodinámico es: p2 , ρ2 , T2 ,
h2 , s2 y su velocidad v2 , la velocidad en el estado 2 es dada por:
v22
= hi − h2 (6.26)
2
Si consideramos ahora una expansión real a partir del estado 1. Los efectos de la irreversibilidad
se manifiestan, y el proceso se aleja de la isentrópica. La entropı́a aumenta y el estado final
20 se posiciona a la derecha de 2. La velocidad final en el estado 20 es dada por:
0
v22 0
= hi − h2 (6.27)
2
Se puede mostrar fácilmente que la velocidad v2 obtenida después una expansión isentropica
0
es superior a v2 que resulta de la expansión real.
Por eso, a partir de:
0
v22 v22 0
= hi − h2 y = hi − h2 (6.28)
2 2
escribimos:
0 0
v22 − v22 = 2(h2 − h2 ) (6.29)
0
Se queda a mostrar que h2 > h2 .
0
Para mostrar que h2 > h2 se usa:
la relación de Gibbs:
1
dh = T ds + dp (6.30)
ρ
sea a presión constante, la variación de entalpı́a h con la entropı́a s es positiva por qué se
escribe: !
∂h
=T (6.31)
∂s p
210 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
Finalmente, mostramos que la diferencia de las energı́a cinéticas correspondiente a los estados
2 y 20 es igual a las sumas de los trabajos disipados (gastados) durante el proceso irreversible.
Tenemos:
δQ = T ds − δWf (6.34)
donde δWf > 0.
Como el proceso es adiabático δQ = 0, tenemos:
T ds = δWf (6.35)
vemos que la diferencia de las energı́as cinéticas es igual a la suma de los trabajos disipados
(gastados) durante el proceso adiabático irreversible.
6.2. FLUJO ISENTRÓPICO UNIDIMENSIONAL DE GAS REALES EN CONDUCTOS DE SECCIONES
dρ dA dv
+ + =0 (6.40)
ρ A v
En flujo incompresible, tenemos:
dρ = 0 (6.41)
y la ecuación 6.40 indica que: el crecimiento de la velocidad es inversamente proporcional a la
disminución de la sección. En efecto tenemos:
dA dv
+ =0 (6.42)
A v
dρ
Cuando el fluido es compresible la ecuación 6.40 no puede ser simplificada se debe expresar ρ
.
Para eso tomamos la ecuación de conservación de la cantidad de movimiento:
dp + ρvdv = 0 (6.43)
podemos escribir:
dp dv
+ v2 =0 (6.44)
ρ v
dρ
Para hacer aparecer la variación de masa especifica del fluido, ρ
, escribimos:
dp dp dρ
= (6.45)
ρ dρ ρ
Por definición: !
∂p
(6.47)
∂ρ s
dp + ρvdv = 0 (6.59)
y dividamos por p
dp ρv 2 dv
+ =0 (6.60)
p p v
Deducimos de la ecuación 6.58, la siguiente expresión:
6.2. FLUJO ISENTRÓPICO UNIDIMENSIONAL DE GAS REALES EN CONDUCTOS DE SECCIONES
dp 1 ρv 2 dA
= (6.61)
p 1 − M2 p A
Los efectos de cambio de sección sobre la velocidad y la presión pueden ser descrito considerando
3 casos:
Por M ≈ 0
• una disminución de la sección produce proporcional de la velocidad
• y también una disminución de la presión.
Por 0 < M < 1
• Es decir a las velocidades subsonicas
• El comportamiento es cualitativamente el mismo
• Una disminución de la sección produce un crecimiento de la velocidad
• La variación de la velocidad es mas importante por que 1 − M 2 < 1
• La presión disminuye es este caso también con la sección.
Por M > 1
• Es decir a las velocidades supersónicas
• 1 − M2 < 0
• un crecimiento de la velocidad es producida por un crecimiento de la sección.
En efecto, cuando el flujo es supersónico, la evolución de la velocidad es diferente de la evolución
que conocemos por el caso incompresible. A velocidades supersónicas, la masa especifica disminuye
mas rápidamente que la velocidad
dρ dv
= −M 2 (6.62)
ρ v
ası́ la sección debe crecer para asegurar la conservación de la masa.
Figura 6.11:
Cuando el flujo desacelera tenemos un difusor para acelerar un gas de una velocidad subsonica a
una velocidad supersónica, se necesita una sección convergente seguida por una sección divergente
(Ver Figura 6.11). Para desacelerar un gas en flujo supersónico, el conducto debe ser primero
convergente (Ver Figura 6.12).
214 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
Figura 6.12:
Figura 6.13:
dv −1 dA
= (6.63)
v 1 − M2 A
Ti v2
=1+ (6.79)
T 2cp T
Ti v2
=1+ (6.80)
T 2cp T
v2
El ratio 2cp T
se puede expresar en función del número de Mach (solamente el número de
Mach):
v2 v2 k − 1 v2 k−1 2
= kr = 2
= M (6.81)
2cp T 2 k−1 T 2 c 2
Ti
Entonces el ratio T
:
Ti k−1 2
=1+ M (6.82)
T 2
A partir de las relaciones isentrópicas anterior (Ec. 6.82) se puede deducir:
6.3. FLUJO ISENTRÓPICO UNIDIMENSIONAL DE GAS IDEALES EN CONDUCTOS DE SECCIONE
para la presión
k ! k
pi Ti k−1 2 k−1
k−1
= = 1+ M (6.83)
p T 2
1 ! 1
ρi Ti k−1 2 k−1
k−1
= = 1+ M (6.84)
ρ T 2
Las 3 relaciones TTi , ppi , ρρi son determinadas en función del número de Mach. La temperatura y la
presión de estancamiento Ti y pi son constantes en flujo isentrópico y por esta razón es cómodo
expresar el flujo másico en función de Ti , pi y la sección A y del número de Mach.
ṁ∗ = ṁ (6.93)
218 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
y
" !# k+1
A 1 2 k−1 2 2(k−1)
= 1+ M (6.94)
A∗ M k+1 2
Se puede mostrar que esta expresión tiene un mı́nimo en M = 1. Para un flujo isentrópico
dado, la sección critica es mı́nima. Esta expresión interviene en númerosos cálculos practicos.
6.3.2.2. Resumen
Las ecuaciones
Ti k−1 2
=1+ M (6.95)
T 2
k ! k
pi Ti k−1 2 k−1
k−1
= = 1+ M (6.96)
p T 2
1 ! 1
ρi Ti k−1 2 k−1
k−1
= = 1+ M (6.97)
ρ T 2
y
" !# k+1
A 1 2 k−1 2 2(k−1)
= 1+ M (6.98)
A∗ M k+1 2
forman un juego completo de ecuaciones que permite determinar las caracterı́sticas de los flujos
isentrópicos de gases ideales en los conductos de sección variable.
v2 v 2 c2 2c
2
M∗2 = = = M (6.100)
c2∗ c2 c2∗ c2∗
c2 T
El ratio c2∗
es igual al ratio T∗
y se puede escribir:
c2 T T Ti
2
= = (6.101)
c∗ T∗ Ti T∗
2
expresamos ahora cc2∗ .
Los 2 ratio presentados a la derecha de la Ec. 6.101 son dados por:
T 1
= k−1 (6.102)
Ti 1 + 2 M2
Ti Ti 1+k
= = (6.103)
T∗ T M =1 2
6.3. FLUJO ISENTRÓPICO UNIDIMENSIONAL DE GAS IDEALES EN CONDUCTOS DE SECCIONE
Sea:
c2 k+1
= (6.104)
2
c∗ 2 + (k − 1)M 2
Finalmente, se deduce:
M 2 (k + 1)
M∗2 = (6.105)
2 + (k − 1)M 2
6.3.4. Resumen
Condiciones de estancamiento para un fluido ideal:
Entalpı́a de estancamiento:
v2
hi = h + (6.106)
2
Temperatura de estancamiento:
!
hi v2 k−1 2
Ti = =T+ =T 1+ M (6.107)
cp 2cp 2
La temperatura
2
T∗ = Ti (6.110)
k+1
La presión
k
2
k−1
p∗ = pi (6.111)
k+1
y masa especifica
1
2
k−1
ρ∗ = ρi (6.112)
k+1
220 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
Figura 6.14:
6.3.7. Ejemplo
Consideramos un túnel de viento supersónico que funcione a un número de Mach M = 2. La
superficie de la sección del cuello es Ac = 0, 1m2 (Ver Figura 6.16). La presión y la temperatura
generadora en el estanque posicionado agua arriba de la tobera son: p0 = 3bar y T0 = 300K
El análisis será suponiendo que el flujo entre el estanque y la sección de ensayo se queda isentró-
pico. Se considera igualmente los parámetros termofı́sicos siguientes: k = 1, 4 y r = 287J.kg−1 .K−1 .
1. Calcular las propiedades del flujo al cuello y en la sección de ensayo (se determinar en
particular ρ, p, T , c, v y la sección A).
2. Determinar el flujo de gas en el túnel de viento.
6.3. FLUJO ISENTRÓPICO UNIDIMENSIONAL DE GAS IDEALES EN CONDUCTOS DE SECCIONE
Figura 6.15:
6.3.7.1. Solución
Determinar las propiedades del fluido en el estanque
p0
p0 , T0 , ρ = rT 0
= 3, 84kg.m−3
√
c0 = krT0 = 347m.s−1
222 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
Para obtener las propiedades del fluido al cuello de la tobera, se usa las tablas de flujo isentropica.
Al cuello tenemos M = 1
T∗ p∗ ρ∗
= 0, 8333 = 0, 5283 = 0, 6339 (6.119)
Ti pi ρi
Aquı́ las condiciones generadoras son también las condiciones de estancamiento isentrópico, por
que el flujo es supuesto también isentropico en toda el túnel de viento. En estas condiciones:
T∗ = 0, 8333Ti = 250 K
p∗ = 0, 5283pi = 1, 58 bar
Consideramos la sección de ensayo del túnel. En esta sección tenemos M = 2. Las tablas dan:
T p ρ
= 0, 555 = 0, 128 = 0, 230 (6.120)
Ti pi ρi
A
= 1, 687 (6.121)
A∗
En estas condiciones:
T = 0, 555Ti = 166 K
v = M c = 517m.s−1
6.3. FLUJO ISENTRÓPICO UNIDIMENSIONAL DE GAS IDEALES EN CONDUCTOS DE SECCIONE
A = 1, 687A∗ = 0, 1687m2
ṁ = ρ∗ A∗ v∗ ≈ 70 kg.s−1
El flujo másico puede ser también calculado en la sección de ensayo ṁ = ρAv ≈ 70 kg.s−1 . Los
dos valores son iguales.
dρ dA dv
+ + =0 (6.122)
ρ A v
dp ρv 2 dv
+ =0 (6.123)
p p v
dv −1 dA
= (6.127)
v 1 − M2 A
Encontramos la expresión que establecimos anteriormente en el caso de los gases reales. Se puede
deducir:
dp kM 2 dA
= (6.128)
p 1 − M2 A
y
dρ M 2 dA
= (6.129)
ρ 1 − M2 A
dT
El ratio T
puede se obtener a partir de:
dT k − 1 dp (k − 1)M 2 dA
= = (6.130)
T k p 1 − M2 A
224 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
Figura 6.17:
Vamos a mirar casos de flujos tı́picos donde se encuentran ondas de choque. Por ejemplo, en
un canal divergente, la onda de choque no es perfectamente derecho. Se observa dos bifurcaciones
asociados al interacción entre el choque y las capa limites parietales (Ver Figura 6.18).
Dibujo a partir obtenidos a partir de strioscopias (Ver Figura 6.19). Para este choque, la razón
de la presión antes y después el choque es del orden de 5. La velocidad aguas abajo es alrededor
de 3 veces mas pequeña que la velocidad aguas arriba.
Figura 6.18: Onda de choque en un canal divergente. La imagen formada es negra en ausencia de
turbulencia y luminosa donde hay turbulencia.
Estructura del flujo alrededor de un proyectil lanzado a grande velocidad (Ver Figura 6.21).
Entrada de aire de un turbo-reactor supersónico (Ver Figura 6.23). La onda de choque cónica
(oblicua) generado por el cuerpo central genera una compresión externa. Una onda de choque
derecha aparece al nivel del entrada de aire.
En comparación con el caso precedente. La intensidad de la onda de choque es mas fuerte
cuando no hay cuerpo central (Ver Figura 6.24).
Sistema de onda de choque en salida de una tobera convergente-divergente (Ver Figura 6.26).
La presión de eyección es inferior a la presión ambiente. Los gases son re-comprimidos por el
226 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
Figura 6.20: Onda de choque despegada cerca de un cuerpo desafilado posicionado en un flujo a
gran velocidad.
Figura 6.21: Estructura del flujo alrededor de un proyectil lanzado a grande velocidad.
intermediario del choque. La interacción de este choque produce un disco de Mach y una sucesión
de ondas de compresión y de expansión.
Onda de choque en un conducto. La capa limite que se desarrolla cerca de las paredes tiene
6.4. ONDAS DE CHOQUE: GENERALIDADES 227
Figura 6.23:
Figura 6.24:
de presión no pueden seguir el rápido movimiento de la fuente que las originó. Todas ellas son
arrastradas detrás de la fuente y son tangentes a una superficie llamada cono de Mach. La frontera
del cono se llama linea o onda de Mach. El medio ángulo α a la cima del cono es (Ver Figura 6.31):
c 1
sin α = = (6.131)
U M
Si M = 2, α = 30◦ .
Si M = 5, α ≈ 11◦ .
Si M = 1, α = 90◦ .
Figura 6.28:
Figura 6.29:
Figura 6.30:
Figura 6.31:
6.5.1.1. hipótesis
Las hipótesis usadas en esta análisis son las siguientes: Las fuerzas de fricción en el volumen
de control que rodea la onda de choque son despreciables. El proceso de choque es realizado a
sección constante, las lineas de corriente se quedan paralelas entre ellas y perpendiculares a la
superficie del choque. El flujo es adiabatico y no hay trabajo exterior sobre el volumen de controle
considerado. Las fuerzas de volumen son despreciables. El flujo se desarrolla en ausencia de
reacciones quı́micas y el fluido se queda en equilibrio termodinámico.
Sin embargo, la onda de choque es una perturbación que hace perder el carácter reversible,
pero conserva el carácter adiabático del flujo.
ṁ = ρ1 V1 A1 = ρ2 V2 A2 = Cte (6.132)
si A1 = A2 = A se puede escribir:
ṁ
= ρ1 V1 = ρ2 V2 (6.133)
A
Figura 6.32:
v12 v22
h1 + = h2 + = hi = Cte (6.136)
2 2
La energı́a interna U y la entalpia h son relacionada por:
p
h=U+ (6.137)
ρ
6.5. ONDAS DE CHOQUE DERECHAS O NORMALES 233
p1 v12 p2 v22
U1 + + = U2 + + = Cte (6.138)
ρ1 2 ρ2 2
6.5.1.6. Resumen
Tenemos:
ṁ = ρ1 V1 A1 = ρ2 V2 A2 = Cte (6.141)
p1 + ρ1 v12 = p2 + ρ2 v22 (6.142)
v12 v2
h1 += h2 + 2 = hi = Cte (6.143)
2 2
P1 P2
s = s(p, ρ) h = h(p, ρ) o = (6.144)
ρ1 T1 ρ2 T2
Tenemos 4 ecuaciones para determinar 4 parámetros no conocidos:
p2 , ρ2 , v2 , h2
Las ecuaciones permiten el estudio de las ondas de choque en los fluidos reales. Sin embargo, es
mas simple de continuar el análisis en el caso de los gases ideales. Esta hipótesis permite
escribir relaciones explicitas entre las caracterı́sticas antes y después el choque.
p = ρrT h = cp T (6.147)
ρ1 V1 = ρ2 V2 (6.148)
v12 v22
cp T1 + = cp T2 + (6.150)
2 2
Es cómodo de hacer parecer en estas 3 relaciones los números de Mach M1 , M2 de los flujos
antes y después el choque. Una vez estas expresiones obtenidas, se podrı́a establecer una relación
entre el numero de Mach M1 (en aguas arriba) y M2 (en aguas abajo).
Después, sera posible presentar todos los ratios entre los parámetros antes y después el
choque como:
p2 T2 ρ2 pi2
(6.151)
p1 T1 ρ1 pi1
en función de los números de Mach (antes o después el choque). Las expresiones deducidas son
muy útiles en practica y permiten la resolución de numerosos problemas de ingenierı́a.
v12 v2
cp T1 + = cp T2 + 2 (6.152)
2 2
El primer miembro de esta ecuación puede se escribir:
!
v2 v2
cp T1 + 1 = cp T1 1 + 2 (6.153)
2 2cp T1
c2
sea como Cp T = k−1
:
!
v2 k−1 2
cp T1 + 1 = cp T1 1 + M1 (6.154)
2 2
Una operación análoga puede ser hecha en el segundo miembro
!
v2 k−1 2
cp T2 + 2 = cp T2 1 + M2 (6.155)
2 2
Sea:
k−1
T2 1+ 2
M12
= k−1 (6.157)
T1 1+ 2
M22
kp1
como c21 = ρ1
p1 + ρ1 v12 = p1 1 + kM12 (6.159)
De la misma manera por el miembro de la derecha y la ecuación ?? se escribe
p1 1 + kM12 = p2 1 + kM22 (6.160)
Sea todavı́a:
p2 1 + kM12
= (6.161)
p1 1 + kM22
y
2
M12 + k−1 (k − 1)M12 + 2
M22 = 2kM12
= (6.166)
−1 2kM12 + 1 − k
k−1
expresa que no hay cambio entre el estado 1 (aguas arriba) y el estado 2 (aguas abajo). Esta raı́z
es trivial, no hay choque
6.5.2.3.2. La segunda raı́z de Ec. 6.164 es asociada a una onda de choque derecha.
(k − 1)M12 + 2
M22 = (6.168)
2kM12 + 1 − k
Vamos a ver mas allá, que el numero de Mach M1 del flujo (aguas arriba) debe ser (es obli-
gatorio) superior a 1. las ondas de choques no pueden formarse en flujo subsonico (M < 1). La
relación 6.168 indica que M2 (aguas abajo) es siempre inferior a 1.
236 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
Figura 6.33:
6.5.3.1. La temperatura
Tenemos:
k−1
T2 1+ 2
M12
= k−1 (6.171)
T1 1+ 2
M22
substituyendo 6.166 tenemos un salto de temperatura a través la onda de choque derecha:
!" #
T2 2k k−1 k−1 2
= M12 − + (6.172)
T1 k+1 k+1 k + 1 (k + 1)M12
6.5. ONDAS DE CHOQUE DERECHAS O NORMALES 237
6.5.3.2. La presión
Substituyendo 6.166 en la expresión
s
p2 M1 T2
= (6.173)
p1 M2 T1
permite determinar el salto de presión a través la onda de choque derecha:
p2 2k k−1 2k 2
= M12 − =1+ M1 − 1 (6.174)
p1 k+1 k+1 k+1
ρ2 p2 T1 (k + 1)M12
= = (6.175)
ρ1 p1 T2 2 + (k − 1)M12
Esta expresión es también el ratio de las velocidades, una consecuencia directa de la ecuación de
conservación de la masa:
v1 (k + 1)M12
= (6.176)
v2 2 + (k − 1)M12
Figura 6.34:
6.5.5. Ejemplo
Un flujo de aire presenta un numero de Mach M1 = 2, una presión de P1 = 1 bar y una
temperatura T1 = 300K. Este flujo pasa a través de un choque derecho estacionario.
Figura 6.35:
p2 ρ2
Esta relación puede ser obtenida a partir de las expresiones p1
y de ρ1
:
p2 2k k−1
= M12 − (6.177)
p1 k+1 k+1
ρ2 (k + 1)M12
= (6.178)
ρ1 2 + (k − 1)M12
Se puede calcular a partir de la Ec. 6.177 el numero de Mach M12 :
(k + 1)p2 + (k − 1)p1
M12 = (6.179)
2kp1
Ponemos la Ec. 6.179 en la Ec. 6.178:
ρ2 (k + 1)p2 + (k − 1)p1
= (6.180)
ρ1 (k + 1)p1 + (k − 1)p2
p2
Se puede ası́ expresar p1
:
p2 (k + 1)ρ2 − (k − 1)ρ1
= (6.181)
p1 (k + 1)ρ1 − (k − 1)ρ2
que podemos escribir:
k+1 ρ2
p2 k−1 ρ1
−1
= (6.182)
p1 k+1
k−1
− ρρ21
Esta ultima relación (Ec. 6.182) es la relación de Rankine-Hugoniot que relaciona las razones de
presión y de masa especifica a través una onda de choque derecha. Esta relación pone en evidencia
el carácter no isentropico del flujo a través la onda de choque. En efecto, en el caso isentropico,
presión y masa especifica de gas ideal son relacionado por la ecuación de Laplace:
!−k
p2 ρ1
= (6.183)
p1 ρ2
Dibujamos las relaciones pp12 = f ρρ12 (Ver Figura 6.36). La divergencia entre las 2 curvas puede
se explicar que por el carácter no reversible del flujo a través la onda de choque, porque
las 2 relaciones se aplican a las evoluciones adiabáticas.
La relación de Hugoniot 6.182 permite mostrar que cuando:
p2
→∞ (6.184)
p1
El ratio de las masas especificas tiende hacia un limite finito. El limite de los grandes números de
Mach:
ρ1 k−1
lı́m = (6.185)
p2
p1
→∞ ρ2 k+1
ρ1
Para un gas caracterizada por k = 1, 4, este ratio no puede ser mas grande que el valor ρ2
= 0, 16.
La linea de Hugoniot y la isentropica son tangentes al punto:
p2 ρ2
=1 y =1 (6.186)
p1 ρ1
Ası́, el caso de onda de choque débil, es decir:
p2
1 y M1 ∼ 1 (6.187)
p1
240 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
Si remplazamos en 6.188
p2 2k 2
=1+ M1 − 1 (6.189)
p1 k+1
ρ2 (k + 1)M12
= (6.190)
ρ1 2 + (k − 1)M12
Se obtiene la variación de entropı́a a través el choque derecho, en función del número de Mach
aguas arriba M1
!k
2kM12
+1−k 2 + (k − 1)M12
s2 − s1 = Cv ln (6.191)
k+1 (k + 1)M12
6.5. ONDAS DE CHOQUE DERECHAS O NORMALES 241
6.5.8.1. La presión
La razón entre la presión de estancamiento pi2 y la presión estática p1 interviene cuando se
hace mediciones de presión en flujo supersónico con un tubo de Pitot.
El tubo mide la presión de estancamiento isentrópico después el choque derecho que aparece en
aguas arriba de la sonda (Ver Figura 6.37). El tubo de Pitot medı́ pi2 y no pi1 . En estas condiciones
es interesante determinar:
pi2 pi2 p2
= (6.192)
p1 p2 p1
Figura 6.37:
pi2
La razón p2
es dado por la relación isentrópica:
! k
p2i k−1 2 k−1
= 1+ M2 (6.193)
p2 2
p2
y p1
es determinado en función de M1 por:
p2 2k k−1
= M12 − (6.194)
p1 k+1 k+1
Se queda remplazar M2 en función de M1 en la Ec. 6.193 y obtenemos:
1
k+1
M12
! " #
k−1
pi2 k+1
= M12 2
(6.195)
p1 2 2k
k+1
M12 − k−1
k+1
242 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
Una otra relación útil es la razón entre las presiones de estancamiento isentrópico en aguas arriba
y en aguas abajo:
pi2 pi2 p2 p1
= (6.196)
pi1 p2 p1 pi1
Se obtiene fácilmente:
" −1
#k ! k−1
pi2 k − 1 2 2k 2 k−1
= + M − (6.197)
pi1 k + 1 (k + 1)M12 k+1 1 k+1
6.5.8.2. La temperatura
La ecuación de la energı́a se escribe:
v12 v2
Cp T1 + = Cp T2 + 2 (6.198)
2 2
sea como tenemos condicion de estancamiento v1 = v2 = 0
muestra que es una función monótona creciente del número de Mach M1 . Sea el gráfico cerca de
1 (ver Figura 6.38).
La única variación de entropı́a compatible con el segundo principio de la termodinámica, es
decir, s2 > s1 , es para M > 1.
No se puede tener onda de choque normal en régimen subsonico.
Figura 6.38:
Figura 6.39: subı́ndice 0 aguas arriba (antes el choque). subı́ndice 1 aguas abajo (después el
choque).
244 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
6.6. Resumen
6.6.1. Ecuaciones fundamentales para las ondas derechas
Ecuación de conservación de la masa
ṁ
= ρ1 V1 = ρ2 V2 (6.202)
A
Ecuación de conservación de la cantidad de movimiento
v12 v2
h1 + = h2 + 2 = hi (6.204)
2 2
Ecuación del estado
6.6.2. Ecuaciones para las propiedades del flujo a través una onda de
choque derecha (caso gas ideal)
2
M12 + k−1
M22 = 2kM12
(6.206)
k−1
−1
p2 2k k−1
= M12 − (6.207)
p1 k+1 k+1
ρ2 v1 (k + 1)M12
= = (6.208)
ρ1 v2 2 + (k − 1)M12
!" #
T2 2k k−1 k−1 2
= M12 − + (6.209)
T1 k+1 k+1 k + 1 (k + 1)M12
" −1
#k ! k−1
pi2 k − 1 2 2k 2 k−1
= + M − (6.210)
pi1 k + 1 (k + 1)M12 k+1 1 k+1
1
k+1
M12
! " #
k−1
pi2 k+1
= M12 2
(6.211)
p1 2 2k
k+1
M12 − k−1
k+1
k+1 ρ2
p2 k−1 ρ1
−1
= (6.213)
p1 k+1
k−1
− ρρ21
6.7. ONDAS DE CHOQUE OBLICUAS 245
Figura 6.40:
Para facilitar, los cálculos se supone que el flujo pasa a través una onda de choque oblicua plana.
ε es el ángulo de inclinación entre el plano de la onda de choque y la dirección del flujo aguas
arriba V1 .
Vn1 y Vt1 los componentes normal y tangencial de la velocidad aguas arriba del choque.
Vn2 y Vt2 los componentes normal y tangencial de la velocidad aguas abajo del choque.
V12 = Vn1
2
+ Vt12 Vt1 = V1 cos ε Vn1 = V1 sin ε (6.218)
V22 = Vn2
2
+ Vt22 Vt2 = V2 cos β Vn2 = V2 sin β (6.219)
Se puede ahora escribir las ecuaciones fundamentales.
Figura 6.41:
2 2
p1 + ρ1 Vn1 = p2 + ρ2 Vn2 (6.221)
2
Vn1 + Vt12 V 2 + Vt22
h1 + = h2 + n2 (6.224)
2 2
Como las velocidades tangencial son iguales:
6.7.2.6. Resumen
Ecuación de conservación de la masa
2
vn1 v2
h1 + = h2 + n2 (6.231)
2 2
Ecuación de estado
Relaciones geométricas
Como antes vamos analizar el flujo por el caso de los gases ideales.
p = ρrT h = cp T (6.235)
V1 = Vn1
V2 = Vn2
De la misma manera, se puede remplazar los números de Mach M1 por Mn1 y M2 por Mn2 en
las expresiones establecidas en el caso de una onda de choque derecha.
Los componentes normales del numero de Mach son dados por:
Vn1 Vn2
Mn1 = ≡ M1 sin ε Mn2 = ≡ M2 sin (ε − δ) (6.239)
c1 c2
Entonces se puede escribir directamente:
p2 2k 2 k−1
= Mn1 − (6.240)
p1 k+1 k+1
6.7. ONDAS DE CHOQUE OBLICUAS 249
2
ρ2 (k + 1)Mn1
= 2
(6.241)
ρ1 2 + (k − 1)Mn1
!" #
T2 2k 2 k−1 k−1 2
= Mn1 − + 2
(6.242)
T1 k+1 k+1 k + 1 (k + 1)Mn1
Sea:
p2 2k k−1
= (M1 sin ε)2 − (6.243)
p1 k+1 k+1
ρ2 (k + 1)(M1 sin ε)2
= (6.244)
ρ1 2 + (k − 1)(M1 sin ε)2
!" #
T2 2k k−1 k−1 2
= (M1 sin ε)2 − + (6.245)
T1 k+1 k+1 k + 1 (k + 1)(M1 sin ε)2
Estas relaciones ponen una diferencia entre onda de choque derecha y onda de choque oblicua:
Choque derecho: depende solamente del número de Mach aguas arriba.
Choque oblicuo depende de 2 parámetros: del número de Mach aguas arriba y del ángulo de
choque ε.
Además tenemos:
vt1 = vt2 (6.248)
Poniendo, las relaciones anteriores en la ecuación de conservación de masa, se obtiene:
tan ε ρ2
= (6.249)
tan(ε − δ) ρ1
Ası́, tenemos una relación que permite determinar el ángulo de choque ε en función del numero
de Mach M1 y de la desviación δ:
ρ2 (k + 1)(M1 sin ε)2 tan ε
= = (6.250)
ρ1 2 + (k − 1)(M1 sin ε)2 tan(ε − δ)
Se puede hacer aparecer la deflexión δ de manera mas explicita en la ecuación 6.250, realizando
los cálculos siguientes:
tan ε tan ε(1 + tan ε tan δ) (k + 1)(M1 sin ε)2
= = (6.251)
tan(ε − δ) tan ε − tan δ 2 + (k − 1)(M1 sin ε)2
250 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
Escribamos esta relación con el mismo denominador y reagrupamos los términos en tan δ en el
primero miembro:
h i
tan δ (k + 1)M12 sin2 ε + (k − 1)M12 sin2 ε tan2 ε + 2 tan2 ε
= 2 tan ε M12 sin2 ε − 1 (6.252)
Algunos cálculos permiten escribir esta relación en una forma final:
!
1 k+1 M12
= − 1 tan ε (6.253)
tan δ 2 M12 sin2 ε − 1
Ası́, para determinar el el ángulo de choque ε en función del numero de Mach M1 y de la desviación
δ podemos usar:
ρ2 (k + 1)(M1 sin ε)2 tan ε
= = (6.254)
ρ1 2 + (k − 1)(M1 sin ε) 2 tan(ε − δ)
o !
1 k+1 M12
= − 1 tan ε (6.255)
tan δ 2 M12 sin2 ε − 1
Si fijamos el número de Mach aguas arriba a M1 = 1, 8 en la Ec. 6.255 tenemos el gráfico de la
Figura 6.42.
Figura 6.42:
M1 sin ε ≥ 1 (6.258)
Figura 6.43:
252 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
Si dibujamos el ángulo de choque ε en función del número de Mach a desviación fijada (ver
Figura 6.44). Cuando:
M1 sin ε = 1 (6.259)
la relación 6.8 muestra que tan δ = 0, es decir de la deflexión del flujo es infinitamente pequeña.
La onda de choque oblicua es llamada onda de Mach. Su dirección es dado por:
1
sin ε = (6.260)
M1
y las presiones, densidad, velocidad, temperatura tienen variaciones de amplitud infinitamente
pequeña:
p2 ρ2 v1 T2
=1 =1 =1 =1 (6.261)
p1 ρ1 v2 T1
Vamos a volver mas allá sobre las ondas de Mach y sobre las ondas de choque oblicuas débiles
Figura 6.44:
El ángulo del choque ε es función de M1 y δ , es entonces posible de eliminar ε usando la Ec. 6.8.
Se obtiene una red de curvas dando el valor del número de Mach aguas abajo M2 en función del
flujo aguas arriba (M1 ) para un ángulo de desviación dado. Sea el gráfico siguiente (Ver Figura
6.45): Para distinguir mejor las propiedades de estas 2 soluciones, hacemos un zoom cerca del
Figura 6.45:
atrás un choque débil, el flujo es mayoritariamente supersónico, excepto en una región limi-
tada.
Además vemos para un ángulo de desviación δ dado, la onda de choque no existe bajo un valor
limite M1,limite del número de Mach aguas arriba.
Existe 2 regı́menes de choques oblicuos según que el flujo aguas arriba es supersónico (M2 > 1)
o subsonico (M2 < 1). Se puede mostrar que el caso donde el flujo se queda supersónico después
el choque es asociado a la más pequeña de los 2 valores posibles del ángulo de choque
254 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
Figura 6.46:
ε. La onda de choque débil aparece en la mayorı́a de los casos. La onda de choque fuerte es poca
comuna.
Por ejemplo (Ver Figura 6.47):
1. La onda de choque fuerte se forma a la punta de un cuerpo central de una entrada de aire
cuando la entrada de aire es completamente tapada
2. Entrada de aire en funcionamiento normal. La onda de choque es del tipo débil.
La figure siguiente da los valores de la razón pp12 en función del numero de Mach M1 y de la
deflexión δ. Si la razón de compresión es fijada, se puede leer sobre la figura la deflexión δ necesaria
para asegurar una compresión.
Figura 6.47:
Figura 6.48:
Se puede entonces leer en las tablas de choque derechos los valores de las razones de propiedades
correspondientes a Mn1 .
256 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
Figura 6.49:
6.7. ONDAS DE CHOQUE OBLICUAS 257
Si la deflexión δ y el numero de Mach aguas arriba son dados, pero no se conoce el ángulo de
choque ε. Se puede usar por ejemplo la Ec.
!
1 k+1 M12
= − 1 tan ε (6.267)
tan δ 2 M12 sin2 ε − 1
Esta ecuación permite obtener una expresión de M12 sin2 ε−1. Por eso, notamos que M12 sin2 ε−1
es del orden de δ cuando δ → 0. En este condición el segundo miembro de la derecha es del orden
δ 2 . Se puede escribir al orden δ:
!
k+1 2
M12 2
sin ε − 1 ≈ M1 tan ε δ (6.272)
2
Finalmente,
k + 1 M12
M12 sin2 ε − 1 ≈ q δ (6.274)
2 M12 − 1
La variación de presión relativa a través el choque es fácilmente deducida de este expresión. En
efecto, tenemos:
∆p p2 − p1 2k 2 2
= = M1 sin ε − 1 (6.275)
p1 p1 k+1
258 CAPÍTULO 6. DINÁMICA DE LOS GASES: FLUJOS COMPRESIBLES
y entonces tenemos:
∆p kM12
≈q δ (6.276)
p1 M12 − 1
Se puede también determinar la variación relativa de la velocidad asociada al choque:
∆v v2 − v1 v2
= = −1 (6.277)
v1 v1 v1
v2
La razón v1
se expresa en función de los ángulos ε y ε − δ
v2 cos ε
= (6.278)
v1 cos(ε − δ)
Ası́:
∆v cos ε − cos(ε − δ)
= (6.279)
v1 cos(ε − δ)
Cuando δ → 0 se puede escribir:
∆v sin ε
≈− δ (6.280)
v1 cos ε
Sea:
∆v 1
≈ −q δ (6.281)
v1 M12 − 1
La variación relativa de la velocidad es ası́ proporcional a la deflexión δ.
6.7.9. Resumen
Principales relaciones por las ondas de choque oblicuas:
p2 2k k−1
= (M1 sin ε)2 − (6.282)
p1 k+1 k+1
ρ2 vn1 (k + 1)(M1 sin ε)2
= = (6.283)
ρ1 vn2 2 + (k − 1)(M1 sin ε)2
!" #
T2 2k k − 1 k − 1 2
= (M1 sin ε)2 − + (6.284)
T1 k+1 k+1 k + 1 (k + 1)(M1 sin ε)2
!
1 k+1 M12
= − 1 tan ε (6.285)
tan δ 2 M12 sin2 ε − 1
2
M12 sin2 ε + k−1
M22 2
sin (ε − δ) = 2kM12 sin2 ε
(6.286)
k−1
−1
Principales relaciones por las ondas de choque oblicuas débiles con desviación muy pequeña:
1
tan ε ≈ q (6.287)
M12 − 1
k + 1 M12
M12 sin2 ε − 1 ≈ q δ (6.288)
2 M12 − 1
∆p kM12
≈q δ (6.289)
p1 M12 − 1
∆v 1
≈ −q δ (6.290)
v1 M12 − 1
Capı́tulo 7
Ayudantias
7.1. Capitulo 1
7.1.1. Problema 1: Flujo de Poiseuille
Queremos estudiar el flujo estacionario de un fluido incompresible de masa especifica ρ, de
viscosidad dinámica µ, en un conducto cilı́ndrico de radio R pequeño, de eje horizontal Oz y de
longitud L. Un punto en el conducto es localizado por sus coordenadas cilı́ndricas.
Supongamos que el campo de velocidad es ~v = vr (r, z)~ez y que el campo de presión es
P = P (r, z).
Supongamos que la influencia de la gravedad sobre el flujo es despreciable.
1. Mostrar con las condiciones evocadas, que si el fluido es incompresible, entonces tenemos:
−−→
~v .grad ~v = 0 (7.1)
259
260 CAPÍTULO 7. AYUDANTIAS
Supongamos que el flujo inducido sobre el aire por el movimiento del liquido es despreciable.
U0 = 0, 20m.s−1 .
P1 = 1200Pa y P2 = 1350Pa.
3. Con esta hipótesis escribir el sistema de ecuación y las condiciones a las limites que representa
el flujo en el espacio entre las dos cintas. Justificar su repuesta.
6. ¿Cual es el valor del gradiente de presión para que se produzca un flujo de retorno?
Preguntas
1. Se supone que el campo de velocidad es invariante por traslación según los ejes x y z
a) Mostrar usando la ecuación de conservación de masa que el vector velocidad se reduce
a 2 componentes no nulo: U (y) y W (y)
b) Deducir, usando la ecuación de conservación de momento según y, que la presión es
una función de x y de z.
c) Deducir, usando la ecuación de conservación de momento según x y z, que el gradiente
de presión en estas 2 direcciones es constante. Se notara:
1 dP
= Cte = A (7.4)
µ deβ
un flujo estacionario (permanente) de aire con un eje de simetrı́a de rotación vertical Oz.
un campo de velocidad ortoradial invariante por rotación y translación según 0z. Se nota:
~v = v(r)e~θ (7.6)
Por supuesto, un verdadero tornado tiene una zona de transición donde el rotacional decrece
progresivamente, pero esta zona tiene un espesor pequeña frente el radio a, eso valida nuestro
modelo.
1. Encontrar la expresión de v(r), si posible sin utilizar las formulas del rotacional en coorde-
nadas cilı́ndricas. Se tratará 2 casos: r < a y r > a
2. Supongamos que el aire sea un gas incompresible, calcular la presión en todo punto exterior
al cilı́ndrico de radio a.
7.1.6. Problema 6
Q
Cuatro fuentes de la misma intensidad m = 2πb (b es la longitud sobre z, aquı́ tomamos b = 1)
se sitúan en las cuatro esquinas de un cuadrado:
(x, y) = (a, a)
(x, y) = (−a, a)
7.1.7.1. Preguntas
Arriba de la astenosfera, se considera una placa tectónica solida llamada placa litosférica. La
placa litosferica tiene un espesor hl y un longitud L muy superior a hl y ha . Es decir: L hl
y L ha . Ademas, la placa litosferica se mueve según la dirección del eje horizontal Ox con
una velocidad constante y uniforme v0 . Eso genera en la astenosfera un flujo estacionario:
~v = u(z)~ex (7.12)
7.1. CAPITULO 1 265
7.1.8.1. Parte 1
Supongamos que las masas especificas de la litosfera y de la astenosfera son iquales: ρl ≈ ρa .
Sea P (x, z) la presión en la astenosfera.
6. Dar una nueva expresión del componente del tensor de esfuerzo viscoso que actúa sobre la
placa en movimiento en función de ha , hl , v0 y µ. Comparar este resultado con el resultado
obtenido a la parte 7.1.8.1.1 pregunta 4. Comentar.
7. Para alumnos con ánimo: Mostrar que la existencia de un gradiente horizontal de presión
conduce a una inclinación de la placa litosferica. Se supone que la inclinación es suficiente
débil para no perturbar el campo de velocidad en la astenosfera. Determinar la diferencia de
altura ∆h entre las 2 extremidades de la placa en función de ha , hl , v0 , L, ρ, g y µ.
El movimiento horizontal de la placa litosferica es provocada por una fuerza de tracción por
unidad de longitud F , constante.
7.1.8.2.1. Preguntas
1. Supongamos que las únicas fuerzas horizontales aplicadas a la placa en movimiento son la
fricción viscosa (calculada en la parte 7.1.8.1.2) y la tracción de la zona de subducción.
7.2. Capitulo 2
7.2.1. Problema 1
Sea la escala de velocidad u y de longitud L de un flujo turbulento. Calcular las escalas de
Kolmogorov, uk y η para los 2 movimientos turbulentos siguientes:
1. Turbulencia de baja capa atmosférica, con u = 0, 5m.s−1 y L = 100m.
3. Dar un orden de magnitud del número de Reynolds de los más pequeños vórtices.
En los 2 casos tomar ν = 15 × 10−6 m2 .s−1 .
Corrección:
1/4
Tenemos: uk = (ν × ε)1/4 y lk = η = (ν 3 /ε)
0,53
1. Calculamos la tasa de disipación: ε = 100
= 1, 25 × 10−3 m2 .s−3
lk = 1, 3mm y vk = 1, 2cm.s−1
23
2. Calculamos la tasa de disipación: ε = 0,02
= 4 × 102 m2 .s−3
lk = 0, 05mm y vk = 28cm.s−1
uk lk
3. En los 2 casos Rek = ν
≈1
7.2.2. Problema 2
Se considera un volumen cubico, lleno de aire, con una dimensión L = 3m. Una porción de
una de la caja es constituida por una superficie caliente de altura h = 0, 1m. La temperatura de
la placa asegura una diferencia de temperatura de ∆T = 10K. La difusividad térmica del fluido
es a = 2 × 10−5 m2 .s−1 .
1. Calcular el tiempo necesario a la difusión del calor en todo la pieza (es decir sobre una
distancia L), cuando tenemos:
268 CAPÍTULO 7. AYUDANTIAS
Corrección:
1. El proceso es solamente difusivo entonces tenemos V~ = ~0, la ecuación del calor se escribe:
∂T ∂T ∂T ∂T
ρCp =k =⇒ =a (7.15)
∂t ∂t ∂t ∂t
El orden de magnitud de la ecuación es:
∂T ∂T
=a (7.16)
∂t
|{z} ∂t}
| {z
∆T
t
a ∆T
2
L
sea si t = tdif f :
∆T ∆T L2 32
∼ a 2 =⇒ tdif f ∼ = = 125h (7.17)
tdif f L a 2 × 10−5
2. Supongamos que la pieza es a temperatura de 20◦ C, es decir T = 293K. Primero verificamos
que el flujo de convección natural es bien turbulento. Se calcula el número de Rayleigh:
1
gβ∆T L3 9, 81 293 10 × 33
Ra = = ≈ 9 × 109 (7.18)
ν2 10−5
El flujo es bien turbulento (ver clase de transferencia de calor).
Sea un proceso de difusión turbulenta donde el tiempo caracterı́stico es tturb . La velocidad
ahora no es nula. Notamos U esta velocidad. El tiempo tturb depende de la velocidad del
fluido U y del largo de la pieza L sea: tturb ∼ UL . Se debe estimar U . U es una velocidad de
convección natural que se desarrolla a lo largo de la placa de altura h. Se debe estimar esta
velocidad U . Hagamos un balance de fuerza sobre una partı́cula de fluido, cerca de la placa,
el movimiento de convección natural es provocado por un fuerza de Arquimedes:
sea si t = tturb :
U ∆ρ
ρ ∼ g∆ρ =⇒ U ∼ gt (7.22)
t ρ
Aquı́ t es el tiempo de calentamiento de la partı́cula frente la placa de altura h. Este tiempo
es t ∼ Uh . Tenemos entonces:
s
∆ρ h ∆ρ
U∼ g =⇒ U ∼ gh (7.23)
ρ U ρ
7.2. CAPITULO 2 269
∆ρ ∆T
∼ (7.24)
ρ T
U
Uconv ∼ ∼ 0, 02 (7.26)
10
Finalmente, el tiempo caracterı́stico tturb se estima:
L 3
tturb = ∼ = 150s (7.27)
Uconv 0, 02
7.2.3. Problema 3
Se considera un flujo turbulento que tiene una velocidad promedio caracterı́stica de Up =
30m.s−1 y una escala caracterı́stica promedio Lp = 10cm. Las fluctuaciones turbulentas es decir
los vórtices de más grande tamaño tienen una velocidad, u, del orden de 20 % de la velocidad
promedio. La viscosidad cinemática del flujo es ν = 1,5 × 10−5 m2 .s−1 . Se supone que podemos dar
un orden de magnitud de los parámetros que describen el flujo turbulento usando la cascada de
Kolmogorov.
Corrección:
1. Queremos determinar las macro-escalas L y u ası́ que las escalas caracterı́sticas de Kolmo-
gorov lk y uk .
Se sabe que la transferencia de energı́a del flujo promedio hacia los grandes vortices se
realiza en el mismo tiempo caracterı́stico:
Lp L
tp = ∼ = tmacro escala (7.29)
Up u
encontramos que:
L = 2cm (7.30)
270 CAPÍTULO 7. AYUDANTIAS
7.2.4. Problema 4
Consideramos un flujo promedio atrás un objeto de dimensión tı́pica Lp = 1m que se desplaza
a la velocidad Up = 10m.s−1 en el agua.
1. ¿Es necesario tomar en cuenta la turbulencia para modelar este flujo?
2. Se estima que la macro-escala de la turbulencia, en la zona situada atrás el objeto, es del
orden de 10 % del tamaño del objeto:
a) Estimar, en esta zona, la escala la más fina del flujo, es decir la escala de Kolmogorov
(escala donde actúa la viscosidad).
b) Estimar la disipación másica en esta zona. Comentar.
3. Queremos simular este flujo resolviendo correctamente las escalas disipativas.
¿Cuentos puntos necesita la malla para modelar el flujo en un cubo de tamaño L3 ?
¿Cuál es el espacio de memoria que se necesita para almacenar los campos a los puntos
de la malla?
¿Este calculo es posible? Si no, ¿qué solución se puede proponer?
Corrección:
1. Rep = 107 es turbulento
2. a) Tenemos una estimación del tamaño de los grande vórtices: L = 0, 1Lp =10cm. Se puede
encontrar su velocidad:
Lp L
tp = ∼ = tmacro escala =⇒ u = 1m.s−1 (7.37)
Up u
Las escalas de Kolmogorov se determine de la misma manera que los ejercicios anterio-
res. Ud. debe elegir el método que necesita el calculo del número de Reynolds basado
sobre el tamaño de los grandes vórtices ReL . No se conoce ε0 en esta pregunta.
lk = 18µm uk =?? (7.38)
7.2. CAPITULO 2 271
u3
ε0 = = 10m2 .s−3 = 10W.kg−1 (7.39)
L
3. Si realizamos una simulación directa de un flujo de fluido debemos describir todas las escalas
del flujo. En una dirección debemos tener a lo menos el numero de punto siguientes
L 3/4
Np ∼ = ReL (7.40)
lk
para una malla 3D tenemos como número de puntos:
3
L
3/4 3
Np3 ∼ = ReL (7.41)
lk
uL 1 × 0, 1
ReL = = ≈ 105 (7.42)
ν 10−6
es decir 9/4
Np3 ∼ 105 ≈ 178 × 109 puntos (7.43)
Si ReL = 107 tenemos Np ≈ 1016 puntos es decir Noctet = 8 × 1017 ≈ 1018 octets, Noctet ≈
109 Go. En realidad se necesita más memoria:
Para bien representar los pequeños vórtices debemos tener más de una malla por vórtice,
debemos tomar a lo menos alrededor de 10 mallas.
DWp
El número de Reynolds del jet es, Re = ν
, donde la viscosidad cinemática, ν = 10−6 m2 .s−1 .
Experimentalmente, sabemos, que si Re > 104 , el jet es turbulento, y en este caso:
• el comportamiento del flujo promedio y el tamaño de los grandes vórtices, Lg , son
independientes de la viscosidad y pueden ser estimado usando la siguiente formula:
3/2
Lg = k ε .
1. Con las hipótesis presentadas anteriormente, ¿podemos decir que los más pequeños vórtices
son independientes de la viscosidad? ¿Por qué?
2. Cuando Wp1 = 1m.s−1 , al punto P situado sobre el eje (del conducto) a 12m de la descarga,
se mide una velocidad promedio de los grandes vórtices Ug1 = 0, 6m.s−1 y una energı́a ciné-
tica turbulenta k = 0, 02m2 .s−2 . Además, se estima que el tamaño de los grandes vórtices es
aproximadamente de Lg1 = 10cm.
3. De la misma manera que en los ejercicios anteriores podemos calcular lk . Ud. puede elegir el
método que necesita el calculo de la tasa de disipación o el método que necesita el calculo
del número de Reynolds basado sobre el tamaño de los grandes vórtices ReL .
lk1 ≈ 75µm
lk2 ≈ 46µm
7.2.6. Problema 6
Una estructura coherente puede ser modelada con un vórtice rectilı́neo con un núcleo viscoso
de radio R = 10cm. Este vórtice se desplaza a la velocidad constante U = 1cm.s−1 . El flujo de
fluido es agua, la viscosidad cinemática es ν = 10−6 m2 .s−1 . Consideramos que la coherencia de la
estructura es mantenida cuando la difusión viscosa no actúa sobre una distancia superior a 10 %
del núcleo del vortice, es decir que la estructura puede ser considera todavı́a como coherente.
2. Consideramos ahora un flujo de fluido turbulento. Las fluctuaciones de velocidad tienen una
escala u = 0, 1cm.s−1 y una escala de tamaño l = 1mm.
Corrección:
Durante este tiempo, la estructura se movió (traslado por el fluido, transporte por convec-
ción) de una distancia:
L = tdif f × U = 100 × 1 = 1m (7.50)
7.2.7. Problema 7
Sea un flujo turbulento donde el campo de velocidad promedio se reduce a la única componente
U (y), con 0 ≤ y ≤ h. Se supone que para y → 0, este campo se escribe:
U (y) = U0 + ay 2 (7.54)
7.3. Capitulo 3
7.3.1. Problema 1: Hidróptero
Queremos propulsar a la velocidad de 72 [km/hr] un bote de 10 [ton] con un ala submarina
donde el coeficiente de sustentación es CL = 0, 7, y el coeficiente de arrastre es CD = 0, 06.
Calcular:
1. Determinar la fuerza de fricción que se ejerce sobre una de las caras de la placa.
7.3.3. Problema 3
Una avión tiene una superficie de ala de 80m2 , un coeficiente de sustentación CL = 0, 80 y
un coeficiente de arrastre de Cd = 0, 04. El avión vuela a una velocidad de 360km.h−1 . La masa
especifica del aire es ρ = 1, 2kg.m3 y su viscosidad cinemática es ν = 1, 5 × 10−5 m2 .s−1 .
2. Si consideramos que 50 % del arrastre es producido por la estructura del avión (carlinga,
cabina, aparato de dirección). Calcular la potencia necesaria para la propulsión del avión.
3. Si el avión vuela en caı́da ahora a 400km.h−1 con una incidencia de −4◦ (respecto al hori-
zontal). Calcular la fuerza de sustentación.
4. Conclusión.
7.3. CAPITULO 3 275
∂u ∂u 1 ∂p ∂ 2u
u + v = − + ν (7.55)
∂x2 ∂y ρ ∂x ∂y 2
∂u ∂v
+ = 0 (7.56)
∂x ∂y
2
Con p + ρ u2e = cte para y = δ.
2. Espesor de desplazamiento.
4. Factor de forma.
276 CAPÍTULO 7. AYUDANTIAS
6. Expresión del esfuerzo de corte en la pared, a partir del espesor de cantidad de movimiento.
10. Fuerza de arrastre, desde el borde de ataque hasta una posición longitudinal x.
7.3.7. Problema 7
Una corriente uniforme U∞ de un fluido de densidad ρ y viscosidad µ que incide sobre una
esfera de radio R. Supongamos que número de Reynolds es grande, Re 1, y que se desarrolla
una capa limite laminar alrededor de la esfera.
−U0 L0
ue (x) = (7.58)
x
donde U0 y L0 son dos parámetros homogéneos a una velocidad y a una longitud respectivamente.
Notamos c0 = U0 L0 .
1. Por qué, tenemos un gradiente longitudinal de presión al largo de la placa? Determinar una
expresión de este gradiente dPdx(x) . Este gradiente es favorable o desfavorable?
2. ¿ Hay en este caso una solución de afinidad para la capa limite sobre este obstáculo?
7.3.9. Problema 9
Sea el perfil de velocidad sugerido por Kármán-Polhausen en 1921 para modelar una capa
limite: 3 4
u y y y
≈2 −2 + (7.62)
Ue δ δ δ
1. Calcular las expresiones de xθ , xδ , Cf y H.
2. ¿Satisface el perfil de la Ec. 7.63 la condición de frontera de flujo laminar sobre la placa
plana?
278 CAPÍTULO 7. AYUDANTIAS
7.3.10. Problema 10
Sea el perfil de velocidad sugerido por Kármán-Polhausen en 1921 para modelar una capa
limite:
u y y2
≈2 − 2 (7.63)
Ue δ δ
2. Sobre una placa plana afilada fluye aire a 20◦ C y 105 Pa con una velocidad de 2m.s−1 .
Calcular:
a) la velocidad u
b) y el esfuerzo de cortadura local τ
2
u y m y
= l + P ara 0 ≤ y ≤ δ (7.64)
ue θ 2 θ
u
= 1 P ara y > δ (7.65)
ue
7.3. CAPITULO 3 279
Calculando θ, mostrar que l y m son necesariamente relacionadas por una función l = l(m),
que no explicitaremos.
Figura 7.7: Resultados experimentales de Fransson & Alfredsson (2003), que muestra la transición
progresiva de la capa lı́mite de Blasius
q a una capa lı́mite aspirada. El eje de ordenadas es la variable
y
adimensional ζ = δ(x) con δ(x) = Uνxe . La aspiración empieza a una distancia xa = 360mm del
borde de ataque.
Figura 7.8: Esquema de la capa lı́mite que se desarrolla sobre una pared impermeable (arriba) y
sobre una pared porosa (abajo)
Propiedades de la pared
La pared puede ser totalmente impermeable, o al contrario, permeable a partir de una distancia
xa (distancia tomada desde el borde de ataque). La parte porosa es conectada a un sistema de
aspiración. Los poros del material son, en primera aproximación, cilindros microscópicos de eje de
revolución paralelo al eje Oy. Ası́, la aspiración se realiza en la dirección normal a la pared.
En los 2 casos (con y sin aspiración) tenemos la condición de adherencia para la velocidad
tangencial:
∀x > 0, u(x, 0) = 0 (7.69)
2. Si el fluido es ideal (sin viscosidad), el caudal a través de la superficie serı́a, qi , que Ud. debe
calcular. ¿Qué se puede decir el signo de q − qi ? Comentar.
τp es el esfuerzo de corte tangencial que ejerce el fluido sobre la pared, a una distancia
x del borde de ataque, ρ es la masa especifica del fluido.
4. Calcular numéricamente este espesor en el caso de una capa lı́mite de aire a temperatura
ambiente (ν = 1, 5 × 10−5 m2 .s−1 ), con un flujo exterior de velocidad Ue = 10m.s−1 y una
Ue
velocidad de aspiración V = 1000 .
5. Expresar analı́ticamente el coeficiente de fricción parietal local de esta capa lı́mite. Comentar
la expresión encontrada.
3. En el caso de una capa lı́mite de aire a temperatura ambiente, con un flujo exterior de
velocidad Ue = 10m.s−1 . Dar el valor de δ ∗ a una distancia del borde de ataque x = 1m.
Comentar y comparar con el valor encontrado anteriormente (Parte 2).
4. Calcular el coeficiente de fricción parietal local de esta capa lı́mite. Dar una expresión ana-
lı́tica donde aparece una derivada de la función f . Calcular numéricamente esta derivada.
6. Comentario final de los resultados de las partes 2 y 3: Comentar los fenómenos que podrı́an
impedir la realización experimentales de estos flujos. En particular, la transición progresiva
de la capa lı́mite de Blasius a una capa lı́mite asintótica.
Para tener una buena calidad en las experiencias que serán realizadas, se desea establecer un flujo
lo más uniforme posible sobre toda la longitud del túnel de viento (Ver Figura 7.9).
Las mediciones de velocidad permitieron constatar una leve aceleración del flujo en la
sección de salida del túnel de viento respecto a la velocidad de entrada. La velocidad a la
entrada es Ue (0) y a la salida es Ue (L).
Figura 7.9: Esquema del flujo uniforme entre las capas lı́mites que se desarrollan sobre las paredes
del túnel de viento. La velocidad a la entrada es Ue (0) y a la salida es Ue (L). El espesor δ de las
capas lı́mites representadas sobre el dibujo es bastante superior al espesor observado en la realidad.
2. a) Explicar por qué el hecho de dar una leve inclinación a las paredes del túnel de viento
(Ver Figura 7.10) permite conservar una velocidad Ue constante sobre todo el largo del
túnel de viento.
b) Notamos la altura de la sección de salida h0 . Expresar h0 en función de α y L
Figura 7.10: Esquema del túnel de viento con sus paredes inclinadas. La sección de salida tiene
ahora una altura h0 . La inclinación representada sobre el dibujo es bien superior a la realidad.
3. Consideramos que el desarrollo de las capas lı́mites sobre las paredes del túnel de viento
esté bien representado por las formulas de Blasius. Es decir que se puede considerar una
capa lı́mite laminar sobre una placa plana en ausencia de gradiente de presión exterior
(Ue (x) = Cte para ∀x). Tenemos entonces:
284 CAPÍTULO 7. AYUDANTIAS
El espesor de desplazamiento:
1, 72x
δ∗ = √ (7.73)
Rex
El espesor de déficit cantidad de movimiento:
x
θ = 0, 664 √ (7.74)
Rex
donde se notó Rex = uνe x el número de Reynolds basado sobre el largo caracterı́stico x y ν
la viscosidad cinemática del aire, ν = 1, 5 × 10−5 m2 .s−1 .
a) Verificar que la capa lı́mite que se desarrolla sobre las paredes del túnel de viento es
esencialmente turbulenta.
1,768
Comentario: Ud. podrı́a usar el gráfico de la función f (z) = z 0,04
+ 39603,2
z 1,04
represen-
tada sobre la Figure 7.11.
1,768 39603,2
Figura 7.11: Gráfico de la función f (z) = z 0,04
+ z 1,04
.
7.4. CAPITULO 4 285
2. Se supone que la evolución del espesor de las capas lı́mites turbulentas que se desarrollan
sobre las paredes del túnel de viento está bien representada por la formula:
0, 37x
δ(x) = (7.76)
Rex 1/5
Si suponemos que el perfil de velocidad en una capa lı́mite turbulenta es de la forma:
1/7
u(x, y) y
= (7.77)
Ue δ
donde y es la distancia a la pared sobre la cual se desarrolla la capa lı́mite.
dθ 1 dUe Cf
+ (δ ∗ + 2θ) = (7.78)
dx Ue dx 2
7.4. Capitulo 4
7.4.1. Problema 1: Modelo simplificado de turbulencia cerca de la pa-
red. Ley de pared
Consideramos un flujo a gran numero de Reynolds paralelamente a una pared plana definida
en y = 0 en un sistema de referencia cartesiano. Estudiamos la capa limite que se desarrolla cerca
de esta pared y a la influencia de la turbulencia sobre esta.
Notamos el vector velocidad promedio V̄ = ū~ex + v̄~ey + w̄~ez y las fluctuaciones de velocidad
u0 , v 0 , w0 según x, y, z respectivamente.
V̄ = ū(y)~ex (7.79)
y que todas las cantidades promedio depende solamente en primera aproximación que de la
coordenada y. Se desprecia los efectos de la gravedad y se considera la presión como presión
motriz: pm = p + ρgy.
1. Escribir las ecuaciones de Reynolds. Mostrar que, en la dirección del flujo tenemos un equili-
brio entre los efectos viscosos y los efectos de turbulencia. Mostrar que cuando nos alejamos
de la pared, la presión a tendencia a disminuir porqué tenemos fluctuaciones turbulencia.
286 CAPÍTULO 7. AYUDANTIAS
2. Integrar una vez (respecto a y) la ecuación de Reynolds en la dirección del flujo. Mostrar
que la constante de integración que aparece puede se expresar simplemente en función del
esfuerzo parietal promedio:
3. a) Construir por análisis adimensional una longitud lp y una velocidad uτ dicha velocidad
de fricción, a partir de ρ, τp , ν viscosidad cinemática del fluido. verificar que lp es
relacionado a uτ por una relación simple.
b) Adimensionar la ecuación diferencial que define Vx (y). Por eso se realizará el cambio de
variable:
y ū(y)
y+ = u+ = = u+ (y + ) (7.81)
lp uτ
4. Justificar fı́sicamente la existencia de una subcapa viscosa cerca (muy) de la pared en la cuál
el esfuerzo de Reynolds es despreciable, y el perfil de velocidad es lineal:
u+ = y + (7.83)
lm = κy (7.84)
u0 v 0
− (7.86)
k
7.5. Capitulo 5
7.5.1. Problema 1: Cuerpo semi-infinito de Rankine
Se considera un flujo potencial compuesto de un flujo uniforme (según el eje x) más de una
Q
fuente de intensidad m = 2πb (b es la longitud sobre z tomada, aquı́ tomamos b = 1) posicionada
en el origen.
Figura 7.13: Esquema del invernadero en el flujo de aire en el sistema de referencia (x, y) y (r, θ)
1. Mostrar que la función φ es la suma de dos flujos elementales φ1 y φ2 . Dar las expresiones
de φ1 y φ2 .
3. Dar una expresión de la presión en todo punto del espacio: P (r, θ). Ud. dará las hipótesis
que permiten aplicar el teorema.
1. El aire es un fluido viscoso newtoniano. Que podemos deducir respecto a la velocidad del
aire en contacto a la pared del cilindro? Determinar el valor de la circulación del vector
velocidad, Γ, sobre el contorno del cilindro.
Figura 7.14: Esquema del invernadero en el flujo de aire en el sistema de referencia (x, y) y (r, θ)
Preguntas conceptuales
Responder en algunas lineas:
1. ¿Por qué se puede considerar que el flujo exterior es un flujo potencial?
2. Se supone que el flujo es bidimensional, horizontal, en un plano xOy. ¿ Para qué condiciones?
7.5. CAPITULO 5 291
Preguntas
1. Escribir el potencial complejo f del flujo como la superposición del vórtice-tornado y del
vórtice virtual.
3. Mostrar que si la frontera del disco (la torre) es una linea de corriente, es decir:
4. Deducir de la pregunta anterior el valor de x2 función de R y x1 que se debe elegir para que
el flujo no penetre en el disco (condición de impermeabilidad).
a) Calcular su circulación.
b) Estimar numéricamente (Pint − Pext ) en z = R cuando R = 10m, x1 = 3R y ρ =
1, 2kg.m−3 . Calcular la fuerza ejercida sobre 1m2 de torre. Explicar fı́sicamente qué va
a pasar.
Preguntas
2. a) Expresar el potencial de las velocidades φ(r, θ) y la función de corriente Ψ(r, θ) del flujo
en función de las coordenadas cilı́ndricas (r, θ).
b) Expresar el potencial de las velocidades φ(x, y) y la función de corriente Ψ(x, y) del
flujo en función de las coordenadas cartesianas (x, y).
3. A partir de la definición de φ(r, θ), expresar las componentes radial y tangencial, Vr y Vθ del
vector velocidad en un punto (r, θ). Deducir el modulo q de la velocidad en este punto.
4. Dar las componentes cartesianas del vector velocidad, Vx y Vy , en un punto (x, y) a partir
de las expresiones de Vr y Vθ .
7. Encontrar el resultado anterior (es decir, el obtenido en la pregunta 6) a partir del conoci-
miento del campo de velocidad en coordenadas cilı́ndricas.
7.5.7. Problema 7
En 1927, un hombre llamado Flettner hizo construir un barco con dos cilindros giratorios para
que actuaran como velas (Ver Figura 7.16). Si la altura de los cilindros es de 15m y el diámetro
es 2, 75m,
1. Encontrar el empuje propulsor máximo posible causado por los cilindros sobre el barco.
Datos
Los cilindros rotan a una velocidad de 750rpm mediante la acción de una máquina de vapor
localizada por debajo de la cubierta.
7.5.8. Problema 8
Se superpone, sobre el mismo eje origen, una fuente lineal de caudal Q por unidad de longitud
y un vórtice de circulación Γ, se obtiene un flujo bidimensional encontrado en las turbo-maquinas.
1. Determinar las ecuaciones del potencial complejo, del potencial de las velocidades y de la
función de corriente de este flujo.
4. Mostrar que las lineas de corriente de este flujo son espirales logarı́tmicas.
x = 3 cm y y = 0, 2 cm
7.5.10. Problema 10
Un viento de velocidad uniforme al infinito V sopla horizontalmente en la dirección Ox sobre
una chimenea vertical de sección circular, de eje perpendicular a Ox y de diámetro 2a = 4m (Ver
Figura 7.17).
1. Determinar:
7.6. Capitulo 6
7.6.1. Problema 1: Tobera de Laval
Un flujo de aire atraviesa el conducto de la Figura 7.18 donde A1 = 24 cm2 , A2 = 18 cm2 y
A3 = 32 cm2 . En la sección 1 tenemos M a1 = 2, 5, P1 = 40 kPa y T1 = 30◦ C. En la sección 2
existe una onda de choque normal. Calcule:
1. El gasto másico.
Se supone por este problema, una entrada en la atmósfera a la altura z = 20 km. A esta
altura la temperatura del aire es de 217 K.
1. Calcular la temperatura de estancamiento por los 2 valores del numero de Mach de la nave,
M1 = 3 y M1 = 5.
2. Se supone que la onda de choque que se forma frente la nave es una onda de choque plana.
Determinar la temperatura T2 del aire entre la onda y la pared de la nave.
296 CAPÍTULO 7. AYUDANTIAS
Preguntas
1. Con las hipótesis ya enunciadas, mostrar que se puede escribir en cada sección:
p
c2 = k (7.96)
ρ
donde c es la velocidad del sonido, p la presión y ρ la masa especifica en una sección.
c2
h= (7.98)
k−1
b) Si en una sección de la tobera tenemos una velocidad v, una presión p y una masa
especifica ρ. Mostrar que se puede escribir la relación siguiente:
v2 p k p0 k
+ = (7.99)
2 ρk−1 ρ0 k − 1
Se pide escribir las formulas utilizadas y detallar los cálculos (No utilizar tablas).
4. Entre una sección x y una sección x + dx, tenemos las relaciones siguientes:
dS 2 dv
= M −1 (7.100)
S v
!−1
dv k−1 2 dM
= 1+ M (7.101)
v 2 M
a) Dar una expresión entre las superficies S1 y S2 y los números de Mach M1 y M2
correspondiente a dos secciones de abscisa x1 y x2 .
b) Calcular la superficie de la sección al cuello, si tenemos un número de Mach, M = 2,
en una sección de superficie S = 1000cm2 .
Preguntas
5. Para cual valor de la presión estática pa una onda de choque derecha se posicionará a la
salida de la tobera.
6. Ahora tenemos una presión estática en el estanque igual a pa = 1, 84bar. Tenemos una onda
de choque en el divergente. Calcular el diámetro de la sección donde se produce el choque.
S2
2. Calcular la razón de las secciones: S1
.
Condiciones estándares
Figura 7.22: Flujo en una tobera. Tenemos una onda de choque a la salida de la tobera.
Appendices
301
Apéndice A
Operadores diferenciales
303
304 APÉNDICE A. OPERADORES DIFERENCIALES
Apéndice B
305
306 APÉNDICE B. OPERADORES EN VARIOS SISTEMAS DE COORDENADAS
Ecuación de Navier-Stokes
volumen (son aquellas que afectan a todo el volumen de la parcela de fluido, son de acción a dis-
tancia.). Generalmente, en mecánica de fluidos consideramos unicamente las fuerzas gravitacional,
es decir, f~m = ~g .
∂u ∂v ∂w
+ + =0 (C.1)
∂x ∂y ∂z
Conservación de la cantidad de movimiento
sobre eje x
" # " #
∂u ∂u ∂u ∂u ∂p ∂ 2u ∂ 2u ∂ 2u
ρ +u +v +w =− + fx + µ + + (C.2)
∂t ∂x ∂y ∂z ∂x ∂x2 ∂y 2 ∂z 2
sobre eje y
" # " #
∂v ∂v ∂v ∂v ∂p ∂ 2v ∂ 2v ∂ 2v
ρ +u +v +w =− + fy + µ + + (C.3)
∂t ∂x ∂y ∂z ∂y ∂x2 ∂y 2 ∂z 2
309
310 APÉNDICE C. ECUACIÓN DE NAVIER-STOKES
sobre eje z
" # " #
∂w ∂w ∂w ∂w ∂p ∂ 2w ∂ 2w ∂ 2w
ρ +u +v +w = − + fz + µ + + (C.4)
∂t ∂x ∂y ∂z ∂z ∂x2 ∂y 2 ∂z 2
El tensor de esfuerzo
Conservación de la masa
sobre eje r
C.4. COORDENADAS ESFÉRICAS 311
" #
∂ur ∂ur uθ ∂ur ∂ur u2θ ∂p
ρ + ur + + uz − = − + fr
∂t ∂r r ∂θ ∂z r ∂r
" #
∂ 2 ur 1 ∂ur ur 1 ∂ 2 ur ∂ 2 ur 2 ∂uθ
+ µ + − 2 + 2 + − 2
∂r2 r ∂r r r ∂θ2 ∂z 2 r ∂θ
(C.6)
sobre eje θ
" #
∂uθ ∂uθ ur uθ uθ ∂uθ ∂uθ 1 ∂p
ρ + ur + + + uz = − + fθ
∂t ∂r r r ∂θ ∂z r ∂θ
" #
∂ 2 uθ 1 ∂uθ uθ 1 ∂ 2 uθ ∂ 2 uθ 2 ∂ur
+ µ + − 2 + 2 + + 2
∂r2 r ∂r r r ∂θ2 ∂z 2 r ∂θ
(C.7)
sobre eje z
" #
∂uz ∂uz uθ ∂uz ∂uz ∂p
ρ + ur + + uz = − + fz
∂t ∂r r ∂θ ∂z ∂z
" #
∂ 2 uz 1 ∂uz 1 ∂ 2 uz ∂ 2 uz
+ µ + + 2 +
∂r2 r ∂r r ∂θ2 ∂z 2
(C.8)
El tensor de esfuerzo
sobre eje r
sobre eje θ
sobre eje ϕ
Apéndice D
∂φ ∂φ
ux = ; uy = (D.1)
∂x ∂y
∂ψ ∂ψ
ux = ; uy = − (D.2)
∂y ∂x
∂φ 1 ∂φ
ur = ; uθ = (D.3)
∂r r ∂θ
1 ∂ψ ∂ψ
ur = ; uθ = − (D.4)
r ∂θ ∂r
∂φ 1 ∂φ ∂φ
ur = ; uθ = ; ux = (D.5)
∂r r ∂θ ∂x
Coordenadas esféricas (r,θ,ϕ)
∂φ 1 ∂φ 1 ∂φ
ur = ; uθ = ; uϕ = (D.6)
∂r r sin ϕ ∂θ r ∂ϕ
1 ∂ψ 1 ∂ψ
ur = − ; ux = (D.7)
r ∂x r ∂r
Coordenadas esféricas (r,θ,ϕ)
313
314APÉNDICE D. RELACIÓN ENTRE VELOCIDAD, POTENCIAL Y FUNCIÓN DE CORRIENTE
1 ∂ψ −1 ∂ψ
ur = ; uϕ = (D.8)
r2 sin ϕ ∂ϕ r sin ϕ ∂r