|
|
|
| EFFECT OF CAPILLARY HYSTERESIS ON TWO-PHASE FLOW IN POROUS MEDIA |
(1. State Key Laboratory of Geomechanics and Geotechnical Engineering，Institute of Rock and Soil Mechanics，
Chinese Academy of Sciences，Wuhan，Hubei 430071，China；2. Central Southern China Electric Power Design Institute，
Wuhan，Hubei 430071，China) |
|
|
|
|
Abstract The seepage process in porous media is significantly influenced by water content and drying/wetting history that the material has experienced. Based on the continuum theory of porous media and an internal state variable-based model of capillary hysteresis，a two-phase flow model of porous media is created and implemented into the U-DYSAC2 finite element code. Comparing with the experimental data，the results show that the proposed model can be reliably and efficiently used to simulate the unsaturated seepage problems of porous media under complex conditions. The new procedure is used to analyze the seepage processes in a soil slope subjected to drying/wetting cycles. The results show that the effect of capillary hysteresis on unsaturated seepage process is significant；and the hydraulic state of soil is not only related to the current value of water content or matric suction，but also depends on the drying/wetting history that the soil has experienced. In particular，it is shown that if the main drying soil-water characteristic curve is adopted in the analysis，the matric suction of soil will be over-predicted；hence the shear strength of soil and the factor of safety of slope stability can be over-estimated based on the traditional methods of stability analysis. Therefore，it is necessary to take into account the effect of capillary hysteresis in modelling unsaturated seepage problems.
|
|
Received: 02 June 2010
|
|
|
|
| [1] IVERSON R M. Landslide triggering by rain infiltration[J]. Water Resources Research，2000，36(7)：1 897–1 910.
[2] BUJANG B K H，FAISAL H J A，LOW T H. Water infiltration characteristics of unsaturated soil slope and its effect on suction and stability[J]. Geotechnical and Geological Engineering，2006，24(5)，1 293–1 306.
[3] TALEBI A，UIJLENHOET R，TROCH P A. A low-dimensional physically based model of hydrologic control of shallow landsliding on complex hillslopes[J]. Earth Surface Processes and Landforms，2008，33(13)：1 964–1 976.
[4] MUALEM Y. Modified approach to capillary hysteresis based on a similarity hypothesis[J]. Water Resources Research，1973，9(5)：1 324–1 331.
[5] POULOVASSILIS A. Hysteresis of pore water，an application of the concept of independent domains[J]. Soil Science，1962，93(6)：405–412.
[6] MUALEM Y. Conceptual model of hysteresis[J]. Water Resources Research，1974，10(3)：514–520.
[7] HAVERKAMP R，REGGIANI P，ROSS P J，et al. Soil water hysteresis prediction model based on theory and geometric scaling[C]// Proceedings of Environmental Mechanics：Water，Mass and Energy Transfer in the Biosphere. Washington：[s.n.]， 2002：213–246.
[8] WEI C F，DEWOOLKAR M M. Formulation of capillary hysteresis with internal state variables[J]. Water Resources Research，2006，42：1–16
[9] GRAY W G. Thermodynamics and constitutive theory for multiphase porous-media flow considering internal geometric constraints[J]. Advances in Water Resources，1999，22(5)：521–547.
[10] WEI C F. Static and dynamic behavior of multiphase porous media：governing equations and finite element implementation[Ph. D. Thesis][D]. Norman，Oklahoma：University of oklahoma，2001.
[11] ZIENKIEWICZ O C，SHIOMI T. Dynamic behavior of saturated porous media：the generalized Biot formulation and its numerical solution[J]. International Journal for Numerical and Analytical Methods in Geomechanics，1984，8(1)：71–96.
[12] MUALEM Y. A new model for predicting hydraulic conductivity of unsaturated porous-media[J]. Water Resources Research，1976，12(3)：513–522.
[13] WEI C F，MURALEETHARAN K K. Linear viscoelastic behavior of porous media with non-uniform saturation[J]. International Journal of Engineering Science，2007，45(2–8)：698–715.
[14] FENG M，FREDLUND D G. Hysteretic influence associated with thermal conductivity sensor measurements[C]// Proceedings of the 52nd Canada Geotechnical Conference and Unsaturated Soil Group from Theory to the Practice of Unsaturated Soil Mechanics，Regina，
Saskatchewan，Canada：[s. n.]，1999：651–657
[15] 徐炎兵，韦昌富，李 幻，等. 非饱和土渗流与变形耦合问题的有限元分析[J]. 岩土力学，2009，30(5)：1 490–1 496.(XU Yanbing，WEI Changfu，LI Huan，et al. Finite element analysis of coupling seepage and deformation in unsaturated soils[J]. Rock and Soil Mechanics，2009，30(5)：1 490–1 496.(in Chinese))
[16] HOA N T，GAUDU R，THIRRIOT C. Influence of hysteresis effect on transient flows in saturated-unsaturated porous-media[J]. Water Resources Research，1977，13(6)：992–996.
[17] HAVERKAMP R，ARRUE J L，SOET M. Soil physical properties within the root zone of the vine area of tomelloso[R]. Brussels，Spain： Local and Spatial Standpoint，1997： 144–185.
[18] RUBIN J. Numerical method for analyzing hysteresis-affected post-infiltration redistribution of soil moisture[J]. Soil Science Society of America Journal，1967，31：13–20.
[19] CHO S E，LEE S R. Instability of unsaturated soil slopes due to infiltration[J]. Computers and Geotechnics，2001，28(3)：185–208.
[20] NG C W W，SHI Q. A numerical investigation of the stability of unsaturated soil slopes subjected to transient seepage[J]. Computers and Geotechnics，1998，22(1)：1–28. |
|
|
|
2010, Vol. 29 
