1.

ECUACIONES FUNDAMENTALES DE LA HIDRAULICA

1.1. Aspectos generales
1.1.1. Principios básicos en el análisis
En la mecánica de fluidos los métodos de análisis consideran la
capacidad de un flujo pata transportar materia y el mecanismo por
el que cambia sus propiedades de un lugar a otro, para lo cual se
establece como axioma que en los fluidos se satisfagan los
principios básicos de la mecánica del medio continuo, a saber:
a) Conservación de la materia (principio de continuidad).
b) Segunda ley de Newton (impulso y cantidad de movimiento).
c) Conservación de la energía (primera ley de la termodinámica).
d) Segunda ley de la termodinámica.

El principio de la conservación de la materia o del transporte de

masa permite derivar la primera ecuación fundamental o de
continuidad, que admite diferentes simplificaciones de acuerdo con
el tipo de flujo de que se trate o de las hipótesis que se deseen
considerar.

La segunda ley de Newton establece la relación fundamental entre

la resultante de las fuerzas que actúan sobre una partícula y la
variación en el tiempo de la cantidad de movimiento.

De acuerdo con la forma en que se aplique, puede conducir a dos

ecuaciones:

la primera (componente escalar según el flujo) llamada dé la

energía, permite calcular las diferentes transformaciones de la
energía mecánica dentro del flujo y las cantidades disipadas en
energía calorífica que, en el caso de los líquidos, no se aprovecha.

La segunda, de tipo vectorial llamada del impulso y cantidad de

movimiento, permite determinar alguna de las fuerzas que
 producen el flujo si se conoce el cambio en la cantidad de
movimiento y las restantes fuerzas.

En la dinámica de fluidos (especialmente en el flujo de gases) el

análisis requiere, además, la inclusión de leyes termodinámicas
referentes al transporte de calor debido al flujo y, para ello, el
principio de la conservación de la energía permite derivar una
ecuación que relaciona la presión, densidad, temperatura;
velocidad, elevación, trabajo mecánico y la cantidad de calor
comunicado al flujo (o el que éste cede).

1.1.2. Flujo con potencial

Para este tipo de flujo la hipótesis consiste en tratarlo como
irrotácional lo que constituye la base de la hidrodinámica clásica,
una rama de la mecánica de fluidos que ocupó la atención de
eminentes matemáticos como Stokes, Rayleigh, Rankine, Kelvin y
Lamb.
En una gran cantidad de problemas prácticos de interés en la
hidráulica, ésta suposición puede ser aceptada debido a que el
agua posee u na viscosidad muy pequeña y se acerca a la
condición de fluido ideal. En otros problemas, es necesario
considerar los efectos viscosos y estudiar las fuerzas de fricción
originadas por la turbulencia que acompaña al movimiento.
En una gran cantidad de problemas prácticos de interés en la
hidráulica, ésta suposición puede ser aceptada debido a que el
agua posee u na viscosidad muy pequeña y se acerca a la
condición de fluido ideal. En otros problemas, es necesario
considerar los efectos viscosos y estudiar las fuerzas de fricción
originadas por la turbulencia que acompaña al movimiento.
1.1.3. Método experimental
El planteamiento racional de un experimento permite continuar,
 complementar o substituir el análisis en aquellos puntos en qué la
