Cap Itulo 1 Fundamentos
Cap Itulo 1 Fundamentos
Fundamentos
1
2 CAPÍTULO 1. FUNDAMENTOS
1.1 Viscosidad
La viscosidad de un fluido es una propiedad importante en el estudio del flujo de fluidos.
La viscosidad es aquella propiedad del fluido mediante la cual éste ofrece resistencia al
esfuerzo cortante. La melaza y la brea, son ejemplos de los lı́quidos altamente viscosos; el
agua y el aire tienen viscosidades muy pequeñas. Para presiones ordinarias, la viscosidad
es independiente de la presión y depende únicamente de la temperatura. Para presiones
muy grandes, los gases y la mayorı́a de los lı́quidos muestran variaciones erráticas de la
viscosidad con la presión. Resolviendo para la viscosidad, de la Ecuación (1.1),
τ
µ=
du/dy
Estática de Fluidos
En este capı́tulo estudiaremos el caso en el que la velocidad del flujo es constante, es decir,
de acuerdo a la segunda ley de Newton, el caso cuando d(mv)/dt = 0. La suposición
fundamental para alcanzar esta condición es que no existe movimiento relativo entre capas
de fluido adyacentes y, consecuentemente, los esfuerzos cortantes son cero. Por consiguiente
se considera que fuerzas normales o fuerzas de presión son las únicas que actúan sobre las
superficies fluidas.
3
4 CAPÍTULO 2. ESTÁTICA DE FLUIDOS
respecto a un plano definido por los vectores direccionales unitarios s1 y s2 que son tangente
al punto de contacto. Por consiguiente se puede definir un sistema de coordenadas locales
en el plano local y el vector fuerza, F6 , se descompone en tres direcciones ortogonales,
∆Fn , ∆Fs1 y ∆Fs2 . La primera es la fuerza normal mientras que la segunda y tercera son
las fuerzas tangenciales.
Dado que el área es una cantidad vectorial que las áreas superficiales están rotando y
cambiando con el tiempo, es extremadamente difı́cil seguir la trayectoria del vector área.
Por consiguiente, conviene definir una representación intensiva conocida como esfuerzo.
Ocuerren dos clases generales de esfuerzos, normales y cortantes. Ambos se definen en un
sentido de lı́mites a medida que el área de contacto incremental tiende a cero, es decir,
∆Fn
σn = lim (2.1)
A→0 ∆A
∆Fs1
τss1 = lim (2.2)
A→0 ∆A
∆Fs2
τss2 = lim (2.3)
A→0 ∆A
Aquı́, el área especı́fica para al ejemplo es ∆A6 , σn , τss1 y τss2 son repectivamente, el
esfuerzo normal y los esfuerzos cortantes, tal como se define en coordenadas locales. La
dimensión de esfuerzo es fuerza por unidad de área.
Existe un número muy grande si no infinito de planos (áreas) que pueden aparecer en
un punto en el continuo fluido y, tal como se anotó, es aconsejable eliminar el proceso
analı́tico de guardar la información requerida para seguir la trayectoria de las coordenadas
locales del área. Por consiguiente es necesario definir el esfuerzo en términos de un sis-
tema ortogonal de fuerzas, referenciado a coordenadas globales y a una serie de planos
ortogonales asociados que pasen a través del origen. Al hacer esto, el estado de esfuerzos
puede definirse sin ambigüedades con un número mı́nimo de componentes. Para ello se re-
quiere de una derivación geométrica extensa, de la cual analizamos solamente los conceptos
básicos.
Con referencia a la Figura ??, el campo de esfuerzos del ejemplo en coordenadas locales
se proyecta a coordenadas globales siguiendo dos pasos. En primer lugar, el plano local
definido por s1 y s2 se proyecta en sus tres componentes ortogonales que pasan a través
del origen. Luego, cada uno de los tres esfuerzos locales σn , τss1 y τss2 se proecta en tres
componentes, en cada uno de los planos ortogonales que pasan a través del origen. Dado
que existen tres planos con tres sistemas de esfuerzos acompañantes, se requiere un total
de nueve componentes de esfuerzo para describir completamente el estado de esfuerzo en
un punto sobre una superficie arbitraria, en términos de un sistema global de coordenadas
2.1. FUERZA, ESUERZO Y PRESIÓN DE PUNTO 5
fijo. La ecuación ?? contiene el tensor total que contiene los nueve esfuerzos requeridos
σxx τyx τzx
τxy σyy τzy = τ̄¯ (2.4)
La notación aquı́ es como sigue: τyx es el esfuerzo cortante que actúa en la dirección x,
en un plano perpendicular al eje y ( es decir, el plano xz) y σyy es el esfuerzo normal que
actúa en la dirección y, en un plano perpendicular al eje y. Todos los vectores esfuerzos
en la Figura ?? están definidos como positivos. El esfuerzo normal positivo se dirige hacia
fuera del plano mientras que la regla del tornillo de la mano derecha se utiliza para definir
los valores positivos para los esfuerzos cortantes.
El promedio de los esfuerzos normales se conoce como esfuerzo volumétrico, σ, el cual
a su vez, se utiliza para definir la presión, p.
1
p = −σ̄ = − (σxx + σyy + σzz )
3
Mientras que los esfuerzos normales positivos se definen como positivos hacia afuera de la
superficie, la presión es positiva hacia el centro de masa de la superficie sobre la que actúa.
Un punto de un fluido en equilibrio tiene la misma presión en todas las direcciones.
Esto significa que un elemento δA de área pequeña, libre para rotar alrededor de su centro,
cuando se sumerge en un fluido en reposo, tendrá una fuerza de magnitud constante que
actúa en cualquiera de sus lados sin importar su orientación.
Para demostrar lo anterior, se toma un cuerpo libre en forma de cuña, de ancho uni-
tarioen el punto (x, y) de un fluido en reposo (ver Figura ??). Como no pueden existir
las fuerzas gravitacionales; entonces las ecuaciones de movimiento en las direcciones x y y
son, respectivamente,
X δx δy
Fx = px δy − ps δs sin θ = ρax = 0
2
y
X δx δy δx δy
Fy = py δx − ps δs cos θ − γ = ρay = 0
2 2
donde px , py y ps son las presiones promedio en las tres caras; γ es el peso especı́fico
del fluido; ρ es su densidad; y ax y ay son las aceleraciones. En el lı́mite que el tamaño del
cuerpo libre tienda a 0, permitiendo que la cara inclinada se aproxime a (x, y) mantieniendo
el mismo ángulo θ y utilizando las relaciones geométricas
δs sin θ = δy δs cos θ = δx
6 CAPÍTULO 2. ESTÁTICA DE FLUIDOS
la ecuación se simplifica a
δx δy
px δy − ps δy = 0 py δx − ps δx − γ =0
2
El último término de la segunda ecuación es un infinitesimal de orden superior que
puede ser despreciado. Si las ecuaciones se dividen por δy y δx, respectivamente, se
obtiene:
ps = px = py
Como θ es un ángulo arbitrario, esta ecuación prueba que la presión es igual en todas
las direcciones en un punto, en un fluido estático. A pesar de que la prueba fue desarrollada
para un caso en dos dimensiones, para uno en tres dimensiones se puede demosstrar que
las mismas ecuaciones de equilibrio para un pequeño tetraedro de fluido con tres caras en
los planos coordenados y una cuarta cara inclinada arbitrariamente.
Si el fluido se encuentra en movimiento, de tal manera que una capa se mueve respecto a
una adyacente, ocurren esfuerzos cortantes, y los esfuerzos normales, en general, no tienen
la misma magnitud en todas las direcciones en un punto. en este caso la presión se define
como el promedio de los esfuerzos de compresión normales, mutuamente perpendiculares
en un punto. En un fluido ficticio de viscosidad cero, es decir, un fluido sin fricción, no
ocurre ningún esfuerzo cortante para cualquier movimiento del fluido, de tal manera que
la presión es la misma en todas las direcciones.
donde δy
2
es la distancia desde el centro a una cara perpendicular a y. Sumando las fuerzas
que actúan en el elemento en la dirección y se obtiene
δp
δFy = − δx δy δz − γδx δy δz
δy
Para las direcciones x y z, debido a que no actúa ninguna fuerza de cuerpo,
δp δp
δFx = − δx δy δz δFx = − δx δy δz
δx δx
El vector de fuerza elemental δF está dado por
!
δp δp δp
δF = −iδFx + jδFy + kδFz = i + j + k δx δy δz − jγδx δy δz
δx δy δz
Si el elemento se reduce a un tamaño cero, después de dividir toda la ecuación por δx δy δz
la expresión se reescribe como
!
δF δ δ δ
=− i +i +i p − jγ lim δx δy δz → 0
δx δy δz δx δx δx
Ésta es la fuerza por unidad de volumen resultante en un punto, que es igual a cero para
un fluido en reposo. La cantidad entre paréntesis es el gradiente, conocido como ∇ (nabla).
δ δ δ
∇=i +j +k
δx δy δz
y el gradiente negativo de p es −∇p, es el vector de campo f de la fuerza de presión
superficial por unidad de volumen,
f = −∇p
Por consiguiente, la ley de fluido estático de variación de presión es
f − jγ = 0 (2.5)
δp δp δp
=0 = −γ =0
δx δy δz
Las ecuaciones diferenciales parciales, para variaciones en las direcciones horizontales, son
una forma de la ley de Pascal; éstas establecen que dos puntos en la misma elevación, en
la misma masa continua de fluido en reposo, tienen la misma presión. Como p es función
únicamente de y,
dp = −γ dy
Esta ecuación diferencial simple relaciona el cambio de presión con el peso especı́fico y el
cambio de elevación, y es válida tanto para fluidos comprensibles como imcompresibles.
Capı́tulo 3
Conservación de la Energı́a
En los problemas del flujo en los conductos hay tres tipos de energı́a que siempre hay
que tomar en consideracón. Tome un elemento de fluido, como el que se muestra en
la Figura ??, el cual es un ducto que puede estar a una cierta elevación z, tener cierta
velocidad v y una presión p. El elemento de fluido tendrá las siguientes formas de energı́a:
Energı́a potencial. Debido a su elevación la energı́a potencial del elemento con respecto
a algún nivel de referencia es:
PE = wz (3.1)
KE = wv 2 /2g (3.2)
FE = wp/γ (3.3)
Trabajo = pAL = pV
9
10 CAPÍTULO 3. CONSERVACIÓN DE LA ENERGÍA
w = γV
Trabajo = pV = pw/γ
wp2 wv22
E2 = + wz2 +
γ 2g
Si no se agrega energı́a al fluido o se pierde entre las secciones 1 y 2, entonces el principio
de conservación de la energı́a requiere que
A1 v1 = A2 v2
v2 = v1 (A1 /A2 )
Puesto que A1 < A2 , v2 debe ser menor que v1 . Y como la velocidad está al cuadrado en
el término correspondiente a la cabeza de velocidad, v22 /2g es mucho menor que v12 /2g.
Tı́picamente, cuando el tamaño de la sección se expande como lo hace en la Figura ??,
la cabeza de presión aumenta debido a que disminuye la cabeza de velocidad. Sin embargo,
el cabio real también se ve afectado por el cambio en la cabeza de elevación.
Cuando se escribe la ecuación de Bernoulli, es escencial que la presión en los dos puntos
de referencia se expresen ambas como presiones absolutas o como presiones manométricas.
Es decir, deben tener las dos la misma presión de referencia.
Aunque la ecuación de Bernoulli es aplicable a una gran cantidad de problemas prácticos,
existen algunas limitaciones que deben tenerse en cuenta con el fin de plicar la ecuación
de manera correcta.
1. Es válida solamente para fluidos incompresibles, puesto que el peso especı́fico del
fluido se tomó como el mismo en las dos secciones de interés.
2. No puede haber dispositivos mecánicos entre las dos secciones de interés que pudieran
agregar o eliminar energı́a del sistema, ya que la ecuación establece que la energı́a
total del fluido es constante.
hR = Pérdidas de energı́a por parte del sistema, debidas a fricción en los conductos, o
pérdidas menores debidas a la presencia de válvulas y conectores.
Los efectos de la transferencia de calor hacia o desde los ductos es también una fuente
de variación de la energı́a que no será considerada en este documento. El balance final de
las energı́as de la ecuación de Bernoulli está dado por,
0p v2
E = +z+ (3.5)
γ 2g
PA = hA W = hA γQ
en la que PA representa la potencia añadida al fluido, γ es el peso especı́fico del fluido que
fluye por la bomba y Q es la rapidez del flujo de volumen del fluido.
La unidad de potencia en el sistema Británico de Unidades es la lb − pie/s. Puesto
que es una práctica común referirse a la potencia en caballos de fuerza (horsepower- hp),
el factor de conversión requerido es
en la Ecuación (3.2) la energı́a agregada, hA , está expresada en pies de fluido que fluye en
el sistema. Entonces, expresando el peso especı́fico del fluido en lb/pie3 , y la rapidez del
flujo en pie3 /s obtendremos la rapidez de flujo del peso, γQ en lb/s. Finalmente, en la
ecuación de la potencia, PA = hA γQ, la potencia estará expresada en lb-pie/s.
El término de eficiencia se utiliza para denotar el cociente de la potencia transmitida
por la bomba al fluido entre la potencia suministrada a la bomba. Debido a las pérdidas
de energı́a ocasionadas por la friccón mecánica en los componentes de la bomba, la fricción
del fluido en la misma y la excesiva turbulencia del fluido que se forma en ella, no toda
la potencia suministrada a la bomba es transmitida al fluido. Entonces, el sı́mbolo eM es
utilizado para representar la eficiencia mecánica, tenemos:
PA
eM =
PI
donde PA es la potencia transmitida al fluido y PI es la potencia puesta en la bomba. El
valor de eM será siempre menor que 1.0.
flujo laminar, mientras que una rápido y caótico, se le conoce como flujo turbulento. Los
métodos que se utilizan para calcular la pérdida de energı́a es diferente para cada tipo de
flujo.
El tipo de flujo se puede predecir mediante el cálculo de un número sin dimensiones,
el número de Reynolds, que relaciona las variables más importantes que describen un
flujo: velocidad, v, diametro del conducto,D, densidad del fluido, ρ y la viscosidad, µ.
Osborne Reynolds fue el primero en demostrar que un flujo laminar o turbulento puede
ser predicho si se conoce la magnitud de un número adimensional, conocido ahora como
número de Reynolds (NR ), el cual se define como
vDρ
NR = (3.6)
µ
Los números de Reynolds grande, tı́picamente debido a una alta velocidad o a una baja
viscosidad, o a ambas, tienden a ser turbulentos. Aquellos fluidos que poseen una alta
viscosidad y/o que se mueven a bajas velocidades tendrán un número de Reynolds pequeño
y tenderán a ser laminares. Para aplicaciones prácticas en flujos de conductos, se tiene
que si el número de Reynolds para el flujo es menor que 2000, el flujo será laminar. Si
es mayor de 4000, se puede asumir que el flujo es turbulento. En el intervalo entre 2000
y 4000, es imposible predecir qué tipo de flujo existe, por consiguiente s a esta región se
le conoce como región crı́tica. Las aplicaciones tı́picas involucran flujos que se encuentran
bien colocados en el intervalo de los flujos laminares o en el de los turbulentos, por lo que
de esta manera evitamos la región de incertidumbre y ası́ cualquier dificultad que pudiera
ocasionar.
32µLv
hL = (3.7)
γD2
3.4. PÉRDIDAS POR FRICCIÓN 15
Diagrama de Moody
El diagrama de Moody se utiliza como una ayuda para determinar el valor del factor de
fricción, f para el flujo turbulento. Deben conocerse los valores del número de Reynolds
y de la rugosidad relativa. Por consiguiente, los datos básicos requeridos son el diámetro
interior del conducto, el material con que el conducto está hecho, la velocidad del flujo y
el tipo de fluido y su temperatura, con los cuales puede encontrar la viscosidad.
Los ejemplos siguientes ilustran el procedimiento para encontrar f .
Ejemplo 1 Determine el factor de fricción, f , si el agua a 160 deg F está fluyendo a 30.0
pies/s en un conducto de hierro forjado no recubierto cuyo diámetro interior es de 1 pulg.
Primero, se debe calcular el número de Reynolds para determinar si el flujo es laminar
o turbulento:
vD
NR =
ν
con D = 11plg=0.0833 pies, pero ν = 4.38 × 10−6 pies2 /s. Ahora tenemos:
(30.0)(0.0833)
NR = = 5.70 × 105
4.38 × 10−6
16 CAPÍTULO 3. CONSERVACIÓN DE LA ENERGÍA
As´ pues, el flujo es turbulento. Ahora se debe calcular la rugosidad relativa. Se considera
el número = 8 × 10−4 pies. Entonces la rugosidad relativa es
D
= 104
Observe que, con el fin de que D/ sea un cociente sin dimensiones, tanto D como deben
estar en las mismas unidades.
Los pasos finales del procedimiento son:
1. Localice el número de Reynolds sobre la abscisa del diagrama de Moody: NR =
5.70 × 105 .
El siguiente es un ejemplo tı́pico ilustrativo programado que ilustra una situación tı́pica
de tuberı́a para fluido. La pérdida de energı́a debido a la fricción debe calcularse como
parte de la solución.
Es cero. Puesto que los puntos A y B están ambos en el mismo conducto con un flujo
estable, las velocidades son iguales.
Esto nos lleva a hL , la pérdida de energı́a debido a la fricción entre A y B. ¿Cuál es el
primer paso?.
La obtención del número de Reynolds es el primer paso. Se debe determinar el tipo de
flujo, laminar o turbulento. Termine el cálculo del número de Reynolds antes de pasar al
siguiente panel.
El valor correcto es NR = 9.54 × 104 . En seguida se presenta cómo se encontró:
NR = vDρ/µ
1m3 /s
!
Q = (110L/min) = 1.83 × 10−3 m3 /s
60 000 L/min
Q 1.83 × 10−3 m3 /s
v = = = 0.932 m/s
A 1.963 × 10−3 m2
Para benceno a 50 ◦C, con una gravedad especı́fica de 0.86, encontramos que:
Entonces tenemos:
(0.932)(0.050)(860)
NR = = 9.54 × 104
4.2 × 10−4
Por lo tanto, el flujo es turbulento. ¿Qué relación debe utilizarse para calcular hL ?.
Para flujo turbulento debe usarse la ecuación de Darcy:
L v2
hL = f × ×
D 2g
Con el fin de utilizar el diagrama de Moody y encontrar el valor de f , ¿Se necesita el valor
de rugosidad relativa, D/?.
En este caso no. Debido a que el conducto es de plástico, la superficie interior es
lisa y podemos utilizar la curva del diagrama de Moody etiquetada como conductos lisos.
Calcular ahora el valor de f.
El resultado es f = 0.018. Ahora se puede terminar el cálculo de hL .
18 CAPÍTULO 3. CONSERVACIÓN DE LA ENERGÍA
f = 64/NR
Esta relación, está graficada en el diagrama de Moody como una lı́nea recta en el lado
izquierdo del Diagrama.
Desde luego, para los números de Reynolds desde 2000 hasta 4000, el flujo se encuentra
en la región crı́tica y es imposible predecir el valor de f .
Por encima del número de Reynolds de 4000, por lo general el flujo se conoce como
turbulento. Sin embargo, en esencia existen dos zonas de interés en este punto. Hacia
el lado derecho del diagrama, el flujo está en la zona de completa turbulencia. Se puede
observar que el valor de f no depende del número de Reynolds, sino solo de la rugosidad
relativa D/. En este intervalo, se aplica la siguiente fórmula:
1
√ = 2 log10 (3.7 D/) (3.8)
f
3.5. ECUACIONES DEL FACTOR DE FRICCIÓN 19
La frontera de esta zona es la lı́nea punteada que corre, por lo general, de la parte superior
izquierda a la parte inferior derecha del diagrama de Moody. La ecuación de esta lı́nea es:
1 NR
√ = (3.9)
f 200(D/)
La tercera zona del diagrama de Moody, que se conoce como zona de transición, se
encuentra entre la zona de completa turbulencia y la lı́nea que se identifica como conductos
lisos. Esta última lı́nea tiene la ecuación
√ !
1 NR f
√ = 2 log10 (3.10)
f 2.51
Siendo lisos, estos conductos no presentan irregularidades superficiales al flujo, de modo
que el factor de fricción sólo es función del número de Reynolds. Los conductos hechos
de vidrio o de cobre tienen un valor de rugosidad relativa que los acerca a la lı́nea de
conductos lisos.
En la zona de transición, el factor de fricción es función tanto del número de Reynolds
como de la rugosidad relativa. Una relación para el factor de fricción en esta zona:
!
1 1 2.51
√ = −2 log10 + √ (3.11)
f 3.7(D/) NR f
La Ecuación (3.12) produce valores para f que se encuentran entre el ±1.0% del valor de
los correspondientes a la ecuación de Colebrook, (3.11), dentro del intervalo de rugosidad
20 CAPÍTULO 3. CONSERVACIÓN DE LA ENERGÍA
relativa, D/, comprendido entre 1000 y 1 × 106 , y para números de Reynolds que van de
5 × 103 hasta 1 × 108 . Ésta es virtualmente la zona de turbulencia completa del diagrama
de Moody.
Bibliography
[1] Goldstein, S., Modern Development in Fluid Dynamics, Vol. II, pp. 676-680, Oxford
University Press, London. 1938.
[2] Moody, L.F. Friction Factors for Pipe Flow. Transactions of ASME 66(8):671-684.
Nueva York: American Society fo Mechanical Engineers. 1944.
[4] Streeter, V.L., E.B. Wylie and K.W. Bedford. Mecánica de Fluidos, 9a Ed. McGraw-
Hill
[5] Swamee, P.K., and A.K. Jain. ”Explicits Equations for Pipe-flow Problems.” Journal
of the Hydraulics Division, 102(HY5):657-664. Nueva York: American Society of Civil
Engineers. 1976.
21