Calculation of Capillary Rise Height of Soils by SWCC Model

*e maximum capillary rise height of soil is a complex system which is mainly determined by the distribution characteristics of soil pores.*e tests of the rising height of capillary water on 8 kinds of soils by the method of vertical tube are widely conducted to measure the maximum capillary rise height. Based on the BCC model and principles of thermodynamics, the soil-water characteristic curve test is designed for the purpose of calculating the pore distribution of soil samples. A new method for calculating the maximum capillary rise height of soil is proposed by the author by using the distribution function of the soil pore. *e coefficient βwhich reflects the relationship between the maximum capillary rise and the average pore radius of soils is utilized during the calculation process, and then the reference range of β for different soils is obtained according to series of experiments corresponding.*e proposed calculation method offers an effective way to calculate the maximum capillary rise height, which can be applied to analyze the capillary effect area of relevant engineering problems.


Introduction
Capillary-driven liquid flow is the main transport mechanism in the soil system of which the water erodes continuously by capillary rise from a lower elevation to higher elevation.Such capillary rise phenomena lead to an increase of the saturation of the soil, which will not only decrease the strength of the soil but also alter the elastic modulus of substructure soils, thereby leading to the corresponding changes in stress and strain response under the external load, for example, the traffic load.
erefore, the defect of the roadbed is closely related to capillary rise erosion.Seasonally frozen ground has always been an important problem in highway construction and channel slope at high latitudes.In the research of frost heaving zone, it is inevitable to determine the height of capillary rise of the substructure.For channel slope and embankment close to the riverside, capillary action path is shortened as well as the effect of capillary action on the supply of water is accelerated.As a result, these infrastructures are more prone to defect and frost heave in such particular areas because of capillary rise.To sum up, study of capillary rise, particularly the maximum capillary rise height, is of great significance to the design of the substructure and channel because the maximum capillary rise height is tightly connected to the strength reduction region and frost area.
A series of studies have been carried out on the capillary rise; LU offered a complete analytical solution for the relationship between the rate and time of capillary rise in soils [1].However, in the equation, the maximum height of capillary rise h c is a known soil parameter.In other words, the accurate determination of the maximum height of capillary rise is the prerequisite to calculate the rate of capillary rise.e maximum height of capillary rise has an important influence on the overall engineering behavior of unsaturated soils and is a highly complex system of both the soil and pore water properties.It is extremely difficult to calculate the maximum capillary height accurately in real soils.In order to overcome this, scholars have done a lot of studies and have put forward some empirical formulas which established the correlation between the maximum capillary height and certain measured soil parameters.e earliest formula was proposed by Lane and Washburn [2,3], after conducting the capillary rising test for 8 kinds of different soils, and the result shows that the maximum capillary height is linear with ln D 10 : h c � −990 ln D 10  − 1540. (1) On the basis of (1), Peak and Hansen put forward another empirical formula [4]: In this formula, 10% particle size d 10 , void ratio e, and the coefficient C(C � 3 × 10 −5 m 2 ∼ 8 × 10 −5 m 2 ) are used to calculate the maximum capillary rise height.However, the coefficient C is difficult to choose because the range of C is relatively large, and the value of C depends hugely on the real condition of soil.Kumar and Malik by carrying out indoor tests summarized as follows [5]: where h a is the height corresponding to air-entry pressure value of soils and r is the equivalent of capillary radius of soils.In this formula, the unit of r is μm, so it has little influence on the result of the calculation.erefore, h a is the only key computational parameter which can be determined by the soil-water characteristic curve.However, it requires an instrument to control the matric suction very precisely especially when the suction is less than 1 kPa [6].
We summarized that according to previous research results, using the micropore distribution of soils and using the pore radius parameter instead of the gradation parameter and void ratio are effective methods to accurately calculate the maximum capillary rise height.In this paper, a new method is proposed to calculate the maximum capillary rising height by using the radius of microscopic pore distribution.A large number of capillary rise tests as well as soil-water characteristic (SWCC) tests have been done for different kinds of soils, from which the parameters of different soils required for the calculation are obtained, and the feasibility of this method will be verified.

