An erratum for this article has been published. To view the erratum, please click here.

Advances in Civil Engineering

Volume 2016, Article ID 8176728, 7 pages

http://dx.doi.org/10.1155/2016/8176728

## Modeling of Hydrophysical Properties of the Soil as Capillary-Porous Media and Improvement of Mualem-Van Genuchten Method as a Part of Foundation Arrangement Research

^{1}Peter the Great St. Petersburg Polytechnic University, Polytechnicheskaya 29, Saint Petersburg 195251, Russia^{2}Leibniz-Centre for Agricultural Landscape Research, Eberswalder Straße 84, 15374 Müncheberg, Germany

Received 30 November 2015; Revised 22 February 2016; Accepted 9 March 2016

Academic Editor: Hossein Moayedi

Copyright © 2016 Vitaly Terleev et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Within the concepts about the capillarity and the lognormal distribution of effective pore radii, a theoretical justification for function of differential water capacity and its antiderivative (function of water-retention capacity in form of a dependence of the soil volumetric water content on capillary pressure of the soil moisture) is presented. Using these functions, the ratio of soil hydraulic conductivity function to the filter coefficient is calculated. Approximations to functions describing the water-retention capacity and relative hydraulic conductivity of the soil have been suggested. Parameters of these functions have been interpreted and estimated with applying the physical and statistical indices of the soil.

#### 1. Introduction

In laboratories usually do not measure the differential water capacity (DWC) function but its antiderivative-water sorption equilibrium isotherm as a dependence of the volumetric soil moisture (cm^{3}·cm^{−3}) on capillary pressure of the soil moisture (cm H_{2}O). The dependence is an indicator of water-retention capacity (WRC) of soil. Calculating the values of DWC is usually limited to interpolating the measured WRC data and to selection of the approximating function, which is then used in the computation by the formula . However, differentiation of approximations is by definition an unstable operation, which can lead to physically absurd results. As a negative example, the degree function was previously widely used for the interpolation of modified WRC points for the subsequent calculation of the DWC values. However, there are some positive examples of the use of WRC approximations, which include the model:where is the relative volumetric water content (effective saturation); is the volumetric water content under conditions of moisture-saturated soil (cm^{3}·cm^{−3}); is the volumetric water content, which corresponds to the minimum of specific volume of water as liquid in soil (cm^{3}·cm^{−3}); (cm ), , and are parameters.

When , model (1) is the WRC proposed by Van Genuchten [1]. When , relation (1) has the form of the precursor’s WRC [2–4]. Van Genuchten uses model (1) to calculate the ratio of the soil hydraulic conductivity function (HCF) to filter coefficients of soil moisture by Mualem method [5]:

Using the variable , (2) reduces to

In the case , (3) permits a simple analytical calculation, and the result of integration is as follows:

Description of hydrophysical properties of the soil using WRC (1) and HCF (4), having common parameters, is called Mualem-Van Genuchten method.

The purpose of this paper is as follows:(i)theoretical justification and formation of a system of functions with common and adequately interpreted parameters, which combine DWС, its primitive in WRC form, and the ratio of HCF to the filter coefficient of soil moisture;(ii)approximation of and in the class of elementary functions and parameter estimation of approximating functions by the physical parameters of the soil.

#### 2. Materials and Methods

##### 2.1. Differential Water Capacity and Water-Retention Capacity of the Soil

It is known that pores in natural structured soils are preferably capillary. In cross section capillaries are very different by the configuration and area, which is caused by the random combination of contacting soil particles of varying shape and size. As a basis for hydrophysical properties of the soil modeling view of the system of cylindrical pores of circular cross section is adopted. This system is equivalent to the actual pore space of the soil by its capillary properties. The model of random lognormal value—the effective radius of the pores—is used to describe the distribution of the volumes of soil capillary [6–9]. By analogy with Kosugi work, is accepted as an effective radius of the pores of the soil, where is radius of the pore; is radius of the smallest pore; is radius of the widest pore [10]. The relations for calculating the proportion of pore volume are recorded, taking into account the random nature of the cross section of soil pores. This proportion accounted for the capillaries, starting with the smallest and ending with effective pore radius :where is lognormal distribution density of the random variable ; and are the most probable value and standard deviation of the logarithms of the effective soil pores radii, respectively. The variable is introduced and formula (5) is rewritten in the waywhere is error function ( and , when ); .

From formulas (6) the relation should be:

The approximation by Winitzki [11] was used , . We reduced it to a simpler form: . The approximate equality follows from formulas (6) and (7): ; it obtains an analogue of the Verhulst equation. Decision of the resulting nonlinear ordinary differential equation of the first order with the boundary condition is logistic curve, which is approximately the ratio in a class of elementary functions:

Given the connection between the values of and , reduce formula (8) in the form

