Ecuacion de Continuidad

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 50

Unidad 4: Dinámica y teoremas

de conservación de la dinámica.
Ing. Nahuel Castello
Mecánica de los fluidos - Departamento de Ingenieria Industrial
Universidad Tecnológica Nacional FRH

2017

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Contenido de La Unidad
• Ecuación de continuidad (Modelos de Euler y de Lagrange)
• Deducción de la Ecuación de Bernoulli para flujo estacionario incompresible
• Equilibrio de fuerzas en traslación
• Teorema de conservación de la cantidad de movimiento
• Ecuación de Cauchy
• Ecuaciones Indefinidas
• Ecuación de Euler
• Equilibrio de momentos en rotación
• Ecuaciones de Navier Stokes
• Ecuaciones de Lamé para fluidos Newtonianos
• Ecuación de Euler y de la estática
• Integración y sistema de Bernoulli para flujo compresible
• Casos particulares del sistema de Bernoulli
• Diagrama Piezométrico
• Corrientes relativas (Bomba centrífuga)
Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuación de Continuidad
• La cantidad de masa de fluido que
fluye por una determinada superficie
por unidad de tiempo se conoce como
Caudal Másico, y se nombra como m
• El punto sobre la m se utiliza para
indicar que corresponde a la masa por
unidad de tiempo.
• El caudal másico a través de un área es
obtenido mediante la siguiente
integración:

m    m   Vn dAc
Ac Ac
m
• Aunque esta expresión es exacta, no
siempre es conveniente para análisis
de Ingenieria.

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuación de Continuidad
• La integral puede reemplazarse en muchos
flujos por la siguiente expresión haciendo uso
de la velocidad y densidad promedio.
1
Ac Ac
Vavg  Vn dAc

m  Vavg Ac

• El caudal volumétrico, V , se define como:

V   V dA
Ac
n c  Vavg Ac  VAc

• Nota: Muchas veces Q es utilizado para


denotar al caudal volumétrico.
• El caudal másico y volumétrico están
relacionados por la siguiente expresión:
m  V
Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuación de Continuidad
• El Principio de Conservación de
la Masa puede expresarse
como:

dmCV
min  mout 
dt
• Donde min y mout son el caudal
másico entrante y saliente en el
volumen de control, y dmCV/dt
es la relación de cambio de
masa dentro del volumen de
control.

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuación de Continuidad
• Para un volumen de control, CV, de
forma arbitraria.
– La tasa de cambio de masa dentro del VC
es:
dmCV d
dt
 
dt CV
 dV

– El flujo neto de masa a través del VC es:


mnet    m   Vn dA 
CS CS
  V n  dA
CS

• Entonces la expresión general de


P.C.M esta dada por:
d

dt CV
 dV    V n dA  0 
CS

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuación de Continuidad
Análisis Diferencial a través del uso de un
volumen de control de tamaño infinitesimal

Teorema de la Divergencia (T. Gauss)

Ecuación de continuidad en
coordenadas Eulerianas

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuación de Continuidad

𝜕𝜌
+ 𝛻 ∙ (𝜌𝑉) = 0
𝜕𝑡
Ecuación de continuidad en
coordenadas Eulerianas

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuación de Continuidad
Nota: En el centro del elemento cubico asumo que el punto P
tienen la densidad ρ. Para determinar cuanto vale la densidad
en el centro de la cara izquierda y derecha usamos una
expansión en serie de Taylor.

𝜕𝜌
+ 𝛻 ∙ (𝜌𝑉) = 0
𝜕𝑡
Ecuación de continuidad en
coordenadas Eulerianas

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuación de Continuidad
𝜕𝜌 Se observa que si 𝛻 ∙ 𝑉 = 0 el flujo es incompresible ya
+ 𝛻 ∙ (𝜌𝑉) = 0
𝜕𝑡 que esto implica que el caudal másico entrante es igual al
Ecuación de continuidad en saliente, y la densidad no varia con respecto del tiempo
coordenadas Eulerianas
dentro del volumen de control.

