Amd apuntes-fluidos

Agustin Martin Domingo
Apuntes de
Mecánica de Fluidos
Agustín Martín Domingo

Copyright

Esta obra “Apuntes de Mecánica de Fluidos” (texto y figuras) es:
Copyright (C) 1997-2011 Agustín Martín Domingo
con las siguientes excepciones:

• La figura 5–2 es Copyright 2006 Isaka Yoji, y se distribuye bajo licencia Creative Commons Attribution-Share Alike 2.1
Japan.
• La figura 5–9 pertenece al dominio público y proviene del informe de Henry Darcy de 1856 Les fontaines publiques de la
ville de Dijon[1].

Algunos derechos reservados.

Licencia de distribución

Este trabajo se distribuye bajo una licencia Creative Commons Reconocimiento-CompartirIgual 3.0 España (CC-
BY-SA). Para ver una copia de esta licencia, visite la página de la licencia
http://creativecommons.org/licenses/by-sa/3.0/es
o envíe una carta a Creative Commons, 171 Second Street, Suite 300, San Francisco, California, 94105, EEUU.
This work is licensed under the Creative Commons Attribution-ShareAlike 3.0 Spain License. To view a copy of
this license, visit
http://creativecommons.org/licenses/by-sa/3.0/es/
or send a letter to Creative Commons, 171 Second Street, Suite 300, San Francisco, California, 94105, USA.

Resumen de la licencia:

Está permitido. . .

• Copiar, distribuir y comunicar públicamente la obra
• Hacer obras derivadas

Bajo las siguientes condiciones

Reconocimiento: Se deben reconocer los créditos de la obra de la manera especificada por el autor o el licenciador.

Compartir bajo la misma licencia: Si se altera o se transforma esta obra, o se genera una obra derivada, sólo se
puede distribuir la obra generada bajo una licencia similar a ésta.

Indice

1. Introducción a los fluidos. La hipótesis del continuo.
1.1.
1.2.
1.3.

2.1.
Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Descripción de un fluido. Hipótesis del continuo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Propiedades de los fluidos. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.3.1.
1.3.2.
1.3.3.
1.3.4.
1.3.5.
1.3.6.
1.3.7.
Densidad ρ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Peso específico γ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Volumen específico v . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Viscosidad. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Presión. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Compresibilidad. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Dilatación térmica. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

2. Estática de fluidos.
Ecuación general de la estática de fluidos. Principio de Pascal. . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.1.1. Expresión diferencial de la misma. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5
5
6
7
7
8
8
8
8
8
9

11
11
11
2.1.2. Caso particular: fluido incompresible. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.1.3. Principio de Pascal. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.1.4. Paradoja de Pascal. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.2. Tubos en U y manómetros. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.2.1. Fluidos miscibles. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.2.2. Fluidos no miscibles. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.2.3. Manómetros. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.2.4. Manómetro diferencial. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.3. Variación de la presión con la altura en un gas perfecto en reposo. . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.3.1. Atmósfera isoterma . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.3.2. Dependencia lineal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.3.3. Atmósfera tipo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.4. Fuerzas hidrostáticas sobre superficies. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.4.1. Fuerza sobre superficies planas. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.4.2. Fuerza sobre superficies curvas. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.5. Equilibrio de un cuerpo sumergido. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.5.1. Flotación. Principio de Arquímedes. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.5.2. Estabilidad de un cuerpo semisumergido . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22

3. Dinámica de fluidos perfectos. 25
3.1. Movimiento de un fluido perfecto. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25

1

2 Indice

3.2.

3.3.

3.4.
3.5.
3.1.1.
3.1.2.
3.1.3.

3.2.2.

3.3.1.
3.3.2.
3.3.3.
3.3.4.
3.3.5.

3.5.1.
3.5.2.
3.5.3.
3.5.4.
Aproximaciones clásicas al estudio de los fluidos. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Tipos de flujo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Senda, línea de corriente y de traza. Tubo de flujo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Ecuación de continuidad. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.2.1. Forma integral. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Forma diferencial. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Fuerza y aceleración en un elemento de fluido. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Fuerza neta sobre un elemento de un tubo de flujo. Variación de la cantidad de movimiento. . . . . . . . .
Ecuaciones de movimiento microscópicas. Aceleración. Ecuación de Euler. . . . . . . . . . . . . . . . .
Aplicación a un fluido en reposo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Aplicación a un fluido en reposo en un sistema de referencia no inercial . . . . . . . . . . . . . . . . . .
Vorticidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Conservación de la energía. Ecuación de Bernoulli. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Aplicaciones de las ecuaciones anteriores. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Velocidad de salida por un orificio (teorema de Torricelli). . . . . . . . . . . . . . . . . . . . . . . . . .
Tiempo de vaciado de un depósito. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Presión dinámica y presión estática. Tubo de Pitot . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Efecto Venturi. Venturímetro. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

4. Dinámica de fluidos reales. Flujo de fluidos en tuberías.
25
26
26
28
28
29
31
31
32
33
33
36
37
40
40
41
41
42

45
4.1. Introducción. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
4.1.1. Concepto de viscosidad. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
4.1.2. Régimen laminar y régimen turbulento. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
4.2. Efecto de la viscosidad en los fluidos reales. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
4.2.1. Fluidos newtonianos. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
4.2.2. De qué depende la viscosidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
4.2.3. La condición de no deslizamiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
4.2.4. Fluidos no newtonianos. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
4.3. Ecuaciones de movimiento de los fluidos reales. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
4.3.1. Fluido perfecto o flujo no viscoso. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
4.3.2. Fluido newtoniano. Las ecuaciones de Navier-Stokes. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
4.4. Distribución de velocidades y tensiones cortantes. Fórmula de Poiseuille. . . . . . . . . . . . . . . . . . . . . . . 54
4.5. Número de Reynolds. Concepto de capa límite. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
4.6. Pérdidas de carga. Generalidades. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
4.6.1. Pérdidas de carga lineales. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
4.6.2. Pérdidas de carga singulares. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
4.6.3. Variación de altura debida a una máquina. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
4.6.4. Representación gráfica de la pérdida de carga. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
4.7. Redes de tuberías . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
4.8. Cavitación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65

Indice 3

4.9.

5.1.

5.2.
Golpe de ariete o choque hidráulico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

5. Mecánica de fluidos en medios permeables.
Fenómenos de superficie: tensión superficial y capilaridad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.1.1.
5.1.2.
5.1.3.

5.2.3.
5.2.4.
5.2.5.
5.2.6.
5.2.7.
5.2.8.
5.2.9.
Tensión superficial. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Ángulo de contacto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Capilaridad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Fluidos en medios permeables. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.2.1.
5.2.2.
Velocidad de descarga y velocidad microscópica promedio. . . . . . . . . . . . . . . . . . . . . . . . . .
La ley de Darcy en una dimensión. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Potencial hidráulico y sus componentes. El campo hidráulico. . . . . . . . . . . . . . . . . . . . . . . . .
La ley de Darcy en tres dimensiones. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Flujo a través de terrenos saturados homogéneos e isótropos. . . . . . . . . . . . . . . . . . . . . . . . .
La ecuación de Laplace y su resolución. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Flujo permeable a través de medios anisótropos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
La red de flujo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Ejemplos de redes de flujo en distintos sistemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.2.10. Subpresiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

Referencias
65

69
69
69
73
74
76
78
79
83
84
85
85
87
90
90
92

95

Indice
4


Capítulo 1

Introducción a los fluidos. La hipótesis del
continuo.

Índice del capítulo
1.1.
1.2.
1.3.
Introducción . . . . . . . . . . . . . . . . . . . .
Descripción de un fluido. Hipótesis del continuo.
Propiedades de los fluidos. . . . . . . . . . . . .
1.3.1. Densidad ρ . . . . . . . . . . . . . . . . .
1.3.2. Peso específico γ . . . . . . . . . . . . . .
1.3.3. Volumen específico v . . . . . . . . . . . .
1.3.4. Viscosidad. . . . . . . . . . . . . . . . . .
1.3.5. Presión. . . . . . . . . . . . . . . . . . . .
1.3.6. Compresibilidad. . . . . . . . . . . . . . .
1.3.7. Dilatación térmica. . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
5
6
7
7
8
8
8
8
8
9

1.1. Introducción

La Mecánica de Fluidos estudia las leyes del movimiento de los fluidos y sus procesos de interacción con los cuerpos
sólidos. La Mecánica de Fluidos como hoy la conocemos es una mezcla de teoría y experimento que proviene por un
lado de los trabajos iniciales de los ingenieros hidráulicos, de carácter fundamentalmente empírico, y por el otro del
trabajo de básicamente matemáticos, que abordaban el problema desde un enfoque analítico. Al integrar en una única
disciplina las experiencias de ambos colectivos, se evita la falta de generalidad derivada de un enfoque estrictamente
empírico, válido únicamente para cada caso concreto, y al mismo tiempo se permite que los desarrollos analíticos
matemáticos aprovechen adecuadamente la información experimental y eviten basarse en simplificaciones artificiales
alejadas de la realidad.
La característica fundamental de los fluidos es la denominada fluidez. Un fluido cambia de forma de manera continua
cuando está sometido a un esfuerzo cortante, por muy pequeño que sea éste, es decir, un fluido no es capaz de soportar
un esfuerzo cortante sin moverse durante ningún intervalo de tiempo. Unos líquidos se moverán más lentamente que
otros, pero ante un esfuerzo cortante se moverán siempre. La medida de la facilidad con que se mueve vendrá dada
por la viscosidad que se trata más adelante, relacionada con la acción de fuerzas de rozamiento. Por el contrario en
un sólido se produce un cambio fijo γ para cada valor de la fuerza cortante aplicada. En realidad algunos sólidos
pueden presentar en cierto modo ambos comportamientos, cuando la tensión aplicada está por debajo de un cierto

F F

γ

Sólido Líquido

Figura 1–1 Diferencias entre el comportamiento de un líquido y de un sólido frente a una fuerza cortante aplicada.

5

6 Capítulo 1. Introducción a los fluidos. La hipótesis del continuo.

umbral presenta el comportamiento habitual, mientras que por encima de un cierto umbral el sólido puede plastificar,
produciéndose una deformación más continua para una fuerza fija, de forma parecida a como ocurre en un fluido. Ésto
es precisamente lo que ocurre en la zona de fluencia. Si la fuerza persiste, se llega a la rotura del sólido.
Así, mientras que un sólido experimenta un desplazamiento definido (o se rompe por completo) bajo la acción de
una fuerza cortante, en los fluidos pequeñas fuerzas producen grandes deformaciones no elásticas (en general no se
recupera la forma) a volumen constante, que se realizan de forma continua. Mientras que para un sólido bajo una
fuerza cortante constante se alcanza un ángulo de deformación determinado y constante, en un fluido debemos hablar
de una velocidad de deformación constante o no, ya que la deformación se produce de forma continua.
Dentro de los fluidos, la principal diferencia entre líquidos y gases estriba en las distintas compresibilidades de los
mismos.

Gases. Los gases presentan una gran compresibilidad, que influye sobre las características del flujo, ya que tanto el
volumen como la densidad varían con facilidad. En el caso de los gases el movimiento térmico vence a las
fuerzas atractivas y, por tanto tienden a ocupar todo el volumen del recipiente que los contiene.

Líquidos. En el caso de los líquidos, por el contrario, la compresibilidad es muy débil. Ésto es debido a que las
fuerzas atractivas entre las moléculas del líquido vencen al movimiento térmico de las mismas, colapsando las
moléculas y formando el líquido. Al contrario que en el caso de los gases, que tendían a ocupar todo el volumen
que los contiene, los líquidos tienden a formar una superficie libre.

La noción de compresibilidad dada es la correspondiente a la estática de fluidos. En dinámica de fluidos, hay casos
en los que la densidad no varía a lo largo del flujo, incluso en un fluido compresible, por lo que a ese flujo le pode-
mos aplicar las leyes de los fluidos incompresibles. Tiene en este caso más sentido hablar de flujo compresible o
incompresible.
Los fluidos no conservan la forma. Al situarlos en un recipiente toman la forma del mismo (si lo llenan) o de parte
del mismo.

1.2. Descripción de un fluido. Hipótesis del continuo.

Para la descripción del movimiento de un fluido recurriremos a las leyes generales de la Mecánica (leyes de Newton,
leyes de conservación de la cantidad de movimiento y de la energía), junto con relaciones específicas condicionadas
por la fluidez.
A escala microscópica la materia, y en particular un fluido está compuesta de moléculas a cierta distancia promedio
con espacio vacío entre ellas. Estas moléculas están continuamente moviéndose y colisionando entre sí. Un análisis
exacto del problema debiera tener en cuenta la acción de cada molécula o grupo de moléculas en el fluido. Este proce-
dimiento, con algunas simplificaciones importantes es el que se adopta en Teoría Cinética y en Mecánica Estadística,
pero es aún demasiado complejo para utilizarlo en el trabajo diario de hidráulica.
En la mayor parte de los cálculos hidráulicos, el interés está realmente centrado en manifestaciones macroscópicas
promedio que resultan de la acción conjunta de una gran cantidad de moléculas, manifestaciones como la densidad,
la presión o la temperatura. En la práctica es posible hacer una simplificación importante, suponer que todas estas
manifestaciones son el resultado de la acción de una hipotética distribución continua de materia, a la que denomi-
naremos el continuo, o el medio continuo, en lugar de estudiar el conglomerado real de las moléculas discretas, de
mucha mayor complejidad. De esta forma a la hora de estudiar nuestros problemas sustituiremos la materia real por
este medio continuo ficticio, cuyas propiedades varían de forma continua y reflejan las propiedades macroscópicas del
medio real.
Este concepto del medio continuo permite una gran simplificación en el análisis. Por supuesto, este enfoque debe
utilizarse únicamente cuando arroje resultados razonablemente correctos. Por ejemplo, no puede utilizarse cuando el
recorrido libre medio de las moléculas es del orden de las magnitudes características del problema. En estas condicio-
nes, la acción de cada molécula individual es significativa y debe estudiarse individualmente.

1.3. Propiedades de los fluidos. 7

Por ejemplo, consideremos la acción sobre una superficie de la pared en el caso de un depósito cerrado que contiene
un gas a una cierta presión, en un estado estacionario. Incluso a baja presión, la gran cantidad de colisiones de molé-
culas sobre la superficie da lugar a una fuerza global que en la práctica puede considerarse independiente del tiempo,
comportamiento que será correctamente simulado por nuestro hipotético medio continuo. Ahora bien, si la presión
fuera tan baja que únicamente quedaran en el tanque unas pocas moléculas de forma que el recorrido libre medio de
las mismas es del orden de magnitud del elemento considerado, se observará una actividad errática según las moléculas
individuales o los grupos de moléculas bombardean la superficie y no se podrá hablar de una fuerza constante, sino de
una serie de choques aleatorios contra la superficie. Este comportamiento no podría ser reflejado por nuestro medio
continuo. Lo mismo ocurriría si considerando el gas discreto real, tomamos una superficie muy pequeña, de forma que
su tamaño es del orden del recorrido libre medio de las moléculas.
Sin embargo, si ya estamos trabajando con el medio continuo y con magnitudes macroscópicas, un elemento de
volumen infinitesimal será un elemento de volumen del medio continuo y no de la materia real discreta, con lo que
trabajaremos con la confianza de que a partir del mismo obtendremos las magnitudes macroscópicas.

Hipótesis del continuo:

La materia y las propiedades físicas asociadas a la misma se consideran dispersas de forma continua
en ella, y no concentradas en pequeñas fracciones (átomos y moléculas) de la misma.

De este modo, sustituimos la materia real de carácter discreto por una materia ficticia continua cuyas propiedades
en cada punto vienen dadas por las propiedades promedio de la materia real en el entorno de ese punto.

1.3. Propiedades de los fluidos.

1.3.1. Densidad ρ

Se define como la masa por unidad de volumen. Sus unidades en el sistema internacional son [kg/m3 ].
Para un fluido homogéneo, la densidad no varía de un punto a otro y puede definirse simplemente mediante

V
ρ= (1–1)
m

Por el contrario, para un fluido inhomogéneo, la densidad ρ varía de un punto a otro. Por tanto tenemos que definir
la densidad en un punto como la masa por unidad de volumen en un elemento diferencial de volumen† en torno a ese
punto:
dm
ρ = ρ(x, y, z, t) = (1–2)
dV

Esto es posible gracias a la continuidad. En los líquidos, al tener baja compresibilidad, la densidad depende de la
temperatura, pero apenas depende de la presión, ρ = ρ(T ). Para los fluidos compresibles, la densidad depende en
general tanto de la presión como de la temperatura, ρ = ρ(p, T ). Para el caso concreto de un gas ideal, con una
ecuación de estado pV = nRT , la densidad tiene la forma concreta:

Mp
ρ(p, T ) = (1–3)
RT

† Nótese que, aunque sea un elemento infinitesimal de volumen, contiene una gran cantidad de partículas y podemos definir propiedades

promedio.


1.3.2. Peso específico γ

El peso específico se define como el peso por unidad de volumen. En el sistema internacional sus unidades son
[N/m3 ]. Para un fluido homogéneo γ = mg/V = ρg, mientras que para un fluido inhomogéneo,

donde g es la aceleración de la gravedad.

1.3.3. Volumen específico v
γ = γ(x, y, z, t) = g

v = v(x, y, z, t) =
dV
dm
dm
dV

= .
= ρg

Se denomina volumen específico al volumen ocupado por la unidad de masa. Para un fluido homogéneo se define
como v = V /m = 1/ρ, mientras que en el caso general de un fluido inhomogéneo tendremos que hablar de su valor
en un punto,
1
ρ
En todos los casos, v = 1/ρ. Sus unidades en el sistema internacional son [m3 /kg].

1.3.4. Viscosidad.
(1–5)
(1–4)

Como se ha dicho en la introducción, la viscosidad refleja la resistencia al movimiento del fluido y tiene un papel
análogo al del rozamiento en el movimiento de los sólidos. La viscosidad está siempre presente en mayor o menor
medida tanto en fluidos compresibles como incompresibles, pero no siempre es necesario tenerla en cuenta. En el caso
de los fluidos perfectos o no viscosos su efecto es muy pequeño y no se tiene en cuenta, mientras que en el caso de los
fluidos reales o viscosos su efecto es importante y no es posible despreciarlo. En el caso del agua a veces se habla del
flujo del agua seca para el flujo no viscoso del agua y del flujo del agua mojada para el flujo viscoso.

1.3.5. Presión.

La presión en un punto se define como el valor absoluto de la fuerza por unidad de superficie a través de una pequeña
superficie que pasa por ese punto y en el sistema internacional su unidad es el Pascal (1 Pa=1 N/m2 ). Mientras que
en el caso de los sólidos en reposo, las fuerzas sobre una superficie pueden tener cualquier dirección, en el caso de
los fluidos en reposo la fuerza ejercida sobre una superficie debe ser siempre perpendicular a la superficie, ya que
si hubiera una componente tangencial, el fluido fluiría. En el caso de un fluido en movimiento, si éste es no viscoso
tampoco aparecen componentes tangenciales de la fuerza, pero si se trata de un fluido viscoso sí que aparecen fuerzas
tangenciales de rozamiento.
De este modo, un fluido en reposo a una presión p ejerce una fuerza −pdS sobre cualquier superficie plana arbitraria
en contacto con el fluido en el punto, definida por un vector unitario dS, perpendicular a la superficie. En general, la
presión en un fluido depende del punto, p = p(x, y, z). Así, para un fluido en reposo la presión se define como la
fuerza normal por unidad de superficie.

1.3.6. Compresibilidad.

Se caracteriza por el coeficiente de compresibilidad, κ, definido como
1 dV
κ=− (1–6)
V dp
que representa la disminución relativa del volumen por unidad de aumento de presión. Sus unidades son de inversa de
presión, en el sistema S.I. [m2 /N]. Su inversa,
1
K= (1–7)
κ
es el módulo de compresibilidad [N/m2 ]. Tanto κ como K dependen de la forma en que se realiza el proceso.

1.3. Propiedades de los fluidos. 9

1.3.7. Dilatación térmica.

Se caracteriza por el coeficiente de dilatación de volumen, que representa el aumento relativo del volumen producido
por un aumento de la temperatura, y está definida como

αV =
1 dV
V dT
(1–8)

donde V es el volumen inicial del líquido. Sus unidades son de inversa de grados [K−1 ] o [◦ C−1 ] y depende de la
forma en que realiza el proceso.

Capítulo 2

Estática de fluidos.

2.1.

2.2.

2.3.

2.4.
Ecuación general de la estática de fluidos. Principio de Pascal. . . .
2.1.1. Expresión diferencial de la misma. . . . . . . . . . . . . . . .
2.1.2. Caso particular: fluido incompresible. . . . . . . . . . . . . .
2.1.3. Principio de Pascal. . . . . . . . . . . . . . . . . . . . . . . .
2.1.4. Paradoja de Pascal. . . . . . . . . . . . . . . . . . . . . . . .
Tubos en U y manómetros. . . . . . . . . . . . . . . . . . . . . . .
2.2.1. Fluidos miscibles. . . . . . . . . . . . . . . . . . . . . . . .
2.2.2. Fluidos no miscibles. . . . . . . . . . . . . . . . . . . . . . .
2.2.3. Manómetros. . . . . . . . . . . . . . . . . . . . . . . . . . .
2.2.4. Manómetro diferencial. . . . . . . . . . . . . . . . . . . . . .
Variación de la presión con la altura en un gas perfecto en reposo. .
2.3.1. Atmósfera isoterma . . . . . . . . . . . . . . . . . . . . . . .
2.3.2. Dependencia lineal . . . . . . . . . . . . . . . . . . . . . . .
2.3.3. Atmósfera tipo . . . . . . . . . . . . . . . . . . . . . . . . .
Fuerzas hidrostáticas sobre superficies. . . . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
11
11
13
13
13
14
14
14
15
15
16
17
17
17
18
2.4.1. Fuerza sobre superficies planas. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.4.2. Fuerza sobre superficies curvas. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.5. Equilibrio de un cuerpo sumergido. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.5.1. Flotación. Principio de Arquímedes. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.5.2. Estabilidad de un cuerpo semisumergido . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22

Comenzaremos estudiando el caso del fluido estático, en el cual todas las partículas tienen (en movimiento prome-
dio) una velocidad nula.† En este fluido en reposo no aparecen fuerzas tangenciales que darían lugar a un movimiento
tangencial. Así, las fuerzas de contacto ejercidas por un fluido en reposo sobre la pared del recipiente que lo contiene
son siempre perpendiculares a la misma. Asimismo, si introducimos un cuerpo cualquiera en el interior de un fluido,
se ejercerá una fuerza sobre toda la superficie del cuerpo en contacto con el fluido y esta fuerza será normal a la
superficie en cada punto. La fuerza por unidad de superficie dependerá del punto de la superficie, pero para un punto
dado no dependerá de la orientación de la superficie del cuerpo en ese punto. Denominaremos presión ejercida sobre
la superficie a la magnitud de esta fuerza normal por unidad de superficie. La presión es una magnitud escalar.

2.1. Ecuación general de la estática de fluidos. Principio de Pascal.

2.1.1. Expresión diferencial de la misma.

El caso más general

En el caso más general, la fuerza de presión por unidad de volumen ejercida por el fluido exterior a un pequeño
elemento de volumen sobre el contenido del mismo, vendrá dada por el gradiente de la presión, ∇p (véase la sección
3.3.2). Si no hay otras fuerzas como la de la gravedad, la presión debe compensarse para que haya equilibrio.
† Rigurosamente, se mueven (en movimiento promedio) con una velocidad constante con respecto de un sistema de referencia inercial.

11

12 Capítulo 2. Estática de fluidos.

Agustin Martin Domingo x
z

z
dz
dy

dx

y

Figura 2–1 Obtención de la ecuación general de la estática de fluidos.

Consideremos una situación en la que estas fuerzas adicionales pueden describirse mediante una energía potencial,
con Φ el potencial o energía potencial por unidad de masa. La fuerza por unidad de masa vendrá dada en esas circuns-
tancias por −∇Φ y si ρ es la densidad del fluido, la fuerza por unidad de volumen por −ρ∇Φ. La suma de las fuerzas
que actúan sobre el elemento de volumen debe ser cero en el equilibrio estático, y por tanto se cumplirá la ecuación

−∇p − ρ∇Φ = 0 (2–1)

que es la ecuación general de la hidrostática. Esta ecuación en el caso general no tiene solución. Si la densidad varía
en el fluido de forma arbitraria no es posible que las fuerzas estén equilibradas y el fluido no puede estar en equilibrio
estático, estableciéndose corrientes de convección. Esto es así ya que ∇p es un gradiente puro mientras que el término
ρ∇Φ no lo es para una densidad ρ variable en general. El equilibrio hidrostático es posible únicamente en el caso de
ρ constante o cuando ρ es función únicamente de p. Veremos a continuación algunos casos en los que esto es así.

La ecuación de la hidrostática para ρ = ρ(z)

Calculemos ahora la dependencia de la presión con la altura en un fluido sometido a la atracción gravitatoria. Para
ello consideraremos un pequeño elemento prismático infinitesimal de volumen dV cuyas caras son perpendiculares a
los ejes coordenados y tienen dimensiones dx, dy y dz, con su cara inferior a una altura z. La fuerza que actúa sobre
la cara inferior del elemento prismático es obviamente

dF (z) = p(z)dxdyuz

mientras que la fuerza que actúa sobre la cara superior es

dF (z + dz) = p(z + dz)dxdy(−uz ) = − [p(z) + dp] dxdyuz

La fuerza de presión neta hacia arriba en la dirección vertical es:

dFp = dF (z) + dF (z − dz) = −dpdxdyuz

2.1. Ecuación general de la estática de fluidos. Principio de Pascal. 13

que a su vez debe ser igual al peso del fluido contenido dentro del elemento prismático ρgdxdydz, de donde se obtiene
la ecuación general de la estática de fluidos en forma diferencial,

que se escribe como
0
dp = −ρgdz

donde el signo negativo refleja el hecho de que la presión disminuye al aumentar la altura.

2.1.2. Caso particular: fluido incompresible.

En el caso de un fluido incompresible la densidad no depende de la presión, y además se tiene una superficie libre.
Es posible describir el problema en función de la profundidad respecto de la superficie libre h, que aumenta según z
disminuye. La ecuación general de la estática quedaría entonces, en forma diferencial, como

dp = −ρgdz = ρgdh

Esta ecuación se puede integrar con facilidad entre la superficie libre y un punto a una profundidad h, quedando
h
dp = p − patm =
0
h

p(h) = patm + ρgh = patm + γh
ρgdh = ρgh

y que es la ecuación general de la estática en forma integral para un fluido incompresible. En concreto, entre dos
puntos del fluido a distinta profundidad, la diferencia de presiones es

p2 − p1 = ρg(h2 − h1 ) = −ρg(z2 − z1 )
(2–2)

(2–3)

Estas expresiones subrayan el importante resultado de que la presión de un fluido en reposo en un punto del mismo
depende exclusivamente de la profundidad a la que se encuentra dicho punto, y de que la diferencia de presiones entre
dos puntos depende exclusivamente de la diferencia de profundidades de los mismos (y, obviamente, de la densidad
del fluido).

2.1.3. Principio de Pascal.

El resultado anterior nos lleva de forma inmediata al principio de Pascal,

En un fluido incompresible, las variaciones locales de presión se transmiten íntegramente a todos los
puntos del fluido y en todos los sentidos, así como a las superficies en contacto con el fluido.

o dicho con otras palabras, si la presión en un punto aumenta, por ejemplo por la aplicación de una fuerza externa,
aumenta por igual en los demás puntos del fluido, ya que según (2–3) la diferencia de presiones entre dos puntos
determinados depende únicamente de ∆h para un fluido estático.
El hecho de que los líquidos sean compresibles (aunque muy poco) hace que este cambio de presión se transmita
como una onda hasta que se restablece el equilibrio mecánico, y entonces ya se cumpliría el principio de Pascal. Este
mismo principio se cumple también para fluidos compresibles, una vez alcanzado el equilibrio.

2.1.4. Paradoja de Pascal.

A primera vista puede resultar algo sorprendente el hecho de que para un fluido dado, la presión dependa exclusi-
vamente de la profundidad y no de otras cosas como el tamaño y forma del recipiente, o pueden resultar algo extrañas
situaciones como la siguiente, conocida como paradoja de Pascal.
Consideremos los tres recipientes de idéntica base que se presentan en la figura 2–2. Los tres están llenos de agua
hasta el mismo nivel, pero sus formas son muy distintas, uno tiene la parte superior muy cerrada, uno es cilíndrico y


Agustin Martin DomingoS

(a)
S

