UNIVERSIDAD NACIONAL DE CAJAMARCA

FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

INDICE

1. Introducción………………………………………………………… Pág. 2

2. Objetivos ……………………………………………………….…… Pág. 2

3. Justificación…………………………………………………………. Pág. 3

4. Métodos Numéricos Aplicados A La Mecánica De Fluidos…........ Pág. 3

5. Marco Teórico ……………………………………………………… Pág. 4

5.1. Ley de Darcy …………………………………………………… Pág. 4

5.2. Ley de Laplace …………………………………………………. Pág. 7

5.2.1 Solución de la Ecuación de Laplace.……………………... Pág. 10

5.2.2 Ecuación de Laplace en Mecánica de fluidos…………… Pág. 12

Ejemplo de Aplicación …………………..……………… Pág. 16

6. Redes de flujo…………………..……………………………………. Pág. 17

6.1. Propiedades de las líneas de corriente y equipotenciales…….. Pág. 18

6.2. Trazado de la red de flujo……………………………………… Pág. 19

6.3. Condiciones hidráulicas de borde………………………………Pág. 20

6.4. Recomendaciones para el trazado de la red de flujo…………..Pág. 21

6.5. Cálculo del gasto o caudal en suelos isótropos…………………Pág. 22

Ejemplo de Aplicación……………………………………Pág. 24

7. Conclusiones…………………………………………………………..Pág. 27

8. Referencias Bibliográficas……………………………………………Pág. 27

MECANICA DE FLUIDOS I 1
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

Métodos Numéricos Aplicados A La Mecánica De Fluidos – Redes de flujo

I. Introducción
La presión intersticial que existe en un suelo con frecuencia no es la que corresponde a
las condiciones hidrostáticas, sino aquella creada por el flujo de agua a través de los poros
de suelo, como ocurre en las presas de hormigón sobre un estrato de suelo u otros
escenarios de obras civiles.
El agua se almacena a una altura determinada aguas arriba, la diferencia de nivel del agua
entre los lados de la presa creará una filtración a través de la cimentación desde el lado de
aguas arriba hasta el lado de aguas abajo. Cuando el flujo comienza, la presión intersticial
en el suelo pasa de los valores iniciales a unos valores finales que deberán ser compatibles
con las nuevas condiciones de frontera hidráulicas, y además con los cambios de volumen
que se producen en la masa de suelo. Durante este periodo el flujo varía en función del
tiempo y se denomina flujo transitorio. Cuando la presión intersticial en toda la masa de
suelo se equilibra con las nuevas condiciones de frontera, el flujo se vuelve independiente
del tiempo, y en ese caso se denomina flujo estacionario.
Las presas de hormigón suelen ser muy extensas en una dirección como para poder asumir
un comportamiento 2-D, que permite definir todas las variables del problema por unidad
de longitud (o espesor) en la dirección seleccionada, sin perjuicio de la desviación de las
soluciones en los extremos de la estructura.

En el presente informe se estudiara la importancia de los métodos numéricos para la

solución de problemas de mecánica de fluidos especialmente sobre redes de flujo que son
necesarios de ser estudiados debido a la importancia e influencia que tienen en presas
hidráulicas.

II. Objetivos
2.1 Objetivo Específico:
 Determinar la importancia de las redes de flujo y los métodos numéricos
aplicados a la mecánica de fluidos.
2.2 Objetivo General.
 Analizar y estudiar las redes de flujo y su influencia en las presas
hidráulicas.
 Estudiar y determinar las leyes de Darcy y de Laplace y la aplicación de
los métodos numéricos en las presas hidráulicas.

MECANICA DE FLUIDOS I 2
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

III. Justificación
 El presente informe se desarrolló para estudiar cual es aplicación de los
métodos numéricos en la Mecánica de Fluidos, además de la influencia y
el estudio de las redes de flujo en las obras hidráulicas.

IV. Métodos Numéricos Aplicados a la Mecánica de Fluidos

Los métodos numéricos permiten resolver de manera aproximada problemas definidos

mediante ecuaciones diferenciales. Es decir, aquellos en los que sabemos cómo cambian
las variables en el espacio y en el tiempo. Son precisamente las ecuaciones diferenciales
del problema las que expresan dichos cambios.

El objetivo es la determinación del valor de las variables en el dominio del problema. En

esta porción del espacio han de conocerse, por un lado, las ecuaciones diferenciales, que
establecen las pautas de variación de las variables en el espacio y en el tiempo. Y por el
otro lado, han de conocerse las condiciones de contorno y la condición inicial. Las
condiciones de contorno definen la relación de los valores de las variables en el interior
del dominio con los valores en el exterior de éste. La condición inicial define el estado de
partida del problema.

Para cada diferente condición inicial, el sistema evolucionará de manera distinta. Más
adelante se hablará de solución tentativa inicial, que es un concepto distinto. Por ejemplo,
la condición inicial podría estimarse a partir de medidas piezométricas en un determinado
instante, y a partir de esta, se podría calcular numéricamente la evolución del nivel
piezométrico en el un acuífero, o región de éste.