Elemento material
moviéndose a través
de una campo fluido.
𝐷() 𝜕() Su densidad cambia de
𝑆𝑒𝑎 𝑒𝑙 𝑜𝑝𝑒𝑟𝑎𝑑𝑜𝑟 = + 𝑉 ∙ 𝛻()
𝐷𝑡 𝜕𝑡 acuerdo a la ecuación
𝐷𝜌
de continuidad.
+𝜌𝛻∙𝑉 =0
𝐷𝑡
Ecuación de continuidad en
coordenadas Lagrangeanas

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Principio de Conservación de la Energía
Fuerza que actúa sobre una partícula de Considerando el flujo estacionario, y las
fluido a lo largo de una línea de corriente fuerzas actuantes sobre la partícula.
𝐷𝑉
P dA – (P + dP) dA – W sin θ = 𝑚 𝐷𝑡
𝐹 = 𝑚𝑎𝑠
𝜕𝑉
P dA – (P + dP) dA – W sin θ = 𝑚 𝜕𝑡 + 𝑚 𝑉 ∙ 𝛻 𝑉
𝑑𝑉
P dA – (P + dP) dA – W sin θ = 𝑚 𝑉 ∙ 𝑑𝑆
𝑑𝑉
P dA – (P + dP) dA – W sin θ = ρ dA ds 𝑉 ∙ 𝑑𝑆
𝑊 sin 𝜃
-dP – dA = ρ 𝑉 ∙ 𝑑𝑣
𝑑𝑧
Siendo sin θ = 𝑑𝑠 y w = ρ g dA ds
𝜌 𝑔 𝑑𝐴 𝑑𝑠 𝑑𝑧
-dP – 𝑑𝐴 𝑑𝑠 = ρ 𝑉 𝑑𝑣
-dP – ρ g dz= ρ 𝑉 𝑑𝑣
𝑑𝑃 𝑑 (𝑉 2 )
- 𝜌 – g dz= 2

𝑑𝑃 𝑑 (𝑉 2 ) Ecuación diferencial
+ g dz + =0
𝜌 2 de Euler

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Principio de Conservación de la Energía
𝑑𝑃 𝑑 (𝑉 2 )
+ g dz + =0 Ecuación diferencial de Euler
𝜌 2
Integrando los diferenciales exactos
𝑑𝑃 𝑉2
+gz+ = Constante a lo largo de una línea de corriente”
𝜌 2
𝑑𝑃
Resolviendo la integral para flujo estacionario incompresible
𝜌

𝑃 𝑉2 Ecuación
+gz+ = Constante a lo largo de una línea de corriente”
𝜌 2 de Bernoulli

𝑃1 𝑃2
+ 𝑔 𝑍1 + 𝑉12 = + 𝑔 𝑍2 + 𝑉22 = 𝐶𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒
𝜌 𝜌

”La suma de la energía de flujo, cinética y potencial de una


partícula de fluido a lo
largo de una línea de corriente es CONSTANTE, cuando el
flujo es estacionario,
incompresible y los efectos de fricción son despreciables”

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Equilibrio de fuerzas en traslación
Las fuerzas que actúan sobre un volumen de
control se dividen de la siguiente forma:
• Fuerzas Másicas (Fbody)
• Fuerzas Superficiales
• Normales
• Viscosas

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Equilibrio de fuerzas en traslación
Las Fuerzas Superficiales Normales y
Viscosas surgen a partir de integrar en el
área al tensor de tensiones. Este tensor
tiene como componentes en su diagonal
a las tensiones normales (Presiones) +
las tensiones viscosas. Fuera de su
diagonal se encuentran las tensiones de
corte.
Fluido en Movimiento

Fluido en Reposo Tensor Viscoso


(Desviador)

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Equilibrio de fuerzas en traslación

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Teorema de conservación de Cantidad de Movimiento.

(𝑒𝑥𝑡)
𝐷𝑣
𝐹 =𝑚
𝐷𝑡
Tomando un volumen diferencial