The possibility of moving from to and from to is considered for conversion of formula (9). The difference between the absolute pressures under the curved boundary between “air-capillary moisture” and a flat free water surface is called a capillary moisture pressure, which is calculated by the law of Laplace: , where is soil capillary radius; ; is the surface tension coefficient of soil moisture at the interface with the air; is contact angle of the surface of soil particles with water; is acceleration of gravity; is density of water. Under the physical modelling of soil drainage using a pneumatic press, water displacement is achieved by exposure of the excess pressure of atmospheric gases. Dissolution of the air in water occurs with the increase of gas pressure above the water in the pores of the soil, from the normal external atmospheric pressure to the value . Excess gas pressure is transferred to the water; the absolute pressure is also increased. When excess gas pressure reaches a value , strength of the interaction between water molecules and the surface of soil particles weakens so that the water goes into the category of free gravitational water which follows from the soil through the membrane into a tray; the air occupies the vacant pore volume. Excess gas pressure in the pneumatic press with respect to the external atmospheric pressure corresponds to the value of the capillary pressure of moisture in the soil of the largest pore . This value is called “bubbling pressure”; it will be considered as the “initial point” of the capillary pressure. With use of , for case (accept ) the formula is derived:where ; is the radius of the soil pore, which corresponds to the most probable value of of the random variable .

Applying formulas (10), enter the following values:

Using law of Laplace and formula (11), relationship (9) can be represented as

Moistening of the soil begins with water filling the smallest pores, and only after they are saturated with moisture, does water start to occupy the larger pores. Drainage of the soil begins with air displacing the water from the largest pores, and only after they are emptied, does air start to occupy the smallest pores. The capillary pressure is the same at all points of topologically closed water space (law of Pascale). The biggest radius among all the moisture-filled pores determines the pressure. These factors make it possible to identify the proportion of the pores volume of with an effective radius not more than with the value of the relative volumetric water content (effective saturation) at the capillary pressure , which uniquely corresponds to . So, mathematical identity is physically reasonable. Relation is a WRC by the definition, so it allows transferring relation (12) to the following form:

Relation is an antiderivative with respect to the relative DWC function :

In case take , and from formulas (13a) and (14) we arrive at the relations, previously obtained by Kosugi [10].

##### 2.2. The Ratio of Hydraulic Conductivity Function to Filter Coefficients of Soil Moisture

Here are not any empirical relationships for application in the calculation by Mualem method, which interpolates the measured WRC data, by analogy with Kosugi [10]. Instead of this, the relative DWC function is used. Moreover, this function is theoretically justified in the frames of some concepts about capillarity and geometry of pore space of the soil. According to formulas (13a) and (13b), the value under the capillary pressure reaches 1. So, using the identity , formula (2) transforms into

Using the transition from to , the final formula, which describes the ratio of HCF to filter coefficient of soil moisture, as well as an approximation of this ratio, has been obtained:

In the particular case (), the ratio (16a) leads to the Kosugi formula [10]. Models (13a), (14), and (16a) form a parameter-closed system of soil and hydrophysical functions. Relation of (13b) and (16b) approximates WRC and the values, relatively, in a class of elementary functions.

##### 2.3. Interpretation of Parameters of the Soil Hydrophysical Functions

Theoretically justified models (13a), (14), and (16a) have the common parameters, for which physical and statistical interpretation is proposed. Parameter is determined: (1) by multiplication of the widest pore radius and the effective pore radius , in which the random variable reaches the most probable value; (2) by coefficient in formula of Laplace, which describes the capillary properties of the soil. If , than . Parameter is a value, inversely proportional to —the standard deviation of the logarithms of the effective pore radii. Parameter is the bubbling pressure on the drainage isotherm of the soil, initially saturated with moisture. Coupling of the parameters and with the certain DWC value can be revealed from the relation , where is the value of function calculated using formula (14) with , where is capillary pressure that corresponds to the effective pore radius (at the effective saturation ). The DWC function reaches its maximum at the inflection point of WRC curve. Authors consider that it corresponds to the maximum capillary-sorption moisture capacity of the soil [12, 13].

#### 3. Results and Discussion

The measured WRC data are usually used to design the curves with respect to formulas (1) and (4). Here we made simulation of dependence applying the data, for which calculating the statistics of lognormal distributed effective soil pore radii and was given, as well as some characteristics of the capillary properties of the soil pore space cm; = −100 cm Н_{2}O and cm^{2} [2, 8, 14, 15] were given. Using formulas (11) and given indices, the parameters of DWC function cm and were calculated. Then by formula (13a) with usage of , , and parameters the curve (white circles on Figure 1), which simulates data, was depicted. With that, according to formula (16a), the curve (black circles on Figure 1), which simulates the ratio, was depicted. Formulas (13b) and (16b) have been used to depict curves (1) and (2) (see Figure 1), which represent the approximation to WRC and relative HCF.