Flux au sol
Storyboard 
Dans le cas du sol, il est possible de le modéliser en supposant qu'il contient de multiples pores formant de petits capillaires qui le traversent. Sur cette base, il est possible d'appliquer les équations pour le flux laminaire à travers des tubes et de calculer les résistances hydrauliques des réseaux capillaires, qui dépendent de la porosité et donc de la proportion des différents composants.
ID:(370, 0)
Flux au sol
Storyboard 
Dans le cas du sol, il est possible de le modéliser en supposant qu'il contient de multiples pores formant de petits capillaires qui le traversent. Sur cette base, il est possible d'appliquer les équations pour le flux laminaire à travers des tubes et de calculer les résistances hydrauliques des réseaux capillaires, qui dépendent de la porosité et donc de la proportion des différents composants.
Variables
Calculs
Calculs
Équations
Si nous consid rons la zone dans la section qui ne contient pas de pores, en soustrayant a coupe des pores ($S_p$) de a coupe des pores ($S$) et en la divisant par la surface d'un grain g n rique avec un rayon de ERROR:10129,0, nous obtenons le nombre de grains visibles dans la section :
$\displaystyle\frac{S-S_p}{\pi r_0^2}=\displaystyle\frac{(1-f)S}{\pi r_0^2}$
o nous utilisons la relation pour a porosité ($f$) :
Si le nombre de grains est gal ERROR:6039,0 avec l'expression :
o le rayon est de le rayon du tube ($R$). Avec cela, nous obtenons la relation :
$\displaystyle\frac{(1-f)S}{\pi r_0^2}=\displaystyle\frac{fS}{\pi R^2}$
aboutissant :
Avec la hauteur a distance infinitésimale ($ds$), le volume de a coupe des pores ($S$) est
$S ds$
et le volume des pores avec a coupe des pores ($S_p$) est
$S_p ds$
Par cons quent, a porosité ($f$) est calcul comme
$f = \displaystyle\frac{S_p ds}{S ds} = \displaystyle\frac{S_p}{S}$
ce qui donne l' quation suivante :
Le volume du capillaire peut tre calcul partir de le rayon du tube ($R$) et a longueur capillaire ($l$), ce qui est gal au volume d'une cha ne de grains de le rayon d'un grain générique ($r_0$) et le longueur du tube ($\Delta L$) multipli par a porosité propre générique ($q_0$) :
$\pi R^2 l = q_0 \pi r_0^2 \Delta L$
Cela, en conjonction avec a porosité ($f$) dans la relation
donne la relation suivante :
Comme a porosité ($f$), calcul avec a coupe des pores ($S_p$) et a coupe des pores ($S$) en utilisant
ainsi que l' quation pour le calcul de a coupe des pores ($S_p$) en fonction de le nombre de capillaires ($N_p$) et le rayon du tube ($R$) travers
donne
$f = \displaystyle\frac{N_p\pi R^2}{S}$
peut tre r solu pour le nombre de capillaires ($N_p$), ce qui donne
Pour calculer le flux total ($J_{Vt}$) en utilisant le nombre de capillaires ($N_p$) et le volumique flux ($J_V$) pour chaque capillaire travers
nous obtenons le nombre de capillaires ($N_p$) avec a porosité ($f$), a coupe des pores ($S$) et le rayon du tube ($R$) travers
et la loi de Hagen-Poiseuille en utilisant a viscosité ($\eta$), a différence de pression ($\Delta p$) et a longueur capillaire ($l$) est calcul e avec
$J_v = - \displaystyle\frac{\pi R^4}{8\eta}\displaystyle\frac{\Delta p}{l}$
En utilisant la relation pour le rayon du tube ($R$) en termes de le rayon d'un grain générique ($r_0$)
et pour a longueur capillaire ($l$), a porosité propre générique ($q_0$) et le longueur de l'échantillon ($\Delta L$)
nous obtenons
Dans le cas des capillaires travers lesquels un liquide avec a densité du liquide ($\rho_w$) s' coule en raison de a différence de pression ($\Delta p$) g n r par une différence de hauteur ($\Delta h$) sous l'influence de la gravit repr sent e par a accélération gravitationnelle ($g$) et calcul avec l' quation :
cela peut tre appliqu dans l' quation de Hagen-Poiseuille, en termes de le flux total ($J_{Vt}$), qui d pend son tour de le rayon d'un grain générique ($r_0$), a porosité propre générique ($q_0$), a porosité ($f$), a viscosité ($\eta$), a coupe ou surface ($S$) et le longueur de l'échantillon ($\Delta L$) comme d crit dans l' quation :
Avec la d finition de a densité de flux ($j_s$) :
$j_s = \displaystyle\frac{J_{Vt}}{S}$
Nous avons :
$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 }$
ce qui donne :
tant donn que a densité de flux ($j_s$) est li le rayon d'un grain générique ($r_0$), a porosité ($f$), a densité du liquide ($\rho_w$), a accélération gravitationnelle ($g$), a viscosité ($\eta$), a porosité propre générique ($q_0$), a différence de hauteur ($\Delta h$) et le longueur de l'échantillon ($\Delta L$) travers l' quation :
Nous pouvons d finir un facteur que nous appellerons a conductivité hydraulique ($K_s$) comme suit :
Ce facteur englobe tous les l ments li s aux propri t s du sol et du liquide qui s' coule travers lui.
Avec la loi de Darcy, o a différence de pression ($\Delta p$) est gal a résistance hydraulique ($R_h$) et le flux total ($J_{Vt}$) :
Ainsi, avec l' quation pour le sol avec a coupe des pores ($S$), le rayon d'un grain générique ($r_0$), a viscosité ($\eta$), a porosité propre générique ($q_0$), a porosité ($f$), a différence de pression ($\Delta p$) et le longueur de l'échantillon ($\Delta L$) :
Par cons quent, a résistance hydraulique ($R_h$) est :
Exemples
Les pores individuels se regroupent pour former des cha nes qui cr ent des capillaires travers lesquels l'eau s' coule.
Pour mod liser ce ph nom ne, il est n cessaire d'estimer la fois le rayon de ces capillaires et leur longueur, en tenant compte du fait qu'ils ne sont g n ralement pas droits.
A coupe des pores ($S$) inclut a coupe des pores ($S_p$) g n r par le nombre de capillaires ($N_p$) :
Si nous observons une coupe transversale du sol, nous remarquerons que les capillaires passent travers les espaces entre les grains. En le faisant, leur nombre est similaire celui des grains eux-m mes, nous pouvons donc supposer que le nombre de capillaires ($N_p$) est similaire au nombre de grains pr sents dans la coupe :
Le d bit total est calcul comme la somme des d bits individuels travers les diff rents pores :
Si nous supposons que tous les pores sont identiques, nous pouvons obtenir le flux total ($J_{Vt}$) en multipliant le volumique flux ($J_V$) individuellement par le nombre de capillaires ($N_p$).
Le flux de liquide dans un milieu poreux comme le sol est mesur l'aide de la variable a densité de flux ($j_s$), qui repr sente la vitesse moyenne laquelle le liquide se d place travers celui-ci. Lors de la mod lisation du sol et de la mani re dont le liquide le traverse, on constate que ce processus est influenc par des facteurs tels que a porosité ($f$) et le rayon d'un grain générique ($r_0$), qui, lorsqu'ils sont plus lev s, facilitent le flux, tandis que a viscosité ($\eta$) entrave le passage travers les capillaires, r duisant ainsi la vitesse d' coulement.
Le mod le int gre finalement ce que nous appellerons a conductivité hydraulique ($K_s$), une variable qui d pend des interactions entre le rayon d'un grain générique ($r_0$), a porosité ($f$), a densité du liquide ($\rho_w$), a accélération gravitationnelle ($g$), a viscosité ($\eta$) et a porosité propre générique ($q_0$) :
a conductivité hydraulique ($K_s$) exprime la facilit avec laquelle le liquide est conduit travers le milieu poreux. En fait, a conductivité hydraulique ($K_s$) augmente avec a porosité ($f$) et le rayon d'un grain générique ($r_0$), et diminue avec a porosité propre générique ($q_0$) et a viscosité ($\eta$).
Si nous examinons la litt rature, nous pouvons trouver des estimations de la conductivit hydraulique pour diff rentes textures de sol, qui sont pr sent es ici en fonction de leur exposant (c'est- -dire, -7 est indiqu pour une conductivit hydraulique de 1E-7 m/s) :
Les r sultats sont r sum s dans le tableau suivant :
| Texture | $g_a$ [%] | $g_i$ [%] | $g_c$ [%] | $f$ [%] | $K$ [m/s] |
| Argile | 0-45 | 0-40 | 55-100 | 40-50 | 1E-9 - 1E-8 [1] |
| Limono-sableux | 23-52 | 28-50 | 8-27 | 40-50 | 1E-7 - 1E-5 [2] |
| Sable | 85-100 | 0-15 | 0-10 | 25-35 | 1E-4 - 1E-2 [3] |
| Limon | 0-20 | 80-100 | 0-13 | 35-45 | 1E-7 - 1E-5 [4] |
| Argilo-limoneux | 0-20 | 40-60 | 40-60 | 40-50 | 1E-9 - 1E-8 [1] |
| Argilo-sableux | 45-65 | 0-20 | 35-55 | 35-45 | 1E-7 - 1E-5 [5] |
| Limono-argileux | 20-45 | 15-53 | 28-40 | 40-50 | 1E-7 - 1E-5 [2] |
| Argilo-limoneux | 0-20 | 40-73 | 28-40 | 40-50 | 1E-8 - 1E-6 [6] |
| Sable limono-argileux | 45-80 | 0-33 | 20-35 | 35-45 | 1E-6 - 1E-4 [1] |
| Limono-argileux | 0-50 | 50-88 | 0-28 | 35-45 | 1E-7 - 1E-5 [4] |
| Sable limoneux | 43-85 | 0-50 | 0-20 | 30-40 | 1E-5 - 1E-3 [2] |
| Sable limono-sableux | 70-90 | 0-30 | 0-15 | 25-35 | 1E-4 - 1E-2 [4] |
Ces donn es ont t obtenues partir de la litt rature suivante, qui est r f renc e dans la colonne de la conductivit hydraulique :
[1] "Principes et Pratiques de G otechnique" de Donald P. Coduto et al., Prentice Hall (1999)
[2] "M canique des Sols et des Fondations" de Muni Budhu, John Wiley & Sons. (2011)
[3] "Introduction l'Ing nierie Environnementale" de Mackenzie Davis et David Cornwell, McGraw Hill (2022)
[4] "Principes de G otechnique" de Braja M. Das, CL-Engineering (2009)
[5] "M canique des Sols dans la Pratique de l'Ing nierie" de Karl Terzaghi et Ralph B. Peck, John Wiley & Sons. (1996)
[6] "M canique des Sols : Concepts et Applications" de William Powrie, CRC Press (2013)
A porosité ($f$) peut tre calcul partir de a coupe des pores ($S_p$) et a coupe des pores ($S$) en utilisant la formule suivante :
tant donn que la section transversale d'un pore de ERROR:5417,0 est $\pi R^2$ et que le nombre de capillaires ($N_p$) est li cela, nous pouvons calculer a coupe des pores ($S_p$) comme suit :
Avec a porosité ($f$) et a coupe des pores ($S$), nous obtenons a coupe des pores ($S_p$), que lorsqu'on divise par la section calcul e du capillaire de le rayon du tube ($R$), donne le nombre de capillaires ($N_p$) comme suit :
Si nous supposons que le nombre de capillaires est gal au nombre de grains visibles dans une section, nous pouvons d montrer que pour un rayon de grain de le rayon d'un grain générique ($r_0$) et une porosité ($f$), le rayon du tube ($R$) sera gal :
Si nous galons le volume d'un capillaire au volume d'une cha ne de grains multipli par le longueur de l'échantillon ($\Delta L$), nous obtenons une relation entre a longueur capillaire ($l$) et a porosité ($f$) donn e par :
Le flux total ($J_{Vt}$) est calcul en multipliant ERROR:6039.1 par la valeur de le volumique flux ($J_V$) dans chaque capillaire, comme suit :
Si nous appliquons l' quation de Hagen-Poiseuille le volumique flux ($J_V$) pour le cas des capillaires avec le rayon du tube ($R$) exprim s en termes de a porosité ($f$) et a longueur capillaire ($l$) en fonction de le longueur de l'échantillon ($\Delta L$), nous pouvons calculer le flux total ($J_{Vt}$) en utilisant
Le r sultat peut s'exprimer en termes de a coupe des pores ($S$), a viscosité ($\eta$), a porosité propre générique ($q_0$), le rayon d'un grain générique ($r_0$) et a différence de pression ($\Delta p$) :
Cette quation correspond l' quation de Kozeny-Carman, qui a t d velopp e par Kozeny et Carman pour mod liser l' coulement d'un liquide travers un milieu poreux et a t publi e dans:
• Sur la conduction capillaire de l'eau dans le sol, ("Ueber kapillare Leitung des Wassers im Boden"), J. Kozeny, Sitzungsber Akad. Wiss., Wien, 136(2a) : 271-306 (1927)
• coulement des fluides travers des lits granulaires, ("Fluid flow through granular beds"), P.C. Carman, Transactions, Institution of Chemical Engineers, London, 15 : 150-166 (1937)
• coulement des gaz travers des milieux poreux, ("Flow of gases through porous media"), P.C. Carman, Butterworths, Londres (1956)
Dans le cas d'un tube travers lequel un liquide avec a densité de l'eau ($\rho_w$) s' coule en raison de a différence de pression ($\Delta p$) g n r par une une différence de hauteur ($\Delta h$) sous l'influence de la gravit repr sent e par a accélération gravitationnelle ($g$) et calcul avec l' quation :
cela peut tre utilis dans l' quation de Hagen-Poiseuille, ainsi qu'avec la d finition de a densité de flux ($j_s$) en termes de le flux total ($J_{Vt}$), qui d pend son tour de le rayon d'un grain générique ($r_0$), a porosité propre générique ($q_0$), a porosité ($f$), a viscosité ($\eta$), et le longueur de l'échantillon ($\Delta L$) :
Pour tablir une relation entre la conductivit hydraulique et les facteurs de masse, nous introduisons le facteur d'échelle capillaire ($\gamma$) avec le rayon d'un grain générique ($r_0$), le rayon du grain de sable ($r_a$), a porosité propre générique ($q_0$) et a porosité du sable ($q_a$) comme
partir de la d finition de a conductivité hydraulique ($K_s$) et le facteur d'échelle capillaire ($\gamma$), nous pouvons exprimer a conductivité hydraulique ($K_s$) en fonction de le rayon du grain de sable ($r_a$), a porosité du sable ($q_a$), a porosité ($f$), a densité du liquide ($\rho_w$), a accélération gravitationnelle ($g$) et a viscosité ($\eta$) comme suit :
Le calcul de le facteur d'échelle capillaire ($\gamma$) est d riv de a porosité ($f$), a densité du liquide ($\rho_w$), a accélération gravitationnelle ($g$) et a viscosité ($\eta$), en excluant le rayon du grain de sable ($r_a$) et a porosité du sable ($q_a$), travers l' quation suivante :
$\gamma=\displaystyle\frac{8K_sq_a}{r_a^2}\displaystyle\frac{(1-f)^2}{f^3}\displaystyle\frac{\eta}{\rho_wg}$
Si nous voulons relier le facteur d'échelle capillaire ($\gamma$) A fraction massique de sable dans l'échantillon ($g_a$), a fraction massique de limon dans l'échantillon ($g_i$) et a fraction massique d'argile dans l'échantillon ($g_c$), nous observons que tandis que le premier varie de 6 ordres de grandeur, les derniers ne varient que de 2 ordres de grandeur. Il est donc logique d' tablir une relation avec le logarithme naturel de $\gamma$. Ainsi, nous effectuons une r gression en utilisant l' quation suivante :
avec le facteur de section capillaire dans le sable ($s_a$), le facteur de section capillaire dans le limon ($s_i$) et le facteur de section capillaire d'argile ($s_c$).
Les donn es moyennes pour chaque intervalle sont les suivantes :
| Type | $g_a$ [-] | $g_i$ [-] | $g_c$ [-] | $f$ [-] | $K_s$ [m/s] | $\gamma$ [-] | $\ln \gamma$ [-] |
| Argile | 0.20 | 0.20 | 0.60 | 0.45 | 1.00E-09 | 6.45E-10 | -21.16 |
| Limon | 0.40 | 0.40 | 0.20 | 0.45 | 1.00E-07 | 6.45E-08 | -16.56 |
| Sable | 0.93 | 0.03 | 0.04 | 0.30 | 1.00E-04 | 3.52E-04 | -7.95 |
| Silt | 0.10 | 0.85 | 0.05 | 0.40 | 1.00E-07 | 1.09E-07 | -16.03 |
| Argile limoneuse | 0.10 | 0.50 | 0.40 | 0.45 | 1.00E-09 | 6.45E-10 | -21.16 |
| Argile sableuse | 0.50 | 0.05 | 0.45 | 0.40 | 1.00E-07 | 1.09E-07 | -16.03 |
| Limono-argileux | 0.30 | 0.35 | 0.35 | 0.45 | 1.00E-07 | 6.45E-08 | -16.56 |
| Limono-argileux limoneux | 0.10 | 0.55 | 0.35 | 0.45 | 1.00E-08 | 6.45E-09 | -18.86 |
| Limono-argileux sableux | 0.60 | 0.13 | 0.27 | 0.40 | 1.00E-06 | 1.09E-06 | -13.73 |
| Limoneux | 0.20 | 0.65 | 0.15 | 0.40 | 1.00E-07 | 1.09E-07 | -16.03 |
| Sableux | 0.65 | 0.25 | 0.10 | 0.35 | 1.00E-05 | 1.92E-05 | -10.86 |
| Argilo-sableux | 0.82 | 0.10 | 0.08 | 0.30 | 1.00E-04 | 3.52E-04 | -7.95 |
La r gression donne une relation lin aire avec une valeur de R carr de 0,9975 et les coefficients suivants, ainsi que les valeurs pour valuer la qualit des coefficients :
| Type | $s$ [-] | p-test |
| Sable (a) | -6.208 | 4.31E-6 |
| Silt (i) | -16.845 | 5.82E-9 |
| Argile (c) | -27.652 | 2.41E-9 |
tant donn les valeurs du test de p, nous pouvons supposer que tous les coefficients sont hautement pertinents.
ID:(370, 0)