(b)
S

(c)

Figura 2–2 Paradoja de Pascal. A pesar de que la cantidad de fluido en los recipientes sea claramente distinta, la presión en la base (que es la
misma para los tres) es la misma.

el otro tiene la parte superior muy abierta, pero todos ellos tienen la misma base.† Puede resultar algo sorprendente el
hecho de que en los tres casos la fuerza ejercida sobre la base sea la misma, como se desprende de la ecuación general
de la estática de fluidos. Sin embargo, no es tan sorprendente si se tiene en cuenta que las paredes del recipiente ejercen
sobre el líquido una fuerza perpendicular a las mismas, que puede tener una componente vertical neta bien hacia abajo
(primer caso), nula (segundo caso) o bien hacia arriba (tercer caso), componente vertical que es necesario tener en
cuenta, y que en el primer caso añade un término adicional al peso del fluido, en el segundo caso no afecta y en el
tercer caso aminora el efecto del peso del fluido.

2.2. Tubos en U y manómetros.

2.2.1. Fluidos miscibles.

Consideremos el caso de un sistema de vasos comunicantes (o un tubo en U) que contiene un fluido de densidad ρ.
Si denominamos 0A y 0B a dos puntos en la base que se encuentran a la misma altura y sin ningún obstáculo entre
medias, está claro que la presión en ambos puntos es la misma, p0A = p0B . Sin embargo queda la duda acerca de
lo que ocurre con las presiones en los puntos M y N , que se encuentran a una altura d por encima de 0A y 0B y
que no están comunicados de una forma tan clara y directa como lo están los puntos 0A y 0B . Sin embargo, se ve
enseguida que si la diferencia de presiones entre M y 0A es pM − p0A = −ρgd y la diferencia de presiones entre M
y 0B es pM − p0B = −ρgd, al tener en cuenta que p0A = p0B , se obtiene directamente que para un fluido en reposo,
pN = pM .

2.2.2. Fluidos no miscibles.

Las cosas no son tan claras sin embargo, si se tienen fluidos no miscibles. Consideremos el caso de la figura en el
que se tienen dos fluidos no miscibles, de forma que el fluido en la columna B es más denso que el fluido en la columna
A. Sea C la línea que pasa por la interfase entre los dos fluidos y une dos puntos a la misma altura. Sean M y N dos
puntos a la altura del nivel superior de la columna B. Si las dos columnas están abiertas a la atmósfera, las presiones
† Lo mismo ocurre en el caso de la fuerza sobre una pequeña superficie en la base de una presa llena de agua hasta un cierto nivel y la fuerza

sobre una superficie similar en un tubo vertical de sección la de la superficie, lleno de agua hasta la misma altura de la presa.

replacementsen U y manómetros.
2.2. Tubos 15

d
M

0A
A
pA

h

ρ
pB
B

N

0B
CA
hA

0A
A
p

M

ρA
dA

ρB
p

respectivas en M y N serán pN = patm y pM = patm + γA dA , por lo que claramente pM > pN . En líquidos no
miscibles las presiones a la misma altura son en general distintas.
B

N

hB

0B
CB

Figura 2–3 Presión en un tubo en U. a) caso de fluidos miscibles. b) Caso de fluidos no miscibles.

El esquema de la figura nos permite obtener la densidad relativa de los dos fluidos. Como las presiones en la línea
C son las mismas para la columna A que para la B, al corresponder a un mismo fluido, se tiene que de pCA =
patm + γA hA y pCB = patm + γB hB se obtiene γA hA = γB hB , lo que permite obtener la densidad relativa en
función de la altura de las distintas columnas como
γA ρA hA
= = (2–4)
γB ρB hB

2.2.3. Manómetros.

Consideremos ahora el caso de un manómetro que queremos utilizar para medir la presión en un recipiente A en la
forma que se muestra en la figura. Como a la altura MN hay un único fluido, pN = pM . Como la presión en M viene
dada por pM = pA + γA d y la presión en N por pN = patm + γh, se tiene que la presión en el depósito A vendrá
dada por
pA = patm + γh − γA d (2–5)
que en el caso habitual de un fluido manométrico mucho más denso que el fluido en A se reduce a

pA = patm + γA h (2–6)

Esta aproximación es en general válida si γ ≫ γA . En caso contrario es necesario utilizar la primera de las expresio-
nes. Normalmente con el mercurio como fluido manométrico la aproximación va bastante bien, a no ser que se quiera
medir la presión de un fluido muy denso, ya que la densidad relativa del mercurio respecto del agua es alta, de 13,6.

2.2.4. Manómetro diferencial.

Un razonamiento análogo se puede hacer para el caso del manómetro diferencial de la figura, con una presión
pM = pA + γA dA + γh en M, y una presión pN = pB + γB en N, lo que da una diferencia de presiones entre A y B

pB − pA = γh + γA dA − γB dB (2–7)

que para el caso habitual γ ≫ γA , γB queda simplemente como pB − pA = γh


Agustin Martin Domingo A
A

pA

h
pA

M
dA

dB
M

N
d
N

pB
B
h

Figura 2–4 Medida de la presión con un manómetro. a) Manómetro normal b) Manómetro diferencial.

2.3. Variación de la presión con la altura en un gas perfecto en reposo.

Hemos visto que la variación diferencial de la presión con la profundidad viene dada por dp = ρgdh = −ρgdz =
−γdz y hemos resuelto el problema para un fluido incompresible (γ y ρ no varían con la presión).
Consideremos ahora el caso de un fluido compresible, y nos restringiremos al caso en que este fluido compresible es
un gas que se comporta como un gas ideal (ésto ocurre con el aire atmosférico y con buena parte de sus componentes
para un rango de presión y temperatura relativamente alto).
Obtengamos en primer lugar la dependencia de la densidad del gas ideal con presión y temperatura. Si tenemos en
¯ ¯
cuenta que ρ = m/V y que n = m/M , con M el peso molecular promedio del gas ideal,
m m
pV = nRT = p = ¯
ρ M
por lo que la densidad queda:
¯
Mp
ρ(p, T ) = . (2–8)
RT
La variación infinitesimal de la presión con la altura se puede entonces escribir como
¯
Mp
dp = −ρgdz = −ρ(p, T )gdz = − gdz (2–9)
RT
Agrupando variables,
dp ¯
M g dz
=−
p R T

2.3. Variación de la presión con la altura en un gas perfecto en reposo. 17

que, al integrar queda

ln
p
p0
ln

=−
p
p0

¯
=−
¯
Mg
R

M g z − z0
R T0
que da una variación exponencial de la presión con la altura

p = p0 exp −

que es una buena aproximación para pequeñas alturas.

2.3.2. Dependencia lineal
¯
Mg
RT0
0
z
dz
T

(z − z0 )
.

(2–11)
(2–10)

Ahora bien, en general T = T (z), es una función de z y por tanto, para poder obtener la variación de la presión con la
elevación es necesario conocer la función T (z).

2.3.1. Atmósfera isoterma

Consideremos primero el caso más sencillo, el de una atmósfera isoterma con T = T0 = cte. En este caso, la
expresión anterior queda

(2–12)

Para grandes alturas, la temperatura de la atmósfera disminuye de forma aproximadamente lineal hasta unos 11000 m,
expresándose esta dependencia en la forma
T (z) = T0 − βter z (2–13)
con T0 la temperatura absoluta al nivel del mar y βter el denominado gradiente térmico. Estas magnitudes varían de
un día a otro, aunque como ejemplo, unos valores típicos serían T0 = 288, 16 K∼ 15◦ C y βter ∼ 6, 50 K km−1 . La
integral queda entonces
¯ z ¯ z ¯
p Mg d M g − ln(T0 − βter z) Mg T0 − βter z
ln =− =− = ln (2–14)
p0 R 0 T0 − βter z R βter 0 Rβter T0

que se puede escribir como
¯
Mg
p T0 − βter z Rβter
ln = ln ,
p0 T0
quedando la dependencia de la presión con la altura, para una dependencia lineal de la temperatura con la altura, en la
forma
¯
Mg
T0 − βter z Rβter
p = p0 , (2–15)
T0

2.3.3. Atmósfera tipo

De entre los distintos casos particulares considerados, la atmósfera de densidad homogénea tiene un interés esencial-
mente teórico, mientras que la aproximación de atmósfera isoterma es aplicable únicamente a capas de la atmósfera
con densidad constante. Esto último es aplicable a la baja estratosfera, especialmente a altas alturas, pero, salvo que se
trate de capas delgadas, es difícilmente aplicable a la troposfera. La aproximación de atmósfera adiabática tiene una
aplicación muy interesante, ya que da un límite superior para el gradiente de una atmósfera verticalmente estable y
tiene la distribución de temperatura de una capa mezclada verticalmente.
La atmósfera real tiene, en general, valores inferiores del gradiente térmico, y aunque en general éste no es constante
con la altura, es posible deﬁnir una atmósfera de gradiente constante (lineal) que se aproxime al caso real.


h = ξ sen θ
hCG = ξCG sen θ
hCP = ξCP sen θ xcp
hcp

CP
hcg

O
CG

x

ycp
h

y

y
ξ
θ

ξcg

x

dA = dxdy
ξcp

Figura 2–5 Cálculo de la fuerza sobre una superficie plana. Consideramos una superficie plana que forma un ángulo θ con la superficie libre del
líquido.

Se han definido distintas atmósferas tipo que se utilizan como referencia para una aproximación a la atmósfera real
y para calibrar y utilizar altímetros. Las condiciones en que se define la atmósfera tipo son las siguientes,

• Atmósfera de aire seco puro con composición química constante en sentido vertical y peso molecular promedio
¯
M = 28, 9644 g.

• El comportamiento del gas puede aproximarse por el de un gas ideal. Ésto es razonable para p ≈1 atm.

• Se utiliza como valor estándar de la aceleración el de la gravedad a nivel del mar, g ≈ 9, 80 ms−2

• El sistema se encuentra en equilibrio hidrostático.

• Se toma como referencia la temperatura T0 = 15◦ C a nivel del mar y la presión p0 = 1 atm.

• Para alturas hasta los 11000 m (tropopausa), se considera constante el gradiente térmico y de valor βter =
−dT /dz = 6, 5 K·km−1

• Para alturas mayores de 11000 m (estratosfera) y hasta 20000 m se considera la temperatura constante e igual
a −56, 5◦C. A partir de esa altura se supone un gradiente térmico constante de valor −1, 0 K·km−1 hasta los
32000 m.

2.4. Fuerzas hidrostáticas sobre superficies.

2.4.1. Fuerza sobre superficies planas.

Consideremos una placa plana de forma arbitraria que se encuentra sumergida completamente en el líquido, como
se muestra en la figura.

2.4. Fuerzas hidrostáticas sobre superficies. 19

Valor de la fuerza resultante.

Tomamos como sistema de referencia el sistema de coordenadas [x, y] centrado en el centro de gravedad de la
superficie. La fuerza total que actúa sobre la placa se expresa como

que es simplemente
FA =
A
pdA =
A
(patm + ρgh)dA = patm A + ρg

A

Ahora bien, si recordamos la definición del centro de gravedad Aξcg =

FA = pcg A
A
A

A
hdA

Si tenemos en cuenta que h = ξ sen θ, y que θ se mantiene constante sobre la placa, la fuerza queda como

FA = Apatm + ρg ξ sen θdA = Apatm + ρg sen θ ξdA

ξdA, la expresión anterior se convierte en

FA = Apatm + ρg sen θAξcg = Apatm + ρgAhcg = (patm + ρghcg )A

La fuerza total que actúa sobre una superficie plana cualquiera sumergida en un fluido uniforme es
igual a la presión que hay en el centro de gravedad de dicha superficie multiplicada por su área, con
independencia de la forma de la superficie plana o de su inclinación θ.

Punto de aplicación de la fuerza resultante.

Imaginemos un recipiente que contiene un líquido y sobre éste una superficie A libre que encaja perfectamente en
la pared del recipiente. Al poder moverse la superficie, ésta será empujada por el líquido y se saldrá. Para evitarlo
(2–16)

habría que aplicar sobre la superficie una fuerza normal a la misma de magnitud la fuerza que ejerce el líquido sobre
la superficie. Para que además la superficie no gire, esta fuerza debe aplicarse en un punto determinado de forma que
el momento total del sistema de fuerzas ejercido por el líquido sobre la superficie se compense con el momento de la
fuerza equivalente aplicada en ese punto. Este punto es el centro de presiones.
Cuanto mayor es la profundidad, mayor es la presión. Por tanto, el punto de actuación de la fuerza total resultante
(centro de presiones) no coincide con el centro de gravedad, sino que debe encontrarse más abajo. La línea de acción
de esta fuerza debe pasar precisamente por este centro de presiones. Para calcular la posición (xcp , ycp ) del centro de
presiones se suman los momentos de las fuerzas elementales pdA respecto del centro de gravedad y se igualan con
el momento (respecto del centro de gravedad) de la fuerza resultante aplicada en el centro de presiones. Obtengamos
ésto para cada una de las componentes,

Eje x: La componente del momento en la dirección x viene dada por

F · ycp = ypdA (2–17)
A

Como p = patm + ρgξ sen θ se tiene a su vez que

F · ycp = ypdA = y(patm + ρgξ sen θ)dA = patm ydA + ρg sen θ yξdA
A A A

El primer sumando es cero por la definición del centro de gravedad, y como y = ξcg − ξ se tiene, con ξ = ξCG − y

F · ycp = ρg sen θ ξcg ydA − y 2 dA
A A

De nuevo el primer sumando es cero por la definición de centro de gravedad de una superficie, por lo que la componente
x del momento es
F · ycp = −ρ sen θIxx (2–18)


con

Agustin Martin Domingo ycp = −
Ixx =

ρg sen θ
F
A

Ixx = −
y 2 dA

el denominado momento de inercia de la sección plana (del área de la placa) respecto de su eje horizontal x, calculado

ρg sen θ
pcg A
Ixx

El signo − indica que el centro de presiones está por debajo del centro de gravedad, y como se puede ver, esta posición
depende de la inclinación θ.
Al aumentar la profundidad a la que se encuentra la placa, ycp se acerca al centro de gravedad, ya que todos los
factores que intervienen en 2–20 permanecen constantes excepto pcg que aumenta.

Eje y: La componente del momento en la dirección y viene dada por

que, al igual que antes se convierte en

F · xcp =
A
xpdA =
A
F · xcp =
A
xpdA

x(patm + ρgξ sen θ)dA = patm
A
xdA + ρg sen θ xξdA.

El primer sumando también es cero por la deﬁnición de centro de gravedad y teniendo de nuevo en cuenta que y =
ξcg − ξ se tiene, al sustituir ξ = ξcg − y
(2–19)

en el plano de la placa. Como F = pcg A, la distancia en el plano de la placa a la que se encuentra el centro de presiones
respecto del centro de gravedad, viene dada por

(2–20)

(2–21)

F · xcp = ρg sen θ ξcg xdA − xydA ,
A A

lo que da, para la componente x de la posición del centro de presiones

F · xcp = −ρ sen θIxy (2–22)

con
Ixy = xydA (2–23)
A
conocido como el producto de inercia. Dividiendo por la fuerza F = pcg , la componente x de la posición del centro
de presiones respecto del centro de gravedad en el sistema de coordenadas xy situado sobre la placa es:
ρg sen θ ρg sen θ
xcp = − Ixy = − Ixy (2–24)
F pcg A
El producto de inercia puede tener cualquier signo, al contrario que el momento de inercia, que puede ser sólo positivo.
Dependiendo de este signo el centro de presiones se encontrará a un lado u otro del centro de gravedad. Si Ixy es < 0,
xcp < 0, mientras que si Ixy es negativo, xcp > 0. En el caso particular en el que Ixy = 0, se tiene que xcp = 0 y, por
tanto, el centro de presiones está directamente debajo del centro de gravedad sobre el eje y en el plano de la placa.
Como normalmente la presión atmosférica pa no se tiene en cuenta ya que actúa sobre ambos lados de la placa (cara
interna de un barco o cara seca de una compuerta o de una presa), se puede escribir pcg = ρghcg y la fuerza neta como
F = ρghcg A. Las expresiones 2–20 y 2–24 se pueden escribir como
Ixx sen θ
ycp = − (2–25)
hcg A
Ixy sen θ
xcp =− , (2–26)
hcg A
válidas únicamente cuando la superﬁcie superior del líquido está a la presión atmosférica.

2.5. Equilibrio de un cuerpo sumergido. 21

Empuje

2.4.2. Fuerza sobre superficies curvas.
Empuje

Figura 2–6 Flotación y principio de Arquímedes.
ρl´ gVdes

mg
ßq

Para calcular la fuerza ejercida por el agua sobre una superficie curva hay que tener en cuenta la combinación de dos
componentes, una horizontal y otra vertical. La componente horizontal se calcula obteniendo la fuerza que actuaría
sobre la proyección de la superficie curva en el plano vertical, con su valor y su línea de aplicación a través del
correspondiente centro de presiones.
La componente vertical se obtiene directamente a partir del peso del agua sobre la superficie curva, o, en su caso del
empuje vertical del agua sobre la misma, actuando a través del dentro de gravedad del líquido o del líquido desalojado,
dependiendo del caso.

2.5. Equilibrio de un cuerpo sumergido.

2.5.1. Flotación. Principio de Arquímedes.

Consideremos un cuerpo parcial o totalmente sumergido en un fluido en reposo (líquido o gas). Este fluido ejerce
presiones sobre todas las partes de la superficie del cuerpo, mayores cuanto mayor es la profundidad. Se denomina
empuje sobre el cuerpo sumergido a la fuerza total hacia arriba ejercida por el agua.
Calculemos el valor de esta fuerza. Si se considera la superficie cerrada que delimita el cuerpo sumergido, la presión
en cada punto de esa superficie dará lugar a una distribución de fuerzas cuya resultante es precisamente el empuje
que estamos calculando. Consideremos ahora que en el recipiente hubiera el mismo fluido y hasta el mismo nivel
que cuando estaba el cuerpo y dentro de este fluido una superficie cerrada imaginaria que coincida con la superficie
exterior del cuerpo que está sumergido en el caso real. Como la presión depende únicamente de la profundidad, la
presión en todos los puntos de esta superficie imaginaria es la misma que había en los puntos correspondientes de la
superficie real, por lo que la fuerza resultante que ejerce el fluido de fuera de la superficie sobre el fluido de dentro
de la misma será la misma que ejerce el fluido sobre el cuerpo. Ahora bien, ahora es fácil hacer el cálculo, ya que el
fluido se encuentra en equilibrio mecánico y por tanto la fuerza que compensa el empuje es el peso del propio fluido
dentro de la superficie imaginaria.
¿Dónde estará aplicada esa fuerza? Volvamos de nuevo al caso del fluido con la superficie cerrada imaginaria. Si el
sistema está en equilibrio mecánico no habrá rotación y para ello es necesario que los puntos de aplicación del empuje
y del peso del fluido desalojado sean el mismo. Por tanto, se tiene el principio de Arquímedes,

Un cuerpo total o parcialmente sumergido en un fluido es empujado hacia arriba con una fuerza igual
al peso del fluido desalojado, que actúa verticalmente a través del centro de gravedad del fluido antes de
ser desplazado.


Agustin Martin Domingo V
ρg2V

mg
V

Figura 2–7 Contraejemplo de porqué el centro de presiones está a lo largo de la línea que pasa por el centro de gravedad del fluido desalojado.

Frecuentemente se tiene la idea de que este empuje actúa a través del centro de gravedad del cuerpo sumergido. Esto
es en general incorrecto y para verlo más gráficamente consideremos un contraejemplo. Sea un cuerpo con forma de
pesa de gimnasia con la particularidad de que cada bola está hecha de un material distinto, una de hierro, la otra de
corcho blanco (poliestireno expandido), ambas del mismo volumen, con una masa y volumen de la barra prácticamente
despreciables. El centro de gravedad de la pesa está a efectos prácticos en el centro de la bola pesada. ¿Qué ocurriría
si el centro de aplicación del empuje estuviera también en el centro de gravedad del cuerpo? Ambas fuerzas estarían
aplicadas sobre el mismo punto por lo que no se produciría rotación de la pesa hasta que la bola pesada quede abajo y
la ligera arriba, contra lo que muestra la experiencia. Por el contrario, en el caso real en el que el empuje actúa sobre
el centro de gravedad del fluido desalojado, este par de fuerzas se produce en el sentido adecuado (como debe ser).

2.5.2. Estabilidad de un cuerpo semisumergido

Al estudiar un cuerpo semisumergido no solamente es importante conocer el valor del empuje, sino cómo actúa
éste. Así, es de gran importancia conocer la capacidad del cuerpo flotante de recuperar el estado de equilibrio una vez
sacado del mismo. No nos sirve sólo que el barco flote, sino la posición en que debe hacerlo.
La estabilidad longitudinal (a lo largo) de los barcos es normalmente muy considerable, por lo que nos centraremos
en la estabilidad transversal (frente a giros a izquierza o derecha). Consideremos el barco cuya sección transversal
se esquematiza en la figura 2–8 y denominaremos eje de flotación OO′ al eje de simetría del barco cuando éste se
encuentra en la posición de equilibrio. En ésta, tanto el centro de gravedad como el de empuje están sobre este eje.
Cuando como consecuencia de las fuerzas que actúan sobre el barco, éste se inclina un cierto ángulo α, una parte del
barco ha salido del agua y otra parte que antes estaba fuera ha entrado. La posición del centro de gravedad del barco
no ha variado en este giro, pero el centro de gravedad del líquido desplazado sí, quedando más hacia la izquierda. Así,
el eje de empuje se desplaza, en este caso hacia la derecha. Si este eje de empuje queda a la derecha del centro de
gravedad se produce un par recuperador que tiende a devolver el barco a su posición original, mientras que si queda
a la izquierda el par fomentará el vuelco del barco. Denominamos metacentro M al punto de cruce del eje de empuje
y el eje de flotación del barco, y altura metacéntrica h a la distancia entre el metacentro y el centro de gravedad a lo
largo del eje de flotación.
La altura metacéntrica es una magnitud característica de la sección transversal del cuerpo para un peso dado y su
valor un indicador de la estabilidad del cuerpo. Valores típicos para barcos comerciales están entre 0, 3 y 0, 8 m. El
cálculo de la altura metacéntrica es, en general, un trabajo pesado.
Dependiendo de la posición relativa del metacentro y el centro de gravedad se tienen 3 casos:

Equilibrio estable
El metacentro se encuentra por encima del centro de gravedad. De esta forma el par de fuerzas es un par
restaurador que lleva al barco a su posición inicial.

2.5. Equilibrio de un cuerpo sumergido. 23

Equilibrio indiferente
G

E
G

Figura 2–8 Estabilidad de un cuerpo semisumergido.

El metacentro y el centro de gravedad coinciden. No hay par de fuerzas.
Equilibrio inestable
M

E

El metacentro se encuentra por debajo del centro de gravedad. El par de fuerzas es un par de vuelco.

Así, la estabilidad del barco será tanto mayor cuanto más bajo se encuentre el centro de gravedad y cuanto mayor sea
la altura metacéntrica.
En condiciones de equilibrio estable el barco oscilará con una cierta frecuencia. El problema es similar a un péndulo
con el metacentro como punto ﬁjo. Así, si MM es el momento de inercia respecto del eje perpendicular al plano de
oscilación que pasa por el metacentro, el periodo será

MM
T = 2π (2–27)
mgh

donde h es la altura metacéntrica. Al disminuir el producto del peso mg por la altura metacéntrica h el periodo aumenta
y por tanto el balanceo es algo menos desagradable. Sin embargo en la práctica ésto implica una disminución de la
estabilidad.

Capítulo 3

Dinámica de fluidos perfectos.

3.1.

3.2.

3.3.

3.4.
Movimiento de un fluido perfecto. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.1.1. Aproximaciones clásicas al estudio de los fluidos. . . . . . . . . . . . . . . . . . . . . . .
3.1.2. Tipos de flujo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.1.3. Senda, línea de corriente y de traza. Tubo de flujo. . . . . . . . . . . . . . . . . . . . . . .
Ecuación de continuidad. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.2.1. Forma integral. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.2.2. Forma diferencial. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Fuerza y aceleración en un elemento de fluido. . . . . . . . . . . . . . . . . . . . . . . . . . .
3.3.1. Fuerza neta sobre un elemento de un tubo de flujo. Variación de la cantidad de movimiento.
3.3.2. Ecuaciones de movimiento microscópicas. Aceleración. Ecuación de Euler. . . . . . . . .
3.3.3. Aplicación a un fluido en reposo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.3.4. Aplicación a un fluido en reposo en un sistema de referencia no inercial . . . . . . . . . .
3.3.5. Vorticidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Conservación de la energía. Ecuación de Bernoulli. . . . . . . . . . . . . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
25
25
26
26
28
28
29
31
31
32
33
33
36
37
3.5. Aplicaciones de las ecuaciones anteriores. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
3.5.1. Velocidad de salida por un orificio (teorema de Torricelli). . . . . . . . . . . . . . . . . . . 40
3.5.2. Tiempo de vaciado de un depósito. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.5.3. Presión dinámica y presión estática. Tubo de Pitot . . . . . . . . . . . . . . . . . . . . . . . 41
3.5.4. Efecto Venturi. Venturímetro. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42

3.1. Movimiento de un fluido perfecto.

3.1.1. Aproximaciones clásicas al estudio de los fluidos.

A la hora de abordar el estudio de los fluidos, surgen dos aproximaciones clásicas, la de Euler y la de Lagrange.

Método de Lagrange. El método de Lagrange constituye una generalización directa de la mecánica del punto mate-
rial. Se estudia un volumen pequeño del fluido y se sigue el movimiento de cada una de las partículas (de coordenadas
x, y, z) en función del tiempo t, a través de la ecuación que describe la trayectoria de cada una de las partículas. El prin-
cipal inconveniente de este sistema es que hacen falta una gran cantidad de ecuaciones para describir el movimiento
del sistema, por lo que en la práctica no es útil.

Método de Euler. En la aproximación de Euler se desiste de describir el movimiento del fluido mediante la historia
de cada una de las partículas individuales. En su lugar se especifica el movimiento del fluido por la densidad ρ(x, y, z)
y la velocidad v(r, t) de las partículas del mismo en ese punto, como una función del tiempo y del espacio. En otras
palabras, en este método se estudia un punto del espacio y como es el movimiento del fluido en ese punto en función
del tiempo. Las herramientas de trabajo serán las típicas de la Teoría de Campos, con un campo de presiones, un
campo de velocidades y un campo de densidades.

25

26 Capítulo 3. Dinámica de fluidos perfectos.

3.1.2. Tipos de flujo.

Atendiendo a las características del flujo, éste puede clasificarse de acuerdo con distintos criterios.

Flujo estacionario/no estacionario. Se dice que el flujo es estacionario si la velocidad v(r) y la densidad ρ(r) del
flujo en un punto no dependen del tiempo y no estacionario en caso contrario. Esto no quiere decir que la velocidad
y la densidad deban ser las mismas en dos puntos distintos, sino sólo que en un mismo punto no deben variar con el
tiempo.

Flujo irrotacional/rotacional. Se dice que el flujo es irrotacional cuando el elemento del fluido en un punto dado
no tiene una velocidad angular neta alrededor de dicho punto y que es rotacional en caso contrario. Un fluido que
circula a través de una tubería recta de sección uniforme sería un ejemplo simple de flujo irrotacional, mientras que un
remolino de un río sería un ejemplo de flujo rotacional.

Flujo compresible/incompresible. Se dice que el flujo es compresible si la densidad ρ en el mismo varía, como
por ejemplo ocurre en los gases en el caso más general, mientras que se dice que el flujo es incompresible cuando
la densidad ρ apenas varía como es el caso de los líquidos. Nótese que es posible tener un flujo aproximadamente
incompresible aunque el fluido en movimiento en sí sea un fluido compresible siempre que a lo largo del flujo en la
región considerada la densidad ρ sea prácticamente la misma en todos los puntos.

Flujo viscoso/no viscoso. Se dice que el flujo es viscoso cuando aparecen en él importantes fuerzas de rozamiento
que no se pueden despreciar. Como consecuencia de estas fuerzas de rozamiento aparecen unas fuerzas tangenciales
entre las capas del fluido en movimiento relativo y hay una disipación de energía mecánica. Por el contrario se dice
que el flujo es no viscoso cuando estas fuerzas de rozamiento son muy pequeñas o bien no se tienen en cuenta.

3.1.3. Senda, línea de corriente y de traza. Tubo de flujo.

Vamos a definir ahora una serie de líneas que nos permitirán describir el movimiento de un fluido. Estas son la
senda, la traza o línea de traza y la línea de corriente o línea de flujo. En general, cada una de éstas líneas es distinta
de las otras, pero en el caso de un flujo estacionario las tres coinciden.

