CAPÍTULO 2

Ecuación de Continuidad

Conservación de Momento

Universidad de Atacama - Departamento de Metalurgia

INTRODUCCION

En el Capítulo 1 se obtuvo la distribución del Esfuerzo de Corte , y la

distribución de velocidades para diferentes situaciones donde
ocurre flujo de fluido, considerando un balance de momentun sobre
un volumen de control escogido.

En este capítulo desarrollaremos Ecuaciones Generales, que podrán

ser aplicadas a cualquier situación.

Estas ecuaciones corresponden:

 Ecuación de continuidad.

 Ecuaciones de conservación del momento.

Universidad de Atacama - Departamento de Metalurgia

ECUACION DE CONTINUIDAD

El principio de la conservación de la masa es un enunciado de que la

masa no se puede crear ni destruir durante un proceso y que,
durante un análisis, debe tomarse en cuenta en su totalidad.

En un flujo estacionario, la cantidad de masa dentro del volumen de

control permanece constante, y como consecuencia, la conservación
de la masa se puede expresar como:

Razón de cambio del  Razón del flujo de  Razón del flujo de 

     
flujo de masa dentro   masa que entra al   masa que sale del 
del volumen de control  volumen de control  volumen de control 
     

Universidad de Atacama - Departamento de Metalurgia

Esta ecuación se deduce aplicando un balance de masa a un volumen
de control de dimensiones x y z, en estado estacionario, a través del
cual está circulando un fluido.

Razón del flujo de masa  Razón del flujo de masa 

 
  
que entra al volumen de control  que sale del volumen de control 

Volumen de Control en coordenadas cartesianas

Universidad de Atacama - Departamento de Metalurgia

La velocidad de masa a la entrada de V.C. a través de la cara x (y en forma
general para las otras caras)

 kg  m
 
densidad  3  X componente x velocidad   X area m 2
m  s
Esto es:
 v x x y z

Y la velocidad a la cual la masa sale del volumen de control a través de la

cara x + x es:

 vx x  x
y z

Universidad de Atacama - Departamento de Metalurgia

Velocidad a la cual la masa ingresa al V de C por  vx|x y z
la cara x en x
Velocidad a la cual la masa sale del V de C por la  vx|x+x y z
cara x en x+x
Velocidad a la cual la masa ingresa al V de C por  vy|y z x
la cara y en y
Velocidad a la cual la masa sale del V de C por la  vy|y+y z x
cara y en y+y
Velocidad a la cual la masa entra al V de C a  vz|z y x
través de la cara z en z
Velocidad a la cual la masa sale del V de C a  vz|z+z y x
través de la cara z en z+z

Universidad de Atacama - Departamento de Metalurgia

Para un V de C, de volumen constante, la velocidad de acumulación de masa se
verá reflejada en la variación de densidad al interior del volumen x y z, por
lo tanto reemplazando en el balance de masa tenemos:

∂
∂t

x y z   v x x y z   v y y x y   v z z x y 

-  v x x x y z   v y yy x y   v z zz x y 
Si dividimos por x y z y luego aplicando el límite cuando x  0 y  0
z  0, la ecuación de continuidad resulta ser:

   
 ( v x )  ( v y )  ( v z )
t x y z

Universidad de Atacama - Departamento de Metalurgia

Vectorialmente la ecuación de continuidad se puede expresar como:

   v 
t
Si  del fluido es constante, es decir, no existe acumulación de masa en el
volumen de control (estado estacionario)

∂ ∂ ∂ 
 ∂x x ∂y y ∂z z   0
v  v  v
 

_
v  0

Universidad de Atacama - Departamento de Metalurgia

 CONSERVACIÓN DEL MOMENTO

Este balance se obtiene al aplicar la segunda ley de Newton del

movimiento a un elemento diferencial del volumen de control, y se puede
enunciar como: la fuerza neta que actúa sobre el volumen de control es
igual a la masa multiplicada por la aceleración del elemento de fluido
dentro del volumen de control, lo cual también es igual a la razón neta
de la transferencia de la cantidad de movimiento de flujo hacia fuera del
volumen de control.

