Flujo en Régimen Estacionario
Flujo en Régimen Estacionario
Flujo en Régimen Estacionario
+ = (6-6)
Donde:
h : carga total
z : carga de posicin o elevacin, altura por encima de un nivel de referencia arbitrario.
w
h
u
: carga de presin, altura de agua en un tubo vertical o piezmetro.
Ingeniera Geotcnica. 6. Flujo en rgimen estacionario
6-5
Figura 6.5 Carga de presin, carga de elevacin y cargas totales para el flujo de agua a travs de un
suelo. Tomado de: Braja Das, Fundamentos de Ingeniera Geotcnica, Thomson Learning, p.80.
Los tubos piezomtricos se utilizan para medir presiones estticas moderadas de un lquido que fluye
dentro de una tubera.
Otra ecuacin que es muy conocida de la mecnica de fluidos es la de Bernoulli para la energa o carga
total de un flujo estacionario incompresible:
g
u
z h
w
2
2
+ + = (6.7)
El ltimo trmino de la ecuacin es la carga de velocidad, que se desprecia porque la velocidad que tiene
el agua a travs de un medio poroso es muy pequea.
- Velocidad de infiltracin.
La velocidad de infiltracin, v
s
, es la tasa de movimiento de un elemento de agua a travs de la masa de
suelo. Esta tasa es importante en problemas geoambientales porque ayuda a determinar cuan
rpidamente viajan a travs del terreno los contaminantes.
En un principio, parecera que v
s
podra ser calculada con la ecuacin 6.3, pero esto sera incorrecto.
Dicha ecuacin describe el flujo a travs de tuberas, donde A es el rea de la seccin transversal de
flujo. Sin embargo la ecuacin 6.4 describe el flujo a travs de un medio poroso, donde A es el rea total
de la seccin transversal, de los vacos y de los slidos. El agua en los suelos solamente fluye a travs de
los vacos, as que estos dos valores de A son incompatibles (ver figura 6.6).
Esta diferencia puede ajustarse calculando la porosidad efectiva,
e
, la cual es el porcentaje de A en la
ecuacin 6.4 que contribuye al flujo, as:
e e
s
ki
= = (6.8)
Donde:
v
s
: velocidad de infiltracin.
v : velocidad de descarga.
k : coeficiente de permeabilidad o conductividad hidrulica.
i : gradiente hidrulico.
Ingeniera Geotcnica Patricia de Hasbun
6-6
Figura 6.6. Diagrama de fase para velocidad de flujo y velocidad superficial (el flujo es perpendicular a la
pgina). Tomado de: Holtz, R.D. y W.D. Kovacs (1981) An Introduction to Geotechnical Engineering.
Prentice Hall, p. 204
Sin embargo, el agua fluyendo a travs de un suelo toma rutas alrededor de las partculas de suelo, la
velocidad de infiltracin est basada en el movimiento sobre una lnea recta equivalente como se muestra
en la figura 6.7. No obstante, este mtodo produce una velocidad que es menor que la verdadera
velocidad (cuando es vista en una escala microscpica), y es ms conveniente para resolver problemas
de transporte de contaminantes.
Figura 6.7 La velocidad de infiltracin est basada en el movimiento, sobre una lnea recta equivalente,
del agua o contaminantes a travs del suelo.Tomado de: Coduto, D., Geotechinical Engineering. Prentice Hall.
p.233.
En suelos arenosos,
e
, es igual a la porosidad como se defini en el captulo 3. Sin embargo, los
suelos arcillosos contienen una capa esttica de agua alrededor de las partculas, de modo el rea de
flujo es menor que la el rea de vacos. Los valores de diseo de
e
para arcillas es mejor determinarlos
usando ensayos especiales de laboratorio (Kim, Edil, y Park, 1997). Si no hay disponibles datos de
ensayos, la Agencia de Proteccin Ambiental de los Estados Unidos usa
e
= 0.10 en arcillas.
Estudiar en el libro Fundamentos de Ingeniera Geotcnica de Braja Das, las pginas 82 y 83 donde se
demuestra que:
=
s
(6.9)
Ingeniera Geotcnica. 6. Flujo en rgimen estacionario
6-7
6.3 Coeficiente de permeabilidad.
El coeficiente de permeabilidad o conductividad hidrulica, k, para un lquido dado y un slido dependen
de muchos factores tales como:
- Propiedades del suelo:
o Tamao de vacos (depende del tamao de partcula, distribucin granulomtrica,
relacin de vacos y otros factores.
o Estructura del suelo.
o Continuidad de vacos.
o Forma y rugosidad de las partculas minerales.
- Propiedades del lquido.
o Densidad.
o Viscosidad.
La concentracin inica y el espesor de las capas de agua adheridas a las partculas, afectan a la
permeabilidad de las arcillas.
La tabla 6.1 presenta valores tpicos del coeficiente de permeabilidad para suelos saturados.
Tabla 6.1. Valores tpicos de k para suelos saturados. Tomada de Coduto, D., Geotechinical Engineering.
Prentice Hall. p.222.
La mayora de problemas prcticos tratan con agua limpia o agua contaminada con cantidades pequeas
de otras substancias. Cualquier variacin en la densidad y viscosidad es pequea y puede se ignorada
usualmente. De este modo, normalmente se piensa en que el coeficiente de permeabilidad depende
solamente del suelo.
Comnmente la unidad de medida de k es en cm/s. Sin embargo, otras unidades tambin son usadas,
incluyendo pies/min, pies/ao, etc. Aunque, las unidades de k, longitud/tiempo, son las mismas que las
que se utilizan para describir la velocidad, esta no es una medida de la velocidad.
Observe el extremadamente amplio rango de valores de k en la tabla 6.1, para arcillas generalmente se
tiene valores de k que son 1,000,000 veces ms pequeos que para las arenas, por tanto, basados en la
ecuacin 6.4, los valores de q tambin deberan ser 1,000,000 veces ms pequeos. El bajo valor de k en
las arcillas es debido al pequeo tamao de sus partculas (y por consiguiente al pequeo tamao de sus
vacos).
Ingeniera Geotcnica Patricia de Hasbun
6-8
6.4 Mtodos de medicin del coeficiente de permeabilidad.
Existen dos tipos de mtodos: directos e indirectos.
6.4.1 Permemetro de carga constante.
til para determinar k de suelos de granos gruesos, con valores de k mayores a 10
-4
m/s.
La muestra de suelo se introduce en un cilindro de plstico transparente de metacrilato de metilo con
filtros de malla de alambre y grava por encima y por debajo, ver figura 6.7. En la pared lateral del cilindro
existen varios puntos de conexin de manmetros para tomar diferentes lecturas de pares de cargas de
presin (en el diagrama slo se muestra uno de estos pares). EI agua que f1uye a travs de la muestra
proviene de un tanque o depsito diseado para mantener una carga constante, y la cantidad de agua se
mide pesando el recipiente recolector.
Puesto que la presencia de burbujas puede afectar mucho los resultados, es necesario asegurarse de
que el aire se haya eliminado del sistema en la medida de lo posible. Esto se lleva a cabo, por una parte,
suministrando agua desairada al tanque de carga constante y, por la otra, aplicando un vaco a la muestra
antes de iniciar la prueba. La prueba se inicia con las vlvulas A y B abiertas y la C cerrada, usndose la
vlvula A para controlar la velocidad de flujo. EI f1ujo se contina hasta lograr un estado estable o de flujo
establecido, esto es, hasta que los niveles en los tubos de los manmetros sean constantes. Una vez que
se haya alcanzado el estado estable, se mide la cantidad que fluye durante un tiempo dado y se registran
las lecturas de los dos niveles manomtricos. Despus se modifica la velocidad de flujo y se repite el
procedimiento. Se deben efectuar varias pruebas con velocidades de flujo y cargas diferentes para
establecer el valor promedio de k.
Para obtener el valor de k se aplica la ley de Darcy: kiA q =
t h A
L V
A h
L
t
V
Ai
q
k = = =
1
* * (6.10)
t
V
q
L
h
i = = ;
Donde:
V = cantidad de agua recolectada en el tiempo T (s), 1 ml = 10
3
mm
3
A = rea de la seccin transversal de la muestra (mm
2
).
h = diferencia de niveles de los manmetros (mm).
L = distancia entre los puntos de conexin de los manmetros (mm).
Ingeniera Geotcnica. 6. Flujo en rgimen estacionario
6-9
Figura 6.7 Permemetro de carga constante.
6.4.2 Permemetro de carga variable.
La prueba de carga variable se usa para determinar el coeficiente de permeabilidad de suelos de grano
fino, tales como arenas finas, limos y arcillas. Para estos suelos, el ritmo de flujo que los atraviesa es
demasiado pequeo para permitir mediciones precisas con el permemetro de carga constante.
En un cilindro de 100 mm de dimetro (generalmente) se introduce una muestra inalterada del suelo; el
cilindro puede ser un tubo U100 de muestreo, o bien el tubo corta ncleos que se usa en los ensayos de
densidad de campo. Las muestras tambin pueden prepararse por compactacin en un molde estndar.
Los extremos superior e inferior de la muestra se protegen con un filtro de malla de alambre y grava. La
base del cilindro queda sumergida en un recipiente con agua que cuenta con vertedor de nivel constante,
y el extremo superior se conecta a un tubo piezomtrico abierto de dimetro conocido. (Ver figura 6.8).
La prueba se lleva a cabo llenando la columna con agua desaireada y permitiendo que se produzca la
infiltracin a travs de la muestra. Se registra la altura de agua en la columna a diferentes tiempos
durante la prueba, y esta se repite con tubos piezomtricos de diferente dimetro. Despus de calcular
los resultados se determina el valor de k. Por lo general se reportan tambin los pesos unitarios inicial y
final y el contenido de humedad de la muestra.
Ingeniera Geotcnica Patricia de Hasbun
6-10
Figura 6.8 Permemetro de carga variable.
Partiendo de la ley de Darcy: kiA q =
Con respecto a la figura 6.8, si el nivel de la columna se reduce dh en un tiempo dt, se tiene:
dt
dh
a q = (6.11)
Y el gradiente hidrulico:
L
h
i =
Por lo que:
L
h
A k
dt
dh
a q = = (6.12)
Donde:
a = rea de la seccin transversal de la columna de agua.
A = rea de la seccin transversal de la muestra.
Reordenando e integrando:
=
2
1
2
1
t
t
h
h
dt
L
A
a
k
h
dh
Ingeniera Geotcnica. 6. Flujo en rgimen estacionario
6-11
( )
1 2
1
2
log t t
L
A
a
k
h
h
e
=
Despejando para k:
( )
( )
1 2
2 1
/ log
t t A
h h L a
k
e
= (6.13)
( )
( )
1 2
2 1 10
/ log 3 . 2
t t A
h h L a
k
= (6.14)
Estudiar el apartado 4.5 del libro Fundamentos de Ingeniera Geotcnica, de Braja Das, sobre
Relaciones empricas para la permeabilidad.
6.5 Infiltracin y redes de flujo: flujo bidimensional.
La infiltracin que se presenta alrededor de "tablestacas", presas y otras estructuras para la retencin de
agua, as como a travs de terraplenes y diques es bidimensional. Es decir, las componentes horizontal y
vertical de la velocidad varan de un punto a otro en la seccin transversal de la masa de suelo.
En el curso se considerar el caso general del flujo bidimensional en una masa de suelo homognea e
isotrpica (esto es k
v
= k
h
) para despus pasar ala representacin grfica conocida como red de flujo.
Cuando el agua fluye a travs de un medio poroso, se pierde energa o carga debido a la friccin, similar
a lo que pasa en flujo a travs de tuberas o canales abiertos.
En el ensayo de permeabilidad del laboratorio al igual que en una presa de tierra o en un tablestacado, se
generan prdidas de carga o energa.
"
Obsrvese como se reducen los niveles de agua en cada piezmetro a medida que fluye el agua a travs
del tablestacado hasta la lnea base agua abajo, ver figura 6.9.b).
Podramos representar el flujo a travs de la presa de la figura 6.10 por medio de lneas de flujo, que
seran la trayectoria promedio del flujo de una partcula de agua que fluye del reservorio aguas arriba al
nivel de descarga.
Podramos tambin representar la energa de flujo usando lneas de igual potencial, llamadas
naturalmente lneas equipotenciales. A lo largo de cada lnea equipotencial la energa disponible para el
flujo es la misma; a la inversa, la energa prdida por el agua que llega a dicha lnea es la misma a lo
largo de la lnea.
La red de lneas de flujo y lneas equipotenciales es llamada red de flujo, ilustran grficamente cmo la
carga o energa se pierde a medida que el agua fluye a travs de un medio poroso.
Ingeniera Geotcnica Patricia de Hasbun
6-12
Figura 6.9.a) Tablestacado para retencin de agua, b) Piezmetros y red de flujo.
Se considerar en el curso el caso de flujo confinado, es decir, donde la infiltracin est confinada entre
dos superficies impermeables. (Ver figuras 6.9.a y 6.10).
Casos de flujos no confinados son el de presas de tierra como la mostrada en la figura 6.11.b y la
infiltracin hacia pozos.
Las redes de flujo son utilizadas para:
Estimacin de prdidas por infiltracin de reservorios.
Determinacin de subpresiones bajo presas.
Revisin de puntos de potencial erosivo, donde i tienda a i
cr.
.
Una red de flujo es realmente la solucin grfica de la ecuacin de Laplace en dos dimensiones:
0
2
2
2
2
=
c
c
+
c
c
z
h
x
h
(6.15)
Donde x y z son las dos direcciones coordenadas.
Ingeniera Geotcnica. 6. Flujo en rgimen estacionario
6-13
La ecuacin de Laplace es una ecuacin muy importante en la fsica y la matemtica, porque representa
la prdida de energa a travs de un medio resistivo. Por ejemplo: puede describir el flujo de electrones, el
flujo de personas hacia hospitales, etc.
Si las condiciones de frontera (geometra, condiciones de flujo y condiciones de energa (carga) en la
frontera) son simples, es posible resolver la ecuacin de manera exacta (cerrada).
Figura 6.10 Red de flujo de la infiltracin por debajo de un vertedero.
Tomada de Coduto, D., Geotechinical Engineering. Prentice Hall. p.244
Reglas de construccin y condiciones lmites de una red de flujo.
Las redes de flujo son construcciones diafragmticas representativas de condiciones de un flujo
bidimensional. Una vez que se ha trazado la red de flujo, es posible evaluar cantidades como el grado de
infiltracin y la carga de presin. Para poder construir correctamente una red de flujo es necesario aplicar
varias reglas (ver figura 6.11):
Figura 6.11 Redes de flujo en: a) tablestacas, y b) presas de tierra.
Campos cuadrados: las reas limitadas por equipotenciales y lneas de corriente deben ser tan cuadradas
como resulte posible.
Intersecciones a ngulos rectos: la interseccin de una equipotencial con una lnea de corriente debe ser
de 90.
Ingeniera Geotcnica Patricia de Hasbun
6-14
Frontera impermeable: puesto que no puede haber flujo a travs de una superficie impermeable, este
lmite es una lnea de corriente (AB y DEF).
Frontera permeable: un lmite permeable sumergido en el cual la carga sea constante constituye una
equipotencial (CD y FG).
Superficie fretica: a lo largo de una superficie fretica (PQ), que tambin se conoce como lnea de flujo
superior, la presin de poro u es igual a cero, de modo que es una lnea de f1ujo.
Superficie de infiltracin: la superficie de infiltracin aparece cuando la superficie piezomtrica intersecta
tangencialmente a la superficie del terreno (QR); tiene las mismas caractersticas de frontera que la
superficie piezomtrica.
Pasos para construir una red de flujo:
1. Trazar la seccin transversal a escala reducida, definiendo todas las fronteras del lugar, la estructura,
etc.
2. Dibujar a lpiz unas cuantas lneas de corriente y equipotenciales de prueba, siguiendo las reglas
explicadas anteriormente. A medida que se aaden ms campos "cuadrados" comienza a tomar forma
la red de flujo. Es importante realizar una serie de ensayos aplicando un buen criterio, un lpiz y un
borrador, con el objetivo final de obtener una red de flujo en que todas las partes cumplan con las
reglas. Se requiere bastante prctica, por lo que conviene trazar tantas redes de flujo como sea posible
para lograr la experiencia adecuada. En teora se pueden trazar tantas lneas de flujo como se desee y
mientras mayor sea su nmero, mejor ser la precisin de los clculos subsecuentes. Sin embargo,
desde un punto de vista prctico, la tarea se simplifica al trazar slo unas cuantas lneas de flujo; casi
nunca se requieren ms de cinco o seis.
No todos los cuadrados deben ser del mismo tamao y no son "cuadrados" en el sentido estricto de la
palabra, se consideran "cuadrados" cuando pueden subdividirse en formas que sean verdaderamente
equilteras.
Ni el nmero de canales de flujo (canales entre lneas de flujo), ni el nmero de cadas de potencial (una
cada es la reduccin en energa o carga h de una lnea equipotencial a la siguiente) tienen que ser un
nmero entero, tambin se permiten cuadrados fraccionales.
En la figura 6.12 se muestra una seccin de una red de flujo, en ella se puede observar cmo entre dos
lneas equipotenciales separadas una distancia l se produce una cada de potencial h. Por otra parte,
se advierte que entre dos lneas de flujo se produce un gasto q.
De la figura 6.12 se puede plantear que el gradiente hidrulico es:
b
N
h
b
h
L
h
i
d
l
=
A
=
A
A
= (6.16)
Donde : b = longitud de trayectoria.
h = cada de potencial entre dos equipotenciales.
h
l
= energa (carga) total prdida en el sistema.
N
d
= nmero total de cadas de potencial.
Ingeniera Geotcnica. 6. Flujo en rgimen estacionario
6-15
Figura 6.12 Interseccin de lneas de flujo y lneas equipotenciales.
De la ley de Darcy:
a
b
N
h
k A
l
h
k q A i k q
d
l
=
A
= A = (6.17)
Y la descarga total q por unidad de profundidad (perpendicular al papel) es:
d
f
l f
N
N
b
a
h k N q q = A = (6.18)
Donde N
f
es nmero total de canales de flujo en la red.
Si esbozamos o dibujamos cuadrados a = b:
d
f
l
N
N
h k q = (6.19)
Con una red de flujo gruesa puede hacerse una estimacin bastante exacta de las cantidades de flujo.
Si a
1
/b
1
= a
2
/b
2
= a
3
/b
3
= a
n
/b
n
(los elementos no son cuadrados).
n
N
N
h k q
d
f
l
= (6.20)
Las subpresiones son las presiones en el agua en los puntos debajo de una estructura. Su clculo es
importante en el diseo de la estabilidad de una estructura como un conjunto.
Ingeniera Geotcnica Patricia de Hasbun
6-16
Figura 6.13 Subpresiones en una presa.
6.6 Gradiente hidrulico crtico.
Cuando el agua fluye a travs de una masa de suelo, la resistencia debida a la viscosidad en los canales
formados por los poros produce unas fuerzas de filtracin que el agua transmite a las partculas de suelo.
En los puntos donde predomina el flujo ascendente, estas fuerzas de filtracin tienden a disminuir el
esfuerzo efectivo entre las partculas del suelo, y por tanto tienden a reducir la resistencia al esfuerzo
cortante de la masa de suelo.
Esto puede demostrarse con el siguiente ejemplo, considere una muestra de arena sometida a una carga
h, (ver figura 6.14):
Figura 6.14 Muestra de arena sometida a una carga Ah.
Ingeniera Geotcnica. 6. Flujo en rgimen estacionario
6-17
El gradiente hidrulico a travs de la arena est dado por:
D
h
i
A
A
=
En la base de la arena el esfuerzo vertical total est dado por:
D
m
A =
La presin intersticial u asociada es:
( ) h D u
w
A + A =
Si se aplica el principio de los esfuerzos efectivos:
( ) h D D u
w m
A + A A = = '
h D D
w w m
A A A = '
Como = =
w m
' peso volumtrico sumergido al sustituir en la ecuacin anterior se obtiene:
h D
w m
A A = ' '
|
|
.
|
\
|
A
A
A =
D
h
D
m
w
m
'
1 ' '
|
|
.
|
\
|
A = i D
m
w
m
'
1 ' '
(6.21)
Al elevar el recipiente con agua, el gradiente hidrulico i a travs de la arena aumenta, y en la ecuacin
6.21 se observa que el esfuerzo vertical efectivo disminuye.
Esta ecuacin indica que el esfuerzo vertical se reduce a cero cuando el gradiente hidrulico alcanza el
valor
m
/
w
. Este valor se denomina gradiente hidrulico crtico, i
c
.
w
m
c
i
'
=
(6.22)
En la tabla 6.2 se presentan rangos valores de gradientes hidrulicos crticos en arenas, en ella puede
observarse que el valor del gradiente hidrulico es aproximadamente igual a uno.
Tabla 6.2 Gradientes hidrulicos en arenas
Suelta Densa
Arena uniforme 0.89 1.08
Arena bien graduada 0.99 1.15
Gradacin
Densidad
Ingeniera Geotcnica Patricia de Hasbun
6-18
En la condicin crtica cuando el esfuerzo efectivo es nulo, las partculas de arena se separan unas de
otras y se presentan como en una suspensin en el agua intersticial; la masa de suelo pareciera en ese
instante que se encontrara en ebullicin. Se dice entonces que la arena est en una condicin de
licuefaccin en la cual su resistencia al esfuerzo cortante es nula y por tanto en altamente inestable.
La licuefaccin generalmente se produce en los limos y en las arenas finas y medias. En los suelos
arcillosos la adherencia entre las partculas de mineral de arcilla ayuda a evitar su separacin; y en los
suelos de alta permeabilidad como arenas gruesas y gravas es improbable disponer de grandes
volmenes de agua necesarios para mantener la condicin de licuefaccin que pudiera producirse.