Facultad de Ciencias Fisicomatemáticas e Ingeniería

Cátedra: Mecánica de Fluidos Práctica N 8 – Flujo laminar

Docentes: Dra. Miralles / Ing. Rubén Galdeano / Ing. Eduardo Contento Revisión: 02

1. Determinar cuál será la velocidad crítica para que un flujo de (a) fuel oíl que fluye a 20°C a través de una tubería
de 15 cm de diámetro (b) agua a 20°C circulando por una tubería de 15 cm de diámetro
Datos: Recrítico = 2300
ρ fuel oíl : 906 Kg/m3 – μ fuel oíl = 3,56 10-3 Kg/ m seg
ρ agua : 998 Kg/m3 – μ agua = 1,005 10-3 Kg/ m seg
Para el fuel oil:
VρLc V Lc 906 × 0 , 15× V
ℜ¿ ó ⇒ 2300= → Vcr í tica=0,060 m/s
μ υ 3,561 0
−3

Para el agua
VρLc V Lc 998 × 0 ,15 × V
ℜ¿ ó ⇒ 2300= → Vcr í tica=0,015 m/s
μ υ 1,0051 0
−3

2.- Un ejemplo de aplicación: Para calcular la viscosidad de un fluido de diseño en un laboratorio del cual se conoce su
densidad, (995 Kg/m3) se utiliza un tubo muy delgado de diámetro interior 1mm, y longitud 1,2 m conectado a un depósito
como se indica en la figura. Al cronometrar el caudal a la salida, se observa que se derraman 6,6 ml en los primeros 10
segundos.

Se pide determinar cuál será la viscosidad buscada para el fluido.

( )
2 2
p2− p1 V 2−V 1
q−w Superficie −g H L = + + g ( z2 −z1 )
ρ 2
Planteamos la ecuación de energía entre la superficie del tanque y la salida del tubo
Eliminando los términos que no son necesarios y escribiendo la ecuación en alturas, tenemos:
2
V2
−H L= + ( z2− z1 )
2g
Eliminamos el termino de las presiones por que tanto 1 y 2 están abiertos a la atmósfera
Por las características del flujo, asumimos que es un flujo laminar y luego lo comprobaremos,
entonces
2 2
V2 −32 μLV m V2
−H L= + ( z2− z1 ) → = + ( z 2−z 1 )
⏟
2
2g γd 2g
ECUACIÓN DE HAGEN

En el único lugar donde tenemos una pérdida de carga es en la conducción de fluido por el tubo,
calculamos la velocidad media con el dato del caudal  6,6 ml / 10 seg por lo que es caudal es 6,6
10-7 m3/seg.
La velocidad entonces se puede calcular como V = 4 Q /  D2  V = 0,84 m/s

2
−32 μ 1 , 2 m0 , 84 m/ s ( 0 , 84 m/ s )
= + ( 0−2 m )
( 995 Kg/m ⋅ 9 ,8 m/ s ) ( 0,001 m ) 29 , 8 m/s 2
3 2 2

2
m s −4 Kg
−2.756 , 64 μ=−1,964 m⇒ μ=7,124610
Kg m⋅ s

Resta comprobar si la hipótesis de flujo laminar se cumple

VρD 0 , 84 × 995× 0,001
ℜ¿ = =1.173 , 1→ La min a r
μ 0,00071246

Flujo laminar incompresible entre placas paralelas:

En el experimento de la figura siguiente Fig.3.2.1, la placa superior se mueve paralela a la
superficie libre y a la dirección del flujo. El flujo se toma unidimensional, es decir, todos los filetes y
sus velocidades en las posiciones y respecto de la base, se repiten iguales para toda la profundidad
z., y no se consideran los efectos del borde de la placa (los valores significativos del experimento se
toman lejos de los bordes).
Como el fluido tiene una componente gravitatoria en la dirección de la longitud l , hay una
variación de presión hidrostática en este sentido. El análisis se hace considerando una partícula

libre de fluido, como la dibujada en el seno del flujo, de dimensiones (  l  y 1 ) es decir de