En definitiva, los métodos numéricos permiten resolver de manera aproximada problemas

definidos mediante ecuaciones diferenciales. La determinación del error de la solución es
esencial es este tipo de problemas.

MECANICA DE FLUIDOS I 3
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

V. Marco teórico.
5.1 Ley de Darcy.
Ecuación empírica que relaciona la pérdida de carga hidráulica (o pérdida de
presión) debido a la fricción a lo largo de una tubería dada con la velocidad media
del flujo del fluido.
Permite el cálculo de la pérdida de carga debida a la fricción dentro una tubería
llena.
Imagen N°1: Esquema del expe

rimento de Bernoulli

Fuente: Bejan A. (2005)

Daniel Bernoulli en 1856 describe el balance de energías a lo largo de una línea de

corriente, que para un fluido incompresible (densidad constante) y no viscoso, en régimen
laminar, puede escribirse mediante la ecuación:
𝑢 𝑉2
ℎ=𝑧+ + … … . (1.1)
𝛾𝑊 2𝑔

Donde:
g: aceleración de la gravedad.
z: Altura de posición (asociada a la posición del punto dentro de la línea de
corriente)
u/𝛾w: Altura asociada a la presión del fluido.
V2/2g: Altura asociada a la velocidad

MECANICA DE FLUIDOS I 4
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

Ésta es una ecuación homogénea en términos dimensionales ya que la ecuación de

dimensión de todos los términos es la de una longitud (m). En el flujo de aguas
subterráneas, los cambios de velocidad son despreciables por lo que la pérdida de energía
cinética también lo es y la ecuación anterior puede simplificarse a la forma:
𝑢
ℎ=𝑧+ … … … (1.2)
𝛾𝑊

La pérdida de carga total entre dos puntos cualesquiera debido al flujo de agua en suelo
saturado viene dada por:
∆ℎ = ℎ1 − ℎ2 … … . (1.3)

Donde:
∆ℎ : Pérdida de carga hidráulica, debido a la viscosidad del agua.
Con lo que, haciendo uso de la ecuación de Bernoulli, podemos escribir:
1
∆ℎ = 𝑧1 − 𝑧2 + (𝑢 − 𝑢1 ) … … … (1.4)
𝛾 1

Por otro lado, Darcy estableció empíricamente, en París alrededor de 1850, con un
dispositivo semejante a la siguiente relación constitutiva lineal entre el flujo de fluido o
caudal, q (m3/s), y el gradiente hidráulico, i (adimensional), que representa la tasa de
pérdida de carga.
𝑞 = 𝑘 𝑖 𝐴 … … . . (1.5)

Donde:
k: es la permeabilidad (m/s), llamada en muchos textos conductividad hidráulica del
medio poroso y A (m2) el área considerada. Se trata de una dependencia aproximada,
para flujos laminares y determinadas condiciones (Bejan, 2004). A través de la relación
entre el caudal, q, y la velocidad del fluido v (m/s), q/A=v, la expresión anterior puede
escribirse en la forma:
𝑣 = 𝑘 𝑖 … … . (1.6)

En términos del nivel piezométrico, el gradiente hidráulico se define como:

∆ℎ
𝑖= … … … (1.7)
∆𝑙

MECANICA DE FLUIDOS I 5
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

Así, en definitiva, la ley de Darcy afirma que la velocidad de descarga es directamente

proporcional al gradiente hidráulico, para medios de permeabilidad constante (medios
saturados en general) y régimen laminar o de pequeña velocidad.
∆ℎ
𝑣=𝑘 … … . . (1.8)
∆𝑙

Imagen N°2: Esquema del experimento de Darcy

Fuente: Bejan A. (2005)

MECANICA DE FLUIDOS I 6
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

5.2 Ley de Laplace.

Imagen N°3: Región de flujo (un elemento de suelo a través del que fluye agua)

Fuente: Bejan A. (2005)

Se supone que la velocidad v con que el agua pasa por el elemento posee tres componentes
Vx , Vy , Vz , y que éstas son sólo función de x, y, z, pero no del tiempo (puesto que por
hipótesis se trata de un régimen estacionario) y dichas componentes son funciones
continuas que admiten cualquier orden de derivación necesario al razonamiento expuesto.
En estas condiciones, si en las caras I de la Imagen, las componentes de la velocidad del
agua son Vx , Vy , Vz; en las caras II estas mismas componentes son:
𝜕𝑣𝑥
𝑣𝑥 + 𝑑𝑥 … … . . (2.1)
𝜕𝑥

𝜕𝑣𝑦
𝑣𝑦 + 𝑑𝑦 … … . . (2.2)
𝜕𝑦

𝜕𝑣𝑧
𝑣𝑥 + 𝑑𝑧 … … . . (2.3)
𝜕𝑧

MECANICA DE FLUIDOS I 7
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

Durante el flujo, la cantidad de agua que entra al elemento, en un régimen establecido,

