Apuntes del curso ME446

Mecánica de los fluidos

Cristóbal E. Castro, Ricardo Fuentes

31 de mayo de 2022
2

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

Índice general

1. Propiedades de los fluidos 5

1.1. El fluido como un medio continuo . . . . . . . . . . . . . . . . 6
1.2. Propiedades termodinámica de los fluidos . . . . . . . . . . . 6

2. Estática de los fluidos 9

2.1. Presión es un escalar . . . . . . . . . . . . . . . . . . . . . . . 10
2.2. Derivación de la variación de presión en un fluido en reposo . 12
2.2.1. Hidrostática en un gas ideal . . . . . . . . . . . . . . . 13
2.2.2. Hidrostática en un líquido . . . . . . . . . . . . . . . . 14
2.2.3. Presión bajo una superficie libre . . . . . . . . . . . . 15
2.3. Presión manométrica . . . . . . . . . . . . . . . . . . . . . . . 16
2.4. Presión absoluta y presión relativa . . . . . . . . . . . . . . . 17
2.5. Fuerza sobre superficies sumergidas . . . . . . . . . . . . . . . 18
2.5.1. Superficies planas . . . . . . . . . . . . . . . . . . . . . 18
2.5.2. Superficies curvas . . . . . . . . . . . . . . . . . . . . . 24
2.6. Fuerza de empuje . . . . . . . . . . . . . . . . . . . . . . . . . 26
2.6.1. Estabilidad de un objeto sumergido . . . . . . . . . . 27
2.7. Distribución de presión en movimiento de cuerpo rígido . . . 27
2.7.1. Aceleración lineal uniforme . . . . . . . . . . . . . . . 28
2.7.2. Rotación de cuerpo rígido . . . . . . . . . . . . . . . . 30

3. Dinámica de los fluidos 31

3.1. Campo de velocidad y líneas de corriente . . . . . . . . . . . 32
3.2. Teorema de Transporte de Reynolds . . . . . . . . . . . . . . 32
3.3. Conservación de la masa . . . . . . . . . . . . . . . . . . . . . 33
3.4. Conservación del momento lineal . . . . . . . . . . . . . . . . 35
3.5. Conservación de la energía . . . . . . . . . . . . . . . . . . . . 35
3.6. Ecuación de Bernoulli . . . . . . . . . . . . . . . . . . . . . . 36
3.6.1. Formas de la ecuación de Bernoulli . . . . . . . . . . . 37
3.6.2. Representación gráfica de la ecuación de Bernoulli . . 40
3.6.3. Aplicaciones de la ecuación de Bernoulli . . . . . . . . 41

3
4 ÍNDICE GENERAL

4. Flujos en tuberías 47
4.1. Pérdidas de carga debido a la fricción . . . . . . . . . . . . . 48
4.1.1. Flujos laminares . . . . . . . . . . . . . . . . . . . . . 49
4.1.2. Flujos turbulentos . . . . . . . . . . . . . . . . . . . . 49
4.1.3. Ecuación de Colebrook . . . . . . . . . . . . . . . . . . 50
4.2. Típicos problemas en flujos en tuberías . . . . . . . . . . . . . 50
4.2.1. Cálculo de las pérdidas de carga . . . . . . . . . . . . 50
4.2.2. Cálculo del caudal . . . . . . . . . . . . . . . . . . . . 50
4.2.3. Cálculo del diámetro . . . . . . . . . . . . . . . . . . . 51
4.2.4. Cálculo mediante herramientas numéricas . . . . . . . 51
4.3. Pérdidas de carga por singularidades . . . . . . . . . . . . . . 51
4.4. Sistemas de tuberías . . . . . . . . . . . . . . . . . . . . . . . 51
4.4.1. Sistemas de tuberías en serie . . . . . . . . . . . . . . 52
4.4.2. Sistemas de tuberías en paralelo . . . . . . . . . . . . 52
4.4.3. Estanques conectados . . . . . . . . . . . . . . . . . . 53
4.4.4. Redes cerradas . . . . . . . . . . . . . . . . . . . . . . 54
4.5. Altura Neta Positiva en la Aspiración (NPSH) . . . . . . . . 55

5. Análisis dimensional 57
5.1. Teorema de Π Buckingham . . . . . . . . . . . . . . . . . . . 58

6. Flujo en canales abiertos 59

6.1. Flujo uniforme en canales de pendiente y sección constante . 61

7. Ejercicios 65
7.1. Flotabilidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
7.1.1. Problema 1 . . . . . . . . . . . . . . . . . . . . . . . . 65
7.1.2. Problema 2 . . . . . . . . . . . . . . . . . . . . . . . . 67
7.2. Fuerzas sobre superficies . . . . . . . . . . . . . . . . . . . . . 69
7.2.1. Problema 1 . . . . . . . . . . . . . . . . . . . . . . . . 69

8. Resolución numérica a problemas de mecánica de fluidos 71

8.1. Colebrook . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
8.1.1. Python . . . . . . . . . . . . . . . . . . . . . . . . . . 71
8.1.2. Excel . . . . . . . . . . . . . . . . . . . . . . . . . . . 72

9. Constantes 73
9.0.1. Ejemplos . . . . . . . . . . . . . . . . . . . . . . . . . 73

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

Capítulo 1

Propiedades de los fluidos

El concepto de fluido en mecánica abarca los gases y líquidos, separando

en esta clasificación a los sólidos. Esta diferenciación es clara para la gran
mayoría de las personas, sin embargo es muy difícil de poner en palabras
para quien no se ha formado en el estudio de los fluidos. Un ejercicio muy
interesante es preguntar a un niño si alguna de las substancias que obser-
vamos en casa es o no es un fluido. Este mismo ejercicio se puede hacer
con los alumnos de pregrado y, con mucha probabilidad, lograrán acertar la
respuesta. Un elemento central en la distinción entre un fluido y un sólido es
la capacidad de resistir un esfuerzo de corte. Por ejemplo, un sólido podrá
resistir un esfuerzo de corte con una deformación estática, por el contrario,
un fluido no.

Un fluido se puede definir formalmente como: “Es una sustancia que

se deforma continuamente cuando es sometido a un esfuerzo de
corte, sin importar la magnitud del esfuerzo”.

Como se mencionó anteriormente, en mecánica distinguimos claramente

los sólidos de los fluidos. En la categoría fluidos podemos encontrar los líqui-
dos y los gases. Nuevamente una persona común podrá distinguir los líquidos
de los gases pero una persona entrenada en la mecánica de fluidos debe po-
der formalizar esta distinción en palabras. Un líquido está formado por un
conjunto de moléculas que se encuentran muy cercanas y con fuerzas cohe-
sivas fuertes. Estas fuerzas mantienen las moléculas cercanas unas a otras
y tienden a formar una superficie libre cuando están sometidas a un campo
gravitacional. Por el contrario, un gas tiene sus moléculas más separadas,
debilitando las fuerzas cohesivas, y se expandirá hasta llenar el volumen dis-
ponible. Si dejamos un gas libre para moverse, este formará una atmósfera
en equilibrio hidrostático. En general, gases no formarán una superficie libre
pero si será relevante las fuerzas de flotabilidad que se desarrollen.

5
6 CAPÍTULO 1. PROPIEDADES DE LOS FLUIDOS

1.1. El fluido como un medio continuo

Si observamos un fluidos, gas o líquido, a escala molecular podemos ver
que consta de un conjunto de moléculas, las cuales están más cerca en los
líquidos y más separadas en los gases. En ambos casos, la distancia entre mo-
léculas es grande comparada al tamaño de las moléculas. Esta distribución
espacial nos obliga a identificar escalas mínimas necesarias para hablar de
una sustancia continua, la cual puede ser definido con las técnicas clásicas de
la física. Para fluidos esa escala mínima está en el orden de 10−9 mm3 [4]. En
este curso consideraremos dimensiones físicas mucho mayores a éste volumen
mínimo lo que nos permitirá hablar de densidad en un punto y aceptar la
continuidad de las propiedades de los fluidos.

1.2. Propiedades termodinámica de los fluidos

Los fluidos, al igual que todas las sustancias que conocemos se comportan
según las Leyes de la Termodinámica. De acuerdo a esto, todo lo aprendido
en el ramo Termodinámica 1 se aplica a este curso. La variables termodiná-
micas de mayor interés en este curso son: densidad ρ, viscosidad µ, presión
p y temperatura T .
Densidad: Identificada por ρ la letra griega minúscula rho, representa
la masa por unidad de volumen y es la inversa del volumen específico v.
Las unidades de densidad son kg/m3 . La densidad es dependiente de la
temperatura y de la presión: en los gases es altamente de pendiente de ambas
mientras que en los líquidos la dependencia de la temperatura es menor y
prácticamente independiente de la presión. Más adelante se tratará el tema
de la incompresibilidad de los líquidos y el módulo a la compresibilidad. En
condiciones normales (20◦ C y 1 atm) el gas más liviano es el hidrógeno ρH =
0,0838 kg/m3 y el líquido más pesado es el mercurio ρHg = 13580 kg/m3 .
Presión: Identificada por p, representa el esfuerzo de compresión a la
que está sometido el fluido. En termodinámica el valor de la presión es
fundamental para caracterizar un estado, mientras que en mecánica de los
fluidos es más relevante el gradiente de la presión, el cual induce un flujo.
Las unidades de la presión es Pascal (Pa), fuerza dividido por área.
Temperatura: La temperatura es la medida del nivel de energía inter-
na del fluido. Si existen grandes diferencias de temperatura en un fluido los
procesos de transferencia de calor pueden ser importantes. En este curso
no veremos la influencia de la transferencia de calor en los flujos. La tem-
peratura se representa con la escala absoluta Kelvin (K) o relativa Celsius
(◦ C).
Viscosidad: La viscosidad relaciona los esfuerzos de corte que actúan
sobre el fluido y su rapidez de deformación. La viscosidad dinámica µ se
identificada por la letra minúscula griega mu y sus unidades son N s/m2 =

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 1. PROPIEDADES DE LOS FLUIDOS 7

kg/m s. La viscosidad cinemática ν = µ/ρ se identifica por la letra griega

nu y sus unidades son m2 /s. La viscosidad juega un rol fundamental en la
descripción del balance de energía para la dinámica de los fluidos. Uno de
los experimentos más largos de la historia tienen a la viscosidad como un
elemento central en el Pitch Drop experiment.
La viscosidad es la variable de proporcionalidad entre el esfuerzo de
corte τ al cual está sometido el fluido y la rapidez de deformación del fluido
du/dy. Para fluidos Newtonianos la viscosidad no depende del esfuerzo de
corte luego se dice que es contante. Una de las ecuaciones más conocidas de
la mecánica de fluidos es la Ley de Viscosidad de Newton

du
τ =µ (1.1)
dy
Peso específico: El peso específico es el peso por unidad de volumen y
se representa por la letra griega minúscula gamma (γ) y se evalúa γ = ρ g.
Las unidades son N/m3 o fuerza por unidad de volumen. Notar que se utiliza
la gravedad para calcular γ.
Gravedad específica: La gravedad específica SG es una variable adi-
mensional que relaciona la densidad del fluido de interés y uno de refe-
rencia. En general para líquidos se utiliza el agua y para gases el aire.
SGgas = ρgas /ρaire , SGliq = ρliq /ρagua .
Tensión superficial: Los líquidos al interactuar con otro fluido genera-
rán una interfaz en su superficie que encontrará un equilibrio de fuerzas y
una geometría particular. La tensión superficial es la fuerza responsable de
este equilibrio y generá las gotas y meniscos. Las moléculas de los líquidos
tienen una distribución espacial que varía según se encuentren al interior
del líquido o en su superficie. Las moléculas cercanas a la superficie se en-
cuentran en un estado de tensiones y el origen de la fuerza generadora de la
tensión superficial radica en complejas reacciones químicas[1].
Coeficiente de compresibilidad: El coeficiente de compresibilidad K
de un fluido nos indica el el cambio de volumen respecto a un cambio de pre-
sión. En los gases, un aumento de presión genera una inmediata disminución
de volumen variando la densidad acordemente. En los líquidos la variación
de volumen es extremadamente pequeña, considerándose así como un fluido
incompresible. El coeficiente de compresibilidad o también llamado modulo
elástico volumétrico del fluido se calcula como K = -V (dp / dV)T donde V
es el volumen, p la presión y la derivada de presión con respecto al volumen
es a temperatura constante[3]. En forma discreta podemos escribir K = -V
∆p / ∆V. Ejemplos de K son KH2O = 2,2GPa, Kaire = 0,14MPa.
Condición de no-deslizamiento: La condición de no-deslizamiento o
non-slip condition en ingles, no es una propiedad de los fluidos pero sí es
consecuencia de las viscosidad. Cuando un fluido se encuentra en contacto
con una superficie, la velocidad de las moléculas de fluido en contacto con
la superficie es cero con respecto a esa superficie. Esta condición se utilizará

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

8 CAPÍTULO 1. PROPIEDADES DE LOS FLUIDOS

en numerosos modelos que describirán el comportamiento de los fluidos. Un

ejemplo de este fenómeno se puede ver en el siguiente video.

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

Capítulo 2

Estática de los fluidos

Muchos casos de estudio involucran fluidos en reposo donde la variable de

mayor interés es presión que se traduce en una distribución de fuerzas que se
ejercen sobre superficies. Producto de estas fuerzas encontramos las fuerzas
de flotación. La condición de fluido en reposos la indicamos con V ~ = ~0
y se llama hidrostática. Esta indicación hace referencia a un sistema de
referencia inercial, por ejemplo un avión que viaja a velocidad constante. El
lector puede profundizar en sistemas de referencia no inercial en [2].
En la condición hidrostática, los cambios de presión se deben al peso
del fluido debido a la gravedad así como también cuando existe una fuerza
externa debido a una aceleración. La hidrostática es fundamental en inge-
niería en el área de la hidráulica entendiendo las presiones y fuerzas sobre
superficies, como también en geofísica y astrofísica. Nos permite estudiar
el manto terrestre y la atmósfera, la presión de la sangre y muchos otros
fenómenos de interés para la sociedad.
Para el estudio de la hidrostática desarrollaremos expresión que se sim-
plificarán a medida que apliquemos algunas hipótesis fundamentales. En ge-
neral utilizaremos las herramientas de la mecánica clásica para desarrollar
estas ecuaciones las cuales son:

Conservación de la masa.

Conservación del momento. Segunda ley de Newton.

P
~ = m~
F a.

Estos principios fundamentales se aplicarán a un elemento de fluido el

cual puede tener variadas formas según el análisis y tiene dimensiones dis-
cretas. Por ejemplo podríamos estar interesados en entender la fuerzas sobre
un elemento de fluido sumergido como se observa en la Fig. 2.1. En esta fi-
gura observamos dos fluidos, los cuales son muy comunes en nuestras vidas,
y que usaremos regularmente como ejemplo sin reducir la generalidad del
análisis. Al existir dos fluidos simultáneamente y encontrarse bajo los efectos
de un campo de gravedad, estos se distribuyen de modo estratificado donde

9
10 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS

Figura 2.1: Se observa un objeto sumergido en agua.

a la interfaz entre ellos se le designa la superficie libre y se identifica con el

triangulo invertido de la figura. Esta terminología se refiere a que el aire de-
fine la atmósfera de nuestro planeta y la superficie expuesta está libremente
conectada con la atmósfera y no confinada por una pared. Todos los fenó-
menos de oleaje en el mar se caracterizan por ser situaciones en las cuales la
superficie libre juega un rol fundamental. En general hablaremos que en la
superficie libre tenemos una presión igual a la presión atmosférica, a menos
que se diga lo contrario. Una propiedad de la superficie libre cuando está en
equilibrio hidrostático es que es perpendicular al vector de aceleración, en
el caso de la Fig. 2.1, es perpendicular a la gravedad y por lo tanto define
una superficie horizontal.
Debido a la aceleración de gravedad y consecuentemente el peso del flui-
do, el objeto sumergido estará sometido a un campo de presión y a fuerzas
sobre sus superficies, las cuales podremos caracterizar completamente con
lo que aprendamos en este capítulo. Para efecto del análisis, podemos con-
siderar que el objeto sumergido es un Elemento de Fluido.

2.1. Presión es un escalar

Imaginemos que el objeto sumergido o elemento de fluido de la Fig. 2.1
tiene forma de cuña como se aprecia en la Fig. 2.2 donde la dirección x
es perpendicular al plano y − z. Este elemento de fluido tiene dimensiones
características dx , dy , dz , ds y un volumen total dV . En esta figura vemos
el elemento de fluido sometido a fuerzas en sus tres caras producto de la
presión, Fy , Fz y Fs y una fuerza de volumen W producto de su propio
peso. Si el elemento no se está acelerando, de acuerdo a la Segunda ley de
Newton
P
~ = m~0 podemos realizar un balance de fuerzas en la dirección
F
z obteniendo

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS 11

Figura 2.2: Elemento de fluido con forma de cuña que nos permite entender
el concepto de presión como un escalar.

Fz − Fs cos θ − W = 0, (2.1)
donde

Fz = pz dx dy, es la fuerza en la dirección z,

Fs = ps dx ds, es la fuerza en la dirección s y

W = ρ g dV , es el peso debido a la densidad, gravedad y volumen del

elemento de fluido.

Hasta ahora hemos asumido que la presión en las distintas caras es dis-
tinta y llamamos pz a la presión en la cara inferior y ps a la presión en la
cara diagonal. Reemplazando las fuerzas en la ecuación anterior obtenemos

pz dx dy − ps dx ds cos θ − ρ g dV = 0 . (2.2)
Reemplazando el volumen del elemento de fluido dV = dx dy2dz , utilizan-
do la relación geométrica dy = ds cos θ y luego de algunos pasos algebraicos
concluimos

dz
pz − ps − ρg = 0. (2.3)
2
Siguiendo los mismos pasos pero en la dirección y, el balance de fuerzas
nos entrega

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

12 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS

py − ps = 0 . (2.4)

A partir de (2.3) y (2.4) podemos concluir que no hay variación de presión

en un plano horizontal py = ps y que la variación de presión en la dirección
vertical depende del tamaño dz del elemento pz = ps + ρg dz 2 .
Ahora, si consideramos que el elemento de fluido es un elemento pequeño,
donde dx = dy = dz → 0 para cualquier angulo θ, este resultado se puede
expandir en palabras indicando que la presión en un punto es un escalar, no
tiene dirección, y es la misma en cualquier dirección que miremos, es decir
py = pz = ps , para cualquier ángulo θ, p es un escalar.

2.2. Derivación de la variación de presión en un

fluido en reposo
Consideremos el elemento de fluido de la Fig 2.3. Asumamos que co-
nocemos la presión en el centro de este elemento y la identificamos con la
letra p. Siendo la presión una función continua y diferenciable en el espa-
cio p = p(x, y, z) podemos estimar la presión en el costado derecho pR e
izquierdo pL de nuestro elemento de fluido, utilizando las series de Taylor
despreciando términos de orden 2 o superiores, de la siguiente forma

dy ∂p
pR ≈ p + (2.5)
2 ∂y
dy ∂p
pL ≈ p− (2.6)
2 ∂y

La presiones sobre superficie generan fuerzas luego podemos escribir que

Figura 2.3: Elemento de fluido con forma de hexaedro que nos permite de-
rivar la ecuación de la hidrostática. El elemento de fluido tiene dimensiones
dx, dy, dz y está sometido a la fuerza de gravedad generando un peso W .

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS 13

dy ∂p
 
FR = pR dx dz = p+ dx dz (2.7)
2 ∂y
dy ∂p
 
FL = pL dx dz = p− dx dz (2.8)
2 ∂y
Ahora podemos realizar un balance de fuerzas en la dirección y y obtener

FL − FR = m ay , (2.9)
dy ∂p dy ∂p
   
p− dx dz − p + dx dz = ρ dV ay . (2.10)
2 ∂y 2 ∂y
Es importante entender que el término m ay se refiere a las fuerzas ex-
ternas que percibe el elemento de fluido en su conjunto, por ejemplo la
aceleración de un camión aljibe.
Reemplazando dV = dx dy dz y repitiendo el análisis en las tres direc-
ciones considerando el peso en la dirección z obtenemos

∂p
− = ρ ax (2.11)
∂x
∂p
− = ρ ay (2.12)
∂y
∂p
− − ρ g = ρ az (2.13)
∂z
que se puede escribir en forma vectorial

−∇p − ρ g k̂ = ρ~
a, (2.14)
donde el operador diferencial nabla ∇ es el vector ∇ = [∂/∂x, ∂/∂y, ∂/∂z] ⊺ ,
k̂ es el vector unitario en la dirección z con k̂ = [0, 0, 1]⊺ y ~ a es el vector
de aceleraciones ~a = [ax , ay , az ]⊺ . Estos tres vectores son vectores columnas
y permiten escribir de forma vectorial el balance por unidad de volumen
(2.14).
Es importante reconocer los supuestos utilizados en la derivación de
(2.14) y las unidades que considera. En este desarrollo se han despreciado las
fuerzas viscosas las cuales no son importantes en un análisis de hidrostática
y las unidades consideran fuerzas divididos por unidad de volumen.
Ejercicio: Analizar las unidades de cada término de la ecuación (2.14)

2.2.1. Hidrostática en un gas ideal

Para entender la hidrostática en un gas ideal necesitamos retroceder al
balance de fuerza por unidad de volumen planteado en (2.14) y si asumimos
que no hay aceleración o fuerza externa escribimos

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

14 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS

−∇p − ρ g k̂ = ~0 (2.15)
Planteando esta ecuación en la dirección vertical tenemos

∂p
= −ρg . (2.16)
∂z
Si consideramos que el fluido es un gas ideal [4] podemos utilizar la ecuación
de gas ideal pV = nRT donde p es la presión absoluta, V el volumen, n el
número de moles, R la constante universal de los gases y T la temperatura
absoluta. Otro modo de escribir esta ecuación es p = ρRT considerando
que R = R/M y M es el peso o masa molecular del gas. Reemplazando la
ecuación de gas ideal en (2.16) tenemos

∂p p
=− g. (2.17)
∂z RT
Si consideramos que la temperatura es una constante, utilizando separación
de variables e integrando, obtenemos
g
ln p = − z + ln c , (2.18)
RT
donde c es una constante de integración que al evaluarla en z = 0 obtenemos
c = po . Finalmente la ecuación de hidrostática aplicada a un gas ideal y
temperatura constante se puede escribir como

gz
 
p(z) = po exp − . (2.19)
RT
Ejercicio: Desarrollar las ecuaciones de hidrostática para un gas ideal
donde la temperatura varía linealmente con la altura a una razón de 0,65
grados Celsius cada 100 metros. Solución

2.2.2. Hidrostática en un líquido

Si consideramos la ecuación de hidrostática (2.14) aplicada a un líqui-
do donde asumimos que ρ = constante y no hay aceleraciones o fuerzas
externas, podemos plantear en una dirección horizontal

∂p
− = 0, (2.20)
∂x
∂p
− = 0. (2.21)
∂y
Ambas ecuaciones implican que la presión no es función de la posición
horizontal luego p = p(z) varía solamente en la dirección vertical producto
de la fuerza de gravedad y consecuentemente del peso del líquido. La tercera
componente de (2.14) confirma esta dependencia y se escribe

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS 15

∂p
− − ρg = 0 (2.22)
∂z
que también podemos escribir como
Z
p(z) = −ρ g dz + C , (2.23)

donde los límites y la constante C de integración dependerán de las con-

diciones particulares del problema. Ecuación (2.23) es bastante general en
el sentido que aún considera la gravedad variable. A continuación veremos
un caso más específico.

2.2.3. Presión bajo una superficie libre

Imaginemos la situación de una piscina donde se verifica una superficie
libre como se muestra en la Fig. 2.4.

Figura 2.4: Contenedor de un fluido líquido donde deseamos calcular la pre-

sión a profundidad z = −h considerando una superficie libre expuesta a la
presión atmosférica.

Si buscamos evaluar la presión a una profundidad h podemos evaluar

(2.23) entre los límites z = 0 que corresponde a presión atmosférica y z = −h
que es el valor a encontrar, asumiendo que la gravedad es constante, situación
válida en la mayoría de los casos de interés en ingeniería.

Z z=−h
p = −ρ g dz
z=0
p = −ρ g (−h − 0)
p = ρgh

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

16 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS

En estricto rigor matemático, la expresión correcta es

Z p(z=−h) Z z=−h
dp = −ρ g dz (2.24)
p(z=0) z=0

lo que significa

p(z = −h) − p(z = 0) = −ρ g (−h − 0) (2.25)

que, de acuerdo al problema presentado en Fig. 2.4 donde la presión en z = 0

es la presión atmosférica p(z = 0) = patm , finalmente se escribe

p(z = −h) = patm + ρ g h (2.26)

Ecuación (2.26) es fundamental para analizar problemas de manome-

tría, donde la presión aumenta si descendemos en el fluido y disminuye si
ascendemos en el fluido. Es importante que el estudiante reconozca que esta
ecuación se puede aplicar entre dos puntos cualquiera de un fluido con la res-
tricción que la densidad permanezca constante en ese tramo. La elección del
punto de referencia, que en (2.26) se consideró la superficie libre con z = 0,
puede cambiar de forma arbitraria, solamente es necesario ser consecuente
con las elección.
Ejercicio: Aplicar (2.26) considerando el punto de referencia el fondo
de la piscina donde asumo una presión conocida p0 y evaluar la presión un
metro más arriba.

2.3. Presión manométrica

En esta sección veremos los problemas típicos de manometría, muy uti-
lizados en el estudio de la mecánica de los fluidos.
Ecuación (2.26) nos muestra como se verifica un cambio de presión frente
a un cambio de altura. Esta relación nos indica que, usando uno o más
líquidos para los cuales conocemos su densidad, podemos medir diferencia de
presión entre dos puntos. Un aparato que logra este objetivo se le denomina
un manómetro.
Los manómetros más comunes son tubos que contienen fluidos de den-
sidad conocida, y mediante la medición de la altura del fluido, es posible
estimar una presión. Por ejemplo en la Fig. 2.5 vemos un recipiente lleno de
un fluido de densidad conocida ρ que, por medio de un tubo, esta expuesto a
la presión atmosférica. La idea de este instrumento es monitorear la presión
en A mediante la columna de fluido de dimensión h. Si h aumenta quiere de-
cir que la presión en A aumenta, mientras que, si h disminuye, la presión en
A disminuye. Esto se formaliza mediante la ecuación de hidrostática (2.26)
donde

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS 17

Figura 2.5: Manómetro de tubo que permite el monitoreo de la presión en

punto A con respecto a la presión de referencia externa que en este caso es
la presión atmosférica.

PA = Patm + ρ g h . (2.27)
Para que este mecanismo funcione, es fundamental que la densidad del
fluido de trabajo sea mayor a la densidad del aire a la cual está expuesto, es
decir ρ > ρaire . Si esta condición no se cumple, el fluido de trabajo dejaría
el tubo y se dispersaría en el ambiente.
Otro instrumento comúnmente usado para medir presiones manométri-
cas es el manómetro de Bourdon.

2.4. Presión absoluta y presión relativa

