Capitulo 1 Balance Momento

Descargar como ppt, pdf o txt
Descargar como ppt, pdf o txt
Está en la página 1de 37

UNIDAD 1 TRANSFERENCIA DE MOMENTO

El estudio del flujo de fluidos en un ducto o una conducción


requiere considerar:

- La geometría del ducto o de la


conducción que contiene al fluido
(ejemplo: longitud, diámetro,
inclinación , etc.)
- La magnitud de la fuerza que causa
el flujo
- Las propiedades físicas del fluido
El esfuerzo de corte () o de
cizalle entre dos capas
contiguos de fluido que se
desplazan a velocidades
diferentes representa la
resistencia que opone el
fluido al movimiento.
La unidad tradicional para la viscosidad es el poise (P)

1 poise = 1(g /cm∙seg) = 1P

Para líquidos y gases se prefiere el centi poise o 1cp.

Ejemplo:

H2O pura a 20° C = 1 cp

Aceite motor 20°C = 8 poise

Escoria líquida = 14 poise


dv x
Si graficamos yx en función de 
dy
tendremos:

 Fluidos Newtonianos que


responden en forma lineal,
pasando por el origen del
sistema de coordenadas y cuya
pendiente corresponde a la
viscosidad del fluido para una
temperatura y presión dada.

 Fluidos no newtonianos que


no se comportan según la
ecuación del esfuerzo de corte.
El estudio de estos fluidos
constituye una ciencia amplia
sobre la deformación de estos,
estudios Reológicos.
Líquido Viscosidad (kg/ms)
Alcohol etílico 1,2 x 10 -3
Gasolina 2,9 x 10 -4
Mercurio 1,5 x 10 -3
Aceite lubricante 0,26
Agua 1,005 x 10 -3

Viscosidad y temperatura
Conservación de momento en flujos en estado estacionario

Para efectuar un balance de momentum al flujo de un fluido, que


permanece en régimen estacionario, debemos considerar un volumen de
control (V.C) escogido convenientemente sobre el cual realizar el
balance. Para tal efecto debemos considerar dos tipos de transporte de
momento:

1. Transporte de momento debido a la existencia de un gradiente de


velocidad en la dirección normal a la dirección del flujo, llamado
Transporte de Momento Viscoso.

2. Transporte de momento debido al movimiento global del fluido en la


dirección de este, llamado Transporte de Momento Convectivo
Considere flujo laminar y estacionario de un fluido newtoniano incompresible
entre dos placas paralelas infinitas. La placa superior se mueve con velocidad
V y la placa inferior está en reposo. La distancia entre estas placas es Y y la
gravedad actúa en dirección z negativa. No hay otra presión en el sistema.

Y
Se define el V.C. de dimensiones x y z, fijo en el espacio con relación a la
placa inferior

Balance de momentun en estado estacionario

velocidad de entrada  velocidad de salida  suma de las 


     
 de 
  de 
 fuerzas que actúan 0
cantidad de momentun cantidad de momentun sobre el sistema 
     
Esfuerzo de corte
dv x

dy
Dimensionalm
ente:

Fuerza Masa  Longitud 1 Masa longitud 1


   x = x x
2 tiempo tiempo Area
Area Tiempo Area

masa  velocidad 1
= x 
tiempo Area

momento
=
tiempo x area
Donde
 = velocidad de transporte de momento por área unitaria
(densidad de flujo de cantidad de movimiento)
velocidad de transporte de momento
convectivo:

masa velocidad kg m m3 kg m m 2 2
  m  v A
tiempo s s m3 m3 s s

ρ = densidad
A = Área transversal donde ocurre el flujo
de fluido.
El transporte de momento
convectivo debido al flujo en la
dirección x, entra al V.C. por la
cara yz en x y sale del V.C. por
la cara yz en x + x.

El transporte de momento
viscoso va desde regiones de
mayor momento hacia las
regiones de menor momento;
luego el transporte entra por la
cara xz en y + y y sale por la
cara xz en y

En el V.C. la única fuerza que


actúa es la fuerza de gravedad
sobre la masa del fluido. Esta
fuerza actúa en la dirección y,
normal a la dirección del fluido,
no influye en el flujo de fluido.
2
A: momento convectivo que entra en la cara yz en x a la velocidad  vx y z
x
2
B: momento convectivo que sale en la cara yz en x+x a la velocidad  v x y z
x  x