(𝑒𝑥𝑡)
𝐷𝑣
𝑑𝐹 = 𝜌 𝑑𝑥 𝑑𝑦 𝑑𝑧
𝐷𝑡

(𝑒𝑥𝑡)
𝐷𝑣
𝑑𝐹 =𝜌 dx dy dz
𝐷𝑡
Tomando la componente x
𝐷𝑢
𝑑𝐹𝑥 = 𝜌 dx dy dz
𝐷𝑡

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Teorema de conservación de Cantidad de Movimiento.

𝜕𝑢 𝜕𝑢 𝜕𝑢 𝜕𝑢
dF𝑥 = 𝜌 +𝑢 +𝑣 +𝑤 dxdydz
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧

Descomponiendo cada término


𝜕𝑢 𝜕 ρu 𝜕𝜌
𝜌 = −𝑢
𝜕𝑡 𝜕𝑡 𝜕𝑡
𝜕𝑢 𝜕 ρu2 𝜕 ρu
ρu = −𝑢
𝜕𝑥 𝜕𝑥 𝜕𝑥
𝜕𝑢 𝜕 ρuv 𝜕 ρv
ρv = −𝑢
𝜕𝑦 𝜕𝑦 𝜕𝑦

𝜕𝑢 𝜕 ρuw 𝜕 ρw
ρw = −𝑢
𝜕𝑧 𝜕𝑧 𝜕𝑧

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Teorema de conservación de Cantidad de Movimiento.

Teorema de conservación de la cantidad de


movimiento.
𝜕 ρu 𝜕𝜌 𝜕 ρu2 𝜕 ρu 𝜕 ρuv 𝜕 ρv 𝜕 ρuw 𝜕 ρw
dF𝑥 = −𝑢 + −𝑢 + −𝑢 + −𝑢 dxdydz
𝜕𝑡 𝜕𝑡 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧
𝜕 ρu 𝜕 ρu2 𝜕 ρuv 𝜕 ρuw 𝜕𝜌 𝜕 ρu 𝜕 ρv 𝜕 ρw
dF𝑥 = + + + −𝑢 + + + dxdydz
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧

𝜕 ρu 𝜕 ρu2 𝜕 ρuv 𝜕 ρuw


dF𝑥 = + + + dxdydz
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧

𝜕 ρu _
dF𝑥 = + 𝛻. ρu𝑉 dxdydz
𝜕𝑡
𝜕 ρu _ 𝜕 ρu _
𝐹𝑥 = + 𝛻. ρu𝑉 dvol = dvol + 𝛻. ρu𝑉 dvol
𝜕𝑡 𝜕𝑡 vol
vol

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Teorema de conservación de Cantidad de Movimiento.
Aplicando el Teorema de Gauss (Divergencia)
𝜕 ρu _ _
𝐹𝑥 = dvol + u ρ(𝑉 . 𝑛) dsup
𝜕𝑡 sup
vol
Generalizando para x, y, z:
_
𝜕 𝜌𝑉 _ _ _
𝐹= dvol + 𝑉 𝜌𝑉 . 𝑛 dsup
𝜕𝑡 sup
vol

_
𝜕 𝜌𝑉 _
𝐹= dvol + 𝑉 d𝐶𝜌
𝜕𝑡 sup
vol

La suma de todas La razón de cambio El flujo neto del momento


Las fuerzas externas Respecto al tiempo lineal hacia afuera de la
= del momento lineal +
que actúan sobre superficie de control por
un VC del contenido del VC el flujo de masa

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuación de Cauchy
Ecuación de Cauchy a partir de un volumen diferencial
En la cara izquierda del elemento cubico asumo que la
momento lineal es m.v., y su flujo m.v.v Para determinar
cuanto vale el momento lineal en el centro de la cara de la
derecha usamos una expansión en serie de Taylor. En las
graficas solo se muestra las componente del momento
lineal en X y de la Fsup X a través del V.C.

