UNIVERSIDAD NACIONAL DEL CENTRO DEL PERÚ

FACULTAD DE INGENIERÍA QUÍMICA

ECUACIONES DE VARIACIÓN PARA

SISTEMAS NO ISOTÉRMICAS

✓ Ecuación de continuidad.
✓ Ecuación de movimiento.
✓ Ecuación de energía

Dr. Orlando A. VILCA MORENO

 INTRODUCCIÓN

Estudiaremos las ecuaciones de continuidad, las que se obtienen del

Principio de la Conservación de la Masa aplicada al escurrimiento de
fluidos, a través de un “volumen de control”.
En efecto, considerando un volumen arbitrario, fijo en el espacio e
inmerso en un medio continuo en movimiento que lo ocupa en cada
punto y en todo instante es evidente que el balance entre la masa
entrante y saliente a través de la superficie del mismo y en un instante
dado, más la variación de la masa en su interior y con la variable
tiempo tendiendo a cero, da inexorablemente una masa resultante nula,
puesto que ésta no puede crearse ni desaparecer.
 OBJETIVOS

• Proporcionar conceptos básicos de ecuación de continuidad y su

aplicación a problemas de interés práctico.
• Desarrollar las ecuaciones de continuidad y cantidad de movimiento,
mediante la aplicación de la ley de la conservación de la materia a un
pequeño elemento de volumen situado en el seno del fluido en movimiento.
• Proporcionar las ecuaciones de continuidad y cantidad de movimiento
en los diferentes sistemas coordenados, con la finalidad de hacer más
sencillo su manejo.
 ECUACION DE LA CONTINUIDAD

La ecuación de continuidad se
obtiene aplicando un balance de
masa a un elemento diferencial de
volumen (V), a través de la cual
está circulando el fluido.

Figura 1. Región de volumen x y z fija en el espacio,

a través de la cual esta circulando el fluido
 Aplicando el balance de masa:
 velocidad de acumulación   velocidad de entrada   velocidad de salida 
 = −  (1)
 de masa   de masa   de masa 

Considerando el par
de caras
perpendiculares al eje
x, y y z, se obtiene que
la velocidad de
entrada de materia a
través de las caras x, y
y z, son:

xyz

t
  (
= yz (v x ) x − ( v x ) x + x + xz  v y

) y − (v y ) y + y 

+ xy (v z ) z − (v z ) z + z  (2)
Dividiendo toda la ecuación por x y z, y tomando el límite cuando x, y,
z tiende a cero, se tiene:

     
t
= −
 ( v x )+ (
v y )+ ( v z ) (3)
 x y z 

Ésta es la ecuación de continuidad y puede escribirse en una forma más

conveniente utilizando notación generalizada:

= −( •  v )
(4)
t

El término se denomina divergencia de v, en donde la divergencia () y la

velocidad (v) son vectores.
Se presentan a continuación la ecuación de continuidad en los demás sistemas
coordenados.
 ECUACION DE CONTINUIDAD EN DISTINTOS
SISTEMAS COORDENADOS

Coordenadas rectangulares (x, y, z):

   
t
+
x
(v x ) +
y
(
v y +
z
)
(v z ) = 0

Coordenadas cilíndricas (r, , z):

 1  1 
+ (rv r ) + (v  ) + (v z ) = 0
t r r r  z

Coordenadas esféricas (r, , ):

 1  2
(
+ 2 r v r +
t r r
)1 
r sen  
(v  sen ) + 1 
r sen  
(v  ) = 0
 ECUACIONES FUNDAMENTALES DE UN FLUJO

Son tres las ecuaciones fundamentales de un flujo:

✓ Ecuación de continuidad (conservación de masa)

✓ Ecuación de la energía (conservación de la energía )