solución matemática se toma imposible o muy compleja, a tal grado
que, para obtenerla, sea necesario conceder hipótesis
simplificatorias; éstas, además de restar generalidad a la misma,
pueden llegar a falsear resultados al punto en que ellos no tengan
semblanza alguna con la situación real del problema.
1.2. Métodos de análisis
los métodos de análisis en la mecánica de fluidos se basan en; una
extensión (de los puntos de vista lagrangiano y euleriano, para describir
un flujo, referidos ahora a regiones dentro del mismo sobre las cuales se
satisfacen los principios fundamentales
En el análisis lagrangiano los principios básicos se aplican a una cantidad
definida de materia que ocupa cierta región del flujo y que recibe el
nombre de sistema. Éste puede cambiar de forma, posición y condición
térmica dentro del flujo, pero debe contener siempre la misma cantidad
de masa en cualquier instante que se considere.
El segundo método de análisis tiene aquí mayor aplicación; se llama
euleriano y estudia el flujo con base en el análisis de un volumen
adecuado de fluido llamado volumen de control fijo respecto de un
sistema coordenado y de forma y magnitud constantes. El contorno de
dicho volumen se llama superficie de control.
Al aplicar la ley de la conservación de la materia, al volumen de control
diferencial, se obtiene la ecuación diferencial de continuidad; si se aplica
la segunda ley de Newton, se obtiene la ecuación diferencial de Navier-
Stokes
1.3. Ecuación de continuidad
1.3.1. Principio de la conservación de la materia
De acuerdo con éste, de la masa de fluido que en la unida d de
tiempo entra a un volumen especificado dentro del flujo, una parte
se queda almacenada en su interior y el resto sale del volumen.
 Matemáticamente es preferible tratar con la cantidad neta de masa
que sale y que entra, sumadas algebraicamente; así, el principio
de la conservación de la materia, aplicado a un volumen de control
fijo completamente arbitrario dentro del flujo, se expresa en la
forma siguiente:
𝑐𝑎𝑛𝑡𝑖𝑑𝑎𝑑 𝑛𝑒𝑡𝑎 𝑑𝑒 𝑙𝑎 𝑚𝑎𝑠𝑎 𝑟𝑎𝑝𝑖𝑑𝑒𝑧 𝑑𝑒
𝑞𝑢𝑒 𝑎𝑡𝑟𝑎𝑣𝑖𝑒𝑠𝑎 𝑙𝑎 𝑠𝑢𝑝𝑒𝑟𝑓𝑖𝑐𝑖𝑒 𝑣𝑎𝑟𝑖𝑎𝑐𝑖𝑜𝑛 𝑑𝑒 𝑙𝑎 𝑚𝑎𝑠𝑎] = 0
[ ]+[
𝑑𝑒 𝑓𝑟𝑜𝑛𝑡𝑒𝑟 𝑑𝑒𝑙 𝑣𝑜𝑙𝑢𝑚𝑒𝑛 𝑐𝑜𝑛𝑡𝑒𝑛𝑖𝑑𝑎 𝑒𝑛 𝑒𝑙
𝑒𝑛 𝑙𝑎 𝑢𝑛𝑖𝑑𝑎𝑑 𝑑𝑒 𝑡𝑖𝑒𝑚𝑝𝑜 𝑣𝑜𝑙𝑢𝑚𝑒𝑛
Este principio se aplica lo mismo a un volumen de control de
tamaño diferencial que a uno finito, de lo cual se deriva la llamada
ecuación de continuidad.
1.3.2. Ecuación diferencial de continuidad

Figura 1.1. derivación de la ecuación diferencial de continuidad.

Si bien esta ecuación no tiene mucha aplicación en los problemas

de flujo unidimensional en hidráulica, aquí se presenta su
derivación para ser utilizada en los problemas de flujo con
potencial. Para obtenerla se aplica el principio de conservación de
la materia al volumen de control diferencial, mostrado en la Fig. 1.1
(de lados dx, dy, d z).
1 𝜕𝜌𝑉𝑥
𝑑𝑥,
2 𝜕𝑥
Donde él subíndice x indica la componente de la función 𝜌𝑉 según
x. De este modo, considerando positiva la masa que sale del
volumen y negativa la que entra, la cantidad neta de masa que
atraviesa estas caras es:
 1 𝜕𝜌𝑉𝑥 1 𝜕𝜌𝑉𝑥
(𝜌𝑉𝑥 + 𝑑𝑥) 𝑑𝑦 𝑑𝑧 − (𝜌𝑉𝑥 − 𝑑𝑥) 𝑑𝑦 𝑑𝑧 =
2 𝜕𝑥 2 𝜕𝑥
𝜕𝜌𝑉𝑥
= 𝑑𝑥 𝑑𝑦 𝑑𝑧
𝜕𝑥

por un razonamiento semejante, la cantidad neta de masa que

atraviesa las caras normales al eje y es:

𝜕𝜌𝑉𝑦
𝑑𝑥 𝑑𝑦 𝑑𝑧
𝜕𝑦

Y, la que atraviesa a las normales al eje Z:

𝜕𝜌𝑉𝑧
𝑑𝑥 𝑑𝑦 𝑑𝑧
𝜕𝑧

Finalmente, la rapidez de variación de la masa contenida en el

volumen elemental es

𝜕
(𝜌 𝑑𝑥 𝑑𝑦 𝑑𝑧);
𝜕𝑡

de tal manera que el principio de conservación de la masa

establece lo siguiente:

𝜕𝜌𝑉𝑥 𝜕𝜌𝑉𝑦 𝜕𝜌𝑉𝑧 𝜕

𝑑𝑥 𝑑𝑦 𝑑𝑧 + 𝑑𝑥 𝑑𝑦 𝑑𝑧 + 𝑑𝑥 𝑑𝑦 𝑑𝑧 + (𝜌 𝑑𝑥 𝑑𝑦 𝑑𝑧)
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡
=0

y, puesto que el volumen elemental escogido no cambia con el

tiempo, la ecuación anterior se puede simplificar y resulta:

𝜕𝜌𝑉𝑥 𝜕𝜌𝑉𝑦 𝜕𝜌𝑉𝑧 𝜕𝜌

+ + + = 0 … … … … . (1.1𝑎)
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡

O bien, recordando que

 𝜕𝜌𝑉𝑥 𝜕𝜌𝑉𝑦 𝜕𝜌𝑉𝑧
𝑑𝑖𝑣(𝜌𝑉) = + +
𝜕𝑥 𝜕𝑦 𝜕𝑧

la ecuación anterior también se expresa en la forma

𝜕𝜌
𝑑𝑖𝑣(𝜌𝑉) + = 0 … … … . . (1.1𝑏)
𝜕𝑡

Las Ecs. (1.1a y b) son dos formas de expresar la ecuación

diferencial de continuidad, que es la más general para un flujo
compresible no permanente; admite las, siguientes simplificaciones:

a) Flujo incomprensible no permanente

𝜕𝜌
( )=0
𝜕𝑡
𝑑𝑖𝑣(𝜌𝑉) = 0 … … … … (1.2)
b) Flujo incompresible no permanente
(𝜌 = 𝑐𝑜𝑛𝑡𝑎𝑛𝑡𝑒)
𝑑𝑖𝑣(𝑉) = 0 … … … … (1.3)
c) Flujo incompresible permanente
𝜕𝜌
(𝜌 = 𝑐𝑜𝑛𝑡𝑎𝑛𝑡𝑒, =0)
𝜕𝑡
𝑑𝑖𝑣(𝑉) = 0

problema 1.1. Un flujo incompresible permanente, con simetría axial respecto del
eje Z (Fig. 1.2), está limitado por una superficie sólida ( con la misma simetría) cuya
forma está definida por la ecuación 𝑧𝑟 2 =b ( r, radio medido desde el eje z, y b una
constante) y tiene un campo de velocidades dado por las componentes en
coordenadas cilíndricas: 𝑣𝑟 = 𝑎𝑟; 𝑣𝜃 = 0; 𝑣𝑧 = −2𝑎𝑧.

a) demostrar que se satisface la ecuación diferencial de continuidad.