tiene que ser igual a la que sale. Por lo tanto, teniendo en cuenta que el caudal que pasa
por una sección puede expresarse como el producto del área de la sección por la velocidad
del flujo puede escribirse:
𝑣𝑥 . 𝑑𝑦. 𝑑𝑧 + 𝑣𝑦 . 𝑑𝑥. 𝑑𝑧 + 𝑣𝑧 . 𝑑𝑥. 𝑑𝑦 =

𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧

(𝑣𝑥 + 𝑑𝑥) . 𝑑𝑦. 𝑑𝑧 + (𝑣𝑦 + 𝑑𝑦) . 𝑑𝑧. 𝑑𝑥 + (𝑣𝑧 + 𝑑𝑧) . 𝑑𝑦. 𝑑𝑥 … … . .2.4)
𝜕𝑥 𝜕𝑦 𝜕𝑧

En la expresión anterior, el primer miembro representa el caudal que entra al elemento y

el segundo el que sale. Reduciendo términos semejantes y reordenando:
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
( ) 𝑑𝑥. 𝑑𝑦. 𝑑𝑧 + ( ) 𝑑𝑦. 𝑑𝑧. 𝑑𝑥 + ( ) 𝑑𝑧. 𝑑𝑦. 𝑑𝑥 = 0 … … . . (2.5)
𝜕𝑥 𝜕𝑦 𝜕𝑧

De donde:
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
+ + = 0 … … . . (2.6)
𝜕𝑥 𝜕𝑦 𝜕𝑧

Esta Ec. (2.6) juega un papel importante en la teoría de flujo de agua y se conoce con el
nombre de ECUACIÓN DE CONTINUIDAD. Si también se supone válida la Ley de
Darcy, la velocidad de descarga a través del elemento es:
𝜕ℎ
𝑣 = 𝑘 ∗ 𝑖 = −𝑘 … … . . (2.7)
𝜕𝑆

Lo cual, expresando al gradiente hidráulico a través de sus tres componentes:

𝜕ℎ
𝑣𝑥 = −𝑘 … … . . (2.8)
𝜕𝑥

𝜕ℎ
𝑣𝑦 = −𝑘 … … . . (2.9)
𝜕𝑦

𝜕ℎ
𝑣𝑧 = −𝑘 … … . . (2.10)
𝜕𝑧

En este caso, el suelo se considera anisótropo en lo referente a su permeabilidad, en las

direcciones x, y, z. Introduciendo estas ecuaciones en la Ec. 2.6, se tiene:
𝜕 2ℎ 𝜕 2ℎ 𝜕 2ℎ
𝑘𝑥 + 𝑘𝑦 + 𝑘𝑧 = 0 … … . . (2.11)
𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2

Esta ecuación describe matemáticamente al flujo en la región considerada.

MECANICA DE FLUIDOS I 8
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

Una característica especial de los problemas en suelos donde se encuentra involucrada el

agua es la importancia que cobra los más pequeños detalles de la agua es la importancia
que cobra los más pequeños detalles de la estructura geológica. Así es como muchas veces
se dan como impermeables macizos rocosos o arcillosos que, sin embargo, luego dejan
pasar caudales considerables, por fisuras, grietas, canalículos o capas limosas difícilmente
discernibles para un observador descuidado. La trascendencia de estas heterogeneidades
se denota por la gran variedad de coeficientes k (parámetro que tiene mayor importancia
entre todos los que se manejan).
Si además, la masa de suelo en estudio, a través de la cual filtra el agua, es isótropa en lo
referente a permeabilidad, entonces:

k 𝑥 = k 𝑦 = k 𝑧 = 𝑘 … … . . (2.12)

Considerando esas hipótesis de la Ec. (2.11), puede simplificarse:

𝜕 2ℎ 𝜕 2ℎ 𝜕 2ℎ
+ + = ∇2 ℎ = 0 … … . . (2.13)
𝜕𝑋 2 𝜕𝑌 2 𝜕𝑍 2

La Ec. 2.13 se la conoce como la ECUACIÓN DE LAPLACE. En los problemas

prácticos de la Mecánica de Suelos, es muy frecuente que el flujo en una sección de la
región considerada, transversal a su eje longitudinal, sea idéntico al que se tiene en
cualquier otra sección. Éste es el caso, por ejemplo, en presas de tierra de eje largo, en
comparación con la altura. Así, los efectos en los bordes de la región de flujo pueden
ignorarse y, de esa manera, el problema de flujo puede estudiarse bidimensionalmente
como contenido todo él en el plano XY. En estas condiciones la Ec. 2.11 puede escribirse
en forma más simplificada como:

𝜕 2ℎ 𝜕 2ℎ
𝑘𝑥 + 𝑘𝑦 = 0 … … . . (2.14)
𝜕𝑥 2 𝜕𝑦 2

Que es la ecuación fundamental para el análisis de flujo bidimensional en una región dada.
La Ec. 2.14 puede simplificarse obteniéndose la Ec. 2.15 para representar
matemáticamente el problema suponiendo el suelo isótropo.

𝜕 2ℎ 𝜕 2ℎ
+ = ∇2 ℎ = 0 … … . . (2.15)
𝜕𝑥 2 𝜕𝑦 2