Capillary Rise Equilibrium Equation in Soil System.
In a single capillary when water column reaches the maximum height, the gravity of the water column is balanced by the surface tension along water-solid interface, as shown in Figure 1. e balance equation is as follows: where α is the liquid-solid contact angle; ρ w is the density of water; r is the radius of the capillary channel (r � d/2); and T S is the surface tension of water.
For soils, we regard it as homogenous porous media, which means the pores inside are similar in diameter and the dispersion is small.In the interior of soils, tiny pores which are connected with each other fulfill the precondition of capillary rise [7,8].e total capillary gas-liquid area A in each cross section can be expressed as follows: where ϕ A is the number of gas-liquid interfaces in total capillary gas-liquid area A of each cross section and f(r) is pore size distribution (PSD) function in the porous system of soils [9,10].When we established the gravity-tension balance equations assuming that soil pores are interconnected with each other, the capillary migration path of water is not like a single capillary tube, meaning that channel interleaving was bent and most of the gas-liquid interfaces on the same section are connected.As a simplified model, we argue that the total capillary route is equal to the maximum capillary rise height h c multiplied by a path retardation coefficient β; thereafter, the capillary equilibrium equation for any gas-liquid interface can be modified from (4) to When ρ w � 1 g/cm 3 , g � 980 cm/m 2 , T S � 72 mN/m, and α � 0 °(dehydration process), (7) can be simplified as follows:

Testing Method of PSD f(r).
One effective way to study the pore characteristics of homogenous porous media is to establish the microscopic pore structure model.Many conventional models for pore characteristics, liquid distribution, flow, and transport in porous media are based on representing pore space geometry as a bundle of cylindrical capillaries (BCC model) [11,12].In the literature, there have been some standardized experimental methods for directly measuring pore size distribution (MIP, mercury intrusion porosimetry; BM, bubble method), and these methods are based on the basic theory of the BCC model.In BCC model, according to the Young-Laplace equation, the relationship between the pore radius r and the suction S (or capillary pressure) is represented as follows: Based on the assumption of pore local equilibrium, the balance of soil is achieved in a given suction S * inside the pore while one part of the pore with the pore radius greater than r * (r * � 2T s cos α/S * ) is filled with air and the other one with pore radius less than r * is filled with water.By defining g(r) as pore size distribution function, g(r) dr represents the percentage of void volume with pore radius ranging from r to r + dr under the unit soil volume.When the suction of the soil reaches S, the maximum pore radius which is filled with water becomes C/S(C � 2T S cos α).erefore, the relationship between soil volumetric water content and pore distribution can be expressed as follows: 2 Advances in Civil Engineering dθ � g(r) dr. ( In this formula, θ indicates the volumetric water content of soil when soil suction equals S and dθ represents the ratio of pore water volume with soil particle radius between r and r + dr to the total pore water volume under a well-saturated condition.In fact, in this formula, the volume of water being evacuated is used as an alternative to that of void in the soil.When the calculated radius is in the range of the whole pore medium, g(r) can be equal to f(r).
In this paper, the pore size distribution of soil was measured, and the methods mentioned above have some limitations when being applied.Taking MIP as an example, sample preparation for MIP is difficult, and under high pressure, liquid mercury will jeopardize soil particles, resulting in the deformation of soil skeleton and pore during the test.Another important reason which restricts the widespread use of conventional methods is the low popularity of test instruments as well as high test costs, which forced us to look for other convenient testing methods.In this case, mercury is replaced by water, and the soil-water characteristic curve is utilized instead of the mercury injection curve.Pores among soil particles rather than the internal pore of the particle are mainly considered with regard to the soil capillary phenomenon and the pore water flow problem.When the soil-water characteristic curve is used, the pressure can be increased to 2 MPa, which means the pore radius can be attached to 1000 Å correspondingly.
In other words, it can obtain the exact pore distribution of the above 1000 Å, allowing us to have enough pore radius distribution data to analyze soil capillary phenomenon.

Basic eories and Calculation Method.
e pore of soil can be regarded as a capillary in the Kelvin formula without violating original characteristics of soil system.Radius, contact angle between soil particles and pore water, and surface tension are denoted by r, α, and T S , respectively.When using measurable physical quantities, Kelvin formula is expressed as where form of matric suction is described as where R � 8.314 J/(mol•K) is the universal gas constant; T � 298 K is the thermodynamic temperature; μ vo (kPa) is the free water saturated vapor pressure in equilibrium state under the condition of temperature T. μ v1 (kPa) is the vapor pressure value of solution at the state of balance; v m � 18 × 10 −6 m 3 /mol is the partial molar volume of water vapor.u a − u w is the matric suction, u v1 /u v0 is defined as relative humidity (RH) and can be obtained from ( 13) which is another expression of ( 11): When the relative humidity or matric suction is in the increment of Step i, the volume change of the pore volume of soil mass which is replaced by the gas can be calculated as Kelvin radius can be estimated by (12), while the actual pore radius is expressed as where the thickness of adsorbed water film t i corresponding to RH or matric suction at Step i can be expressed by Halsey equation as [13] where τ is the effective diameter of adsorbed water molecules.Assuming that the cross-sectional area occupied by the water molecules is approximately 10.8 Å and the Avogadro constant N A is 6.02 × 10 23 /mol, the value of τ can be calculated as

Advances in Civil Engineering
After terminating the calculation of pore radius which corresponds to different suctions, points can be selected in the SWCC to calculate the pore size distribution curve.Considering the nonuniform pore size and inhomogeneous pore distribution as well as the data of pore distribution curve obtained by point-selecting method, r 0 as average pore radius of soil by using weighted average of the statistical theory is defined as Similar to the parameters of d 10 , d 60 , and C u found in the particle size distribution curve, the average capillary pore radius can be obtained by pore size distribution curve.According to the definition, r 0 represents the average radius within the calculating range.For example, if suction value ranges from 1 kPa to 10000 kPa, then the r 0 is the average radius of the pore which falls into the corresponding range.
When the calculated aperture range relates to the whole soil (suction value ranges from the air-entry pressure value S S to the residual suction S r ), we can approximate r 0 as r 0 �  n i�1 r i /n, by  n i�1 r i ≈ nE(r) where We can easily argue that By conducting substitution of (20) to ( 8), ( 8) can be simplified as follows: As a parameter, average pore radius r 0 characterizing the pore state can be obtained by SWCC, and path retardation coefficient β as an empirical parameter will be measured by the following tests.

Test Materials.
e test soils are distinguished as two different kinds of soils: fine-grained soils and coarse-grained soils.Fine-grained soils include low liquid limit clay (CL), high liquid limit clay (CH), low liquid limit silt (ML), and high liquid limit expansive clay (CHE), while coarse-grained soils include silt sand (SS), fine sand (FS), medium sand (MS), and coarse sand (CS).
e parameters of the soil samples are obtained by the conventional tests and are shown in Tables 1 and 2.

Capillary Rise Test.
e method of vertical tube (a transparent plexiglass tube with diameter of 4.5 cm and height of 500 cm) is applied when performing soil capillary rise test.20 capillary plexiglass tubes with 8 kinds of soil samples are tested and the compaction degree of soils is controlled by the total density of soils during the sample loading process.e test results of the steady rising height of capillary water in different soil samples are shown in Table 3.

SWCC Test.
SWCC is mainly determined by pressure plate method, in which four samples of each kinds of soil and the as-fabricated samples were placed in a vacuum saturator for 48 h to ensure the full saturation.PF moisture characteristic curve instrument (DIK-4303) is used for measuring soil-water characteristic curve.First of all, samples are weighted by using electronic balance to ensure that the mass is recorded.And then, samples are put into a determinator (4 samples at a time) which connects a water channel between the instrument and the ceramic plate.
e upper cover should be well tightened in case of air leakage.e maximum pressure is set to 1200 kPa, which is controlled by the increment of each step of the 100 kPa, and before applying the pressure, each sample's gravimetric mass is recorded.e suction balance will be achieved if no water flows out from the pipe within 48 hours.

Analysis and Calculation
e results of the SWCC test are shown in Figures 2 and 3, but due to the design restriction in pressure plate method which means that a pressure over 1200 kPa cannot be exerted and maintained in the experiment, points exceeding 1200 kPa in SWCC are obtained by fitting with VG model [14].To be precise, the datum acquired from the test is injected into Van-Genuchten model in ORIGIN to obtain soil-water characteristic curve in a wider range of suction.
e fitting equation is where θ is the volumetric water content; θ s is the saturated volumetric water content; θ r is the residual volumetric water content; and a, m, and n are independent coefficients (satisfy the condition m � 1 − 1/n as much as possible).e fitting results from (22) for the soil samples are shown in Table 4.When selecting suction calculation range, the value range should be considered to represent the entity of soil as well as the measurement range of suction.After many attempts, we propose that fine-grained soils matric suction range be controlled at 200 kPa∼10000 kPa and that the range of coarse-grained soils be controlled at 50 kPa∼4000 kPa.( e external values of the test range are calculated by fitting the results of VG model.) According to the theory introduced above, soil-water characteristic curve is used to calculate the radius r 0 .Calculation results of soil sample of CL with compaction degree 0.9 are listed in Table 5 as an example while the numerical integration calculation steps are as follows: (1) Get the SWCC of the soil and select the data of 200 kPa∼20000 kPa (the data outside the range of test are obtained by fitting with VG model).(2) Calculate RH (relative humidity) of matrix suction according to (13).(3) Translate mass water content into per unit mass of soil-water filled pore volume (volumetric water content).( 4) Use (12) to calculate the Kelvin radius r i k under the suction of different substrates with which the contact angle α is 0 (soil in dehydration process).

Advances in Civil Engineering
(5) Calculate the thickness of the water film t i under the suction of different substrates by ( 16). ( 6) Calculate actual pore radius r i p by results ( 4) and ( 5) applying (15).( 7) Calculate the reduction of pore volume water Δv p in per unit mass of soil when the matrix suction is changed form one substrate to another.(Here, Δv p is the pore volume corresponding to the water that is discharged.)(8) Calculate the average pore radius (r p )avg in the process of evacuation of pore water by averaging the result of ( 6) for per step.(Here (r p )avg is the average pore radius corresponding to the pore volume v p .) (9) Add the pore volume to calculate the cumulative pore volume  (V p ) of per unit mass, and draw relationship diagram between the pore volume Δv p and the average pore radius (r p )avg.(10) Calculate the r 0 by (18).e average pore radius r 0 of soil sample is calculated by the method above.After determining the radius r 0 , the retardation coefficient β of different soil samples will be obtained.And the results of coefficient β are shown in Table 6.
As clearly observed from the above chart, the coefficient of the fine-grained soils was about 21 cm −2 while the range barely fluctuates with value relatively fixed.But for coarsegrained soils, β changes from 23 cm −2 to 26 cm −2 .One main reason we argue for coarse-grained soils is that the proportion of large pore size increased significantly, contributing greatly to a maximum capillary rise.As a result, for   Advances in Civil Engineering 7 (21) the numerical value of the coefficient β is considered to be 21 cm −2 and 25 cm −2 for fine-grained soils and coarsegrained soils, respectively.

Verification and Discussion
In this study, in order to verify the accuracy of the maximum capillary rise height formula proposed by the author, capillary rise tests are carried out on 44 kinds of soil samples, and thereafter, a comparison is made between the capillary rise estimated from our equation (21) and the estimations from (1)-(3).Here, the coefficient C for (2) takes 5 × 10 −5 m 2 e test results and calculation results of the maximum capillary rise height for 44 kinds of soil samples are shown in Table 7.
It can be seen from the table that the estimating error of the maximum height from the formula proposed by the author is always less than 10% and the maximum error is 22 cm (204 cm), indicating its relatively higher accuracy compared with other formulas mentioned in this paper.is  Advances in Civil Engineering also means that the results coincide better with the actual situation and have higher practical value.Accurately determining the maximum capillary rise height of the soils is a complicated issue.In this paper, we calculated the average pore radius by method of SWCC and then determined the maximum capillary rise height using the pore radius.Obtaining the pore distribution data from the SWCC data is the core of the calculation method.We can easily get the pore distribution curve through the soil-water characteristic curve.However, similar to the particle size curve, we need to find a set of parameters to evaluate this curve as well as calculate the maximum capillary rise height.At first, we tried to use r 50 as the evaluation parameter but later found when fine-soil samples are at different degrees of compaction, the r 50 values are very close.According to the definition, as long as the soils have the same matrix suction at the saturation of 50%, the values of r 50 still remain the same.It is difficult to use a specific pore radius to represent the equivalent capillary radius in the calculation of the maximum capillary rise height of the soil; therefore, we chose the weighted average of the pore radius r 0 as the calculation parameter.Whether there is a better way to evaluate the distribution of soil pore and to analyze the behavior of soils is a project that we still need to explore.
Another important parameter is β which is an empirical one obtained from the test.We use r 0 to distinguish different soil samples, so we hope different soil samples would obtain the same values of β as much as possible.Test results show that β of fine-grained soil is very close and can be directly taken to 21, but β also has the tendency to augment with the increase of average pore radius r 0 and this trend is particularly evident in coarse-grained soils.erefore, there is an issue about how to pick the value of β in the actual calculation process and we argue that β of fine-grained soils can be fixed to 21, but for coarse-grained soil the value of β is not a settled figure, ranging from 24 to 26.Moreover, the value rises with the increase of the content of the coarse particles.Here the author suggests taking an average of 25.

Conclusion
After theoretical calculation and large numbers of indoor experiments, the following conclusions can be summarized as follows: (1) e average pore radius formula is deduced on summary of experiments and statistical theory, and the calculation theory can be used to distinguish the pore distribution characteristics of different compaction degree of soil samples with different particle size distribution.erefore, it can offer references for similar research and also can be used for estimating the maximum height of capillary.In particular, for fine-grained soils the calculation results for the maximum capillary height are relatively accurate.(2) Using the soil pore radius to calculate the maximum capillary rise height is an accurate and feasible approach.During the calculation process, the coefficient β of the fine-grained soil was about 21 to 22, while β of the coarse-grained soil was about 24 to 26. e greater the number of coarse particles is, the greater the retardation coefficient β is.
(3) e maximum capillary rise height of the soil is mainly determined by the distribution of the large size pore, and meanwhile, the effect of the small size pore is quite insignificant.Fine-grain diameter of ten percent of total soil mass e:

Notations
Void ratio of soils h a : Height corresponding to air-entry pressure value of soils α: Liquid-solid contact angle r: Radius of capillary channel ρ w : Density of water T S : Surface tension of water A: Total capillary gas-liquid area in each cross section ϕ A : Number of gas-liquid interface in total capillary gas-liquid area A f(r): Pore size distribution (PSD) function of soil β: Path retardation coefficient S: Matric suction R: Universal gas constant T: ermodynamic temperature μ v0 : Free water saturated vapor pressure in equilibrium state under the condition of temperature T μ v1 : Vapor pressure value of solution at the state of balance v m : Partial molar volume of water vapor u a : Pore air pressure u w : Pore water pressure RH: Relative humidity v p : Pore volume of soil r k : Kelvin pore radius r p : Actual pore radius t: ickness of adsorbed water film τ: Effective diameter of adsorbed water molecules N A : Avogadro constant 6.02 × 10 23 /mol d 60 : Fine-grain diameter of 60 percent of total soil mass C u : Nonuniform coefficient of soil r 0 : Average pore radius θ: Volumetric water content θ S : Saturated volumetric water content θ r : Residual volumetric water content a, m, n: Independent coefficients of VG model.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Figure 1 :
Figure 1: Mechanical equilibrium of capillary rise in the small diameter tube.
h c : Maximum capillary rise height of soils D 10 : Particle size d 10 of soils d 10 :

Table 1 :
Main physical indexes of fine-grained soils.

Table 2 :
Main physical indexes of coarse-grained soils.

Table 4 :
Fitting results of SWCC parameters for soil samples.

Table 5 :
Equivalent radius calculated form of soil sample of CL with compaction degree 0.9.

Table 6 :
Calculation result of coefficient β for soil samples.

Table 7 :
Test and calculation results of maximum capillary rise height for 44 kinds of soils.