Mecanica Fluidos
Mecanica Fluidos
Mecanica Fluidos
CAMPUS TECNOLGICO DE LA UNIVERSIDAD DE NAVARRA. NAFARROAKO UNIBERTSITATEKO CAMPUS TEKNOLOGIKOA Paseo de Manuel Lardizbal 13. 20018 Donostia-San Sebastin. Tel.: 943 219 877 Fax: 943 311 442 www.tecnun.es [email protected]
ISBN
Reservado todos los derechos. Queda prohibida la reproduccin total o parcial sin autorizacin previa. Primera Edicin: 2008
MECNICA DE FLUIDOS
INDICE (I)
1 CONCEPTOS INTRODUCTORIOS
1.1 Prlogo 1.2 La Hiptesis del Continuo 1.2.1 Partcula Material 1.2.2 Sistema Material 1.3 Definicin de Fluido 1.4 Enfoques Lagrangiano y Euleriano 1.4.1 Enfoque Lagrangiano 1.4.2 Enfoque Euleriano 1.4.3 Ejemplo Final 1.5 Leyes Fundamentales 1.5.1 Consideraciones Finales
3
4 6 9 10 11 15 16 21 24 27 28
30
31 39 41 45 55 58 64 66 67 70 74 75 76 78 81 89 97
MECNICA DE FLUIDOS
INDICE (II)
5 INSTALACIONES HIDRULICAS
5.1 Generalidades 5.1.1 Definicin y Modelado de una instalacin hidrulica 5.1.2 Elementos de una instalacin hidrulica 5.2 Prdidas de carga en tuberas 5.2.1 Ecuacin de Darcy-Weisbach 5.2.2 Secciones no circulares. Dimetro Hidrulico 5.2.3 Problemas Bsicos en tuberas 5.3 Vlvulas 5.3.1 Funciones y Tipos 5.3.2 Prdidas de carga en vlvulas 5.4 Modelo Matemtico de una instalacin hidrulica 5.4.1 Ecuaciones Fundamentales 5.4.2 Condiciones de Contorno 5.4.3 Resolucin 5.4.4 Formulacin por Caudales 5.4.5 Formulacin por Alturas
147
148 149 150 154 155 160 163 166 167 169 171 172 173 174 175 177
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
1.1 PRLOGO
MECNICA DE FLUIDOS
1.1 Prlogo
MECNICA DE FLUIDOS: Estudia las fuerzas y energas que los fluidos generan cuando se encuentran en reposo (fluidoesttica) y en movimiento (fluidodinmica). Naturaleza de los fluidos (Qu son?). Magnitudes que se emplean para analizar su comportamiento. Leyes Fsicas que gobiernan su comportamiento. Modelacin matemtica y resolucin de problemas de Mecnica de Fluidos. PREGUNTA: Por qu es necesario para un ingeniero estudiar Mecnica de Fluidos?. RESPUESTA: Porque son amplsimos los campos de la ingeniera donde aparecen fluidos en movimiento o reposo. Transporte de Fluidos. Generacin de Energa. Control Ambiental Transporte
Pocos son los ingenieros que pueden desempear su funcin de manera efectiva sin, por lo menos, un conocimiento rudimentario de la Mecnica de Fluidos. . P. M. Gerhart, R. J. Groos and J.
I. Hochstein. Fundamentos de Mecnica de Fluidos. Ed. Addison-Wesley
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
SOLUCIN: Adoptar un modelo (aproximacin) de la materia que se denomina CONTINUO. La materia es continua, llena todo el espacio que ocupa. (macroscpicamente lo parece). unas propiedades macroscpicas que definen el comportamiento de la materia y que varan continuamente en el espacio ocupado esta.
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
(VP)
una
superficie
(SP)
Su comportamiento se analiza a travs de la evolucin de una serie de magnitudes: Cinemticas: Movimiento y deformacin. Ejemplos r P, v P, a P y D P. Dinmicas: Fuerzas que actan sobre la partcula. Ejemplos g, pP y Tp. Termodinmicas: Transferencia de energa entre la partcula y su entorno. Ejemplos: pP, TP y P
MECNICA DE FLUIDOS
Cualquier magnitud asociada al sistema material puede expresarse a partir de las propiedades de las partculas materiales que lo constituyen (siempre las mismas). Ejemplo: Cantidad de Movimiento M (t ) = P (t ) vP (t ) dVP
V (t )
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
CONCLUSIN: Un slido adquiere una deformacin esttica ante una solicitacin cortante.
MECNICA DE FLUIDOS
CONCLUSIN: Un fluido NO adquiere una deformacin esttica ante una solicitacin cortante, por el contrario se deforma continuamente (velocidad de deformacin).
& F
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
Ejemplo: Posicin r(X0,t). Fijando t : r(X0,t) proporciona la posicin de todas las partculas en ese instante. Fijando X0 : r(X0,t) proporciona la evolucin temporal de la posicin de la partcula que en t=0 tiene las coordenadas X0.
Y X
rP(X0P,t)
vP(X0P,t=0)
vP(X0P,t)=dr/dt
Partcula Material P en t
Sistema Material en t
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
t (+)
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
S.C. V.C.
S.C. e V.C. s
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
BECARIO (1)
4
INGENIERO Jr. (2)
6
INGENIERO Sr. (3)
8
INGENIERO Sr. (3)
10
GERENTE (4)
14 Ao
16
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
Conservacin de la masa: La rapidez de variacin en el tiempo de la masa de un sistema material es nula. La masa de un sistema material es constante. 2 Ley de Newton: La rapidez de variacin en el tiempo de la cantidad de movimiento de un sistema material es igual a la fuerzas externas que actan sobre l. 1 Ley de la Termodinmica: La rapidez de variacin en el tiempo de la energa de un sistema material es igual a la velocidad de transferencia neta de energa entre el sistema y su entorno. 2 Ley de la Termodinmica: (Poco empleada en este curso de Mecnica de Fluidos. No enunciada)
MECNICA DE FLUIDOS
v P (t ) = v (x P , t )
t 0
v P (t + t ) v (x P , t + t )
a (x, t ) Dv Dt v t lim
t 0
v (x , t + t ) v (x , t )
MECNICA DE FLUIDOS
= lim
M (t + t ) M (t )
t 0
DM (t ) Dt
dMVC dt
= lim
MVC (t + t ) MVC (t )
t 0
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
div (v ) = 0 v v = p + 2 v + fV r
+ CONDICIONES DE CONTORNO
X
V.C.
r(x)
e
Y
MECNICA DE FLUIDOS
CONCLUSIN: Adoptar un mtodo de anlisis donde las incgnitas fueran magnitudes integrales en las que el ingeniero est interesado y matemticamente ms simple.
Los nicos flujos que se analizarn por el mtodo diferencial sern flujos en Rgimen Laminar (no turbulencia) y en Geometra muy simples (unidireccionales-unidimensionales)
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
qe = qs = q pe 1 q + 2 g Ae FW K q 2 = ps + 1 q Q 2 g As 1 1 = ps As pe Ae + q 2 A s Ae
2 2
MECNICA DE FLUIDOS
METODO DIFERENCIAL
METODO INTEGRAL
qe = qs = q pe 1 q + 2 g Ae FW K q 2 = ps + 1 q Q 2 g As 1 1 = ps As pe Ae + q 2 A s Ae
2 2
div (v ) = 0 v v = p + 2 v + fV r
+ CONDICIONES DE CONTORNO
PREGUNTA: Qu mtodo parece ms sencillo en su formulacin y proporciona respuestas prcticas de forma ms inmediata en su resolucin?
MECNICA DE FLUIDOS
ANLISIS INTEGRAL
Experimentacin o Anlisis Diferencial (C.F.D)
qe = qs = q
2 2
pe 1 q + 2 g Ae FW
K q 2 = ps + 1 q Q 2 g As 1 2 1 = ps As pe Ae + q A s Ae
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
Es la propiedad ms importante en el anlisis de los flujos: Es la principal magnitud cinemtica. Todas las dems mag. cinemticas se definen a partir de ella. En flujos incompresible si se conoce v(x,t) el flujo esta resuelto.
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
r 2 v (r ) = u z (r ) k = U 0 1 k R
MECNICA DE FLUIDOS
Z m
A B
2 2 z 3 y z 3 v( y , z ) = u (y , z ) i = 6 3 U m i a 2 a a
MECNICA DE FLUIDOS
d XX
D (X ) = Dt X
MECNICA DE FLUIDOS
N vN. t
dYY
D (Y ) = Dt Y
Y(t+ t)
E(e) S vS. t
MECNICA DE FLUIDOS
d XX (x, t ) =
d rr =
u r u u 1 u ; d = + r ; d zz = z r z r r
d rr =
u u r u u r u 1 u 1 ; d = + r ; d = + + cot r r r r sen r r
MECNICA DE FLUIDOS
A =
(v E v O ) t
X
B =
(uN uS ) t
Y
EO =
A (v E v O ) v = = t X x
SN =
(u uS ) = u B = N t Y y
MECNICA DE FLUIDOS
D R
1 R = (A B ) 2 1 D = (A + B ) 2
Z =
R 1 A B 1 1 v u = = ( OE + SN ) = 2 t t 2 2 x y t D 1 A B 1 1 v u = + + = ( OE SN ) = 2 t t 2 2 x y t
d XY = dYX =
MECNICA DE FLUIDOS
d XY = d YX =
1 v u 1 w u 1 w v ; d XZ = d ZX = + + + ; d YZ = d ZY = 2 x y 2 x z 2 y z
Todas las velocidades de deformacin se agrupan en una sola magnitud (tensor o matriz) denominada matriz de velocidad de deformacin. Para un flujo tridimensional y tridireccional en coordenadas cartesianas:
dXX dXY dXZ
dYY
u 1 u v + x 2 y x v D y SIM
1 u w + 2 z x 1 v w z + y 2 w z
dZY
dZZ
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
; Y = X = ; Z = 2 y z 2 z x 2 x y
La velocidad de rotacin se puede expresar de forma independiente al sistema de coordenadas como:
1 w
1 u
1 v
= 2 = rot [v (x ,t )] v (x ,t )
MECNICA DE FLUIDOS
1 D (V & V = V Dt
En coordenadas cartesianas
D (V ) D (X Y Z ) = = Dt Dt D (X ) D (Y ) D (Z ) Y Z + X Z + X Y Dt Dt Dt D (V ) = (d xx + d YY + d ZZ ) X Y Z Dt 1 D (V ) & V = = d xx + d YY + d ZZ Tr (D ) V Dt
V = X Y Z
MECNICA DE FLUIDOS
& V (x ,t ) = Tr [D (x ,t )] div [v (x ,t )]
En coordenadas cartesianas:
u v w & V (x, t ) = + + x y z
En coordenadas cilndricas y esfricas respectivamente (Tablas Apuntes):
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
m (t ) =
P (t ) dV
VC
(x, t ) dV
mVC (t )
MECNICA DE FLUIDOS
B (t ) =
VC
(x, t ) (x, t ) dV
BVC (t )
Particularizando para las propiedades masa (m), cantidad de movimiento (M) y energa (E) se tiene que:
m (t ) = M (t ) = E (t ) =
VC
(x, t ) dV
VC
(x, t ) v (x, t ) dV
VC
(x, t ) e (x, t ) dV
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
t
. .
n .
V = v t S cos
Por lo que el caudal volumtrico q en S es:
q =
V = v S cos t
El signo q (definido por cos) indica si el volumen de fluido est entrando (-) saliendo (+) del V.C. Ambas posibilidades se contemplan utilizando:
h' h
q =
V = v n S = v S t
MECNICA DE FLUIDOS
q = dq =
{ vv n dS = v d S
n
Si la integral se realiza para toda la S.C. se obtiene el caudal neto que atraviesa S.C.
q SC =
SC
dq = v n dS = v (x ,t ) d S
SC SC
Considerando nicamente las zonas de la S.C. donde el fluido est entrando (e) o saliendo (s) la expresin de qSC puede escribirse como:
qSC =
v d S + v d S = qs qe s e s e
SC s
1 4 42 3
>0
SC e
1 4 42 3
<0
MECNICA DE FLUIDOS
n .
B = v t S cos
Siendo la propiedad extensiva expresada por unidad de masa (m~1, M~v y E~e). El flujo a travs de S se define como la cantidad de propiedad B que atraviesa S en la unidad de tiempo por lo tanto: & B = v S B = t
> 0 v S > 0 ( Entrada ) 0 < 0 v S < 0 ( Salida ) < 0 v S > 0 ( Entrada ) 0 > 0 v S < 0 ( Salida )
MECNICA DE FLUIDOS
wm S .C .
& & B = d B = v d S
& Bs2
r P (X 0 P ,t)
e1
s1
Z
s2
& m SC =
V .C .
X Y
SC
(v d S )
& Bs1
MECNICA DE FLUIDOS
q SC =
& B SC =
SC
dq = (v v SUP ) d S = v RS d S
SC SC SC
SC
[(v v SUP ) d S ] = (v RS d S )
En las expresiones de los flujos en lugar de la velocidad del fluido v(x,t) aparece la velocidad del fluido relativa a la superficie. Si vSUP(x,t)=0 el V.C. Es fijo e indeformable y se obtienen las expresiones anteriores
MECNICA DE FLUIDOS
v =
q A
El caudal neto que atraviesa la S.C. puede escribirse en funcin de las velocidades medias como: q SC = v A s v A e
A partir del flujo msico y del caudal es posible definir en una superficie la Densidad Promedio como: & m = q De esta forma el flujo msico neto a travs de la S.C. Es:
& m SC =
( q )s ( q )e
s e
MECNICA DE FLUIDOS
& M SC =
(v q ) (v q )
s s e
Flujo Incompresible
& M SC = v q s
) (v q )
s e
(e q ) (e q )
s s e
Flujo Incompresible
& E SC = e q s
) (e q )
e e
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
PREGUNTA: Qu expresin proporcionar en funcin de propiedades integrales la rapidez del cambio en el tiempo de la propiedad B del sistema? RESPUESTA: La rapidez del cambio en el tiempo de la propiedad B de un sistema se define como:
DB (t ) B (t + t ) B (t ) = lim t 0 Dt t
MECNICA DE FLUIDOS
Sistema en t V (t)=VVC(t)
s2
s2
e1
S.C.
e1
S.C.
B (t ) = BVC (t )
MECNICA DE FLUIDOS
DB (t ) dB VC (t ) & = + BSC (t ) Dt dt
Particularizando para: Masa: Cantidad de movimiento:
d M VC dt
x
& + M SC
[M
i + My j + Mz k
d dt
[M
i + My j + Mz k
VC
SC
Energa:
DE dE VC & = + E SC Dt dt
MECNICA DE FLUIDOS
aP =
v (t + t ) v P (t ) Dv P = lim P t 0 Dt t
En general la rapidez del cambio en el tiempo de una propiedad de una partcula de fluido se expresa como:
(t + t ) P (t ) D P = lim P t 0 Dt t
A la rapidez del cambio en el tiempo de una propiedad de una partcula se le denomina Derivada material de la propiedad .
MECNICA DE FLUIDOS
Derivada Material y T.T. de Reynolds/2.4.1 Derivada Material (II) 2.4 T.T. de Reynolds y Derivada Material-2.4.2 Derivada Material (II)
v P (t ) = v (x ,t ) v P (t + t ) = v (x ,t + t ) + v v P (t + t ) = v (x ,t + t ) +
r = v (x ,t ) t v r r ( x ,t + t )
Sustituyendo en el lmite:
a (x ,t ) =
Dv v (x ,t + t ) v (x ,t ) v (x ,t ) = lim + r v (x ,t ) Dt t 0 t
a (x ,t ) =
Dv = Dt
ACELERACI N LOCAL
v (x ,t ) v (x ,t ) + v (x ,t ) r 1 t4 42 3 14 24 4 4 3
ACELERACI N CONVECTIVA
MECNICA DE FLUIDOS
D = Dt
DERIVADA LOCAL
t {
DERIVADA CONVECTIVA
v 2 r3 1
r
Magnitud vectorial (velocidad) /r es una matriz. Si el flujo es estacionario y los campos de propiedades no dependen del tiempo, (x), la parte local de la derivada material es nula y slo existe parte convectiva:
D (x ) v (x ) = Dt r
MECNICA DE FLUIDOS
u y v y w y
Coordenadas Cartesianas:
u z v z w z
ur z u z uZ z
1 ur u r r 1 u ur r + r 1 uZ r
1 ur u r r 1 u ur r + r 1 u r
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
DM = Fext Dt
Si est ocupando el VC en t
Dv p Dt
m p = Fext
Si p est ocupando la posicin x en t
a (x , t ) (x , t ) V = Fext a (x , t ) = v ( x , t ) v (x , t ) v (x , t ) + t r
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
( FV )p = (fV )p V p
Siendo fV la fuerza volumtrica por unidad de volumen que acta sobre una partcula del sistema.
FV =
d FV
fV dV
MECNICA DE FLUIDOS
( FV )p = (fV
V )p = fV (x, t ) V
fV ser funcin de (x,t) expresa la fuerza volumtrica por unidad de volumen que acta sobre las partculas que en el instante t estn ocupando el V.C.
FV
VC
fV (x, t ) dV
Resultante de las fuerzas de volumen que actan sobre el sistema que en el instante t est ocupando el Volumen de Control.
MECNICA DE FLUIDOS
Las fuerzas de volumen son normalmente conocidas siendo la ms habitual la debida a un campo gravitatorio:
=g es el peso especfico del fluido. []=FL-3=ML-2T-2 Otras fuerzas que pueden aparecer tambin son: Fuerzas de inercia (Sistema de Referencia No Inercial):
fVI = (a arr + a c )
MECNICA DE FLUIDOS
(E)
ENTORNO F S
(E)
(E)
(E)
E F (E ) = F ( ) t ( ) = t (E )
Tercera Ley de Newton t es el vector tensin que expresa la fuerza por unidad de superficie.
MECNICA DE FLUIDOS
t P ( 1 ) t P ( 2 )
t = f ( punto , n )
MECNICA DE FLUIDOS
t = f ( punto, n)
F( ) = T n( ) S
F(E ) = T n(E ) S
Como n()=-n(E)
t = T(punto ) n
F( ) = F(E )
( )
F( ) = F(E ) =
T n( ) dS
(S )
(E )
T n(E ) dS = T n(S ) dS
Resultantes de las fuerzas de superficie que actan sobre el sistema (la ejerce el entorno) y sobre su entorno (la ejerce el sistema).
MECNICA DE FLUIDOS
F ( ) F (E )
SC
T (x, t ) n dS T (x, t ) n dS
SC
Resultantes de las fuerzas de superficie que el Entorno ejerce sobre el sistema que en el instante t est ocupando el Volumen de Control y la que este sistema ejerce sobre el Entorno.
MECNICA DE FLUIDOS
F ( )
SC
T (x, t ) n dS
T n dS
e
F ( )
SC
T (x, t ) n dS =
SC = s + e + w + w m
s
F ( ) = Fs + s
Fe +Fw + Fw m ( ) e
F (E ) = F ( )
MECNICA DE FLUIDOS
F = t E S E + tO S O + t N S N + t S S S
t E = TE n E ; t O = TO nO ; t N = TN n N ; t S = TS n S
S E = S O = Y Z ; S N = S S = X Z
n E = i ; n O = -i ; n N = j ; n S = j
F = (TE TO ) i Y Z + (TN TS ) j X Z
MECNICA DE FLUIDOS
TE = T (x e , t ) = T (x, t ) +
t t t t t t F = xx i + yx j + xy i + yy j X Y Z = xx + xy y x y y x x
t t i + yx + yy x y
j V
Resultante de las fuerzas de superficie que actan sobre una partcula de fluido que en el instante t se encuentra en x.
F = div [T (x, t )] V
MECNICA DE FLUIDOS
1 (r t rr ) 1 t r t rz t r r + r + z r 2 t 1 t 1 r tr div [T (x ,t )]= 2 + + z r r z r 1 (r t zr ) + 1 t z + t zz r r r z
t r t + t 1 r 2 t rr (t sen ) 1 1 + r + 2 r r sen r r sen r 3 t t cot (t sen ) 1 1 1 r tr + + div [T (x ,t )]= 3 r sen r sen r r r 3 t sen t t cot 1 1 1 r tr + + + r 3 r sen r sen r r
MECNICA DE FLUIDOS
DM = Fext Dt
Si est ocupando el VC en t
Fext
Dv p Dt m p = Fext
= FV + F ( )
VC
fV dV +
SC
T n dS
a V = Fext
Si p est ocupando la posicin x en t
Fext
= FV + F = [fV + div (T )] V
MECNICA DE FLUIDOS
Fext = [fV
+ div (T )] V
Dv V = [fV + div (T )] V Dt
PREGUNTA: En un flujo incompresible las incgnitas eran p(x,t) y v(x,t) Dnde aparece la presin?. Como se eliminan las tensiones (T)
RESPUESTA: Es necesario introducir una RELACIN CONSTITUTIVA RESPUESTA: Es necesario introducir una RELACIN CONSTITUTIVA
MECNICA DE FLUIDOS
Propiedades Dinmicas/2.5.3 Relacin Constitutiva un Fluido Newtoniano (III) 2.5 Magnitudes Dinmicas-2.5.5 Relacin Constitutiva dede un Fluido Newtoniano (II)
En general las relaciones constitutivas son expresiones entre las variables de un problema que permiten modelar el comportamiento del sistema estudiado. Se deducen de la experiencia no de las leyes fundamentales (sin violarlas). Ejemplo: Viga sometida a esfuerzo axial se desea relacionar la deformacin con la solicitacin N: Relacin Constitutiva
A N = 0
= E
=
L
N E A
MECNICA DE FLUIDOS
La ecuacin constitutiva de un fluido relaciona su matriz de tensiones (T) con : Presin (p) Velocidad de deformacin (D): La matriz de tensiones se descompone en una Parte Hidrosttica (TH) y en otra Parte Desviadora (TD). Cumplindose que tr(TD)=0.
T = TH + TD
1 TH = tr (T ) I 3 4 123 4
TD = T TH = T I
MECNICA DE FLUIDOS
F ( ) = (TH ) n( ) S ( ) = I n( ) S ( ) = n( ) S ( )
1. En un fluido en reposo slo existe hidrosttica de la matriz de tensiones. parte
2. La parte desviadora estar asociada al movimiento y por tanto a la velocidad. Evidencias experimentales establecen que: El valor coincide con la presin termodinmica ( =-p). TD se relaciona con la matriz de velocidades de deformacin D. La relacin ms sencilla es:
1 T (x, t ) = p I + 2 D Tr (D) I 3
MECNICA DE FLUIDOS
1 T (x, t ) = p I + 2 D Tr (D) I 3
Al coeficiente se le denomina viscosidad dinmica del fluido o simplemente viscosidad y es una propiedad termodinmica (depende de la presin y la temperatura). Los fluidos que poseen la relacin constitutiva se les denomina FLUIDOS NEWTONIANOS. Los fluidos ms comunes son newtonianos (ejemplos: agua, aire, mercurio, aceites) En ocasiones se trabaja con la denominada viscosidad cinemtica =/. []=ML-1.T-1 en el S.I. sus unidades son kg/(ms) Ns/m2. []=L2T-1 en el S.I. sus unidades son m2/s.
MECNICA DE FLUIDOS
F
Siendo:
= [ p + 2 div (D )] V
( p + 2 v ) V =
2w 2w 2w j+ x 2 + y 2 + z 2 k
2u 2u 2u 2v 2v 2v i + 2 + v = 2 + + + 2 2 2 x x y z 2 y z
2
Fext = p + 2 v + fV V
MECNICA DE FLUIDOS
u 1 u r u r u z + + r r r z r 1 u z u 1 u u r + + p + 2 z r r r u z p + 2 z
u 1 u r + r r r r sen u u 1 1 u u r + p + 2 + sen r sen r r r u u r u 1 p + 2 r sen + r + r cot 1 u u r +r r r r sen
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
As por ejemplo el mercurio (Hg) tiene una densidad relativa al agua swHg=13.6. Los fluidos de referencia suelen ser fluidos comunes cuya densidad es bien conocida, como por ejemplo el agua en el caso de lquidos y el aire en el de gases. En Mecnica de Fluidos, sobre todo en Esttica de Fluidos, suele emplearse mucho, el Peso Especfico de un fluido ()
= g
[]=M.L-2.T-2=F.L-3 y sus unidades en el S.I. son kg/(m2.s2) N/m3
MECNICA DE FLUIDOS
Presin absoluta (pabs). Valor siempre (+). Rango: 0<pabs< Presin Manomtrica (pman=pabs-patm). Valor (+) (-). Rango: -patm<pman< El valor de diseo de la presin atmosfrica local es de 101.3 kPa=1003 mBar
MECNICA DE FLUIDOS
p = ha a = hb b
Si se elige otro fluido (c) para expresar como una altura la presin p, las diferentes alturas que expresan la misma presin se relacionan mediante:
p = ha a = hb b = hc c
ha = hb
[ha]=L de a. En el S.I. metros de columna de a.
b
a {
a sb
= hc
c
a {
a sc
MECNICA DE FLUIDOS
pa = p0 + ha
pa p0 = (ph )a = ha
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
Dm =0 Dt
& ms 2
& m e1
Si el sistema que se considera es aquel que en el instante t est ocupando el volumen de control:
MECNICA DE FLUIDOS
dmVC + ( q )s ( q )e = 0 dt s e
En el caso de flujo incompresible (=cte) entonces mVC=VVC la ecuacin queda como:
dVVC + qs qe = 0 dt s e
dVVC dt
+ (v A )s (v A )e = 0
s e
MECNICA DE FLUIDOS
& & ms me = 0
s e
q s q e = 0
s e
En el caso de un VC con una entrada y una salida: Flujo compresible estacionario =(x)
& me & ms
qs = qe = q
(v A )s = (v A )e = q
MECNICA DE FLUIDOS
D(m )
Dt
=0
m = V
& m e1
& ms 2
D( V ) D D(V ) = V + =0 Dt Dt Dt
Introduciendo la Volumtrica: Velocidad de Deformacin
& ms 1
D( V ) D & = V + V V Dt Dt
Considerando una partcula que en el instante t est ocupando una posicin x en el V.C. r=r(x), v=v(x,t) y (x,t)
MECNICA DE FLUIDOS
div (v ) = 0
Un caso particular muy interesante es el flujo incompresible, completamente desarrollado y en rgimen laminar en conductos (conductos rectos de gran longitud) de cualquier tipo de seccin.
u = u (y,z ) ; v = w 0
Se satisface la ecuacin de continuidad
u v w + + =0 x y z { { {
=0 =0 =0
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
DM = Fext Dt
& ms 2
& m e1
Si el sistema que se considera es aquel que en el instante t est ocupando el volumen de control:
( )
SC
= (Fx )ext
MECNICA DE FLUIDOS
Utilizando la cantidad de movimiento por unidad de masa promedio en las superficies de entrada y salida se puede escribir:
& M SC = ( q v )s ( q v )e
s e
Fext = FV + F( )
+ F ( )
e
d MVC + dt
( q v )s ( q v )e = FV
s
MECNICA DE FLUIDOS
( q v )s ( q v )e = FV
s e
+ F
F ( ) =
Fe ( ) + Fs ( ) + +FW ( ) + FW
e s
( )
En las entradas y en las salidas las fuerzas de superficie se descomponen en suma de una debida a las presiones y otra debida a la viscosidad.
F ( ) =
(Fp + F )e ( ) + (Fp + F )s ( ) + FW ( ) + FW
e s
( )
Despreciando las fuerzas viscosas en las entradas y las salidas y considerando superficies planas:
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
q v s = (v n )s + (v t )s ( v n + v SUP )s = n + v SUP { A s
0
q v e = (v n )e + (v t )e ( v n + v SUP )e = n + v SUP 13 2 A e
0
MECNICA DE FLUIDOS
q [( v n )s + ( v n )e ] = FV + FW ( ) + FW m ( ) (p A n )s (p A n )e
MECNICA DE FLUIDOS
Considerando una partcula que en el instante t est ocupando una posicin x en el V.C. r=r(x) y v=v(x,t). Para un fluido newtoniano y suponiendo flujo incompresible:
Fext
= FV + F = [fV + div (T )] V T = p I + 2 D
F = div ( T ) V = p + 2 v V
MECNICA DE FLUIDOS
2u 2u 2u v = u i + v j + w k = 2 + 2 + 2 x y z
2 2 2 2
2v 2v 2v i + + + x 2 y 2 z 2
2w 2w 2w j+ + + x 2 y 2 z 2
En coordenadas cilndricas o esfricas las expresiones se pueden encontrar en las Tablas de los Apuntes:
MECNICA DE FLUIDOS
v v + v = p + 2 v + fV 2 t r { 1 3 f
p
Dos incgnitas de flujo, v(x,t) y p(x,t). Ecuaciones Diferenciales de Cant. De Mov. y Continuidad rigen cualquier flujo incompresible (Ecuaciones de Navier-Stokes) Ecuaciones EN DERIVADAS PARCIALES NO LINEALES. Las ms complejas de la fsica. Poqusimos flujos poseen solucin analtica (geometras sencillas y rgimen laminar). Flujo incompresible y completamente desarrollado en un conducto recto y de seccin arbitraria constante en rgimen laminar es uno de ellos.
MECNICA DE FLUIDOS
h h h k j+ fV = g [sen ( )i cos ( ) j + 0 k ]= i + z y x
(X) (Y) (Z)
2u 2u p 0 = + h + 2 + 2 y z x 0 = 0 = p +h y p +h z
MECNICA DE FLUIDOS
(Y) (Z)
0 =
p +h y
p 0 = + h z
H H (y ,z )
(X)
2u 2u dH + 2 + 2 0 = y dx z
H He dH = cte = s dx L
2 u 2u dH + 2 + 2 0 = y dx z
MECNICA DE FLUIDOS
dH 1 d du z = r dz r dr dr
CONDICIONES DE CONTORNO uZ(r=R)=Uwi y H/z Uw= 0 y H/z0 UW0 y H/z=0 UW0 y H/z0 Flujo Poiseuille Flujo de Couette Flujo Poiseuille+Couette
MECNICA DE FLUIDOS
er
dH 1 d du z = r dz r dr dr
CONDICIONES DE CONTORNO uZ(r=Ri)=Uwi, uZ(r=Re)=Uwe y H/z Uwi=Uwe=0 y H/z0 Uwi Uwe 0 y H/z=0 Uwi Uwe 0 y H/z0 Flujo Poiseuille Flujo de Couette Flujo de Poiseuille +Couette
Re
Ri
MECNICA DE FLUIDOS
dH d 2u = 2 dx dy
CONDICIONES DE CONTORNO u (y=0)=Uwi, u(y=h)=Uws y H/x Uwi=Uws=0 y H/z0 Uwi Uws 0 y H/z=0 Uwi Uws 0 y H/z0 Flujo Poiseuille Flujo de Couette Flujo de Poiseuille +Couette
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
DE & & = Q +W Dt
Siendo E la energa total del sistema, suma de su energa interna, cintica y potencial (Ek+Ep=Em).
E=(+Ek+Ep) =(+Em)
& Q potencia neta en forma de Calor (T)
& W potencia neta en forma de Trabajo (Fuerzas)
& & Q > 0,W > 0 Energa Entorno Sistema & & Q < 0,W < 0 Energa Sistema Entorno
Si el sistema que se considera es aquel que en el instante t est ocupando el V.C.
MECNICA DE FLUIDOS
VC
fV v dV + t v dS
SC
En las entradas y salidas de fluido al V.C. la potencia asociada a las fuerzas de superficie es:
En una entrada o en una salida la potencia asociada a las fuerzas viscosas se desprecia y la asociada a las fuerzas de presin: & & & W p s p v dS = p (v v SUP )dS p v SUP dS WF +WD Quedando: s s s
( )
)e
MECNICA DE FLUIDOS
( )
Cuando el flujo es incompresible (=cte) la 2a ley de la Termodinmica: ~ dUVC ~ & & & +USC Q =WL 0 dt 14243
( )
DU Dt
& + WF
& ) + (W
s s
& + WF
& ) W
s
( )
( )
MECNICA DE FLUIDOS
( )
( )
Caso Particular V.C. Fijo e Indeformable (vSUP 0) Flujo Incompresible (=cte) V.C. Con una entrada y una salida (qe=qs=q) Flujo Estacionario (v=v(x) d(Em)VC/dt0)
p p ek + ep + ek + ep + = w m w L s e
MECNICA DE FLUIDOS
p p hk + h + hk + h + = H m hL s e
Ecuacin de Bernoulli Ecuacin de Bernoulli
A la k+h+(p/) se le denomina Bernoulli del fluido (B) en la superficie y est expresado como una altura de columna de fluido [B]=L. Normalmente la altura de energa cintica promedio en una superficie k puede expresarse como:
v hk 2g
Siendo el coeficiente de correccin de la energa cintica.
MECNICA DE FLUIDOS
hL = Be Bs + H m 0 123
B
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
La constatacin de la existencia de los distintos regmenes de un flujo proviene de antiguo: Leonardo da Vinci (Estudio sobre el Agua). En el siglo XIX comenzaron los primeros estudios cientficos sobre el tema:
p p~v 1.75
G. H. L. Hagen (1839). Primeros indicios experimentales. Caida de presin en conductos largos de latn. Osborne Reynolds (1883). Pionero en el estudio de los regmenes de flujo.
p~v v
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
D v Re =
Siendo: D Dimetro de la tubera. v Velocidad media. Densidad del fluido. la viscosidad del fluido.
En cualquier flujo existen dos regmenes y la existencia de uno u otro depende de su nmero de Reynolds que viene dado por:
Re =
L U
MECNICA DE FLUIDOS
Re =
v v v = p + 2 v + fV fp + f + fv + fi = 0 + t r
Nmeros de Reynolds elevados (Reg. Turbulento): fi>>f. En el flujo predominan las fuerzas de inercia. Nmeros de Reynolds bajos (Reg. Laminar): fi<<f. En el flujo predominan las fuerzas viscosas.
MECNICA DE FLUIDOS
REGIMEN LAMINAR
Re=2300
TRANSICIN
Re=4000 Re=
REGIMEN TURBULENTO
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
div (v ) = 0 Dv = p + 2 v + fV Dt Las ecuaciones de Navier-Stokes rigen el flujo de un fluido en rgimen laminar y turbulento.
Las capacidades de clculo actuales no son capaces de resolverlas para cualquier tipo de flujo en rgimen turbulento. RESPUESTA:
uP(t)
Desde el punto de vista ingenieril NO es interesante conocer los valores instantneos de las variables de flujo sino su valor medio temporal. Las variables de flujo se descomponen en en un valor promedio y en una fluctuacin:
v = V + v p = P + p
UP t
u'P(t)
MECNICA DE FLUIDOS
(x, t ) = (x, t )
1 N = lim (n ) (x, t ) N N n =1
uP(t)
Siendo T un perodo de tiempo mayor que cualquier perodo significativo de las fluctuaciones y N un nmero de experimentos. Los promedios cumplen ciertas reglas como:
= ; = 0; 0 = x x
(t)
'(t)
UP t
1 (x ) = (x, t ) = (x, t ) dt T 0
u'P(t)
MECNICA DE FLUIDOS
Incluso en flujos sencillos (i.e.: flujo completamente desarrollado en un conducto) NO puede obtenerse analticamente el perfil de velocidades promedio
xy
dU = + ( xy )T y
Y 0
w
H d dU = y u v x dy
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
5.1 GENERALIDADES
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
0 (m) V 0 a b q01
1 c C1
q13
d B1
MECNICA DE FLUIDOS
j qij
Elementos Pasivos: En la relacin entre prdidas y caudal hL=hL(q) suele intervenir adems del caudal tambin otros parmetros caractersticos del fluido ( y ), la geometra y el material (rugosidad ) del elemento:
hL =h L (q, , , ,Geometra )
En el caso de una tubera hL=hL(q,,,,L,D).
MECNICA DE FLUIDOS
K=
hL hK
Siendo hk una altura de energa cintica caracterstica del elemento (entrada o salida). En el caso que existan dos velocidades medias es posible definir dos K segn la que se considere. Ambos estn relacionados (viAi= vjAj =q). La relacin hL=hL(q,,,,Geometra) es similar a la relacin entre K y unos parmetros adimensionales 1, 2,...,k obtenidos a partir de los dimensionales dependientes (q,,,,Geometra)
Tubera:
i q ij
j D, L y
K = K (Re , L D, D )
K = K ()
Vlvulas:
Grado de Apertura
MECNICA DE FLUIDOS
v2 K = q2 hL = K hK = K 2g 2g A 2 123 4 4
R
Pregunta: cmo se determina la relacin entre el coeficiente adimensional de prdidas o la ecuacin caracterstica de un elemento y el resto de parmetros. Respuesta: Es necesario resolver el flujo en el elemento (v y p): Anlisis Diferencial: Analticamente. Escasos casos en rgimen laminar. Numricamente. Mecnica de Fluidos Computacional (CFD). Experimentacin. Normalmente se combinan resultados experimentales y numricos.
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
B i (hf Bi = H i + (hk )i
)ij
= Bj
v i2 = Hi + i 2g v2 j 2g
B j = H j + (hk ) j = H j + j
H i H j = (hf )ij
MECNICA DE FLUIDOS
q ( j v j i v i ) =
Pw
= L A sen L tW dP + (pi p j ) A
vi=vj=q/A (Continuidad)
A (hi h j ) L tW dPW + (p i p j ) A = 0
Pw
Por la simetra la tensin cortante en la pared es igual en todo el permetro mojado Pw=D
(hi h j ) + pi
p j L Pw = tW A
Hi H j =
L Pw L 4 tW H i H j = tW A D
(hf )ij = L 4 tW
D
Para relacionar las prdidas de carga con el caudal (velocidad media) es necesario obtener una relacin entre ste ltimo y tw
MECNICA DE FLUIDOS
8 tW f = v2
(hf )ij
L 4 L v2 = tW = f D D 2g
Regimen Laminar (Re<2300): La relacin entre q y tw f(Re,/D) mediante la resolucin de la ecuacin diferencial del flujo:
1 d du H = r x r dr dr u (R ) = 0
r u (r ) = R 1 4 x R R H q = u (r ) 2 dr = 8 x 0 H
2
2 2q r u (r ) = 1 R 2 R H i H j H 8 R 4 = = q 4 x L R
tW =
du dr
=
r =R
4 8 q = v R R2 D
f (Re ) = 64 Re
(hf )ij
128 L q 4 g D
MECNICA DE FLUIDOS
Tuberas Lisas
MECNICA DE FLUIDOS
D f = 0.25 log 3 .7
MECNICA DE FLUIDOS
pi p j L Pw (hi h j ) + = tW A
Las prdidas de carga quedara como:
(hf )ij
L Pw tW A
DH =
4 A P w
Esta expresin semejante a la obtenida para una tubera circular.
(hf )ij
L 4 tW DH
MECNICA DE FLUIDOS
f =
8 tW v2
PREGUNTA: Son los resultados obtenidos para tuberas de seccin circular tiles para las no circulares, sustituyendo el dimetro por el dimetro hidrulico? RESPUESTA: Regimen Laminar (Re<2300): El factor de friccin de Darcy para tuberas de seccin circular no sigue la relacin f=64/ReH, Siendo ReH =v.DH/. En general en rgimen laminar la relacin del factor de friccin de Darcy es de la forma:
f (ReH ) = C ReH
C es un coeficiente (no tiene porque ser constante) diferente para cada tipo de seccin y que puede obtenerse integrando las ecuaciones diferenciales. (e.g.: seccin anular C=C(Ri/Re)).
Regimen Turbulento (Re>4000): Para tuberas de seccin no circular, el baco de Moody es vlido simplemente tomando en lugar del dimetro DH
MECNICA DE FLUIDOS
(hf )ij
Seccin No Circular
L v2 8f L 2 = 2 5 q ij = f D ij 2g D g ij 1 24 4 3 14 244 4 3
K ij Rij
(hf )ij
L v2 f L 2 = q ij = f D 2g 2g D A 2 H H 1 24 ij 4 3 1442443 4ij
K ij Rij
El factor de friccin de Darcy viene dado por: Rgimen Laminar (Re<2300): f=C/ReH (i.e.: seccin Circular C=64). Rgimen Turbulento (Re>4000): f=f(ReH,/DH). baco de Moody.
MECNICA DE FLUIDOS
(1) Conocido el caudal qij y el dimetro Dij, calcular la prdida de carga (hf)ij I. Calcular el nmero de Reynolds (Re)ij. II. Calcular el factor de friccin f=f(Re,/D)ij. III. Calcular la prdida de carga con la expresin de Darcy (hf)ij=Rij.q2ij.
(hf)ij=Hi-Hj (?)
MECNICA DE FLUIDOS
(hf)ij=Hi-Hj (C)
MECNICA DE FLUIDOS
(hf)ij=Hi-Hj (C)
MECNICA DE FLUIDOS
5.3 VLVULAS
MECNICA DE FLUIDOS
Q2
RED POLGONO
A
0
MECNICA DE FLUIDOS
Vlvula de mariposa
Vlvula de compuerta
MECNICA DE FLUIDOS
K () =
(hL )V (hK )V
(hL )V =
K () q 2 = KQ () q 2 2 2g 4 1 2A 4 3
KQ
caudal.
KQ() (dimensional [KQ]=Altura/Caudal2) se denomina Coeficiente de Prdidas referido al KQ() igual que K(), tiene un valor mnimo (KQ)0 K0 cuando la vlvula se halla
completamente abierta (=100%). A medida que se cierra ( disminuye) va aumentando hasta hacerse infinito cuando la vlvula se halla completamente cerrada (=0%). Para evitar trabajar con KQ() y K() , que toman valores tan elevados cuando la vlvula se halla casi cerrada, se introduce otro coeficiente (dimensional) denominado Coeficiente de Flujo KV():
KV () =
(hL )V
(hL )V
q2 = 2 KV ()
[KV]=Caudal/(Presin)1/2
MECNICA DE FLUIDOS
Se suele trabajar con otro coeficiente adimensional denominado Coeficiente de descarga Cd() que se define como: v Cd () = 2g (hL )V + v 2
Cd () =
1 + K ()
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
Ejemplo:
(hL )ij = (hL )ia + (hf )bc + (hf )dj (Hm )ij = (HB )ab (HT )cd
Nota: En las instalaciones hidrulicas suele despreciarse los trminos de energa cintica del Bernoulli.
qij + Qi = 0
j
Ejemplo:
q ij + q im q ki Q i = 0
MECNICA DE FLUIDOS
Nincgnitas=NL+NN=13
Q2 Q4 4 q41
A
q24
2 q12
hA (m)
V 0 a b
q01 1
Q0
MECNICA DE FLUIDOS
MECNICA DE FLUIDOS
[(h ) (H ) ] = 0
L ij m ij
qij + Qi = 0
j
MECNICA DE FLUIDOS
POBLACIN
RED 2
H 0 + H B (hL )V (hf
Q2 Q4 4 q41
A
q24
2 q12
hA (m)
V 0 a b
q01 1
Q0
q 01 + q 13 + q15 + q 12 q 41 = 0 q 24 + q 41 + Q 4 = 0 q 24 + q 25 q12 Q 2 = 0
MECNICA DE FLUIDOS
Ecuaciones formulacin en alturas: Despejando los caudales de NL Ecs. Bernoulli en las lneas de la red.
q ij
j
+ Qi = 0
Fij (H i , H j ) + Q i
j
=0
MECNICA DE FLUIDOS
Q2 Q4 4 q41
A
q15 = F15 (H1 ) q 41 = F41 (H1, H 4 ) q12 = F12 (H1, H 2 ) q24 = F24 (H 2 , H 4 ) q25 = F25 (H 2 )
3 Ecs. Continuidad en nodos de altura desconocida
q24
2 q12
hA (m)
V 0 a b
q01 1
Q0
F01 (H 1 ) + F13 (H 1 ) + F15 (H 1 ) + F12 (H 1, H 2 ) F41 (H 1, H 4 ) = 0 F24 (H 2 , H 4 ) + F41 (H 1, H 4 ) + Q 4 = 0 F24 (H 2 , H 4 ) + F25 (H 2 ) F12 (H 1, H 2 ) Q 2 = 0
ISBN