profundidad unitaria, (por supuesto las escala con que se dibujó la partícula ha sido muy ampliada).
La lámina se mueve con velocidad constante U hacia abajo, , inicialmente el fluido está en reposo.
La ecuación de movimiento a partir del balance de fuerzas, para la placa en la dirección L es:
dp dτ
p δy−( p δy+ δl δy)−τδl+(τδl+ δy δl)+ γ (δl δy 1)⋅ senθ=0
dl dy
Dividimos m.a.m. por el volumen del elemento y simplificamos quedando:

−dp dτ
+ + γ senθ=0
dl dy
Fig.3.2.1

−dh
Como de la figura auxiliar: senθ= →
dl
dτ d d
= ( p+ γh) → dτ=[ ( p+ γ h)]dy
dy dl dl
Como no hay aceleración respecto de y, ya que se trata de un deslizamiento de una placa sobre
un seno fluido, todo el paréntesis es constante para la integración, quedando:
d
τ=y ( p+ γh)+ C1
dl
Como asumimos el comportamiento como flujo laminar, sustituimos acuerdo a la ley de viscosidad
para fluidos:
du
τ =μ
dy
du 1 d C1
Por lo tanto: = ( p+ γh) y +
dy μ dl μ
Que integrando respecto de y nos da el perfil de velocidades entre la base y la placa.
1 d 2 C1
u= (p +γh) y + y +C 2 (3.2.1)
2 μ dl μ
C1 y C2 son constantes de integración que pueden calcularse a partir de las condiciones de
contorno siguientes:
Para y=0  u = 0
Para y=a  u = U
1 d 2 C1
Por lo tanto: C2 = 0 ; ( p+ γh)a + a=U
2 μ dl μ
O sea, de la segunda:
C1 U 1 d
= − ( p+γh )a
μ a 2 μ dl

que reemplazando en 3.2.1 queda:

Uy 1 d 2
u= − [ ( p+ γ h)](a y− y )
a 2 μ dl (3.2.2)

Si el esquema de la figura fuese horizontal, no habrá variación de presión hidrostática y

d U y
( p  h) 0  u 
dl a
es decir daría una distribución lineal, cuidado que en
este caso supone también que no hay variación de
presión, P / L 0 , es decir el conducto plano tiene
áreas normales al movimiento constantes.

Fig.3.2.2
3.- Ejemplo de aplicación
En el esquema de la figura siguiente, una placa móvil se mueve sobre la base fija según se indica:
 0,80 poise

 93,5 UTM / m 3 UTM Unidad técnica

de masa
Determinar:
a) La distribución de velocidades.
b) El caudal.
c) La tensión de corte sobre la placa
móvil.

Fig.3.2.3

Resolución:
a.- la distribución de velocidades:

En el punto superior calculamos ( p   h) :

 2
KgF 4 cm 4 KgF UTM m KgF
p=1 , 5 2
⋅1 0 2
=1 ,5 ×1 0 2
γ = ρ g=93 , 5 3
⋅ 9 , 81 2
=917 , 2 3 ¿
cm m m m se g m
KgF KgF
( p+ γ h)|1=1 , 1 0
4
2
=10000 2
m m
Con los datos anteriores
Δ( p+ γ h) d ( p+γ h) 10000−17750
= = =−1825
Δl dl 3 √2
KgF
21825
m 3 también como datos tenemos:

a 5 10  3 m
m
V  1
seg
Los cuales reemplazados en la ecuación (3.2.2) queda:
−y 1825 2
u= −3
+ (0,005 y− y )
5.1 0 (2 , 08/98)
KgF .seg
( )
En la cual se hizo la conversión de poise a m2 : como poise es la unidad “cgs” de viscosidad:
2 KgF .seg
dyn . seg 1 KgF cm 0,8 ( )
1 poise= 1 2
⋅ ( )⋅10000 ( 2 )= m2
cm 980000 dyn m 98

Operando sobre el resultado final de “u” queda para la ecuación de distribución de velocidades:
u 359y  111800 y 2