b) Determinar la expresión para el gasto a través de la sección horizontal A-A y
de la sección cilíndrica B-B.
c) Determinar la velocidad en el punto P (r = z = 1.5 m ) cuando Q =10.64 m
3/seg (Ref. 20).
 Figura 1.2. flujo de problema 1.1.

Solución a). El campo de velocidades, definido en coordenadas cilíndricas,

equivale a las siguientes expresiones en coordenadas cartesianas
𝑣𝑥 = 𝑎𝑥
𝑣𝑦 = 𝑎𝑦
𝑣𝑧 = −2𝑎𝑧
Resulta entonces que
𝑑𝑖𝑣 𝑉 = 𝑎 + 𝑎 − 2𝑎 = 0
esto es, se satisface la ecuación de continuidad (1.3) y se verifica que el flujo
es incompresible.
Para los restantes puntos conviene más utilizar las coordenadas polares.
Solución b). Para la sección horizontal A-A, el gasto es
√𝑏/𝑧

√𝑏/𝑧
𝑄 = −∫ 2𝜋𝑟(−2𝑎𝑧)𝑑𝑟
0

√𝑏/𝑧
𝑟2
𝑄 = 4𝜋𝑎𝑧 [ ] = 2𝜋𝑎𝑏
2 0

Para la sección cilíndrica B - B se tiene:

𝑏/𝑟 2
𝑏/𝑟 2
𝑄=∫ 2𝜋𝑟(𝑎𝑟)𝑑𝑧 = 2𝜋𝑎𝑟 2 [𝑧]0
0
 𝑄 = 2𝜋𝑎𝑏

Solución c). para el punto P

𝑏 = 𝑧𝑟 2 = 1.5 ∗ 2.25 = 3.375
y, considerando el valor de Q, se tiene entonces que
𝑄 10.64
𝑎= = = 0.502 𝑠𝑒𝑔−1
2𝜋𝑏 2 ∗ 3.1416 ∗ 3.375
por tanto, la magnitud de la velocidad en el punto P, es:

𝑣 = √𝑣𝑟2 + 𝑣𝑧2 = 𝑎√𝑟 2 + 4𝑧 2 =

= 0.502√2.25 + 4 ∗ 2.25
𝑣 = 1.684𝑚/𝑠𝑒𝑔

1.3.3. Ecuación de continuidad para una vena líquida La vena líquida

mostrada en la Fig. 1.3 está limitada por la superficie 3 (que
generalmente coincide con una frontera sólida, o por ésta y una
superficie libre) y por las secciones transversales 1 y 2, normales
al eje que une los centros de gravedad de todas las secciones. Las
velocidades en Cada punto de una misma sección transversal
poseen un valor medio V, que se considera representativo de toda
la sección y de dirección tangencial al eje de la vena.
Se considera el volumen elemental de líquido “mostrado en la Fig.
1.3” limitado lateralmente por la superficie que envuelve a la vena
líquida, así como por dos secciones transversales normales al eje
de la vena, separadas la distancia de, donde s representa la
coordenada curvilínea siguiendo el eje de la vena. La cantidad neta
de masa que atraviesa la superficie de frontera, del volumen
elemental en estudio, es:
𝜕(𝜌𝑉𝐴) 𝜕(𝜌𝑉𝐴)
[𝜌𝑉𝐴 + 𝑑𝑠] − 𝜌𝑉𝐴 = 𝑑𝑠
𝜕𝑠 𝜕𝑠
y, la rapidez con que varía la masa dentro del mismo, es
𝜕(𝜌𝐴𝑑𝑠)/𝜕𝑡. Por tanto, el principio de conservación de la masa
establece que
 𝜕(𝜌𝑉𝐴) 𝜕
𝑑𝑠 + (𝜌𝐴𝑑𝑠) = 0 … … … … … … (1.4. )
𝜕𝑠 𝜕𝑡

Figura 1.3. ecuación de continuidad para una vena liquida.

Sin cometer prácticamente error se puede aceptar, en la mayoría

de los problemas, que la longitud de del elemento de volumen
considerado, no depende del tiempo. Éste puede salir de la
derivada del segundo término de la ecuación anterior y simplificarse
con; el que aparece en el primero, de lo cual resulta:

𝜕(𝜌𝑉𝐴) 𝜕(𝜌𝐴)
+ = 0 … … … … (1.5𝑎)
𝜕𝑠 𝜕𝑡

Recordando que 𝜌, V, A son funciones de s y t, al desarrollar las

derivadas parciales indicadas se obtiene:

𝜕𝑉 𝜕𝐴 𝜕𝜌 𝜕𝐴 𝜕𝜌
𝜌𝐴 + 𝜌𝑉 + 𝑉𝐴 +𝜌 +𝐴 = 0 … … … . (1.5𝑏)
𝜕𝑠 𝜕𝑠 𝜕𝑠 𝜕𝑡 𝜕𝑡

O bien, con V=ds/dt:

𝜕𝑉 𝜕𝐴 𝑑𝑠 𝜕𝐴 𝜕𝜌 𝑑𝑠 𝜕𝜌
𝜌𝐴 +𝜌( + )+𝐴( + ) = 0 … … … (1.5𝑐)
𝜕𝑠 𝜕𝑠 𝑑𝑡 𝜕𝑡 𝜕𝑠 𝑑𝑡 𝜕𝑡

Dividiendo la Ec. (1.5c) entre p A y recordando el desarrollo de la

derivada total, resulta entonces:

𝜕𝑉 1 𝑑𝐴 1 𝑑𝜌
+ + = 0 … … … … … (1.5𝑑)
𝜕𝑠 𝐴 𝑑𝑡 𝜌 𝑑𝑡
que es la ecuación de continuidad para una vena líquida donde se
produce un flujo no permanente y compresible. Un ejemplo clásico
de su aplicación lo constituye el problema de golpe de ariete. En
problemas de flujo no permanente a superficie libre (tránsito de
ondas de avenida en canales y de mareas en estuarios), donde se
considera que el líquido es incompresible, desaparece el último
término de la Ec. (1. 5d).

Si el escurrimiento es permanente las derivadas con respecto a t

que aparecen en la Ec. (1. 5a) se eliminan y esta ecuación ‘resulta:

𝜕(𝜌𝑉𝐴)
… … … … … … … (1.6𝑎)
𝜕𝑠

O bien,

𝜌𝑉𝐴 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 … … … … … … … . (1.6𝑏)

Si, además, el fluido es incomprensible:

𝑉𝐴 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 … … … … … … . . (1.7𝑎)

Esto significa que es constante el gasto que circula por cada

sección de la vena líquida en un flujo permanente; o bien, que para
dos secciones transversales 1 y 2 de la misma, se cumple lo
siguiente:

𝑄 = 𝑉1 𝐴1 = 𝑉2 𝐴2 … … … … … … … … … . (1.7𝑏)

Problema 1.2. En la Fig. 1.4 se muestra la bifurcación de un tubo

circular que tiene los diámetros indicados. El agua que escurre
dentro del tubo, entra en A y sale en C y D. Si la velocidad media
en B es de 0.60 m/seg, y en C es de 2.70 m/seg, calcular las
velocidades medias en A y D; el gasto total; y el gastó en cada rama
de la tubería.
Figura. 1.4. flujo en la bifurcación

Solución. La ecuación de continuidad (1. 7) aplicada a la vena

líquida, considerada en la Fig. 1.4, conduce a que:

𝜋𝐷𝐴2 𝜋𝐷𝐵2
𝑉𝐴 = 𝑉𝐵
4 4

De donde

0.30 2
𝑉𝐴 = 0.60( ) = 2.40𝑚/𝑠𝑒𝑔
0.15

En forma análoga:

𝜋𝐷𝐵2 𝜋𝐷𝐶2 𝜋𝐷𝐷2

𝑉𝐵 = 𝑉𝐶 + 𝑉𝐷
4 4 4

0.30 2 0.10 2
𝑉𝐷 = 0.60 ( ) − 2.7 ( )
0.05 0.05

𝑉𝐷 = 21.6 − 10.8

𝑉𝐷 = 10.8𝑚/𝑠𝑒𝑔

El gasto total es

𝜋𝐷𝐴2 𝜋𝐷𝐶2 𝜋𝐷𝐷2

𝑄 = 𝑉𝐴 = 𝑉𝐶 + 𝑉𝐷
4 4 4

𝑄 = 2.4 ∗ 0.785 ∗ 0.0225

𝑄 = 0.042𝑚3 /𝑠𝑒𝑔

El gasto por el tubo C es entonces:

 𝜋𝐷𝐶2
𝑄𝐶 = 𝑉𝐶 = 2.70 ∗ 0.785 ∗ 0.01
4

𝑄𝐶 = 0.021𝑚3 /𝑠𝑒𝑔

Y, el gasto por el tubo D, el siguiente:

𝜋𝐷𝐷2
𝑄𝐷 = 𝑉𝐷 = 10.8 ∗ 0.785 ∗ 0.0025
4

𝑄𝐷 = 0.021𝑚3 /𝑠𝑒𝑔

Esto es el gasto total vale

𝑄 = 𝑄𝐶 + 𝑄𝐷 = 0.021 + 0.021

𝑄 = 0.042𝑚3 /𝑠𝑒𝑔

1.4. Ecuación de la energía

1.4.1. Ecuaciones del movimiento
Si no se incluyen los efectos termodinámicos en el flujo ni la adición
o extracción de energía mecánica desde el exterior (bomba o
turbina), es posible derivar las ecuaciones del; movimiento
aplicable; al flujo de líquidos a partir de la segunda ley de Newton.
Para ello es necesario considerar las fuerzas que se oponen al
movimiento, las cuales desarrollan un trabajo mecánico
equivalente a la energía disipada al vencer dichas fuerzas.
Cuando se aplica la segunda ley de Newton a un elemento
diferencial de masa de líquido, en la forma dF= dma, se obtienen
las ecuaciones del movimiento a lo largo de una línea de corriente
para el flujo de un líquido real, no permanente;
Para el planteo de las ecuaciones es necesario establecer el
equilibrio dinámico de las fuerzas en las direcciones tangencial,
normal y binormal, que actúan sobre el elemento líquido (mostrado
en las figuras 1.6), con la fuerza de peso como única fuerza de
cuerpo. Dicho elemento encierra al punto ρ, en el cual existen los
valores v, p, ρ, 𝜏 (velocidad, presión, densidad, esfuerzo de
fricción). Las. componentes de las fuerzas que actúan sobre el
elemento en la dirección + s son las siguientes:
a) La fuerza de superficie resultante de un gradiente de presiones
en la dirección del movimiento; para la dirección positiva de la
coordenada curvilíneas (Fig. 1.6b) es:
1 𝜕𝑝 1 𝜕𝑝 𝜕𝑝
(𝑝 − 𝑑𝑠) 𝑑𝑛 𝑑𝑏 − (𝑝 + 𝑑𝑠) 𝑑𝑛 𝑑𝑏 = − 𝑑𝑠 𝑑𝑛 𝑑𝑏
2 𝜕𝑠 2 𝜕𝑠 𝜕𝑠

Figura 1.6 a. Elemento de líquido en un campo de flujo.

b) La fuerza de superficie, debida a la resistencia al movimiento,

