Tema 5. Distribuciones de Velocidad 2023-24

FENÓMENOS DE TRANSPORTE
AVANZADOS.
Tema 5. Distribuciones de velocidad
Fenómenos de transporte avanzados
BLOQUE A. PRINCIPIOS MATEMÁTICOS Y BALANCES

Tema 1. Introducción a los fenómenos de transporte
Tema 2. Balance microscópico de materia
Tema 3. Balance microscópico de cantidad de movimiento
Tema 4. Balance microscópico de energía
BLOQUE B. DISTRIBUCIONES DE PROPIEDAD Y ANALOGÍAS

Tema 5. Distribuciones de velocidad
Tema 6. Distribuciones de temperatura
Tema 7. Distribuciones de concentración
Tema 8. Transporte de interfase y analogías
Indice
• 1. Introducción
• 2. Aplicación de las ecuaciones de variación a la obtención de perfiles de

presión y velocidad
• 2.1. Flujo en una película descendente en un plano inclinado
• 2.2. Flujo de Coutte totalmente desarrollado
• 2.3. Flujo interno en un tubo circular
• 2.4. Flujo descendente a lo largo de la pared exterior de un tubo
vertical
• 2.5. Forma de la superficie de un líquido que gira
Bibliografía
BIBLIOGRAFÍA FUNDAMENTAL
• Bird, R.B.; Stewart, W.E.; Lightfoot, E.N. “Fenómenos de Transporte”.

2ª Ed. Limusa Wiley (2006).
BIBLIOGRAFÍA COMPLEMENTARIA
• Costa, E.; et al. “Ingeniería Química, volumen II: Fenómenos de

Transporte”. Ed. Alhambra (1984).
• Welty, J.R.; Wicks, C.E.; Wilson, R.E. & Rorrer, G.L. “Fundamentals of
Momentum, Heat, and Mass Transfer”. 5ª Ed. John Wiley & Sons
(2008)
1. Introducción
• La obtención de los perfiles de velocidad y/o presiones para flujo laminar en un
sistema isotérmico y de un único componente (fluido puro) puede realizarse mediante
la resolución de la ecuaciones de variación:
• Ecuación de continuidad
• Ecuación de movimiento
Procedimiento:
1. Se simplifican las ecuaciones de variación para el caso bajo estudio.

2. Se sustituye la expresión adecuada de la densidad de flujo de cantidad de movimiento
(ley de Newton para fluidos newtonianos) y se obtiene la ecuación diferencial que
describe la distribución de velocidad.
3. Se integran las ecuaciones diferenciales obtenidas previamente, considerando las
adecuadas condiciones límite (o de frontera). De esta forma se obtienen,
explícitamente las expresiones de la distribución de presiones, de la densidad de flujo
de cantidad de movimiento y de velocidad en el sistema.
4. A partir de la información anterior se procede al cálculo de las magnitudes de interés:
velocidad media y máxima, caudal, fuerzas que actúan sobre las superficies límite.
1. Introducción
• La determinación de los valores de las constantes de integración se realiza a partir de
la formulación de las correspondientes CONDICIONES LÍMITE:
CONDICIÓN LÍMITE DE NO-DESLIZAMIENTO:
Para un fluido en contacto con una pared sólida, la velocidad del fluido debe ser igual a la
de la pared:
En otras palabras, no hay deslizamiento entre el fluido y la pared. Las partículas de fluido
adyacentes a la pared se adhieren a la superficie de ésta y se mueven a la misma
velocidad que la pared.
En el caso especial de una pared en reposo, el fluido

adyacente a ésta tiene también velocidad cero.
Para los casos en los que también se consideren los

efectos de temperatura, la temperatura del fluido debe
ser igual a la de la pared; es decir, Tfluido = Tpared.
1. Introducción
CONDICIÓN LÍMITE DE INTERFASE:
Cuando dos fluidos A y B se encuentran en una interfase, las condiciones frontera de la

interfase son:
donde además de que las velocidades de los dos fluidos han de ser iguales, los esfuerzos
cortantes que actúan sobre una partícula de fluido adyacente a la interfase, en la
dirección paralela a la interfase, también han de ser iguales entre los dos fluidos.
Esta condición se aplica en las interfases

entre dos líquidos, o entre un líquido y un
gas cuya interfase tiene una curvatura
pronunciada, como en los meniscos de
los líquidos que se elevan en un tubo
capilar.
Además, si se consideran los efectos de

temperatura, TA = TB en la interfase.
1. Introducción
Cuando dos fluidos A y B se encuentran en una interfase, las condiciones frontera de la

interfase son:
Ejemplo: agua líquida en contacto con aire. La interfase es plana y el agua se desplaza
horizontalmente.
Las velocidades del aire y del agua deben de ser
iguales en la interfase y el esfuerzo cortante que
actúa sobre una partícula de agua en la interfase
debe ser igual al que actúa sobre una partícula de
aire justo en la interfase.
Como magua es casi 50 veces mayor que maire, para

que los esfuerzos cortantes sean iguales:
v x v x

y agua
y
aire
Por tanto, es razonable considerar el esfuerzo cortante que actúa en la superficie del agua
como despreciable (nulo) frente a los esfuerzos cortantes en cualquier otra parte del agua.
1. Introducción
Otra manera de decir esto es que el agua en movimiento arrastra al aire junto a ella con
poca resistencia del aire; o bien, el aire no frena al agua de manera significativa. En
resumen:
Para el caso de un líquido en contacto con un gas, y con efectos de tensión superficial
despreciables, las condiciones frontera de superficie son:
Interfase líquido  gas Plíquido  Pgas  
ij líquido 0
2. Aplicación de las ecuaciones de variación a la obtención de perfiles
2.1. Flujo en una película descendente en un plano inclinado
• Aplicaremos las ecuaciones de variación para el caso del flujo de un líquido que
desciende por una lámina plana inclinada de longitud L y ancho W (flujo viscoso
isotérmico bajo la influencia de la gravedad).
• Este estudio tiene interés para los siguientes casos:

• Torres de pared mojada.
• Experimentos de evaporación y absorción de gases.
• Aplicaciones de recubrimientos.
Actividad propuesta: determine la distribución de presión y de velocidad, la densidad de flujo de

cantidad de movimiento, la velocidad máxima y media, el caudal volumétrico, el espesor de la película y
la fuerza que ejerce el fluido sobre la superficie sólida.
HIPÓTESIS DE PARTIDA:
• Régimen Laminar.
• Fluido puro.
• Flujo estacionario e isotérmico.
• Densidad y viscosidad constantes.
• Se considerará que L y W son grandes en comparación con el espesor de la
película δ, para poder despreciar las perturbaciones de entrada y salida
• La anchura de la placa W se considera suficientemente grande como para
considerar flujo bidimensional en las direcciones “x” y “z”.
x=0
dx
gz x=δ
gx
gx = g cos α
gz = g sen α
z
α
x
ECUACIONES DE VARIACIÓN: x=0
dx
Ecuación de continuidad (ρ = cte ): gz x=δ
gx
gx = g cos α
z gz = g sen α
α
x
Ecuación de movimiento
( ρ = cte, μ = cte y fluido newtoniano. Ecuación de Navier-Stokes)