El valor de la presión no es siempre un indicativo para identificar fuerzas
totales ejercidas sobre objetos o superficies sumergidas. Como vimos ante-
riormente en (2.14), es el gradiente de presión quién genera una fuerza. Es
común utilizar valores de referencia de la presión para entender la situa-
ción en que nos encontramos. El valor de referencia natural para nosotros
es la presión atmosférica. En la Fig. 2.6 vemos una representación gráfica
del concepto de presión absoluta y presión relativa. Si consideramos la pre-
sión atmosférica como la referencia, podemos encontrar presiones positivas
y negativas con respecto a la presión atmosféricas. Siguiendo este valor de
referencia, una presión relativa igual a 0 Pa es equivalente a un valor 101,3
kPa de presión absoluta.
La unidades de presión y sus equivalencias más comunes se pueden ob-
servar en la tabla
En el caso de referir la presión como altura de un líquido es necesario
considerar las condiciones estándar debido a la fuerte dependencia de la
densidad y la temperatura. En esta tabla

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

18 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS

Figura 2.6: Representación gráfica del concepto de presión absoluta y rela-

tiva.

Cuadro 2.1: Unidades de presión y equivalencias.

atm Pa kPa psi mmHg inHg mca

1,0 1,013 × 105 101,3 14,7 760 29,9 10,33

2.5. Fuerza sobre superficies sumergidas

Cuando tenemos una superficie sumergida en un fluido, la presión ejer-
cida por el fluido produce fuerzas distribuidas sobre la superficie. Conocer y
entender estas fuerzas nos permite calcular estanques de agua, compuertas,
represas y otras obras hidráulicas. Una estrategia para entender estas fuer-
zas es estudiarlas en casos simples para luego extender el análisis a casos
más complejos. Iniciaremos por estudiar el caso de una superficie plana su-
mergida y estudiaremos las fuerzas distribuidas y la resultante equivalente.
Con esta información extenderemos el análisis a una superficie curva.

2.5.1. Superficies planas

Imaginemos una superficie plana sumergida que forma un angulo θ con
la superficie libre como se observa en la Fig. 2.7. En esta figura podemos
ver el eje y en la dirección del plano de la superficie y el eje x en la in-
tersección de la proyección del plano de la superficie con la superficie libre.
Adicinalmente se muestra una vista del pano x − y abatida para observar
la forma de la superficie sumergida. Producto del peso del fluido y conse-
cuente presión hidrostática, que aumenta con la profundidad, se genera una

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS 19

Figura 2.7: Esquema de una superficie plana sumergida en un fluido. El eje y

tiene origen en la superficie libre a un ángulo θ con respecto a esta. El eje x
sale del papel. En la parte derecha-abajo, se muestra la superficie sumergida
en el plano abatido x − y.

fuerza distribuida sobre esta superficie. Nuestro interés es evaluar cual es la

fuerza distribuida F que actúa sobre la superficie, la fuerza resultante FR y
la posición de acción de esta fuerza resultante yR .

Fuerza resultante FR
Para encontrar este resultado utilizaremos la ecuación de hidrostática
para evaluar la presión a profundidad h, la cual es p = p0 + ρ g h donde p0 es
la presión en la superficie libre. Si despreciamos p0 = 0 obtenemos p = ρ g h.
Por geometría podemos decir que sin(θ) = h/y con y la distancia desde el
origen o hasta el diferencial de área de interés. Luego podemos escribir

p = ρ g y sin(θ) . (2.28)
Esta presión actuando sobre el elemento de área dA generará un diferen-
cial de fuerza

dF = ρ g y sin(θ) dA . (2.29)
Integrando este diferencial de fuerza sobre el área total de la superficie se
obtiene la fuerza total actuando
Z Z
F = ρ g y sin(θ) dA = ρ g sin(θ) y dA (2.30)
A A

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

20 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS

donde hemos asumido que la densidad del fluido no cambia con la profun-
didad, la superficie es plana y estática y la gravedad es una constante. Si
recordamos ahora el concepto de primer momento de área con respecto al
eje x podemos evaluar la posición del centroide yc de la superficie sumergida
como
1
Z
yc =
y dA . (2.31)
A A
Combinando esta ecuación con (2.30) obtenemos

F = ρ g sin(θ) A yc . (2.32)
Nuevamente por argumentos geométricos podemos escribir sin(θ) = hc /yc
lo que nos permite definir la fuerza total en base a la presión en el centroide
de la superficie pc = ρ g hc obteniendo

F = pc A . (2.33)
Ecuación (2.33) nos permite cuantificar la magnitud de la fuerza dis-
tribuida actuando sobre la superficie sumergida. La siguiente pregunta es
evaluar la posición de acción de la fuerza equivalente FR = pc A de tal modo
que genere el mismo momento sobre el eje x. En la Fig. 2.8 vemos el esquema
de la fuerza resultante FR y su posición de acción yR que calcularemos de
la siguiente manera.

Posición de acción yR de fuerza resultante

Para encontrar yR debemos igualar el momento generado por la fuerza
distribuida F con el momento generado por la fuerza resultante FR . El mo-
mento producido por la fuerza distribuida se calcula integrando el diferencial
de momento dM = y dF de la siguiente forma

Z Z Z Z
2
M= dM = y dF = ρ g sin(θ) y dA = ρ g sin(θ) y 2 dA (2.34)
A A

Si recordamos ahora el concepto de momento de inercia1 Ix = 2 dA

R
Ay
con respecto a un eje x podemos escribir

M = ρ g sin(θ)Ix , (2.35)
Considerando que la superficie podría estar a una profundidad arbitraria
utilizaremos el Teorema de los Ejes Paralelos o Teorema de Steiner el que
nos indica

Ix = Ixc + A yc2 . (2.36)

1
Cuán difícil es romper la inercia a la rotación

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS 21

Figura 2.8: La fuerza distribuida es equivalente en magnitud a la fuerza

resultante FR la cual actúa en la posición yR . Ese punto se llama centro de
presión de la superficie.

Este teorema nos permite calcular el momento de inercia con respecto

a un aje arbitrario x a partir de el momento de inercia con respecto al eje
que pasa por el centroide de la superficie más la distancia al cuadrado del
centroide al eje de interés como se observa en la Fig. 2.9. Los valores de Ixc
están tabulados para las principales geometrías.

Figura 2.9: El Teorema de Steiner nos permite calcular el momento de inercia

con respecto a un eje paralelo al eje que pasa por el centroide de la superficie
a distancia yc .

Retomando el momento de la fuerza distribuida (2.35) más el teorema

de los ejes paralelos (2.36) e igualando al momento producido por la fuerza
resultante podemos escribir

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

22 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS

M = FR yR
ρ g sin(θ)(Ixc + A yc2 ) = ρ g sin(θ) A yc yR
(Ixc + A yc2 ) = A yc yR
Ixc
yR = yc + (2.37)
A yc

Ecuación (2.37) puede ser utilizada bajo el supuesto que la presión en

superficie libre p0 actúa en ambos lados de la superficie sumergida o cuando
es despreciable. En caso contrario la posición de acción es

Ixc
yR = yc + 
p0
 (2.38)
A yc + ρ g sin θ

Desarrollo alternativo

Ecuación (2.37) es válida de acuerdo al sistema coordenado mostrado

en la Fig. 2.7. La definición de un sistema coordenado es arbitraria luego
podemos repetir el análisis pero considerando una posición alternativa de
ejes x−y. En particular nos interesa ubicarlos en el centroide de la superficie
a estudiar como se observa en Fig. 2.10
En este desarrollo, la presión sobre la superficie es correctamente calcu-
lada como F = pc A. Sin embargo, debido al cambio del origen del sistema
coordenado debemos desarrollar nuevamente la posición de acción de la fuer-
za resultante, definiendo el centro de presión ycp . La fuerza resultante debe
generar un momento de igual magnitud que la fuerza distribuida, luego po-
demos escribir

Z
ycp FR = y dF
Z
ycp pc A = y p dA
Z
= y (p0 + ρgξ sin(θ)) dA
Z Z
= p0 y dA + ρg sin(θ) y ξ dA (2.39)

La primera integral de (2.39) es cero debido al origen del sistema coor-

denado. Definiendo la variable auxiliar ξ = ξCG − y podemos reemplazar en
(2.39) obteniendo

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS 23

Figura 2.10: Esquema de una superficie plana sumergida en un fluido. Los

ejes x − y tienen origen en el centroide de la superficie sumergida, siento
el eje x paralelo a la intersección entre el plano de la superficie sumergida
y la superficie libre. En la parte derecha-abajo, se muestra la superficie
sumergida en el plano abatido x − y. Se define una variable auxiliar xi que
mide la distancia desde la superficie libre sobre el plano de la superficie
sumergida.

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

24 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS

 Z Z 
2
ycp pc A = ρg sin(θ) ξCG y dA − y dA
Z
ycp pc A = −ρg sin(θ) y 2 dA
Ixc
ycp = −ρg sin(θ) (2.40)
pCG A

De esta forma, la posición del centro de presión está definida a partir de

la posición del centroide y siempre se encontrará a una profundidad mayor
indicado por el signo menos en (2.40).
Ecuaciones (2.37), (2.38) y (2.40) son todas correctas pero difieren entre
ellas las suposiciones utilizadas como considerar o no la presión en superficie
libre p0 o el origen del sistema coordenado y. El estudiante debe reconocer
estos detalles para poder hacer correcta aplicación de las formulas.

2.5.2. Superficies curvas

En el caso que se deba estudiar las fuerzas sobre una superficie sumergi-
da curva, es posible descomponer este problema en dos subproblemas más
simples. Veamos la Fig. 2.11 donde se observa la superficie curva y la fuerza
ejercida sobre ella producto de la presión. Analizando lo que sucede a una
porción del áreas curva dA podemos identificar una fuerza dF . Esta fuerza es
siempre normal a la superficie luego sobre la superficie curva la dirección de
la fuerza normal cambia. Esta fuerza normal se puede descomponer en una
fuerza horizontal dFH y una componente vertical dFV del siguiente modo

dF = dFH ı̂ + dFV ̂ (2.41)

donde ı̂ y ̂ son los vectores unitarios en la dirección horizontal y vertical

respectivamente. Del mismo modo, la superficie dA se puede descomponer
en sus proyecciones horizontal dAH = dA sin θ y vertical dAV = dA cos θ.
Ahora nuestro problema lo podemos enfrentar en dos pasos.

Fuerzas horizontales

El valor de la fuerza horizontal dFH se puede obtener con lo aprendido

en la sección 2.5.1, es decir

dFH = pc dAv . (2.42)

Generalizando esta idea, la fuerza horizontal FH = dFH será la fuerza

R

actuando sobre la proyección vertical de la superficie curva.

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS 25

Super
cie curva

Figura 2.11: Descomposición de fuerzas y áreas para el calculo de fuerzas

sobre superficies curvas.

Fuerzas verticales
Para el caso de la fuerza vertical, es interesante observar que

dFV = p dAH
dFV = ρgh dAH
dFV = ρg dV
(2.43)

donde p = ρgh es la presión a la profundidad h en que se encuentra la

superficie dA. El producto h dAH = dV esR el volumen de fluido sobre la
superficie dA. Finalmente, la fuerza FV = ρg dV lo que equivale al peso
del fluido que se encuentra sobre la superficie curva.

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

26 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS

2.6. Fuerza de empuje

La fuerza de empuje es una fuerza vertical hacia arriba responsable de la
flotabilidad de un cuerpo sumergido. La fuerza de empuje se opone al peso
del cuerpo. Cuando estas dos fuerzas son iguales se habla de flotabilidad
neutra. Cuando los cuerpos se encuentran parcial o totalmente sumergidos
desplazan un volumen de fluido. El empuje se puede evaluar fácilmente con-
siderando el peso del fluido desplazado

E = ρf Vo g (2.44)

donde ρf es la densidad del fluido, Vo es el volumen sumergido del objeto

y g es la constante de gravedad. Esta ecuación también se le conoce como
Principio de Arquímedes quien en el año 250 AC propuso su trabajo sobre
flotabilidad de los cuerpos, fundando los conceptos de hidrostática.
La fuerza de empuje puede derivarse a partir de los conocimientos de
fuerzas sobre superficies sumergidas. La fuerza de empuje resulta del gra-
diente de presión encontrado entre la superficie superior e inferior del cuer-
po sumergido. Si consideramos un cuerpo sumergido como en la Fig. 2.12.
podemos subdividirlo en pequeños volúmenes dV con forma cilíndrica. En
cada uno de estos cilindros las fuerzas laterales se anularán mientras que
las fuerzas verticales no. Si conocemos la presión en la cara inferior de estos
cilindros y la denominamos p0 , la presión sobre la cara superior a distancia
h será p0 − ρgh. Ambas presiones generarán una fuerza vertical opuesta y
de magnitud diferente. Esta diferencia de fuerzas es el empuje.

Figura 2.12: Diagrama de fuerzas en un cilindro sumergido en un fluido de

densidad ρf

La fuerza inferior esta dada por Finf = p0 A siendo A el área de la base

del cilindro. La fuerza en la cara superior es Fsup = (p0 −ρf gh) A. Realizando
la sumatoria de fuerzas externas sobre la superficie del cilindro obtenemos

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS 27

E = Finf − Fsup
E = p0 A − (p0 − ρf gh) A
E = ρf g h A
E = ρf g Vo (2.45)

donde Vo = h A es el volumen del cilindro. Como se observa en 2.45 la

fuerza de empuje existe mientras el volumen sumergido sea mayor a cero y
la densidad del fluido sea mayor a cero. Esta fuerza actúa en el centro del
volumen sumergido y en dirección vertical hacia arriba.

2.6.1. Estabilidad de un objeto sumergido