Senda.

Se denomina al camino seguido realmente por una partícula de fluido. La senda se define para una partícula a
lo largo del tiempo. Su carácter es fundamentalmente experimental y se obtendría experimentalmente soltando una
partícula marcador y haciendo una fotografía a obturador abierto durante el tiempo del estudio.
Para obtener analíticamente la senda integraríamos el campo de velocidades para obtener las ecuaciones paramétri-
cas (x(t), y(t), z(t)) de la senda para un elemento genérico del flujo que en un cierto instante de referencia pasa por
un punto de referencia dado (lo que nos permite obtener las constantes de integración). Si eliminamos el tiempo entre
las distintas ecuaciones obtenemos la ecuación e la senda.

Línea de traza.

Se denomina línea de traza al lugar geométrico de las partículas que en instantes sucesivos pasaron por un punto
dado. Su carácter es también fundamentalmente experimental y en la práctica se obtendría por ejemplo inyectando de
forma continua en un punto fijo del flujo una serie de partículas marcadas y tomando una fotografía.
Para obtener analíticamente la ecuación de la traza debemos previamente obtener la familia de constantes de inte-
gración correspondientes a las sendas de las partículas que en instantes tt anteriores al instante genérico t considerado

3.1. Movimiento de un fluido perfecto. 27
y

3,5

2,5
4

3

2

1,5

1
1
Líneas de campo

Senda

1,5
Traza
2 2,5
Senda(s)

3
Traza

Líneas de campo

3,5
x

Figura 3–1 Representación espacial (x, y) de un campo de velocidades y de algunas líneas de flujo junto con varias sendas y una línea de traza. El
campo de velocidades representado es el campo (x/(1 + t), y, 0) en el instante t = 1 s. La senda representada en trazo más grueso corresponde a
la seguida por un elemento del fluido que en t = 1 s pasa por el punto (2,7, 2,7) y la traza al lugar geométrico de los elementos del fluido que en
instantes sucesivos pasan por dicho punto. Junto a la senda anterior se muestran las sendas de las partículas que pasarán por el citado punto en los
instantes t = 1,5 s, t = 2 s y t = 3 s, representadas con líneas finas de trazos.

pasaron por el punto (x0 , y0 , z0 ) de referencia. Sustituyendo estos valores en las ecuaciones paramétricas de la senda
tendríamos las sendas de cada una de estas partículas y eliminando tt tendríamos la ecuación de la traza, en el instante
t, de las partículas que en instantes previos pasaron por (x0 , y0 , z0 ).

Línea de corriente o línea de flujo.

Se define la línea de corriente como una línea que en un instante dado es tangente al vector velocidad en cada punto.
Es importante recalcar que la línea de corriente está definida para un instante dado, mientras que la senda y la línea de
traza contienen información de otros instantes.
Al contrario de las dos anteriores, de carácter básicamente experimental, la línea de corriente tiene un profundo
sustrato matemático. Para su obtención experimental se extenderían en el flujo un conjunto de partículas marcadoras
y se tomaría una fotografía del conjunto “casi” instantánea (obturador abierto durante un breve periodo de tiempo),
de modo que en ésta se aprecian para cada partícula líneas pequeñas que marcan la dirección (ángulo de la línea) y el
módulo de la velocidad (relacionado con la longitud de la línea).
Como la velocidad debe ser tangente en cada punto a la línea de flujo, debe cumplirse v × dr = 0 con dr =
(dx, dy, dz). Así, éstas deben cumplir
dx dy dz dr
= = = (3–1)
vx vy vz v


compleja.

Tubo de flujo.
1

A1
vdt
2

A2

Figura 3–2 Tubo de corriente utilizado para la obtención de la ecuación de continuidad en forma integral.

que, si vx , vy vz son conocidas en función de la posición y el tiempo, pueden integrarse para dar la ecuación de la
línea de corriente que en un determinado instante pasa por un punto dado. Nótese que esta integración puede ser muy

Un tubo de flujo o de corriente está formado por un haz de líneas de flujo. Como las líneas de flujo son tangentes
al vector velocidad en cada punto, las líneas de flujo no atraviesan las paredes del tubo de flujo. Esto es así porque
las paredes están a su vez formadas por líneas de flujo y si las atravesaran dos líneas de flujo se cortarían, estando la
velocidad indeterminada en el punto de corte. Por tanto, a pesar de estar limitado por la superficie ficticia que envuelve
el haz de líneas de flujo, éste se comporta a todos los efectos como una superficie impermeable.

3.2. Ecuación de continuidad.

3.2.1. Forma integral.

Consideremos un tubo de corriente estrecho, de forma que se pueda considerar uniforme la velocidad en cualquier
sección del tubo perpendicular al flujo. En el interior del tubo la velocidad del flujo es paralela a la línea de corriente
en cada punto, pudiendo ser estas velocidades distintas en cada punto.
Sea v1 la velocidad de la partícula en el punto 1, y v2 la velocidad de la partícula en el punto 2, con A1 y A2
las secciones transversales de los tubos, perpendiculares a las líneas de corriente. Si el tubo es estrecho v1 y v2 son
uniformes en A1 y A2 respectivamente.
En un intervalo de tiempo dt, un elemento del fluido recorrerá una distancia v dt, por lo que en el tiempo dt pasará
por A1 la masa de fluido
dm1 = ρ1 · A1 · v1 dt
donde ρ1 es la densidad del fluido al pasar por la sección 1. El flujo de masa o caudal másico se define como la masa
que atraviesa una sección en la unidad de tiempo, y viene dado por
dm1
Qm = = ρ1 A1 v1 , (3–2)
dt
donde se considera implícitamente que en ese intervalo infinitesimal de tiempo ni A ni v varían apreciablemente en el
recorrido del fluido v dt. El caudal másico a través de la sección A1 es ρ1 A1 v1 y a través de la sección A2 es ρ2 A2 v2 .

3.2. Ecuación de continuidad. 29

Agustin Martin Domingo Figura 3–3

⇒
S2

Como las partículas del flujo no pueden atravesar las paredes del tubo de flujo debe cumplirse que, si el régimen es
permanente (o estacionario) y no hay fuentes ni sumideros de partículas, ambos caudales másicos han de ser iguales

Qm = ρ1 A1 v1 = ρ2 A2 v2 Qm = ρAv = cte

y análogamente para cualquier otra sección A perpendicular al tubo de flujo, por lo que esta ley de conservación de la
masa o ecuación de continuidad se puede escribir simplemente como

ρAv = cte

a través de cualquier sección del tubo de flujo perpendicular al mismo en régimen estacionario.
(3–3)

Para el caso particular de flujo incompresible ρ no depende del punto y esta ecuación de continuidad puede escribirse
como
A1 v1 = A2 v2 ⇒ Q = Av = cte (3–4)
donde Q es el caudal o volumen que atraviesa la sección en la unidad de tiempo.
Por ejemplo, en una canalización por la que circula un fluido incompresible, se tiene la sencilla relación S1 v1 =
S2 v2 , que da, para la relación entre velocidades
S1
v2 = v1 (3–5)
S2
En todo el cálculo anterior hemos considerado implícitamente que la velocidad v es uniforme en cada sección. Esto
no es cierto en el caso general, pero la ecuación de continuidad sigue siendo válida en las mismas condiciones si la
densidad es uniforme en la sección y en vez de la velocidad se utiliza la velocidad promedio en la sección,
1
v=
¯ v dA (3–6)
A A

3.2.2. Forma diferencial.

Estudiaremos la conservación de la masa en un prisma infinitesimal de aristas dx, dy y dz fijo en el espacio. El
caudal másico que atraviesa la cara perpendicular al eje y de la izquierda es

dQmy = ρvy dAy = ρvy dxdz

mientras que el caudal másico que atraviesa la cara perpendicular al eje y de la derecha (a y + dy) es
∂ ∂
dQmy+dy = ρvy dAy + [ρvy dAy ]dy = [ρvy + (ρvy )dy]dxdz
∂y ∂y
De este modo, el caudal másico neto a través de las caras perpendiculares al eje y es
∂
Qneto = Qmy − Qmy+dy = −
my (ρvy )dxdydz (3–7a)
∂y



y
z
z

dy
dx
dz

y

Figura 3–4 Obtención de la ecuación de continuidad en forma diferencial.
x

Una expresión similar se obtiene para los caudales másicos netos a través de las caras perpendiculares al eje x,

∂
Qneto = Qmx − Qmx+dx = −
mx (ρvx )dxdydz, (3–7b)
∂x
y otra para los caudales másicos netos a través de las caras perpendiculares al eje z,

∂
Qneto = Qmz − Qmz+dz = −
mz (ρvz )dxdydz. (3–7c)
∂z
Si tenemos en cuenta que la masa total contenida en el prisma es dm = ρdV = ρdxdydz, el caudal másico total
también debe ser
∂
dQneto = (ρdV )
total
∂t
que representa la velocidad con que varía la masa en el punto ﬁjo del espacio en el que estamos considerando el
elemento de volumen.† Este caudal másico debe ser a su vez la suma de los caudales másicos a través de las caras
perpendiculares a cada una de las direcciones x, y y z.

∂ ∂ρ ∂ ∂ ∂
dQneto =
total (ρdV ) = dxdydz = − (ρvx ) − (ρvy ) − (ρvz ) dxdydz. (3–8)
∂t ∂t ∂x ∂y ∂z

La parte entre corchetes no es más que menos la divergencia de ρv por lo que la expresión anterior se puede escribir
como
∂ρ
= −∇(ρv) (3–9)
∂t
o
∂ρ
+ ∇(ρv) = 0. (3–10)
∂t
† Nótese que estamos estudiando un elemento del espacio ﬁjo en el mismo y por ello aparecen derivadas parciales. Véase la discusión sobre la

velocidad en torno a la ecuación (3–23).

3.3. Fuerza y aceleración en un elemento de fluido. 31

Agustin Martin Domingo dm1
S1
11
′

v1 dt
v1
2

Figura 3–5 Variación de la cantidad de movimiento.
2′

v2 dt
dm2
S2
v2

que es la expresión de la ecuación de continuidad en forma diferencial. Para el caso particular de un fluido incompre-
sible, la densidad ρ no depende de la posición ni del tiempo y la ecuación de continuidad queda como
∇v = 0. (3–11)

En régimen estacionario, la densidad en un punto no cambia con el tiempo, por lo que ∂ρ/∂t = 0 y la ecuación de
continuidad queda ∇(ρv) = 0.

3.3. Fuerza y aceleración en un elemento de fluido.

3.3.1. Fuerza neta sobre un elemento de un tubo de flujo. Variación de la cantidad de movimiento.

Calculemos la fuerza total que actúa sobre un elemento del fluido que circula por un tubo de flujo estrecho en
régimen estacionario. Para ello, calcularemos la variación de la cantidad de movimiento de dicho elemento por unidad
de tiempo, lo que nos dará la fuerza neta que actúa sobre el mismo.
Consideremos ahora el elemento de fluido que se muestra en la figura 3–5 y que inicialmente se encuentra entre las
secciones 1 y 2. Un intervalo de tiempo dt después, el fluido habrá avanzado hasta las secciones 1′ y 2′ . La cantidad
de movimiento inicial del elemento será la suma de las cantidades de movimiento de los subelementos 1 − 1′ y 1′ − 2,
(mv )inicial = dm1 v1 + (mv )1′ −2
mientras que la cantidad de movimiento del elemento pasado un intervalo de tiempo dt será la suma de las cantidades
de movimiento de los subelementos 1′ − 2 y 2 − 2′ ,
(mv )f inal = (mv )1′ −2 + dm2 v2 ,
siendo la cantidad de movimiento de la parte 1′ −2 la misma en ambos casos al ser el régimen permanente. La variación
de la cantidad de movimiento en dicho intervalo es entonces
d(mv ) = (mv )f inal − (mv )inicial = dm2 v2 − dm1 v1 .
La ecuación de continuidad (3–3) nos dice que las masas dm1 y dm2 deben ser iguales, por lo que la variación de
cantidad de movimiento por unidad de tiempo puede escribirse como
d dm
(mv ) = Ftot = (v2 − v1 ) = Qm (v2 − v1 ) (3–12)
dt dt
donde Qm es el caudal másico (3–2), y esta variación de cantidad de movimiento por unidad de tiempo debe ser igual
a la fuerza neta que actúa sobre el elemento.



y
−j

i
z

dy

dy
k

−k

3.3.2. Ecuaciones de movimiento microscópicas. Aceleración. Ecuación de Euler.
dx

dx
−i
dz

y

x
j

Figura 3–6 Volumen de control utilizado para la obtención de la ecuación de Euler.

Seguimos considerando el caso de un fluido perfecto (no viscoso). Para resolver el problema del movimiento de
un fluido perfecto hay que escribir las ecuaciones de movimiento que relacionan las velocidades con las fuerzas que
actúan en el fluido.
Recordemos que, para un fluido perfecto, las fuerzas en el seno del fluido en movimiento son siempre perpendicu-
lares a la superficie, son las fuerzas de presión. El fluido que se encuentra del lado de la superficie a donde apunta n
ejerce una fuerza F = −p dS n sobre la superficie dS.
Consideremos el elemento de volumen diferencial de la figura,† de dimensiones dx dy dz y caras perpendiculares a
los ejes coordenados. En la dirección y actúa, sobre la cara que se encuentra a una distancia y del origen, una fuerza

dF (y) = p(y) dx dz (3–13)

mientras que sobre la cara que se encuentra a una coordenada y + dy actúa una fuerza

∂p
dF (y + dy) = −p(y + dy) dx dz = − p(y) + dy dx dz. (3–14)
∂y

La componente y de la fuerza de presión neta será la suma de las dos fuerzas anteriores,

∂p ∂p
dFy = −p(y) dx dz − dy dx dz + p(y) dx dz = − dy dx dz (3–15a)
∂y ∂y

obteniéndose de la misma manera las expresiones para las componentes netas de la fuerza en las direcciones x y z,

∂p
dFx = − dy dx dz (3–15b)
∂y
∂p
dFz = − dy dx dz (3–15c)
∂y
† Ahora, este elemento de volumen es un elemento de volumen del fluido que se desplaza con el tiempo y por tanto debemos utilizar una

derivada total.


La suma vectorial de estas componentes da la fuerza de presión total que actúa sobre el elemento de volumen, de masa
ρ dV = ρ dx dy dz.
dFpres. = −
∂p
∂x
i+
∂p
∂y
j+
∂p
∂z

la ecuación de movimiento para los fluidos perfectos queda

dFtot = dFpres. + fext dm =

que se puede escribir como
dv
dt
dv
dt

1
ρ
k dx dy dz = −∇p dx dy dz.

Si recordamos la segunda de las leyes del movimiento de Newton,

F = ma = m
dv
dt

dm = −∇p dx dy dz + fext ρ dx dy dz

= − ∇p + fext

La ecuación anterior no es más que la ecuación de movimiento F = ma en la que, junto a la contribución de la
(3–16)

(3–17)

donde fext representa la aceleración debida a una fuerza externa que actúa sobre el elemento de fluido. Normalmente,
fext = g = −guz .

fuerza de masa a través de fext aparece un término debido a las diferencias de presión entre dos puntos. Como ∇p está
dirigido hacia el sentido de presiones crecientes y el movimiento tiene lugar empujado por las altas presiones hacia
las bajas presiones, aparece el signo negativo delante del gradiente.

3.3.3. Aplicación a un fluido en reposo

Si además de que la única aceleración de masa a que está sometido el fluido es la aceleración gravitatoria, éste está
dv
en reposo, v es cero en todos los puntos del fluido así que = 0 y la ecuación anterior queda
dt

−∇p + ρg = 0 (3–18)

de componentes

∂p
− =0
∂x
∂p
− =0
∂y
∂p
− − ρg = 0.
∂z
De las dos primeras expresiones de obtiene directamente que la presión en este caso no puede depender ni de la
variable x ni de la variable y, es decir, debe ser uniforme en un plano horizontal. Como en este caso, p no depende
ni de x ni de y, dp/dz = ∂p/∂z y de la última expresión se obtiene la conocida ecuación general de la estática de
fluidos, dp = −ρgdz. Esta ecuación queda, para un fluido incompresible (cuya densidad ρ no varía), en la conocida
forma p = p0 − ρgz.

3.3.4. Aplicación a un fluido en reposo en un sistema de referencia no inercial

Superficie libre en un depósito que se mueve con movimiento rectilíneo acelerado.

Estudiemos ahora el caso que se presenta en la figura 3–7, donde un depósito que contiene un líquido de densidad ρ
se mueve con movimiento acelerado.


(a)

una fuerza ficticia igual a −a además de la fuerza gravitatoria.

−
dv
dt

ρ ∂x
=a
1
a

= 0 = − ∇p + fext .
ρ
(b) z

a
x

En el sistema de referencia no inercial, fext = −aux − guz e igualando las distintas componentes se tiene

1 ∂p
⇒
∂p
∂x
= −ρa
g

Figura 3–7 Se tiene un depósito que se mueve con una aceleración a (a) y se quiere saber la distribución de presiones en el líquido y la ecuación
de la superficie libre. Visto el problema desde el sistema de referencia no inercial que se mueve con el depósito (b), sobre éste y su contenido actúa

En vez de estudiar el problema desde el sistema de referencia inercial lo haremos desde un sistema de referencia que
se mueve con el depósito. Visto desde este sistema de referencia, aparece una fuerza ficticia con aceleración −a en el
sentido contrario a a, y el sistema se encuentra en reposo. Las ecuaciones de Euler aplicadas a este problema nos dan

1 ∂p ∂p
− =0 ⇒ =0
ρ ∂y ∂y
1 ∂p ∂p
− =g ⇒ = −ρg
ρ ∂z ∂z

La variación infinitesimal de presión entre dos puntos muy próximos vendrá dada por

∂p ∂p ∂p
dp = dx + dy + dz = −ρadx − ρgdz,
∂x ∂y ∂z

que, tomando el origen de coordenadas como se indica en la figura da, al integrar,

p(x, z) = −ρax − ρgz + p0 (3–19)

donde p0 es la presión en el origen de coordenadas, en este caso la presión atmosférica. Para un valor de x dado, según
aumenta la profundidad (z cada vez más negativo), la presión aumenta linealmente con la profundidad, mientras que
para un z dado, como la superficie libre no es horizontal en este sistema de referencia, la presión disminuye según
aumenta x.
La superficie libre será el lugar geométrico de los puntos que se encuentran a la presión atmosférica, es decir,
aquéllos para los que al ser p = p0 se cumple

p − p0 = 0 = −ρax − ρgz

lo que da como ecuación de la superficie libre
a
z(x) = − x, (3–20)
g
que es un plano inclinado un ángulo arctan a/g.


(a) (b)

Agustin Martin Domingo ω

Superficie libre en un depósito circular que gira sobre su eje.
z

r

Figura 3–8 Superficie libre en un depósito circular que gira sobre su eje.
r

g
ωr2

Estudiaremos ahora otro caso de sistema acelerado, pero en este caso con aceleración angular. Consideremos un
depósito cilíndrico colocado en posición vertical que gira alrededor de su eje con una velocidad angular ω.
En el sistema de referencia del cilindro que rota aparece una fuerza ficticia de masa, la fuerza centrífuga, como
consecuencia de que este sistema de referencia no es un sistema de referencia inercial. Así, la aceleración total debida
a las fuerzas de masa en el sistema de referencia no inercial es fext = ω 2 rur − guz y las distintas componentes de las
ecuaciones de Euler (en coordenadas cilíndricas) son:

−
1 ∂p
ρ ∂r
+ ω2r = 0 ⇒
∂p
∂r
= ρω 2 r

1 1 ∂p ∂p
− =0⇒ =0
ρ r ∂θ ∂θ
1 ∂p ∂p
− −g =0⇒ = −ρg
ρ ∂z ∂z

Una variación infinitesimal de p en el entorno de un punto se representaría por

∂p ∂p ∂p
dp = dr + dθ + dz = ρω 2 rdr − ρgdz
∂r ∂θ ∂z
que, integrando en el sistema de ejes de la figura, da

1
p(r, z) = p0 + ρω 2 r2 − ρgz (3–21)
2
donde p0 es la presión en el origen de coordenadas, en este caso la atmosférica. Para un z dado la presión aumenta
cuadráticamente con r, ya que cada vez hay más fluido encima del punto, mientras que para un r dado, la presión
aumenta linealmente con la profundidad.
La superficie libre será el lugar geométrico de los puntos que se encuentren a la presión p0 , en este caso la atmosfé-
rica. Para estos puntos se tiene
1
p − p0 = 0 = ρω 2 r2 − ρgz
2
lo que da, para la ecuación de la superficie libre,

ω 2 r2
z(r) = (3–22)
2g

que es un paraboloide (y en el plano de la figura se ve como una parábola).


3.3.5. Vorticidad

Volvamos sobre la ecuación de Euler (3–17) y discutamos sobre la aceleración dv/dt. En principio podría pensarse
que la expresión para la aceleración es de la sencilla forma ∂v/∂t, pero esto no es así. ∂v/∂t representa la rapidez con
que la velocidad cambia en un punto fijo del espacio, pero lo que necesitamos es la variación de la velocidad de un
elemento particular del fluido. Si la trayectoria tiene una curvatura o el fluido se acelera de un punto a otro, aparece
una componente de aceleración del elemento aunque la velocidad en un punto se mantenga constante. Así, la variación
de velocidad de un punto a otro infinitamente próximo viene dada por

dv =

dt
=
∂v
∂t
de forma que la aceleración dv/dt se escribe como
dv ∂v
∂t
+
dt +

dv
dt
∂v
∂x
∂v
∂x

vx +
dx +

∂v
∂y
∂v
∂y

vy +

que puede escribirse, representando ∇ como un vector, en la forma

= (v · ∇)v +
dy +

∂v
∂t
∂v
∂z

.
∂v
∂z

vz ,
dz

Debe observarse que es posible una aceleración aunque ∂v/∂t sea 0, es decir, aunque la velocidad del fluido en un
punto dado no esté variando. Por ejemplo, un elemento de agua que fluye en círculo a una velocidad constante está
acelerándose sometido a una aceleración centrípeta aunque la velocidad en un punto dado no cambie.
Se denomina aceleración convectiva o aceleración de transporte al primero de los sumandos que aparecen en la
(3–23)

ecuación (3–23) y corresponde al ritmo de cambio de la velocidad debido al cambio en la posición en que se mide. El
segundo sumando es la denominada aceleración local y corresponde al ritmo de cambio de la velocidad en un punto
concreto del flujo.
Utilizaremos la siguiente identidad del análisis vectorial,
1
(v · ∇)v = (∇ × v) × v + ∇(v · v)
2
y definiremos un nuevo campo vectorial Ω como el rotacional de v,

Ω = ∇ × v, (3–24)

lo que nos permite escribir la identidad vectorial anterior como
1
(v · ∇)v = Ω × v + ∇v 2 .
2
Si la fuerza mfext admite un potencial, la aceleración fext también admitirá uno. Así, la aceleración fext puede
escribirse en función de su potencial Φ como fext = −∇Φ. Si la única fuerza de masa que actúa es la gravitatoria, el
potencial es simplemente Φ = gz. En función del potencial, la ecuación de movimiento (3–17) queda en la forma
dv ∂v 1 1
= + Ω × v + ∇v 2 = − ∇p − ∇Φ (3–25)
dt ∂t 2 ρ

Se denomina vorticidad al campo vectorial Ω. Si la vorticidad es cero en cualquier punto (es decir, si ∇ × v = 0) se
dice que el flujo es irrotacional.
Retomemos ahora un concepto habitual en electromagnetismo, la circulación de un campo vectorial. Ésta se define
como la integral de línea, a lo largo de cualquier camino cerrado del campo vectorial correspondiente multiplicado
por el elemento de longitud a lo largo de ese camino. En particular, la circulación del vector velocidad v a lo largo de
cualquier camino cerrado en un fluido se define como

v dl.

3.4. Conservación de la energía. Ecuación de Bernoulli. 37

Utilizando el teorema de Stokes esta circulación por unidad de área es igual a ∇ × v. Así, la vorticidad Ω es la
circulación del vector velocidad alrededor de un área unidad (perpendicular a la dirección de Ω). También puede
deducirse que un elemento de fluido en un punto del líquido rotará con una velocidad angular Ω/2.
Particularicemos ahora al caso de un fluido incompresible en el que nuestro interés se centra en el campo de velo-
cidades ya que en estas condiciones es posible eliminar la presión de las ecuaciones anteriores. Para ello se toma el
rotacional de los dos miembros de la ecuación (3–25). Como el rotacional de cualquier gradiente es cero y en un
líquido la densidad es constante, se tiene, utilizando la ecuación de continuidad

para un fluido incompresible,

ecuación que, junto con las ecuaciones

y
∂Ω
∂t
∇v = 0

+ ∇ × (Ω × v) = 0,

∇v = 0

Ω=∇×v
determina completamente el campo de velocidades v. Así, si conocemos la vorticidad Ω en un determinado instante
se tiene el rotacional del vector velocidad y como además conocemos que la divergencia es cero, dada esa situación
física es posible determinar la velocidad v en cualquier instante y lugar. Conociendo v la ecuación (3–26) permite
obtener la vorticidad Ω en un instante de tiempo posterior. Utilizando de nuevo (3–24) se tiene la nueva velocidad y
vuelta a empezar. Hay que notar que siguiendo este procedimiento se obtiene información únicamente sobre el campo
de velocidades y se ha perdido por completo la información sobre la presión.
(3–11)

(3–26)

(3–11)

(3–24)

Veamos un caso particular de lo anterior. Si Ω es cero en todo el fluido en un determinado instante de tiempo, ∂ Ω/∂t
también lo será, según la ecuación (3–26), y por tanto Ω será cero en cualquier instante de tiempo inmediatamente
posterior. Así, el flujo es permanentemente irrotacional. Si la rotación del flujo es inicialmente cero, siempre será
irrotacional y las ecuaciones a resolver serán

∇v = 0
∇×v =0

similares a las ecuaciones de los campos electrostático y magnetostático.

3.4. Conservación de la energía. Ecuación de Bernoulli.

Hemos obtenido la ecuación de continuidad para flujo estacionario y hemos demostrado que, para un flujo incompre-
sible vA es una constante. Si además, el trabajo realizado por las fuerzas no conservativas es despreciable, la energía
mecánica se conserva.
Estudiaremos esto de dos formas, en primer lugar utilizando la ecuación de movimiento del fluido perfecto escrita
en función de la vorticidad, (3–25) y en segundo lugar mediante el análisis de las variaciones de energía en el fluido
contenido dentro de un tubo de flujo estrecho y del trabajo mecánico realizado sobre el mismo,

Obtención a partir de las ecuaciones de Euler

Consideremos la ecuación (3–25) para el caso particular de flujo estacionario en el que ∂v/∂t = 0 y recordemos
que en flujo estacionario, las trayectorias reales de las partículas (la senda) tienen lugar a lo largo de las líneas de flujo.
Multiplicando escalarmente la ecuación (3–25) por v el término v(Ω × v) se anula y queda

p 1 2
v·∇ + v +Φ =0 (3–27)
ρ 2


o
Si consideramos un sistema de ejes coordenados que tiene uno de ellos (denominémosle xv ) dirigido a lo largo de una
línea de flujo y los otros dos en direcciones perpendiculares, la única componente del gradiente que no se anula en el
producto escalar es precisamente la componente en la dirección del flujo, quedando la ecuación anterior en la forma

p 1 2
v
∂
∂xv

∂
∂xv

+ v + Φ = cte
ρ 2

∇
p 1 2
+ v +Φ
ρ 2

p 1 2
+ v +Φ
ρ 2

p 1 2
+ v +Φ
ρ 2
= 0.

=0

Esta ecuación nos dice que a lo largo de una línea de corriente el término entre corchetes debe ser constante,

en una línea de flujo

conociéndose a esta expresión como ecuación de Bernoulli. Esta ecuación es válida a lo largo de una línea de flujo

de donde se deduce que en este caso particular de flujo irrotacional,
p 1 2
+ v + Φ = cte
ρ 2
= 0.

en el flujo
(3–28)

(3–29)

para un flujo estacionario e incompresible que conserva la energía. La constante será en general distinta para distintas
líneas de flujo.
Un caso particular interesante se tiene cuando el flujo es además irrotacional. En este caso la vorticidad es nula en
cualquier punto del flujo y la ecuación (3–25) queda como

(3–30)

