An Efficient Numerical Integration Method for the Capillary Hysteresis Internal Variable Model
-
摘要: 毛细滞回内变量模型是基于热动力学基本原理推导出来的土水特征本构模型,能够有效地描述干湿循环情况下非饱和土中的毛细滞回现象.研究了该模型的形式特点和数值积分方法,就经典的Euler法、4阶Runge-Kutta法和4阶Adams-Bashforth法的积分效果进行了对比.计算结果表明,Euler法精度很低,累计误差较大;4阶Adams-Bashforth法精度最高,且运算效率高于同阶的Runge-Kutta法,适合对该模型进行求解和参数标定.将基于Adams-Bashforth法的子程序集成到有限元程序中,两相渗流模拟结果的精度有了较大提高.Abstract: Derived from thermo dynamic principle, the capillary hysteresis internal variable model is capable of describing the capillary hysteretic phenomena in unsaturated soils. The mathematical characteristics of this model were studied, followed by a numerical experimentation with classical Euler method, fourth-order Runge-Kutta method and fourth-order Adams-Bashforth method. The results show that Euler method has lower accuracy and larger accumulated error, whereas Adams-Bashforth method holds the upper most accuracy and the better efficiency compared with the same order Runge-Kutta method and is suitable for the solution and calibration of the internal variable model. Moreover, Adams-Bashforth method is implemented into the finite element programme, leading to more accurate results in simulation of two-phase flow.
-
Key words:
- capillary hysteresis /
- internal variable /
- numerical integration /
- two-phase flow /
- unsaturated soil
-
[1] Hillel D. Fundamentals of Soil Physics [M]. New York: Academic Press Inc, 1980. [2] Huang H C, Tan Y C, Liu C W, Chen C H. A novel hysteresis model in unsaturated soil[J]. Hydrological Processes,2005,19(8): 1653-1665. [3] Mualem Y A. Conceptual model of hysteresis[J]. Water Resources Research,1974, 10(3): 514-520. [4] Parlange J Y. Capillary hysteresis and the relationship between drying and wetting curves[J]. Water Resources Research,1976, 12(2): 224-228. [5] Wei C F, Dewoolkar M M. Formulation of capillary hysteresis with internal state variables[J]. Water Resources Research,2006, 42(7): W07405. [6] 陈盼, 魏厚振, 徐炎兵, 李幻, 韦昌富. 多孔介质中两相流动过程的毛细滞洄效应[J]. 岩石力学与工程学报, 2010,29(10): 21482158.(CHEN Pan, WEI Hou-zhen, XU Yan-bing, LI Huan, WEI Chang-fu. Effect of capillary hysteresis on twophase flow in porous media[J]. Chinese Journal of Rock Mechanics and Engineering,2010, 29(10): 2148-2158.(in Chinese)) [7] Muraleetharan K K, Liu C, Wei C, Kibbey T C G, Chen L. An elastoplatic framework for coupling hydraulic and mechanical behavior of unsaturated soils[J]. International Journal of Plasticity,2009, 25(3): 473-490. [8] Liu C,Muraleetharan K K.Coupled hydroMmechanical elastoplastic constitutive model for unsaturated sands and silts—I: formulation[J]. Int J Geomech,2012, 12(3): 239-247. [9] Dafalias Y F. Bounding surface of plasticity—I:mathematical foundation and hypoplasticity[J]. J Eng Mech,1986, 112(9): 967-987. [10] Vogel T, van Genuchten M T, Cislerova M.Effect of the shape of the soil hydraulic functions near saturation on variablysaturated flow predictions[J].Advances in Water Resources,2000, 24(2): 133-144. [11] 李庆扬, 王能超, 易大义.数值分析[M]. 北京:清华大学出版社, 2008.(LI Qing-yang, WANG Neng-chao, YI Da-yi. Numerical Analysis [M]. Beijing: Tsinghua University Press, 2008.(in Chinese)) [12] Hoa N T, Gaudu R, Thirriot C. Influence of the hysteresis effect on transient flows in saturatedunsaturated porous media[J]. Water Resources Research,1977, 13(6): 992-996. [13] 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. [14] Wei C F. Static and dynamic behavior of multiphase porous media: governing equations and finite element implementation[D]. Norman: University of Oklahoma, 2001.
点击查看大图
计量
- 文章访问数: 1592
- HTML全文浏览量: 47
- PDF下载量: 1077
- 被引次数: 0