La velocidad máxima se presenta donde se obtienen los máximos y mínimos de la función u (y), o
sea para
du
0
dy
−3
259−2⋅111800 y=0→ y=359 /(2 ⋅111800)=1 ,60 ⋅1 0 m=1 , 6 mm

Las velocidades mínimas, es decir donde el líquido tiene un punto de estagnación, las obtenemos
para:
y=0  u=0
y=5mm  u=359⋅ 0,005−111800 ⋅¿

y la velocidad máxima:
y=5mm  u=359⋅ 0,0016−111800 ⋅¿
Y el perfil de velocidades sería:
Fig.3.2.4

b.-Cálculo del caudal:

Como el diferencial de área es dA = dy . 1
a m2
Va 1 d [ ]
Q=∫ udy = −
3
( p+ γ h)a seg
0 2 12 μ dl

2 2 3
m m ⋅m m 1
Pero la unidad [ ]= =[ ] es decir caudal por unidad de dimensión en la
seg seg ⋅m seg m
3
m 1
profundidad de la placa. Reemplazando valores Q=−0 ,17 ⋅1 0−3 [ ] da en un sentido negativo,
seg m
o sea hacia arriba, es decir que es mayor el flujo ascendente que el descendente.

c.- Cálculo de la tensión de corte:

Para calcular la tensión de corte que en definitiva nos permitiría calcular la fuerza necesaria para
mover la placa corriente arriba, calculamos primero (du/dy) o sea el gradiente de velocidad para esa
posición:

(
du
|
)
dy y=0,005
=(
0,8
98
)(−759)
kgf . seg 1
m 2
.
seg
kgf
=−6 , 19 2
m
Como da negativo la tensión para mover la placa debe ser en sentido negativo o sea oponiendo al
movimiento de la placa. Observe que para calcular la fuerza total deberíamos multiplicar esta
tensión de corte por el área de la placa y agregar la componente de peso de la placa según la
dirección del movimiento.

Análisis de Flujo Laminar Bidimensional:

Consideremos la condición de flujo incompresible, permanente y viscoso (que produce la
condición de flujo laminar) desplazándose por un tubo de sección circular inclinado, de radio interior
r = a ,como se muestra en la figura 4.4.1 siguiente, que es la condición más general posible.

Sobre el elemento másico de fluido de forma anular, concéntrico con el tubo , tomamos todas las
fuerzas de campo y superficiales que nos aparecen.
Aplicando la ecuación de movimiento en la dirección longitudinal l con aceleración nula tenemos 

∑ Fi=0 (Movimiento.Uniforme) 
dp d
2 πrδr p−(2 π rδr p+ 2 π rδr ( )δl)+2 πrδl τ−[2 πr δl τ + (2 πr δl τ )δr ]+ γ 2 πr δr δl senθ=0
dl dr
−d h
Como =sen θ ; sustituyendo y dividiendo por el volumen del cuerpo libre (2 πrδrδl ), queda:
dl
dp d d d 1 d (4.4.1)
−(2 πrδr δl) − (2 πrδlδr )τ−γ (2 πr δr δl)=0→ ( p+ γ h)+ (τ ⋅r )=0
dl dr dl dl r dr
Multiplicando ambos miembros por (r . dr ) queda:

{ dld ( p+ γ h)}r ⋅dr + d (τ r )=0

Fig.4.4.1

Como la llave no es una

función de r, y τ es constante para el tubo de corriente r ,entonces integrando queda:
2
d r d
( p+ γ h)∫ r ⋅ dr +τ ∫ dr =0 → (p +γ h)+ τ r +C 1=0
dl 2 dl

Como esta ecuación debe ser válida para todo el campo del flujo, aun cuando r = 0  C1 = 0 
r d
( p+ γ h)+τ =0 (simplifico una r)
2 dl