se puede evaluar en términos del esfuerzo tangencial de
fricción 𝜏, el cual varía únicamente en la dirección n dado que
en la inmediata vecindad del punto P no hay variación de la
velocidad en la dirección b. Esta fuerza es:
1 𝜕𝜏 1 𝜕𝜏 𝜕𝜏
(𝜏 + 𝑑𝑛) 𝑑𝑠 𝑑𝑏 − (𝜏 − 𝑑𝑛) 𝑑𝑠 𝑑𝑏 = − 𝑑𝑛 𝑑𝑠 𝑑𝑏
2 𝜕𝑛 2 𝜕𝑛 𝜕𝑛
c) La componente de la fuerza de cuerpo, debida al propio peso
del elemento. Con 𝑐𝑜𝑠𝜃 = 𝜕𝑧/𝜕𝑠, vale :
𝜕𝑧
−𝜌 𝑔 𝑑𝑠 𝑑𝑛 𝑑𝑏 𝑐𝑜𝑠𝜃 = −𝜌 𝑔 𝑑𝑠 𝑑𝑛 𝑑𝑏
𝜕𝑠
Puesto que en todos los términos que representan fuerzas aparece
el volumen del elemento ds dn db, resulta entonces:
𝜕𝑝 𝜕𝜏 𝜕𝑧
[− + − 𝜌𝑔 ] 𝑑𝑠 𝑑𝑛 𝑑𝑏 =
𝜕𝑠 𝜕𝑛 𝜕𝑠
𝜕 𝑣2 𝜕𝑣
= 𝜌 [ ( ) + ] 𝑑𝑠 𝑑𝑛 𝑑𝑏
𝜕𝑠 2 𝜕𝑡
dado que 𝜌 ds dn db representa la masa del elemento, si los
términos de la ecuación anterior se dividen entre aquella, cada
término representará una fuerza por unidad de masa.

Figura 1.6 b. Componentes de las fuerzas que actúan sobre el

elemento
1 𝜕𝑝 1 𝜕𝜏 𝜕𝑧
− + −𝑔 =
𝜌 𝜕𝑠 𝜌 𝜕𝑛 𝜕𝑠
𝜕 𝑣2 𝜕𝑣
= ( )+ … … … … … (1.8𝑎)
𝜕𝑠 2 𝜕𝑡
ésta es la primera ecuación diferencial del movimiento. El primer
término es debido al gradiente de presiones en la dirección de la
línea de corriente; el segundo, la fuerza de resistencia causada por
la fricción interna y que induce la disipación de energía; el tercero,
la fuerza de peso (todas estas fuerzas son por unidad de masa);
finalmente, el cuarto término (segundo miembro) es el cambio de
energía cinética (aceleración conectiva) que experimenta la unidad
de masa a lo largo de la línea de corriente; y, el último, la
aceleración local de la misma.
 1 𝜕𝑝 𝜕𝑧 𝑣2
− −𝑔 = − … … … … … … . (1.8𝑏)
𝜌 𝜕𝑛 𝜕𝑛 𝑟
La Ec. (1. 8b) permite determinar la distribución de la presión en la
dirección de la normal principal de la línea de corriente, si se
conoce la distribución de v sobre la misma.
En el caso de que la línea de corriente sea de curvatura
despreciable (𝑟 = ∞)el segundo término de la Éc. (1. 8b) vale cero.
Finalmente, del equilibrio dinámico según la dirección de la
binormal, resultaría:
1 𝜕𝑝 𝜕𝑧
− −𝑔 = 0 … … … … … … . (1.8𝑐)
𝜌 𝜕𝑏 𝜕𝑏
La ecuación (1.8c) es válida para el flujo permanente o. no
permanente y sus términos también representan a fuerzas por
unidad de masa.
Si se trata del flujo de líquidos los efectos térmicos no tienen
influencia en 𝜌 y, además, es común que los cambios de p y 𝜏, con
la posición del punto, sean más importantes que los que pueda
experimentar 𝜌 ( aun en golpe de ariete). Por tanto, las Ecs. (4. 8)
para el flujo de líquidos se pueden escribir en la forma:
𝜕 𝑝 𝜕𝑧 𝜕 𝜏
− ( )−𝑔 + ( )=
𝜕𝑠 𝜌 𝜕𝑠 𝜕𝑛 𝜌
𝜕 𝑝 𝜕𝑣
= ( )+ … … … … … … … . . (1.9𝑎)
𝜕𝑠 𝜌 𝜕𝑡
𝜕 𝑝 𝜕𝑧 𝑣2
− ( )−𝑔 = − … … … . . (1.9𝑏)
𝜕𝑛 𝜌 𝜕𝑛 𝑟
𝜕 𝑝 𝜕𝑧
− ( )−𝑔 = 0 … … … … … (1.9𝑐)
𝜕𝑏 𝜌 𝜕𝑏
Todavía la forma vectorial de las ecuaciones del movimiento (1.9a,
b, c)
𝑝 𝜕 𝜏
−𝑔𝑟𝑎𝑑 ( + 𝑔𝑧) + ( )𝑠 =
𝜌 𝜕𝑛 𝜌
𝑣2 𝜕𝑉
= 𝑔𝑟𝑎𝑑 ( ) + 𝑟𝑜𝑡𝑉 ∗ 𝑉 + … … … … … (1.9)
2 𝜕𝑡
1.4.2. Ecuaciones del movimiento sobre una línea de corriente
Es importante el poder efectuar la integración de la Ec. (4.9a) a lo
largo de una línea de corriente. Sin embargo, debido al carácter
tensorial del esfuerzo de fricción 𝜏, dicha integración es compleja
si no se hacen consideraciones simplificatorias.
Puesto que los términos de la Ec. (4. 9a) representan fuerzas por
unidad de masa, al dividir la misma entre g dichos términos
expresarán ahora fuerzas por unidad de peso. Haciendo esta
operación con 𝛾 = 𝜌𝑔, y ordenando, resulta:
𝜕 𝑝 𝑣2 𝜕 𝜏 1 𝜕𝑣
− (𝑧 + + ) + ( )= … … … … . (1.10𝑎)
𝜕𝑠 𝛾 2𝑔 𝜕𝑛 𝛾 𝑔 𝜕𝑡
Si, además, los términos de la ecuación anterior se multiplican por
ds, los resultantes expresarán los trabajos mecánicos realizados
por. las fuerzas (por unidad de peso) a lo largo de la línea de
corriente, o bien, las energías equivalentes (también por unidad de
peso):
𝜕 𝑝 𝑣2 𝜕 𝜏 1 𝜕𝑣
− (𝑧 + + ) 𝑑𝑠 + ( ) 𝑑𝑠 = 𝑑𝑠 … … … … … (1.10𝑏)
𝜕𝑠 𝛾 2𝑔 𝜕𝑛 𝛾 𝑔 𝜕𝑡
La integración de esta ecuación sobre una línea de corriente
conduce a que:
𝑝 𝑣2 𝜕 𝜏
𝑧+ + −∫ ( ) 𝑑𝑠 =
𝛾 2𝑔 𝑠 𝜕𝑛 𝛾