Existen al menos dos formas de estabilidad de un cuerpo flotante. La
primera tiene que ver con la magnitud del empuje en relación a la magnitud
del peso. Cuando ambas fuerzas son iguales (siempre opuestas en dirección)
se verifica la flotabilidad neutra. Si el empuje es menor al peso, el objeto
bajará hasta lograr un nuevo equilibrio el cual puede ser llegar al fondo del
contenedor del fluido o aumentar el volumen desplazado aumentando conse-
cuentemente el empuje. Si el empuje es mayor al peso, el objeto subirá hasta
encontrar un nuevo equilibrio lo que se verificará con el objeto parcialmente
sumergido.
El segundo modo de estabilidad de cuerpos flotantes se refiera a la ro-
tación que se puede generar sobre el cuerpo cuando la fuerza de empuje y
el peso no se encuentran alineadas. En este caso se puede generar pares de
fuerzas (momento) restauradoras del equilibrio o pares de fuerzas desestabili-
zadoras. En este segundo caso el cuerpo encontrará una nueva configuración
estable, como es el caso de un iceberg que pierde parte de su masa.

2.7. Distribución de presión en movimiento de cuer-

po rígido
Ciando un fluido, en su conjunto, se mueve de forma acelerada, las leyes
de la hidrostática también son válidas. Este tipo de movimiento se llama mo-
vimiento de cuerpo rígido. Algunos casos típicos de este tipo de movimiento
son:

Un camión que transporta un líquido acelera cuando la luz del semá-

foro es verde.

Un contenedor cilíndrico con un líquido gira sobre su eje de simetría

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

28 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS

Estos dos casos son solo algunas de las posibilidades de movimiento de

cuerpo rígido y las más comunes en general. Existen más caso pero no los
estudiaremos en estos apuntes. El lector puede rpofundizar más sobre este
tema en la Sección 2.9 Pressure Distribution in Rigid-Body Motion de [4].
La ecuación general de hidrostática 2.14 desarrollada en este capítulo
nos servirá para estudiar la distribución de presión en los siguientes casos.
Escribamos 2.14 considerando el vector gravedad ~ g = −g k̂ de la siguiente
forma −∇p + ρ ~ g = ρ~a , luego podemos despejar el gradiente de la presión
como

∇p = ρ (~ a)
g −~ (2.46)
con ~
g= [0, 0, −g]T y~
a = [ax , ay , az ]T .

2.7.1. Aceleración lineal uniforme

En Fig. 2.13 vemos un contenedor de agua bajo los efectos de la gra-
vedad y con aceleración horizontal ax . En este caso la ecuacion genera de
hidrostática es ∇p = ρ (~
g −~ a = [ax , 0, 0]T . Bajo estas condiciones
a) con ~
podemos escribir

∂p
= −ρ ax ,
∂x
∂p
= 0, (2.47)
∂y
∂p
= −ρ g .
∂z

Figura 2.13: Diagrama de un contenedor de agua acelerado uniformemente.

Podemos entender la variación de presión en estas condiciones calculando

la derivada dp(x, y, z)

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS 29

∂p ∂p ∂p
dp = dx + dy + dz . (2.48)
∂x ∂y ∂z
Considerando que no hay variación de presión en la dirección y tenemos

∂p ∂p
dp = dx + dz . (2.49)
∂x ∂z
Una superficie a presión constante (isobárica), como la superficie libre,
cumple dp = 0, luego

∂p ∂p
dp =dx + dz = 0 . (2.50)
∂x ∂z
Si reemplazamos 2.47 obtenemos

−ρ ax dx − ρ g dz = 0 . (2.51)
or
ax
dz = − dx (2.52)
ρ
Integrando podemos calcular la ecuación que describe una superficie iso-
bárica
ax
z=− x + cte (2.53)
ρ
La constante de integración se puede evaluar para la superficie libre con
x = 0 y z = z0 . Le flecha azul en la Fig. 2.14 corresponde a la dirección de
máxima variación de presión.

Figura 2.14: Diagrama de un contenedor de agua acelerado uniformemente.

Para encontrar la distribución de presión podemos integrar (2.49) obte-

niendo

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

30 CAPÍTULO 2. ESTÁTICA DE LOS FLUIDOS

p(x, z) = −ρ ax x − ρ g z + cte (2.54)

La constante la podemos evaluar en un punto conocido, por ejemplo en
p(x = 0, z = z0 ) = p0 luego

p(x = 0, z = z0 ) = −ρ g z0 + cte = p0
cte = p0 + ρ g z0 (2.55)

Finalmente la distribución de presión es

p(x, z) = p0 + ρ g z0 − ρ (ax x + g z) (2.56)

2.7.2. Rotación de cuerpo rígido

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

Capítulo 3

Dinámica de los fluidos

La dinámica de los fluidos es uno de los fundamentos más usados de

la ingeniería aplicada. Desde el agua potable que es un circuito de flujo a
presión en tuberías, las descargas de aguas docimilciliarias que es un flujo
con superficie libre, hasta el motor Merlin que impulsa el primer cohete
espacial reutilizable Falcon9.
Podemos afrontar el estudio de la dinámica de los fluidos desde, al menos,
don formulaciones. La primera modela la naturaleza como una sustancia
continua, donde se plantean ecuaciones que describen el comportamiento de
cada pequeña unidad de volumen, este enfoque se llama análisis diferencial.
Podemos escribir por ejemplo el principio de conservación de la masa de la
siguiente forma

∂ρ ∂ρu ∂ρv ∂ρw