Ecuación de Cauchy

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuación de Cauchy
Ecuación de Cauchy a partir de un volumen diferencial
En el centro del elemento cubico asumo que la momento
lineal es m.v., y su flujo m.v.v Para determinar cuanto vale
el momento lineal en el centro de la cara izquierda y
derecha usamos una expansión en serie de Taylor. En las
graficas solo se muestra las componente del momento
lineal en X y de la Fsup X a través del V.C.

Ecuación de Cauchy

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuación Indefinidas de la Dinámica
Aplicando la regla de Producto:

Ecuación de Cauchy

Forma alternativa de la ecuación de Cauchy

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuación Indefinidas de la Dinámica
Ecuación de Cauchy

𝜕𝑃 𝜕𝜏𝑥𝑥 𝜕𝜏𝑦𝑥 𝜕𝜏𝑧𝑥


− + + +
𝜏𝑥𝑥 𝜏𝑥𝑦 𝜏𝑥𝑧 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧
−𝑃 0 0 𝜕𝑃 𝜕𝜏𝑥𝑦 𝜕𝜏𝑦𝑦 𝜕𝜏𝑧𝑦
𝜕 𝜕 𝜕 𝜏𝑦𝑥 𝜏𝑦𝑦 𝜏𝑦𝑧
𝛻 ∙ σ𝑖𝑗 = , , ∙ 0 −𝑃 0 + = − 𝜕𝑦 + 𝜕𝑥 + 𝜕𝑦 + 𝜕𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧
0 0 −𝑃 𝜏𝑧𝑥 𝜏𝑧𝑦 𝜏𝑧𝑧 𝜕𝑃 𝜕𝜏𝑥𝑧 𝜕𝜏𝑦𝑧 𝜕𝜏𝑧𝑧
− + + +
𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝑃 𝜕𝜏𝑥𝑥 𝜕𝜏𝑦𝑥 𝜕𝜏𝑧𝑥 𝐷𝑢
𝜌𝑔𝑥 − + + + = 𝜌
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝐷𝑡
𝜕𝑃 𝜕𝜏𝑥𝑦 𝜕𝜏𝑦𝑦 𝜕𝜏𝑧𝑦 𝐷𝑣
𝜌𝑔𝑦 − + + + = 𝜌
𝜕𝑦 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝐷𝑡
𝜕𝑃 𝜕𝜏𝑥𝑧 𝜕𝜏𝑦𝑧 𝜕𝜏𝑧𝑧 𝐷𝑤
𝜌𝑔𝑧 − + + + = 𝜌
𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝐷𝑡

Ecuaciones Indefinidas de la
Se introduce un Nuevo tensor Dinámica de Fluidos
con las tensiones de corte 𝐷𝑉
ρ𝑔 − 𝛻𝑝 + 𝛻 ∙ τ = ρ
como valores desconocidos 𝐷𝑡

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuación de Euler
Considerando zona potencial, se tiene:

Ecuación de Euler
𝐷𝑉
ρ𝑔 − 𝛻𝑝 = ρ
𝐷𝑡

Ecuaciones
de Euler

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Equilibrio de Momentos en Rotación

Para Flujo Estacionario: Si r y el flujo de momento están en el mismo plano:

Ecuación de Euler para


Turbo máquinas

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Equilibrio de Momentos en Rotación

Despreciando los términos de orden superior,

Análogamente, en los otros planos,

NO EXISTE ECUACION DIFERENCIAL DEL MOMENTO CINETICO. La aplicación de la


ecuación en forma integral a un volumen de control infinitesimal proporciona el
resultado bien conocido de que las tensiones de corte son simétricas.

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuaciones de Lamé
Una componente i-j del tensor de tensiones, para un fluido en movimiento (caso contrario
seria nula), se puede escribir como:

Donde, sij es una componente del tensor de tensiones,


p es la presión, ij es una componente de la matriz
identidad y tij es una componente del tensor desviador
de tensiones (contiene las tensiones de corte).

El tensor desviador de tensiones puede escribirse como:

Donde, Kijpq es una matriz que describe la relación