✓ Ecuación de cantidad de movimiento (conservación de

cantidad de movimiento
 ECUACIONES FUNDAMENTALES DE UN FLUJO
ECUACIÓN DE CONTINUIDAD
En un flujo permanente, m 1 = m 2 = m
𝜌1 𝑣1 𝑆1 = 𝜌2 𝑣2 𝑆2 = 𝜌𝑣𝑆
1  Q1 =  2  Q2 =   Q
1. Por el siguiente conducto circula un caudal de 1250 m³ cada hora ,
conociendo que el diámetro de la sección “1” es de 8 cm y en la sección
“2” es de 12 cm . Calcular la velocidad del fluido en cada una de las
secciones. El fluido es agua.

Figura 1. Volumen de control para flujo a través de tuberías en serie.

 Por el principio de continuidad tenemos: 𝑄1 = 𝑄2 = 𝑄 1

𝑣1 𝑆1 = 𝑄 = 𝑣2 𝑆2 2
𝑄=𝑣𝑆
𝑚3
Datos: 𝑄 = 1250 ; 𝐷1 = 8 𝑐𝑚; 𝐷1 = 12 𝑐𝑚
ℎ

De la ecuación 2:

𝑚3 1ℎ
𝑄 1250 .
ℎ 3600 𝑠
𝑣1 = = = 69,08 𝑚/𝑠
𝑆1 𝜋 4𝑥10−2 𝑚 2

𝑣2 𝑆2 = 𝑣1 𝑆1

𝜋 2 2 2 2
𝑆1 𝐷1 𝐷1 𝑅1 𝑚 4 𝑚
𝑣2 = 𝑣1 4
= 𝑣1 𝜋 = 𝑣1 = 𝑣1 = 69,08 = 30,702
𝑆2 𝐷22 𝐷2 𝑅2 𝑠 6 𝑠
4
2. Por la tubería horizontal mostrada en la figura circula agua a 68 °F. El
diámetro de la sección “1” y “3” son de 20 cm, reduciéndose a la mitad
en la sección “2”. Si la velocidad en la sección “2” es de 12 m/s, calcular:
a. La velocidad en las secciones “1” y “3”
b. La velocidad de flujo volumétrico.

Solución:
𝑣1 = 3𝑚/𝑠 𝑄 = 0.09425 𝑚3 /𝑠

𝑣3 = 3𝑚/𝑠
3. Un jardinero usa una manguera para llenar una cubeta de 30 litros, el jardinero
observa que tarda 2 minutos en llenar la cubeta. A la manguera se le conecta
una boquilla con abertura de 0.5 cm² de área de sección transversal. Determine
la velocidad de salida del chorro de agua.
Haciendo uso de la ecuación de la continuidad:
LEY DE NEWTON:

La ecuación de movimiento, para un fluido Newtoniano:

Formas simplificadas de la ecuación de movimiento:

Fluidos de densidad y viscosidad constantes: (Ecuación de Navier - Stokes

Sistemas de flujo sin efectos viscosos: (Ecuación de Euler)

Fluido en reposo:
La ecuación de movimiento en coordenadas rectangulares:
En función de esfuerzo de corte
La ecuación de movimiento en coordenadas rectangulares:
Para fluidos Newtonianos de 𝝆 𝒚 𝜇 constantes
La ecuación de movimiento en coordenadas cilíndricas:
En función de esfuerzo de corte
La ecuación de movimiento en coordenadas cilíndricas:
Para fluidos Newtonianos de 𝝆 𝒚 𝜇 constantes
La ecuación de movimiento en coordenadas esféricas:
En función de esfuerzo de corte
La ecuación de movimiento en coordenadas esféricas:
Para fluidos Newtonianos de 𝝆 𝒚 𝜇 constantes
 Componentes del Tensor Esfuerzo Cortante:

Coordenadas rectangulares Coordenadas cilíndricas Coordenadas esféricas

 Ecuación de la Continuidad

Coordenadas rectangulares (x, y, z):

Coordenadas cilíndricas 𝒓, 𝜽, 𝒛 :

Coordenadas esféricas 𝒓, 𝜽, ϕ :
Ejemplo

Un Líquido newtoniano e incompresible se

encuentra retenido en la zona anular comprendida
entre 2 cilindros concéntricos (ver Figura 1). El
cilindro externo (radio 𝑅2 ), tiene un movimiento
oscilatorio simple dado por: 𝑣 = 𝑣𝑜𝑠𝑒𝑛(𝑤𝑡) . El
cilindro interno (radio = kR; k<1) esta quieto. Asumir
que el fenómeno se desarrolla en forma estable.
Partiendo de las ecuaciones de movimiento y
continuidad obtener el perfil de velocidad y el perfil
del esfuerzo cortante.

Solución
La ecuación de movimiento en coordenadas cilíndricas:
 Solución
𝜕𝑣𝑟 𝜕𝑣𝑟 𝜕𝑣𝑟
𝑣𝑟 ≠ 0 = 0; = 0; ≠0
ቐ𝑣𝜃 = 0 𝜕𝑧 𝜕𝜃 𝜕𝑟
1
𝑣𝑧 ≠ 0 𝜕𝑣𝑧 𝜕𝑣𝑧 𝜕𝑣𝑧
≠ 0; = 0; ≠0
𝜕𝑟 𝜕𝜃 𝜕𝑧

La ecuación de movimiento en coordenadas cilíndricas, componente r:

Ecuación de movimiento en coordenadas cilíndricas, componente z:

3
Solución
La ecuación de continuidad en coordenadas cilíndricas:

Si: 𝜕𝑣𝑧
𝐿≫𝑅 𝑣𝑟 ≈ 0 =0 5
𝜕𝑧

El perfil del esfuerzo cortante:

 Energía de un líquido en movimiento

𝑊𝐹1 + 𝐸𝐶1 + 𝐸𝑃1 = 𝑊𝐹2 + 𝐸𝐶2 + 𝐸𝑃2 = 𝐶𝑡𝑒. 1

1 1
𝐹1 ∆𝑥1 + 𝑚 𝑣12 + 𝑚𝑔ℎ1 = 𝐹2 ∆𝑥2 + 𝑚 𝑣22 + 𝑚𝑔ℎ2 2
2 2
Dividiendo la ecuación (2) entre su volumen: (V=S ∆𝑥)
𝐹1 ∆𝑥1 1 𝑚 𝑣12 𝑚𝑔ℎ1 𝐹2 ∆𝑥2 1 𝑚 𝑣22 𝑚𝑔ℎ2
+ + = + +
𝑆1 ∆𝑥1 2 𝑉1 𝑉1 𝑆2 ∆𝑥2 2 𝑉2 𝑉2

1 1
𝑃1 + 𝜌1 𝑣12 + 𝜌1 𝑔ℎ1 = 𝑃2 + 𝜌2 𝑣22 + 𝜌2 𝑔ℎ2
2 2

𝑃1 𝑣12 𝑃2 𝑣22
+ + 𝑔ℎ1 = + + 𝑔ℎ2 = constante 3
𝜌 2 𝜌 2

𝑃 𝑣2 𝑚2 𝑘𝑔 𝑘𝑔. 𝑚 𝑚 𝑁. 𝑚 𝐽
𝐸= + + 𝑔ℎ = 2 . = . = = 4
𝜌 2 𝑠 𝑘𝑔 𝑠 2 𝑘𝑔 𝑘𝑔 𝑘𝑔
Si ha la ecuación (3) lo dividimos entre su gravedad:

𝑃1 𝑣12 𝑃2 𝑣22
+ + 𝑔ℎ1 = + + 𝑔ℎ2 3
𝜌 2 𝜌 2

𝑃1 𝑣12 𝑃2 𝑣22
+ + ℎ1 = + + ℎ2
𝜌𝑔 2𝑔 𝜌𝑔 2𝑔

𝑃1 𝑣12 𝑃2 𝑣22
+ + ℎ1 = + + ℎ2 4
𝛾 2𝑔 𝛾 2𝑔

𝑃 𝑣2 𝑁 𝐽
ℎ= + + ℎ = 𝑚. =
𝛾 2𝑔 𝑁 𝑁

o en metros de columna de líquido,

Ala de un Avión

𝑷𝟏 < 𝑷𝟐
𝟏
𝑷𝟏 + 𝝆 𝒗𝟏 − 𝒗𝟐 = 𝑷𝟐
𝟐
Ala de un Avión

𝑷𝟏
𝒗𝟏

𝒗𝟏 > 𝒗𝟐 𝑷𝟏 < 𝑷𝟐 𝒗𝟏 > 𝒗𝟐 𝑷𝟐

𝑷𝟏 < 𝑷 𝟐

𝒗𝟐
𝑭
En la figura se representa el perfil del ala de un avión, con un determinado ángulo de ataque, dentro de una corriente
de aire laminar.
Si aplicamos de una manera cualitativa la ecuación de continuidad y el teorema de Bernoulli, se llega a las siguientes
conclusiones:
1. La máxima deformación de las líneas de corriente se produce en
la zona superior del borde de ataque, por lo tanto hay un
aumento de velocidad del fluido, consecuentemente, esto lleva
implícito una disminución de presión, muy marcada en el borde
de ataque, disminuyendo hacia el borde de fuga.
2. Justo por debajo del borde de ataque se aprecia una zona que no
hay líneas de corriente (4), la velocidad del fluido en esta zona es
nula, es la denominada zona de remanso. Por el teorema de
Bernoulli la presión aumentará en el borde de ataque,
encontrando una zona de sobrepresión, disminuyendo conforme
se entra en el perfil hacia el borde de fuga.
3. Finalmente por debajo del perfil y cerca del borde de fuga, se produce un pequeño aumento de la
velocidad y por lo tanto una pequeña de presión, que compensará en parte, la producida en la misma zona
por encima del perfil.
PRINCIPIO DE BERNOULLI

Es aquel efecto producido en los fluidos

dinámicos e ideales (viscosos), que explica el
comportamiento de un fluido moviéndose a lo
largo de un tubo de flujo.
El efecto Bernoulli describe, prácticamente, que
la presión de un fluido esta relacionada con la
velocidad de desplazamiento del mismo.
Este principio explica claramente el por qué lo
aeroplanos se pueden mantener en el aire. La
explicación a esto es sencilla.
Gracias a la forma aerodinámica de las alas de
estas maquinas, las líneas de flujo permanecen
juntas, formando en la parte inferior de las alas
una fuerza llamada “sustentación”, que consiste
en que la velocidad del flujo de aire por la parte
de arriba de la aleta va a ser mayor comparada
con la parte de abajo.
En colusión tenemos que a mayor velocidad,
menor presión.
Problema 1.
El agua con una densidad de 1 g/cm³ fluye a través de una tubería horizontal. La sección transversal
del tubo se estrecha de 80 cm² a 40 cm² en un reductor. La presión estática antes del reductor es
de 4 bar y la velocidad de flujo es de 4 m/s. El flujo es incompresible y sin fricción (no viscoso).
Determine la presión estática que se mide después del reductor.

Solución
Por la ecuación de Bernoulli

𝑃1 𝑣12 𝑃2 𝑣22
+ + 𝑔ℎ1 = + + 𝑔ℎ2
𝜌 2 𝜌 2

Datos:
𝑔 𝑘𝑔
𝜌=1 = 1000
𝑐𝑚3 𝑚3
𝑆1 = 80 𝑐𝑚2 ; 𝑆2 = 40 𝑐𝑚2
𝑃1 = 4 𝑏𝑎𝑟
𝑣1 = 4 𝑚/𝑠

𝑃2 𝑃1 𝑣12 − 𝑣22 𝑣12 − 𝑣22

= + 𝑃2 = 𝑃1 + 𝜌
𝜌 𝜌 2 2
 Por la ecuación de la continuidad para un fluido incompresible:

𝑄 = 𝑣1 𝑆1 = 𝑣2 𝑆2

𝑆1 𝑚 80
𝑣2 = 𝑣1 =4 = 8 𝑚/𝑠
𝑆2 𝑠 40

16 − 64 𝑚2 𝑘𝑔
𝑃2 = 4 𝑏𝑎𝑟 + 1000
2 𝑠2 𝑚3

1 𝑏𝑎𝑟
𝑃2 = 4 𝑏𝑎𝑟 − 24000 𝑃𝑎
100000 𝑃𝑎

𝑷𝟐 = 𝟑. 𝟕𝟔 𝒃𝒂𝒓

El aumento de la velocidad del flujo hace que la presión estática caiga de 4,00 bar a 3,76 bar. Esto puede
explicarse por el hecho de que parte de la energía asociada con la presión estática tuvo que utilizarse para
acelerar el agua. El aumento de la energía cinética del agua se produce a expensas de la presión estática.
 Problema 2

Una manguera con una sección transversal interna

de 1,24 cm² está conectada a un grifo de agua. La
manguera conduce a una altura de 6 metros sobre
el suelo, donde el agua fluye por una boquilla y se
recoge en una piscina. La piscina se llena con 30
litros por minuto. Un metro por encima del suelo,
se conecta un manómetro a la manguera para
medir la presión estática. El manómetro indica
una presión de 2 bar. La presión del aire ambiente
es de 1 bar. El flujo es incompresible. Determina la
velocidad que sale el agua por la boquilla.

Solución
Solución

𝑃1 𝑣12 𝑃2 𝑣22
+ + 𝑔ℎ1 = + + 𝑔ℎ2
𝜌 2 𝜌 2

𝑣22 − 𝑣12 𝑃1 − 𝑃2
= + 𝑔 ℎ1 − ℎ2
2 𝜌

𝑣22 − 𝑣12 𝑃1 − 𝑃2
= + 𝑔 ℎ1 − ℎ2
2 𝜌

𝑃1 − 𝑃2
𝑣2 = 𝑣12 + 2 − 2𝑔 ℎ2 − ℎ1
𝜌

𝑣2 =
AHORA ES TU TURNO
 PROBLEMA DE APLICACIÓN

El agua fluye a una velocidad de 200 litros por segundo hacia arriba a
través de una tubería vertical cónica. El diámetro en la parte inferior es
de 240 mm y en la parte superior de 200 mm y la longitud es de 5 m. La
presión en la parte inferior es de 8 bar y la presión en la parte superior es
de 7,3 bar. Determine
• La pérdida de carga a través de la tubería.
• Expresarlo sólo como una función de la cabeza de velocidad existente
 PROBLEMA DE APLICACIÓN
A través de una refinería, el etanol combustible fluye en una tubería
a una velocidad de 1 m/s y una presión de 101300 Pa. La refinería
necesita que el etanol esté a una presión de 2 atm en un nivel más
bajo. Determine la altura que debe tener la tubería para alcanzar
esta presión. Suponga que la velocidad no cambia.
Considere: La densidad del etanol es 789 𝑘𝑔/𝑚3