+ + + = 0
∂t ∂x ∂y ∂z
∂ρ ~
+ ∇ · ρV = 0 (3.1)
∂t
donde ρ = ρ(~ x, t) es la densidad en cada punto x ~ = (x, y, z) y tiempo t y
~
V = (u(~ x, t), v(~ x, t)) es el campo de velocidad en cada posición y
x, t), w(~
tiempo.
La segunda formulación describe los fenómenos de la naturaleza mediante
la definición de un volumen de control que rodea el fenómeno de interés,
planteando ecuaciones que definen lo que entra y lo que sale del volumen de
control. Este enfoque se llama análisis integral y es el que utilizaremos en
este capítulo.
Ambos análisis son equivalentes en el sentido que respetan las leyes de la
naturaleza, sin embargo, cada uno de ellos nos permite entregar información
distinta sobre el fenómeno de interés. La descripción de una bomba para agua
se hace regularmente con el análisis integral donde se plantean ecuaciones
que relacionan la presión y el caudal. Sin embargo, un análisis diferencial
del flujo al interior de la bomba es una tarea de extrema dificultad que en

31
32 CAPÍTULO 3. DINÁMICA DE LOS FLUIDOS

años recientes se a facilitado gracias a la mecánica computacional. Por otro

lado, el desarrollo de la aeronáutica no habría sido posible sin el análisis
diferencial y los flujos potenciales.
Un elemento común a ambos análisis es el concepto de campo de veloci-
dad y la descripción de este.

3.1. Campo de velocidad y líneas de corriente

Una de las variables de mayor interés en dinámica de los fluidos es el cam-
po de velocidad V~ = (u, v, w), el cual describe la velocidad de una partícula
que se ubica en la posición x ~ al tiempo t. Asociado al campo de velocidad
tenemos varias líneas que se pueden construir. Entre ellas encontramos las
lineas de corriente las cuales se definen como la linea tangente al vector ve-
locidad en un instante dado. De su definición se entiende que son lineas que
cambian en el tiempo y que dos lineas de corriente no se pueden intersectar
debido a que un punto no puede tener dos velocidades distintas. Sobre una
línea de corriente se verifica que

dx dy dz d~
r
= = = (3.2)
u v w ~
V
donde ~
r = (x, y, z).
Para un instante dado, las líneas de corriente se pueden dibujar conec-
tando el campo de velocidad. Otro elemento interesante de observar es un
tubo de corriente formado por un conjunto cerrado de lineas de corriente
como se observa en la Fig. 3.1. El campo vectorial de velocidad es grafica-
do en algunos puntos con flechas negras. La lineas de corriente en rojo son
tangentes a estas flechas.

Figura 3.1: Representación gráfica de las lineas de corriente sobre el campo

de velocidad y un tubo de corriente formado por un conjunto cerrado de
líineas de corriente.

3.2. Teorema de Transporte de Reynolds

Antes de desarrollar el Teorema de Transporte de Reynolds, debemos
dejar claro que es un sistema desde el punto de vista de la mecánica de los

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 3. DINÁMICA DE LOS FLUIDOS 33

fluidos, a lo que llamaremos un Volumen de Control (VC). Un Volumen de

Control es un concepto que define una región del espacio que está limitada
por una Superficie de Control (SC) y sobre el cual nosotros decidimos poner
nuestra atención. Si el problema a estudiar es tridimensional el Volumen
de Control será un volumen y la Superficie de Control será un área. Por el
contrario, si el problema de estudio es bidimensional, el VC será un área y
la SC será una curva cerrada en torno a esta área. Un volumen de control
puede permanecer fijo, puede moverse y puede deformarse en el tiempo.
Además, a través de la superficie de control puede fluir masa. En este curso
consideramos que el volumen de control es fijo e indeformable. Este caso es
el más simple y aún válido en muchas condiciones típicas de ingeniería.
El Teorema de Transporte de Reynolds nos da una herramienta matemá-
tica que nos permite formular rigurosamente como cambia una propiedades
extensiva (masa, momento lineal, energía por ejemplo) al interior del Volu-
men de Control.

dBsys d
Z Z  
= ρb dV + ~ · n̂ dA
ρb V (3.3)
dt dt VC SC

Ecuación (3.3) es válida para un Volumen de Control fijo donde el tér-

mino de la izquierda nos indica el cambio temporal de la propiedad B. El
primer término de la derecha nos indica razón de acumulación de B en el
volumen de control y el ultimo término nos indica el flujo neto de B que
entra o sales del VC. La variable b = B/m es una propiedad intensiva que
se deriva de B al dividirla por la masa.
Para flujos estacionarios tenemos

dBsys
Z  
= ~ · n̂ dA
ρb V (3.4)
dt SC

3.3. Conservación de la masa

La conservación de la masa es un principio fundamental de la física y
nos permite describir con precisión aplicaciones de ingeniería de nuestra
vida cotidiana. Caso especial es cuando se aplica la formula propuesta por
Einstein que relaciona directamente masa y energía E = mc2 .
La conservación de la masa aplicada a los fluidos nos permite describir
el comportamiento de flujos de fluidos compresibles e incompresibles. Por
ejemplo si combinamos 2 Kg de Hidrógeno H2 y 16 Kg de Oxigeno O2
obtenemos 18 kg H2 O.
Si se define un volumen de control indeformable, la masa que entra al
volumen menos la masa que sale del volumen tiene que ser igual al cambio
de masa al interior del volumen de control, todo esto en un intervalo de
tiempo definido. En ecuaciones esto se escribe

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

34 CAPÍTULO 3. DINÁMICA DE LOS FLUIDOS

Figura 3.2: Representación gráfica de la ley de conservación de la masa

aplicada a un volumen de control. En rojo se identifica el volumen de control.

d mvc
ṁen − ṁsa = , (3.5)
dt

donde ṁen es la masa que entra al volumen de control por unidad de tiempo,
ṁsa es la masa que sale del volumen de control por unidad de tiempo y
d mvc
dt es la variación de masa por unidad de tiempo al interior del volumen
de control, lo que representamos gráficamente en fig. 3.2
Si consideramos que la entrada al volumen de control en fig. 3.2 tiene
una velocidad V1 y área A1 la masa por unidad de tiempo se expresa como
ṁen = ρ1 V1 A1 . Siguiendo esta representación, para una salida del volumen
de control con velocidad V2 y área A2 la masa por unidad de tiempo se
expresa como ṁsa = ρ2 V2 A2 . Si asumimos un flujo estacionario, es decir, la
cantidad de masa al interior del volumen de control no cambia en el tiempo
podemos expresar la conservación de la masa como

ρ1 V1 A1 = ρ2 V2 A2 . (3.6)

Ahora si limitamos el análisis a un fluido incompresible (3.6) se reduce

a Q = V1 A1 = V2 A2 .

Ejercicio: Analice el flujo en una tubería de sección constante, para un

fluido de densidad constante, entre las secciones x1 y x2 con una inclinación
sobre la linea horizontal arbitraria θ. Utilizando el principio de conservación
de la masa, desarrolle las ecuaciones que describen la velocidad en la sección
x2 en función de la velocidad en la sección x1 .

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 3. DINÁMICA DE LOS FLUIDOS 35

3.4. Conservación del momento lineal

El momento lineal se define como p~ = mV~ el producto de la masa y su
velocidad, por lo tanto es una propiedad vectorial del flujo.
Recordando la segunda ley de Newton,

~
F = m~a
dV~
~
F = m
dt
dmV ~
~
F =
dt
~ d~
p
F = (3.7)
dt
donde hemos asumido que la masa se mantiene constante. En este caso, la
Segunda Ley de Newton nos indica que un cambio en el momento lineal
tiene asociado una fuerza neta actuando sobre el cuerpo.

3.5. Conservación de la energía

Si analizamos lo que sucede sobre una línea de corriente en un flujo esta-
cionario, es decir, un flujo que no cambia en el tiempo, podemos definir un
volumen de control que considera un volumen alrededor de esta línea de co-
rriente como se muestra en la Fig. 3.3. En ese volumen de control definimos
una dirección longitudinal a la línea de corriente s y una dirección normal
a la linea de corriente n. El volumen de control está inclinado de acuerdo a
la línea de corriente un angulo θ con respecto a una línea horizontal. Consi-
derando los efectos de la gravedad, sobre este volumen de control podemos
aplicar el balance de fuerzas por unidad de volumen (2.14) desarrollado en
la sección 2.2.

−∇p − ρ g k̂ = ρ~
a, (3.8)
Esta ecuación es vectorial y se puede aplicar en cualquier dirección de
acuerdo a nuestra definición de sistema coordenado. En particular aplicán-
dola en la dirección s podemos escribir

∂p ∂V
− − ρg sin θ = ρV (3.9)
∂s ∂s
donde ρ es la densidad del fluido y V es la magnitud de la velocidad en la
dirección s. Integrando sobre la línea de corriente y considerando que ρ y g
son constantes obtenemos
1 2
ρV + p + ρgz = C (3.10)
2

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

36 CAPÍTULO 3. DINÁMICA DE LOS FLUIDOS

Figura 3.3: Volumen de control definido sobre una línea de corriente.

siendo C una constante sobre la línea de corriente. Esta ecuación es una

de las más famosas en Mecánica de los Fluidos y fue propuesta por Daniel
Bernoulli en su trabajo Hydrodynamica en 1738 por lo que se le conoce
como el Principio de Bernoulli. Esta es una de las ecuaciones más abusadas
de la ingeniería luego es vital recopilar las suposiciones utilizadas para su
derivación

Sobre una línea de corriente

Flujo estacionario

Densidad constante

No considera efectos viscosos

3.6. Ecuación de Bernoulli

1 2
ρV + p + ρgz = C (3.11)
2
La ecuación de Bernoulli es una de las herramientas mas importantes en
la mecánica de fluidos para entender el comportamiento de la dinámica de
los fluidos. Esta ecuación describe el balance de energía que tiene el fluido
sobre una línea de corriente relacionando los términos de presión p, velocidad
V y cota geométrica z; para un fluido de densidad constante ρ y sobre un
campo gravitacional constante g. El valor de la constante C depende de la
condición del flujo sobre la línea de corriente donde se evalúe (3.11).
Para utilizar (3.11) correctamente es importante recordar siempre co-
mo se derivó a partir de los principios fundamentales. Las condiciones de
aplicabilidad son:

Flujo estacionario: Esta condición es muy restrictiva y nos indica que

(3.11) no puede aplicarse en los inicios y términos de procesos de trans-
porte de fluidos donde aun los transientes dominan el problema.

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 3. DINÁMICA DE LOS FLUIDOS 37

Flujo sobre una línea de corriente: Le ecuación de Bernoulli se derivó

realizando un balance de fuerzas sobre un volumen de control a lo largo
de una línea de corriente. Esta condición implica que la constante C de
la ecuación (3.11) depende de la línea de corriente donde se aplique la
ecuación. Dicho de otro modo, en general C será distinto sobre distintas
líneas de corriente.

Efectos viscosos despreciables: La ecuación de Bernoulli no considera

los efectos viscosos. Cada flujo experimenta en alguna cantidad los
efectos de la fricción debido a la interacción de fluidos viscoso y las
paredes que están en contacto con este. Los efectos viscosos pueden o
no ser despreciables.

Elementos que extraen o aportan energía al flujo: La ecuación de Ber-

noulli se deriva a partir de un balance de fuerzas sobre una línea de
corriente para flujo estacionario. La presencia de bombas, turbinas,
ventiladores o cualquier accesorio que perturbe el flujo estacionario
inhabilita el uso de la ecuación de Bernoulli. Estos elementos se con-
siderarán de manera específica en el balance de energía en el siguiente
capítulo.

3.6.1. Formas de la ecuación de Bernoulli

La ecuación de Bernoulli puede escribirse de diferentes formas al multi-
plicarla por diferentes factores. Esta representaciones están caracterizadas
por las unidades en las cuales se expresa. Si estudiamos las unidades de la
forma descrita en (3.11)

Forma 1:
La primera forma que veremos de esta ecuación es la que derivamos en
3.5 donde

1 2
ρV + p + ρ g z = C1 , (3.12)
2
Analizando las dimensiones de esta forma y considerando el SI de uni-
dades podemos ver que

1 2 kg m 2 kg m 1 N
   
ρV = = 2 = 2 = Pa, (3.13)
2 m 3 s s m 2 m
N
[p] = = Pa,
m2
kg m kg m 1 N
[ρ g z] = 3 2
m= 2 2
= 2 = Pa.
m s s m m

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

38 CAPÍTULO 3. DINÁMICA DE LOS FLUIDOS

Figura 3.4: Ilustración que presenta zonas en las cuales la ecuación de Ber-
noulli se puede aplicar y donde no se puede aplicar. Esta imagen es del libro
Fluid Mechanics [4]

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 3. DINÁMICA DE LOS FLUIDOS 39

De esta forma verificamos la consistencia dimensional de la ecuación y

encontramos que las unidades de esta forma son P a es decir dimesniones de
presión.

Forma 2:
La segunda forma de esta ecuación la obtenemos al dividir 3.12 por el
factor constante ρg obteniendo

V2 p C1
+ + z = C2 = , (3.14)
2g ρg ρg
Analizando las dimensiones de esta forma y considerando el SI de uni-
dades podemos ver que

m 2 m2
" #
V2


s s2
= m = m = m, (3.15)
2g s2 s2
N kg m
p Pa
 
m2 s2 m2
= kg m
= kg m
= kg m
= m,
ρg
m3 s2 m3 s2 m3 s2
[z] = m ,

es decir, la energía aportada por cada término está en dimensiones de lon-

gitud y considerando el SI de unidades, unidades de metros o altura de una
columna de fluido.

Forma 3:
La tercera forma de esta ecuación la obtenemos al multiplicar 3.14 por
el factor constante g obteniendo

1 2 p C1
V + + z g = C3 = C2 g = , (3.16)
2 ρ ρ

2
1 2 m m2
  
V = = , (3.17)
2 s s2
N kg m
p Pa m2
 
m2 s2 m2
= kg
= kg
= kg
= ,
ρ s2
m3 m3 m3
m m2
[z g] = m = .
s2 s2
concluyendo que la energía aportada por cada témino está en dimensiones
de energía dividido por masa lo que se representa en el SI de unidades por
m2 /ss = J/kg. Recordar que J = N m.

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

40 CAPÍTULO 3. DINÁMICA DE LOS FLUIDOS

Es importante realizar que la constante de (3.12), (3.14) y (3.16) tienen

valor distinto. Esto nos dice que las tres expresiones son equivalentes en su
capacidad de representar las condiciones del flujo pero no deben intercam-
biarse sin extrema atención a las condiciones particulares de su formulación.

3.6.2. Representación gráfica de la ecuación de Bernoulli

Si consideramos la segunda forma de la ecuación de Bernoulli 3.14, la
energía disponible en el flujo esta en dimensión de longitud de acuerdo a los
siguientes términos:

V 2 /(2g): energía cinética.

p/(ρ g): energía mecánica o de flujo.

z: energía potencial.

Gracias a las unidades de esta forma, los tres términos pueden agruparse
en forma gráfica de la siguiente forma. Imaginemos una línea de corriente
que aumenta su cota geométrica en forma constante como se ve en la figura
3.5. Sobre una línea de corriente que va desde x1 a x2 la energía total dispo-
nible es representada por la constante C. La energía potencial varía con la
cota geométrica desde z1 a z2 . Si suponemos que la velocidad se mantiene
constante a lo largo de la línea de corriente, la energía cinética representada
por el término V12 /(2g) = V22 /(2g) se mantendrá constante. En consecuen-
cia, para cumplir la conservación de la energía descrita por la ecuación de
Bernoulli la energía de flujo cumplirá que p1 /(ρ g) > p2 /(ρ g). Esta forma
gráfica de representar la condición de flujo es muy útil para entender las
relaciones entre las tres formas de energía. Esta forma gráfica es equivalente
a evaluar la ecuación de Bernoulli en x1 y x2

p1 v2 p2 v2
+ 1 + z1 = + 2 + z2 . (3.18)
ρg 2g ρg 2g

Si agrupamos los términos p/(ρ g) + z se define la Línea de Gradien-

te Hidráulica (LGH). Si a la LGH le agregamos el término de velocidad
p/(ρ g) + V 2 /(2g) + z se define la Línea de Energía (LE). Ver por ejemplo
esta secuencia de elementos.
Adicionalmente a esta representación gráfica, siguiendo la forma 1 de la
ecuación, tenemos los términos p que llamamos presión estática, el término
1 2 1 2
2 ρV que llamamos presión dinámica y el término p + 2 ρV presión total.
En estricto rigor, los términos presión estática y presión son lo mismo y
dependen del estado en el cual se encuentre el fluido y no de su velocidad.

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 3. DINÁMICA DE LOS FLUIDOS 41

Figura 3.5: Representación gráfica de la energía disponible descrita por la

ecuación de Bernoulli.

3.6.3. Aplicaciones de la ecuación de Bernoulli

Presión de estancamiento
Para entender la presión de estancamiento (PS) primero debemos enten-
der que es un punto de estancamiento. Si recordamos el concepto de lineas de
corriente presentado en sección 3.1 y la ecuación de Bernoulli, podemos ver
que la energía que tiene el fluido sobre la linea de corriente debe ser cons-
tante y puede distribuirse “idealmente” en los 3 términos de la ecuación.
Estos son presión, velocidad y altura. Si hipotizamos una linea de corriente
que reduce su velocidad hasta llegar a cero y sin cambiar su altura, como se
presenta en el punto 2 en la Fig. 3.6, veremos que la energía será transferida
al término de presión como se observa a continuación.
En el punto donde se produce el estancamiento, es decir, donde la veloci-
dad se reduce a cero, la energía del fluido es transferido en forma de presión.
Sobre esa línea de corriente podemos plantear la siguiente ecuación

p1 v2 p2 v2
+ 1 + z1 = + 2 + z2 . (3.19)
ρg 2g ρg 2g
Asumiendo que la cota geométrica desde 1 a 2 no cambia y utilizando el
hecho que la velocidad en el punto de estancamiento es cero tenemos

1
p2 = p1 + ρv12 . (3.20)
2

Tubo Pitot
El tubo pitot es un instrumento utilizado para medir la velocidad de
un fluido y puede ser observado en aviones y autos de formula 1. Este ins-

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

42 CAPÍTULO 3. DINÁMICA DE LOS FLUIDOS

1 2
PS

Figura 3.6: Representación gráfica de la línea de corriente y el punto de

estancamiento.

0
2

4
3

Figura 3.7: Representación de un tubo de Pitot.

trumento fue inventado por el ingeniero francés Henri Pitot en 1732. Este
instrumento hace uso del concepto de punto de estancamiento al enfrentar
un tubo a la dirección del fluido. Cuando este fluido es forzado a detenerse
al interior del tubo la energía cinética del flujo es transformada en presión
midiéndose la presión de estancamiento. Comparando la presión de estan-
camiento con la presión estática podemos derivar una expresión para la
velocidad del fluido con respecto al instrumento. De acuerdo a la figura 3.7
el fluido tiene una presión estática p0 = p1 . Siguiendo el razonamiento de la
presión de estancamiento, la presión p2 = p0 + 21 ρv02 . El tubo pitot registra
estas dos presiones en p4 = p1 y p3 = p2 . Con esta información podemos
concluir que comparando la presión de estancamiento con la presión estática
tenemos

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 3. DINÁMICA DE LOS FLUIDOS 43

s
2(p3 − p4 )
v0 = (3.21)
ρ

Tubo Venturi

El físico italiano Giovanni Venturi (1746-1822) estudió el flujo en tu-

berías con contracciones y expansiones cónicas. Estudiando esta geometría
Venturi se dio cuenta como la energía que tiene el flujo puede variar entre
sus términos presión, velocidad y cota geométrica siguiendo el balance de
energía propuesto en la ecuación de Bernoulli 3.11. Para entender bien es-
te proceso es necesario primero entender el principio de conservación de la
masa cuando el flujo recorre esta geometría convergente divergente. Obser-
vando la geometría en la Fig. 3.8 podemos analizar el volumen de control al
interior de la tubería desde la sección 1 y 2.
Estudiemos el caso en que el fluido que recorre el tubo venturi es in-
compresible y asumamos un régimen estacionario. Bajo estos supuestos el
balance de masa aplicado al volumen de control (linea segmentada roja en
Fig. 3.8) nos indica

ρ, V1 , A1 = ρ, V2 , A2 (3.22)
V1 , A1 = V2 , A2 (3.23)
A1
V2 = V1 (3.24)
A2

De acuerdo a la geometría el área de la sección 1 y 2 son A1 = πD 2 /4

y A2 = πd2 /4 respectivamente. Si asumimos que D > d entonces A1 > A2
por lo tanto

V2 > V1 . (3.25)

De acuerdo a este análisis, el fluido debe aumentar su velocidad cuando

recorre la sección de menor diámetro de la tubería para cumplir con la ley
de conservación de la masa. Este resultado es muy importante de entender
ya que nos da argumentos para el siguiente análisis aplicando la ley de
conservación de la energía mediante la ecuación de Bernoulli. Asumamos que
existe una linea de corriente que recorre la tubería convergente-divergente a
lo largo de su eje axial y la energía disponible sobre esta línea de corriente
es representativa de la energía del flujo. Planteando el balance de energía
entre las secciones 1 y 2 podemos escribir,

p1 V2 p2 V2
+ 1 + z1 = + 2 + z2 . (3.26)
ρg 2g ρg 2g

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

44 CAPÍTULO 3. DINÁMICA DE LOS FLUIDOS

Figura 3.8: Representación de un tubo de Venturi.

Si asumimos que la energía disponible en los términos de presión y ve-

locidad son mucho mayores a los términos geométricos podemos despreciar
z1 y z2 de la ecuación anterior y escribir,

p1 V2 p2 V2
+ 1 = + 2 . (3.27)
ρg 2g ρg 2g
Ahora podemos despejar la presión en la sección 2 obteniendo
1  
p2 = p1 + ρ V12 − V22 (3.28)
2
De acuerdo al análisis anterior 3.25 con respecto a la conservación de la
masa podemos concluir que el término V12 − V22 < 0 y por lo tanto en la

sección 2 la presión p2 < p1 será menor a la presión en la sección 1. Los

valor numéricos de esta diferencia de presión dependerán de las condiciones
particulares del problema. La conclusión general que podemos obtener es
que siempre en una restricción del área observaremos un aumento de la
velocidad y una disminución de la presión. Esta condición del tubo venturi
es ampliamente utilizada en la industria para introducir mediante la succión
un fluido a la red principal de distribución. Otra aplicación del tubo venturi
es la de medir la velocidad del flujo aunque no es la más recomendada de
acuerdo a criterios de eficiencia que se podrán observar en la experiencia de
laboratorio.

Fenómeno de cavitación
De acuerdo al análisis presentado en la sección anterior, nos dimos cuenta
que cuando el flujo se ve enfrentado a una restricción del área de la sección,
se producirá un aumento de la velocidad y en consecuencia una disminución
de la presión. Esta disminución de la presión puede dar origen al fenómeno
llamado cavitación. La cavitación se produce cuando la presión de un líquido
se reduce bajo la presión de vapor en ese fluido, formado pequeñas burbujas
de la misma sustancia pero en fase gaseosa. Todo esto sucede en lugar muy
acotados de la geometría del problema y en general a consecuencia de vál-
vulas mal cerradas. El mayor problema se produce cuando esta burbujas en

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 3. DINÁMICA DE LOS FLUIDOS 45

fase gaseosa avanzan con el flujo y aumentan su presión forzándolas a im-

plosionar. En ese proceso se liberan ondas de presión que recorren el fluido
y golpean las paredes de las tuberías generando procesos de erosión.

Teorema de Toricelli

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

46 CAPÍTULO 3. DINÁMICA DE LOS FLUIDOS

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

Capítulo 4

Flujos en tuberías

Analizar el flujo en tuberías es una de las áreas más importantes de

las aplicaciones de mecánica de los fluidos. En el capítulo anterior vimos
la derivación de la ecuación de Bernoulli que representa la conservación de
energía sobre una línea de corriente. Una gran limitante a la aplicabilidad de
la ecuación de Bernoulli es que no considera efectos viscosos y por ende no
se hace cargo de las pérdidas por fricción en los flujos en tuberías. Durante
décadas ingenieros y científicos estudiaron el problema y construyeron las
herramientas que hoy tenemos para el estudio y diseño de flujos en tuberías.
En 1839 en Ingeniero Alemán Gotthilf Heinrich Ludwig Hagen propuso
la existencia de dos régimen de flujos viscosos. Midiendo la presión del flujo
al inicio y fin de una tubería de bronce encontró que la pérdida de presión
a lo largo de la tubería se podía estimar como ∆p ∝ LQ R4
, es decir propor-
cional al caudal y largo e inversamente proporcional al radio. Aumentando
el caudal encontró que la pérdida de carga era proporcional a ∆p ∝ Q2 . En
ese momento no realizó el efecto de la viscosidad del fluido pero si encontró
evidencia de dos régimen, los que hoy llamamos laminar y turbulento.
En 1883 el Ingeniero Británico Osborn Reynolds mostró que el cambio
de régimen depende de la relación entre las fuerzas inerciales y las fuerzas
viscosas, proponiendo el parámetro adimencional llamado número de Rey-
nolds

ρvd
Re = (4.1)
µ
donde ρ y µ son la densidad y viscosidad del fluido, v es la velocidad media
en la tubería y d el diámetro de la tubería. Cuando las fuerzas inerciales
tienen mayor influencia que las fuerzas viscosas hablamos de un número de
Reynolds alto y flujo turbulento. Por el contrario, cuando las fuerzas visco-
sas tienen más influencia que las fuerzas inerciales hablamos de un número
de Reynold bajo y flujo laminar. Para efectos de este curso, hablaremos de
Re = 2300 como transición entre ambos régimen luego para Re < 2300 tene-
mos flujo laminar, para Re > 2300 tenemos flujo turbulento. Es importante

47
48 CAPÍTULO 4. FLUJOS EN TUBERÍAS

Figura 4.1: Balance de fuerzas en un flujo en tubería.

reconocer que es el flujo quien es laminar o turbulento y no el fluido. Un

fluido altamente viscoso puede fluir en régimen laminar y turbulento según
la velocidad a la cual se desplace.
Ejercicio: Analizar las unidades de cada término de la ecuación del
número de Reynolds y verificar que es adimensional.
Ejercicio: Calcular la velocidad para la cual se verifica la transición
entre laminar y turbulento para aire y agua en una tubería de 5 cm de
diámetro.

4.1. Pérdidas de carga debido a la fricción

Las pérdidas de presión encontradas por G. H. Ludwig Hagen y carac-
terizadas entre régimen laminar y turbulento por el número de Reynolds
son debidas a los efectos viscosos entre el fluido y las paredes rugosas de
la tubería. Este efecto puede ser modelado mediante un esfuerzo de corte τ
que se opone al flujo de la siguiente forma
Realizando un balance de fuerzas sobre el volumen de control presentado
en la fig. 4.1 vemos como el peso del fluido debido a la gravedad actúa en
el centro de masa, a la entrada y salida la presión por el área aportan con
fuerzas en la dirección longitudinal de la tubería y finalmente el esfuerzo
de corte multiplicado por el área de acción, es decir el manto interior de la
tubería se opone al flujo. Esto se puede escribir de la siguiente forma
X
F = p1 A1 − p2 A2 + γAL sin(θ) − τ LPer = 0 (4.2)

Asumiendo una tubería de sección constante, A1 = A2 = A. El manto

interior de la tubería es el largo multiplicado por el perímetro Per = 2πr y
el área es A = πr 2 donde r es el radio de la tubería. Utilizando la relación
geométrica L sin(θ) = z1 − z2 podemos escribir

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 4. FLUJOS EN TUBERÍAS 49

∆p 2τ L
+ ∆z = (4.3)
ρg ρgr
De acuerdo a los innumerables experimentos realizados en laboratorio se
2
ha propuesto que el esfuerzo de corte τ = λ ρv8 donde λ define el coeficiente
adimensional de fricción. Reemplazando esto en (4.3) obtenemos luego de
algunos pasos algebraicos

∆p Lv 2
+ ∆z = λ (4.4)
ρg d2g
donde d = 2r es el diámetro de la tubería. En 1850 Julius Weisbach propuso
la ecuación que conocemos como Darcy-Weisbach

Lv 2
hf = λ (4.5)
d2g
donde hf son las perdidas de carga debido a los efectos viscosos.

4.1.1. Flujos laminares

Para flujos laminares Darcy propuso la relación

64
λ= (4.6)
Re
lo que nos permite calcula las pérdidas por fricción

Lv 2
hf = λ (4.7)
d2g
64 Lv 2
= (4.8)
Re d2g
32 µ L v
= (4.9)
ρ g d2
128 µ L Q
= (4.10)
π ρ g d4

4.1.2. Flujos turbulentos

Para flujos turbulentos Prandtl en 1935 propuso que el coeficiente de
fricción para tuberías de paredes suaves podía ser evaluado como

1  √ 
√ = 2 log Re λ − 0,8 (4.11)
λ
Para el caso de tuberías con paredes rugosas se propuso la siguiente
relación

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

50 CAPÍTULO 4. FLUJOS EN TUBERÍAS

1 ǫ/d
 
√ = −2 log (4.12)
λ 3,7
donde ǫ es la rugosidad de las paredes de la tubería con d el diámetro de la
tubería. El término ǫ/d es la rugosidad relativa.

4.1.3. Ecuación de Colebrook

En 1939 Colebrook propuso una ecuación unificada para calcular el co-
eficiente de fricción en flujos laminares y turbulentos.

1 ǫ/d 2,51
 
√ = −2 log + √ (4.13)
λ 3,7 Re λ
En 1944 Moody graficó esa ecuación para que sea más fácil de utilizar en
los cálculos. Este gráfico es conocido como el diagrama de Moody y puede
ser consultado en el siguiente enlace donde se observa también la rugosidad
de diferentes materiales.
La ecuación de Colebrook es una función de tres parámetros del problema
los cuales son el coeficiente de fricción, la rugosidad relativa y el número de
Reynolds. Esto lo podemos escribir como λ = λ(ǫ/d, Re).
Con esta información podemos nuevamente escribir la ecuación de ener-
gía pero ahora considerando las pérdidas por fricción siendo esta una expre-
sión más realista del fenómeno de interés en comparación a la ecuación de
Bernoulli (3.11).

p1 V2 p2 V2
+ 1 + z1 = + 2 + z2 + hf (4.14)
ρg 2g ρg 2g

4.2. Típicos problemas en flujos en tuberías

Cuando se estudia el flujo en una tubería, se deben plantear diferentes
ecuaciones para describir el problema. En particular, las ecuaciones funda-
mentales son el balance de energía, balance de masa, número de Reynolds
y la ecuación de Colebrook. Según la pregunta que deseamos responder se
reconocen tres estrategias de cálculo

4.2.1. Cálculo de las pérdidas de carga

pdf

4.2.2. Cálculo del caudal

pdf

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 4. FLUJOS EN TUBERÍAS 51

4.2.3. Cálculo del diámetro

pdf

4.2.4. Cálculo mediante herramientas numéricas

Una alternativa moderna a las estrategias propuestas es la implementa-
ción de un algoritmo computacional que resuelva las ecuaciones pertinentes.

4.3. Pérdidas de carga por singularidades

Para poder implementar una red de tuberías será necesario utilizar ele-
mentos adicionales a las tuberías, como por ejemplo, codos, te, válvulas.
Todos estos elementos generan pérdidas de energía en el flujo en ubicaciones
precisas las que llamaremos singularidades. Estas pérdidas están asociadas a
los cambios bruscos de dirección y recirculaciones locales las cuales no se han
estudiado profundamente, de hecho lo que se usa actualmente es mediciones
experimentales. Normalmente el fabricante de uno de estos elementos decla-
ra la relación entre el cambio de presión ∆p/ρg y la carga cinética v 2 /2g
cuando el flujo recorre el elemento (la singularidad).
∆p
ρg
K= V2
. (4.15)
2g

Escribiendo
∆p V2
= hm = K , (4.16)
ρg 2g
podemos cuantificar las pérdidas asociadas a estos elementos lo que nos
permitirá mejorar la precisión del balance de energía

p1 V2 p2 V2
+ 1 + z1 = + 2 + z2 + hf +
X
hm (4.17)
ρg 2g ρg 2g
El término hm significa la suma de todos los elementos. Especial aten-
P

ción se debe tener cuando hay cambios de diámetro en las tuberías. Los
valor de K pueden ser consultados en los libros clásicos de mecánica de flui-
dos [3, 4] y directamente por los fabricantes de cada elemento. Para efectos
de este curso utilizaremos los valores presentados en [4] y que pueden ser
consultados en el siguiente enlace.

4.4. Sistemas de tuberías

En general el flujo en tuberías considera más de una tubería luego po-
demos clasificar los casos principales como tuberías en seria y tuberías en
paralelo

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

52 CAPÍTULO 4. FLUJOS EN TUBERÍAS

A 3 B
2
1

Figura 4.2: Sistema de tuberías en serie. El flujo entra en A y sale en B. Las

tuberías 1, 2 y 3 tienen áreas A1 , A2 y A3 respectivamente.

4.4.1. Sistemas de tuberías en serie

Cuando nos enfrentamos a un flujo en tuberías en serie tenemos que
identificar los elementos fundamentales que nos permitan plantear las ecua-
ciones pertinentes. Imaginemos tres tuberías de distinta sección conectadas
en serie como se observa en la Fig. 4.2.
El balance de masa nos indica que si la tubería es indeformable, el flujo
másico que entra ṁA = ṁB . Si consideramos un fluido incompresible tene-
mos que el caudal que entra por A debe ser igual al caudal que sale por B,
es decir QA = QB = Q, luego podemos plantear

v1 A1 = v2 A2 = v3 A3 = Q (4.18)
Ahora si aplicamos el balance de energía podemos concluir que las pér-
didas entre A y B serán la suma de las pérdidas de carga en la tuberia 1
más la 2 más la 3.

hf = hf 1 + hf 2 + hf 3 (4.19)

L1 v12 L2 v22 L3 v32

hf = λ1 + λ2 + λ3 (4.20)
d1 2 g d2 2 g d3 2 g

4.4.2. Sistemas de tuberías en paralelo

Cuando nos encontramos con una geometría donde el flujo puede dividir-
se en dos o mas tuberías es necesario poner especial atención en reconocer los
principios de conservación. El primer elemento fundamental es plantear de
buena forma el balance de masa. Si observamos la Fig. 4.3 y considerando el
flujo que entra en A y sale en B debemos considerar la suma de los caudales
que viaja por cada tubería. Cada punto de bifurcación lo definiremos como
un nodo donde planteamos el balance de masa, es decir en Fig. 4.3 decimos

QA = Q1 + Q2 + Q3 = QB (4.21)

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 4. FLUJOS EN TUBERÍAS 53

1
A B

Figura 4.3: Sistema de tuberías en paralelo. El flujo entra en A y sale en B.

Las tuberías 1, 2 y 3 tienen áreas A1 , A2 y A3 respectivamente.

Desde el punto de vista del balance de energía, las pérdidas de carga

totales que se encuentren en cada tubería 1, 2 y 3 deberán ser iguales. Esta
condición se puede derivar al plantear el balance de energía en cada tubería.
Por simplicidad asumiremos que solo hay pérdidas por fricción lo que no
quita generalidad a las ideas propuestas. Evaluemos una linea de corriente
que inicia en A, termina en B y recorre la tubería 1.

pA vA 2 pB v2 L1 v12
+ + zA = + B + zB + λ1 (4.22)
ρg 2g ρg 2g d1 2g
Del mismo modo para las tuberías 2 y 3 tenemos

pA vA 2 pB v2 L2 v22
+ + zA = + B + zB + λ2 (4.23)
ρg 2g ρg 2g d2 2g
pA vA 2 pB 2
vB L3 v32
+ + zA = + + zB + λ3 (4.24)
ρg 2g ρg 2g d3 2g
Esto nos permite concluir

L1 v12 L2 v22 L3 v32

λ1 = λ2 = λ3 (4.25)
d1 2g d2 2g d3 2g
Si esto no se cumpliera, la presión en B sería distinta si se recorren
distintos caminos, lo que contradice la suposición de una condición de flujo
en régimen estacionario.
Las velocidades de cada sección dependerán de los diámetros de las tu-
berías y las pérdidas que encuentre el flujo.

4.4.3. Estanques conectados

Esta tipología de problemas se refiere a los casos donde hay poca clari-
dad cual es la condición de equilibrio. Por ejemplo 3 estanques conectados

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

54 CAPÍTULO 4. FLUJOS EN TUBERÍAS

Qin,A Qou,B
A 1 B

Circuito 1 2

3 C
Qou,C

Figura 4.4: Sistema de tuberías en redes cerradas.

mutuamente por tuberías. Ver por ejemplo [3] [4]

4.4.4. Redes cerradas

Una combinación de tuberías en serie y en paralelo puede ser un modelo
más complejo de analizar. Un ejemplo de esto es la metodología propuesta
por el Ingeniero Estructural Hardy Cross en 1936 que permite aproximar la
solución a una red como la que se muestra en la Figura 4.4.
Para resolver este tipo de redes el método Hardy Cross asume las si-
guientes condiciones: La energía del flujo atribuida al término de velocidad
y las pérdidas por singularidades son despreciables frente a las pérdidas por
fricción y presión del sistema. Adicionalmente se asume que conocemos los
caudales de entrada y salida de la red, los diámetros y longitudes de las
tuberías y que las pérdidas por fricción responde a una ley de acuerdo a
hf = rQn con r y n números reales.
Bajo estos supuestos que definen una aproximación a la realidad, el mé-
todo si respeta los prinicpios de conservación de la masa donde los caudales
de entrada y salida en cada bifurcación son iguales; y el principio que indi-
ca que en tuberías en paralelo las pérdidas de carga son iguales lo que se
traduce que en un circuito cerrado las pérdidas son cero.
Bajo estas condiciones Hardy Cross propone dos métodos:

Método de balance de carga: Se impone continuidad de masa y se

corrige evaluando la continuidad de las pérdidas de carga.

Métodos de balance de caudal: Se impone continuidad en las pérdidas

de carga y se corrige evaluando la continuidad de masa.

Si revisamos la suposición del cálculo de las pérdidas de carga hf = rQn ,

podemos confrontarlo con dos modelos clásicos en la literatura.
El primero es la ley de Darcy-Weisbach que plantea

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 4. FLUJOS EN TUBERÍAS 55

LV 2
hf = λ (4.26)
d2g
LQ2 16
hf = λ (4.27)
d2gπ 2 d4
8λL 2
hf = Q (4.28)
gπ 2 d5
8λL
donde n = 2 y r = gπ 2 d5 .
Una segunda alternativa es la ley de Hazen-Williams que propone

L10,67
hf = Q1,85 (4.29)
C 1,85 d4,87
L10,67
donde n = 1,85 y r = C 1,85 d4,87 donde C es un coeficiente de fricción que
depende del material de las tuberías.
La teoría que sustenta este métodos se basa en una aproximación cons-
truida siguiendo las Series de Taylor del siguiente modo. En la Fig. 4.4 vemos
un circuito cerrado con tres nodos donde se bifurca el caudal. Asumimos que
conocemos los caudales de entrada Qin,A y de salida Qou,B Qou,C . Iniciamos
proponiendo un caudal para cada tubería del sistema Qi,0 . Si recorremos el
circuito en la dirección de la flecha roja podemos decir que

ri Qni i = r1 Qn1,0 + r2 Qn2,0 + r3 Qn3,0

X
1 2 3
(4.30)
i

Si los caudales propuestos inicialmente para Qi están correctos, i ri Qni i =

P

0 esta sumatoria de pérdidas por fricción debe ser 0 como lo plantea la teo-
ría. Por el contrario, si los Qi no están correctos, la sumatoria de pérdidas
en el circuito cerrado no será cero. En este caso, se propone una corrección
dado por ∆Q para las tres tuberías obteniéndose los caudales nuevos en
el proceso iterativo Qi,1 = Qi,0 + ∆Q. En este caso, se espera que para el
caudal corregido se cumpla i ri Qni,1i = 0
P

La corrección ∆Q se obtiene truncando una expansión de Taylor defi-

niéndose así
P ni
i ri Qi
∆Q = − P ni −1 (4.31)
i ri ni Qi

4.5. Altura Neta Positiva en la Aspiración (NPSH)

La Altura Neta Positiva de la Aspiración o Net Positive Suction Head
(NPSH) es fundamental para el diseño de sistemas de bombeo y nos indica
una condición mínima para evitar cavitación en la entrada de la bomba.

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

56 CAPÍTULO 4. FLUJOS EN TUBERÍAS

Figura 4.5: Succión de la bomba y Altura Neta Positiva de la Aspiración.

!
p2 V2 pv
NPSHr = + 2 − , (4.32)
ρg 2g ρg
donde p2 es la presión y V2 la velocidad a la entrada de la bomba como se
muestra en la Fig. 4.5; y pv es la presión de vapor del fluido el cual depende
fuertemente de la temperatura. El valor de NPSHr es un dato del fabricante
de la bomba y resumen la dinámica interna de la bomba entre el punto de
succión y el punto más desfavorable al interior.
Haciendo un balance energético entre 1 y 2 en la Fig. 4.5 podemos escri-
bir,

p1 V2 p2 V2
+ 1 + z1 = + 2 + z2 + hf , (4.33)
ρg 2g ρg 2g
donde hf son las pérdidas de carga en el tramo 1 - 2. Los casos típicos de
succión nos permiten asumir que p1 = Patm = 101,3 kPa y V1 ≈ 0 m/s para
luego despejar
!
p2 V2 Patm
+ 2 = + z1 − z2 − hf (4.34)
ρg 2g ρg
que podemos reemplazar en 4.32 para obtener NPSHd .

Patm − pv
NPSHd = + z1 − z2 − hf . (4.35)
ρg
El valor de NPSHd calculado en (4.35) corresponde al Net Positive Suc-
tion Head Disponible y para evitar el efecto de cavitación es condición ne-
cesaria y suficiente que NPSHd > NPSHd .

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

Capítulo 5

Análisis dimensional

El análisis dimensional es una muy buena herramienta en el análisis de

datos experimentales y busca resumir la información aportada por un set
de datos en los elementos fundamentales que lo componen, construyendo
curvas o conjunto de curvas que modela todo el proceso.
Para entender mejor el proceso debemos precisar los conceptos de unida-
des y dimensiones. Las unidades utilizadas en cualquier medición dependerán
del sistema de unidades en que trabajemos, ya sea el S.I. o el sistema Inglés.
En ambos casos hay una equivalencia entre las unidades de ambos sistemas.
Por otro lado, las dimensiones son las dimensiones físicas que presenta el
fenómeno de estudio y estas se reducen a siete dimensiones primarias las
cuales son Masa, Longitud, Tiempo, Temperatura, Corriente eléctrica, Can-
tidad de luz y Cantidad de materia. En la siguiente table se puede ver la
relación entre dimensiones y unidades.
Junto con las dimensiones y unidades debemos entender que son las
variables dimensionales, es decir las cantidades que varían. Las constantes
dimensionales que son la las cantidades que pueden cambiar de un experi-
mento a otro pero son constantes durante el experimento, y las constantes

Cuadro 5.1: Relación entre dimensiones y unidades

Dimensiones Símbolo Unidades S.I. Unidades Imperiales

Masa m kg lbm
Longitud L m ft
Tiempo t s s
Temperatura T K R
Corriente eléctrica I A A
Cantidad de luz C cd cd
Cantidad de materia N mol mol

57
58 CAPÍTULO 5. ANÁLISIS DIMENSIONAL

puras que no tienen dimensiones y siempre son constantes. Estas aparecen

de la manipulación matemática. Consideremos dos ejemplos de ecuaciones
donde caracterizaremos cada uno de sus elementos.

gt2
S = S 0 + V0 t + (5.1)
2
En esta expresión tenemos las variables dimensionales S y t, las cons-
tantes dimensionales S0 , V0 , g y las constantes puras 12 .
Otro ejemplo es la ecuación de Bernoulli donde

p V2
+ + gz = CT E (5.2)
ρ 2
En esta ecuación tenemos variables dimensionales p, V , z, constantes dimen-
sionales ρ, g, CT E, y las constantes puras 21 y el exponente 2.
Ambas ecuaciones (5.1), (5.2) deben respetar la ley de homogeneidad
dimensional, es decir cada témico que se suma o multiplica debe contener
las mismas dimensiones.
La derivación e integración puede cambiar las dimensiones pero no debe
cambiar la condición de homogeneidad.

5.1. Teorema de Π Buckingham

Si observamos un fenómeno que depende de n variables dimensionales,
el análisis dimensional nos permitirá reducir su descripción a m variables
adimensionales donde n − m = 1, 2, 3, 4, dependiendo de la complejidad del
problema. Una de las técnicas más famosas

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

Capítulo 6

Flujo en canales abiertos

Los canales abiertos son una de las formas más antiguas de ingeniería y
muchas de sus soluciones aun funcionan después de miles de años como por
ejemplo el acueducto que alimenta la Fontana di Trevi en Roma.
A nosotros nos interesa conocer en más detalle la dinámica de los canales
abiertos y como calcular su caudal. Los flujos en canales abiertos se caracte-
rizan por tener una superficie libre expuesta a la atmósfera, en general son
problemas tridimensionales los cuales pueden ser aproximados por modelos
en dos dimensiones para estuarios y lagos, y en una dimensión para canales
y ríos.
Es posible clasificar el flujo de acuerdo a las siguientes definiciones:

Continuo: Es flujo es independiente del tiempo.

Discontinuo: El flujo depende del tiempo.

Uniforme: El flujo es independiente de la posición.

No-Uniforme: El flujo depende de la posición.

En general nos interesa modelar y estudiar los flujos continuos no-uniformes,

es decir, donde la velocidad media del canal V = V (x) y la altura de agua
sobre la base del canal y = y(x) dependen de la posición a lo largo del canal.
Uno de los primeros elementos que se deben entender en el estudio de
los canales abiertos es la relación entre la gravedad y el flujo. Dicho de otra
manera, la gravedad es la fuerza que genera el flujo en un canal abierto. En
un canal abierto horizontal y donde la superficie libre es también horizontal,
no se generará un flujo. Adicional a la clasificación del flujo mencionado
anteriormente, se define el parámetro adimensional que relaciona las fuerzas
inerciales con la fuerza de gravedad, definiendo el número de Froude (F r).
Para el caso de canales abierto, el número F r se puede evaluar como

V
Fr = √ (6.1)
gy

59
60 CAPÍTULO 6. FLUJO EN CANALES ABIERTOS

donde V = V (x) y y = y(x) son la velocidad media y la profundidad del flujo

respectivamente en la posición x. Esta definición nos indica que el número
de Froude F r = F r(x) varía a lo largo del canal.
De acuerdo al numero F r, podemos clasificar el flujo en 3 condiciones.
Si F r > 1 las fuerzas inerciales tienen mayor peso que la fuerza de gravedad
y el flujo se llama Supercrítico. En el caso en que F r = 1 el flujo se llama
Crítico. El caso en que la fuerza de gravedad dominan sobre las fuerzas
inerciales donde F r < 1 se denomina Subcrítico.
Las leyes de conservación de la masa y energía también aplican para
canales abiertos. En particular, si estudiamos un flujo continuo no-uniforme
de un fluido de densidad constante, la conservación de la masa nos dice que
Q = V (x)A(x) = cte.
Para el balance de energía podemos plantear

p1 V2 p2 V2
+ 1 + h1 = + 2 + h2 + hf (6.2)
ρg 2g ρg 2g

donde hemos imaginado una linea de corriente sobre la superficie libre

con h1 = z1 + y1 y h2 = z2 + y2 como se observa en la Fig. 6.1. Como men-
cionamos anteriormente, un canal abierto siempre requiere que su superficie
libre esté expuesta a la atmósfera o otro ambiente de presión constante. De
este modo, p1 = p2 = patm , luego la ecuación de balance de energía se reduce
a

V12 V2
+ h1 = 2 + h2 + hf (6.3)
2g 2g

Las pérdidas por fricción se calculan de manera similar al las pérdidas

en tuberías donde

LV 2
hf = λ (6.4)
Dh 2g

donde Dh es el diámetro hidráulico que se calcula a partir del perímetro

mojado

4A
Dh = (6.5)
P
con P el perímetro mojado, perímetro de la sección del flujo que no considera
la superficie libre. El radio hidráulico se define como Rh = Dh /4 = A/P .
Ejemplo: Calcular perímetro mojado de un canal de sección recta de
base b y altura de las paredes a donde fluye agua con altura y. Respuesta
P = b + 2 y.

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 6. FLUJO EN CANALES ABIERTOS 61

Super cie libre

Fondo del canal

Figura 6.1: Perfil longitudinal de un canal abierto a la atmósfera.

6.1. Flujo uniforme en canales de pendiente y sec-

ción constante
Un canal de pendiente y sección constante significa que el fluido que
recorre el canal no observa ni percibe cambio en su comportamiento, es decir
la velocidad y altura del agua son constantes y no dependen de la posición.
En la Fig. 6.1 se observa este tipo de canal con pendiente constante dada
por θ ente los puntos x1 y x2 .
Aplicando la ley de conservación de la masa podemos decir que

Q1 = Q2 = Q (6.6)
y reconociendo que la velocidad media del flujo V y la sección de canal
donde se encuentra el fluido A no dependen de la posición, luego para todo
x1 < x < x2 ,

Q=V A (6.7)
Como el canal es de sección constante y el área es constante, la altura
del fluido y = y1 = y2 también será constante. Con estos resultado más el
balance de energía obtenemos que

h1 − h2 = hf (6.8)
z1 − z2 = hf (6.9)
(6.10)
Utilizando un argumento geométrico, podemos decir que
z1 − z2
S0 = tan (θ) = (6.11)
|x1 − x2 |

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

62 CAPÍTULO 6. FLUJO EN CANALES ABIERTOS

donde el nombre S0 es generalmente utilizado en el estudio de canales y se

considera positivo cuando el flujo desciende. De () y () escribimos

hf = S0 L (6.12)
L V02
λ = S0 |x1 − x2 | (6.13)
Dh 2g
V02
λ = S0 (6.14)
Dh 2g
V02
λ = S0 (6.15)
4Rh 2g

En esta derivación hay pequeñas aproximaciones como cos(θ) ≈ 1 si θ es

muy pequeño luego L ≈ |x1 − x2 |. Finalmente, la velocidad del flujo en el
canal se puede evaluar como
1/2 p
8g

V0 = Rh S 0 (6.16)
λ
Generalmente en un canal la geometría y rugosidad son constantes, a
 1/2
menos en secciones, luego se utiliza el coeficiente de Chèzy C = 8g
λ
obteniendo la fórmula de Chèzy, Ingeniero Francés que alrededor de 1769
planteó

V0 = C Rh S0 (6.17)
p

√
El coeficiente de Chèzy varía en el rango 30 ∼ 90 m/s. Por muchos años
se trabajo en estimar el coeficiente C para las distintas condiciones de los
canales. Algunos de los que propusieron valores son Ganguillet y Kutter en
1869, Manning en 1889, Bazin en 1897 y Powel en 1950. Equivalente mente a
el flujo en tuberías, es posible evaluar el coeficiente de fricción directamente
considerando que los canales son grandes y rugosos

14,8 Rh
 −2
λ = 2,0 log (6.18)
ǫ
La ecuación para el coeficiente de fricción en canales se deriva directa-
mente de (4.12). En casos especiales donde el fondo del canal o ríos tienen
principalmente rocas, se considera la ecuación

Rh
 −2
λ = 1,2 + 2,03 log (6.19)
d84 %
donde d84 % es el tamaño para el cual el 84 % de las rocas son menores. Es
muy común en ingeniería usar el coeficiente de Manning que en 1891 propuso

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 6. FLUJO EN CANALES ABIERTOS 63

1/2 1/6
8g R

C= ≈α h (6.20)
λ n
donde n es el coeficiente de rugosidad

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

64 CAPÍTULO 6. FLUJO EN CANALES ABIERTOS

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

Capítulo 7

Ejercicios

7.1. Flotabilidad

7.1.1. Problema 1

Un instrumento científico fue instalado en el fondo del mar para medir

las corrientes del agua. Considere el instrumento con un volumen de 10 cm3
y una masa de 100 kg. Luego de dos años de operación es necesario recuperar
el instrumento y para eso se propone usar una boya rígida esférica de masa
total 1kg. Considere que el agua de mar tiene densidad 1000 kg/m3 y la
gravedad es 9,81 m/s2.
¿Cuál es el volumen de la boya necesario para lograr una flotabilidad
neutra?

Figura 7.1: Instrumento al fondo del mar que debe ser recuperado mediante
una boya de volumen desconocido.

65
66 CAPÍTULO 7. EJERCICIOS

Solución:

Llamemos con el subíndice M las variables del instrumento sumergido,

luego tenemos MM = 100 kg la masa del instrumento, WM = MM g el peso
del instrumento y VM = 10cm3 = 1 × 10−5 m3 .
Llamemos con el subíndice O las variables de la boya usada para lograr
la flotabilidad neutra, luego tenemos MO = 1 kg la masa de la boya, WO =
MO g el peso de la boya.
La flotabilidad neutra se logra cuando el peso es igual a la fuerza de
empuje luego queremos resolver para la siguiente condición:

WM + WO = EM + EO (7.1)

El empuje calculado a partir del peso desplazado por un cuerpo y que

genera el instrumento EM y la boya EO corresponde a

EM = ρf VM g (7.2)
EO = ρf VO g (7.3)

siento VM el volumen del instrumento y VO el volumen de la boya, ambos

completamente sumergidos. De estas variables, debemos calcular VO que es
nuestra incógnita del problema. A partir de (7.1) podemos despejar EO

EO = WM + WO − EM (7.4)

Usando (7.3) y despejando VO tenemos

EO WM + WO − EM
VO = = (7.5)
ρf g ρf g
MM g + MO g − ρf VM g
VO = (7.6)
ρf g
MM + MO − ρf VM
VO = (7.7)
ρf
100 + 1 − 0,01
VO = (7.8)
1000

VO = 0,10099 m3 = 100,99 lt (7.9)

Se sugiere usar el resultado obtenido y calcular (7.1) para verificar la

flotabilidad neutra.

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 7. EJERCICIOS 67

7.1.2. Problema 2
Para mantener un instrumento oceanográfico flotando se utiliza una boya
la cual tiene la geometría de un paraboloide de revolución de altura H = 150
cm y diámetro D = 100 cm como se ve en la figura 7.2 En el punto inferior del
paraboloide se conecta el instrumento de masa M y volumen despreciable.
Considere despreciable el peso de la boya.
a) ¿Cuál es el peso máximo del instrumento de masa M que puede sostener
la boya?
b) ¿Cuál es la posición de flotación de la boya si el instrumento tiene masa
M = 250 kg?