MECANICA DE FLUIDOS I 9
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

5.2.1 Solución de la Ecuación de Laplace.

Atendiéndose al caso del flujo bidimensional, se plantean dos grupos de funciones que
son solución de dicha ecuación y que tienen una interpretación geométrica muy útil:
 Función Potencial (Ф)
 Función de flujo (Ψ)
En primer caso se define una función Ф:

Ф = −k ∗ h + c … … . . (2.16)

La Ec. (2.16) satisface la ecuación de Laplace, por lo tanto se cumple:

𝜕 2Ф 𝜕 2Ф
+ = 0 … … . . (2.17)
𝜕𝑥 2 𝜕𝑦 2

Así, la función Ф(x, y) = 𝑐𝑡𝑒, es solución de la ecuación de Laplace. Esta solución

representa una infinidad de funciones según sea el valor de la constante c que intervenga.
De inmediato, puede darse una interpretación geométrica a esta solución, pues la
expresión Ф(x, y) = 𝑐𝑡𝑒 puede representar a una familia de curvas que se desarrollan en
la región plana en la que ocurre el flujo, obteniéndose una curva específica de la familia
para cada valor de la constante que se tome.
En cuanto a la interpretación física, se sigue que si una curva une puntos en que Φ es
constante, en esos puntos h también es constante. Es decir, en la curva Φ igual a una
constante, todos los puntos tienen la misma carga hidráulica h. Por esta razón, estas curvas
que unen puntos de igual carga hidráulica reciben el nombre de líneas equipotenciales.
Considérese ahora una función 𝛹(𝑥, 𝑦) = 𝑐𝑡𝑒, llamada función de flujo y definida de
modo que:
𝜕𝛹
𝑉𝑥 = … … . . (2.18)
𝜕𝑌

𝜕𝛹
𝑉𝑌 = … … . . (2.19)
𝜕𝑋

Puede demostrarse que una función ψ también cumple con la ecuación de Laplace y por
lo tanto es solución de la misma. Analizando el sentido físico de la curvas ψ = cte.
según la Fig. N°3:

MECANICA DE FLUIDOS I 10
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

Imagen N° 4: Trayectoria del Agua que pasa por un punto.

Fuente: Ms. Sc. Ing. Jorge B. Gutiérrez (2011)

Si se tiene en cuenta la trayectoria del agua que pasa por el punto P(x; y); el agua posee
una velocidad v, que es, naturalmente, tangente a su trayectoria. Si se trata ahora de
encontrar la ecuación de esa trayectoria a lo largo de la curva se tiene:
𝑉𝑦 𝑑𝑌
𝑡𝑔𝜃 = = … … . . (2.20)
𝑉𝑥 𝑑𝑋
De aquí:
𝑉𝑦 𝑑𝑋 − 𝑉𝑥 𝑑𝑌 = 0 … … . . (2.21)
De la expresión anterior y de las Ec. 2.18 y Ec. 2.19, se obtiene:
𝜕𝛹 𝜕𝛹
𝑑𝑋 + 𝑑𝑌 = 0 … … . . (2.22)
𝜕𝑋 𝜕𝑌
Finalmente tenemos:
𝜕 2𝛹 𝜕 2𝛹
+ =0
𝜕𝑥 2 𝜕𝑦 2

La solución de estas ecuaciones permite dibujar las iso-líneas de las correspondientes

magnitudes potenciales Φ y ψ a intervalos arbitrarios o, como es costumbre en general
para una inspección inmediata de los resultados, a intervalos constantes. Ambas funciones
potenciales se cortan formando siempre ángulos rectos en medios isótropos y
determinando sobre el plano un mallado formado por cuadrados curvilíneos denominado
red de flujo, como el de la Imagen N° 4.

MECANICA DE FLUIDOS I 11
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

Imagen N° 5: Líneas Equipotenciales y Líneas de Flujo

Fuente: Ms. Sc. Ing. Jorge B. Gutiérrez (2011)

5.2.2 Ecuación de Laplace en Mecánica de fluidos

En la mecánica de fluidos tenemos ciertas condiciones de borde usuales para la solución
de las ecuaciones que rigen el flujo subterráneo.
a) Condición de Potencial Impuesto
Es el caso que se presenta cuando el acuífero está en contacto con una masa libre de agua
tal como un rio o el mar. En esta situación la carga potencial es constante en todos los
puntos de la superficie de contacto entre el acuífero y el rio o entre el acuífero y el mar, y
está definida por la altura del agua en el rio o en el mar. En estas masas de agua las
pérdidas de carga son prácticamente despreciables pues si bien es cierto que la carga
pueda tener ciertas variaciones con el tiempo, dichas variaciones no dependen del
funcionamiento del acuífero sino de condiciones externas a del como lo son las
precipitaciones, por ejemplo. La condición se expresa por lo tanto como h = cte.
b) Condiciones de Flujo Impuesto
𝜕ℎ
Son equivalentes a la condición de Neumann, ya que si se impone un valor a (el
𝜕𝑛
𝜕𝐻 𝜕ℎ
gradiente en dirección n), se tiene a partir de la ley de Darcy: 𝑉𝑁 = −𝑘 𝜕𝑁 y como 𝜕𝑛 =

