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)



Mécanismes

Iframe

>Top



Code
Concept

Mécanismes

ID:(15203, 0)



Couler à travers le sol

Concept

>Top


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)



Rapport volumique de porosité

Concept

>Top


ID:(2284, 0)



Échantillon de coupe de pores

Concept

>Top


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

Concept

>Top


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

Concept

>Top


ID:(2291, 0)



S'écoule à travers les pores individuels

Concept

>Top


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 du sol

Équation

>Top, >Modèle


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

$ K_s \equiv \displaystyle\frac{ r_0 ^2 }{8 q_0 }\displaystyle\frac{ f ^3 }{(1- f )^2}\displaystyle\frac{ \rho_w g }{ \eta }$

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

$ j_s =-\displaystyle\frac{ r_0 ^2 }{8 q_0 }\displaystyle\frac{ f ^3 }{(1- f )^2}\displaystyle\frac{ \rho_w g }{ \eta }\displaystyle\frac{ \Delta h }{ \Delta L }$



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

$ K_s \equiv \displaystyle\frac{ r_0 ^2 }{8 q_0 }\displaystyle\frac{ f ^3 }{(1- f )^2}\displaystyle\frac{ \rho_w g }{ \eta }$

Ce facteur englobe tous les éléments liés aux propriétés du sol et du liquide qui s'écoule à travers lui.



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

ID:(4739, 0)



Conductivité hydraulique pour différents sols

Image

>Top


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

Top

>Top



Paramètres

Symbole
Texte
Variable
Valeur
Unités
Calculer
Valor MKS
Unités MKS
$\rho_w$
rho_w
Densité du liquide
kg/m^3
$\pi$
pi
Pi
rad
$\eta_w$
eta_w
Viscosité de l'eau
Pa s

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

Calculs


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

Calculs

Symbole
Équation
Résolu
Traduit

Calculs

Symbole
Équation
Résolu
Traduit

Variable Donnée Calculer Cible : Équation À utiliser




Équations

#
Équation

$ f = \displaystyle\frac{ S_p }{ S }$

f = S_p / S


$ \gamma = \displaystyle\frac{ r_0 ^2}{ r_a ^2}\displaystyle\frac{ q_a }{ q_0 }$

g = r_0 ^2 * q_a /( r_a ^2* q_0 )


$ j_s =-\displaystyle\frac{ r_0 ^2 }{8 q_0 }\displaystyle\frac{ f ^3 }{(1- f )^2}\displaystyle\frac{ \rho_w g }{ \eta }\displaystyle\frac{ \Delta h }{ \Delta L }$

j_s = - r_0 ^2* rho_w * g * f ^3 * Dh /(8* q_0 * eta * DL )


$ J_{Vt} =-\displaystyle\frac{ r_0 ^2}{8 \eta q_0 }\displaystyle\frac{ f ^3}{(1- f )^2}\displaystyle\frac{ S }{ \Delta L } \Delta p $

J_Vt = - S * r_0 ^2* f ^3 * Dp /(8* eta * q_0 *(1- f )^2* DL )


$ J_{Vt} = N_p J_V $

J_{Vt} = N_p * J_V


$ K_s \equiv \displaystyle\frac{ r_0 ^2 }{8 q_0 }\displaystyle\frac{ f ^3 }{(1- f )^2}\displaystyle\frac{ \rho_w g }{ \eta }$

K_s = r_0 ^2 * f ^3 * rho_w * g /(8* q_0 *(1- f )^2* eta )


$ l = q_0 \displaystyle\frac{(1- f )}{ f } \Delta L $

l = q_0 *(1 - f )* DL / f


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

N_p = f * S /( pi * R ^2)


$ R = \sqrt{\displaystyle\frac{ f }{1- f }} r_0 $

R = sqrt( f /(1- f ))* r_0


$ R_h = \displaystyle\frac{8 \eta q_0 }{ r_0 ^2}\displaystyle\frac{(1- f )^2 }{f ^3}\displaystyle\frac{ \Delta L }{ S }$

R_h = 8* eta * q_0 * (1- f )^2* DL /( r_0 ^2* f ^3 * S )


$ S_p = N_p \pi R ^2$

S_p = N_p * pi * R ^2


$ K_s = \displaystyle\frac{ r_a ^2}{8 q_a }\displaystyle\frac{ f ^3}{(1- f )^2}\displaystyle\frac{ \rho_w g }{ \eta } \gamma $

K_s = r_a ^2* f ^3* rho_w * g * gamma /(8* q_a *(1- f )^2* eta )


$\ln \gamma = s_a g_a + s_i g_i + s_c g_c $

\ln \gamma = s_a * g_a + s_i * g_i + s_c * g_c

ID:(15222, 0)



Porosité basée sur la section

Équation

>Top, >Modèle


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 :

ID:(938, 0)



Coupe des pores

Équation

>Top, >Modèle


Étant donné que la section transversale d'un pore de rayon du tube ($R$) 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 :

ID:(4362, 0)



Nombre de pores

Équation

>Top, >Modèle


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 :

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