Figura 7.2: Boya con geometría paraboloide de revolución soporta un ins-

trumento de masa M .

Solución:

La primera parte de esta pregunta se puede resolver al saber que el

volumen de un paraboloide de revolución es Vtot = π2 r 2 H siendo H la altura
del paraboloide y r = D/2 el radio en su punto más amplio. En este caso, el
empuje máximo se produce cuando la boya está completamente sumergida
luego desplazando un Vtot = 0,589 m3 . Si consideramos la densidad del fluido
desplazado ρ = 1000 kg/m3 y la gravedad g = 9,81 m/s2 tenemos un empuje
máximo de Emax = Vtot ρ g = 5778,57 N. De esta información respondemas
la pregunta a) indicando que la masa M máxima que puedo soportar la boya
es Mmax = Emax /g = 589 kg.

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

68 CAPÍTULO 7. EJERCICIOS

La segunda parte de esta pregunta requiere un poco más de desarrollo.

Veamos la forma de la parábola que describe la boya. Si definimos el origen
del sistema coordenado en el punto más bajo de la boya podemos decir que
z = α r 2 + cte siendo α y la cte parámetros por definir. Evaluado la ecuación
de la parábola en el origen tenemos que cte = 0. Evaluando la ecuación de la
parábola en su punto superior derecho podemos decir que H = α (D/2)2 lo
que nos da α = 4 H/D 2 . Ahora sabemos que en cualquier punto de flotación
de la boya debe respetar la condición z = α r 2 . Si consideramos que la masa
del instrumento es M = 250 kg, tendremos equilibrio cuando el empuje
generado por la boya y el peso generado por la masa M sean iguales, luego