C: momento viscoso que entra en la cara xz en y+y a la velocidad  yx x z


y  y

D: momento viscoso que entra en la cara xz en y a la velocidad  yx x z


y

2 2
 v x x y z   v x y z   yx x z   yx x z  0
x  x y  y y
En régimen estacionario la velocidad es la misma en todos los puntos. (Notar
que vx solamente depende de y, no depende de x y del tiempo)

v x | x  v x | x  x

(  yx | y  y   yx | y ) x z  0 |: Δx.Δy.Δz

Resolviendo y tomando el límite a esta expresión cuando Δy → 0


tenemos:
d yx
0
dy

Integrando se tiene τyx = C1


dv x
Reemplazando τyx por la ecuación de Newton:   C1
dy
C1
Integrando nos queda: vx   y  C2
μ

Las constantes C1 y C2 son evaluadas según las siguientes condiciones de


borde:
CB 1: vx = 0 si y = 0;
CB 2: vx = V si y = Y;
V
C1  
De donde C2 = 0 y Y

El perfil de velocidades está dado por V


vx  y
Y
Notar que:

 vx es proporcional al gradiente de velocidad e independiente de la


viscosidad del fluido;

 vx es una función lineal de y.

Reemplazando el valor de C1 , podemos calcular yx

V
 yx   
Y

Luego, yx es independiente de y; pero es proporcional a la viscosidad del


fluido
El procedimiento estándar para calcular el esfuerzo cortante y la distribución de
velocidades a partir de un balance de momento es el siguiente.

 Aplicar un balance de momento a una delgada envoltura de fluido


(volumen de control) y reducido a una ecuación diferencial que contenga a
.
 Integrar la ecuación diferencial obtenida.

 Obtener una expresión para  y sustituir la ecuación de Newton de la


viscosidad, integrar una segunda vez para obtener una ecuación para la
velocidad.

 Evaluar las dos constantes de integración a partir de dos condiciones de


Borde.
En las integraciones mencionadas aparecen constantes de integración que son
evaluadas considerando condiciones límites. Es decir, determinaciones de
hechos físicos para valores bien concretos de la variable independiente.

Algunas situaciones límites:

 En las interfases Sólido-Fluido, la velocidad del fluido es igual a la velocidad


con que se mueve la superficie misma. (El fluido se considera adherido a la
superficie de la placa).

 Interfase Líquido-Gas, la viscosidad de los gases es muy pequeña comparada


con la viscosidad de los líquidos. Cuando un líquido se desplaza, la fase gaseosa
de la interfase no ofrece resistencia a este desplazamiento, por lo tanto la
densidad de flujo de cantidad de momento y el gradiente de velocidad en la fase
líquida es muy pequeño, lo que significa que en la mayoría de los cálculos puede
suponerse igual a cero.
 En las interfases líquido-líquido, tanto la densidad de flujo de cantidad de
movimiento como la velocidad son continuas a través de la interfase; es
decir, que son iguales a ambos lados de la interfase.

La resolución de este balance permite calcular : velocidad media,


velocidad máxima, velocidad volumétrica de flujo, pérdida de presión
y fuerzas que actúan sobre las superficies límites.
1.1 Flujo de Fluido entre dos Placas Planas Paralelas

Considere flujo laminar y estacionario de un fluido newtoniano


incompresible, contenido entre dos placas planas paralelas
infinitas, sometido a un gradiente de presión en la dirección x
Consideraciones previas:

 Este flujo requiere una fuerza externa que impulsa al fluido en la dirección
de las placas planas paralelas.

La velocidad del fluido en contacto con las placas, es cero.

vx | y=  =0

 La velocidad del fluido aumenta hacia la línea central del fluido.


dv x
 La máxima velocidad ocurre en el plano medio de las placas ( dy  0 ), por
tanto, según la ecuación de Newton, el esfuerzo de corte yx es igual a cero.

 Estas consideraciones se utilizan como condiciones de borde en el balance


de cantidad de movimiento, o balance de momento, aplicado a un volumen
de control
Aspectos sobre la simetría del sistema:

 Es conveniente que y = 0 este ubicado en el plano central de las placas