constitutiva del material en un fluido newtoniano. Kijpq =
Es la matriz de elasticidad, que originalmente tiene 81
componentes, y que por simetría de los tensores puede
reducirse a 36 y considerándose isotrópico
(elásticamente equivalente en todas las direcciones),
puede escribirse hacienda uso de solo dos constantes
de elasticidad, las constantes de Lame.

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuaciones de Lamé
Las constantes l y m son coeficientes de viscosidad’, análogos a los parámetros
de Lamé. Reformulando la expresión del tensor de tensiones se tiene que:

𝜎𝑖𝑗 = −𝑝𝛿𝑖𝑗 − λ∗ 𝛿𝑖𝑗 𝐷𝑘𝑘 + 2𝜇∗ 𝐷𝑖𝑗

Donde, sij es una componente del tensor de tensiones, p es la presión, ij es una
componente de la matriz identidad, Dkk es la traza de la matriz Dpq y Dij es una
componente de la matriz Dpq.

La nulidad de la viscosidad volumétrica se conoce como condición de stokes.

1 2 2 ∗
𝑘 ∗ = 3 λ∗ + 3 𝜇 ∗ = 0 λ∗ = 𝜇
3

Luego, el tensor de tensiones, para un fluido newtoniano, aceptando la


condición de Stokes resulta:
2 ∗
𝜎𝑖𝑗 = −𝑝𝛿𝑖𝑗 − 𝜇 𝛿𝑖𝑗 𝐷𝑘𝑘 + 2𝜇∗ 𝐷𝑖𝑗
3

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuaciones de Lamé
Ecuaciones de Navier Stokes
Insertando la expresión anterior en las ecuaciones indefinidas se tiene que:

𝑑𝑉
𝜌 𝑔− = 𝛻𝑃 − 13𝜇𝛻(𝛻 ∙ 𝑉) − 𝜇𝛻 ∙ 𝛻𝑉
𝑑𝑡

Esta expresión se conoce como Ecuación de Navier – Stokes.

En el caso de fluido incompresible (𝛻 ∙ V = 0) la expresión anterior puede


escribirse como:

𝑑𝑉
𝜌 𝑔− = 𝛻𝑃 − 𝜇𝛻 ∙ 𝛻𝑉
𝑑𝑡

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuación de Navier Stokes
Fluido en Movimiento

Tensor Viscoso
Fluido en Reposo
(Desviador)

Se introduce un Nuevo tensor


con las tensiones de corte
como valores desconocidos

Se define el tensor de deformación viscoso. Se considera Flujo


Incompresible y adiabático (esto implica viscosidad constante)

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuación de Navier Stokes
Partiendo de la Ecuación de Cauchy:
El producto puede escribirse de la siguiente forma:

Consideramos el flujo incompresible

Analizamos solo la dirección en X

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuación de Navier Stokes

Ecuaciones de Navier Stokes

Ecuación de Navier Stokes

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Ecuaciones de Euler y de la estática
Partiendo de la Ecuación de Navier – Stokes, si se considera un fluido ideal:

𝐷𝑉
𝜌 𝑔− = 𝛻𝑃
𝐷𝑡

Esta expresión se conoce como ecuación de Navier Sotkes para fluidos ideales.

En el caso de la estática de los fluidos (V = 0), puede reducirse a:

𝜌𝑔 = 𝛻𝑃

Esta expresión da lugar a la ecuación diferencial fundamental de la estática,


que en un espacio unidimensional puede reducirse a:

d𝑃 = 𝜌 𝑔 dℎ
Ecuación General de la Estatica

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Principio de Conservación de la Energía
Partiendo de la ecuación de Navier Stokes para fluidos ideales :
𝐷𝑉
𝜌 𝑔− = 𝛻𝑃
𝐷𝑡
Multiplicando escalarmente ambos miembros por el vector ds = dx i + dy j + dz k
y pasando la densidad al término de la derecha:

𝐷𝑉 𝛻𝑃
𝑔 ∙ 𝑑𝑠 − ∙ 𝑑𝑠 = ∙ 𝑑𝑠
𝐷𝑡 𝜌

I II III

Se analiza cada miembro por separado.