1 𝜕𝑣
= 𝐶(𝑇) − ∫ 𝑑𝑠
𝑔 𝑠 𝜕𝑡
0 𝜕 𝜏
El término − ∫𝑠 ( ) 𝑑𝑠 se interpreta como la energía, por unidad
𝜕𝑛 𝛾

de peso, utilizada para vencer las fuerzas de fricción y que se

transforma en energía calorífica no aprovechable en el
movimiento. Por esta razón se considera una pérdida de energía
que se designará por hr. De esta manera, la Ec. (1.10) será:
𝑝 𝑣2 1 𝜕𝑣
𝑧+ + + ℎ𝑟 = 𝐶(𝑇) − ∫ 𝑑𝑠 … … … … . . (1.11)
𝛾 2𝑔 𝑔 𝑠 𝜕𝑡
 donde C(t) es una constante de integración que es función
únicamente del tiempo. Esta es la ecuación del movimiento pata
una línea de corriente en un flujo de un líquido real (rotacional) no
permanente;
Su forma diferencial, equivalente a la Ec. (1.10b), es
𝜕 𝑝 𝑣2 1 𝜕𝑣
(𝑧 + + + ℎ𝑟) = − … … … … … (1.12)
𝜕𝑠 𝛾 2𝑔 𝑔 𝜕𝑡
La Ec. 4.11 admite las siguientes simplificaciones:
a) Si el flujo es permanente, la integral de la Ec. (1.11) desaparece
y C(t)=𝐶1 (constante).
𝑝 𝑣2
𝑧+ + + ℎ𝑟 = 𝐶1 … … … … … … (1.13)
𝛾 2𝑔
b) Si en el flujo, además, no hay fricción la Ec. (1.13) toma la
expresión:
𝑝 𝑣2
𝑧+ + = 𝐶2 … … … … … … (1.14)
𝛾 2𝑔
Que es la ecuación de Bernoulli para una línea de corriente.
Por lo que respecta a la componente dada por la Ec. (1.9 b), es
interesante la integración para el caso en que las líneas de
corriente fuesen rectas o de curvatura despreciable, en un flujo
permanente. Para este caso, r =∞ muy grande y dicha ecuación
es:
𝜕 𝑝
( + 𝑔𝑧) = 0
𝜕𝑛 𝜌
Por tanto, la integración en la dirección de la normal a la línea de
corriente conduce a:
𝑝
+ 𝑧 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒𝑠 … … … … … … . . (4.15)
𝛾
lo cual significa que la presión se distribuye de manera hidrostática
en la dirección de la normal principal.
1.4.3. Ecuación de la energía para una vena líquida al existir una
distribución de velocidades en la sección, que además se aparta
del valor medio V (Fig. 1.7), se comete un error en el cálculo de
dicho valor medio. Puesto que en las ecuaciones (1.11) y (1.12) el
término 𝑣 2 /2𝑔 representa la energía .cinética que posee la unidad
de peso,

la que corresponde al peso del líquido qué atraviesa el área dA en

la unidad de tiempo será: 𝛾𝑣𝑑𝐴𝑣 2 /2𝑔. En la misma forma, la
energía cinética que posee todo el peso del líquido que fluye a
través de una sección de la vena líquida, en la unidad de tiempo,
es 𝛾𝑉𝐴𝛼𝑣 2 /2𝑔, donde a corrige el error de considerar el valor
medio de la velocidad. Se debe entonces satisfacer lo siguiente:
𝑣2 𝑣2
∝ 𝛾𝑉𝐴 = ∬ 𝛾𝑣𝑑𝐴
2𝑔 𝐴 2𝑔

Puesto que 𝛾 representa el valor medio del peso específico en toda

la sección, resulta que
1 𝑣 8
∝= ∬ ( ) 𝑑𝐴 … … … … … … … … (1.16)
𝐴 𝐴 𝑉
Por un razonamiento análogo con el último término de la Ec. (4.12),
se tiene

𝛽𝑉𝜌𝑉𝐴 = ∬ 𝑣𝜌𝑣𝑑𝐴
𝐴

1 𝑣 2
𝛽= ∬ ( ) 𝑑𝐴 … … … … … … … (1.17)
𝐴 𝐴 𝑉
Los coeficientes 𝛼 𝑌 𝛽 se conocen como coeficientes de Coriolis y
de BousSinesq, respectivamente. Con estas correcciones la Ec.
(4.12) resulta así:
 𝜕 𝑝 𝑉2