para puntos en cualquiera de las líneas de flujo del fluido. Si como es habitual la única fuerza de masa que actúa sobre
el fluido es la atracción gravitatoria, Φ = gz y la ecuación de Bernoulli queda
p 1 2
+ v + gz = cte (3–31)
ρ 2
o, escrito en la forma habitual,
1
p + ρgz + ρv 2 = cte, (3–32)
2
o también como
p v2
+z+ = cte, (3–33)
ρg 2g
En esta última forma cada uno de los elementos de la ecuación de Bernouilli tiene dimensiones de longitud. Física-
mente eso representa energías por unidad de peso del fluido. Además, cada uno de los sumandos de (3–33) recibe un
nombre especial,

p
es la altura de presión. Representa la energía por unidad de peso debida a las fuerzas de presión.
ρg
z es la altura geométrica y representa la energía potencial gravitatoria del fluido por unidad de peso.
v2
2g es la altura de velocidad o altura cinética y representa la energía cinética por unidad de peso del fluido en un punto.

p
La suma de las alturas de presión y geométrica + z recibe el nombre de altura piezométrica y corresponde a la
ρg
altura que alcanzaría el líquido en un tubo vertical abierto a la atmósfera (piezómetro) colocado en ese lugar del flujo.
Asimismo se denomina altura total a la suma de los tres elementos de la ecuación de Bernoulli, correspondiendo a la
energía total por unidad de peso del fluido.

3.4. Conservación de la energía. Ecuación de Bernoulli. 39

p1 A1
1

v1 dt

A1
v1

A′
1
v2
A2
2

v2 dt
A′
2

z2 − z1
p2 A2

Figura 3–9 Obtención de la ecuación de Bernoulli para una línea de flujo. En la figura se representa un tubo de flujo estrecho de forma que la
velocidad pueda considerarse uniforme en cada sección.

Obtención a través de la conservación de la energía

Veamos cómo la ley de conservación de la energía relaciona la altura del fluido, el módulo de la velocidad del flujo
y la presión para puntos que se encuentran a lo largo de una línea de flujo. Para ello consideraremos un tubo de flujo
estrecho alrededor de la línea de flujo considerada, como se muestra en la figura 3–9. De este modo la velocidad en
las secciones transversales del tubo de flujo puede ser considerada uniforme.
En primer lugar hallaremos el trabajo realizado en un intervalo de tiempo dt sobre el fluido que está en la región
limitada por las secciones A1 y A2 y el tubo de corriente. El fluido que se encuentra a la izquierda de la superficie A1
ejerce sobre ésta una fuerza p1 A1 , perpendicular a la superficie. En el intervalo de tiempo dt esta fuerza realizará un
trabajo
F1 dx1 = p1 A1 v1 dt
sobre el volumen de fluido considerado, aumentando en esa cantidad la energía del mismo. Análogamente, el fluido
que se encuentra a la derecha de la superficie A2 ejerce sobre ésta una fuerza p2 A2 , por lo que el trabajo realizado por
esta fuerza sobre el elemento de volumen considerado es

F2 dx2 = −p2 A2 v2 dt,

por lo que el trabajo total realizado sobre el fluido es

δW = p1 A1 v1 dt − p2 A2 v2 dt, (3–34)

expresión válida para cualquier tipo de fluido. Para el caso concreto de un fluido incompresible, se cumple v1 A1 =
v2 A2 , por lo que la expresión anterior puede escribirse como

δW = (p1 − p2 )v1 A1 dt = (p1 − p2 )v2 A2 dt

o simplemente como
δW = (p1 − p2 ) dV (3–35)


donde dV es el volumen del elemento infinitesimal.
Este trabajo realizado sobre el fluido debe traducirse en un aumento de su energía total. Así, hay un cambio de
la energía cinética en la región A1 -A′ que se acaba de abandonar y en la región A2 -A′ que se acaba de ocupar, y

2

que, al simplificar queda

o, escrito en la forma habitual,

o también como
1

1 2

p1
ρg
2

1 2

v2
+ z1 + 1 =
2g
1 2

p2
ρg
2

1 2

v2
1 2
p1 − p2 = ρgz2 + ρv2 − ρgz1 − ρv1
2 2

p1 + ρgz1 + ρv1 = p2 + ρgz2 + ρv2 = cte,
2 2

+ z2 + 2 = cte,
2g
2
también un cambio en la energía potencial. En el volumen dV , que al ser el fluido incompresible, es el mismo en las

cambio en la energía total queda entonces como

1 2 1 2
dU = gz2 dm + v2 dm − gz1 dm + v1 dm
1
dos regiones la masa es ρdV , la energía potencial es gzdm (con z la altura del punto) y la energía cinética 2 v 2 dm. El

y debe ser igual al trabajo realizado sobre el fluido, dU = δW , ya que no se produce ni disipación ni otro tipo de
intercambio de calor.
1 2
dU = gz2 dm + v2 dm − gz1 dm + v1 dm = δW = (p1 − p2 )dV = (p1 − p2 )
2
dm
ρ

expresión conocida como ecuación de Bernoulli, válida a lo largo de una línea de flujo para un flujo que conserva la
(3–36)

(3–37)

(3–38)

(3–39)

energía, estacionario en un flujo incompresible.
Ya se ha comentado en la sección anterior el significado físico de cada uno de los términos que aparecen en la
Ecuación de Bernoulli (3–33) o (3–39).
Cuando el flujo es irrotacional, esta expresión es válida entre dos puntos de un mismo tubo de flujo, y no solamente
entre puntos de la misma línea de flujo.

3.5. Aplicaciones de las ecuaciones anteriores.

3.5.1. Velocidad de salida por un orificio (teorema de Torricelli).

Aplicando la ecuación de Bernoulli entre los puntos 1
1 y2
1 2 1 2
ρv + pa + ρgh = ρv2 + pa + ρgh
2 1 2
Como S1 ≫ S2 , de v1 S1 = v2 S2 se tiene que v1 ≃ 0,
h por lo que la velocidad de salida por un orificio en un
depósito abierto a la atmósfera es:

v2 = 2gh (3–40)
2
Obsérvese que la expresión obtenida es válida única-
mente cuando el depósito está abierto a la atmósfera
h0 y no cuando el depósito está a presión o hay una capa
de otro líquido encima de la del líquido que está desaguando.

3.5. Aplicaciones de las ecuaciones anteriores. 41

3.5.2. Tiempo de vaciado de un depósito.

dy
A

a
y y0

3.5.3. Presión dinámica y presión estática. Tubo de Pitot
En un intervalo de tiempo dt, el nivel del líquido en el
depósito desciende en dy. La variación de volumen será
entonces:

y0
dy
y 1/2

de donde queda
−Ady = Q dt =

Agrupando e integrando,
y
=−

t=
√
a 2g
A

A
a

tvac =
0
t

2 1/2
y
g 0
− y 1/2
2gy a dt

dt = 2y 1/2
y

y0

Para que el depósito se vacíe completamente se necesi-
taría un tiempo

A
a
2
g
y0
=−
at
A
2g

(3–41)

(3–42)

El tubo de Pitot se utiliza para medir la velocidad de un fluido. Para ello se utilizan dos aberturas, A y B, perpendi-
cular y paralela a la dirección del flujo. Si pC es la presión en el gas sobre el fluido (que puede regularse mediante la
llave), claramente se tiene
pA = ρg(H + h) + pC
⇒ pA − pB = ρgh (3–43)
pB = ρgH + pC

En B, la presión será la presión estática de la corriente fluida, pB , mientras que en A el fluido está estancado con
velocidad 0 y la presión es la presión dinámica. Como ambos puntos están a la misma altura, se tiene

1
pA + 0 + 0 = pB + ρv 2 + 0
2
(dicho en otras palabras, a la hora de calcular la altura total en A hay que tener en cuenta el aumento de altura debido
al fluido que intenta entrar y no puede). Así, se tiene

2(pA − pB ) 2(pA − pB )
v2 = ⇒ v= . (3–44)
ρ ρ

Esto se puede expresar en función de la diferencia de alturas h en la forma

v= 2gh (3–45)

Los tubos de Pitot para gases tienen la forma que se muestra en la figura, y en ellos se utiliza un fluido manométrico
de densidad ρ′ . En este caso, la velocidad del fluido (de densidad ρ) toma la forma

2ρ′ gh
v= (3–46)
ρ


v

A
v
A
h

B
H
C

B

Figura 3–10 Tubo de Pitot.

h

Figura 3–11 Tubo de Pitot para gases.

3.5.4. Efecto Venturi. Venturímetro.

El efecto Venturi es una consecuencia inmediata de la ecuación de continuidad y de la ecuación de Bernoulli.
Consideremos dos puntos A y B a la misma altura, pero con distintas secciones SA y SB . La ecuación de Bernoulli
nos da
1 2 1 2
pA + ρvA = pB + ρvB
2 2
que implica, al ser SA > SB

1 1 2 S2
pB − pC = ρ(vB − vC ) = ρvB 1 − B
2 2
2 > 0. (3–47)
2 2 SA

Como se ve, la presión es más débil donde la velocidad es mayor, es decir, donde la sección es menor.
Una aplicación inmediata de este efecto es el denominado venturímetro. Consiste en un calibrador colocado en una
tubería para medir la velocidad de ﬂujo de un líquido. Despejando la velocidad de la ecuación (3–47) se tiene
2
2 2(pB − pC ) SA
vB = 2 − S2
ρ SA a

que da
2
SB 2 S 2 2(pB − pC ) SA2
2
vA = 2 vB = B2 2 2
SA SA ρ SA − Sa

3.5. Aplicaciones de las ecuaciones anteriores. 43

Agustin Martin Domingo 1
pA

v

v′
pB

Figura 3–12 Efecto Venturi en un estrechamiento.

A

h
2
a

ρ′
ρ
v
pC

Figura 3–13 Relación entre la velocidad y el estrechamiento.

resultando ﬁnalmente, para la velocidad del ﬂujo

2
2SB (pA − pB ) ρ′ gh2
vA = 2 2 = SB 2 2 = (3–48)
ρ(SA − SB ) ρ(SA − SB )

Otra aplicación de este efecto es la denominada trompa de agua, que se utiliza para hacer bajos vacíos de forma limpia.
Este mismo efecto se utiliza también en pulverizadores.

Capítulo 4

Dinámica de fluidos reales. Flujo de fluidos
en tuberías.

4.1.

4.2.

4.3.

4.4.
Introducción. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.1.1. Concepto de viscosidad. . . . . . . . . . . . . . . . . . . . . . . .
4.1.2. Régimen laminar y régimen turbulento. . . . . . . . . . . . . . . .
Efecto de la viscosidad en los fluidos reales. . . . . . . . . . . . . . . . .
4.2.1. Fluidos newtonianos. . . . . . . . . . . . . . . . . . . . . . . . . .
4.2.2. De qué depende la viscosidad . . . . . . . . . . . . . . . . . . . .
4.2.3. La condición de no deslizamiento . . . . . . . . . . . . . . . . . .
4.2.4. Fluidos no newtonianos. . . . . . . . . . . . . . . . . . . . . . . .
Ecuaciones de movimiento de los fluidos reales. . . . . . . . . . . . . . .
4.3.1. Fluido perfecto o flujo no viscoso. . . . . . . . . . . . . . . . . . .
4.3.2. Fluido newtoniano. Las ecuaciones de Navier-Stokes. . . . . . . . .
Distribución de velocidades y tensiones cortantes. Fórmula de Poiseuille.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
45
45
46
47
47
48
49
50
51
53
53
54
4.5. Número de Reynolds. Concepto de capa límite. . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
4.6. Pérdidas de carga. Generalidades. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
4.6.1. Pérdidas de carga lineales. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
4.6.2. Pérdidas de carga singulares. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
4.6.3. Variación de altura debida a una máquina. . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
4.6.4. Representación gráfica de la pérdida de carga. . . . . . . . . . . . . . . . . . . . . . . . . . 63
4.7. Redes de tuberías . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
4.8. Cavitación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
4.9. Golpe de ariete o choque hidráulico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65

4.1. Introducción.

A menudo, las fuerzas no conservativas en un fluido no pueden despreciarse. Estas fuerzas dan lugar a una disi-
pación de la energía mecánica del fluido en energía interna del mismo, igual que las fuerzas de rozamiento disipan
energía de un bloque deslizante en energía interna del bloque y de la superficie sobre la que se desliza (aumentan sus
temperaturas). Como consecuencia de ésto, cuando un fluido viscoso circula por una tubería horizontal uniforme, la
presión a lo largo de una línea de flujo disminuye, al contrario de lo que ocurre en un fluido perfecto (fig. 4–1). Debe
notarse que al no variar la sección de la tubería, la velocidad del flujo en la misma no varía, al contrario de la presión.

4.1.1. Concepto de viscosidad.

Todo fluido real posee la propiedad de la viscosidad. La viscosidad es la causa de la resistencia que aparece en el flujo
de un fluido en tuberías y canales y está relacionada con la oposición de un fluido a las deformaciones tangenciales,
en las que se produce rozamiento entre distintas partes del fluido. Por ejemplo, a causa de las fuerzas viscosas de

45

46 Capítulo 4. Dinámica de fluidos reales. Flujo de fluidos en tuberías.

Agustin Martin Domingo Fluido perfecto

turbulento
∆h
Fluido real

Figura 4–1 La presión disminuye a lo largo del flujo en un fluido real, al contrario de lo que ocure en un fluido perfecto

laminar

Figura 4–2 Régimen laminar y régimen turbulento. En régimen laminar el movimiento de las capas del fluido es ordenado, en régimen turbulento
no.

rozamiento en una tubería, las capas de fluido más próximas a la pared (más lentas) frenan a las capas más interiores
(más rápidas), apareciendo una distribución de velocidades.
El efecto de estas fuerzas de rozamiento se traduce en una disipación de energía que debe tenerse en cuenta en la
ecuación de Bernoulli (3–33) en la forma siguiente, si el flujo va de 1 a 2,
p1 v2 p2 v2
+ z1 + 1 − ∆h12 = + z2 + 2 , (4–1)
ρg 2g ρg 2g
donde ∆h12 es la energía disipada entre 1 y 2 por unidad de peso del fluido, a la que se denomina pérdida de carga.

4.1.2. Régimen laminar y régimen turbulento.

Dependiendo de la importancia relativa entre las fuerzas viscosas y las fuerzas de inercia, el régimen de movimiento
de un fluido puede ser de dos tipos,

Régimen laminar o de Poiseuille: El régimen del flujo es laminar cuando el flujo tiene un movimiento ordenado, en
el que las partículas del fluido se mueven en líneas paralelas, sin que se produzca mezcla de materia entre las
distintas capas.
Régimen turbulento o de Venturi: El régimen de movimiento de un fluido es turbulento cuando el fluido presenta
un movimiento desordenado con mezcla intensiva entre las distintas capas.

Para distinguir conceptualmente ambos régimenes es conveniente recordar la forma en que experimentalmente se
observa el paso de un régimen de movimiento a otro. Consideremos una tubería cilíndrica de sección constante por la
que circula un fluido a una cierta velocidad. Si introducimos en el flujo una aguja fina por la que lentamente inyectamos
tinta como se muestra en la figura 4–2 podemos tener dos comportamientos:

• Cuando el régimen del flujo es laminar, la tinta avanza con el agua como un fino hilo, mezclándose lentamente
con ésta por procesos de difusión. La mezcla entre el agua y la tinta es muy lenta. Esta situación se tiene para
diámetros pequeños y velocidades bajas del fluido.
• A medida que aumenta la velocidad del flujo, llega un momento en que el régimen de movimiento pasa a ser
turbulento. Cuando el régimen de movimiento del flujo es turbulento, la tinta se mezcla de forma inmediata con
el agua, diluyéndose en ésta.

4.2. Efecto de la viscosidad en los fluidos reales. 47

dn

n
n

τ + dτ

τ
v + dv

v

pared
v(n)

β

x, v

Figura 4–3 En un fluido real en movimiento aparecen fuerzas viscosas de rozamiento Froz entre las distintas capas del fluido, con una fuerza por
unidad de superficie τ . Esto hace que la velocidad en una sección del fluido no sea uniforme.

Dependiendo del tipo de flujo, laminar o turbulento, las fuerzas de rozamiento que aparecen serán de distintos tipos.
En el caso de régimen laminar, las fuerzas que se ejercen entre las distintas capas del fluido son tangentes a la dirección
del movimiento. Por el contrario, cuando se está en régimen turbulento aparece una segunda contribución a la fuerza
de rozamiento debida a la mezcla entre las distintas capas.

4.2. Efecto de la viscosidad en los fluidos reales.

4.2.1. Fluidos newtonianos.

Comencemos considerando el caso en el que el movimiento de un fluido tiene lugar de forma bien ordenada,
moviéndose éste en líneas paralelas, en lo que acabamos de denominar régimen laminar, teniéndose que las fuerzas de
rozamiento son tangentes a la dirección del movimiento. Estas fuerzas de rozamiento presentan en general distintos
comportamientos dependiendo de la rapidez con que varía la velocidad en la dirección perpendicular a la superfi-
∂v
cie, rapidez representada por . Una situación muy frecuente se tiene cuando esta relación es de proporcionalidad,
∂n
denominándose fluidos newtonianos a los fluidos que satisfacen esta condición.
Así, en los fluidos newtonianos, las fuerzas de rozamiento que aparecen entre las capas de fluido en movimiento una
con respecto a la otra, para una superficie de contacto dS, vienen dadas por la fórmula de Newton,

∂v
dFroz = −µ dS, (4–2)
∂n
que en función de la tensión cortante (fuerza por unidad de superficie) queda en la forma

∂Froz ∂v
τ =− =µ , (4–3)
∂S ∂n
∂v
conocida como ley de Newton. En esta expresión es la derivada que caracteriza la variación de la velocidad en
∂n
la dirección transversal al movimiento y la constante de proporcionalidad µ es la viscosidad dinámica o coeficiente
de viscosidad dinámica que depende de las propiedades del fluido, indicando los signos negativos que se trata de una


fuerza de rozamiento que se opone al movimiento. Las unidades de la viscosidad dinámica en el sistema internacional
son [kg m−1 s−1 ] o [N s m−2 ≡Pa s]. La viscosidad dinámica µ depende mucho de la temperatura y poco de la presión,
como se trata en más detalle en el apartado 4.2.2. Solamente los fluidos newtonianos (por ejemplo agua, aire, Hg,
productos petrolíferos típicos, etc...) cumplen esta ley de Newton. Algunos ejemplos de fluidos no newtonianos son

• Fluidos biológicos.
• Pinturas al aceite.
• Soluciones coloidales.
• Soluciones acuosas de combinados de gran masa molecular.

La ley de Newton, que en principio fue propuesta como hipótesis para una primera aproximación, está ampliamente
comprobada experimentalmente en buen número de fluidos y de hecho es la más usual en Mecánica de Fluidos cuando
se quieren calcular tensiones tangenciales en el movimiento laminar de un fluido.
Se define la viscosidad cinemática ν en la forma ν = µ/ρ, y sus dimensiones en el sistema internacional son [ν] =
m2 s−1 . La viscosidad cinemática no depende de las características dinámicas del fluido.
En los sistemas CGS e internacional las unidades de las viscosidades dinámica y cinemática son, respectivamente,

µ
ν
CGS
1 poise = 1 dy·s·cm−2 = 1 g/cm·s = 10−1 kg/m·s
1 stokes = 1 cm2·s−1 = 10−4 m2 s−1

4.2.2. De qué depende la viscosidad
MKS
1 kg/m·s = 1 Pa·s = 10 poises
1 m2 s−1 = 104 stokes

En general la viscosidad no depende apenas de la presión. Sin embargo, sí depende de la temperatura, y de una forma
distinta para líquidos y gases. La viscosidad de un líquido disminuye con un aumento de la temperatura, mientras que
la viscosidad de un gas tiene el comportamiento contrario, aumenta con la temperatura. Esto es debido a los distintos
orígenes de la viscosidad en ambos casos.
La movilidad de las moléculas de un líquido es más limitada, con una presencia importante de fuerzas de cohe-
sión, relacionada con las fuerzas de rozamiento y por tanto con la viscosidad, ya que son estas fuerzas de cohesión
las que hacen que durante el movimiento del líquido unas capas frenen a otras. Cuando la temperatura aumenta la
cohesión entre las moléculas disminuye, al aumentar la energía cinética de las mismas. Esto da como consecuencia
una disminución de la viscosidad.
La viscosidad de los líquidos se aproxima por la siguiente ley empírica,

µ=A e −BT
(4–4)

donde tanto A como B son magnitudes que se obtienen al ajustar la expresión anterior con los resultados experimen-
tales de un líquido en particular.
En el caso de los gases las moléculas tienen una movilidad mucho mayor y el origen de la viscosidad está relacionado
con el paso de moléculas de diferente energía entre distintas capas. Imaginemos dos capas adyacentes, cada una de
ellas con una distinta velocidad de movimiento. En promedio las moléculas de la capa más rápida se mueven a mayor
velocidad que las moléculas de la capa más lenta. Como en el sentido del flujo las moléculas de la capa más rápida se
mueven en promedio a una velocidad mayor que las de la capa más lenta, el paso de moléculas de la capa más rápida
a la más lenta se traducirá en una aceleración de la misma, mientras que el paso de moléculas de la capa más lenta a
la más rápida se traducirá en un frenado de estas últimas, dando este movimiento aleatorio lugar al comportamiento
viscoso del gas. Como superpuesto a este movimiento se tiene un movimiento aleatorio de las moléculas relacionado
con la temperatura del gas, pasarán moléculas de la capa más rápida a la más lenta y viceversa. Este paso será tanto
más rápido cuanto mayor sea la temperatura, y esto dará lugar a una mayor intensidad de la aceleración y el frenado a
que hacíamos mención en el párrafo anterior, aumentando por tanto la viscosidad con la temperatura.

4.2. Efecto de la viscosidad en los fluidos reales. 49

Tabla 4–1 Viscosidades dinámica y cinemática del agua a distintas temperaturas con sus correspondientes densidades. La viscosidad dinámica se
da en unidades de 10−6 Pa s=10−5 poise, la cinemática en unidades de 10−6 m2 /s y la densidad en kg/m3 .

T (◦ C)

0
10
20
30
40
50
Viscosidad
dinámica
1793
1309
1006
800
657
550

Fluido

Mercurio

Etanol
Metanol
Benceno
Gasolina
Glicerina
Densidad

999,841
999,7
998,2
995,65
992,21

Aceite de oliva
Aceite lubricante
Aire
Viscosidad
cinemática
1,793
1,309
1,008
0,803
0,662

T (◦ C)

0
50
100
20
20
20
20
20
20
20
20
T (◦ C)

60
70
80
90
100

Viscosidad
dinámica
1690
1410
1220
1192
591
649
290
83×104
84×104
26×104
18
Viscosidad
dinámica
469
406
357
315
284

Densidad

13595
13472
13351
791
810
899
680-780
1260
920
900-920
1,293
Densidad

983,21
977,79
971,8
965,31
958,35

Viscosidad
cinemática
0,124
0,105
0,091
1,506
0,730
0,722
0,427
658,7
913,0
279
13,92
Viscosidad
cinemática
4,770
4,152
3,673
3,263
2,963

Tabla 4–2 Viscosidades dinámica y cinemática de varios fluidos junto con su densidad. La viscosidad dinámica se da en unidades de
10−6 Pa s=10−5 poise, la cinemática en unidades de 10−6 m2 /s y la densidad en kg/m3 .

Dióxido de carbono 20 16 1,977 8,09
Nitrógeno 20 18,4 1,251 14,71
Oxígeno 20 20,9 1,429 14,62

La dependencia de la viscosidad de los gases con la temperatura se aproxima por una de las siguientes leyes empí-
ricas, conocidas como ley de Sutherland

µ0 (T /T0 )3/2 (T0 + S)
µ= (4–5)
T +S
o por la ley de la potencia
T n
µ = µ0 (4–6)
T0
En ambos casos µ0 es una viscosidad dinámica conocida a una cierta temperatura, mientras que los valores S y n se
ajustan con los datos experimentales de un líquido concreto.

4.2.3. La condición de no deslizamiento

Resumiendo, el esfuerzo cortante en un fluido viscoso newtoniano en régimen laminar es proporcional a la pendiente
de la velocidad y es máximo en la pared. Además, justo en la pared, la velocidad del fluido es cero con respecto a la
pared. A este hecho se le denomina condición de no deslizamiento y es una característica de los fluidos viscosos. Un
ejemplo sencillo de ésto se tiene en las aspas de un ventilador, que tienen polvo fino aunque el ventilador se encuentre
en movimiento, ya que justo en la superficie del ventilador el aire se encuentra en reposo. Al girar las aspas, éstas
pueden mover las partículas grandes de polvo pero no las pequeñas. Cuanto más potente sea el ventilador podrá mover
más partículas, pero siempre quedarán algunas pequeñas que no se mueven.


τ

viscoplástico

4.2.4. Fluidos no newtonianos.
α
dilatantes

pseudoplástico
newtonianos

dv/dn
τ

Figura 4–4 Reogramas para distintos tipos de fluidos.
β

Acabamos de ver que denominamos fluidos newtonianos a aquéllos en los que la dependencia entre las fuerzas de
rozamiento y el perfil de velocidades es lineal. Sin embargo, no todos los fluidos tienen el mismo comportamiento,
aunque la mayor parte de los gases y líquidos simples se comportan como fluidos newtonianos. Por ejemplo, pastas,
lodos y polímeros de alta densidad son ejemplos de fluidos que no tienen un comportamiento newtoniano. Deno-
minaremos fluidos no newtonianos o anómalos a aquellos fluidos que no siguen la ley de Newton del rozamiento
interno.
β

dv/dn

El comportamiento de los fluidos no newtonianos se representa mediante las curvas de corriente, curvas reológicas†
o reogramas, como se muestra en la figura 4–4, donde se representa la tensión cortante producida por las fuerzas de
rozamiento frente a la tasa de variación de la velocidad en una sección dv/dn. Estas curvas se construyen a partir de
datos experimentales obtenidos con viscosímetros de distintos tipos en distintas condiciones de funcionamiento.
Las curvas de corriente correspondientes a los fluidos newtonianos son lineales y se representan en el diagrama por
rectas que pasan por el origen de coordenadas, con la viscosidad igual a la pendiente de la recta,
τ
µ = tan α =
dv/dn
que caracteriza completamente las propiedades reológicas del fluido a temperatura y presión dadas, con independencia
del gradiente de velocidad.
Las curvas de corriente de los líquidos no newtonianos son muy diversas y en el caso general no lineales. Depen-
diendo del comportamiento de estos líquidos en el diagrama reológico se tienen distintos comportamientos no newto-
nianos:

Fluidos pseudoplásticos: En éstos, al aumentar la velocidad de variación de la velocidad la fuerza cortante viscosa
aumenta menos que con los newtonianos (véase la figura). Un ejemplo de materiales pseudoplásticos son
muchos lodos no newtonianos o fundidos de polímeros.
Fluidos dilatantes: Los fluidos dilatantes tienen el comportamiento contrario al de los fluidos pseudoplásticos. Muchas
pinturas son fluidos dilatantes.

Las curvas de corriente de los fluidos pseudoplásticos y dilatantes se describen muy bien en régimen laminar mediante
la expresión
n
dv
τ =κ (4–7)
dn
† Precisamente se denomina Reología a la rama de la Física de Materiales que estudia la respuesta de los fluidos a los esfuerzos.

4.3. Ecuaciones de movimiento de los fluidos reales. 51

donde κ y n son magnitudes constantes para un líquido dado. κ es una medida del grado de consistencia del líquido
(teniéndose mayor κ a mayor viscosidad) mientras que n representa el grado de comportamiento no newtoniano. Los
fluidos newtonianos corresponden al caso particular κ = µ y n = 1, los fluidos pseudoplásticos al caso n < 1 y los
dilatantes a n > 1.

Materiales no lineales de Bingham o viscoso-plásticos Otro tipo de fluido no newtoniano corresponde no descrito
por la ecuación anterior corresponde a los materiales lineales de Bingham o fluidos viscoplásticos. En estos
materiales el comportamiento es similar al de un sólido para tensiones cortantes por debajo de un determinado
umbral, es decir presenta una deformación fija para un cierto valor de la tensión cortante. Sin embargo, cuando
el esfuerzo cortante es mayor que ese umbral, el material presenta el comportamiento de un fluido viscoso
newtoniano, en la forma
τ = τ0 + µ
dv
dn
donde se conoce a τ como tensión de desplazamiento, tensión inicial de desplazamiento o límite de fluidez.

En el trabajo común con fluidos no newtonianos se introduce el concepto de viscosidad efectiva, que permite con

coordenadas con los puntos de la curva de corriente para las condiciones del flujo,

µe = tan β =
τ
dv/dn
(4–8)