I 𝑔 ∙ 𝑑𝑠 = 𝕏, 𝕐, ℤ ∙ 𝑑𝑥, 𝑑𝑦, 𝑑𝑧 = 𝕏 𝑑𝑥 + 𝕐 𝑑𝑦 + ℤ 𝑑𝑧

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Principio de Conservación de la Energía
Dado que el campo vectorial g deviene del potencial gravitatorio es posible
expresarlo como el gradiente de una función potencial 𝑈𝑚 :
𝜕𝑈𝑚 𝜕𝑈𝑚 𝜕𝑈𝑚
𝑔 = 𝕏, 𝕐, ℤ = , , = 𝛻𝑈𝑚
𝜕𝑥 𝜕𝑦 𝜕𝑧
Aplicando el producto escalar por un diferencial de línea ds :
𝜕𝑈𝑚 𝜕𝑈𝑚 𝜕𝑈𝑚
𝑔 ∙ 𝑑𝑠 = 𝛻𝑈𝑚 ∙ 𝑑𝑠 = 𝑑𝑥 + 𝑑𝑦 + 𝑑𝑧 = 𝑑𝑈𝑚
𝜕𝑥 𝜕𝑦 𝜕𝑧
II 𝐷𝑉 𝐷𝑢 𝐷𝑣 𝐷𝑤
∙ 𝑑𝑠 = , , ∙ 𝑑𝑥, 𝑑𝑦, 𝑑𝑧
𝐷𝑡 𝐷𝑡 𝐷𝑡 𝐷𝑡
𝐷𝑉 𝐷𝑢 𝐷𝑣 𝐷𝑤 𝑑𝑥 𝑑𝑦 𝑑𝑧
∙ 𝑑𝑠 = 𝑑𝑥 + 𝑑𝑦 + 𝑑𝑧 = 𝑑𝑢 + 𝑑𝑣 + 𝑑𝑤
𝐷𝑡 𝐷𝑡 𝐷𝑡 𝐷𝑡 𝑑𝑡 𝑑𝑡 𝑑𝑡
𝐷𝑉 𝑑𝑥 𝑑𝑦 𝑑𝑧
∙ 𝑑𝑠 = 𝑑𝑢 + 𝑑𝑣 + 𝑑𝑤 = 𝑢 𝑑𝑢 + 𝑣 𝑑𝑣 + 𝑤 𝑑𝑤
𝐷𝑡 𝑑𝑡 𝑑𝑡 𝑑𝑡
𝐷𝑉 𝑑 𝑢2 𝑑 𝑣2 𝑑 𝑤2 𝑑(𝑉 2 )
∙ 𝑑𝑠 = + + =
𝐷𝑡 2 2 2 2
Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Principio de Conservación de la Energía
III
𝛻𝑃 1 𝜕𝑃 𝜕𝑃 𝜕𝑃 1 𝜕𝑃 𝜕𝑃 𝜕𝑃
∙ 𝑑𝑠 = , , ∙ 𝑑𝑥, 𝑑𝑦, 𝑑𝑧 = ( 𝑑𝑥 + 𝑑𝑦 + 𝑑𝑧)
𝜌 𝜌 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜌 𝜕𝑥 𝜕𝑦 𝜕𝑧

𝛻𝑃 𝑑𝑃
∙ 𝑑𝑠 =
𝜌 𝜌

Finalmente, sumando los tres términos se tiene la ecuación diferencial de Euler:

𝐷𝑉 𝛻𝑃 𝑑𝑃 𝑑(𝑉 2 )
𝑔 ∙ 𝑑𝑠 − ∙ 𝑑𝑠 = ∙ 𝑑𝑠 + − 𝑑𝑈𝑚 = 0
𝐷𝑡 𝜌 𝜌 2
Ecuación diferencial de Euler

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Principio de Conservación de la Energía
Integración y Sistema de Bernoulli

Integrando la expresión anterior se tiene:

𝑑𝑃 𝑑(𝑉 2 )
+ − 𝑑𝑈𝑚 = ℂ
𝜌 2
𝑑𝑃 𝑉 2
+ − 𝑈𝑚 = ℂ
𝜌 2
Integral Generalizada de Bernoulli

La integral que involucra dP y  requiere un análisis adicional.

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Principio de Conservación de la Energía
Integración y Sistema de Bernoulli

Casos particulares del sistema de Bernoulli


Los casos de interes para la integración del sistema de Bernoulli son los
siguientes:
1. Isopínica (Incompresible,  es constante).
2. Politrópica. Simplificación adiabática, expresión the Saint Venant -
Wenczel.
𝑑𝑃 𝑃
= ℂ
𝜌 𝜌𝑛

1. Incompresible:
𝑑𝑃 1 𝑃
𝜌= ℂ = 𝑑𝑃 =
𝜌 𝜌 𝜌

𝑃 𝑉2
+ − 𝑈𝑚 = ℂ Ecuación de Bernoulli para fluidos incompresibles
𝜌 2

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Principio de Conservación de la Energía
Integración y Sistema de Bernoulli
2.Poli trópica: dada la complejidad de la poli trópica, se toma un caso particular
de la misma que es la adiabática:
𝑃
= ℂ
𝜌𝛾
1 𝛾−1
−𝛾+1
𝑃 𝛾1 𝑃 −𝛾1 𝑃 𝑃 𝛾
𝜌= ( ) → ( ) 𝑑𝑝 → ℂ1/𝛾 𝑃−1/𝛾 𝑑𝑃 → ℂ1/𝛾 → ℂ1/𝛾
ℂ ℂ 1 𝛾−1
−𝛾 + 1 𝛾
1
𝛾 𝛾−1
𝑃 𝛾 𝛾 𝑃
→ 𝑃 𝛾 →
𝛾−1 𝜌 𝛾−1𝜌

𝛾 𝑃 𝑉2
+ − 𝑈𝑚 = ℂ Ecuación de Bernoulli para flujo adiabático.
𝛾−1𝜌 2

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Principio de Conservación de la Energía
Integración y Sistema de Bernoulli

Casos particulares del sistema de Bernoulli .


Ecuación de Saint Venant - Wenczel
Partiendo de la Ecuación de Bernoulli para flujo adiabático y recordando que
para gases ideales 𝑃 = 𝜌𝑅𝑇, 𝑅 = 𝐶𝑝 − 𝐶𝑣 y 𝛾 = 𝐶𝑝 /𝐶𝑣 se tiene:
𝐶𝑝 /𝐶𝑣 𝜌𝑅𝑇 𝑉 2 𝐶𝑝 𝑉2
+ − 𝑈𝑚 = 𝑅𝑇 + − 𝑈𝑚 = ℂ
𝐶𝑝 /𝐶𝑣 − 1 𝜌 2 𝐶𝑝 − 𝐶𝑣 2

𝐶𝑝 𝑉2
𝑅𝑇 + − 𝑈𝑚 = ℂ
𝑅 2

Si se considera despreciable el termino geodésico (Um)

𝑉2
𝐶𝑝 𝑇 + − 𝑈𝑚 = ℂ Ecuación de Saint Venant - Wenczel.
2

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Principio de Conservación de la Energía
Integración y Sistema de Bernoulli
Diagrama Piezómetro
Dividiendo la ecuación de Bernoulli para fluidos incompresibles por la gravedad
se tiene que:
1 𝑃 𝑉2 𝑃 𝑉 2 𝑈𝑚
+ − 𝑈𝑚 = + − =ℂ
𝑔 𝜌 2 𝜌𝑔 2𝑔 𝑔
Donde: 𝛾 = 𝜌𝑔 es el peso específico y 𝑈𝑚 = −g z , con lo se tiene:
𝑃 𝑉2
+ +𝑧 =ℂ
𝜌𝑔 2𝑔
𝑃 𝑉2
EGL: Línea de energía (𝜌𝑔 + 2𝑔 + 𝑧)
𝑃
HGL: Línea Hidráulica (𝜌𝑔 + 𝑧)
P
= Carga de Presion (columna de presion Estatica)
ρg
V2
2g
= Carga de Velocidad (columna necesaria para lograr V)
z = Carga de altura (Energía Potencial del fluido)

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Corrientes relativas. Bomba centrífuga.