Dv  
r  P + m  2v + r g
Dt
æ v v v x v x ö æ 2v 2v x 2v x ö
p
rç x + vx x + vy + vz ÷   + m çç 2 + 2 + 2 ÷÷ + r g x
x
è t x y z ø x è x y z ø
æ v v v y v y ö æ 2v 2v y 2v y ö
p
÷   +mç ÷+ rg
r çç y + v x y + vy + vz ÷
y
+ +
è t x y z ø y ç x 2
y 2
z 2 ÷ y
è ø
æ v v z v z v z ö æ 2v 2v z 2v z ö

p
rç z + vx + vy + vz ÷   + m çç 2 + 2 + 2 ÷÷ + r gz
z
è t x y z ø z è x y z ø
SIMPLIFICACIONES: x=0
dx
• Movimiento unidireccional (eje z): x=δ
gz
vx = 0 vy = 0
gx
gx = g cos α
z gz = g sen α
α
• Régimen estacionario: x
• Campo de velocidad bidimensional:
• La componente “y” de la aceleración de la gravedad es nula: gy = 0
• La presión permanece constante para cada valor de “z”, es P

decir, no existe gradiente de presión aplicado que empuje al 0
z
fluido en la dirección z; el flujo es debido únicamente a la acción
de la gravedad.
x=0
dx
ECUACIONES DE VARIACIÓN SIMPLIFICADAS:
gz x=δ
Ecuación de continuidad (ρ = cte ): gx

gx = g cos α
z gz = g sen α
v z x
α
0
z
vx = 0 vy = 0 vz = vz(x)
Ecuación de movimiento:
Componente X: vx = 0
æ v v x v x v x ö æ 2v 2v x 2v x ö
p p
rç x + vx + vy + vz ÷   + m çç 2 + 2 + 2 ÷÷ + r g x
x
 r gx
è t x y z ø x è x y z ø x
estacionario vx = 0 vx = 0 vx = 0 vx = 0 vx = 0
Componente Y:
æ v v y v y v y ö æ 2v 2v y 2v y ö p
p
÷   +mç ÷ + rg  0
r çç y
+ vx + vy + vz ÷
y
+ +
ç x ÷ y y
è t x y z ø y 2 2 2
è y z ø
estacionario vy = 0 vy = 0 vy = 0 vy = 0 vy = 0 vy = 0 gy = 0
p = p(x)
2.1. Flujo en una película descendente en un plano inclinado x=0
dx
gz x=δ
Ecuación de continuidad (ρ = cte ):

gx
v z vz = vz(x) gx = g cos α
0 gz = g sen α
z z
α
x
Componente X: p
 r gx
x
Componente Y: p p = p(x)
 0
y
Componente Z: p = cte para cada z continuidad
æ v v z v z v z ö æ 2v 2v z 2v z ö æ  2v z ö

p
z 0  mç 2
ç x ÷÷ + r g z
è t x y z ø z è x y z ø è ø
estacionario vx = 0 vy = 0 continuidad f. bidimensional
x=0
dx
gz x=δ
Ecuación de continuidad:
v z gx
0 vz = vz(x) gx = g cos α
z z gz = g sen α
α
x
Componente X: Componente Y: Componente Z:
p p æ  2v z ö
 r gx  0 p = p(x) 0  m ç 2 ÷ + r gz
x y ç x ÷
è ø
Distribución de presiones
• A partir de la componente “x”, de la ecuación de movimiento simplificada, separando
variables e integrando:
dp
dx
 r gx dp  r g x dx  dp   r g dx  dp   r g cos dx
x p  r g x cos  + C
• Teniendo en cuenta la condición de contorno: x = 0, p = patm se obtiene:
p  r g x cos  + patm
ecuación que permite conocer el perfil de presión que se establece en la dirección “x”.
x=0
dx
gz x=δ
v z gx
0 vz = vz(x) gx = g cos α
z z gz = g sen α
α
x
Componente X: Componente Y: Componente Z:
p p æ  2v z ö
 r gx  0 p = p(x) 0  m ç 2 ÷ + r gz
x y ç x ÷
è ø
Distribución de velocidades
• A partir de la componente “z” de la ecuación de movimiento simplificada, separando
variables:
æ d 2v ö æ d 2v ö æ d 2v ö
0  mç z ÷+ rg mç z ÷  r g  rg sen ç z ÷   r g sen
ç dx 2 ÷ z ç dx 2 ÷ z ç dx 2 ÷ m
è ø è ø è ø
• e integrando:
dv z r g sen rg sen 2
 x + C1 v z   x + C1x + C2
dx m 2m
x=0
dx
dv z r g sen rg sen 2 gz x=δ

 x + C1 v z   x + C1x + C2
dx m 2m
gx
Evaluación de la constante C1 gx = g cos α
z gz = g sen α
• Condición de contorno: x
α
Interfase líquido  gas Plíquido  Pgas  

ij líquido 0
La condición de contorno en este caso es que el esfuerzo en la superficie de la fase

líquida es extremadamente pequeño frente a los esfuerzos en cualquier otra parte del
líquido:
dv z
x 0 0
dx
dv z r g sen r g sen
 x + C1 0 0 + C1 C1  0
dx m m
x=0
dx
dv z r g sen r g sen 2 gz x=δ

 x + C1 vz   x + C1x + C2
dx m 2m
gx
gx = g cos α
Evaluación de la constante C2 gz = g sen α
z
α
• Condición de contorno: x
No hay deslizamiento entre el fluido y la pared. En el caso de una pared en reposo, el

fluido adyacente a ésta tiene también velocidad cero.
r g sen 2 r g sen 2 r g sen 2
vz   x + C2 x  vz  0 0  + C2 C2  
2m 2m 2m
• Por consiguiente, la ecuación que permite conocer el perfil de velocidades es:
r g sen 2 r g sen 2 r g sen 2

vz  
2m
x +
2m
 
2m
  x2   r g 2sen æ æ x ö ö
2
vz  ç1  ç ÷ ÷
2m ç èø ÷
è ø
x=0
Densidad de flujo de cantidad de movimiento
dx
• La única componente no nula de la velocidad x=δ
gz
es vz y sólo se produce variación de la
misma en la dirección x, por consiguiente: gx
gx = g cos α
z gz = g sen α
α
x
A partir del perfil de velocidad:
r g 2sen æ æ x ö ö
2
dv z r g sen
vz  ç1  ç ÷ ÷  x
2m ç èø ÷ dx m
è ø
dv z
 xz  m  rg x sen
dx
• La densidad de flujo de cantidad de movimiento es nula en la superficie libre (x = 0) y

