Utilisateur:


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.

>Modèle

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

Symbole
Texte
Variable
Valeur
Unités
Calculer
Valor MKS
Unités MKS
$J_{Vt}$
J_Vt
Circulation à travers l'échantillon
m^3/s
$S_t$
S_t
Coupe de sol
m^2
$S_p$
S_p
Coupe des pores
m^2
$S$
S
Coupe ou surface
m^2
$J_{Vt}$
J_Vt
Débit selon la loi de Hagen Poiseuille pour les sols
m^3/s
$\rho_w$
rho_w
Densité du liquide
kg/m^3
$\Delta p_2$
Dp_2
Différence de pression 2
Pa
$J_V$
J_V
Flux Hagen-Poiseuille
m^3/s
$J_{Vt}$
J_Vt
Flux total
m^3/s
$\Delta h$
Dh
Hauteur de la colonne de liquide
m
$L$
L
Longueur capillaire
m
$N_p$
N_p
Nombre de capillaires
-
$f$
f
Porosité
-
$R$
R
Rayon de porosité
m
$S_1$
S_1
Section au point 1
m^2
$S_2$
S_2
Section au point 2
m^2
$\eta_w$
eta_w
Viscosité de l'eau
Pa s

Calculs


D'abord, sélectionnez l'équation:   à ,  puis, sélectionnez la variable:   à 

Symbole
Équation
Résolu
Traduit

Calculs

Symbole
Équation
Résolu
Traduit

 Variable   Donnée   Calculer   Cible :   Équation   À utiliser



É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$) :

equation=938

Si le nombre de grains est gal ERROR:6039,0 avec l'expression :

equation=4363

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 :

equation

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 :

equation

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

quation=109

donne la relation suivante :

quation

Comme a porosité ($f$), calcul avec a coupe des pores ($S_p$) et a coupe des pores ($S$) en utilisant

equation=938

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

equation=4362

donne

$f = \displaystyle\frac{N_p\pi R^2}{S}$



peut tre r solu pour le nombre de capillaires ($N_p$), ce qui donne

equation

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

equation=4364

nous obtenons le nombre de capillaires ($N_p$) avec a porosité ($f$), a coupe des pores ($S$) et le rayon du tube ($R$) travers

equation=4363

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$)

equation=109

et pour a longueur capillaire ($l$), a porosité propre générique ($q_0$) et le longueur de l'échantillon ($\Delta L$)

equation=2215

nous obtenons

equation

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 :

equation=4345

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 :

equation=4365

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 :

equation

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 :

equation=4366

Nous pouvons d finir un facteur que nous appellerons a conductivité hydraulique ($K_s$) comme suit :

equation

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}$) :

equation=3179

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$) :

equation=4365

Par cons quent, a résistance hydraulique ($R_h$) est :

kyon


Exemples


mechanisms

Les pores individuels se regroupent pour former des cha nes qui cr ent des capillaires travers lesquels l'eau s' coule.

image

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$) :

image

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 :

image

Le d bit total est calcul comme la somme des d bits individuels travers les diff rents pores :

image

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$) :

kyon

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) :

image

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)


model

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 :

image

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 :

image

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 :

kyon

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 :

kyon

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 :

kyon

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 :

kyon

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

equation=4364

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$) :

kyon

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 :

equation=4345

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$) :

kyon

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

kyon.

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 :

kyon.

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 :

kyon

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.


>Modèle

ID:(370, 0)



Mécanismes

Définition


ID:(15203, 0)



Couler à travers le sol

Image

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.

ID:(937, 0)



Échantillon de coupe de pores

Citation

A coupe des pores ($S$) inclut a coupe des pores ($S_p$) généré par le nombre de capillaires ($N_p$) :

ID:(2285, 0)



Relation entre le nombre de grains et de pores

Exercer

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 :

ID:(2283, 0)



Propre porosité et capillaires

Équation

ID:(2291, 0)



S'écoule à travers les pores individuels

Script

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$).

ID:(2286, 0)



Conductivité hydraulique pour différents sols

Variable

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)

ID:(4740, 0)



Modèle

Audio


ID:(15222, 0)