(𝑧 + + 𝛼 + ℎ𝑟) =
𝜕𝑠 𝛾 2𝑔
1 𝜕𝛽
=− … … … … … … … … . (1.18)
𝑔 𝜕𝑡
que es la ecuación diferencial de la energía para una vena líquida,
llamada también ecuación dinámica. Si esta ecuación se integra
entre dos secciones, 1 y 2 de la vena líquida, se obtiene:
2
𝑝1 𝑉12 𝑝2 𝑉22 1 2 𝜕(𝛽𝑉)
𝑧1 + + 𝛼1 = 𝑧2 + + 𝛼2 + ∑ ℎ𝑟 + ∫ 𝑑𝑠. . (1.19)
𝛾 2𝑔 𝛾 2𝑔 𝑔 1 𝜕𝑡
1

es decir, la ecuación general de la energía para una vena líquida,

donde ∑21 ℎ𝑟 representa la disipación de energía interna del flujo,
entre las secciones 1 y 2, que además, incluye la constante de
integración C(t).
1.5. Ecuación de la cantidad de movimiento
La ecuación de la cantidad de movimiento en un cuerpo libre o volumen
de control se deriva de la segunda ley de Newton. Se conoce como la
cantidad de movimiento de un elemento; de masa Mal producto de ésta
por su velocidad. Por tanto, la segunda ley de Newton establece lo que
sigue.
La suma vectorial de todas las fuerzas F que actúan sobre u na masa de
fluido es igual a la rapidez del cambio del vector lineal cantidad de
movimiento de la masa de fluido, es decir:
𝑑(𝑀𝑉)
𝐹= … … … … … … … … … . . (1.20)
𝑑𝑡
Las fuerzas externas son de dos tipos:
a) Fuerzas de superficie que actúan sobre la masa de fluido
b) Fuerzas de cuerpo Fc, generalmente las de peso propio.
 Figura 1.8. derivación de la ecuación de la cantidad de movimiento para
un volumen de control
La variación en el tiempo, de la cantidad de movimiento a través del
elemento dA, será entonces.

𝜌𝑉(𝑉 ∗ 𝑑𝐴)

En cualquier instante la masa de un elemento diferencial es 𝜌𝑑𝑣, donde la

densidad del elemento depende del instante que se considere y de la
posición del mismo dentro del volumen de control. La cantidad de
movimiento de dicho elemento de volumen será entonces: 𝑉𝜌𝑑𝑣.

El cambio total de la cantidad de movimiento en el tiempo, en todo el

volumen de control, será entonces:

𝑑(𝑀𝑉) 𝜕
= ∬ 𝜌𝑉(𝑉 ∗ 𝑑𝐴) + ∭ 𝑉𝜌𝑑𝑣 … … … . (1.21)
𝑑𝑡 𝑆𝐶 𝜕𝑡 𝑉𝐶

la Ec. (4.29) aplicada al volumen de fluido de la Fig. 4.12 fijo respecto de

un marco de referencia, conduce a que

𝜕
𝐹𝑝 + 𝐹𝑡 + 𝐹𝑐 = ∬ 𝑉𝜌(𝑉. 𝑑𝐴) + ∭ 𝑉𝑑𝑣 … … … . (1.22)
𝑆𝐶 𝜕𝑡 𝑉𝐶

sea, la ecuación de la cantidad de movimiento para un volumen de control

fijo.

si en esta ecuación se considera que el flujo ocurre únicamente a través

de por acciones de la superficie SC, siendo los vectores velocidad
 aproximadamente normales a la sección (con valores medios para (𝑉 𝑦 𝜌),
la primera integral de la Ec. (1.22) para cada porción de la SC, es de la
formas siguiente:

𝜕
𝐹𝑝 + 𝐹𝑡 + 𝐹𝑐 = ∑(𝜌𝑄𝛽𝑉) + ∭ 𝑉𝜌𝑑𝑣 … … … … … . (1.23)
𝜕𝑡 𝑉𝐶

Si el flujo fuese unidimensional el cuerpo libre estudiado sería como el que

se muestra en la Fig. 1.3 y la integral de la Ec. (1.23) se podría calcular
como sigue:

𝜕 𝜕 𝜕 𝜕
∭ 𝑉𝜌𝑑𝑣 = ∭ 𝑉𝜌𝑑𝐴𝑑𝑠 = ∫ 𝜌𝑑𝑠 ∬ 𝑉𝑑𝐴 = ∫ 𝜌𝑄𝑑𝑠
𝜕𝑡 𝑉𝐶 𝜕𝑡 𝑉𝐶 𝜕𝑡 𝑠 𝐴 𝜕𝑡 𝑆

y la Ec. (4.31) para el flujo unidimensional sería

𝜕
𝐹𝑝 + 𝐹𝑡 + 𝐹𝑐 = ∑(𝜌𝑄𝛽𝑉) + ∫ 𝜌𝑄𝑑𝑠 … … … … (1.24)
𝜕𝑡 𝑆

Si el flujo es permanente la integral en las Ecs. (1.23) y (1.24) vale cero.

Si además de permanente es incompresible, p es constante y la Ec. (1.24)
resulta:

𝐹𝑝 + 𝐹𝑡 + 𝐹𝑐 = 𝜌 ∑(𝑄𝛽𝑉) … … … … … … … … . (1.25)

1.6. Sobre la aplicación de las ecuaciones de la energía y de la

cantidad de movimiento.
Las ecuaciones de la energía y de la cantidad de movimiento se aplican
de manera diferente y, si se hace correctamente, ellas describirán un flujo
con idénticos grados de exactitud. Sus principales diferencias se
encuentran en su estructura: mientras la ecuación de la cantidad de
movimiento es Vectorial y engloba fuerzas totales y condiciones externas
sin tomar en cuenta los cambies internos de energía la ecuación de la
energía es por el contrario escalar y toma en cuenta los cambios internos
de energía y no las fuerzas totales y condiciones externas.
En general, cualquiera que sea el sistema de ecuaciones por usar, éste
se deberá plantear entre secciones finales con condiciones de frontera
perfectamente definidas, es decir, entre aquellas secciones de la
conducción en las que se conozcan con exactitud los valores de la
energía de posición, de presión y de velocidad y, por lo misma, la energía
total.
Estas secciones son las siguientes.
a) La superficie libre del líquido, en un recipiente al cual se conecta el
conducto.
b) La sección final de un chorro descargado por un chiflón a las
condiciones atmosféricas (o dentro de un espacio lleno de gas a
presión constante).
c) Secciones intermedias de una conducción a las cuales confluyen o se
bifurcan ramales, donde la energía sea común para todas las ramas.
si se considera que la velocidad en un punto cualquiera de la sección se
puede determinar con el valor de la media, más una fracción de la misma,
se puede escribir que:
𝑣 = 𝑉 + 𝑘 𝑉 = (1 + 𝑘 ) 𝑉
En que −1 < 𝑘 ≤ 1, siendo k una función de punto.