va aumentando a medida que nos acercamos a la superficie sólida, donde alcanza su
valor máximo.
dv z
dx
Velocidad máxima x=0
dx
• Vz debe satisfacer la condición de máximo dvz/dx:
gz x=δ
dv z r g sen
 x
dx m gx
gx = g cos α
dv z r g sen z gz = g sen α
 x 0 x 0 x
α
dx m
Es decir, el valor máximo de la componente z de la velocidad se alcanza en x = 0

Comprobemos que se trata de un máximo:
æ d 2v ö
ç z ÷   rg sen  0
ç dx 2 ÷ m
è ø
• De acuerdo con la distribución de velocidad obtenida:
r g 2sen æ æ x ö ö
2
r g 2sen
vz  ç1 ç ÷ ÷ vz  v z max 
2m ç èø ÷ x 0
2m
è ø
La velocidad varía entre este valor máximo en x = 0 y vz = 0 para x = 
x=0
Velocidad media
dx
• Esta velocidad media corresponde con la
observada o determinada a nivel gz x=δ
macroscópico.
gx
gx = g cos α
• Esta velocidad puede determinarse, cuando gz = g sen α
z
se dispone del perfil de velocidades, α
x
calculando la media referida a la sección
transversal de la película
vz 
1
S
v z
dS
S
• La sección transversal de la película “S” viene dada por el producto de la anchura (en
la dirección y) y de la altura (en la dirección x).
  
 W v dx z W  v z dx  v dx z
vz  0
 0
 0
W W 
dx
Velocidad media
gz x=δ
• Sustituyendo vz:

gx
 v dx
z
vz 
æ
1 ç ÷
2ö
r g 2sen ç æ x ö ÷
z
gx = g cos α
gz = g sen α
vz  0
2m ç èø ÷ x
α
 è ø
 é æ ö2 ù
rg  2 x


sen ê1 ç ÷ ú dx
v z
dx 2m ê èø ú
ë û rg   é æxö ù
2
sen  ê1 ç ÷ ú dx 
0
vz  0
 
  2m ê èø ú
0 ë û
é ù

é ù
rg  ê  x3 ú  r g  sen ê   ú  rg  sen 2  rg  sen
2
sen x  2
2m ê 0 3 ú 2m ë 3 û 2m 3 3m
ë 0û
• Es decir:
r g 2
vz  sen
3m
dx
Velocidad media
gz x=δ
rg  2
vz  sen gx
3m gx = g cos α
gz = g sen α
z
• Por tanto, la relación entre la velocidad x
α
máxima y la media es la siguiente:
rg 2sen vz 3 2
vz  max
 v z  vz
max 2m vz 2 3 max
Caudal volumétrico
• El caudal volumétrico correspondiente al nivel macroscópico se determina como el
producto de la velocidad media por la sección de paso:
r g 2 rgW 3

Qv  Sv z  W  
3m
sen 
3m
sen
dx
Espesor de la película 
gz x=δ
• La determinación del espesor de la película
en función de las variables determinables a gx
escala macroscópica responde a las gx = g cos α
gz = g sen α
siguientes expresiones: z
α
x
• En función de la velocidad media:
r g 2 3mv z
vz  sen 
3m rg sen
• En función del valor del caudal volumétrico:
r gW 3 3mQv
Qv  sen  3
3m rgWsen
dx
Fuerza que ejerce el fluido sobre la superficie sólida
gz x=δ
• La componente z de la fuerza (Fz) sobre la
superficie se obtiene evaluando la densidad gx
de flujo de cantidad de movimiento sobre la gx = g cos α
gz = g sen α
interfase sólido-fluido: z
α
x
dv z
dx
Fz  L W  xz x 
 LW r g sen
• Como puede observarse el valor obtenido coindice exactamente con la componente z

de la fuerza peso del fluido contenido en la película (m·gz)
• Todos los resultados obtenidos son válidos únicamente cuando la película desciende
en flujo laminar en líneas de corriente rectas. Estas condiciones se cumplen para el
flujo lento de películas viscosas delgadas.
2.2. Flujo de Couette totalmente desarrollado
• Considere el caso del flujo laminar y estacionario de un fluido newtoniano
incompresible en la rendija estrecha formada por dos placas paralelas infinitas.
• La placa superior se mueve con velocidad uniforme vo, y la placa inferior está en
reposo. La distancia entre estas dos placas es h y la gravedad actúa en la dirección z
negativa (normal hacia el plano de la pantalla).
• No existe gradiente de presión aplicado que empuje al flujo en la dirección x.; el flujo
se establece debido a los esfuerzos viscosos que se provocan por el movimiento de la
placa superior. Este flujo se llama flujo de Couette.
vo Placa móvil HIPÓTESIS DE PARTIDA:
h • Régimen Laminar
Fluido
• Fluido puro
y
• Flujo estacionario e isotérmico
Placa fija x • Densidad y viscosidad constantes
• Las placas son infinitas en “x” y “z”
Actividad propuesta: determine la distribución de presión y de velocidad. Estime la densidad de flujo de

cantidad de movimiento y la fuerza de corte por unidad de área que actúa sobre la placa inferior
vo Placa móvil
ECUACIONES DE VARIACIÓN:
h
Ecuación de continuidad (ρ = cte ): Fluido
y
Placa fija x

Dv  
r  m 2 v  P + r g
Dt
æ v v v x v x ö æ 2v 2v x 2v x ö
p
rç x + vx x + vy + vz ÷   + m çç 2 + 2 + 2 ÷÷ + r g x
x
è t x y z ø x è x y z ø
æ v v v y v y ö æ 2v 2v y 2v y ö
p
÷   +mç ÷+ rg
r çç y + v x y + vy + vz ÷
y
+ +
è t x y z ø y ç x 2
y 2
z 2 ÷ y
è ø

p
z
è t x y z ø z è x y z ø
vo Placa móvil
SIMPLIFICACIONES:
• Flujo paralelo a las placas: h

Fluido
vy = 0 vz = 0 y
Placa fija x
• Régimen estacionario:
P
• La presión permanece constante para cada valor de “x”: 0
x
v x
• El campo de velocidad es bidimensional: 0
z
• La gravedad actúa en la dirección “z” negativa: gx = gy = 0 y gz= -g

Ecuación de continuidad (ρ = cte ):
v y = 0 y vz = 0
v x
0 vx = vx(y)
x
vy = 0 vz = 0
vo Placa móvil
Ecuación de movimiento: h
Fluido
y
Componente X: Placa fija x
Continuidad P = cte respecto a x gx = 0
æ v v x v x v x ö æ 2v 2v x 2v x ö 2v x

p
rç x + vx + vy + vz ÷   + m çç 2 + 2 + 2 ÷÷ + r g x
x
0
è t x y z ø x è x y z ø y 2
f. estacionario Vy = 0 Vz = 0 Continuidad Campo velocidad