𝑐𝑡𝑒, entonces 𝑉𝑁 = 𝑐𝑡𝑒 y por consiguiente el caudal o flujo es también constante.

MECANICA DE FLUIDOS I 12
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

c) Condiciones de Fourier
Supónganse un rio que drena o alimenta un acuífero, que tiene un fondo colmatado por
un material poco permeable.
Imagen N° 6 Drenaje de un rio con suelo poco permeable

Fuente: Burden & Faires, (2008)

Imagen N°7. Drenaje de un acuífero

Fuente: Burden & Faires, (2008)

d) Condiciones de flujo a superficie libre

Dos condiciones definen una superficie libre:
- La presión sobre todo el punto M de la superficie libre es la presión atmosférica.
Se puede escribir entonces h = z.
- Además la superficie a flujo impuesto, que puede ser nulo si el acuífero no es
𝜕ℎ
alimentado por su superficie, o sea = 0, y si la napa es recargada por su
𝜕𝑛
𝜕ℎ
superficie, entonces = 𝑎. Esta “alimentación” puede ser también negativa,
𝜕𝑛

como en el caso en que haya evaporación.

MECANICA DE FLUIDOS I 13
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

- Aparece aquí entonces una doble determinación. EL problema principal reside en

el hecho de que la posición de la superficie libre no es conocida sino porque por
el contrario debe ser determinada y a su vez dicha superficie constituye una
condición de borde de flujo
Imagen. N°8: Superficie libre de un acuífero

Fuente: Burden & Faires, (2008)

Se trata entonces de una superficie que cumpla simultáneamente las dos

𝜕ℎ
ecuaciones h = z, y = 𝑐𝑡𝑒. Lo que se hace en la práctica es determinar la
𝜕𝑛

posición de la superficie por aproximaciones sucesivas. Primero se supone la

posición de la superficie, limitando asi el domino de integración, luego se fija la
𝜕ℎ
carga para dicho dominio h = z y se verifica que el caudal calculado 𝑘 𝜕𝑛, sea

correcto. Si dicho flujo no es correcto, se varía la posición de la superficie libre.

Hay condiciones de límites con flujo a superficie libre, por ejemplo en los
acuíferos libres en los cuales la superficie piezometrica es la misma superficie
freática. También en el caso del flujo a través de una presa de tierra, la línea de
saturación constituye un límite de flujo a superficie libre. En muchos casos la
superficie libre es cortada por una superficie que está en contacto con la atmosfera,
y aparece lo que se denomina una línea de emergencia del fluido, dejando de
existir una continuidad entre la superficie libre y el plano de agua hacia abajo.
Dicha superficie de contacto entre la superficie libre y la atmosfera es llamada
superficie de goteo. Como ejemplos de superficies de goteo se pueden precisar los
mostran flujo a través de una presa de tierra, flujo hacia un pozo, contacto de un
acuífero con una masa libre de agua. En la Imagen el sector AB es la superficie de
goteo.

MECANICA DE FLUIDOS I 14
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

En este caso entonces las condiciones de borde en la superficie de goteo se

expresan por la ecuación:
h=z
Imagen N°9: Línea de Flujo en una presa de tierra y línea de emergencia del fluido

Fuente: Burden & Faires, (2008)

Aquí también se presenta el problema de determinar la extensión de la superficie de goteo,

lo cual se hace también por aproximaciones sucesivas, como en el caso de la posición de
la superficie libre.
En ciertos casos, cuando se supone que el dominio de integración es infinito, es posible
abstraerse de las condiciones de frontera. Esto es muy utilizado cuando se están buscando
soluciones analíticas a la ecuación de difusión. Los métodos numéricos se adaptan mejor
cuando se tienen condiciones de frontera conocidas.
Para los problemas de flujo transitorio es necesario definir las condiciones iniciales del
problema o sea el valor de h en todo el dominio, cuando t=0.

MECANICA DE FLUIDOS I 15
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

 Ejemplo de Aplicación
Encontrar el caudal que fluye debajo de una presa que descansa sobre una fundación
permeable.

Solución:
Considerando el acuífero confinado y el flujo permanente se tiene:
∇2 ℎ = 0
Se considera además, que el flujo es unidimensional:
𝜕 2ℎ
=0
𝜕𝑥 2
Integrando esta última ecuación se tiene:
ℎ = 𝐶1 𝑥 + 𝐶2
Las condiciones de borde son:
Para x = 0, h =H1 y para x = B, h = H2 Esto implica que: C1 = H1 y C2 = H2
La cabeza piezometrica en cualquier punto debajo de la presa será entonces:
𝐻2 − 𝐻1
ℎ= 𝑥 + 𝐻1
𝐵
El caudal total, si L es la longitud de la presa, será:
𝑄 = 𝑉𝐴 → 𝐴 = 𝐿𝑐
𝑑ℎ 𝐻2 − 𝐻1
𝑉 = −𝐾 =𝐾
𝑑𝑥 𝐵
𝐿𝑐
𝑄 = 𝐾 (𝐻2 − 𝐻1 )
𝐵