Fuerzas de cuerpo (fuerzas de

gravedad, eléctricas, térmicas)
Fuerzas que actúan sobre
el volumen de control
Fuerzas de superficie (fuerzas de
presión, esfuerzos de corte)

Universidad de Atacama - Departamento de Metalurgia

CONSERVACIÓN DEL MOMENTO

La diferencia entre la velocidad de entrada de momento al volumen de

control y la velocidad de salida de momento del volumen de control
debe ser igual a la acumulación de momento en el volumen de control.

Para un elemento de volumen de control de dimensiones x y z, se

puede escribir el siguiente balance de cantidad de momento.

velocidad de  velocidad  velocidad  suma de las 

acumulacio n  de entrada  de salida  fuerzas que 
       
    
  
  
 de cantidad   de cantidad   de cantidad   actúan sobre 
de momento  de momento  de momento  el sistema 

Universidad de Atacama - Departamento de Metalurgia

Consideremos el esfuerzo de corte para la componente x del flujo

(a) Dirección del transporte de momento debido a la componente x de la

velocidad (b) Dirección de las fuerzas viscosas debido al transporte de
momento

Universidad de Atacama - Departamento de Metalurgia

Consideremos el esfuerzo de corte para la componente x de flujo en
nuestro sistema en coordenadas cartesianas rectangulares

Notar que xx (transporte de momento viscoso en la dirección x por el flujo

de fluido en la dirección x) se origina cuando la velocidad del flujo en el
sentido de x varía, es decir, vx cambia, ya sea por aceleración o
desaceleración del flujo en la dirección de x.

Luego, todo cambio de velocidad implica la existencia de una fuerza que

actúa en un plano normal a la dirección de vx. Por lo tanto hablaremos de
un Esfuerzo Normal al plano perpendicular al eje x del volumen de
control, en el punto x.

xx x y z xx x  x y z   yx y x z 
yx y  y
 
xz  
 yx y  y



xz 

  xy  zx z  z xy 

zx z

Universidad de Atacama - Departamento de Metalurgia

Para evaluar el trasporte de momento convectivo, debemos considerar el
flujo convectivo de la cantidad de momento x a través de las seis caras;
esto es:

 v x v x x -  vx vx  y z   v x v y y -  v x v y yy x z
x  x

  v x v z z -  v x v z z  z  x y

En la mayoría de las situaciones reales, las únicas fuerzas importantes

serán las procedentes de la presión del fluido y la fuerza de gravedad por
unidad de masa. La resultante de estas fuerzas en la dirección x será:

P x 
 P x  x y z   g x x y z

Universidad de Atacama - Departamento de Metalurgia

Para el Balance de momento, la suma de los términos de las
ecuaciones descritas es igual a la acumulación de momento en la
dirección de x en el Volumen de control, esto es

  v x 
 x y z
t

 v x    xx x  x   xx x  yx y  y   yx y  zx z  z   zx z 
    
t  x y z 
 

 v x v x x  x  v x v x x v x v y y  y  v x v y y v x v z z  z  v x v z z 
   
  x  y z 
 

P x  x  P x 
   g x
 x 
 

Universidad de Atacama - Departamento de Metalurgia

 Tomando límite x  0 y  0 z  0, resulta

v x     yx zx   v x v x  v x v y  v x v z   P

  xx           g x
t  x y z   x y z  x

De manera similar, realizando un balance de momento en la dirección y, tenemos:

v y    xy  yy zy   v y v x  v y v y  v y v z   P

            g y
t  x y z   x y z  y

Y en la dirección z:

v z     yz zz   v z v x  v z v y  v z v z   P

  xz           g z
t  x y z   x y z  z

Universidad de Atacama - Departamento de Metalurgia

Combinando estas tres ecuaciones, se obtiene la ecuación para la
conservación de momento, en notación vectorial:


 v           v v   P   g
t

La ecuación de momento expresada de esta forma, establece que un

pequeño elemento de volumen que se mueve con el fluido es acelerado
por las fuerzas que actúan sobre él

Universidad de Atacama - Departamento de Metalurgia

Para obtener la distribución de velocidades debemos conocer la relación
entre el esfuerzo (tensión) y el gradiente de velocidad. Esto esta
expresado por la ecuación de Newton para la viscosidad, según:

Esfuerzo (tensión) Normal

 xx  2
dv x

2

  v 
x 3

v y
 v 
2
 yy  2 
y 3

v z
  v 
2
 zz  2
z 3

Universidad de Atacama - Departamento de Metalurgia

Esfuerzo cortante
 dv y v x 
 xy   yx    
 dx y 

 v v 
 xz   zx   z  x 
 x z 

 v z v y 
 yz  zy    
 y z 

Universidad de Atacama - Departamento de Metalurgia

ECUACIONES DE NAVIER STOKES PARA FLUIDOS CON DENSIDAD
Y VISCOSIDAD CONSTANTES

Para un fluido con densidad constante la ecuación de continuidad corresponde a:

v x v y v z
v    0
x y z

Y la componente x del término (  v v) puede ser escrita como:

 v x v x v y v x v z v x 
 v x  vy  vx  vy  vx  vz 
 x x y y z z 

 v x v x v x  v x v y v z 
  v x  vy  vz  vx    
x y z
 x y z
  

Universidad de Atacama - Departamento de Metalurgia

Pero el término entre paréntesis es igual a cero de acuerdo a la ecuación de la
continuidad

Si desarrollamos  según la dirección x, debemos considerar solo xx , xy y

xz , ya que los otros esfuerzos de corte no consideran al eje x

v x
 v 
2
 xx  2
x 3

 v v x 
 xy   yx   Y  
 x y 

 v Y v x 
 xz   zx    
 x z 

Universidad de Atacama - Departamento de Metalurgia

   v x    v y v x    v x 
    2   
       
x  x  y  x y  z  z 
 2vx 2vy 2vx  2
v  2
v x 

  2    z

 x 2 y x  y 2
 z x  z 2 
 
 2vx 2vx 2vx   v x v y v z 
         
 x y z x  x y z 
2 2 2

Universidad de Atacama - Departamento de Metalurgia

 El balance de momento según la dirección x para un fluido de densidad () y
viscosidad () constante es:

v x  v x v x v x    2 v x 2 2
 vx  vx 

P
     vx  vy  vz    2     gx
t  x y z   x
2
z  x
y 
o bien:

 v x v x v x v x    2vx 2
 vx
2
 vx  P
  vx  vy  vz     2 
2

2

   gx
 t x y z   x y z  x

Universidad de Atacama - Departamento de Metalurgia

Ecuación de continuidad en distintos sistemas coordenados

Universidad de Atacama - Departamento de Metalurgia

Universidad de Atacama - Departamento de Metalurgia
En función del gradiente de velocidad para un fluido de  y µ constante

Componente x

Componente y

Componente z

Universidad de Atacama - Departamento de Metalurgia

Universidad de Atacama - Departamento de Metalurgia
En función del gradiente de velocidad para un fluido de  y µ constante

Componente r

Componente 

Universidad de Atacama - Departamento de Metalurgia

Componente z

Universidad de Atacama - Departamento de Metalurgia

Universidad de Atacama - Departamento de Metalurgia
En función del gradiente de velocidad para un fluido de  y µ constante

Componente r

Componente 

Universidad de Atacama - Departamento de Metalurgia

Componente 

Universidad de Atacama - Departamento de Metalurgia

Componentes del Esfuerzo de corte en coordenadas cartesianas

Universidad de Atacama - Departamento de Metalurgia

Universidad de Atacama - Departamento de Metalurgia
Componentes del Esfuerzo de corte en coordenadas cilíndricas

Universidad de Atacama - Departamento de Metalurgia

Universidad de Atacama - Departamento de Metalurgia
Componentes del Esfuerzo de corte en coordenadas esféricas

Universidad de Atacama - Departamento de Metalurgia

Universidad de Atacama - Departamento de Metalurgia
Ejemplo

Considere el flujo de aceite lubricante comprendido entre un cilindro

estacionario y una varilla que se mueve a una velocidad de 0,5 m/s. La
varilla y el cilindro son coaxiales

a) Hallar la distribución de velocidad

b) Evaluar la velocidad en r = 0,6 cm

= 3000 kg/m3

 = 10 kg/m-s

Universidad de Atacama - Departamento de Metalurgia