Objetivo: aumentar la presión total de un líquido.


La bomba es una turbo máquina que entrega energía al fluido.

∆𝑃= PB - PA
Anti-retorno
Energía B > Energía A
Tubería de
descarga H1 > H0
Tanque B

H1 B
H0
Tubería de
aspiración
Tanque A
Filtro

• Funcionamiento V P
• Cavitación

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Corrientes relativas. Bomba centrífuga.

H (carga de la Bomba) = hout - hin, energía aportada,


despreciando efectos viscosos y térmicos.

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Corrientes relativas. Bomba centrífuga.

El BEP es el punto de mayor eficiencia.

El punto de operación se corresponde cuando el H


disponible es igual al H requerido por el sistema.

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Corrientes relativas. Bomba centrífuga.
Bomba centrífuga. Cavitación

Para evitar la cavitación es necesario evitar que el fluido este


por debajo de su presión de vapor en cualquier punto dentro
de la bomba.
El criterio de cavitación suele especificarse a la entrada de la
bomba.

“ Si la carga neta de aspiración positiva NPSH requerida es


mayor que la NPSH disponible, existe cavitación.”

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Corrientes relativas. Bomba centrífuga.

Curvatura del álabe:


• b < 90°, curvatura hacia atrás. El ensanchamiento es gradual. Se recomienda
valor de curvatura de 25° < b < 60°.
• b = 90°, sin curvatura o álabe radial.
• b > 90°, curvatura hacia adelante. Mayor rozamiento por la distancia recorrida
por el fluido. Puede haber Cavitación.

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Corrientes relativas. Bomba centrífuga.

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Corrientes relativas. Bomba centrífuga.
Teóría de hilo de corriente.
Hipótesis:
• Los álabes son infinitos.
• Los álabes tienen espesor diferencial.

El torque en el eje es equivalente al cambio de cantidad de


momento entre la salida 2 y la entrada 1. Aplicando la
ecuación de Euler para turbo maquinas.

Considerando un rendimiento de la bomba = 1

Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Corrientes relativas. Bomba centrífuga.
Aplicando Teorema de Coseno
(1)
Según diagrama vectorial
(2)

Sustituyendo (2) en (1)

(3)

Aplicando Euler y considerando un rendimiento = 1

(4)

Sustituyendo (3) (para sub índice 2 y 1) en (4)


(5)
La carga neta es igual al cambio de energía cinética absoluta mas el cambio de la energía
cinética de la puntera del rotor menos la energía cinética relativa entre la entrada y salida.
Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica
Corrientes relativas. Bomba centrífuga.
(5)

(6)

Igualando (5) con (6)

Ecuación de Bernoulli para en un marco rotacional


2 2
𝜔2 𝑟22 − 𝑟12 𝑉1𝑟𝑒𝑙𝑎𝑡𝑖𝑣𝑒 − 𝑉2𝑟𝑒𝑙𝑎𝑡𝑖𝑣𝑒
∆𝑃 = 𝑃2 − 𝑃1 = 𝜌 +
2 2
Nótese que incrementar 𝜔 y 𝑉𝑟𝑒𝑙𝑎𝑡𝑖𝑣𝑒 aumentará el salto de presión de la bomba. Pero
existen límites:
• Cavitación – Erosión.
• La velocidad relativa (𝑉𝑟𝑒𝑙𝑎𝑡𝑖𝑣𝑒 ) limita el ángulo 𝛽 del álabe.
• Agrandar el rotor (>𝑟2 ) incrementa las pérdidas por rozamiento.
• Elevados valores de 𝜔 comprometen estructuralmente la bomba.
Mecánica de Los Fluidos Ingenieria Industrial Unidad 4: Dinámica y teoremas de conservación de la dinámica

También podría gustarte