MECANICA DE FLUIDOS I 16
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

VI. Redes de flujo

El empleo de la representación gráfica para el flujo que traspasa el suelo se lleva a cabo
a través de lo que se denomina como red de flujo o red de corriente. El camino seguido
por una partícula de agua en su recorrido de escurrimiento o filtración a lo largo de una
masa de suelo saturado es llamado línea de flujo o de corriente. Elementos de la red de
flujo:
Imagen N°10: partes de una red de flujo

Fuente: Ms. Sc. Ing. Jorge B. Gutiérrez (2011)

 Línea de Corriente:
También llamada Línea de flujo es la trayectoria seguida por las partículas de agua al fluir
a través del suelo.
 Línea Equipotencial:
Es aquella que une puntos en donde se tiene el mismo potencial hidráulico o carga
hidráulica.
 Tubo de Corriente:
Es el espacio comprendido entre líneas de corriente vecinas.
 Celda de Flujo:
Es el espacio comprendido entre dos líneas equipotenciales vecinas y dos líneas de
corriente vecinas.

MECANICA DE FLUIDOS I 17
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

Cada línea de flujo comienza en un punto localizado sobre la línea AB, con una carga de
presión (h = h1 – h2) que se disipa por la fricción viscosa, y termina en la línea CD, donde
la carga de presión es igual a h2. A todo lo largo de la línea existe un punto donde el agua
ha disipado una porción específica conocida como su potencial. Una línea que conecte
todos esos puntos de igual carga es llamada línea equipotencial.
Dichas líneas se asemejan a curvas de nivel de igual energía y las líneas de flujo las cortan
en ángulo recto, ya que el agua se mueve desde los niveles de mayor a los de menor
energía, siguiendo los caminos de gradiente de máxima energía; de la misma manera que
el agua corre hacia abajo en la ladera de una colina, de los niveles superiores a los
inferiores siguiendo la misma línea de pendiente máxima.
6.1 Propiedades de las líneas de corriente y equipotenciales

La cantidad de agua que pasa entre dos líneas de flujo es constante en cualquier sección
que se tome entre las líneas, dicho espacio se denomina canal de flujo. En
efecto:
Imagen N°11: líneas de equipotencial

Fuente: Orsolini H. (2009)

MECANICA DE FLUIDOS I 18
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

Imagen N°12: canal de flujo

Fuente: Orsolini H. (2009)

Donde q representa el caudal en el canal por unidad de longitud en dirección normal al

papel.
Las líneas de flujo no pueden cortarse dentro de la región de flujo. En efecto, si dos líneas
de flujo convergen en el punto de contacto no hay área para el paso del agua y en ese
punto no se respeta la continuidad del gasto, lo cual es imposible bajo las hipótesis de la
teoría de estudio. Las líneas equipotenciales tampoco pueden cortarse entre sí.

6.2 Trazado de la red de flujo:

El trazado de la red comprende en la práctica los siguientes pasos:

Delimitación de la zona de flujo que se desea estudiar, analizando sus condiciones
específicas de frontera.
Trazado de dos familias de curvas ortogonales entre sí que satisfagan las condiciones de
frontera y que constituyen una solución única de la ecuación de Laplace.

MECANICA DE FLUIDOS I 19
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

Imagen N° 13: Trazado de la red de flujo

Fuente: Orsolini H. (2009)

6.3 Condiciones hidráulicas de borde:

La determinación de las condiciones hidráulicas de borde consiste en establecer las
condiciones de carga hidráulica y circulación de flujo en todos los puntos límite del perfil
que se estudie. Estas condiciones, para el caso de la pantalla de tablestacas de la figura,
son definidas en forma completa por cuatro líneas significativas:
La línea AB, incluyendo su prolongación hacia la izquierda y hasta donde se extienda la
sección, es una línea equipotencial a lo largo de la cual la altura de agua es igual a “h1”
(primer línea equipotencial).
La línea CD incluyendo su prolongación hacia la derecha, es una línea equipotencial a lo
largo de la cual la altura del agua para este problema es igual a “h2“(última línea
equipotencial).
Desde el punto B (el cual toca la tablestaca en la superficie de la interfaz suelo agua); la
línea que sigue la superficie impermeable de la tablestaca hacia abajo hasta el punto E y
después hacia arriba por el otro lado de la tablestaca hasta el punto C (primer línea de
flujo).
La línea FG, incluyendo sus extensiones hacia ambos lados derecho e izquierdo (última
línea de flujo).

MECANICA DE FLUIDOS I 20
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

Imagen N°14: condiciones hidráulicas de borde

Fuente: Orsolini H. (2009)

6.4 Recomendaciones para el trazado de la red de flujo:

Al intentar el trazado de las familias de líneas equipotenciales y de flujo, surge el