E = Vsum ρ g = M g = W (7.10)
Vsum = M/ρ (7.11)
3
Vsum = 0,250 m (7.12)

Para lograr este volumen sumergido la boya debe estar a un nivel dado
por hsum = α r 2 y su volumen sumergido se calculará como

π 2 π 2
Vsum = r hsum = h (7.13)
2 2 α sum
s
Vsum 2 α
hsum = (7.14)
π
hsum = 0,977 m (7.15)

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

 CAPÍTULO 7. EJERCICIOS 69

7.2. Fuerzas sobre superficies

7.2.1. Problema 1
El estanque cilíndrico de la figura contiene un líquido fermentando de
gravedad específica SG = 1, 1 presurizado. La presión en el estanque es me-
dida por un manómetro en su parte superior y actualmente está registrando
P0 = 100 Pa. Debido al proceso de fermentación, la presión P0 puede au-
mentar. En la sección inferior, el estanque tiene una compuerta rectangular
de altura a = 100 cm y de base b = 30 cm que se utiliza para vaciar el
estanque. Esta compuerta se mantiene cerrada gracias a una fuerza externa
F = 10 kN. La altura del líquido es H = 5 m.

Calcule la magnitud de la fuerza distribuida sobre la compuerta.

Calcule la posición de acción de la fuerza equivalente puntual.