ciertas limitaciones realizar los cálculos mediante las fórmulas de fluidos newtonianos. Para un líquido dado ésta
no es una magnitud constante, sino que depende de la tasa de cambio de la velocidad dv/dn y de la tensión de
desplazamiento τ0 . Su valor se calcula en el diagrama reológico como la pendiente de las rectas que unen el origen de

(4–9)

Hay que insistir en que los valores de la viscosidad efectiva se determinan solamente por el estado instantáneo del
flujo.
En los fluidos pseudoplásticos la viscosidad efectiva µe disminuye con el aumento de τ o dv/dn, dando la apariencia
de que se diluyen al fluir (al fluir aparece una τ y un dv/dn). Por el contrario en los fluidos dilatantes µe aumenta con
τ o dv/dn y por tanto parece que se espesan al fluir.
Otro tipo de materiales presentan una combinación de características viscosas y elásticas especial y se les denomina
fluidos viscoelásticos. Estos materiales se comportan como sólidos ante variaciones rápidas del esfuerzo aplicado, pero
se comportan como fluidos ante pequeños esfuerzos continuados. Los plásticos a temperatura ambiente sometidos a
una carga presentan un comportamiento viscoelástico, así como de forma muy espectacular una papilla espesa de
harina de maíz.

4.3. Ecuaciones de movimiento de los fluidos reales.

¿Qué forma tienen las ecuaciones del movimiento para un fluido viscoso? En el caso de un fluido viscoso tendremos
que tener en cuenta que, además de las fuerzas de presión y de las fuerzas de masa como la gravitatoria, sobre un
fluido viscoso en movimiento actúan las fuerzas viscosas de rozamiento.
El estado de tensiones en un punto del fluido puede expresarse mediante el tensor de tensiones

−p + τxx τxy τxz
 

[σij ] =  τyx −p + τyy τyz (4–10)
 

τzx τzy −p + τzz

donde cada elemento ij del tensor representa la fuerza por unidad de superficie que actúa en la dirección i sobre la
cara perpendicular al eje j. Esta tensión es en general inhomogénea, es decir, varía de un punto a otro del flujo. Si
consideramos el elemento de volumen infinitesimal de la figura, las fuerzas que actúan sobre las caras x y x + dx en
la dirección x serán
dFxx (x) = [p(x) − τxx (x)] dydz


−p(x) + τxx (x) −
∂p
∂x
dx +

x
∂τxx
∂x
dx

dy
τxy (y)
z
τxz (z) +
∂τxz
∂z

τxz (z)
dz

−p(x) + τxx (x)

τxy (y) +

y
∂τxy
∂y

dx
dy
dz

Figura 4–5 Componente x de las fuerzas que actúan sobre un ﬂuido real en movimiento en las distintas caras de un elemento de volumen.

sobre la cara que se encuentra a una distancia x del plano yz y

dFxx (x + dx) = [−p(x + dx) + τxx (x + dx)] dydz

sobre la cara que se encuentra a una distancia x + dx del plano yz. Si desarrollamos en serie esta última expresión, se
tiene
∂p ∂τxx
dFxx (x + dx) = −p(x) − dx + τxx (x) + dx dydz
∂x ∂x

quedando la fuerza neta que actúa en la dirección x sobre las caras perpendiculares a esta misma dirección

neta ∂p ∂τxx
dFxx = dFxx (x) + dFxx (x + dx) = − + dxdydz (4–11a)
∂x ∂x

Un cálculo análogo para las fuerzas en la dirección x que actúan sobre las caras perpendiculares a las direcciones y y
z da

neta ∂τxy
dFxy = dFxy (y) + dFxy (y + dy) = dxdydz (4–11b)
∂y

neta ∂τxz
dFxz = dFxz (z) + dFxz (z + dz) = dxdydz (4–11c)
∂z

quedando la fuerza neta en la dirección x como

neta ∂p ∂τxx ∂τxy ∂τxz
dFx = − + + + . (4–12a)
∂x ∂x ∂y ∂z

4.3. Ecuaciones de movimiento de los fluidos reales. 53

Un cálculo análogo da, para las fuerzas netas en las direcciones y y z

o en notación tensorial
neta
neta
dFy

neta
dFz

dFsup = −∇p +
=

=
−

−
∂p ∂τyx
∂y

∂z
+

+

∂τxx
∂x
∂τyx
∂x
∂τzx
∂x

∂ ∂ ∂
,
∂x ∂y ∂z

neta
∂x

∂p ∂τzx
∂x

+

+

+

,
+

+
∂τyy

∂τxy
∂y
∂τyy
∂y
∂τzy
∂y
∂y

∂τzy
∂y

+

+

+
+

+
∂τyz

∂τxz
∂z
∂τyz
∂z
∂τzz
∂z
∂z

∂τzz
∂z

i+

j+

k.

De una forma más compacta esta expresión se puede escribir en la siguiente forma

τxx τyx τzx
.



·  τxy τyy τzy 


τxz τyz τzz

(4–12b)

(4–12c)

Multiplicando estas ecuaciones respectivamente por los vectores unitarios i, j y k y sumando vectorialmente se tiene,
para la fuerza neta que actúa sobre las distintas superficies
neta
dFsup = −
∂p
∂x
∂p ∂p
i + − j + − k+
∂y ∂z

(4–13)

(4–14)

dFsup = −∇p + ∇îj .
τ (4–15)
De esta forma, la ecuación de conservación de la cantidad de movimiento para un fluido viscoso se escribe, en su
forma diferencial, como
dv ∂v
ρ = ρg(−k) − ∇p + ∇îj
τ =ρ + (v · ∇)v (4–16)
dt ∂t
donde en la parte derecha se ha descompuesto la aceleración en sus componentes local y convectiva (3–23).

4.3.1. Fluido perfecto o flujo no viscoso.

Obviamente para el caso de un fluido perfecto, todas las componentes del tensor τij son nulas y se recupera la
ˆ
ecuación de Euler para flujo no viscoso
dv
ρ = ρg(−k) − ∇p (3–17)
dt

4.3.2. Fluido newtoniano. Las ecuaciones de Navier-Stokes.

Para un fluido newtoniano los esfuerzos viscosos son proporcionales a la derivada de la velocidad y al coeficiente
de viscosidad. Para el caso de flujo incompresible, las componentes del esfuerzo viscoso para un fluido newtoniano
quedan
∂vx
τxx = 2µ (4–17a)
∂x
∂vy
τyy = 2µ (4–17b)
∂y
∂vz
τzz = 2µ (4–17c)
∂z


para las componentes normales, y

ρ

ρ

ρ
dvx
dt
dvy
dt
dvz
dt
τxy = τyx = µ

τxz = τzx = µ

τyz = τzy = µ

= ρgx −

= ρgy −

= ρgx −
∂p
∂x
∂p
∂y
∂p
∂z
+µ

+µ

+µ
∂x
∂vx
∂y
∂vz
∂x
∂vy
∂z

2

∂ 2 vy
∂x 2
+

+

+

∂ 2 vz
∂x2
+

+
∂vy
∂x
∂vx
∂z
∂vz
∂y

+
∂ 2 vx
∂y 2

∂ 2 vy
∂y 2

∂ 2 vz
∂y 2
+

+

+
∂ 2 vx
∂z 2
∂ 2 vy
∂z 2
∂ 2 vz
∂z 2
(4–17d)

(4–17e)

(4–17f)

para las cortantes. Sustituyendo estas expresiones en la ecuación anterior se obtienen las ecuaciones de Navier-Stokes

∂ 2 vx
(4–18a)

(4–18b)

(4–18c)

para un fluido newtoniano. Nótese que éste es un sistema de ecuaciones diferenciales en derivadas parciales no lineales
de segundo orden que tiene cuatro incógnitas, vx , vy , vz y p. Para resolverlo, necesitaremos una cuarta ecuación que
no es sino la ecuación de continuidad (3–11) para un flujo incompresible.

4.4. Distribución de velocidades y tensiones cortantes. Fórmula de Poiseuille.

Determinaremos ahora el caudal o gasto Q en una tubería cilíndrica de radio R como la representada en la figura 4–6
en condiciones de flujo laminar. A través de la tubería de sección constante fluye un líquido (homogéneo e incom-
presible) de viscosidad ν y las presiones en dos puntos 1 y 2 separados una distancia L son respectivamente p1 y
p2 .
Si consideramos régimen estacionario, al ser v = v(t) no hay aceleración en la dirección axial y la suma de fuerzas
en esta dirección es cero. Observando el equilibrio de fuerzas en un cilindro imaginario de radio r (distancia al eje del
cilindro) y longitud L dirigido según el eje mostrado en la figura, que forma un ángulo θ con la horizontal se tiene, en
la dirección axial,
F = 0 ⇒ p1 (πr2 ) − p2 (πr2 ) − τ Slat. − mg sen θ = 0.

Como la masa m es ρSL y L sen θ = z, se tiene

τ Slat. = (p1 − p2 )S − ρgzS.

Podemos expresar la diferencia de presiones p1 − p2 en función de la pérdida de carga ∆h12 a partir de la ecuación
(4–1), que teniendo en cuenta que v1 = v2 al ser S1 = S2 y que z = z2 − z1 da

p1 − p2 = ρg∆h12 + ρgz.

Como la sección de la tubería es constante, la pérdida de carga ∆h12 es igual a la variación de la altura piezométrica
entre los dos puntos, ∆hpiezo = ∆h12 , por lo que podemos escribir también
12

p1 − p2 = ρg∆hpiezo + ρgz.
12

Sustituyendo en la ecuación anterior se van los términos en z y se obtiene la distribución de tensiones cortantes

ρg∆hpiezo r
12
τ= (4–19)
2L

4.4. Distribución de velocidades y tensiones cortantes. Fórmula de Poiseuille. 55

v

p1 S

τ = −µ
∂v
∂r
=
p1
τ

ρg∆hpiezo r
12
2L
donde vm´x es la velocidad máxima sobre el eje, quedando
a
⇒−
τ
L

θ

mg

v

vm´ x
a
θ

Figura 4–6 Tubería cilíndrica por la que circula un fluido newtoniano.

dv =
0
r
p2 S

Si el régimen del flujo no fuera laminar habría dos componentes de la tensión cortante, la mostrada arriba y además
una originada en la mezcla de fluido. En régimen laminar sólo aparece ésta y aplicando a este flujo la ley de Newton,†
p2

ρg∆hpiezo rdr
12
2µL
R
r

z

ρgr2
vm´x − v =
a ∆hpiezo
12
4µL
Si tenemos en cuenta que para r = R (es decir sobre la pared) debe cumplirse v = 0 se tiene
ρgR2
vm´ x =
a ∆hpiezo
12
4µL
quedando la ley de distribución de velocidades para un fluido newtoniano en régimen estacionario y laminar en la
forma parabólica
ρg∆hpiezo 2
12
v(r) = (R − r2 ) (4–20)
4µL
coincidiendo el eje de la parábola con el eje de la tubería.
El gasto o caudal a través de una sección elemental dS en la cual la velocidad es constante viene dado por dQ = vdS,
siendo el caudal total la integral de los caudales elementales
R
πR4
vS = Q =
¯ vdS = v2πrdr = ρg∆hpiezo
12 (4–21)
S 0 8µL

Si comparamos con la expresión anterior para la velocidad máxima se observa que Q = 1 πR2 vm´x y si tenemos
2 a
además en cuenta que la velocidad media está relacionada con el caudal a través de Q = v S y que la sección es
¯
precisamente S = πr2 se tiene de forma inmediata que
Q ρgR2 1
v=
¯ = ∆hpiezo = vm´x
12 a (4–22)
S 8µL 2
para esta distribución parabólica de velocidades.
† El signo − aparece al tomarse los ejes respecto de r = 0


zona de fluido libre

v0
v0

v = v(y)

capa límite

Figura 4–7 Consideremos una pequeña parte de la sección de un flujo viscoso alejada de las paredes, en la que podamos considerar la velocidad
uniforme. Si introducimos una cuchilla en el flujo y observamos la distribución de velocidades en las cercanías de la cuchilla se observan dos
comportamientos bien diferenciados. Lejos de la cuchilla el efecto de las fuerzas de inercia predomina sobre el de las fuerzas viscosas y se tiene la
zona de fluido libre, mientras que cerca de la cuchilla el efecto predominante es el de las fuerzas de inercia y se tiene la denominada capa límite.

4.5. Número de Reynolds. Concepto de capa límite.

Resulta de gran interés conocer la influencia relativa de las fuerzas de inercia y viscosas en el movimiento del fluido
y caracterizar esta influencia.
El número de Reynolds es un parámetro adimensional que surge a partir de la relación entre las fuerzas de inercia
y las fuerzas viscosas en las ecuaciones de Navier-Stokes (4–18). El término de inercia ρ(v · ∇)v es el asociado con
la componente convectiva de la aceleración y dimensionalmente va como ρv 2 /L, mientras que el término viscoso
µ∇2 v es el asociado con la parte viscosa, y va dimensionalmente como µv/L2 , siendo L y v una longitud y una
velocidad característica respectivamente. Como ambos términos son dimensionalmente equivalentes, su cociente será
una magnitud adimensional.
Así, se define el número de Reynolds en en la forma
ρv 2 /L ρ¯D
v vD
¯
⇒ Re = = (4–23)
µv/L2 µ ν
donde v es la velocidad media característica del flujo, D es una longitud característica (típicamente el diámetro), ρ es
¯
la densidad del fluido, µ es la viscosidad dinámica y ν = µ/ρ es la viscosidad cinemática.
Se denomina velocidad característica a la velocidad por debajo de la cual toda turbulencia es amortiguada por la
acción de la viscosidad del fluido. Experimentalmente se observa que de forma aproximada para números de Reynolds
menores que 2300† el régimen de movimiento del fluido es laminar (predominan las fuerzas de rozamiento) mientras
que para números de Reynolds mayores que 2300 el régimen de movimiento es turbulento (predominan las fuerzas de
inercia). Experimentalmente se observa el paso de uno a otro régimen introduciendo un colorante en el flujo como se
mencionó anteriormente.
Para fluidos de pequeña viscosidad o en general para valores grandes del número de Reynolds, el movimiento del
fluido suele estudiarse por separado en dos zonas, como se muestra en la figura 4–7.

1. Zona de espesor pequeño (capa límite) en la que la influencia de las fuerzas viscosas es tan importante como la
† En realidad, el cambio de régimen tiene lugar entre 2000 y 2400, dependiendo de las condiciones concretas en que se realiza el proceso.

Algunos autores dan otros valores en ese rango como referncia para el cambio de régimen

4.6. Pérdidas de carga. Generalidades. 57

de las otras fuerzas.
∆h =
v2
2g

2. Zona donde la influencia de la viscosidad es despreciable (zona de fluido libre).

4.6. Pérdidas de carga. Generalidades.
Válvula cerrada

Válvula abierta (fluido perfecto)

Válvula abierta (fluido real)

Figura 4–8 Se comparan los niveles del agua en los piezómetros cuando la válvula está cerrada, y cuando, estando abierta, circula por ella un fluido
perfecto y un fluido real.

4.6.1. Pérdidas de carga lineales.

Las pérdidas de carga lineales son aquellas que tienen lugar de forma continua a lo largo del flujo, por el hecho
de que el fluido avance, aunque no cambien ni la dirección del flujo ni la sección de la tubería. Como un ejemplo,
consideremos un depósito que desagua a través de un tubo de sección constante en el cual se han colocado una serie
de piezómetros, como se muestra en la figura.
Si la válvula de salida está cerrada, el nivel en todos los piezómetros es el mismo e igual al nivel de agua en el
depósito. Al abrir la válvula se observan distintos comportamientos dependiendo de si tenemos un fluido perfecto o un
fluido real. Si lo que tenemos es un fluido perfecto, el nivel de agua en todos los piezómetros baja la misma cantidad,
y ésta es igual a la altura de velocidad
v2
∆h =
2g
Por el contrario, si lo que tenemos es un fluido real, la pérdida de carga se reparte de forma continua a lo largo del
tubo. Esto hace que el nivel del agua en los piezómetros vaya disminuyendo según se avanza en el tubo.
Estas pérdidas de carga son, en general, una función compleja de la geometría del sistema, de las propiedades
del fluido, y de las características del flujo. Experimentalmente se observa que las pérdidas de carga son en general
proporcionales al cuadrado de la velocidad promedio en la mayor parte de los casos de interés técnico. Esta observación
se expresa en función de la altura de velocidad en la conocida ley de Darcy-Weiesbach para las pérdidas de carga
lineales debidas al rozamiento,
L v2
¯
∆hl = f . (4–24)
D 2g
El factor de proporcionalidad f es el coeficiente de fricción que en general varía poco para pequeños cambios del flujo
y es constante para el caso particular de un flujo fuertemente turbulento.
En general el coeficiente de fricción f depende del número de Reynolds, de la rugosidad relativa de la tubería ǫ/D


y de la sección de la tubería,

tubería representadas por el diámetro de la misma D

Rh =
per.

∆hl =
=
ǫ
f = f (Re,
D

ρgπR4

∆hl =
, sección de la tubería)
donde la rugosidad relativa representa la relación entre la rugosidad de un elemento típico y las dimensiones de la

No sólo es importante la magnitud, sino también la forma de la sección de la canalización, ya que diferencias de la
forma circular dan lugar a flujos secundarios que aumentan la caída de presión. Cuando se trabaja con canalizaciones
no circulares se utiliza generalmente el radio hidráulico o el diámetro hidráulico
2A 2 × área de la sección
perímetro mojado

=
8µL¯

32νL¯
2
v
v
ρgR2

=
=
32µL¯
ρgD2
v
=

64 L v 2
¯
Dh = 2Rh =

32νL¯
gD2
v
4A
per.

Si tenemos en cuenta la definición del número de Reynolds, Re = v D/ν, esta expresión queda como
¯

=f
L v2
¯
(4–25)

en vez del radio o el diámetro, y se trata la canalización como si fuera de sección circular, con el número de Reynolds
calculado en función del diámetro hidráulico.

Régimen laminar.

Para tuberías de sección circular uniforme por las que circula un fluido newtoniano en régimen laminar, esta pérdida
de carga que aparece tanto en tuberías lisas como rugosas puede deducirse de la fórmula de Poiseuille. A partir de la
expresión para el caudal (4–21) y teniendo en cuenta que el diámetro es dos veces el radio, se tiene que
8µLQ
(4–26)

(4–27)
gD Re D 2g D 2g
v2
¯
donde f = 64/Re es el coeficiente de fricción (adimensional), L es la longitud, D el diámetro y la altura de
2g
velocidad. Esta expresión coincide con la ley empírica (4–24) de Darcy-Weiesbach
L v2
¯
∆h = f (4–24)
D 2g
que se utiliza para un flujo general, tomando para el caso particular del flujo laminar, la relación f = 64/Re, conocida
también como ecuación de Hagen-Poiseiuille.
Estas relaciones permiten obtener la tensión cortante en la pared de la tubería circular para un flujo laminar. De las
expresiones anteriores (4–19), se tiene
γR γR
τ0 = ∆hl = ∆hl
2L 2L
Al despejar ∆hl se obtiene
2τ0 L 4τ0 L
∆hl = =
γR γD
que, al comparar con la ecuación de Darcy-Weiesbach da la tensión en las paredes de la tubería
γ¯2
v ρ¯2
v
τ0 = f =f (4–28)
8g 8

Cuando las canalizaciones no son cilíndricas la relación entre el número de Reynolds y el coeficiente de fricción es,
para régimen laminar y fluido newtoniano, de la forma
k
f= con 48 ≤ k ≤ 96 (4–29)
Re
en función de un factor geométrico k. Algunos valores típicos del factor geométrico k son los siguientes


Régimen turbulento

viene dado por f = 64/Re.

1
√ = 0,869 ln √
f
Cuadrado
Rectángulo 2:1
Rectángulo 5:1
Láminas paralelas

f = 0,316/Re0,25

Re
f
− 0,8 = 2 log √
56,91
62,19
76,28
76,28

La ley experimental de Darcy-Weiesbach es válida para todo tipo de flujo con el coeficiente de fricción adecuado.
Como acabamos de ver, para un fluido newtoniano en régimen laminar y tuberías circulares, el coeficiente de fricción

Tuberías lisas. En régimen turbulento (Re >2300) y tuberías lisas se utilizan diversas expresiones f (Re) depen-
diendo del rango del número de Reynolds. Para números de Reynolds en el rango (3000<Re<100000) se utiliza la
fórmula de Blausius

Re
f
− 0,8
(4–30)

mientras que para números de Reynolds superiores a éstos, pero menores que 3 × 106 se utiliza la fórmula de von
Karman
(4–31)

Tuberías rugosas En régimen turbulento y tuberías rugosas se utilizan distintas funciones del tipo f = f (ǫ/D, Re)
donde ǫ es el coeficiente de rugosidad y d el diámetro de la conducción.
Para régimen turbulento en tuberías circulares la más utilizada de todas es la fórmula de Colebrook (1939)

1 ǫ/D 2,523 ǫ/D 2,523
√ = −0,869 ln + √ = −2 log + √ . (4–32)
f 3,7 Re f 3,7 Re f

Ésta es una ecuación implícita que converge bien en pocas iteraciones cuando se aborda mediante métodos numéricos.
El diagrama de Moody que se presenta en la figura 4–9 es una representación logarítmica de la ecuación de Colebrook
en ejes coeficiente de fricción/número de Reynolds para distintos valores de la rugosidad relativa, combinado con la
relación f = 64/Re para flujo laminar.
En el intervalo de rugosidades relativas 10−6 ≤ ǫ/D ≤ 10−2 y de números de Reynolds 5000 ≤ Re ≤ 108 se
puede hacer la aproximación explícita

1, 325 0,25
f= 2 = 2 (4–33)
ǫ/D 5,74 ǫ/D 5,74
ln + 0,9 log + 0,9
3,7 Re 3,7 Re

que da valores aproximados de la ecuación de Colebrook en buena parte del rango de interés.
Para rugosidades relativas altas y números de Reynolds altos, f depende fundamentalmente de la rugosidad relativa
y poco del número de Reynolds, como se observa en el diagrama de Moody y puede ajustarse por

1/ f = 1,14 − 0,869 ln(ǫ/D) = 1,14 − 2 log(ǫ/D). (4–34)


0,08

Agustin Martin Domingo 0,05
0,04
0,03
0,02
0,015

Rugosidad relativa ǫ/D
0,04
Coeficiente de fricción f

0,01
0,008
0,006
0,004

0,002
0,02 0,001
0,0008
0,0006
0,0004
0,0002
0,0001
0,00 0,00005
0,01 001
103 104 105 106 107
Número de Reynolds Re

Figura 4–9 Diagrama de Moody para tuberías de sección circular obtenido a partir de la relación f = 64/Re para régimen laminar y la ecuación
de Colebrook (4–32).

Cuando se conocen las características del flujo y de la tubería es sencillo ∆h
obtener el coeficiente de fricción f y la pérdida de carga directamente desde
v
el diagrama de Moody. Cuando la situación es la contraria, es decir se cono-
cen las características de la tubería y de alguna forma la pérdida de carga (p.ej, f f
0
midiendo las presiones entre dos puntos si se conoce el diámetro de la tubería)
y se quiere obtener el caudal, no es posible hacer el proceso inverso salvo que Re
se conozca f para esas características del flujo (f depende tanto de las carac-
terísticas de la tubería como de las características del flujo). Para resolver el D. Moody
problema, se utilizan métodos iterativos. Se parte de un valor de f estimado (p.ej, si no se tiene más información, algo
en torno a 0,02), se calcula la velocidad a partir de la ecuación de Darcy-Weiesbach. De aquí y de las características
de la tubería se obtiene el número de Reynolds y con la rugosidad relativa se extrae del diagrama de Moody un nuevo
valor del coeficiente de fricción. Repitiendo el proceso hasta que converja suficientemente se tiene un buen valor de f
y el caudal como Q = vS.

4.6.2. Pérdidas de carga singulares.

Representan disipaciones adicionales de energía en el flujo, normalmente causadas por los flujos secundarios (turbu-
lencias) inducidos en cambios de sección o de dirección de la tubería. Aunque en algunos sistemas constituyen una
parte importante de la pérdida de carga total, se les suele denominar denominar pérdidas secundarias o menores.
Corresponden a:


•
•
•
•
•
•

Secciones de rejilla.
Ensanchamientos.
Estrechamientos.
Codos.
Válvulas.
Bifurcaciones.
Para obtener las pérdidas de carga locales o singulares, en general se utilizan expresiones experimentales de la forma

∆h = K
v2
2g

donde K es un coeficiente que se obtiene a veces de forma teórica y a veces de forma experimental y que se encuentra
en los manuales de hidráulica. En el caso concreto de los ensanchamientos esta expresión aparece en función de (∆v)2
(que en el fondo es similar a 4–35).
No hay que olvidar que, aunque el coeficiente K parece constante, en general tiene alguna dependencia adicional
de las condiciones del flujo. En particular, los siguientes factores afectan a K

• La geometría exacta del componente en cuestión.

• El número de Reynolds.

• Proximidad a otros elementos, etc... (Los valores tabulados son para componentes “aislados”, con una cierta
longitud de tubería entre el elemento analizado y el siguiente).

Por ejemplo, para los ensanchamientos graduales la dependencia de K con las secciones inicial y final del tubo es tal
que la pérdida de carga viene dada por†
(4–35)

(∆v)2
∆h = K (4–36)
2g
donde en K aún queda incluida una dependencia en el ángulo del cono.
A continuación se da un resumen de algunos de los coeficientes de pérdidas más habituales. En todas las expresiones
que siguen se ha utilizado el subíndice 1 para representar los valores aguas arriba y 2 para representar la situación aguas
abajo. Para enfatizar a que parte corresponden las velocidades se ha intentado utilizar las mayúsculas para los valores
correspondientes a la sección mayor y las minúsculas para los valores correspondientes a la sección menor.

1. De depósito a tubería (pérdida a la entrada)
Conexión a ras de la pared tubería entrante conexión abocinada
v2 v2 v2
0,50 2 1,00 2 0,05 2
2g 2g 2g
2
v1
2. De tubería a depósito (pérdida a la salida) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1,00
2g

(v1 − V2 )2
3. Ensanchamiento brusco . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2g

(v1 − V2 )2
4. Ensanchamiento gradual . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . K
2g
† Esto no es más que un caso particular de 4–35, ya que por la ecuación de continuidad, v1 S1 = v2 S2 y por tanto,
2 2
(∆v)2 (v2 − v1 )2 S2 v2
∆h = K =K =K 1−
2g 2g S1 2g


ángulo | d1 /d2
4◦
10◦
15◦
20◦
30◦
50◦
60◦

D1 /d2
K
1,2
0,08
1,2
0,02
0,04
0,09
0,16
0,25
0,35
0,37

1,4
0,17
1,4
0,03
0,06
0,12
0,23
0,36
0,50
0,53

1,6
0,26
1,6
0,03
0,07
0,14
0,26
0,42
0,57
0,61

1,8
0,34
1,8
0,04
0,07
0,15
0,28
0,44
0,61
0,65

2,0
0,37

codo a 45◦ . . . . . . . . . . . . . . . . . .
codo a 90◦ . . . . . . . . . . . . . . . . . .
bifurcaciones en T . . . . . . . . . . .
válvulas de compuerta (abierta)
2,0
0,04
0,07
0,16
0,29
0,46
0,63
0,68

2,5
0,41
2,5
0,04
0,08
0,16
0,30
0,48
0,65
0,70

3,0
0,43

0,35 − 0,45
0,50 − 0,75
1,50 − 2,00
≃ 0,25
3,0
0,04
0,08
0,16
0,31
0,48
0,66
0,71

4,0
0,45
4,0
0,04
0,08
0,16
0,31
0,49
0,67
0,72

5. Venturímetros, boquillas y orificios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

5,0
0,46
5,0
0,04
0,08
0,16
0,31
0,50
0,67
0,72

1
cr
−1

6. Contracción brusca . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . K

7. Codos, accesorios, válvulas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . K
2
v2
2g
2
v2
2g

v2
2g

válvulas de control (abierta) . . . ≃ 3,0

8. Salida de tubería a la atmósfera (desagüe libre a la atmósfera).
La pérdida en el desagüe libre de una tubería a la atmósfera es generalmente muy pequeña al no producirse
apenas turbulencias en la tubería en ese proceso.

4.6.3. Variación de altura debida a una máquina.

Las máquinas hidráulicas pueden

• Aportar energía al flujo (bombas).

• Extraer energía del flujo (turbinas).