$S_t$
Coupe de sol
$m^2$
6041
$N_p$
Nombre de capillaires
$-$
6039
$\pi$
Pi
3.1415927
$rad$
5057
$f$
Porosité
$-$
5805
$R$
Rayon de porosité
$m$
6021

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

$ f = \displaystyle\frac{ S_p }{ S }$



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

$ S_p = N_p \pi R ^2$



donne

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



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

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

ID:(4363, 0)



Rayon des pores

Équation

>Top, >Modèle


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

$ R = \sqrt{\displaystyle\frac{ f }{1- f }} r_0 $

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 rayon d'un grain générique ($r_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$) :

$ f = \displaystyle\frac{ S_p }{ S }$



Si le nombre de grains est égal à nombre de capillaires ($N_p$) avec l'expression :

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



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

$ R = \sqrt{\displaystyle\frac{ f }{1- f }} r_0 $

ID:(109, 0)



Longueur capillaire

Équation

>Top, >Modèle


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 :

$ l = q_0 \displaystyle\frac{(1- f )}{ f } \Delta L $

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 :

ID:(2215, 0)



Flux poreux total

Équation

>Top, >Modèle


Le flux total ($J_{Vt}$) est calculé en multipliant ($$) par la valeur de le volumique flux ($J_V$) dans chaque capillaire, comme suit :

$ J_{Vt} = N_p J_V $

$J_V$
Flux Hagen-Poiseuille
$m^3/s$
6023
$J_{Vt}$
Flux total
$m^3/s$
6042
$N_p$
Nombre de capillaires
$-$
6039

ID:(4364, 0)



Écoulement à travers un sol poreux (Kozeny-Carman)

Équation

>Top, >Modèle


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

$ J_{Vt} = N_p J_V $



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

$ J_{Vt} =-\displaystyle\frac{ r_0 ^2}{8 \eta q_0 }\displaystyle\frac{ f ^3}{(1- f )^2}\displaystyle\frac{ S }{ \Delta L } \Delta p $

$S_t$
Coupe de sol
$m^2$
6041
$J_{Vt}$
Débit selon la loi de Hagen Poiseuille pour les sols
$m^3/s$
6043
$\Delta p_2$
Différence de pression 2
$Pa$
5820
$L$
Longueur capillaire
$m$
6022
$f$
Porosité
$-$
5805
$R$
Rayon de porosité
$m$
6021
$\eta_w$
Viscosité de l'eau
$Pa s$
6020

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

$ J_{Vt} = N_p J_V $



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

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



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

$ R = \sqrt{\displaystyle\frac{ f }{1- f }} r_0 $



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

$ l = q_0 \displaystyle\frac{(1- f )}{ f } \Delta L $



nous obtenons

$ J_{Vt} =-\displaystyle\frac{ r_0 ^2}{8 \eta q_0 }\displaystyle\frac{ f ^3}{(1- f )^2}\displaystyle\frac{ S }{ \Delta L } \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)

ID:(4365, 0)



Densité de flux entre les colonnes

Équation

>Top, >Modèle


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 :

$ \Delta p = \rho_w g \Delta h $



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

$ j_s =-\displaystyle\frac{ r_0 ^2 }{8 q_0 }\displaystyle\frac{ f ^3 }{(1- f )^2}\displaystyle\frac{ \rho_w g }{ \eta }\displaystyle\frac{ \Delta h }{ \Delta L }$

$J_{Vt}$
Circulation à travers l'échantillon
$m^3/s$
6044
$S_t$
Coupe de sol
$m^2$
6041
$\rho_w$
Densité du liquide
$kg/m^3$
5407
$\Delta h$
Hauteur de la colonne de liquide
$m$
5819
$L$
Longueur capillaire
$m$
6022
$f$
Porosité
$-$
5805
$R$
Rayon de porosité
$m$
6021
$\eta_w$
Viscosité de l'eau
$Pa s$
6020

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 :

$ \Delta p = \rho_w g \Delta h $



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 :

$ J_{Vt} =-\displaystyle\frac{ r_0 ^2}{8 \eta q_0 }\displaystyle\frac{ f ^3}{(1- f )^2}\displaystyle\frac{ S }{ \Delta L } \Delta p $



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 :

$ j_s =-\displaystyle\frac{ r_0 ^2 }{8 q_0 }\displaystyle\frac{ f ^3 }{(1- f )^2}\displaystyle\frac{ \rho_w g }{ \eta }\displaystyle\frac{ \Delta h }{ \Delta L }$

ID:(4366, 0)



Facteur d'échelle capillaire

Équation

>Top, >Modèle


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

$ \gamma = \displaystyle\frac{ r_0 ^2}{ r_a ^2}\displaystyle\frac{ q_a }{ q_0 }$

.

ID:(15101, 0)



Conductivité hydraulique et facteur d'échelle

Équation

>Top, >Modèle


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

$ K_s = \displaystyle\frac{ r_a ^2}{8 q_a }\displaystyle\frac{ f ^3}{(1- f )^2}\displaystyle\frac{ \rho_w g }{ \eta } \gamma $

.

ID:(977, 0)



Régression pour la conductivité hydraulique

Équation

>Top, >Modèle


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 :

$\ln \gamma = s_a g_a + s_i g_i + s_c g_c $



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:(15100, 0)