bidimensional
Componente Y:
Vy = 0 gy = 0
æ v v v v ö æ 2v  2
v  2
vy ö
p ç ÷ + rg p
r çç y
+ vx y
+ vy y
+ vz y
÷   +m y
+ y
+  0
è t x y z ÷ø y ç x 2 y 2 z 2 ÷ y
y
è ø
f. estacionario Vy = 0 Vy = 0 Vy = 0 Vy = 0 Vy = 0
p = p(z)
vo Placa móvil
h
Fluido
Ecuación de movimiento: y
Componente Z: Placa fija x

p
z
p
è t x y z ø z è x y z ø  rg
z
f. estacionario Vz = 0 Vz = 0 Vz = 0 Vz = 0 Vz = 0 Vz = 0
• Por tanto, las ecuaciones de variación simplificadas son:

2v x p
vx = vx(y) 0 p = p(z)  rg
y 2
z
• A partir de la componente “z” de la ecuación de movimiento simplificada, separando
variables e integrando:
dp p = -rgz + C1
 rg siendo C1 una constante arbitraria.
dz
ECUACIONES DE VARIACIÓN SIMPLIFICADAS: vo Placa móvil
vx = vx(y) 2v x p
0 p = p(z)  rg h
y 2
z Fluido
y
Distribución de velocidades Placa fija x
• A partir de la componente “x” de la ecuación de movimiento simplificada, separando

variables:
d 2v x dv x
0  C2 v x  C2 y + C3
dy 2 dy
vo Placa móvil
• Condición de contorno:
h
Fluido
Placa inferior: y = 0 vx = 0 y
Placa superior: y = h v x = vo
Placa fija x
• Por tanto, C3 = 0 y C2 = vo/h: vo
vx  y
h
vo Placa móvil
• La única componente no nula de la velocidad es h
vx y sólo se produce variación de la misma en la Fluido
dirección y, por consiguiente: y
dv v
 yx  m x  m o Placa fija x
dy h
• Representa la fuerza por unidad de área que actúa sobre el elemento diferencial de
volumen de fluido. Actúa en la dirección x negativa, ya que la placa “tira de él” hacia la
derecha.
Fuerza de corte por unidad de área que actúa sobre la placa inferior
• La fuerza por unidad de área que actúa sobre la pared es igual y opuesta a τyx (tercera
ley de Newton) :
F v
m o
A h
• Actúa en la dirección x positiva, ya que el fluido intenta “tirar” de la placa hacia la
izquierda debido a los esfuerzos viscosos (fricción).
Utilidad del flujo de Couette
Cilindro exterior en reposo
• El flujo de Coutte ocurre en el interior de un Fluido: m, r
viscosímetro rotacional, que es un instrumento que
se utiliza para medir la viscosidad. w
Ri
• Consta de dos cilindros concéntricos de longitud L:
un cilindro sólido interior rotacional de radio Ri y un
Ro
cilindro hueco exterior estacionario de radio Ro.
• La rendija formada por los dos cilindros es muy

pequeña (Ro – Ri) << Ro y contiene el fluido cuya
viscosidad se va a medir. Cilindro interior en rotación
¿Cómo se obtiene el valor de la viscosidad?

• En una medición de viscosidad, se mide la velocidad angular del cilindro interior w, así
como el torque (Taplicado) necesario para rotar al cilindro.
• Del ejemplo anterior se sabe que el esfuerzo cortante que actúa sobre un elemento de
fluido adyacente al cilindro interior es aproximadamente igual:
dv x v
 yx  m  m o
dy h
Utilidad del flujo de Couette Cilindro exterior en reposo

Fluido: m, r
=
• donde la velocidad de la placa superior móvil se w

Ri
sustituye por la velocidad de giro, en sentido contrario
a las manecillas del reloj, de la pared rotacional del
cilindro interior. Ro
• En estado estacionario, el torque en sentido de las

manecillas del reloj (torque viscoso), que actúa sobre
la pared del cilindro interior debido a la viscosidad del
fluido, se equilibra con el torque aplicado en sentido Cilindro interior en rotación
contrario a las manecillas del reloj (conocido):
siendo el torque viscoso es igual a este esfuerzo cortante por el área de la pared,
por el brazo de palanca.
w Ri
Tvis cos o  Taplicado  m  2Ri L  Ri
Ro  Ri
2.3. Flujo en un tubo circular vertical
p0 HIPÓTESIS DE PARTIDA:
r • Régimen Laminar
z • Fluido puro
• Densidad y viscosidad constantes
L • Se considerará que la longitud del tubo es muy
R
grande respecto al radio del tubo, de modo que
los “efectos finales" carezcan de importancia a lo
largo de la mayor parte del tubo.
Actividad propuesta: determine la distribución de presión y

de velocidad, velocidad máxima, velocidad media, caudal
pL volumétrico, densidad flujo de cantidad de movimiento y
fuerza que ejerce el fluido sobre la superficie sólida
2.3. Flujo en un tubo circular
• Régimen Laminar
• Fluido puro
• Se considerará que la longitud del tubo es muy grande respecto al radio del tubo,
de modo que los “efectos finales" carezcan de importancia a lo largo de la mayor
parte del tubo.
r 1  1  
+
t r r
 
rr v r +
r q
 
rv q +
z
 
rv z  0
æ v v r v q v r v q2 v r ö
Componente r: rç r + vr +  + vz ÷   p
ç t r r q r z ÷ø r
è
é æ ö 1 2v v q 2v r ù
 1  2
+m ê ç
r
êë è r r
 
rv r ÷ +
ø r q2 2
r 
r 2 q
+
z úû
2
ú + rgr
æ v v q v q v q v r v q v q ö 1 p
rç q + vr + + + vz
Componente θ: ÷
è t r r q r z ø r q
é æ ö 1 2v v r 2v q ù
 1  2
+m ê ç  
êë r è r r
rv q ÷ +
2
ø r q
q
2
+ +
r q z úû
2 2
ú + rgq
æ v v v v v ö p
Componente z: r ç z + vr z + q z + v z z ÷  
è t r r q z ø z
é æ v ö 1 2v 2v ù
1 
+m ê ç r z ÷ + 2 2z + 2z ú + rgz
êë r r è r ø r q z úû
SIMPLIFICACIONES:
• Movimiento unidireccional (eje z):
v r = 0 y vq = 0
• Régimen estacionario: v z r
0 0
t t
• Por simetría: v z Si no se cumpliese esta condición de simetría, el fluido

0
q se deformaría y no tendría geometría cilíndrica
• Las componentes “r” y "q" de la aceleración de la gravedad son nulas:
gr = 0; gq= 0; gz=g
Z
r 1  1  
+
t r r

rr v r +
r q

rv q +
z
 
rv z  0   z
(x, y, z)
(r, q, z)
estacionario vr = 0 vq = 0
 
z
 
rv z  0
z
vz  0 (1) q r
Y
r constante
vz = vz(r)
X
vr = 0 vr = 0 vq = 0 vr = 0
æ v v v 2
v r ö
÷   p
v v
Componente r: r ç r + vr r + q r q
 + vz