Figura 1.9. distribución de velocidades en una sección.

Entonces, el coeficiente de coriolis vale:
1 𝑣 3
𝛼= ∬ ( ) 𝑑𝐴 =
𝐴 𝐴 𝑉
1
= ∬ (1 − 𝐾)3 𝑑𝐴 =
𝐴 𝐴
1
= ∬ (1 + 3𝐾 + 3𝐾 2 + 𝐾 3 )3 𝑑𝐴
𝐴 𝐴
O bien
3 3 1
𝛼 = 1+ ∬ 𝑘 𝑑𝐴 + ∬ 𝑘 2 𝑑𝐴 + ∬ 𝑘 3 𝑑𝐴 … … … . . (1.26)
𝐴 𝐴 𝐴 𝐴 𝐴 𝐴
Por otra parte
1
∝= ∬ 𝑣 𝑑𝐴 =
𝑉 𝐴

= ∬ (1 + 𝑘)𝑑𝐴 = 𝐴 + ∬ 𝑘𝑑𝐴
𝐴 𝐴

Se deduce
3
𝛼 ≈1+ ∬ 𝑘 2 𝑑𝐴 … … … … … … … … (1.27)
𝐴 𝐴
En la misma forma β resulta:
3
𝛽 = 1+ ∬ 𝑘 2 𝑑𝐴
𝐴 𝐴
Combinando estas dos ecuaciones se tiene.
𝛼−1
𝛽 ≈1+ … … … … … … … … (1.28)
3
Problema 1.6. Determinar la velocidad media y los coeficientes 𝛼 𝑦 𝛽 en
un conducto cilíndrico donde se produce:

Figura 1.10. esquema ilustrativo del problema 1.6.

a) unos escurrimientos laminan cuya distribución de velocidades sigue la
ley
𝑟
𝑣 = 𝑣𝑚𝑎𝑥 [1 − ( )2 ]
𝑅
b) Un escurrimiento turbulento cuya distribución de velocidades sigue la
ley
𝑟 1 𝑦
𝑣 = 𝑣𝑚𝑎𝑥 [(1 − )7 = 𝑣𝑚𝑎𝑥 ( )1/7 ]
𝑅 𝑅
En ambos casos 𝑣𝑚𝑎𝑥 es la velocidad en el eje del tubo; R el radio
del mismo; y y =R — r la distancia a la pared de los puntos de radio r
y velocidad v (Fig. 1.10).

Solución a) La velocidad media es

𝑅
1 𝑟 2 𝑣𝑚𝑎𝑥
𝑉= 2
∫ 𝑣𝑚𝑎𝑥 [1 − ( ) ] 2𝜋𝑟𝑑𝑟 =
𝜋𝑅 0 𝑅 2
La ley de distribución de velocidades se escribe en la forma:
𝑣 𝑟 2
= 2 [1 − ( ) ]
𝑉 𝑅
El coeficiente α:
8 𝑟
𝑟 2 3
𝛼= ∫ [1 − ( ) ] 2𝜋𝑟𝑑𝑟 =
𝜋𝑅 2 0 𝑅
𝑟
𝑟 2 3 2𝑟𝑑𝑟
= −8 ∫ [1 − ( ) ] (− 2 ) = 2
0 𝑅 𝑅
De la Ec. (1.28) el valor aproximado de β es
2−1
𝛽 = 1+ = 1.33
3
Solución b) La velocidad media V, resulta de su definición, a saber:
𝑅
𝜋𝑅 2 𝑉 = 2𝜋 ∫ 𝑣𝑟𝑑𝑟
0

donde r = R — y, dr = — dy. Haciendo caso omiso del signo menos,

se tiene que:
𝑅 1
2
𝑦 7
𝜋𝑅 𝑉 = 2𝜋𝑣𝑚𝑎𝑥 ∫ (𝑅 − 𝑦) ( ) 𝑑𝑦 =
0 𝑅
𝑅
6/7 1/7
𝑦 8/7
2𝜋𝑣𝑚𝑎𝑥 ∫ (𝑅 𝑦 − 1/7 ) 𝑑𝑦 =
0 8
Resolviendo la integral resulta así:
 49
𝑉= 𝑣
60 𝑚𝑎𝑥
La ecuación para la distribución de velocidades puede expresarse
como
𝑣 49 𝑦 1/7
= ( )
𝑉 60 𝑅
El coeficiente de Coriolis a resulta de la
Ec. ( 4 . 1 6 ) :