El aporte de energía que realiza una bomba al flujo se representa habitualmente a través de la energía aportada por
unidad de peso del fluido, a la que se denomina altura de la bomba o altura aportada por la bomba HB . Veamos ahora
la relación entre la potencia de la bomba y la altura de la misma. Para ello consideraremos una bomba colocada en
una tubería horizontal de sección constante y estudiaremos el aumento de energía en un elemento dm del fluido que
atraviesa la bomba en un intervalo de tiempo dt. Si denominamos h1 a la energía por unidad de peso a la entrada de
la bomba y h2 a la energía por unidad de peso a la salida de la bomba, la energía de un elemento de masa dm antes
de entrar en la bomba será gh1 dm mientras que la energía del elemento de masa dm al salir de la bomba será gh2 dm,
por lo que la variación de energía del fluido en el intervalo de tiempo dt es g(h2 − h1 )dm = gHB dm donde HB es
la altura de la bomba o energía por unidad de peso aportada por la bomba. La potencia de la bomba será entonces, en
función de la altura de la misma
gHB dm dm
P = = gHB = gHB Qm = ρgQHB (4–37)
dt dt


A
Bomba

pA + ρgHbomba = pB
1

B
Figura 4–10
2

A
h2

h1

Turbina

pA − ρgHturbina = pB
B

Figura 4–11 Expresiones para incluir correctamente el efecto de las bombas y de las turbinas.

La misión de la turbina es extraer energía del fluido, disminuyendo su energía. Análogamente al caso de las bombas,
se denomina altura HT de la turbina a la energía que ésta extrae por unidad de peso del fluido. De forma semejante al
caso de la bomba, la potencia extraída por la turbina está relacionada con la altura HT de ésta mediante
Ptur = ρgHT Q = gHT Qm (4–38)
Para considerar correctamente los signos de las alturas de bomba y turbina, es conveniente seguir el sentido del flujo,
añadiendo a la altura total inicial todas las alturas aportadas por bombas y restándole todas las alturas de pérdidas y
todas las alturas extraídas por turbinas, llegándose a la ecuación de Bernoulli generalizada.
p1 v2 p2 v2
+ 1 + h1 + HB − HT − ∆htot = + 2 + h2 (4–39)
ρg 2g ρg 2g

4.6.4. Representación gráfica de la pérdida de carga.

Consideremos el caso de una tubería alimentada por un gran depósito cuyo nivel se mantiene constante, que descarga
en otro gran depósito cuyo nivel se mantiene también constante. Las líneas con que se trabaja habitualmente son la
línea de alturas totales y la línea de alturas piezométricas.

1. La línea de altura total en cada punto es el resultado de restar a la del depósito (nivel de energía del depósito)
todas las pérdidas de carga encontradas en el camino (en este caso todas las pérdidas de carga a la izquierda del
punto) y sumar todas las ganancias añadidas por las bombas.
v2
2. La línea de altura piezométrica se obtiene restando la altura de velocidad de la línea de altura total
2g


4.7. Redes de tuberías

para su cálculo.
Htotal

Hpiezo
∆Ht =

Figura 4–12 Línea de alturas totales y línea de alturas piezométricas.
nivel

i

Analizaremos ahora algunos ejemplos sencillos de redes de tuberías y los elementos fundamentales a tener en cuenta
∆hi

El caso más sencillo es el de tuberías en serie. En este caso, en régimen permanente el caudal másico Qm que circula
a través de cada tubería es el mismo, Qm1 = Qm2 ... = Qmn , mientras que la pérdida de carga total a lo largo del
sistema de tuberías en serie se obtiene como la suma de las pérdidas de carga en cada uno de los tramos (más las
pérdidas en los elementos de unión de las tuberías),

∆htot = ∆hi
i

En las tuberías en paralelo, la corriente que originalmente circulaba por una tubería se distribuye en un conjunto de
tuberías para posteriormente reunirse de nuevo en una sola. En este caso el caudal en la línea principal es la suma de
los caudales en las tuberías en que se divide, Qprincipal = Qm1 + ... + Qmn . Como los valores de la altura de carga
en los puntos de bifurcación deben ser los mismos para todas las tuberías colocadas en paralelo, las pérdidas de carga
en todas ellas deben ser iguales.
Cuando tenemos redes de tuberías más complejas hay dos reglas fundamentales que deben tenerse en cuenta:

Q1 B

Q3
Q1
Q2

A C
Q2
Q3

Figura 4–13 Reglas fundamentales para trabajar con redes de tuberías

1. En cualquier nodo (unión de dos o más tuberías) debe conservarse la masa. Es decir, el caudal másico neto debe
ser cero en régimen estacionario, todo lo que entra en un intervalo de tiempo debe salir en el mismo intervalo,

4.8. Cavitación 65

no puede producirse acumulación de materia. Para el caso de un fluido incompresible, ésto es cierto también
para el caudal.

en el bucle.

4.8. Cavitación
i
Q mi = 0
i

2. No debe de haber pérdida de carga neta a lo largo de un bucle. Si comenzamos en una posición del bucle y
sumamos las variaciones de carga hidráulica, cada una con su correspondiente signo, al llegar a la posición
original la carga hidráulica debe ser la misma que al principio,

∆hrozi + ∆Hbombas + ∆Hturbinas = 0

Se denomina cavitación a la formación de cavidades en el líquido en movimiento llenas de vapor. Se produce cuando
la presión en algún lugar del fluido se hace menor que la presión de vapor saturante del líquido, es decir la presión que
corresponde a la ebullición del líquido a una temperatura dada.
Por ejemplo, acelerando agua desde el reposo hasta velocidades de 15 m/s la presión desciende en aproximadamente
una atmósfera y produciéndose cavitación. De hecho esto ocurre a menudo en los elementos en los que se producen
pérdidas locales cuando se produce un gran aumento de velocidad al estrecharse el flujo después de dilatarse.
Así, en las partes de baja presión puede comenzar una intensa evaporación de líquido, formándose las cavidades de
cavitación formadas por burbujas de vapor. Al aumentar posteriormente la presión el vapor se condensa, normalmente
con un brusco chasquido y las cavidades de cavitación desaparecen.
La cavitación va acompañada de un notable aumento de las pérdidas locales y se caracteriza a través un parámetro
adimensional denominado número de cavitación
2(p − pv )
Ca = (4–40)
ρv 2
donde p y v son la presión y la velocidad antes de la resistencia local pv la presión de vapor saturante del agua.
Para que en el elemento local no se produzca cavitación el número de cavitación debe ser menor que un valor crítico
que se determina de forma experimental, Ca < Ca cr. .
También se produce cavitación en hélices, bombas y turbinas. La cavitación produce grandes oscilaciones, con
aumento de las pérdidas de carga y además ayuda a la corrosión de los metales. Ésto último es consecuencia de los
procesos casi instantáneos de contracción de las burbujas que producen grandes presiones que inicialmente desmigajan
lentamente el material y posteriormente extienden su efecto al interior del mismo, combinados a menudo con efectos
químicos. Como consecuencia de todo ésto el metal se corroe y llega a adquirir una estructura esponjosa.

4.9. Golpe de ariete o choque hidráulico

Se denomina golpe de ariete o choque hidráulico al brusco aumento de la presión que se produce en las tuberías
cuando se detiene bruscamente el líquido que circula por ellas. Ésto ocurre, por ejemplo, cuando se cierran rápidamente
alguno de los distintos dispositivos (válvulas, grifos) presentes en la tubería o cuando se produce una parada imprevista
de las bombas que impelen el líquido.
El golpe de ariete es especialmente peligroso en las tuberías largas por las que circulan grandes cantidades de líquido
a altas velocidades. En estos casos, el golpe de ariete puede producir graves daños en las juntas de unión, en la propia
tubería o en las máquinas.
Para estudiar este fenómeno consideremos una tubería horizontal de sección circular constante S y radio R por la que
circula un líquido de velocidad v constante. Si se cierra la válvula de forma idealmente instantánea la capa de líquido


Agustin Martin Domingo v
p0

dl
p0 + ∆p

Figura 4–14 Golpe de ariete. El líquido que circula por la tubería se detiene súbitamente.

en contacto con esta se detendrá de forma inmediata produciéndose un aumento de la presión como consecuencia de
la transformación de energía cinética en energía potencial de presión.
Consideremos primero un caso extremo, el caso de un fluido completamente incompresible que circula por una
tubería totalmente rígida. En este caso, si la válvula se cierra en un tiempo ∆t todo el fluido se detiene simultáneamente
en ese tiempo la bajo la acción de la fuerza ∆pπR2 en la superficie en contacto con la válvula y la cantidad de
movimiento de éste cambia de ρLπR2 v a 0 en el intervalo ∆t. El impulso de la fuerza mencionada debe ser igual a la
variación de la cantidad de movimiento del líquido,

∆pπR2 ∆t = ρLπR2

Esta ecuación da un buena estimación del aumento de presión asociado al cierre lento de una válvula. Sin embargo
da un resultado pequeño y no demasiado interesante. También hay que recalcar que da como resultado una presión
infinita si el cierre de la válvula se produce de forma instantánea, pero para presiones muy altas las aproximaciones de
líquido totalmente incompresible y tubería totalmente rígida son malas.
Cuando la válvula se cierra de forma instantánea el fluido se frenará rápidamente y la presión aumentará, pero este
aumento de presión producirá una pequeña compresión en el fluido y una pequeña expansión de la tubería (que espera-
mos elástica). Así, el líquido en contacto con la válvula se detendrá instantáneamente, pero al ser éste compresible, el
resto del líquido en la tubería no se detendrá de forma instantánea. El límite del volumen que incluye el líquido parado
se desplaza alejándose de la válvula con una cierta velocidad c a la que denominamos velocidad de propagación de la
onda de presión.
Examinemos el elemento de volumen del líquido Scdt = Sdl junto a la válvula. En un tiempo dt este volumen se
ha detenido y su cantidad de movimiento ha pasado de ρSdl v a 0. El impulso de la fuerza que actúa en el tiempo dt
deteniendo el fluido es ∆pSdt. Al igualarlo a la variación de la cantidad de movimiento del elemento en ese tiempo
dt tenemos
∆pSdt = ρSdl v
Como c = dl/dt se obtiene el valor del aumento de presión en el golpe de ariete:

∆p = ρcv (4–41)

Después se detiene la siguiente capa y así sucesivamente, y la presión elevada que inicialmente surge en las proxi-
midades de la válvula se propaga poco a poco por toda la tubería con una velocidad c, en contra del flujo original del
líquido.
Si la presión al principio de la tubería no cambia (por ejemplo, cuando el líquido viene de un depósito abierto a la
atmósfera) una vez que la onda de presión alcanza la sección inicial de la tubería, ésta comienza a propagarse en el
sentido contrario con la misma velocidad c, ahora como onda de disminución de presión. Cuando esta onda alcanza la
válvula la presión cae y se hace menor que antes del cierre y de nuevo una onda de presión sale en sentido contrario.
Así, se alternan ciclos de aumento y disminución de la presión en intervalos de tiempo iguales al tiempo empleado por
la onda en recorrer completamente en los dos sentidos el tramo de tubería entre la válvula y el principio.

4.9. Golpe de ariete o choque hidráulico 67

Así, el líquido realiza un movimiento oscilatorio que se va amortiguando como consecuencia de la disipación de
energía que tiene lugar en el proceso.
La velocidad de propagación de la onda de presión depende de las características del líquido y de la tubería y viene
dada por la expresión:†

la tubería.
c=
K
ρ
1+
KD
E∆R
−1

donde K es el módulo de compresibilidad del líquido (inversa del coeﬁciente de compresibilidad), ρ la densidad del
líquido, E el módulo de Young del material de la tubería, D el diámetro de la misma y ∆R el espesor de la pared de

Si se considera totalmente la tubería totalmente rígida (E = ∞) la velocidad de propagación de la onda de presión
queda en la forma

igual a la velocidad de propagación del sonido en un líquido.
c=
K
ρ
,

Para evitar los efectos del golpe de ariete se utilizan diversos sistemas:

1. Dispositivos que aumentan el tiempo de cierre de todo tipo de válvulas aminorando así el efecto.
2. Uso de un diámetro grande en la tubería para que las velocidades sean bajas.
3. Dispositivos de seguridad en las bombas para que en caso de corte de energía reduzca lentamente la velocidad
del motor.
4. Colocación de válvulas de seguridad automáticas, campanas de aire y chimeneas de gran sección.
(4–42)

(4–43)

† que se obtiene de la condición de igualdad entre la energía cinética del líquido que se mueve por la tubería y la suma de los trabajos de

compresión del líquido y alargamiento del tubo.

Capítulo 5

Mecánica de fluidos en medios permeables.

5.1. Fenómenos de superficie: tensión superficial y capilaridad . . . .
5.1.1. Tensión superficial. . . . . . . . . . . . . . . . . . . . . . .
5.1.2. Ángulo de contacto . . . . . . . . . . . . . . . . . . . . . .
5.1.3. Capilaridad . . . . . . . . . . . . . . . . . . . . . . . . . .
5.2. Fluidos en medios permeables. . . . . . . . . . . . . . . . . . . .
5.2.1. Velocidad de descarga y velocidad microscópica promedio. .
5.2.2. La ley de Darcy en una dimensión. . . . . . . . . . . . . . .
5.2.3. Potencial hidráulico y sus componentes. El campo hidráulico.
5.2.4. La ley de Darcy en tres dimensiones. . . . . . . . . . . . . .
5.2.5. Flujo a través de terrenos saturados homogéneos e isótropos.
5.2.6. La ecuación de Laplace y su resolución. . . . . . . . . . . .
5.2.7. Flujo permeable a través de medios anisótropos . . . . . . .
5.2.8. La red de flujo. . . . . . . . . . . . . . . . . . . . . . . . .
5.2.9. Ejemplos de redes de flujo en distintos sistemas . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
69
69
73
74
76
78
79
83
84
85
85
87
90
90
5.2.10. Subpresiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92

5.1. Fenómenos de superficie:
tensión superficial y capilaridad.

5.1.1. Tensión superficial.

Las superficies libres de los líquidos presentan una serie de fenómenos cuyo origen se reduce a la tendencia de la
superficie a permanecer lo más pequeña posible. Este hecho, conectado con la presencia de una tensión en la película
superficial de forma análoga a la tensión en una membrana elástica fina que se encuentra tensa puede observarse de
forma sorprendente mediante una serie de experimentos simples.

Algunos ejemplos de tensión superficial.

Consideremos un bastidor como el de la figura 5–1, que se sumerge en una solución jabonosa, quedando al sacarlo
un fina película de la misma. Si sobre esta película depositamos un hilo cerrado en un lazo, se observa que éste no
adopta una posición especial, sino la posición que haya tomado al caer. Sin embargo, si se pincha la parte de la película
dentro del hilo, éste adopta la forma de una circunferencia por la acción de la tensión superficial. De esta forma, la
película de fluido ha alcanzado su mínima superficie.
Otro efecto de la tensión superficial es el siguiente: ¿Qué ocurre cuando un objeto de mayor densidad que el agua se
deposita sobre la misma? La respuesta más inmediata sería que éste se hundiría. Sin embargo, ésto no siempre es así.
Consideremos el caso de una aguja sobre el agua. Si ésta se deposita con el suficiente cuidado es posible hacerla flotar
sobre la superficie del líquido. Esto es debido a que las fuerzas que soportan la aguja son debidas no sólo al empuje
del agua sobre la misma, que es obviamente menor que su peso, sino también a la tensión superficial.

69

70 Capítulo 5. Mecánica de fluidos en medios permeables.

Agustin Martin Domingo Figura 5–1 Efecto de la tensión superficial en una película delgada.

Otro efecto de la tensión superficial es el hecho de que un líquido que no se encuentre sometido a fuerzas externas
tiende a tener una forma esférica, minimizando de esta forma su superficie para un volumen dado.
La tensión superficial es utilizada por algunos insectos como el zapatero (gerris lacustris) que se muestra en la
figura 5–2. Éste se mantiene muy por encima del agua gracias a la tensión superficial en sus patas, que utiliza para

Figura 5–2 El zapatero se desliza sobre la superficie del agua gracias a la tensión superficial. Fotografía de ISAKA Yoji bajo licencia CC-BY-SA.

5.1. Fenómenos de superficie: tensión superficial y capilaridad 71

Tabla 5–1 Valores de la tensión superficial medidos en dy/cm. En todos los casos se asume que el líquido está en contacto con aire.
(1dy/cm=10−3 Nm−1 ). Salvo indicación en sentido contrario, los datos están tomados a 20◦ C

patas.
Sustancia

Acetona
Metanol
Etanol
Aceite de oliva
Aceite de parafina
Benceno
Glicerina

Principios físicos de la tensión superficial.
tensión
superficial
23,7
22,6
22,3
33
26
28,9
64
Sustancia

Agua a 5◦ C
Agua a 10◦ C
Agua a 15◦ C
Agua a 20◦ C
Agua a 25◦ C
Agua a 30◦ C
Mercurio
tensión
superficial
74,92
74,22
73,49
72,75
71,97
71,18
475

remar (las centrales) o a modo de timón (las posteriores), mientras que caza a sus presas con las delanteras, pero todas
ellas contribuyen a que se mantenga por encima del agua. Obsérvese la curvatura de la superficie del agua bajo las

Cualitativamente ésto tiene su origen en la fuerza de atracción molecular no compensada en las moléculas de la
superficie. Mientras que en el interior de un líquido, cada molécula está rodeada de otras moléculas por todos lados,
y por tanto las fuerzas intermoleculares se compensan unas a otras, en la superficie no hay moléculas por encima de
las moléculas superficiales. Así, las moléculas de la superficie experimentan una fuerza neta dirigida hacia el interior
de modo que éstas no escapan. De este modo la superficie intenta hacerse lo menor posible. Asimismo, si se eleva
ligeramente una molécula superficial, los enlaces moleculares con las moléculas adyacentes se alargan, produciéndose
una fuerza restauradora que tira de la molécula desplazada de nuevo hacia la superficie.
Esto es lo que ocurre cuando se coloca cuidadosamente una aguja sobre la superficie. Las moléculas superficiales
descienden ligeramente, separándose de las demás moléculas de la superficie, que ejercen una fuerza atractiva sobre
las mismas, soportando la aguja. Así, la superficie de un líquido presenta una cierta semejanza con un membrana
elástica tensa.
Este concepto de tensión en la superficie de un fluido, conocido como tensión superficial fue introducido por J. A.
von Segner como la fuerza que actúa en la superficie de un fluido a cada lado de una línea de longitud unidad dibujada
en la superficie.
La fuerza necesaria para romper la superficie se puede medir levantando un alambre fino. La fuerza necesaria para
romper la superficie resulta ser proporcional a la longitud de la superficie que se rompe, que es el doble de la longitud
del alambre, si hay una película de superficie en ambos lados del mismo, o a la longitud del mismo, si sólo hay
película superficial en un lado. Si el alambre posee una masa m y una longitud l, y es necesario ejercer una fuerza F
para separarlo de la superficie (de dos caras), el coeficiente de tensión superficial σ es:
F − mg
σ= (5–1)
2l
Algunos valores de éste coeficiente de tensión superficial se dan en la tabla 5–1.

Relación entre la tensión superficial y la diferencia de presiones a través de una superficie.

Estas no son las únicas consecuencias de la tensión superficial. Hasta el momento se ha dicho que, una vez alcanzado
el estado de equilibrio, las presiones de dos cuerpos en contacto han de ser iguales. Sin embargo, esta afirmación es
cierta únicamente cuando el efecto de los fenómenos de superficie es despreciable. Cuando se tiene en cuenta la tensión
superficial, las presiones de dos medios en contacto pueden ser diferentes.


σdl2
σdl1

dl2
r2

r1
dα1
dα2

σdl2

dl1
σdl1
α/2

r

α/2
α/2

Figura 5–3 Fuerza ejercida por la tensión superficial sobre un elemento infinitesimal rectangular de una superficie curva.

Como ejemplo de lo que se acaba de afirmar, consideremos el caso de una gota de líquido en el aire. La tendencia
de ésta a disminuir su superficie conduce a una compresión de la gota y a un aumento de la presión interna, siendo la
presión del líquido en la gota mayor que la del aire circundante. A la diferencia entre ambas presiones se le denomina
diferencia, variación de presión o presión diferencial y la designaremos por pdif .
Para calcular esta magnitud, debemos tener en cuenta que el trabajo realizado por las fuerzas superficiales al dismi-
nuir la superficie de la gota en la magnitud dS viene dado por la disminución de la energía superficial σdS. Este trabajo
es realizado por las moléculas de la superficie sobre el líquido del interior de la gota, y puede expresarse también como
el trabajo volumétrico realizado sobre el líquido, en función de la variación de volumen dV de la gota, en la forma
(pgota − pext )dV = pdif dV . Ambas expresiones se refieren al mismo trabajo, en un caso realizado sobre el líquido
de la gota y en el otro caso realizado por el líquido de la gota, y por tanto son iguales:

σdS = pdif dV (5–2)

Veamos ahora que ocurre con una gota esférica de radio r, con una superficie S = 4πr2 y un volumen V = 4πr3 /3.
Sustituyendo ésto en la expresión anterior, se tiene:

σ8πrdr = pdif 4πr2 dr

con lo que tenemos que la presión diferencial en el caso de una gota esférica viene dada por:
2σ
pdif = (5–3)
r
Un cálculo análogo puede realizarse para una superficie cilíndrica, como la que aparece cuando entre dos planos
verticales próximos hay agua. En este caso, el volumen es V = πr2 L y la superficie S = 2πrL, con L la altura del
cilindro. Así, para el caso de una superficie cilíndrica, la presión diferencial es:
σ
pdif = (5–4)
r
Para obtener la relación entre la tensión superficial σ y la fuerza resultante perpendicular a la superficie consideremos
el elemento diferencial de superficie rectangular que se observa en la figura 5–3. Si σ es la tensión superficial, r1
y r2 los radios de curvatura de la superficie en las direcciones indicadas y dl1 y dl2 las longitudes de los lados del
rectángulo, actuarán unas fuerzas σdl1 y σdl2 en los lados del elemento de superficie. Si dα1 y dα2 son los ángulos
subtendidos por los lados dl1 y dl2 respectivamente, dα1 /2 y dα2 /2 serán los ángulos que forman las fuerzas σdl1 y


σdl2 con el plano tangente a la superficie, y se cumplirá que dl1 = r1 dα1 y dl2 = r2 dα2 . Por tanto, la componente
perpendicular de la fuerza de tensión superficial que actúa sobre el lado 1 viene dada por

2σdl1 sen dα2 /2 ≈ 2σdl1 dα2 /2 = σdl1

y la componente perpendicular de la fuerza que actúa sobre el lado 2 por

2σdl2 sen dα1 /2 ≈ 2σdl2 dα1 /2 = σdl2

surgida de la diferencia de presiones entre los dos medios de la forma

∆F = ∆pdl1 dl2

pdif = σ
1
r1
+
1
r2
dl2
r2

dl1
r1
Para que se esté en condiciones de equilibrio mecánico, la suma de estas fuerzas debe ser compensada por una fuerza

La diferencia de presiones entre las dos caras de una superficie líquida es proporcional a la suma de las inversas
de los radios de curvatura en dos direcciones mutuamente ortogonales, siendo la presión mayor siempre hacia la
parte cóncava de la superficie. Hay que recalcar que la suma 1/r1 + 1/r2 es un invariante respecto de la rotación de
(5–5)

(5–6)

(5–7)

Así, en el caso general de una superficie de radios principales de curvatura (mínimo y máximo) r1 y r2 , la diferencia
de presiones viene dada por la ecuación de Laplace (o Laplace-Young).

(5–8)

coordenadas, de forma que la ecuación (5–8) es independiente de la dirección en la que se toman los ejes, obviamente
de acuerdo con el punto de vista físico.
La ecuación (5–8) se puede particularizar fácilmente al caso cilíndrico con r1 = ∞ y al caso esférico con r1 = r2

5.1.2. Ángulo de contacto

Se denominan fuerzas de cohesión a las fuerzas de atracción entre una molécula del líquido y otras moléculas del
mismo. La fuerza que tiene lugar entre una molécula del líquido y otra sustancia, como la pared de un tubo delgado,
es una fuerza adhesiva.
Cuando la fuerza adhesiva es menor que la fuerza de cohesión el ángulo se hace grande, mientras que cuando la
fuerza de cohesión es menor que la fuerza adhesiva, el ángulo se hace pequeño.
Estudiemos ahora ésto en función de las tensiones superficiales implicadas en el problema. Por el borde de la
superficie de un líquido que se halla en un recipiente concurren tres medios: (1) sólido, que es la pared del recipiente,
(2) líquido, y (3) gaseoso.
En la línea de contacto están aplicadas tres fuerzas de tensión superficial, estando cada una de ellas dirigida según la
tangente hacia el interior de la superficie de contacto de los dos medios correspondientes. El ángulo entre la superficie
del líquido y la pared sólida se denomina ángulo de contacto θ.
La superficie del líquido adquiere una forma tal que la resultante de las tres fuerzas σSL , σSG y σLG no tenga
componente a lo largo de la pared del recipiente (la reacción de la pared elimina la componente perpendicular a la
misma), quedando la condición de equilibrio:

σSG = σSL + σLG cos θ

donde no se ha incluido la longitud de la línea de contacto, que es la misma en todos los casos. Así, el coseno del
ángulo de contacto queda:
σSG − σSL
cos θ = (5–9)
σLG
Como se ve, el ángulo de contacto depende exclusivamente de la naturaleza de los tres medios en contacto (de las


S
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
σSG

σSL
θ

σLG
G

θ
L

1234567890123456789012345678901212345678901234567
1234567890123456789012345678901212345678901234567
S
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567
1234567

Figura 5–4 Equilibrio de fuerzas en la línea de contacto sólido-líquido-gas

θ
σSG

σSL
σLG
G

θ

L

12345678901234567890123456789012123456789012345678
12345678901234567890123456789012123456789012345678
1234567890123456789012345678901212345678901234567
1234567890123456789012345678901212345678901234567 12345678901234567890123456789012123456789012345678
12345678901234567890123456789012123456789012345678
1234567890123456789012345678901212345678901234567
1234567890123456789012345678901212345678901234567 12345678901234567890123456789012123456789012345678
12345678901234567890123456789012123456789012345678
1234567890123456789012345678901212345678901234567
1234567890123456789012345678901212345678901234567 12345678901234567890123456789012123456789012345678
12345678901234567890123456789012123456789012345678
1234567890123456789012345678901212345678901234567
1234567890123456789012345678901212345678901234567 12345678901234567890123456789012123456789012345678
12345678901234567890123456789012123456789012345678

σSG > σSL σSG < σSL
cos θ > 0◦ cos θ < 0◦
θ < 90◦ θ > 90◦
Figura 5–5 Distintos comportamientos de una gota sobre una superficie.

tensiones superficiales en los límites) y no depende ni de la forma del recipiente ni del valor de g. Estas tensiones
superficiales y las fuerzas asociadas son muy sensibles al estado de las superficies divisorias y a sus purezas.
Si σSG > σSL , es decir si la tensión superficial sólido-gas es mayor que la tensión superficial sólido-líquido, se
tiene que cos θ > 0 y que θ < 90◦ . El ángulo es agudo y el menisco tiene la forma cóncava de la parte izquierda de
las figuras 5–4 y 5–5. En este caso se dice que el líquido “moja” a la superficie sólida, una gota del líquido depositada
sobre la superficie se extiende sobre la misma.
Cuando σSG < σSL , la tensión superficial sólido-gas es menor que la tensión superficial sólido-líquido, teniéndose
que en este caso cos θ < 0◦ y que θ > 90◦ . Este caso corresponde al menisco convexo, y una gota depositada sobre la
superficie se recoge sobre sí misma, como se muestra en la parte derecha de las figuras 5–4 y 5–5. En este caso se dice
que el líquido “no moja” la superficie del sólido. Ejemplos de ángulos de este tipo se tienen en los casos de mercurio
sobre vidrio (150◦ ) y de agua sobre parafina (105◦). Las gotas de estos líquidos se elevan intentando disminuir el área
de contacto con la superficie.

5.1.3. Capilaridad

La elevación de un líquido por un tubo delgado, llamado tubo capilar, es otro fenómeno de superficie íntimamente
ligado con la tensión superficial. En este proceso aparecen fuerzas tanto de cohesión como de adhesión. Si sumergimos


un tubo delgado de vidrio en agua, hemos visto que, ya que el agua moja bastante bien al vidrio, tendremos un ángulo
de contacto muy agudo, de forma que la superficie estará muy curvada si el tubo es lo suficientemente delgado. El
resultado de todo ésto es que se producirá una fuerza ascendente, subiendo el agua por el tubo hasta que el peso de la
columna de agua iguala esta fuerza.

r

h
θ

F

Figura 5–6 Capilaridad en un tubo fino