du
Sustituyendo ahora τ por su valor para flujo viscoso τ =μ( )
dr
r d
2 dl
du −1 d
( p+ γ h)+ μ( )= r ( p+ γ h)→ du=
dn 2 μ dl
−1 d
2 μ dl {
( p+γ h) r dr → }
{ }
2
−1 d
u=
2 μ dl
( p+ γ h) ∫ rdr ⟹ u=−r d
4 μ dl
( p+ γ h)+ C2 (5.4.2)
Aplicamos otra condición de borde para calcular la constante C2.
2
a d
Cuando r = a  u = 0  C 2= ( p+γ h) queda la expresión general de la velocidad del campo
4 μ dl
a2 −r 2 d
u=( )− ( p+γh )
fluido: 4μ dl (4.4.3)

y para la velocidad máxima en r = 0:

a2 d
umx = ( p+γh )
4 μ dl
y para tubo horizontal:
2 2
a −r Δp
d h=0→ h=0 →u=( )( )
4μ Δl
siendo el valor de la velocidad máxima:
2
a Δp
Umx= ( )
4 μ Δl

Δp / Δl es la razón de la presión entre dos puntos del tubo entre los cuales calculamos la velocidad,
o mejor el perfil de velocidades. La representación de la velocidad “u” en función del radio “r” como
se ve en la ecuación es claramente una parábola para cada sección longitudinal y un “paraboloide
de revolución” para la distribución bidimensional de velocidades en el tubo.

Fig. 4.4.2

Como la distribución de velocidades es un paraboloide, es fácil demostrar que la velocidad media

es exactamente la mitad de la velocidad máxima o sea:
2
a d
V= ( p+γ h)
8 μ dl
y, por lo tanto, el caudal es: V ⋅ π a2=Q
O sea:
4 4
πa d D πD d
Q= ( p+ γ h), Como a= →Q= ( p+ γ h)
8 μ dl 2 128 μ dl
Caída de presión para flujo laminar en tubos, fórmula de Hagen-Poiseuille:
A partir de la fórmula de caudal para tubo horizontal, la cual se mostró en análisis dimensional:
4
π D Δp
Q= ( )
128 μ Δl
y como en este caso:
Δp ΔP
( )=
Δl L
Siendo L la distancia entre la ubicación de dos manómetros o dos puntos cualesquiera, queda:
4
π D ΔP 128 μ LQ
Q= → Δ P= 4
128 μ L πD
Siendo:
2
πD
Q= ⋅V
4

Por lo tanto, vemos que la pérdida de presión (energía piezométrica) es directamente proporcional
a la viscosidad, a la longitud y al caudal, e inversamente proporcional a la cuarta potencia del
diámetro.
Obsérvese que la rugosidad del tubo no interviene en la ecuación, simplemente porque la velocidad
es cero en los bordes.
128 μ LQ Δ P 128 μ LQ
Δ P= ⟹ =
πD
4
γ ρgπD
4

Llamaremos pérdida de carga a

2
128 μ LQ 128 μ LV π ( D ) 32 μ L V
H L= 4
= 4
= 2
ρ g π ( D) ρ g π ( D) 4 ρ g ( D)

Esta ecuación es conocida como ecuación Hagen Poiseuille, y fueron obtenidas entre 1839
y 1840. Las ecuaciones son válidas para flujo ya establecido en cañería.
Para el tramo inicial de la cañería, si proviene de un depósito con una entrada abocinada, la
velocidad es casi uniforme al valor de la velocidad media. Después, debido a la ecuación de
continuidad, la velocidad aumenta sobre la línea central y disminuye hacia los bordes.
Cuando estudiamos las propiedades de la capa límite, vemos que dentro del tubo el espesor
de la misma va aumentando desde la entrada del depósito para cada sección longitudinal y en
poco trecho la capa invade todo el volumen del tubo desde la periferia interior de velocidad cero
hasta la línea central de velocidad máxima (recordar lo visto en capa límite).
Se puede establecer analíticamente que la longitud de transición L = 0,058 Re.D y está de
acuerdo a las observaciones experimentales. Podremos estar seguros del establecimiento de este
patrón de flujo laminar como ya dijimos, si el número de Reynolds Re es inferior a 2.300, es
decir, si se cumple:
 ρV D
ℜ¿ <2.300
μ
lo que indica flujo laminar desarrollado.

