Guida Online

Material Models in Flow Analysis

The steady state flow analysis is driven by Darcy’s law specifying the relationship between the flow velocity and the gradient of hydraulic head. The current version of the program assumes constant values of coefficients of permeability independent of pore pressure.

The program also requires specifying the initial void ratio e₀ for the determination of the current porosity n and subsequently the actual velocity of water flowing through pores only v_s = v/n, where v is the average flow velocity through the entire seepage area. Generally, the void ratio e₀ = 1 corresponds to soil porosity of n = 50%.

By introducing the coefficient of relative permeability K_r the program allows for tracking the transition zone between fully saturated (S = 1, K_r = 1) and unsaturated (K_r => 0) region of the soil body. As an example we may consider the problem of unconfined flow. The process of tracking the transition zone is governed by one of the three models of transition zone determining the evolution of coefficient of relative permeability K_r as a function of pore pressure head, see Figures.

(a) Log-linear model [1],  (b) Van Genuchten model [2]

Log-linear model

The Log-linear transition zone model described e.g. in [1] is defined by the following parameters:

The coefficient of relative permeability K_r is given by:

Gardner model

This is an equivalent model depending on a single parameter α [1/m] only. The coefficient of relative permeability K_r is in this case given by [4]:

Van Genuchten model

In this case, the value of coefficient of relative permeability K_r is given by:

where α [1/m], n > 1, m = 1 - 1/n are model parameters. Their values can be obtained from laboratory measurements of retention curves approximated by:

Parameter Θ is in general provided by:

The current degree of saturation S can be expresses as a ratio of the water content θ and porosity n as follows (it is necessary to distinguish between n representing porosity and n, which appears in the van Genuchten model, they are two different variables):

The Log-linear and Gardner models adopt a simplified version of the van Genuchten model according to [5]:

where b > 0 [-] is a fitting parameter allowing for a better approximation of the retention data of a given soil.

We recommend the following tables of parameters which are actually used in the program. These parameters are derived from those given in the original tables.

Optimal values of parameters of the van Genuchten model for various classifications based on USDA and FAO are presented in the following tables.

Table with regression coefficients for grain size USDA according to Van Genuchten (1991)

Table with regression coefficients for grain size FAO according to Van Genuchten (1998)

Table: FAO texture classification system

Literature:

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.