Consideremos el caso de un líquido que ha alcanzado una altura h en un tubo capilar de radio r. El tubo está abierto
en su parte superior y se encuentra, por tanto, a la presión atmosférica. La componente vertical de la fuerza que sostiene
al líquido es F · cos θ que, al ser la longitud de la línea de contacto l = 2πr, queda como σ2πr cos θ.
Si despreciamos la ligera curvatura de la superficie del líquido, el volumen del mismo en la columna es πr2 h
quedando, al hacer el equilibrio de las fuerzas:

σ2πr cos θ = ρghS = ρghπr2

dando para la altura de la columna capilar:

2σ cos θ
h= (5–10)
ρgr


5.2. Fluidos en medios permeables.

En los capítulos anteriores hemos estudiado el movimiento de fluidos “libres”. Sin embargo, los fluidos también
pueden moverse a través de ciertos materiales, como ocurre con las aguas subterráneas que se filtran a través de los
terrenos, o el petróleo que se mueve en las capas petrolíferas a los pozos petrolíferos. Este movimiento es consecuencia
de la acción conjunta de fuerzas gravitatorias, de presión, y otras de diverso origen, que dan lugar a un desplazamiento
que se produce siempre hacia las zonas de menor energía potencial.
Se denominan materiales o medios permeables a aquéllos que permiten el paso de fluidos a través de los mismos, ya
se trate de gases o de líquidos. Los materiales permeables pueden ser de muy distintos tipos, esponjosos como algunas
espumas, fibrosos, granulares como la arena de los terrenos, etc... En todos ellos la característica común es que el
flujo tiene lugar a través de una intrincada red de microcanales, con unas elevadas pérdidas de carga dada la pequeña
sección de éstos. En el caso de líquidos aparecerán también fuerzas de ascensión capilar junto a las fuerzas que dan
lugar al flujo.

Naturaleza de los terrenos como medios permeables.

Los distintos terrenos consisten en una gran cantidad de partículas pequeñas en contacto, con poros o espacios
huecos entre las mismas que pueden estar ocupados por agua, gas o vapor. El fluido puede circular de una parte a otra
del medio a través de multitud de microcanales entre las partículas que comunican los distintos poros, comportándose
como medios permeables. En el caso de un fluido en reposo en un medio permeable, éste asciende por el medio

Tabla 5–2 Valores típicos de la subida capilar para una serie de materiales
Material
Arena gruesa
Subida capilar
2 − 5 cm
Arena normal 12 − 35 cm
Arena fina 35 − 70 cm
Arcillas 2 − 4 m y superior

permeable por encima de la altura de la superficie libre del líquido de fuera del medio permeable. Esto es debido a
lo siguiente: El medio permeable está constituido por un complejo conjunto de microcanales, a través de los cuales
asciende el fluido por capilaridad. En la tabla 5–2 de dan algunos valores típicos de esta subida capilar para una serie
de materiales. A la hora de modelizar el comportamiento de este medio permeable frente a este efecto, se suele sustituir
en los cálculos la red real de capilares por un sistema de capilares verticales de un cierto diámetro con una determinada
densidad de capilares por unidad de superficie. De esta forma se puede tratar el problema, así como el problema que
surge cuando hay dos medios permeables de distintas características superpuestos.
La naturaleza del sistema de poros es variada dependiendo del tipo de terreno, aunque hay dos casos límites intere-
santes, el de las arenas y el de las arcillas.
Las arenas se componen de partículas macroscópicas de forma redondeada o angulosa. Estas partículas drenan
fácilmente, no absorben apenas agua, tienen un pequeño potencial capilar y cuando se secan no encogen de tamaño.
Por el contrario, las arcillas están compuestas de partículas microscópicas de forma laminar. Estas absorben agua con
facilidad, son bastante impermeables una vez mojadas, tienen un elevado potencial capilar y disminuyen de volumen
de forma considerable al secarse.
Las arenas se aproximan más al medio poroso ideal y son representativas del tipo de medio poroso que trataremos.
Aunque ésto parece imponer serias restricciones al estudio que vamos a hacer, en la mayoría de los problemas reales la
baja permeabilidad de las arcillas las hace ser relativamente impermeables en comparación con los terrenos de grano
más grueso.
Denominaremos nivel freático al lugar geométrico de los puntos que se encuentran a la presión atmosférica. Corres-
ponde al nivel que se alcanzaría en cada punto de la superficie del terreno si en el mismo se introdujera un piezómetro
la profundidad suficiente. Si tomamos como referencia de medida de presiones la presión atmosférica, la presión del

5.2. Fluidos en medios permeables. 77

terreno seco
nivel capilar

nivel freático

Figura 5–7 Nivel freático y capilar en un terreno. Por encima del nivel capilar no se encuentran granos de arena húmedos.

agua es positiva por debajo del nivel freático. Por encima del nivel freático pueden encontrarse granos de arena moja-
dos debido al efecto de la capilaridad hasta el nivel capilar. La presión del agua capilar es negativa y dará lugar al
denominado potencial de matriz.

Porosidad.

Asumamos que las partículas del medio poroso son todas de forma esférica uniforme. Si denominamos V al volumen
total y Vh al volumen de los huecos, se deﬁne la porosidad como

Vh
n= (5–11)
V
que corresponde a la fracción del volumen total ocupada por los poros, y el índice de poros como

Vh
e= (5–12)
V − Vh

que corresponde al cociente entre el espacio ocupado por los poros y el espacio ocupado por las partículas sólidas.
Ambos están relacionados por las expresiones
e n
n= e= (5–13)
1+e 1−n

Para un empaquetamiento cúbico de esferas de diámetro d el volumen total (por esfera) es V = d3 , mientras que el
volumen hueco es Vh = d3 − πd3 /6 quedando una porosidad de

π
n = 1− = 0,476.
6

Para un empaquetamiento romboédrico, que representa la forma más compacta de empaquetar una serie de esferas
uniformes, la porosidad es:
√
2
n=1− π = 0,26.
6

En la tabla 5–3 se dan valores para las porosidades de algunos terrenos naturales. Como hemos visto en los ejemplos
de empaquetamiento cúbico y romboédrico de esferas uniformes, incluso en un medio poroso ideal el espacio hueco
no es regular, sino que está compuesto por toda una serie de celdas cavernosas interconectadas por canales estrechos.


El medio poroso como un medio continuo.
Tabla 5–3 Porosidades de algunos terrenos naturales.
Medio
Arena uniforme, suelta
Arena uniforme, densa
Arcilla glacial blanda
Arcilla glacial dura
Arcilla orgánica blanda
Bentonita blanda
Porosidad
46
34
55
37
75
84

Los terrenos naturales contienen partículas que se desvían de forma considerable del comportamiento esférico (como
hemos visto por ejemplo en las arcillas) y que además están lejos de ser de tamaño uniforme. Por tanto la red de
microcanales es extremadamente complicada y de hecho es imposible en la práctica tratar de una forma razonable un
terreno atendiendo de forma detallada a la naturaleza microscópica de su estructura de dichos microcanales. Por esta
razón no fue posible tratar de forma sistemática el flujo permeable hasta el tratamiento de Darcy.
Afortunadamente, en los problemas de flujo a través de medios permeables no es necesario tener en cuenta el
flujo a través de los canales individuales. Nuestro interés principal es el flujo macroscópico, entendido como el flujo
global a través de una sección que contiene multitud de microcanales, y que puede ser considerada uniforme a efectos
macroscópicos, en contraste con la distribución aproximadamente parabólica que se tiene para un fluido newtoniano en
régimen laminar a través de un microcanal individual de sección circular. De este modo, en vez de trabajar con el medio
permeable real con toda su estructura microscópica, trabajaremos con un medio permeable ficticio de propiedades
continuas en cada punto que coinciden con el valor medio de las propiedades microscópicas del mismo en el entorno
de dicho punto. Esto es de nuevo la aproximación del continuo o hipótesis del continuo que se vio para el análisis del
movimiento de los fluidos, pero ahora extendida al medio permeable a través del cuan fluyen éstos.
Matemáticamente el flujo permeable es una situación análoga a la que se tiene en los procesos de conducción del
calor en un cuerpo entre puntos que se encuentran a distintas temperaturas, actuando en ese caso la temperatura como
magnitud potencial. De esta forma, se produce un flujo de calor a través del medio entre puntos que se encuentran
a distintas temperaturas, dependiendo este flujo de las características del medio. En el estudio del problema de las
filtraciones a través de un medio permeable, la magnitud potencial que regula el proceso será la energía por unidad de
peso del fluido, produciéndose entre puntos en los que el fluido tiene distinta energía por unidad de peso una filtración
que depende de las características del medio.

5.2.1. Velocidad de descarga y velocidad microscópica promedio.

La velocidad de descarga o velocidad de flujo se define como la cantidad de fluido que se filtra a través del medio
poroso por unidad de área y de tiempo. La velocidad de descarga no corresponde a la velocidad real del flujo del agua
subterránea ya que en realidad el agua en realidad circula por el medio permeable en una compleja trayectoria no
rectilínea, recorriendo grandes distancias, sino que tiene un significado únicamente macroscópico.
Obviamente, como el flujo puede tener lugar únicamente entre poros de terreno saturado interconectados entre sí,
la velocidad a través de cualquier sección puede únicamente considerarse en un sentido estadístico, lo que hacemos
a través del tratamiento del medio real como un medio continuo. Si denominamos xp a la relación efectiva entre el
área Ap correspondiente a los poros y el área total A de la sección, está claro que xp = Ap /A. El caudal filtrado es
entonces
Q = Sv = xp A¯m ,v (5–14)
donde xp vm es la velocidad de descarga y a vm se le denomina velocidad microscópica promedio† de filtración. Así,
¯ ¯
esta velocidad es la velocidad promedio a través de los canales individuales, mientras que la velocidad de descarga
† En algunos textos en castellano se utiliza la denominación velocidad de filtración para la velocidad de descarga. Como ésto no es generalizado

y en textos sajones es más frecuente lo contrario (denominar velocidad de filtración a la velocidad microscópica promedio), evitaremos en lo posible
el término velocidad de filtración.


es la velocidad promedio a través de todo el medio, incluyendo los granos de terreno. Realmente ni siquiera esta
velocidad es la velocidad real, sino sólo el promedio de las componentes de la velocidad microscópica en la dirección
del flujo. Dependiendo de la importancia de las componentes en las direcciones perpendiculares al flujo la velocidad
microscópica real puede ser aún más alta.
Indaguemos algo más en xp . Si consideramos una sección A del medio que se encuentra a una altura z y denomina-
mos Ap a la superficie de esa sección correspondiente a los poros, se tiene que a la altura z:

Promediando a un cilindro de altura h, el valor promedio de xp será

que, dividiendo por A y operando, da

xp =
¯
1
hA
xp =
¯

0
xp (z) =

h
1
h 0
Ap (z)

h

Axp (z)dz =
A
.

xp (z)dz

1
V 0
h
Ap (z)dz
(5–15)

(5–16)

(5–17)

donde V es el volumen total del cilindro, mientras que la integral es el volumen hueco total en el cilindro. Así, el valor
medio de xp es precisamente la porosidad n.
En general se trabaja con velocidades de descarga, y todo lo que sigue se referirá a velocidades de descarga. Sólo en
algunos casos muy concretos en los que se busca el comportamiento microscópico se trabaja con la velocidad real.

5.2.2. La ley de Darcy en una dimensión.

Como es bien sabido de la mecánica de fluidos libres, para flujo estacionario no viscoso e incompresible, la ecuación
de Bernoulli
p v2
¯
+ + z = Htot = cte (3–33)
ρg 2g
p v2
¯
establece que la suma de la alturas de presión , de velocidad o cinética y geométrica z es constante en cualquier
ρg 2g
punto a lo largo de una línea de flujo.† En el flujo a través de medios permeables como en el caso del permeámetro que
se muestra en la figura 5–8 es necesario tener en cuenta las pérdidas de energía producidas por la resistencia viscosa
en los poros. Por tanto es necesario considerar la ecuación de Bernoulli generalizada en la forma ya conocida
pA v2
¯ pB v2
¯
+ A + zA − ∆h = + B + zB (4–1)
ρg 2g ρg 2g
donde ∆h representa la pérdida total de carga (energía por unidad de peso del fluido) en la distancia ∆l. Se denomina
gradiente hidráulico promedio al cociente
¯ = − ∆h .
i (5–18)
∆l
Tomando distancias ∆l cada vez más pequeñas respecto del punto inicial, se tiene, en el límite infinitesimal, el
gradiente hidráulico en cada punto
∆h dh
i = − l´ ım =− (5–19)
∆l→0 ∆l dl
En la inmensa mayoría de los problemas de filtraciones, la altura de velocidad es tan pequeña que puede despreciarse.
Por ejemplo, una velocidad de filtración de 30 cm/seg, que es grande comparada con las velocidades de filtración
típicas, daría una altura de velocidad de 0,5 cm. Por tanto, la ecuación de Bernoulli para flujo en medios permeables
puede escribirse a efectos prácticos en la forma
pA pB
+ zA − ∆h = + zB (5–20)
ρg ρg
† Si el flujo es no rotacional, ésto es cierto en general entre dos puntos cualesquiera del tubo de flujo.


Agustin Martin DomingopA
ρg

HA
zA
A

∆l
∆h

B

Figura 5–8 Flujo a través de un permeámetro.

a través de las alturas piezométricas, que es lo que realmente podemos medir.
pB
ρg

zB
HB

Por tanto, en los problemas normales de permeación, es indiferente hablar de altura piezométrica o total y se deno-
mina a ésta potencial hidráulico Φ que tiene, en un punto genérico del flujo permeable, la forma
p
Φ = Htot = +z (5–21)
ρg

Debido a la compleja naturaleza de los terrenos, antes de 1856 el estudio del flujo permeable del agua en los mismos
desafiaba a todos los análisis. En ese año, Henry Darcy publicó una relación simple basada en sus experimentos en
“Les fontaines publiques de la ville de Dijon”[1]. Para ello utilizó el permeámetro que se muestra en la figura 5–9 y,
para distintos valores del caudal, analizó la pérdida de carga entre el principio y el final del permeámetro. De esta forma
observó que para caudales no demasiado grandes se satisface una relación de proporcionalidad entre la velocidad de
descarga v = Q/S y el gradiente hidráulico o pérdida de carga por unidad de longitud del medio permeable. Esta
ecuación,
Q ∆h dΦ
=v=k ⇒ v = ki = −k (5–22)
S ∆l dl
conocida como ley de Darcy, muestra una dependencia lineal entre el gradiente hidráulico y la velocidad de descarga v.
En Mecánica del Suelo, se denomina al coeficiente de proporcionalidad k coeficiente de permeabilidad o conductividad
hidráulica que como se ve claramente en la ecuación (5–22) tiene dimensiones de velocidad. Este coeficiente representa
la facilidad que ofrece el medio para dejarse atravesar por el fluido. Su valor depende del tamaño promedio de los poros
y está relacionado también con la distribución de tamaños y formas de las partículas y con la estructura del medio. Así,
éste puede cambiar en un factor de hasta 106 entre el de una grava y el de una arcilla. Incluso una pequeña proporción
de material fino en uno más grueso puede cambiar de forma importante su magnitud.
Aunque en la ecuación (5–22) la ley de Darcy se presenta en forma diferencial es necesario recalcar que en modo
alguno describe la situación microscópica en la red de poros del medio. De forma estricta, la ley de Darcy representa
el equivalente estadístico macroscópico de las ecuaciones del movimiento de Navier-Stokes para el flujo viscoso a
través de un medio permeable. Es precisamente esta equivalencia la que permite el tratamiento del flujo permeable en
el marco de una teoría de campos potenciales. A nuestros efectos, todas las contribuciones viscosas serán consideradas
a través de la ley de Darcy, y el flujo a todos los demás efectos será considerado como no viscoso o perfecto.


Figura 5–9 Permeámetro original utilizado por Darcy. Extraído de su informe Les fontaines publiques de la ville de Dijon, 1856[1].

Rango de validez de la ley de Darcy.

Como se ha visto en la sección 4.1.2 la observación visual de colorantes inyectados en los líquidos llevó a Reynolds
(1883) a concluir que el orden dentro del flujo depende de la velocidad del mismo. A velocidades pequeñas el flujo era
ordenado, en capas, es decir, laminar. Por el contrario, al aumentar la velocidad se produce una mezcla entre el agua y
el colorante, el flujo se hace irregular o turbulento.
Las distintas experiencias realizadas muestran de forma concluyente que para pequeñas velocidades (en el rango
de flujo laminar), la ley de Darcy (5–22) representa de forma adecuada el flujo a través de un medio permeable, con


una relación lineal entre el gradiente hidráulico y la velocidad del flujo. Según el flujo se va haciendo turbulento, el
comportamiento del gradiente hidráulico se aproxima mejor mediante una expresión de la forma

i=−

viscosidad dinámica y ν su viscosidad cinemática.

µ
dΦ
dx

=
= av + bv n

turbulento de agua a través de columnas de perdigones de tamaño uniforme, se obtiene n = 2 (Lindquist).

µ

vD
ν
ν

≤1
(5–23)

donde a y b son constantes y n varía entre 1 y 2, dependiendo de las características del flujo. Por ejemplo, para el flujo

Al igual que en el caso de los fluidos libres (véase la sección 4.5), el criterio para caracterizar el tipo de régimen,
laminar o turbulento, es utilizar el número de Reynolds (número adimensional que relaciona las fuerzas de inercia con
las viscosas), pero adaptado al flujo a través de medios permeables,

Re =
vDρ
=
vD
, (4–23)

donde ahora v es la velocidad de descarga, D el diámetro promedio de las partículas, ρ la densidad del fluido, µ su

El valor crítico del número de Reynolds en el cual el flujo permeable cambia de laminar a turbulento está entre 1 y
12. Sin embargo para nuestros efectos, será suficiente dar por válida la ley de Darcy cuando el número de Reynolds es
menor o igual que la unidad.
vDρ
(5–24)

Sustituyendo los valores de ρ y de µ para el agua en la ecuación anterior y asumiendo un valor de la velocidad de
descarga de 0,25 cm/seg, se obtiene D = 0,4 mm. Este es un valor típico del tamaño medio de partículas para la arena
gruesa. Por tanto parece bastante razonable que para medios permeables normales, en los cuales el tamaño medio de
grano es considerablemente menor, la validez de la ley de Darcy sea clara. De hecho, hay que subrayar que el carácter
laminar del flujo del agua a través de medios permeables en las condiciones dadas es uno de los pocos ejemplos válidos
de este tipo de flujo en situaciones reales de ingeniería hidráulica.

La conductividad hidráulica.

Hemos visto que para flujo laminar a través de un medio permeable se satisface la Ley de Darcy

dΦ
v = −k
dl
donde k es la conductividad hidráulica. La densidad y la viscosidad del fluido ejercen una influencia en la velocidad
con la que éste circula a través del medio permeable, por lo que es interesante aislar la parte de k que depende de
esas propiedades. Para ello, introducimos la permeabilidad física o permeabilidad intrínseca k0 (con unidades de
superficie), que es un parámetro que caracteriza las propiedades del medio y es independiente de las propiedades del
fluido. La relación entre la permeabilidad física y la conductividad hidráulica viene dada por
ρg g
k = k0 = k0 (5–25)
µ ν

Sustituyendo esta expresión en la ecuación de Darcy (5–22), se tiene

ρg dΦ g dΦ
v = −k0 = −k0 (5–26)
µ dl ν dl

que indica que la velocidad de descarga es inversamente proporcional a la viscosidad del fluido. Es interesante utilizar
esta última expresión cuando se trabaja con más de un fluido o cuando las propiedades de éste varían por ejemplo
a causa de cambios de temperatura. Aunque la permeabilidad física es idealmente la misma para distintos fluidos y
depende únicamente del medio permeable, en la práctica éste puede incluso cambiar de dimensiones con distintos
fluidos obteniéndose en ese caso distintos valores de la permeabilidad física


Tabla 5–4 Valores aproximados de la conductividad hidráulica y el tamaño medio de grano para distintos medios porosos.

Limo
Tipo de terreno
Gravas
Arena gruesa
Arena media
Arena fina
Arena limosa

Turba o limo fino
Arcilla meteorizada
Arcilla no meteorizada
Conductividad hidráulica
k (cm/seg)
≥ 1, 0
1,0 − 10−1
10−1 − 10−2
10−2 − 10−3
10−3 − 10−4
10−4 − 10−5
10−5 − 10−6
10−6 − 10−7
≤ 10−7
Tamaño medio
de grano (mm)
≥2
0,6 − 2
0,2 − 0,6
0,06 − 0,2
0,02 − 0,06
0,006 − 0,02
0,0002 − 0,006

En los problemas típicos de flujo a través de medios permeables, en los cuales estaremos interesados en el flujo de
un único fluido relativamente incompresible es más sencillo utilizar la ley de Darcy con k. Algunos valores típicos de
la conductividad hidráulica para diversos medios se dan en la tabla 5–4.
Puede relacionarse la permeabilidad física k0 con la porosidad n del material para flujo en régimen laminar como
el que se tiene en la región de validez de la ley de Darcy a partir de la fórmula de Poiseuille (4–21) y de la relación
(5–14). En la drástica aproximación de tuberías cilíndricas de radio R el correspondiente al tamaño promedio de los
poros, la velocidad de descarga queda en la forma

vm = n
¯
Q
A
=
nQ ρg ∆hpiezo
8 µ L
que al comparar con la ley de Darcy en función de la permeabilidad física (5–26) da un valor para la permeabilidad
física en la aproximación que utilizamos de
nR2
k0 = (5–27)
8

5.2.3. Potencial hidráulico y sus componentes. El campo hidráulico.

Como hemos dicho anteriormente, la resultante de las fuerzas que actúan en cada punto del flujo sobre el fluido
puede obtenerse como el gradiente de la energía potencial en dicho punto del medio poroso. Cuando se realiza un
trabajo sobre el fluido que circula por el medio poroso, esta energía potencial cambia en una cantidad igual al trabajo
realizado. Se denomina potencial hidráulico Φ a esta variación de la energía potencial del fluido respecto de la que
tenía en un estado de referencia, medida por unidad de peso del fluido, es decir, la energía total por unidad de peso del
fluido. Como hemos visto, hay distintas contribuciones a este potencial hidráulico, el potencial de presión o altura de
presión Φp = p/(ρg), el potencial gravitacional o altura geométrica Φz = z, y el potencial de velocidad o altura de
velocidad o cinética Φv = v 2 /(2g), que como hemos visto es despreciable en los problemas habituales de movimiento
de fluidos a través de un medio permeable.
Hay otra contribución al potencial hidráulico de la que aún no hemos hablado de forma explícita, y es el potencial
de matriz o de succión, Φm . Esta contribución está ligada a la presencia de fuerzas de ascensión (o descenso) capilar.
Éstas aparecen como consecuencia de la interacción del agua con el esqueleto del medio poroso o matriz de partículas
sólidas en que está embebida (véase la sección 5.1.3). Estas fuerzas capilares son en general atractivas para los suelos
y cerámicas porosas, pero pueden ser negativas para otros materiales porosos. Cuando las fuerzas son atractivas el
potencial de matriz tiene valores negativos.
Finalmente mencionar el potencial osmótico, que únicamente aparece cuando existen membranas semipermeables
en el medio poroso (por ejemplo, las raíces de las plantas que se encuentran en el suelo) y que no tendremos en cuenta.
Cuando entre dos puntos existe una diferencia de potencial hidráulico aparece una fuerza que intenta producir
movimiento del fluido hacia posiciones de menor energía potencial, y este movimiento se producirá únicamente entre


puntos con distintos potenciales hidráulicos. En general, el potencial hidráulico puede variar con el tiempo, y además
no está distribuido de forma uniforme (en ese caso no habría movimiento). Por tanto, lo que en realidad se tiene es una
distribución espacio-temporal del potencial hidráulico, el denominado campo hidráulico

que puede escribirse como
v = −k
dΦ
Φ = Φ(x, y, z, t)

en el mismo punto en el mismo instante, las superficies equipotenciales no pueden cruzarse.

5.2.4. La ley de Darcy en tres dimensiones.

vy = kiy = −k

vz = kiz = −k

ux +
dΦ
uy +
dΦ
dΦ
dy
dΦ
dz

uz = −k ∇Φ.
(5–28)

Se denomina superficie equipotencial hidráulica al lugar geométrico de los puntos del medio permeable en los cuales
el potencial hidráulico tiene el mismo valor. Obviamente, como el potencial hidráulico no puede tener valores distintos

Si en un punto genérico del flujo permeable a través de un medio isótropo observáramos la experiencia anterior en
tres direcciones mutuamente perpendiculares, x, y y z, se tendrían las expresiones para la ley de Darcy (5–22) en cada
una de estas direcciones como

vx = kix = −k
dΦ
dx

(5–29)
dx dy dz
para un medio permeable isótropo.
Si el medio es anisótropo la situación es más compleja, y en general la conductividad hidráulica viene representada
por un tensor con varias componentes distintas de cero, expresándose la ley de Darcy en un punto en la forma
dΦ
 
     dx 
vx kxx kxy kxz   dΦ 

v =  vy  = −  kyx kyy kyz    ˆ
 = −k ∇Φ (5–30)
 dy 

vz kzx kzy kzz 
dΦ 
dz
ˆ
El tensor k es un tensor simétrico de segundo orden, y por lo tanto existe un sistema de referencia (sistema de ejes
principales) en el cual todos los elementos de fuera de la diagonal son nulos, denominándose a los elementos de
la diagonal en ese sistema elementos principales del tensor o, en este caso concreto, conductividades hidráulicas
principalesconductividad hidráulica!principal. En este sistema de ejes principales, r1 , r2 y r3 , la expresión anterior
toma la forma
dΦ
 
     dr1 
v1 k1 0 0  
 v2  = −  0 k2 0   dΦ  .
 
v=  dr2  (5–31)
v3 0 0 k3   dΦ 


dr3
El vector gradiente es un vector que es siempre perpendicular a las superficies equipotenciales del campo escalar en
cada punto. Por tanto al ser, para un material isótropo v = −k ∇Φ la velocidad de descarga es siempre perpendi-
cular a las superficies de potencial hidráulico constante para un material isótropo. Por el contrario, para un material
anisótropo, aunque el gradiente sigue siendo perpendicular a las equipotenciales, la velocidad no lo será al aparecer


la conductividad hidráulica en forma de tensor en la ecuación (5–31) con distintos valores de k. Sin embargo, sí se
cumplirá para los materiales isótropos ficticios equivalentes a uno anisótropo que se obtienen redimensionando los
ejes y se tratan en la sección 5.2.7.

5.2.5. Flujo a través de terrenos saturados homogéneos e isótropos.

La conductividad hidráulica de un suelo saturado es constante si su estructura es estable, porque todos sus poros
y capilares están ocupados por el agua, igual que cualquier otro material poroso en idénticas condiciones. Por el
contrario, la conductividad hidráulica de un suelo no saturado puede cambiar continuamente en función de los cambios
que se produzcan en el potencial de matriz. Asimismo, en algunas situaciones como durante las lluvias se tienen
cambios temporales del contenido de agua en el suelo. Lo mismo ocurre si se está produciendo una degradación del
medio o del fluido durante el proceso de flujo permeable.
En el caso general en el que puede haber variación temporal del flujo, es necesario tener en cuenta no sólo la ley de
Darcy, sino también la ecuación de continuidad
∂ρ
∂t
= −∇(ρv).

Sin embargo cuando se tiene en cuenta que la compresibilidad de los líquidos es muy pequeña, la densidad tiene una
variación despreciable con el tiempo y la posición, y la expresión anterior se puede escribir como

∇v = 0
(3–9)

(3–11)

Estudiaremos ahora como es la ecuación diferencial que describe el flujo. Para ello nos limitaremos al caso en que el
medio permeable es estable, a régimen estacionario y a un fluido incompresible como el agua.
Si el medio es además isótropo, el tensor conductividad hidráulica es un tensor diagonal con todas las componentes
de la diagonal iguales a un único valor de la conductividad hidráulica, y la ley de Darcy puede escribirse como

v = −k ∇Φ (5–29)

Sustituyendo esta expresión en la ecuación anterior se tiene que la ecuación diferencial que hay que resolver para
obtener la distribución de potencial hidráulico es

∂2Φ ∂2Φ ∂2Φ
∇2 Φ = + + = 0. (5–32)
∂x2 ∂y 2 ∂z 2
Esta ecuación diferencial es la conocida ecuación de Laplace, ampliamente estudiada en el electromagnetismo y en
el tratamiento de algunos problemas de conducción del calor. Como se observa en esta ecuación, la forma del campo
hidráulico no depende de la conductividad hidráulica, sino que depende únicamente de las condiciones de contorno en
los límites, es decir de la forma del recinto y de las características de estos límites. Obviamente, el flujo sí depende de
la conductividad hidráulica a través de la ley de Darcy.