que están separadas a una distancia 2.

 El V.C. puede localizarse en cualquier lado de y = 0 y lo suficientemente


alejado de la entrada del flujo a objeto que se cumple la condición de
Régimen Estacionario.

 Notar que en Régimen Estacionario la velocidad local o instantánea es


independiente de la posición en la dirección del flujo.

 La cantidad de movimiento viscoso se transporta desde y = 0 (máxima


velocidad) hacia las dos placas.
Consideramos un Volumen de Control de dimensiones x; y; z y
apliquemos un balance de cantidad de movimiento:
Velocidad con el cual el momento
convectivo ingresa al V.C.  v 2x |x  y  z 
Velocidad con el cual el momento  v 2x x  x  y  z 
convectivo sale del V.C.
Velocidad con el cual el momento viscoso  yx | y  x  z 
ingresa al V.C.
Velocidad con el cual el momento viscoso  yx |y  y  x  z 
sale del V.C.
Suma de fuerzas que actúan sobre el fluido (Px|x – Px|x + x )(y z)
en el V.C. (Caída de presión en la dirección
del flujo)
Fuerza de gravedad El flujo de fluido es horizontal,
no influye la fuerza de
gravedad.
 v 2
x x  v 2x x  x  y z   yx y   yx
y  y
 x z   P x

 P x  x y z  0

Si consideramos que el flujo se encuentra en régimen estacionario, es decir,


la velocidad no depende de la posición en el eje x, depende solamente de
su posición en el eje y, tenemos:

V x| x = V x| x+x

 yx y   yx  x z   P  P
y  y x x  x
 y z  0
      P  P
yx y yx y  y
x x  x
 0
y x
  yx y   yx y  y

lim   d yx
y  0     dy
 y 
 

P x  P x  x dPx
lim 
x 0 x dx

d yx dPx dP P
  0 pero    = inicial – final
dy dx dx L

d yx P

dy L
P
La primera integral da:  yx  y  C1
L

Condición de Borde 1: Para y = 0 yx = 0, luego C1 = 0;

P
 yx  y
L

Según esta ecuación,  yx es proporcional a la caída de presión por unidad


de longitud a través del flujo, además es una función lineal de y, con un valor
máximo de
P  en la placa.
L
Según Ley de Newton para la viscosidad yx = -  dv x y la ecuación queda como
dy

dv x P
  y
dy L

P 2
Integrando vx   y  C2
2 L

Condición de Borde 2 si y = δ entonces vx = 0 luego

P 2 2
vx  (  y )
2L
P 2
a) Velocidad máxima (y = 0) v x max  
2L

Distribución de velocidades y esfuerzo de corte para un fluido entre


placas planas paralelas
b) El componente de la fuerza F del fluido sobre la superficie se obtiene
mediante el esfuerzo de corte que actúa sobre la placa debido al flujo de un
fluido viscoso, el cual a su vez es causado por la aplicación de una presión
normal al desplazamiento del flujo. Para placas de ancho W y largo L, la fuerza
que actúa sobre la placa es:
F = (Esfuerzo de corte) x
(Área) = ( y x| y = ) (2 L W)

P
=    2 LW 
L

= P  2  W
c) Velocidad media del fluido vx
 v x dA
vx 
 dA
W
  v x dA
vx  0 0
W
  dA
0 0

P  2 
v      y 2  dy
x 2 L  
0

P 2
vx  
3  L
d) Velocidad de flujo volumétrico Q  V 
Q = (velocidad media) x (área de la sección normal al flujo)
2 P 3
Q  vx  2  W   W
3 L

e) Velocidad de flujo Másico M

M = (velocidad de flujo volumétrico) x (densidad del fluido)

2 P 3
M
  W
3 L
1.2 Flujo de fluido de una película descendente
Consideremos una superficie plana inclinada, tal como se indica en la
figura. Se supone que la viscosidad y densidad del fluido son constantes. Se
considera una región de longitud L, suficientemente alejada de los extremos
de la pared, de forma que las perturbaciones de la entrada y la salida no
están incluidas en L.
1.3 Flujo de Fluido en un Cilindro Vertical
Consideremos el flujo laminar en estado estacionario, de un fluido de
densidad constante  en un tubo «muy largo» de longitud L y radio R.
1.4 Flujo de Fluidos en Tubos Concéntricos

También podría gustarte