problema de que por cada punto de la región debe pasar en principio precisamente una
línea de flujo y una equipotencial, pues en cada punto de la región de flujo, el agua tiene
una velocidad y una carga hidráulica. Esto llevaría, de trazar todas las líneas posibles, a
una solución que formaría una mancha uniforme, de modo que, al proceder así no se
tendría valor práctico, pues las soluciones encontradas en los diferentes problemas serían
uniformemente inútiles. Por lo tanto, es preciso no trazar todas las líneas de flujo y
equipotenciales posibles sino sólo algunas seleccionadas. El procedimiento más
conveniente es el siguiente:
 Dibujar las líneas de flujo de manera que el caudal que pase por el canal formado
entre cada dos de ellas sea el mismo (Δq).
 Dibujar las líneas equipotenciales de manera que la caída de carga hidráulica entre
cada dos de ellas sea la misma (Δh).
 Examinar las condiciones hidráulicas de borde del problema y determinar su
efecto sobre la forma de las líneas de corriente.

MECANICA DE FLUIDOS I 21
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

 Trazar todas las líneas de corriente recordando que deben empezar y terminar con
una tangente vertical, es decir, ser normales a las superficies del terreno, tanto
aguas abajo como aguas arriba. Dibujar líneas equipotenciales y de flujo siempre
completas.
 Obtener una primera aproximación dibujando varias curvas suaves que
representen líneas de corriente, que determine un número de canales de flujo entre
cuatro y cinco.
 Examinar cuidadosamente la red de filtración construida para identificar los
efectos más visibles, sin tratar de corregir detalles hasta que toda ella esté
aproximadamente bien trazada.
 Corregir los errores dibujando una nueva red. Proseguir con el procedimiento de
ajuste hasta obtener una red satisfactoria.
6.5 Cálculo del gasto o caudal en suelos isótropos:
Imagen N°15: líneas equipotenciales

Fuente: Orsolini H. (2009)

Como se dijo antes el espacio entre cualquier par de líneas de corriente o flujo se lo
denomina canal de flujo y se cumple que el caudal que circula a través de cada uno de
ellos es el mismo. Si se considera un tubo de corriente definido por dos líneas de corriente
(1-1 y 3-3), situadas en un mismo plano paralelo al movimiento (sin olvidar que el flujo
es bidimensional) y las líneas 2-2 y 4-4, correspondientes a aquellas situadas en un plano

MECANICA DE FLUIDOS I 22
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

paralelo al suyo a una distancia b, el agua entre las cuatro líneas se comporta
esencialmente como si circulara por una cañería limitada por ellas.

Dónde:
hi: altura piezométrica en cada una de las equipotenciales
Li: distancia entre las líneas equipotenciales consecutivas
Cuando la sección del suelo limitado por las cuatro líneas disminuye, como el caudal no
varía, la velocidad de corriente aumenta. Igualando los caudales en dos secciones, de áreas
A1 y A2, y teniendo en cuenta que en un suelo homogéneo el coeficiente de permeabilidad
es constante en toda la masa, tenemos:
𝑖1 ∗ 𝐴1 = 𝑖2 ∗ 𝐴2 … … … … … … … . (14.1)

h0 − h1
i1 = … … … … … … … … . (14.2)
L1

h1 − h2
i2 = … … … … … … … … . . (14.3)
L2

A1 = a1 ∗ b … … … … … … … … … .14.4)

𝐴2 = 𝑎2 ∗ 𝑏 … … … … … … … … … . (14.5)

Si se decide trazar las líneas equipotenciales de tal manera que las pérdidas de carga entre
cada par de líneas contiguas sean iguales:
∆ℎ = ℎ0 − ℎ1 = ℎ1 − ℎ2 … … … … … … … … … . (14.6)

Donde Δh es la caída de potencial hidráulico entre dos líneas equipotenciales,

sustituyendo los valores de las 1 a la 2 en la 6 obtenemos:
Δh ∗ a1 ∗ b Δh ∗ a2 ∗ b
= … … … … … … … … … … . (14.7)
L1 L2

a1 a2
= … … … … … … … … … … … … … … … … … … (14.8)
L1 L2

Esto significa que la relación de los lados de los rectángulos de una red de corriente,
formados por la intersección de líneas de flujo y equipotenciales, es constante. Así, si un
rectángulo de una red de corriente es aproximadamente cuadrado, a1 = L1, todos los otros
rectángulos de la red son aproximadamente cuadrados. El caudal Δq que pasa por el canal
vale, según la Ley de Darcy:
∆𝑞 = 𝑣 ∗ 𝐴 = 𝑣1 ∗ 𝑣2 = 𝑘 ∗ 𝑖 ∗ 𝑎 ∗ 𝑏 … … … . (14.9)

MECANICA DE FLUIDOS I 23
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

Si se supone un ancho unitario (b=1) y se calcula el caudal por ancho unitario, la ecuación
(9) se transforma en:
∆ℎ
∆𝑞 = 𝑘 ∗ 𝑖 ∗ 𝑎 = 𝑘 ∗ 𝑎 … … … … (14.10)
𝐿