ç t r r q r z ÷ø r
è
é æ ö 1 2v v  2 ù
estacionario  1  2 vr
+m ê ç
êë r è r r
 
rv r ÷ +
ø r q2 2
r  q +
r q z úû
2 2
ú + rgr
vr = 0 vr = 0 vq = 0 vr = 0 qr = 0
p
 0 (2)
r
vq = 0 vq = 0 vq = 0 vq = 0
rç q + vr + + + vz
Componente q: ÷
estacionario é æ ö 1 2v v r 2v q ù
 1  2
+m ê ç
êë r è r r
 
rv q ÷ +
2
ø r q 2
q + +
r q z úû
2 2
ú + rgq
1 p
 0 (3) vq = 0 vq = 0 vr = 0 vq = 0 gq = 0
r q
vr = 0 vq = 0
continuidad
æ v v z v q v z v z ö p
Componente z: rç z + vr + + vz ÷
è t r r q z ø z continuidad
estacionario é æ v ö 1 2v 2v ù
1 
+m ê çr z ÷ + 2 z + z ú + rg
z
êë r r è r ø r q2 z 2 úû
gz = g
simetría
p é 1  æ v öù
 + mê ç r z ÷ú + rg  0 (4)
z êë r r è r øúû
2.3. Flujo en un tubo circular p0
r
 p 1 p p é 1  æ v öù z
vz  0  0  0  + mê ç r z ÷ú + rg  0
z r r q z êë r r è r øúû
L
(1) (2) (3) (4) R
• Teniendo en cuenta las relaciones obtenidas para la presión:
p p pL
 0 p  f r , q 
r q
p é 1  æ v öù
 mê ç r z ÷ú + rg p p  f (z)
z êë r r è r øúû
• Integrando resulta:
pL L
é 1  æ  v z öù L
z0 p  po
p dP  0 m êë r r ççè r  r ÷÷øúû dz + 0 r g dz zL p  pL
o
r
é 1  æ  vz öù z
p L  po  m ê çç r ÷÷ú L + r g L
ë r r è  r øû
p  p é 1  æ v öù L
R
• Reordenando la ecuación: L o
 r g  mê ç r z ÷ú
L êë r r è r øúû
• Si se realiza el cambio de variable p = P + rgz
æ presión ö æ presión ö æ presión ö pL

ç ÷ç ÷+ç ÷
ètotal ø è del fluido ø è hidrostática ø
• Se tiene que: z0 po   o + rg 0 p L  po  + rgL   o  o

rg  L rg  L
zL pL   L + rgL L L L
• Y por tanto, se obtiene la siguiente distribución de presiones:
PL  Po é 1  æ v öù
 mê ç r z ÷ú
r
• Teniendo en cuenta que, según las simplificaciones del ejercicio y z
el resultado de la ecuación de continuidad:
v r = 0 y vq = 0 Movimiento unidireccional (eje z) L
R
v z v z Por simetría y a partir de la ecuación de
 0
q z continuidad
• La única componente no nula del tensor de esfuerzo cortante será rz:

pL
dv
 rz  m z
dr
• Y sustituyéndola en la expresión de la distribución de presiones:
PL  Po é 1  æ v öù
 mê ç r z ÷ú
Po  P L é1 
P L  Po é 1  æ  rz öù
 m ê çç r ÷÷ú ê r rz ùú
L ë r r è m øû
L ë r r û
r
Po  PL é 1  ù z
L
ê
ë r r
 
r  rz ú
û
L
• Integrando: R
Po  P L Po  P L r 2
 
 rz  L r r
 r   r  rz 
L 2
+ C3
P  P L r C3
 rz  o + pL
L 2 r
• La condición de contorno es: r 0  rz 0 C3 = 0

r 0
• Por con siguiente:

En consecuencia la distribución de esfuerzo
P  PL
 rz  o r cortante resulta una función lineal respecto del
2L radio, aumentando desde el centro de la conducción
hacia las paredes.
r
z
vz = 0
DISTRIBUCIÓN DE
VELOCIDADES VZ (r)
vz,máx
rz = 0
Po  PL
 rz  r
2L
rz, máx
DISTRIBUCIÓN DE
DENSIDAD DE FLUJO
DE CANTIDAD DE
MOVIMIENTO
r
z
• A partir del perfil del esfuerzo cortante y considerando la
definición de rz puede obtenerse la distribución de velocidades:
L
dv Π  ΠL Πo  ΠL r 2 R
τ rz  μ z  o r vz   + C4
dr 2L 2μL 2
• La correspondiente condición de contorno a la interfase sólido-

fluido permite evaluar la constante de integración:
pL
P  PL 2
r R vz  0 C4  o R
4m L
• Sustituyendo en la ecuación del perfil de velocidades, se obtiene:
æ P  P ö 2 2 æ P P ö 2 é æ r ö2 ù
 
v z  ç 0 L ÷ R r  ç 0 L ÷R ê1 ç ÷ ú
è 4m L ø è 4m L ø êë è R ø úû
r
z
vz = 0
æ P P ö 2 é æ r ö 2 ù
v z  ç 0 L ÷R ê1 ç ÷ ú
è 4m L ø êë è R ø ú
û
DISTRIBUCIÓN DE
VELOCIDADES VZ (r)
vz,máx
rz = 0
rz, máx
DISTRIBUCIÓN DE
DENSIDAD DE FLUJO
DE CANTIDAD DE
MOVIMIENTO
Velocidad máxima
é 2ù
Po  PL 2 ê æ r ö ú
vz  R 1 ç ÷
4m L ê èR ø ú
ë û
• Aplicando la condición de máximo (dvz/dr)=0 y sustituyendo:
dv z Π  ΠL dv z
0 o 2r  0 si r  0
dr 4μL dr
Es decir, el valor máximo de la componente z de la velocidad se alcanza en r = 0

Se trata de un máximo ya que:
d 2v z Po  PL
 2 0
dr 2 4m L
• De acuerdo con la distribución de velocidad obtenida:
Po  PL 2
v z,máx  R
4m L
La velocidad varía entre este valor máximo y vz = 0 para r = R
Velocidad media
• Esta velocidad media corresponde con la observada o determinada a nivel
macroscópico.
• Esta velocidad puede determinarse cuando se dispone del perfil de velocidades,

calculando la media referida a la sección transversal de la película
vz 
1
S
v z dS 
S
• La sección transversal del elemento diferencial de volumen considerado es una corona

circular
R
2   v z rdr R
2
vz  0
 2  v z rdr
 R2 R 0