Otra expresión de la fórmula de Hagen-Poiseville:

ΔP d ΔP d p
Como ( )= (p+ γ h) →( )=γ ( +h)
Δl dl Δl dl γ
d p
Al valor ( +h) se lo llama “caída de presión por unidad de longitud”, se lo indica con la letra “ j”.
dl γ

Entonces: (
ΔP
ΔL
)=γ ⋅ j ⟹ j=( )
ΔP
γ ΔL
=
L [
H L m de columna de líquido
metros de cañería ]
De lo cual también se pueden escribir las ecuaciones anteriores como:
4 2
πD D 32 μ
Q= (γ j); V = (γ j)→ j= 2
V
128 μ 32 μ γD
ΔP
Asimismo a partir de la expresión: ( )=γ j → ΔP=γ j ΔL
ΔL
Y el valor de la fuerza resistente en la sección será:
2 2
πD πD
W= ( ΔP)= γ j ΔL
4 4
Y la resistencia en el contorno de la tubería será:
2
πD
γ j ΔL
4 γDj
τC= =
π D ΔL 4
Tensión tangencial superficial en el contorno interior de la tubería.

4.- Ejemplo de aplicación:

Determinar las características del flujo en el tubo representado en la figura 8.1 siguiente, siendo:

3
γ =800 kg f /m μ=0 , 40 poises

Calcular el caudal en el [litros / min.] y el

n° Re del flujo.

Fig.4.6.1
En la posición 1:
4 kg f kg f kg f
p1 +γ h1=1 ,5 ⋅1 0 2
+ 800 3 ⋅5 m=19.000 2
m m m
En la posición 2:
4 kg f kg f
p2 +γ h2=2 ⋅1 0 2
=20.000 2
m m
Como p2 > p1 suponemos en principio que el flujo se mueve de 2 a 1.
La expresión:
Δ 19.000−20.000 −1.000 kg f
( p+ γ h)= = =−100 3
ΔL √ 5 + 8.6 6
2 2 10 m
que reemplazada en la expresión del caudal:
4 4 4 3
πa Δ πD Δ 3 , 14 ⋅0 , 0 2 −5 m
Q= (p +γ h)= ( p+ γ h)= ⋅(−100)=−9 , 6 ⋅10
8 μ ΔL 128 μ ΔL 128 ⋅10 , 40/98 seg

En litros por minuto, el caudal será:

3
−5 m 3 lit seg
Q=9 ,6 ⋅1 0 ⋅1 0 2 ⋅ 60
seg m l Como el caudal tiene signo negativo será hacia arriba.
min
min
El N° Re será:
V Dρ 0,306 ⋅0 , 02 ⋅800 ⋅98
R= = =122 , 4 ≪2.300
μ 9 , 81 ⋅0 , 40
Este valor obtenido, significa que el régimen es laminar y las fórmulas están bien aplicadas, ahora,
en la expresión anterior tomamos:
2
Q/ π D γ
V= y ρ=
4 g
y el valor superficial de la tensión:
γ D j γ D 32 μ V 8 μ V 8⋅(0 , 40 /98)⋅0,306 kg f
τ= = 2
= = ≅ 0,5 2
4 4 γD D 0 , 02 m

5.- Un eje rotatorio a 300 RPM y diámetro 30 mm está montado a un núcleo metálico sólido en una
máquina a través de un cojinete de latón que está solidarizado al núcleo y en contacto giratorio
suelto con el eje. En el huelgo “h” entre ambos se coloca una grasa lubricante con las siguientes
características a la temperatura estable de funcionamiento.
ρ = 888 Kg/m3
μ = 0.837 Kg/m s
Se pide determinar cual deberá ser el huelgo óptimo de diseño entre el eje y el cojinete para que
dentro del huelgo exista un régimen laminar. Utilizar un Reynolds crítico de 2100
6.- Para un canal abierto a la atmósfera inclinado un ángulo θ respecto al terreno y de ancho b, que
trasporta una película de altura a de un fluido en régimen laminar, se sabe que la velocidad
γa γ
genérica de las láminas para cada cota viene dada por: Vy = y sin θ− y 2 sinθ
2μ 2μ