Teniendo en cuenta las dos primeras recomendaciones que se siguen para construir la red
de flujo según el punto 2-3-2:
𝑞
∆𝑞 = … … … … … … … . . … … … (14.11)
𝑁𝑓

ℎ
∆ℎ = … … … … … … … … … … (14.12)
𝑁𝑐

Siendo:
Nf: cantidad de canales de flujo en la red
Nc: número de caídas de potencial
q: caudal por ancho unitario
h : pérdida de carga total
Entonces la Ec. 14.10 puede escribirse siendo a/l una constante:
𝑁𝑓
𝑞 =𝑘∗ℎ∗
𝑁𝑐

Dicha ecuación es la fórmula más sencilla que permite calcular el gasto ó caudal por
unidad de longitud normal a la sección estudiada que ocurre a través de una región de
flujo en la que se ha dibujado la red correspondiente.

Ejercicio de Aplicación

La presa del dibujo se asienta sobre materiales cuya conductividad hidráulica es 0,3 m/día.
Bajo dichos materiales se encuentra un sustrato impermeable.
Se pide:
a) Dibujar la red de flujo bajo la presa
b) Calcular el flujo por metro de presa (un metro perpendicular al dibujo)

MECANICA DE FLUIDOS I 24
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

Solución:
a) Dibujo de la red de flujo
Se han dibujado tres líneas de flujo, pero eso es subjetivo (red más densa, solución más
precisa, pero más difícil, o imposible, de realizar a mano).
La base de la presa y la formación impermeable inferior funcionan como líneas de flujo,
por tanto, las equipotenciales deben cortarlas perpendicularmente.

Fuente: Orsolini H. (2009)

b) Cálculo del caudal bajo la presa

Vamos a aplicar la Ley de Darcy a un tramo de presa de 1 metro. Calcularemos el caudal
para uno de los cuatro tubos de corriente, por ejemplo el que aparece punteado en la figura.

Fuente: Orsolini H. (2009)

MECANICA DE FLUIDOS I 25
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

Fuente: Orsolini H. (2009)

Consideremos ese tubo aislado y apreciamos que la sección es igual a: Sección a = 1m.
La distancia entre esa equipotencial y la siguiente es x y la diferencia de potencial entre
ellas será:
6 𝑚
∆𝑣 = = 0.4 𝑚𝑒𝑡𝑟𝑜𝑠
15 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑜𝑠

Por tanto, el gradiente hidráulico entre las dos superficies equipotenciales dibujadas en

la última figura será

∆ℎ/∆𝑥 = 0.4/∆𝑥

Finalmente aplicamos la Ley de Darcy:

∆ℎ
𝑄 = 𝑘 ∗ 𝑠𝑒𝑐𝑐𝑖𝑜𝑛 ∗ ( )
∆𝑥
0.3𝑚 0.4
𝑄= ∗ (𝑎 ∗ 1)𝑚2 ∗ ( ) = 0.12𝑚3 /𝑑𝑖𝑎
𝑑𝑖𝑎 ∆𝑥

Aunque desconocemos los valores de a y de x, se cancelan, ya que son iguales al ser la
red cuadrada.
Finalmente, multiplicamos por el número de tubos:
Q total = Q por cada tubo · nº de tubos = 0,12 · 4 = 0,48 m3/día
Este sería el caudal por cada metro, faltaría multiplicar por los metros de longitud de la
presa.

MECANICA DE FLUIDOS I 26
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

VII. Conclusiones.
 Se logró determinar la importancia de los métodos numéricos en la mecánica de
fluidos, simplificando las ecuaciones empíricas en ecuaciones aritméticas las cuales
nos hace más sencillo encontrar la solución de un problema tanto de flujo como el
caudal, la carga total y el gradiente hidráulico de una presa o de un suelo granular.
 Se determinó la importancia de las redes de flujo en los diferentes tipos de casos que
se puedan dar como presas concreto, de tierra, suelos, drenajes, etc.
 Se logró estudiar y determinar la importancia de las leyes de Darcy y Laplace y su
influencia en la determinación de las redes de flujo en las presas hidráulicas, así como
su empleo en el cálculo del caudal de flujo.
VIII. Referencias Bibliografía.
 Ms. Sc. Ing. Briones Jorge, (2011), red de flujo en secciones típicas presas de tierra
homogéneas.
 http://tesis.uson.mx/digital/tesis/docs/8530/Capitulo2.pdf
 Burden & Faires, “Análisis Numérico”. Ed. México: Thomson Internacional., 2008
 Bejan A. “Convection heat transfer”, 3ª ed. John Wiley & Sons, Hoboken, New
Jersey (2005)
 Orsolini, Hugo; Zimmermann, Eric; Basile, Pedro – (2009, 3ª Edición) – Hidrología:
Procesos y Métodos – Rosario: UNR Editora.

MECANICA DE FLUIDOS I 27
 UNIVERSIDAD NACIONAL DE CAJAMARCA
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICO PROFESIONAL DE INGENIERÍA CIVIL

MECANICA DE FLUIDOS I 28