dr
dq ds
r
Arco =r·dq
dS = 2 π r dr
Velocidad media
• Sustituyendo vz (perfil radial de velocidad):
2 Πo  ΠL 2 é æ r ö ù
2
2 é Po  P L 2 Po  P L 2 r3 ù
R R R
vz  2  R ê1  ç ÷ ú rdr  2 ê  R rdr   R 2 dr ú 
R 0 4μL êë è R ø úû R ë 0 4m L 0
4m L R û
2 éP  P 2 R
Po  P L r 4 ù 2 Po  P L é R4 R4 ù 2 Po  P L R4
R
2r
ê o L
R  ú ê  ú 2
R2 êë 4 m L 2 0
4m L 4 0 ú R 2 4m L ë 2
û
4 û R 4m L 4
• Por tanto:
Po  PL 2
vz  R
8m L
• Comparando con el valor obtenido para vz,máx:
Po  PL 2 1
v z,máx  R v z  v z,máx
4m L 2
Caudal volumétrico
• El caudal volumétrico (determinado macroscópicamente) es el producto de la
velocidad media de la corriente de fluido por la sección de paso,
éP  P 2ù 2 R4
Qv  v zS  ê o
ë 8m L
L
R ú R 
û 8m L
P o  PL  
Ecuación de Hagen-Poiseuille
• Teniendo en cuenta que será el producto del área lateral (2RL) por el valor de la
fuerza por unidad de área en la pared (que coincide con el valor del esfuerzo cortante
evaluado en la pared rz|r=R):
Po  P L
Fz  2 R L  rz
r R
 2 R L
2L

R  R 2 Po  P L 
• Deshaciendo el cambio de variable que se realizó anteriormente para la presión:
z0 P o  po  rg0
zL PL  pL  rgL
   
Fz  R 2 po  pL + rgL  R 2 po  pL + R 2rgL
   
Fz  R 2 po  pL + rgL  R 2 po  pL + R 2rgL
• La fuerza neta que actúa en el sentido de la corriente sobre el cilindro de fluido,

debido a la aceleración de la gravedad y a la diferencia de presión, se equilibra
exactamente por la fuerza viscosa Fz, que tiende a oponerse al movimiento del
fluido.
2.4. Flujo descendente a lo largo de la pared exterior de un tubo vertical
• Considere el caso de un tubo vertical, de radio exterior R, que está lleno de un líquido
que asciende por el interior del mismo, para descender después por la parte exterior
del tubo en forma de película de espesor uniforme , tal y como se indica en la figura.
• Se desprecian los efectos finales en la parte superior del tubo y la posible influencia
de la fase gaseosa sobre el líquido descendente.
• Supóngase, asimismo, que la pérdida de presión debida al rozamiento es
despreciable y, por tanto, la presión se mantiene constante en la dirección z.
r
z
R 
Actividad propuesta: determine a) la distribución de velocidades en la película externa del
tubo, b) la tensión de rozamiento en la parte externa del tubo, c) el caudal volumétrico en
la mencionada película.
Z
(x, y, z)
z (r, q, z)
• Fluido puro y newtoniano
• Flujo estacionario e isotérmico Y
q r
• Presión constante en la dirección “z” X
• Se considerará que la longitud del tubo es muy grande respecto al radio del tubo,
de modo que los “efectos finales" carezcan de importancia a lo largo de la mayor
parte del tubo.
r 1  1  
+
t r r

rr v r + 
r q
rv q +  
z
rv z  0  
ç t r r q r z ÷ø r
è
é æ ö 1 2v v q 2v r ù
 1  2
+m ê ç
r
êë è r r
 
rv r ÷ +
ø r q2 2
r 
r 2 q
+
z úû
2
ú + rgr
rç q + vr + + + vz
é æ ö 1 2v v r 2v q ù
 1  2
+m ê ç  
êë r è r r
rv q ÷ +
2
ø r q
q
2
+ +
r q z úû
2 2
ú + rgq
é æ v ö 1 2v 2v ù
1 
+m ê ç r z ÷ + 2 2z + 2z ú + rgz
Z
SIMPLIFICACIONES:
• Movimiento unidireccional (eje z):
v r = 0 y vq = 0
(x, y, z)
z (r, q, z)
• Régimen estacionario: v z r
0 0
t t
Y
q r
v z
• Por simetría: 0
q X

gr = 0; gq= 0; gz=g
• La presión se mantiene constante en la dirección z.

Z
r 1  1  
+
t r r

rr v r + 
r q
rv q +
z
 
rv z  0   z
(x, y, z)
(r, q, z)
estacionario vr = 0 vq = 0
 
z
 
rv z  0
z
vz  0 (1) q r
Y
r constante
vz = vz(r)
X
vr = 0 vr = 0 vq = 0 vr = 0
ç t r r q r z ÷ø r
è
estacionario é æ ö 1 2v v q 2v r ù
 1  2
+m ê ç
êë r è r r
 
rv r ÷ +
ø r q2 2
r  +
r q z úû
2 2
ú + rgr
vr = 0 vr = 0 vq = 0 vr = 0 qr = 0
p
 0 (2)
r
vθ = 0 vθ = 0 vθ = 0
rç q + vr + + + vz
estacionario vθ = 0 é æ ö 1 2v v r 2v q ù
 1  2
+m ê ç
êë r è r r
 
rv q ÷ +
2
ø r q
q
2
+ +
r q z úû
2 2
ú + rgq
vθ = 0 vθ = 0 vr = 0 vθ = 0 gq = 0
1 p
 0 (3)
r q vr = 0 vθ = 0 continuidad p = cte en z
è t r r q z ø z continuidad
1 
+m ê çr z ÷ + 2 z + z ú + rg
z
êë r r è r ø r q2 z 2 úû
simetría
é 1  æ v öù
mê ç r z ÷ú + rg  0 (4)
êë r r è r øúû
ECUACIONES DE VARIACIÓN SIMPLIFICADAS: r
z
 p 1 p é 1  æ v öù
vz  0  0  0 mê ç r z ÷ú + rg  0
z r r q êë r r è r øúû
(1) (2) (3) (4)
• A partir de la integración de la ecuación (4):
é 1  æ v öù R 
 æ v z ö rg d æ dv z ö rg
mê ç r z ÷ú + rg  0 çr ÷ r çr ÷ r
êë r r è r øúû r è r ø m dr è dr ø m
dv z rg r 2 dv z rg r C1 rg r 2
r  + C1  + vz   + C1 Lnr + C2
dr m 2 dr m 2 r 2m 2
• Las condiciones límite son:
r R vz  0
v z
r R+  rz 
r
0 Interfase líquido  gas Plíquido  Pgas  
ij líquido 0
Distribución de velocidades r
z
dv z rg r C1 rg r 2
 + vz   + C1 Lnr + C2
dr m 2 r 2m 2
• Constante C1:
v z dv z rg r C1
r R+  rz  0  +
r dr m 2 r
rg R +  C rg
 
2
0 + 1 C1  R+ R 
m 2 R+ 2m
• Constante C2:
rg r 2
r R vz  0 vz   + C1 Lnr + C2
2m 2
rg R 2 rg r g éR2 ù
   