Siendo: y la cota genérica desde la base en [m], γ el peso específico de la sustancia transportada.
Se pide encontrar el valor que tendrá la velocidad media de la lámina y el caudal volumétrico
transportado por el canal por unidad de ancho y el total, utilizando lo visto en los ejemplos
anteriores.

7.- Un aceite lubricante medio, de densidad relativa 0,86 es bombeado a través de una tubería
horizontal de 5 cm de diámetro y 300 m de longitud. El caudal bombeado es de 1,20 l/seg. Si la
caída de presión es de 2,10 kg/cm2, ¿cuál es la viscosidad absoluta del aceite?
4
π D ΔP 128 μ LQ
Q= → Δ P= 4
128 μ L πD

8.- Un caudal de 44 l/seg de un aceite de viscosidad absoluta 0.0103 kg seg/m 2 y densidad relativa
0,850 está circulando por una tubería de 30 cm de diámetro y 3000 m de longitud. ¿Cuál es la
pérdida de carga en la tubería?

9.- Del punto 1 al punto 2, de una tubería está fluyendo un fuel-oil pesado; la tubería es de acero y
se encuentra horizontal. La longitud de la tubería es de 900 m y tiene 15 cm de diámetro. La
presión en A es de 11,0 kg/cm2 y en B de 0,35 kg/cm2. La viscosidad cinemática es 4,13 x 10 -4
m2/seg y la densidad relativa 0,918. ¿Cuál es el caudal en l/seg que circula por esa tubería?
Resuelva el problema considerando como primer hipótesis que el flujo es laminar y luego
compruebe.
10.- Como se muestra en la figura, un tubo redondo de 0,5 mm de diámetro y 750 mm de longitud
está conectado al pleno de un ventilador. El ventilador produce una presión relativa de – 3,81 cm de
H20 en la cámara impelente y toma aire (20 ° C) en el micro canal. ¿Cuál es la velocidad media del
aire en el micro canal? Suponga que la única pérdida de carga está en el tubo del micro canal.
Datos del aire a 20°C, densidad 1,2047 Kg/m 3, viscosidad cinemática 1,488 10-5 m2/s, viscosidad
absoluta 18,288 10-7 Kg seg/m2.

Adenda de fórmulas para flujo laminar:

Con referencia a la figura (lado derecho abajo) podemos tener diferentes expresiones para las
variables solicitadas
A.- Flujo Unidimensional:

A1.- Placa oblicua a velocidad U, flujo en movimiento por gravedad:

u=
Uy 1 Δ
a
− {
2 μ Δl }
( p+ γh) (ay − y 2 )

du
τ =μ( )|y =a
dy
u=velocidad gnerica en cualquier valor de r
A2.- Horizontal con gradiente presión según la longitud:

Uy 1 Δ
u= − ( p)(ay− y 2 )
a 2 μ Δl
A3.- Horizontal sin gradiente de presión:

Uy
u =
a

B.-Flujo Bidimensional

B1.- Tubo inclinado y gradiente de presión:

Velocidad genérica:

a 2−r 2 Δ
u=( ) ( p+γh)
4μ Δl
2
a Δ
umx= ( p+γh)
4 μ Δl
Velocidad media:

a2 Δ
V= ( p+γ h )
8 μ Δl
Caudal:

πa4 Δ πD 4 Δ
Q= ( p+γh)= ( p+γh)
8 μ Δp 128 μ Δl

b. Tubo horizontal:
2 2 2 2
dh
=0 →h=0 velocidad gen é rica ⟹ u= a −r ( Δp )velocidad media ⟹ V = a ( Δp ) V = a ( Δp )
dl 4μ Δl 8 μ Δl 8 μ Δl
4 4
π a Δp π D Δp
caudal ⟹ Q= ( )= ( )
8 μ Δl 128 μ Δl
2
πD
Δp
4 DΔp
tensi ó n de corte en la superficie de la ca ñ er í a ⟹ τ= =
πD 4