Flujo por el Suelo
Storyboard 
En el caso del suelo, se puede modelar asumiendo que contiene múltiples poros que forman pequeños capilares que lo atraviesan. Con esta base, se pueden aplicar las ecuaciones para el flujo laminar a través de tubos y calcular las resistencias hidráulicas de las redes de capilares, que dependen de la porosidad y, por lo tanto, de la proporción de las distintas componentes.
ID:(370, 0)
Flujo por el Suelo
Storyboard 
En el caso del suelo, se puede modelar asumiendo que contiene múltiples poros que forman pequeños capilares que lo atraviesan. Con esta base, se pueden aplicar las ecuaciones para el flujo laminar a través de tubos y calcular las resistencias hidráulicas de las redes de capilares, que dependen de la porosidad y, por lo tanto, de la proporción de las distintas componentes.
Variables
Cálculos
Cálculos
Ecuaciones
Si consideramos el rea en la secci n que no contiene poros, restando la sección de poros ($S_p$) de la sección del flujo ($S$) y dividiendo por el rea de un grano gen rico de ERROR:10129,0, obtenemos el n mero de granos visibles en el corte:
$\displaystyle\frac{S-S_p}{\pi r_0^2}=\displaystyle\frac{(1-f)S}{\pi r_0^2}$
donde utilizamos la relaci n para la porosidad ($f$):
Si el n mero de granos es igual a ERROR:6039,0 con la expresi n:
donde el radio es el radio del tubo ($R$). Con esto, obtenemos la relaci n:
$\displaystyle\frac{(1-f)S}{\pi r_0^2}=\displaystyle\frac{fS}{\pi R^2}$
lo que resulta en:
Con la altura la distancia infinitesimal ($ds$), el volumen de la sección del flujo ($S$) es
$S ds$
y el de los poros con la sección de poros ($S_p$) es
$S_p ds$
por lo tanto, la porosidad ($f$) se calcula como
$f = \displaystyle\frac{S_p ds}{S ds} = \displaystyle\frac{S_p}{S}$
lo que resulta en la siguiente ecuaci n:
El volumen del capilar se puede calcular a partir de el radio del tubo ($R$) y la largo del capilar ($l$), lo cual es igual al volumen de una cadena de granos de el radio de un grano genérico ($r_0$) y el largo de tubo ($\Delta L$) multiplicado por la porosidad propia genérica ($q_0$):
$\pi R^2 l = q_0 \pi r_0^2 \Delta L$
Esto, en conjunto con la porosidad ($f$) en la relaci n
resulta en la siguiente relaci n:
Como la porosidad ($f$), calculado con la sección de poros ($S_p$) y la sección del flujo ($S$) usando
junto con la ecuaci n para el c lculo de la sección de poros ($S_p$) en funci n de el número de capilares ($N_p$) y el radio del tubo ($R$) mediante
con lo cual se obtiene
$f = \displaystyle\frac{N_p\pi R^2}{S}$
se puede despejar el número de capilares ($N_p$) resultando en
Para calcular el flujo total ($J_{Vt}$) utilizando el número de capilares ($N_p$) y el flujo de volumen ($J_V$) para cada capilar a trav s de
obtenemos el número de capilares ($N_p$) con la porosidad ($f$), la sección del flujo ($S$) y el radio del tubo ($R$) mediante
y la ley de Hagen-Poiseuille utilizando la viscosidad ($\eta$), la diferencia de presión ($\Delta p$) y la largo del capilar ($l$) se calcula con
$J_v = - \displaystyle\frac{\pi R^4}{8\eta}\displaystyle\frac{\Delta p}{l}$
Utilizando la relaci n para el radio del tubo ($R$) en t rminos de el radio de un grano genérico ($r_0$)
y para la largo del capilar ($l$), la porosidad propia genérica ($q_0$), y el largo de la muestra ($\Delta L$)
obtenemos
En el caso de los capilares por los cuales un l quido de la densidad del líquido ($\rho_w$) fluye debido a la diferencia de presión ($\Delta p$) generado por una diferencia de altura ($\Delta h$) bajo la influencia de la gravedad representada por la aceleración gravitacional ($g$) y calculado con la ecuaci n:
esto puede ser aplicado en la ecuaci n de Hagen-Poiseuille, en t rminos de el flujo total ($J_{Vt}$), que a su vez depende de el radio de un grano genérico ($r_0$), la porosidad propia genérica ($q_0$), la porosidad ($f$), la viscosidad ($\eta$), la sección o superficie ($S$) y el largo de la muestra ($\Delta L$) como se describe en la ecuaci n:
Junto con la definici n de la densidad de flujo ($j_s$):
$j_s = \displaystyle\frac{J_{Vt}}{S}$
Tenemos:
$j_s=\displaystyle\frac{J_{Vt}}{S}=\displaystyle\frac{ r_0 ^2 }{8 q_g }\displaystyle\frac{ f ^3 }{(1- f )^2}\displaystyle\frac{ \rho_w g }{ \eta }\displaystyle\frac{ \Delta h }{ L }$
lo que resulta en:
Dado que la densidad de flujo ($j_s$) est relacionado con el radio de un grano genérico ($r_0$), la porosidad ($f$), la densidad del líquido ($\rho_w$), la aceleración gravitacional ($g$), la viscosidad ($\eta$), la porosidad propia genérica ($q_0$), la diferencia de altura ($\Delta h$) y el largo de la muestra ($\Delta L$) a trav s de la ecuaci n
Podemos definir un factor al que llamaremos la conductividad hidráulica ($K_s$) de la siguiente manera:
Este factor incorpora todos los elementos relacionados con las propiedades del suelo y del l quido que fluye a trav s de l.
Usando la Ley de Darcy, donde la diferencia de presión ($\Delta p$) se iguala a la resistencia hidráulica ($R_h$) y el flujo total ($J_{Vt}$):
As , con la ecuaci n para el suelo con la sección del flujo ($S$), el radio de un grano genérico ($r_0$), la viscosidad ($\eta$), la porosidad propia genérica ($q_0$), la porosidad ($f$), la diferencia de presión ($\Delta p$) y el largo de la muestra ($\Delta L$):
Entonces, la resistencia hidráulica ($R_h$) es:
Ejemplos
Los poros individuales se agrupan en cadenas que forman capilares por los cuales el agua circula.
Para modelar este fen meno, es necesario estimar tanto el radio de estos capilares como su longitud, teniendo en cuenta que no suelen ser rectos.
Si se toma una sección del flujo ($S$) de una altura la distancia infinitesimal ($ds$)
se tendra que el volumen ser
$S ds$
meintras los poros ocuparan un volumen
$S_p ds$
lo que permite calcular la porosidad.
La sección del flujo ($S$) incluye la sección de poros ($S_p$) generada por el número de capilares ($N_p$):
Si observamos una secci n transversal del suelo, notaremos que los capilares atraviesan los espacios entre los granos. Al hacerlo, su n mero es similar al de los mismos granos, por lo que podemos suponer que el número de capilares ($N_p$) es similar al n mero de granos presentes en la secci n:
Si se compara el volumen del capilar con el de los granos del suelos se encontrara que en el suelo compactado los capilares son formados por la porosidad propia genérica ($q_0$) de los granos esta en esa misma relaci n con el volumen de los granos.
Por ello se puede establecer una relaci n entre ambos volumenes e igualados a la porosidad propia genérica ($q_0$) lo que permite generar una relaci n que permite calcular la largo del capilar ($l$).
El flujo total se calcula como la suma de los flujos individuales a trav s de los distintos poros:
Si asumimos que todos los poros son id nticos, podemos obtener el flujo total ($J_{Vt}$) multiplicando el flujo de volumen ($J_V$) individualmente por el número de capilares ($N_p$).
El flujo de l quido en un medio poroso como el suelo se mide mediante la variable la densidad de flujo ($j_s$), que representa la velocidad media a la que el l quido se desplaza a trav s de este medio. Al modelar el suelo y el paso del l quido a trav s de l, se descubre que este proceso est influenciado por factores como la porosidad ($f$) y el radio de un grano genérico ($r_0$), que, al ser mayores, facilitan el flujo, mientras que la viscosidad ($\eta$) dificulta el paso a trav s de los capilares, lo que reduce la velocidad de flujo.
El modelo finalmente incorpora lo que llamaremos la conductividad hidráulica ($K_s$), una variable que depende de las interacciones entre el radio de un grano genérico ($r_0$), la porosidad ($f$), la densidad del líquido ($\rho_w$), la aceleración gravitacional ($g$), la viscosidad ($\eta$) y la porosidad propia genérica ($q_0$):
la conductividad hidráulica ($K_s$) expresa la facilidad con la que el l quido se conduce a trav s del medio poroso. De hecho, la conductividad hidráulica ($K_s$) aumenta con la porosidad ($f$) y el radio de un grano genérico ($r_0$), y disminuye con la porosidad propia genérica ($q_0$) y la viscosidad ($\eta$).
Si examinamos la literatura, podemos encontrar estimaciones de la conductividad hidr ulica para diferentes texturas de suelo, que se muestran aqu en funci n de su exponente (es decir, se indica -7 para una conductividad hidr ulica de 1E-7 m/s):
Los resultados se resumen en la siguiente tabla:
| Textura | $g_a$ [%] | $g_i$ [%] | $g_c$ [%] | $f$ [%] | $K$ [m/s] |
| Arcilla | 0-45 | 0-40 | 55-100 | 40-50 | 1E-9 -1E-8 [1] |
| Franco | 23-52 | 28-50 | 8-27 | 40-50 | 1E-7 - 1E-5 [2] |
| Arena | 85-100 | 0-15 | 0-10 | 25-35 | 1E-4 - 1E-2 [3] |
| Limo | 0-20 | 80-100 | 0-13 | 35-45 | 1E-7 - 1E-5 [4] |
| Arcilla Limosa | 0-20 | 40-60 | 40-60 | 40-50 | 1E-9 - 1E-8 [1] |
| Arcilla Arenosa | 45-65 | 0-20 | 35-55 | 35-45 | 1E-7 - 1E-5 [5] |
| Franco Arcilloso | 20-45 | 15-53 | 28-40 | 40-50 | 1E-7 - 1E-5 [2] |
| Limo Arcilloso | 0-20 | 40-73 | 28-40 | 40-50 | 1E-8 - 1E-6 [6] |
| Franco Arcilloso Arenoso | 45-80 | 0-33 | 20-35 | 35-45 | 1E-6 - 1E-4 [1] |
| Franco Limoso | 0-50 | 50-88 | 0-28 | 35-45 | 1E-7 - 1E-5 [4] |
| Franco Arenoso | 43-85 | 0-50 | 0-20 | 30-40 | 1E-5 - 1E-3 [2] |
| Franco Arenoso Limoso | 70-90 | 0-30 | 0-15 | 25-35 | 1E-4 - 1E-2 [4] |
Estos datos se obtuvieron de la siguiente literatura, que se referencia en la columna de conductividad hidr ulica:
[1] "Geotechnical Engineering Principles and Practices" de Donald P. Coduto et al., Prentice Hall (1999)
[2] "Soil Mechanics and Foundations" de Muni Budhu, John Wiley & Sons. (2011)
[3] "Introduction to Environmental Engineering" de Mackenzie Davis y David Cornwell, McGraw Hill (2022)
[4] "Principles of Geotechnical Engineering" de Braja M. Das, CL-Engineering (2009)
[5] "Soil Mechanics in Engineering Practice" de Karl Terzaghi y Ralph B. Peck, John Wiley & Sons. (1996)
[6] "Soil Mechanics: Concepts and Applications" de William Powrie, CRC Press (2013)
La porosidad ($f$) se puede calcular a partir de la sección de poros ($S_p$) y la sección del flujo ($S$) utilizando la siguiente f rmula:
Dado que la secci n de un poro de ERROR:5417,0 es $\pi R^2$ y el número de capilares ($N_p$) est relacionado con eso, podemos calcular la sección de poros ($S_p$) de la siguiente manera:
Usando la porosidad ($f$) y la sección del flujo ($S$), obtenemos la sección de poros ($S_p$), que al dividirlo por la secci n del capilar calculada de el radio del tubo ($R$), resulta en el número de capilares ($N_p$) de acuerdo con:
Si asumimos que el n mero de capilares es igual al n mero de granos visibles en una secci n, podemos demostrar que para un radio de un grano de el radio de un grano genérico ($r_0$) y una porosidad ($f$), el radio del tubo ($R$) ser igual a:
Si igualamos el volumen de un capilar al volumen de una cadena de granos multiplicado por el largo de la muestra ($\Delta L$), obtenemos una relaci n entre la largo del capilar ($l$) y la porosidad ($f$) dada por:
El flujo total ($J_{Vt}$) se calcula multiplicando un número de capilares ($N_p$) por el valor de el flujo de volumen ($J_V$) en cada capilar, de la siguiente manera:
Si aplicamos la ecuaci n de Hagen-Poiseuille a el flujo de volumen ($J_V$) para el caso de capilares con el radio del tubo ($R$) expresados en t rminos de la porosidad ($f$) y la largo del capilar ($l$) como funci n de el largo de la muestra ($\Delta L$), podemos calcular el flujo total ($J_{Vt}$) utilizando
El resultado puede expresarse en t rminos de la sección del flujo ($S$), la viscosidad ($\eta$), la porosidad propia genérica ($q_0$), el radio de un grano genérico ($r_0$) y la diferencia de presión ($\Delta p$):
Esta ecuaci n corresponde a la ecuaci n de Kozeny-Carman, que fue desarrollada por Kozeny y Carman para modelar el flujo de un l quido a trav s de un medio poroso y se public en:
• Sobre la conducci n capilar del agua en el suelo., ("Ueber kapillare Leitung des Wassers im Boden."), J. Kozeny, Sitzungsber Akad. Wiss., Wien, 136(2a): 271-306 (1927)
• Flujo de fluidos a trav s de lechos granulares., ("Fluid flow through granular beds."), P.C. Carman, Transactions, Institution of Chemical Engineers, London, 15: 150-166 (1937)
• Flujo de gases a trav s de medios porosos., ("Flow of gases through porous media."), P.C. Carman, Butterworths, London (1956)
Como la conductancia hidráulica ($G_h$) est relacionado con la sección del flujo ($S$), el radio de un grano genérico ($r_0$), la viscosidad ($\eta$), la porosidad propia genérica ($q_0$), la porosidad ($f$) y el largo de la muestra ($\Delta L$), mediante
y dado que la conductancia hidráulica ($G_h$) es el inverso de la resistencia hidráulica ($R_h$), podemos afirmar que
En el caso de un tubo por el cual un l quido con la densidad del agua ($\rho_w$) fluye debido a la diferencia de presión ($\Delta p$) generado por una una diferencia de altura ($\Delta h$) bajo la influencia de la gravedad representada por la aceleración gravitacional ($g$) y calculado con la ecuaci n:
esto puede emplearse en la ecuaci n de Hagen-Poiseuille, junto con la definici n de la densidad de flujo ($j_s$) en t rminos de el flujo total ($J_{Vt}$), que a su vez depende de el radio de un grano genérico ($r_0$), la porosidad propia genérica ($q_0$), la porosidad ($f$), la viscosidad ($\eta$) y el largo de la muestra ($\Delta L$):
Para relacionar la conductividad hidr ulica con los factores de masa, introducimos el factor de escala de capilares ($\gamma$) junto con el radio de un grano genérico ($r_0$), el radio del grano de arena ($r_a$), la porosidad propia genérica ($q_0$) y la porosidad propia de la arena ($q_a$) como
A partir de la definici n de la conductividad hidráulica ($K_s$) y el factor de escala de capilares ($\gamma$), podemos expresar la conductividad hidráulica ($K_s$) en funci n de el radio del grano de arena ($r_a$), la porosidad propia de la arena ($q_a$), la porosidad ($f$), la densidad del líquido ($\rho_w$), la aceleración gravitacional ($g$) y la viscosidad ($\eta$) de la siguiente manera:
El c lculo de el factor de escala de capilares ($\gamma$) se obtiene a partir de la porosidad ($f$), la densidad del líquido ($\rho_w$), la aceleración gravitacional ($g$), y la viscosidad ($\eta$) eliminando el radio del grano de arena ($r_a$) y la porosidad propia de la arena ($q_a$) mediante la siguiente ecuaci n:
$\gamma=\displaystyle\frac{8K_sq_a}{r_a^2}\displaystyle\frac{(1-f)^2}{f^3}\displaystyle\frac{\eta}{\rho_wg}$
Si deseamos relacionar el factor de escala de capilares ($\gamma$) con la fracción de masa de arena en la muestra ($g_a$), la fracción de masa de limo en la muestra ($g_i$) y la fracción de masa de arcilla en la muestra ($g_c$), notamos que mientras el primero var a en 6 potencias, los segundos var an solo en 2, lo que justifica establecer una relaci n con el logaritmo natural de $\gamma$. Por lo tanto, realizamos una regresi n utilizando la ecuaci n:
con el factor de sección del capilar en arena ($s_a$), el factor de sección del capilar en limo ($s_i$), y el factor de sección del capilar en arcilla ($s_c$).
Los datos medios para cada rango son los siguientes:
| Tipo | $g_a$ [-] | $g_i$ [-] | $g_c$ [-] | $f$ [-] | $K_s$ [m/s] | $\gamma$ [-] | $\ln \gamma$ [-] |
| Arcilla | 0.20 | 0.20 | 0.60 | 0.45 | 1.00E-09 | 6.45E-10 | -21.16 |
| Marga | 0.40 | 0.40 | 0.20 | 0.45 | 1.00E-07 | 6.45E-08 | -16.56 |
| Arena | 0.93 | 0.03 | 0.04 | 0.30 | 1.00E-04 | 3.52E-04 | -7.95 |
| Limo | 0.10 | 0.85 | 0.05 | 0.40 | 1.00E-07 | 1.09E-07 | -16.03 |
| Arcilla limosa | 0.10 | 0.50 | 0.40 | 0.45 | 1.00E-09 | 6.45E-10 | -21.16 |
| Arcilla arenosa | 0.50 | 0.05 | 0.45 | 0.40 | 1.00E-07 | 1.09E-07 | -16.03 |
| Franco arcilloso | 0.30 | 0.35 | 0.35 | 0.45 | 1.00E-07 | 6.45E-08 | -16.56 |
| Franco arcilloso limoso | 0.10 | 0.55 | 0.35 | 0.45 | 1.00E-08 | 6.45E-09 | -18.86 |
| Franco arcilloso arenoso | 0.60 | 0.13 | 0.27 | 0.40 | 1.00E-06 | 1.09E-06 | -13.73 |
| Franco limoso | 0.20 | 0.65 | 0.15 | 0.40 | 1.00E-07 | 1.09E-07 | -16.03 |
| Franco arenoso | 0.65 | 0.25 | 0.10 | 0.35 | 1.00E-05 | 1.92E-05 | -10.86 |
| Arena arcillosa | 0.82 | 0.10 | 0.08 | 0.30 | 1.00E-04 | 3.52E-04 | -7.95 |
La regresi n resulta en una relaci n lineal con un R cuadrado de 0.9975 y los siguientes coeficientes, junto con los valores para evaluar la calidad de los coeficientes:
| Tipo | $s$ [-] | p-test |
| Arena (a) | -6.208 | 4.31E-6 |
| Limo (i) | -16.845 | 5.82E-9 |
| Arcilla (c) | -27.652 | 2.41E-9 |
Dado los valores del p-test, podemos asumir que todos los coeficientes son altamente relevantes.
ID:(370, 0)