2 2
0 + R +  LnR + C2 C2  ê  R +  LnR ú
2m 2 2m 2m ë 2 û
Distribución de velocidades r
z
rg r 2
vz   + C1 Lnr + C2
2m 2
• Sustituyendo en vz los valores de C1 y C2:
rg r g éR2 ù
   
2 2
C1  R+ C2  ê  R +  LnR ú
2m 2m ë 2 û
R 
rg r 2 r g r g éR2 ù
   
2 2
vz   + R +  Ln r + ê  R +  LnR ú
2m 2 2m 2m ë 2 û
rg é r 2 R2 ù
   
2 2
vz  ê + R +  Ln r +  R +  LnR ú
2m ë 2 2 û
r g éR2  r 2 rù
 
2
vz  ê + R +  Ln ú
2m ë 2 Rû
Tensión de rozamiento r
z
dv z
 rz  m
dr
• Conocemos el valor de la derivada anterior a partir del cálculo de
la distribución de velocidad:
rg
 
2
dv z rg r C1 con C1  R+
 + 2m
dr m 2 r
é 2ù R 
 rz 
rg ê
r 

R+  ú
2ê r ú
êë úû
Caudal volumétrico
• El caudal volumétrico (determinado macroscópicamente) es el producto de la
velocidad media de la corriente de fluido por la sección de paso,
é ù
 
2
Qv  v zS   ê R +   R 2 úv z
ë û
2.4. Flujo descendente a lo largo de la pared exterior de un tubo vertical r
z
Caudal volumétrico
é ù
 
2
Qv  v zS   ê R +   R 2 úv z
ë û
• Tendremos que calcular la velocidad media:
vz 
1
S
v z
dS
S
R 
• La sección transversal del elemento diferencial de volumen considerado es una corona
circular
R +
 2  v rdr z
2
R +
vz 

R
 R +    R 2 2
 
R +  
2
R 2  v rdr
R
z
dr
• Sustituyendo la distribución de velocidad en la integral anterior: dq ds

r
Arco =r·dq
dS = 2 π r dr
Caudal volumétrico
2 R+
r g éR2  r 2 rù
  
2
vz  ê + R +  Ln ú r dr
 R + 
2
R 2
R 2m ë 2 Rû
rg éR+ R 2  r 2 R+
r ù
ê 
r dr +  R +  Ln r dr ú 
2
vz 
é ù
  2 R
2
mê R +   R2úë R R û
ë û
rg é R 2 R+ 1 R+ 3 R+

r ù
 r dr  2  r dr + R +    
2
vz  ê Ln r dr ú
é 2 ùë 2 R
  R
2
mê R+ R ú R R û
ë û
rg éR 2 é 2ù 1é 4ù
R+
r ù
      
2 4 2
vz  ê ê R +   R ú  ê R +   R ú + R+ Ln r dr ú
é ù ë û ë û
 
R +   R2úë 4 8 R
2
mê R û
ë û
Integral I
Caudal volumétrico
R+
r
I  Ln
R
r dr
R
v  Ln
r 1
dv  dr  
d uv  udv + vdu  
vdu  d uv  udv
R r
R+ R+ R+ R+
r2
 v du   d uv    u dv  uv  u dv
R+
du  r dr u 
R
2 R R R R
   
R+ R+ R+ 2 2
r2
r R+ 2
r 1 r 2
r r 2 R+ R+ R 2
R R+ R2
I  Ln
2 RR
  2 r
dr  Ln
2 RR

4

2
Ln
R

2
Ln 
R 4
+
4
R R
 R +   
2 2
R+ R + R
2
I Ln 
2 R 4
Caudal volumétrico
rg éR 2 é 2ù 1é 4ù
R+
r ù
      
2 4 2
vz  ê ê R +   R ú  ê R +   R ú + R+ Ln r dr ú
é ù ë û ë û
 
R +   R2úë 4 8 R
2
mê R û
ë û
Integral I
 R +   
2 2
R+ R+ R
2
I Ln 
2 R 4
é 2 2ù
     
4 4
rg êR é R +   R2 ù  1é R +   R4 ù + R +  R+ R+ R2 R +  ú
   
2 4
vz  ê ê ú ê ú Ln  + ú
é ù ë û ë û
 
R +   R 2 ú êë 4 8 2 R 4 4
2
mê úû
ë û
• Desarrollando la ecuación anterior se llega a:
   
4
rg é 2 ù rg R+ R+
 
2
vz  êR  3 R +  ú+ Ln
ë û
 
2
8m 2m R +   R 2 R
Caudal volumétrico
• Por lo tanto, el caudal será:
é ù
 
2
Qv  v zS   ê R +   R 2 úv z
ë û
   
4
rg é 2 2ù r g R+ R+
vz  ê
ë 
R  3 R +  ú+
û  Ln
 
2
8m 2m R +   R 2 R
Qv 
 rg
R +  2  R 2 R 2  3R +  2 +  r g R +  4 Ln R +  
  
8m 2m R
2.5. Forma de la superficie libre de un líquido que gira
• Considere el caso de fluido de densidad y viscosidad constantes que está contenido
en un recipiente cilíndrico vertical de radio R.
• El recipiente gira con velocidad angular W constante alrededor de su eje vertical.
Debido al efecto centrífugo del movimiento, el líquido tiende a tomar un nivel más alto
en el borde del recipiente y más bajo en el centro.
• Suponga que el movimiento ocurre en régimen

laminar, estacionario, con simetría axial.
• Considere que el líquido se mueve en
trayectorias circulares concéntricas con el eje z
Actividad propuesta: determine la distribución de velocidades y de presión y la forma de la

superficie libre una vez alcanzado el estado estacionario.
Z
• Fluido puro y newtoniano (x, y, z)
z (r, q, z)
Y
• Trayectorias circulares concéntricas con z. q r
• No existe gradiente de presión en la dirección q

X
r 1  1  
+
t r r

rr v r + 
r q
rv q +
z
 
rv z  0  
ç t r r q r z ÷ø r
è
é æ ö 1 2v v q 2v r ù
 1  2
+m ê ç
r
êë è r r
 
rv r ÷ +
ø r q2 2
r 
r 2 q
+
z úû
2
ú + rgr
rç q + vr + + + vz
é æ ö 1 2v v r 2v q ù
 1  2
+m ê ç  
êë r è r r
rv q ÷ +
2
ø r q
q
2
+ +
r q z úû
2 2
ú + rgq
é æ v ö 1 2v 2v ù
1 
+m ê ç r z ÷ + 2 2z + 2z ú + rgz
SIMPLIFICACIONES: Z
• Movimiento circular (trayectorias

concéntricas en planos z =cte):
v r = 0 y vz = 0
(x, y, z)
z (r, q, z)
• Régimen estacionario: v q r
0 0
t t Y
q r
• No existe gradiente de presión aplicado en la

dirección θ:
gr = 0; gq= 0; gz=g
• La línea central del recipiente es un eje de simetría, lo que requiere que vθ sea
independiente de z. Si esta condición no se cumpliese, el fluido se deformaría-
Z
r 1  1  
+
t r r