5.2.6. La ecuación de Laplace y su resolución.

Como se decía anteriormente, la ecuación de Laplace ha sido ampliamente estudiada en Física Matemática por
su interés en diversos campos que van desde electromagnetismo a la conducción del calor. Su resolución depende
principalmente de la geometría del recinto espacial en el que se verifica ésta (como por ejemplo un suelo, un material
poroso como la zapata de cimentación de un edificio, un muro de carga o una presa de gravedad) y de las condiciones
de contorno que se cumplen en los límites del recinto. Es decir, la resolución de la ecuación de Laplace depende de
las características o los valores que ha de tener el potencial hidráulico (o su gradiente, directamente relacionado con la
velocidad de descarga) en la superficie límite del recinto en estudio. Estas condiciones de contorno vienen, en general,
impuestas por la naturaleza física de estos límites, como por ejemplo, la presencia de una capa impermeable.
Se dice que estas condiciones de contorno son condiciones de Dirichlet o de primera especie cuando lo que se
especifica son los valores del potencial hidráulico en la superficie límite, mientras que se dice que estas condiciones


son condiciones de von Neumann o de segunda especie si lo que se especifica son condiciones en el gradiente del
potencial hidráulico (en la velocidad). Se dice que se tienen unas condiciones de contorno mixtas cuando sobre una
parte de la superficie límite se conocen los valores del potencial hidráulico y sobre otra parte se conoce información
sobre la velocidad.
Las técnicas de resolución de la ecuación de Laplace son muy conocidas en el ámbito de la Física Matemática y
utilizan fundamentalmente el método de separación de variables en coordenadas cartesianas, esféricas o cilíndricas,
según la simetría del recinto permeable. Estos son los casos más favorables y al mismo tiempo, los menos frecuen-
tes. Si el recinto no tiene una simetría apropiada se utiliza, como técnica alternativa, la transformación conforme en
variable compleja o el método de imágenes. Para la mayor parte de los problemas que son intratables analíticamente,
es necesario recurrir a métodos numéricos de resolución de ecuaciones diferenciales, como los de diferencias fini-
tas o elementos finitos. A menudo se utilizan estos métodos también para abordar problemas que, aunque resolubles
analíticamente, requieren un excesivo esfuerzo.

Condiciones de contorno para la resolución de la ecuación de Laplace.

En un medio permeable saturado aparecen distintos tipos de condiciones de contorno, todas presentes en el ejemplo
de la presa de tierra que se muestra en la figura. Consideremos en detalle las características de cada una de ellas.

h
Superficie filtrante
B
Superficie libre

C
Superficie de escorrentía

A
D

Superficie impermeable

Figura 5–10 Presa de tierra en la que aparecen los cuatro tipos de condiciones de contorno que estudiaremos.

Superficies impermeables.
Una superficie impermeable separa el medio permeable (saturado) de un medio impermeable con el que está en
contacto. El fluido no puede penetrar a través de una superficie impermeable, ya que las características del material
impermeable lo impiden. Por tanto, la componente de la velocidad perpendicular a la superficie debe ser cero en cada
punto. Sin embargo no se impone restricción alguna a sus componentes tangenciales. Si denominamos n y t a las
direcciones normal y tangencial, respectivamente, en un punto de la superficie, se tiene

∂Φ
= 0. (5–33)
∂n
es decir, una superficie impermeable define el lugar geométrico de una serie de líneas de flujo. Al mismo tiempo una
línea de flujo satisface la condición de una superficie impermeable y, de hecho, las paredes de un tubo de flujo, se
consideran a todos los efectos como superficies impermeables.
Hay dos tipos de superficies impermeables típicas. Una es la superficie impermeable que tiene el fluido por encima,
como ocurre cuando se considera el flujo a través de un medio permeable que se encuentra sobre una base de otro
material de conductividad hidráulica mucho menor que la del medio permeable. En este caso la superficie impermeable
está definida por las líneas de corriente más bajas. También puede tenerse una estructura impermeable sobre el terreno,
como una presa de gravedad, o la pared vertical de un muro.


En un material isótropo, al ser las superficies equipotenciales del campo hidráulico perpendiculares a la velocidad
de descarga, éstas deben ser perpendiculares asimismo a la superficie impermeable. Ésto no es cierto para un mate-
rial anisótropo en el que formarán ángulos distintos de 90◦ . Sin embargo sí se cumple para los materiales isótropos
equivalentes a un anisótropo que se tratarán en la sección 5.2.7.

Superficies filtrantes.
Las superficies filtrantes corresponden a los límites del medio permeable en contacto con los límites del fluido libre.
En estas superficies filtrantes el valor del potencial hidráulico será el del fluido libre en contacto con la misma, es
decir,
Φ=
p
ρg
+ z = cte. (5–34)

Las superficies filtrantes en contacto con una masa de fluido libre son por tanto superficies equipotenciales, cuyo
potencial hidráulico es igual a la altura de carga total del fluido en contacto con el medio permeable. Obviamente,
distintas superficies filtrantes del medio permeable en contacto con distintas masas de fluido libre con distintas alturas
de carga totales tendrán distintos potenciales hidráulicos.
Para materiales isótropos la velocidad de descarga será perpendicular a las superficies equipotenciales del campo
hidráulico y por tanto a las superficies filtrantes. Para materiales anisótropos esto no será cierto, a no ser que se trate
de los materiales isótropos equivalentes a uno anisótropo que se tratan en la sección 5.2.7.

Superficie de escorrentía.
La superficie de escorrentía representa un límite en el cual el flujo emergente del medio permeable sale a una zona
libre de líquido y de terreno, es decir, sale a la atmósfera. Como la presión a la salida es la atmosférica, se cumple en
esta superficie que el potencial hidráulico es
Φ=z (5–35)
En principio no hay restricciones en v, salvo que obviamente debe tener una componente perpendicular a la superficie.

Superficie libre o línea de filtración.
La superficie libre está compuesta por las líneas de flujo superiores en el terreno permeable y coincide con el nivel
freático. Esta superficie se encuentra dentro del medio permeable y separa la región de flujo saturado de la región del
medio permeable a través de la cual no hay flujo. Se diferencia de la superficie impermeable en que en este caso son
las características del flujo las que hacen que no se atraviese esta superficie, en vez de las características del material
de la superficie, que para la superficie libre es el mismo que en la zona en la que sí hay flujo. Su determinación es uno
de los principales problemas que surgen a la hora de estudiar de forma analítica el problema del flujo permeable. Por
un lado esta superficie debe comportarse como una superficie impermeable, al estar compuesta por líneas de flujo, y
por otro lado es evidente que la presión en cada punto de esa superficie es constante e igual a la presión atmosférica,
por lo que el potencial hidráulico será
Φ=z (5–36)
lo que muestra que el potencial hidráulico a lo largo de una superficie libre varía linealmente con la altura.

5.2.7. Flujo permeable a través de medios anisótropos

Hemos considerado el caso de un flujo isótropo en el que el coeficiente de permeabilidad del medio en un punto es
independiente de la dirección del flujo en ese punto. Si además el medio tiene el mismo coeficiente de permeabilidad
en todos los puntos del recinto se dice que es homogéneo e isótropo. Cuando el coeficiente de permeabilidad depende
de la dirección y la dependencia es la misma para todos los puntos se dice que el medio es homogéneo y anisótropo.
La mayor parte de los terrenos son anisótropos en alguna medida. Por ejemplo, los terrenos sedimentarios están
habitualmente formados por finas capas alternas y además, la tensión de compresión en la dirección horizontal es
distinta de la que se tiene en la dirección vertical. Generalmente en terrenos sedimentarios homogéneos naturales, el
coeficiente de permeabilidad es mayor en la dirección horizontal que en la vertical. En esta sección se abordará lo

PSfrag replacements

Perpendicular

siguiente:
Paralelo
k1

k2

kn−2
kn−1
kn

y
b

...

L
d1

d2

dn−2
dn−1
dn

Figura 5–11 Flujo paralelo y perpendicular a través de un material permeable formado por capas de distintos materiales.

Flujo paralelo y perpendicular a través de un medio anisótropo compuesto de capas isótropas.

Consideremos una sección vertical de un terreno estratificado formado por n capas isótropas, cada una de ellas de
espesor d1 , d2 , . . . , dn , con unas conductividades hidráulicas k1 , k2 , . . . , kn , respectivamente y analicemos dos casos
particulares interesantes.
x

Flujo paralelo: Consideremos en primer lugar el flujo paralelo a las capas. El caudal filtrado será la suma de los
caudales filtrados a través de cada una de las capas,
n n
hi − hf hi − hf
Qi = di bki ⇒ Q= Qi = b ki di
L i=1
L i=1

donde h1 − h2 es la pérdida de carga en las capas en una distancia L. La velocidad promedio en la dirección x será
n
Q ki di hi − hf
vx = =
S i=1
d L

por lo que se puede escribir una conductividad hidráulica efectiva cuya componente x es
n
ki di
kx = (5–37)
i=1
d

Flujo perpendicular: Consideremos ahora el flujo perpendicular a las capas. Para que se satisfaga la ecuación de
continuidad la componente y de la velocidad debe ser la misma en todas las capas. Como la pérdida total de carga
debe ser igual a la suma de las pérdidas de carga en cada capa, se tiene
n n
di di
∆hi = v ⇒ ∆h = ∆hi = v
ki i=1 i=1
ki
y
hi − hf d hi − hf
v= n = n
di di d
k
i=1 i i=1
ki


que corresponde a una conductividad hidráulica efectiva de componente y

Flujo a través de un medio anisótropo general.

∇v =
ky = n

i=1

v3 = −k3

∂v1
∂r1
+
∂v2
∂r2
∂r2
∂Φ
∂r3
d

di /ki

En el sistema de ejes principales la ecuación de continuidad se escribe como

+
∂v3
∂r3
=0

En un medio homogéneo como el que estamos considerando, [kij ] no depende del punto, y por tanto es constante
frente a las derivadas espaciales, quedando, al sustituir la ley de Darcy en la ecuación de continuidad
(5–38)

Para estudiar lo que ocurre en un medio anisótropo consideraremos el problema descrito en los ejes principales del
tensor conductividad hidráulica. En esos ejes, en función de las conductividades hidráulicas principales k1 , k2 y k3 ,
las componentes de la velocidad de descarga que se obtienen de la Ley de Darcy quedan:

v1 = −k1

v2 = −k2
∂Φ
∂r1
∂Φ
(5–39a)

(5–39b)

(5–39c)

(5–40)

∂2Φ ∂2Φ ∂2Φ
k1 2 + k2 ∂r2 + k3 ∂r2 = 0. (5–41)
∂r1 2 3

Si consideramos el siguiente cambio de variable, en función de una constante arbitraria k con dimensiones de conduc-
tividad hidráulica,
k k1
q1 = r1 r1 = q1
k1 k
k k2
q2 = r2 ⇐⇒ r2 = q2 (5–42)
k2 k
k k3
q3 = r3 r3 = q3
k3 k
la ecuación (5–41) se puede reescribir como

∂2Φ ∂ 2Φ ∂ 2Φ
k 2 + k ∂q 2 + k ∂q 2 = 0. (5–43)
∂q1 2 3

que, eliminando k se puede escribir como

∂2Φ ∂2Φ ∂2Φ 2
2 + ∂q 2 + ∂q 2 = ∇q Φ = 0.
∂q1 2 3

que es de nuevo la ecuación de Laplace, pero ahora en el sistema de coordenadas reescalado. Así, mediante una simple
expansión o contracción de las coordenadas espaciales podemos transformar una región del flujo anisotrópa en una
región isótropa ficticia en la que la ecuación de Laplace es válida y en la que en principio podríamos obtener una
solución para ésta. Una vez se haya resuelto el problema para esta región ficticia, para obtener la solución para el flujo
real que se tiene basta con invertir la transformación.
El hecho de que se cumpla la ley de Darcy (5–31) no puede verse afectado por la elección del sistema de ejes y
en concreto por el redimensionado de los mismos que acabamos de tratar. No solo eso, si en la ley de Darcy (5–31)


realizamos el cambio de variables (5–42) tanto en las componentes dri /dt de la velocidad como en las ∂Φ/∂ri del
gradiente la ley de Darcy queda, para el material anisótropo en el sistema de coordenadas reescalado en la forma
vq = −k ∇q Φ similar a la ley de Darcy para un material isótropo (5–29) en función del parámetro k. Esto implica
que en el material anisótropo analizado en estos ejes como un material isótropo ficticio la velocidad de descarga es
también perpendicular a las superficies equipotenciales del campo hidráulico expresado también en estos ejes.

5.2.8. La red de flujo.

Se denomina red de flujo a una representación gráfica de las líneas de corriente y de sus correspondientes superfi-
cies equipotenciales en una región del flujo. La construcción de la red de flujo es uno de los métodos más valiosos
empleados en el análisis de los problemas de flujo en dos dimensiones, aunque a menudo requiere tediosos ajustes de
prueba y error. Al construir la red de flujo se llena la zona en permeación con una serie de líneas de flujo y otra serie
de superficies equipotenciales que forman un ángulo recto entre unas y otras en los puntos en que se cruzan.
Nótese que la construcción de la red de flujo sólo proporciona una solución aproximada al problema, aunque puede
ser buena incluso en el caso de geometrías complejas a partir de una serie de reglas sencillas (inicialmente debidas
a Phillip Forchheimer en torno a 1900 y después mejoradas por Arthur Casagrande en 1937). El método habitual de
construcción de la red de flujo es el siguiente:

1. Dibujar a escala los límites del recinto de forma que puedan hacerse terminar en dichos límites todas las líneas de
corriente y equipotenciales. En caso de tratarse de un material anisótropo homogéneo debe primero reescalarse
el problema en la forma descrita en la sección anterior.
2. Trazar de forma aproximada unas pocas líneas de flujo, teniendo en cuenta las condiciones en los límites del
recinto. Es decir, entre otras cosas, que las líneas de corriente han de ser tangentes a las superficies impermeables
y perpendiculares a las superficies filtrantes. No debe perderse la noción de que éstas son sólo unas pocas de
entre el continuo de líneas de corriente. Como una ayuda a la hora de espaciar las líneas de corriente, obsérvese
que la distancia entre líneas de corriente consecutivas aumenta en la dirección de mayores radios de curvatura.
Es conveniente empezar en las zonas más sencillas, por ejemplo en aquéllas en las que las líneas de flujo
son paralelas entre sí, lo mismo que las equipotenciales, de forma que los “cuadrados” de la red de flujo son
realmente cuadrados, y después seguir hacia las zonas de mayor complejidad.
3. Trazar las equipotenciales correspondientes, considerando que deben intersectar a ángulos rectos a todas las
líneas de corriente para materiales isótropos o anisótropos reescalados. Al inscribir un círculo en uno de estos
“cuadrados curvilíneos” éste debe tocar cada lado sólo una vez. Estos “cuadrados” deben cambiar de tamaño
gradualmente.
4. Afinar las líneas de corriente y equipotenciales trazadas para que se satisfagan bien las condiciones anteriores.
Asimismo debe tenerse en cuenta la simetría del problema por si es posible simplificar la representación.
5. Se recomienda como una prueba adicional comprobar que las diagonales de los “cuadrados” formados por líneas
de corriente y equipotenciales forman curvas sin variaciones excesivamente bruscas que se cruzan formando
ángulos rectos.

5.2.9. Ejemplos de redes de flujo en distintos sistemas

Presas de tierra

Cuando las presas de tierra están en las condiciones que se muestran en la figura 5–10, con una superficie de
escorrentía de esa forma, puede producirse perforación de la presa como consecuencia de una erosión interna que
poco a poco va arrastrando la tierra y produciendo un canal por el interior de la presa. Esta erosión, denominada
tubificación, que comienza en la parte aguas abajo y va propagándose por el interior de la presa a través de líneas de
conductividad hidráulica máxima puede evidenciarse muchos años después de que ésta se ha construido.
Para minimizar el riesgo de perforación se utilizan distintos sistemas, como paredes impermeables enterradas o la
colocación de un núcleo impermeable en la presa. Uno de los sistemas consiste en controlar el efecto de la filtración
mediante un dren interno como se muestra en la figura 5–12 diseñado de forma que la superficie libre está siempre por


h

base impermeable

Φ10 = h
superficie filtrante

superficie libre
superficie de escorrentía

Figura 5–12 Red de flujo en una presa de tierra típica con un dren interno.

h
dren

Φ0 ≃ 0 m

Φ = 0,9 h Φ = 0,1 h

Φ = 0,8 h Φ = 0,2 h

Φ = 0,7 h Φ = 0,3 h
Φ = 0,6 h Φ = 0,4 h
Φ = 0,5 h

Figura 5–13 Red de flujo en un presa de gravedad con su base al nivel del suelo.Las superficies equipotenciales se muestran con líneas de trazos y
las líneas de flujo con líneas continuas.

h

Φ10 = h Φ0 ≃ 0 m

Φ9
Φ1

Figura 5–14 Red de flujo en una presa de gravedad ligeramente enterrada y con una pantalla impermeable y una base impermeable a una cierta
profundidad.


Agustin Martin Domingo h0
p
ρg

Figura 5–15 Distribución de presiones bajo una presa
h1

debajo de la superficie externa de la presa y por tanto no haya superficie de escorrentía más que en el dren, donde se
controla que la tierra no pase al mismo mediante rejillas y filtros de tamaño de poro adecuado.

Presas de gravedad

En la figura 5–13 se muestra un ejemplo de red de flujo bajo una presa de gravedad sencilla cuya base está al
mismo nivel del terreno, que se supone homogéneo e isótropo. Como se ve en la figura, en estas condiciones el
flujo es simétrico y los potenciales hidráulicos van desde un valor máximo igual a h en la superficie filtrante aguas
arriba a un valor mínimo aproximadamente nulo (suponiendo un caudal muy pequeño) aguas abajo. Las superficies
equipotenciales se muestran con líneas de trazos y las líneas de flujo con líneas continuas.
Puede mejorarse notablemente el comportamiento de esa presa frente a las filtraciones mediante dos sencillas modi-
caciones en el diseño. Por una parte se entierra parcialmente la presa, se forma que su parte inferior queda por debajo
del nivel del terreno y por otra se coloca una pantalla impermeable en su parte inferior como se muestra en la figura
5–12. Estas medidas hacen que el peligro de perforación (tubificación) se reduzca notablemente así como que las
velocidades de salida aguas abajo se reduzcan como consecuencia del alargamiento de los recorridos de las líneas de
flujo, dando como resultado una disminución de las filtraciones bajo la presa.

5.2.10. Subpresiones

La presencia de filtraciones bajo una estructura como la de una presa introduce un elemento adicional a tener en
cuenta en cuanto a su estabilidad, debido a la aparición de una serie de presiones en la base del elemento, denominadas
subpresiones. Estas subpresiones producen por una parte un empuje y por otra parte un momento de vuelco a causa de
la no homogeneidad de su distribución, debiendo ambos efectos ser tenidos en cuenta.
Por ejemplo, en una presa en la que no se han tomado medidas para aminorar el efecto de estas subpresiones la
distribución de presiones en su base puede ser, en la práctica, de la forma que se muestra en la figura 5–15, apro-
ximadamente lineal. Esto hace que, a la hora de estimar el momento de vuelco sea necesario considerar, además de


los momentos debidos a las fuerzas actuando sobre las paredes verticales el efecto de las presiones sobre la base.
Asimismo, si el suelo está por debajo del nivel freático puede ser necesario tomar medidas para evitar el efecto del
empuje sobre éste.
Para amortiguar el efecto de las subpresiones se instalan, cerca de la zona aguas arriba, drenajes verticales dentro de
la presa desde la base de la misma a una galería de desagüe a un nivel superior al nivel aguas abajo. Este mismo sistema
se utiliza en losas grandes, donde es preferible taladrar verticalmente el piso para aliviar el efecto de las subpresiones
y evacuar el agua ﬁltrada. En otros casos se utilizan tirantes anclados bajo el piso para contrarrestar el efecto de las
subpresiones.

B
Indice alfabético

A
aceleración

altura
convectiva, 36
de transporte, 36
local, 36

cinética, 38, 83
de presión, 38, 83
de velocidad, 38, 83
geométrica, 38, 83
metacéntrica, 22
piezométrica, 38
total, 38
ángulo de contacto, 73
aproximación del continuo, 78
Arquímedes, principio de, 21
atmósfera
adiabática, 17
de densidad homogénea, 17
isoterma, 17
tipo, 17

Bernoulli, ecuación de, 37, 38, 40
condición

D
de contorno, 85
de Dirichlet, 85
de primera especie, 85
de segunda especie, 86
de von Neumann, 86
mixta, 86
de no deslizamiento, 49
conductividad hidráulica, 80, 82
contacto
ángulo de, 73
línea de, 73
continuo, 6
hipótesis del, 7, 78
curvas
de corriente, 50
reológicas, 50

Darcy, Henry, 80
Darcy-Weiesbach, ley de, 57
densidad, 7
diagrama de Moody, 59
diámetro hidráulico, 58
Bernoulli, ecuación generalizada de, 63 E
Bingham, material lineal de, 51 ecuación
Blausius, fórmula de, 59 de Bernoulli, 40
generalizada, 63
C de Bernoulli, 37, 38, 79
campo de Bernoulli generalizada, 79
de densidades, 25 de continuidad, 29
de presiones, 25 de continuidad,
de velocidades, 25 forma diferencial, 29, 31
hidráulico, 84 forma integral, 28
capa límite, 56 de Euler, 32
capilaridad, 74 de Hagen-Poiseiuille, 58
carga, pérdida de, 46 de Laplace, 73, 85
caudal, 29 de Laplace-Young, 73
másico, 28 general de la hidrostática, 12
cavitación, 65 ecuaciones de Navier-Stokes, 54
centro efecto Venturi, 42
de gravedad, 19 eje
de presiones, 19 de empuje, 22
choque hidráulico, 65 de flotación, 22
coeficiente equilibrio
de compresibilidad, 8, 67 estable, 22
de dilatación de volumen, 9 indiferente, 23
de fricción, 57 inestable, 23
de permeabilidad, 80 Euler
de viscosidad dinámica, 47 ecuación de, 32
Colebrook, fórmula de, 59 método de, 25
compresibilidad, 8
coeficiente de, 8 F
módulo de, 8 fluidez, 5

95

96 Indice alfabético

fluido

flujo

G
anómalo, 50
dilatante, 50
newtoniano, 47
no newtoniano, 50
pseudoplástico, 50
viscoelástico, 51
viscoplástico, 51

compresible, 26
estacionario, 26
incompresible, 26, 29
irrotacional, 26, 36
no estacionario, 26
no viscoso, 26
rotacional, 26
viscoso, 26
fórmula
de Blausius, 59
de Colebrook, 59
de Newton, 47
de Poiseuille, 54
de von Karman, 59
fuerza hidrostática
sobre superficies curvas, 21
sobre superficies planas, 18

golpe de ariete, 65
M
de corriente, 27
de filtración, 87
de traza, 26

manómetro, 15
diferencial, 15
material lineal de Bingham, 51
medio
anisótropo, 87
continuo, 6
homogéneo, 87
isótropo, 87
permeable, 76
metacentro, 22
método
de Euler, 25
de Lagrange, 25
módulo
de compresibilidad, 8, 67
de Young, 67
Moody, diagrama de, 59

N
Navier-Stokes, ecuaciones de, 54, 80
Newton, ley de, 47, 55
nivel
capilar, 77
freático, 76, 87
gradiente número
hidráulico, 79 de cavitación, 65
hidráulico promedio, 79 de Reynolds, 56, 82

H P
hipótesis del continuo, 7, 78 paradoja de Pascal, 13
Pascal,
I paradoja de, 13
índice de poros, 77 principio de, 13
pérdida de carga, 46
L lineal, 57
Lagrange, método de, 25 local, 60
laminar representación gráfica, 63
flujo, 81 singular, 60
régimen, 46 permeabilidad
Laplace, ecuación de, 73, 85 física, 82
Laplace-Young, ecuación de, 73 intrínseca, 82
ley permeámetro, 79
de Darcy-Weiesbach, 57 peso específico, 8
de conservación de la masa, 29 piezómetro, 38
de Darcy, 80 Poiseuille
de la potencia, 49 fórmula de, 54
de Newton, 47, 55 régimen de, 46
de Sutherland, 49 porosidad, 77
límite de fluidez, 51 potencial
línea de matriz, 77, 83
de flujo, 27 de presión, 83
de alturas de succión, 83
piezométricas, 63 de velocidad, 83
totales, 63 gravitacional, 83
de contacto, 73 hidráulico, 80, 83

Indice alfabético 97

R

red
presa
osmótico, 83

de gravedad, 92
de tierra, 86, 90
presión, 8, 11
dinámica, 41
estática, 41
principio
de Pascal, 13
de Arquímedes, 21
producto de inercia, 20

radio hidráulico, 58

de flujo, 90
de tuberías, 64
régimen
de Poiseuille, 46
de Venturi, 46
laminar, 46
turbulento, 46
reogramas, 50
Reología, 50
Reynolds, número de, 56, 82
rugosidad, 58

S
relativa, 58
V
flujo, 81
régimen, 46

vasos comunicantes, 14
velocidad
característica, 56
de descarga, 78
de flujo, 78
microscópica promedio, 78
Venturi
efecto, 42
régimen de, 46
venturímetro, 42
viscosidad, 5, 8, 45
cinemática, 48
dinámica, 47
efectiva, 51
en gases, 48
en líquidos, 48
volumen específico, 8
von Karman, fórmula de, 59
vorticidad, 36

Z
zona de fluido libre, 57

senda, 26
Stokes, teorema de, 37
subpresiones, 92
superficie
de escorrentía, 87
filtrante, 87
impermeable, 86
libre, 87
superficie libre, 34

T
tensión
cortante en la pared de una tubería, 54, 58
de desplazamiento, 51
inicial de desplazamiento, 51
superficial, 69
teorema
de Stokes, 37
de Torricelli, 40
Torricelli, teorema de, 40
tuberías
en paralelo, 64
en serie, 64
redes complejas, 64
tubificación, 90
tubo
capilar, 74
de corriente, 28
de flujo, 28
en U, 14
turbulento

Indice alfabético
98


Referencias

[1] Darcy, H., Les Fontaines Publiques de la Ville de Dijon. Dalmont, París, 1856.
[2] White, F. M., Mecánica de fluidos. McGraw-Hill, México, 1983.
[3] Shames, I. H., Mecánica de los fluidos. Ed. del Castillo, S.A., 1970.

[4] Streeter, V. L. y E. B. Wylie, Mecánica de los fluidos. McGraw-Hil / Interamericana de México, S.A., México,
1988.
[5] Giles, R. V., Mecánica de los fluidos e hidráulica. McGraw-Hill, México, 1988.
[6] Rabinóvich, Y. Z., Hidráulica. Ed. Rusia-URSS, Moscú, 1994.
[7] Feynmann, R. P., R. B. Leighton y M. Sands, Electromagnetismo y materia, vol. 2 de Física. Addison-Wesley
Iberoamericana, Wilmington, 1987.
[8] Iribarne, J. V. y W. L. Godson, Termodinámica de la atmósfera. Monografías, Ministerio de Medio Ambiente.
Dirección General del Instituto Nacional de Meteorología, Madrid, 1996.
[9] Tipler, P. A., Física. 2 vols, Reverté, Barcelona, 1992.
[10] Isenberg, C., The science of soap films and bubles. Tieto Ltd, Clevedon, 1978.
[11] Prandtl, L. y O. G. Tietjens, Fundamentals of hydro and aeromechanics. Dover, N. Y., 1957.
[12] Landau, L., A. Ajiezer y E. Lifshitz, Curso de Física general: mecánica y Física molecular. Ed. Mir, Moscú,
1979.
[13] Harr, M. E., Groundwater and seepage. Dover, New York, 1991.
[14] Jiménez Salas, J. A. y J. L. de Justo Alpañes, Propiedades de los suelos y de las rocas, vol. I de Geotecnia y
cimientos. Rueda, Madrid, 2a edición, 1975.
[15] Marshall, T. J., Soil Physics. Cambridge University Press, 1988.
[16] González de Posada, F., F. A. González Redondo, M. González Redondo y M. D. Redondo Alvarado, Hidráulica
del medio permeable: teoría físico-matemática de Darcy. Ed. por el autor, ETSAM, Madrid, 1992.

99

Amd apuntes-fluidos

Más contenido relacionado

Similar a Amd apuntes-fluidos (20)

Más de GROVER MORENO (13)

Último (20)

Amd apuntes-fluidos