El operador del estanque necesita saber la máxima presión P0 que

puede contener el estanque antes que se abra la compuerta.

Figura 7.3: Estanque presurizado con una compuerta rectangular en la parte

inferior. La compuerta se mantiene cerrada gracias a la fuerza externa F .

Solución:
Para calcular la magnitud de la fuerza distribuida debemos evaluar la
presión Pc en el centroide de la compuerta. Sabemos que la altura del líquido
es H = 5 m, la presión en la superficie ligre el P0 = 100 Pa y que el centroide
se encuentra a mitad de altura de la compuerta a una distancia a/2 = 50
cm desde la base. Luego la presión en el centroide de la compuerta es

Cristóbal E. Castro 2022 Apuntes de Mecánica de Fluidos

70 CAPÍTULO 7. EJERCICIOS

Pc = P0 + ρ g hc (7.16)
con hc la profundidad del centroide a una distancia hc = H − a/2 = 5 −
0,5 = 4, 5 m desde la superficie libre. La magnitud de la fuerza distribuida
se calcula como FT ot = Pc A siendo A = a b m2 el área de la superficie
sumergida. En este ejercicio podemos evaluar la fuerza como

F = Pc A = (P0 +ρ g hc ) a b = (P0 +SG ρH2O g hc ) a b = 14597,85 N , (7.17)

considerando que Pc = 48659,5 Pa y A = a b = 0,3 m2 .

Para calcular el punto de acción de esta fuerza resultante debemos definir
un origen del sistema coordenado. Definiremos el origen con z = 0 en el
centroide de la superficie sumergido y positivo en la dirección contraria a la
gravedad. Bajo estas condiciones, la posición de acción de la fuerza resultante
se calcula como

Ixc
zeq = −ρ g (7.18)
Pc A
donde Ixc es el momento de inercia de la superficie de la compuerta calculado
con respecto a un eje horizontal que pasa por el centroide de la superficie
sumergida, lo que corresponde a
1 3
Ixc = b a = 0,0250 m4 . (7.19)
12
Con esta información podemos calcular que el punto de acción de la fuerza
puntual equivalente esta a zeq = −1,848 cm, es decir, 1,848 cm por debajo
del centroide de la compuerta.
Para responder a la tercera pregunta es necesario entender que la com-
puerta se mantendrá cerrada mientras el momento generado por la fuerza
externa F sea mayor a el momento generado por la fuerza equivalente ac-
tuando en el punto zeq . Luego podemos plantear

a
 
F a = (Pc A) × + zeq (7.20)
2
a Ixc
 
Fa = (Pc A) × − ρg (7.21)
2 Pc A
a
Fa = Pc A − ρ g Ixc (7.22)
2
a
Fa = (P0 + ρ g hc ) A − ρ g Ixc (7.23)
2
2 (F a + ρ g Ixc )
P0 = − ρ g hc (7.24)
Aa
P0 = 19,91kPa (7.25)

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

Capítulo 8

Resolución numérica a
problemas de mecánica de
fluidos

8.1. Colebrook
La ecuación de Colebrook 4.1.3 nos permite evaluar el coeficiente de
fricción y depende de el número de Reynolds 4.1 y la rugosidad relativa
ǫ/d. Esta es una ecuación no lineal para el coeficiente de fricción λ y his-
tóricamente se ha evaluado de acuerdo al diagrama de Moody que es una
representación gráfica de la ecuación de Colebrook. En tiempos modernos
con el uso de computadores y para cálculos más precisos, es posible evaluar
λ y aquí mostraremos dos alternativas.
Para resolver la ecuación necesitamos escribirla de la siguiente forma

1 ǫ/d 2,51
 
f (λ) = √ + 2 log + √ (8.1)
λ 3,7 Re λ
Si λ es solución de la ecuación, entonces f (λ) = 0 debe ser igual a cero,
luego es lo que se conoce como encontrar la raíz de la función.

8.1.1. Python
Python es un lenguaje de programación que nos permite implementar
algoritmos para resolver numerosos problemas en ingeniería. Pueden leer
más en este enlace. En Python podemos hacer uso de librerías como Numpy
para resolver numéricamente ecuaciones no lineales, por ejemplo, usando el
método de Newthon-Raphson. En el siguiente video podrán ver cómo imple-
mentar un algoritmo para resolver la ecuación de Colebrook numéricamente.
El algoritmo implementado en Python pueden descargarlo desde este enlace.

71
 CAPÍTULO 8. RESOLUCIÓN NUMÉRICA A PROBLEMAS DE
72 MECÁNICA DE FLUIDOS
8.1.2. Excel
En una planilla de cálculo como Excel de Microsoft es posible resolver
ecuaciones no lineales. Para eso primero es necesario activar el módulo. Esto
se logra en el menú Archivo-Opciones-Complementos aquí. Para lograr esto
en Excel introducimos en una celda la función f (λ) y en otra el valor de
λ y ejecutamos el complemento Solver. Pueden ver una explicación en el
siguiente video.

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

Capítulo 9

Constantes

9.0.1. Ejemplos
Ejemplo 1: Un agricultor

Cuadro 9.1: Constantes físicas

Nombre Símbolo Valor Unidades

Constante universal de los gases R 8.314 J / (mol K)
Constante del aire seco R 287.058 J / (kg K)

73
74 CAPÍTULO 9. CONSTANTES

Apuntes de Mecánica de Fluidos Cristóbal E. Castro 2022

Bibliografía

[1] A. W. Adamson. Physical Chemistry of Surfaces,. 5th ed., Interscience,

New York,., 1990.

[2] M. Blokh. Combining laws of mechanics and hydrostatics in non-inertial

reference frames. International Journal of Physical and Mathematical
Sciences, 9(9):513–516, 2015.

[3] Yunus A. Cengel and John M. Cimbala. Fluid Mechanics, Fundamentals

And Applications. McGraw-Hill, third edition, 2014.

[4] Frank M. White. Fluid Mechanics. McGraw-Hill, fourth edition, 1998.

75