rr v r +
r q
rv q +  
z
rv z  0   z
(x, y, z)
(r, q, z)
estacionario vr = 0 vz = 0
 vq vq = vq(r) (1) Y
0 q r
q
X
vr = 0 vr = 0 vz = 0
ç t r r q r z ÷ø r
è continuidad
estacionario é æ ö 1 2v v q 2v r ù
 1  2
+m ê ç
êë r è r r
 
rv r ÷ +
ø r q2 2
r  +
r q z úû
2 2
ú + rgr
vr = 0 vr = 0 vr = 0 qr = 0
2
v p
r q
 (2)
r r
continuidad vr = 0 vz = 0 P es constante en la
dirección θ
rç q + vr + + + vz
estacionario vr = 0 é æ ö 1 2v v r 2v q ù
 1  2
+m ê ç
êë r è r r
 
rv q ÷ +
2
ø r q
q
2
+ +
r q z úû
2 2
ú + rgq
continuidad vr = 0 geometría gq = 0
 æ1  ö
0m ç
r è r r
 
r vq ÷
ø
(3)
vz = 0 vz = 0 vz = 0
1 
+m ê çr z ÷ + 2 z + z ú + rg
z
êë r r è r ø r q 2
z úû
2
vz = 0 vz = 0 vz = 0
p
0 r g (4)
z
v q2 p  æ1  ö p
vq = vq(r) r 
r r
0m ç
r è r r
r 
v q ÷
ø
0
z
r g
(1) (2) (3) (4)

d æ1 d ö 1d d r2
ç
dr è r dr
 
r vq ÷  0
ø r dr
 
r v q  C1
dr
 
r v q  C1r r v q  C1 + C2
2
C C
vq  1 r + 2
2 r
• Las condiciones límite son: • Como vq debe ser finito, para r = 0, C2=0
r 0 vq  0
• Con r = R se obtiene que C1 = 2W. Por
r R v q  WR tanto:
v q  Wr
v q2 p  æ1  ö p
vq = vq(r) r 
r r
0m ç
r è r r
r 
v q ÷
ø
0
z
r g
(1) (2) (3) (4)
v q2 W2r 2 1 2 2
 p   r r r  p   r r r p rW r + f  z 
2
• Derivando la expresión anterior:
p f  z 

z z
• E igualándola con (4):
f  z 
 rg f  z   rgz + C3
z
1 2 2
p rW r  rgz + C3
2
1 2 2
p rW r  rgz + C3
2
• Las condiciones límite son:
r 0 p  po z  zo
• Por tanto C3 viene dado por:
C3  po + rgzo
• Y la distribución de presión en la superficie:

1 2 2
p rW r  rg  z  zo  + p0
2
• Como la superficie libre se define como el lugar geométrico de los puntos en los que
p = po:
1 2 2
z  zo  Wr
2g

Tema 5. Distribuciones de Velocidad 2023-24

Cargado por

Copyright:

Formatos disponibles

Tema 5. Distribuciones de Velocidad 2023-24

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tema 5. Distribuciones de Velocidad 2023-24

Cargado por

Copyright:

Formatos disponibles

FENÓMENOS DE TRANSPORTE

BLOQUE A. PRINCIPIOS MATEMÁTICOS Y BALANCES

BLOQUE B. DISTRIBUCIONES DE PROPIEDAD Y ANALOGÍAS

• 2. Aplicación de las ecuaciones de variación a la obtención de perfiles de

• Bird, R.B.; Stewart, W.E.; Lightfoot, E.N. “Fenómenos de Transporte”.

• Costa, E.; et al. “Ingeniería Química, volumen II: Fenómenos de

1. Se simplifican las ecuaciones de variación para el caso bajo estudio.

CONDICIÓN LÍMITE DE NO-DESLIZAMIENTO:

En el caso especial de una pared en reposo, el fluido

Para los casos en los que también se consideren los

CONDICIÓN LÍMITE DE INTERFASE:

Cuando dos fluidos A y B se encuentran en una interfase, las condiciones frontera de la

Esta condición se aplica en las interfases

Además, si se consideran los efectos de

Cuando dos fluidos A y B se encuentran en una interfase, las condiciones frontera de la

Como magua es casi 50 veces mayor que maire, para

• Este estudio tiene interés para los siguientes casos:

Actividad propuesta: determine la distribución de presión y de velocidad, la densidad de flujo de

æ v v z v z v z ö æ 2v 2v z 2v z ö

• Campo de velocidad bidimensional:

• La componente “y” de la aceleración de la gravedad es nula: gy = 0

• La presión permanece constante para cada valor de “z”, es P

Ecuación de continuidad (ρ = cte ): gx

Ecuación de continuidad (ρ = cte ):

Componente Z: p = cte para cada z continuidad

æ v v z v z v z ö æ 2v 2v z 2v z ö æ  2v z ö

• Teniendo en cuenta la condición de contorno: x = 0, p = patm se obtiene:

dv z r g sen rg sen 2 gz x=δ

Interfase líquido  gas Plíquido  Pgas  

La condición de contorno en este caso es que el esfuerzo en la superficie de la fase

dv z r g sen r g sen 2 gz x=δ

No hay deslizamiento entre el fluido y la pared. En el caso de una pared en reposo, el

• Por consiguiente, la ecuación que permite conocer el perfil de velocidades es:

r g sen 2 r g sen 2 r g sen 2

A partir del perfil de velocidad:

• La densidad de flujo de cantidad de movimiento es nula en la superficie libre (x = 0) y

Densidad de flujo de cantidad de movimiento

Es decir, el valor máximo de la componente z de la velocidad se alcanza en x = 0

• Como puede observarse el valor obtenido coindice exactamente con la componente z

vo Placa móvil HIPÓTESIS DE PARTIDA:

Actividad propuesta: determine la distribución de presión y de velocidad. Estime la densidad de flujo de

æ v v z v z v z ö æ 2v 2v z 2v z ö

• Flujo paralelo a las placas: h

ECUACIONES DE VARIACIÓN SIMPLIFICADAS:

Componente X: Placa fija x

Continuidad P = cte respecto a x gx = 0

æ v v x v x v x ö æ 2v 2v x 2v x ö 2v x

f. estacionario Vy = 0 Vz = 0 Continuidad Campo velocidad

Componente Z: Placa fija x

æ v v z v z v z ö æ 2v 2v z 2v z ö

• Por tanto, las ecuaciones de variación simplificadas son:

Distribución de velocidades Placa fija x

• A partir de la componente “x” de la ecuación de movimiento simplificada, separando

• La rendija formada por los dos cilindros es muy

¿Cómo se obtiene el valor de la viscosidad?

Utilidad del flujo de Couette Cilindro exterior en reposo

• donde la velocidad de la placa superior móvil se w

• En estado estacionario, el torque en sentido de las

Actividad propuesta: determine la distribución de presión y

• Movimiento unidireccional (eje z):