Aide contextuelle

Modèles de matériau dans l'analyse des écoulements

L’analyse des écoulements en régime permanent est régie par la loi de Darcy qui spécifie la relation entre la vitesse d’écoulement et le gradient de la charge hydraulique. La version actuelle du programme suppose des valeurs constantes des coefficients de perméabilité indépendants de la pression interstitielle.

Le programme impose également de valoriser l'indice des vides initial e₀ permettant de déterminer la porosité courante n, puis la vitesse réelle de l’eau traversant uniquement les pores v_s = v/n, où v est la vitesse moyenne d’écoulement dans toute la zone d’infiltration. En règle générale, à l'indice des vides vérifiant e₀ = 1 correspond une porosité du sol valant n = 50%.

En introduisant le coefficient de perméabilité relative de K_r, le programme permet de suivre la zone de transition entre la région complètement saturée (S = 1, K_r = 1) et la région non saturée (K_r → 0) du corps de sol. A titre d'exemple, nous pouvons considérer le problème du flux non confiné. Le processus de suivi de la zone de transition est régi par l’un des trois modèles de zone de transition déterminant l’évolution du coefficient de perméabilité relative de K_r en fonction de la hauteur de pression interstitielle, voir Figures.

(a) modèle Log-linéaire [1],  (b) modèle de Van Genuchten [2]

Modèle Log-linéaire

Le modèle log-linéaire de zone de transition, décrit par exemple dans [1], est défini par les paramètres suivants :

Le coefficient de perméabilité relative K_r est donné par :

Modèle de Gardner

Ceci est un modèle équivalent dépendant uniquement d'un seul paramètre α [1/m]. Le coefficient de perméabilité relative K_r est dans ce cas donné par [4] :

Modèle de Van Genuchten

Dans ce cas, la valeur du coefficient de perméabilité relative K_r est donnée par :

où α [1/m], n > 1 et m = 1 - 1/n sont les paramètres du modèle. Leurs valeurs peuvent être obtenues à partir des courbes de rétention approchées résultant de mesures en laboratoire par :

Le paramètre Θ est en général fourni par :

Le degré actuel de saturation S peut être exprimé comme le rapport de la teneur en eau θ et de la porosité n comme suit (il est nécessaire de distinguer entre n représentant la porosité et n qui apparaît dans le modèle de van Genuchten, ce sont deux variables différentes) :

Les modèles Log-linéaire et de Gardner adoptent une version simplifiée du modèle van Genuchten selon [5] :

où b > 0 [-] est un paramètre d'ajustement permettant une meilleure approximation des données de rétention d'un sol donné.

Nous recommandons les tableaux de paramètres suivants qui sont réellement utilisés dans le programme. Ces paramètres sont dérivés de ceux donnés dans les tableaux d'origine.

Les valeurs optimales des paramètres du modèle de van Genuchten pour diverses classifications basées sur l'USDA et la FAO sont présentées dans les tableaux suivants.

Tableau des coefficients de régression de Van Genuchten pour les classes de taille de grains de l'USDA (1991)

Tableau des coefficients de régression de Van Genuchten pour les classes de taille de grains de la FAO (1998)

Tableau : Système de classification des textures de la FAO

Littérature :

Details can be found in [2].

[1] D.M. Potts, L. Zdravkovič, Finite element analysis in geotechnical engineering - theory, Thomas Telford, London, 1999.

[2] M. Th. Van Genuchten, A closed formulation for predicting the hydraulic conductivity of unsaturated soils, Journal Soil Science Society of America 44, 239-259, 1988..

[3] M. Šejnoha, Finite element analysis in geotechnical design, to appear (2015)

[4] W. R. Gardner, Some steady-state solutions of the unsaturated moisture flow equation to evaporation from a water table, Soil Science 85(4), 228-232, 1958.

[5] M. Šejnoha, T. Janda, H. Pruška, M. Brouček, Metoda konečných prvků v geomechanice: Teoretické základy a inženýrské aplikace, předpokládaný rok vydání (2015)

[6] USDA 1951. Soil Survey Manual. Soil Conservation Service. U.S. Department of Agriculture Handbook No. 18. US Government Printing Office. Washington DC.

[7] Wösten, J.H.M., et. al. 1998. Using existing soil data to derive hydraulic parameters for simulation models in environmental studies and in land use planning. Final Report on the European Union Funded project. DLO Winand Staring Centre. Report 156, Wageningen, NL. p. 106. ISSN 0927-04537.