Sciences in Cold and Arid Regions ›› 2019, Vol. 11 ›› Issue (1): 21-28.doi: 10.3724/SP.J.1226.2019.00021

Previous Articles     Next Articles

The analysis of heat and water fluxes in frozen silty soil

DaHu Rui1,2,*(),Ming Lu1,Kunio Watanabe3,Jun Zhang2,4   

  1. 1. School of Civil Engineering, Henan Polytechnic University, Jiaozuo, Henan 454000, China
    2. State Key Laboratory of Frozen Soil Engineering, Cold and Arid Regions Environmental and Engineering Research Institute, Chinese Academy of Sciences, Lanzhou, Gansu 730000, China
    3. Graduate School of Bioresources, Mie University. Mie 514-8507, Japan
    4. University of Chinese Academy of Sciences, Beijing 100049, China
  • Received:2018-09-16 Accepted:2018-10-09 Online:2019-02-01 Published:2019-03-22
  • Contact: DaHu Rui
  • About author:DaHu Rui, School of Civil Engineering, Henan Polytechnic University, Century Road, No. 2001, Shanyang District, Jiaozuo, Henan 454000, China.


In this paper, based on the basic equations of water flow and heat transfer, the hydrothermal coupling model is established. The numerical model was realized in COMSOL Multiphysics software, and simulation results are compared with the experimental results of Watanabe and Wake (2008) to verify the effectiveness of the model. Through the calculation, we can obtain the dynamic changes of heat and water fluxes, thermal and hydrological properties, matric potential and temperature gradient in unsaturated freezing soil; and these variables are unmeasurable in practice.

Key words: unsaturated freezing soil, hydrothermal coupling, heat and water fluxes, COMSOL

Figure 1

Plots (a) and (b) indicate the experimental results (Watanabe and Wake, 2008) and simulation results of the distribution of the temperature profile, respectively"

Figure 2

Plots (a) and (b) indicate the experimental results and simulation results of the distribution of the liquid-water content and total-water content in the unidirectional freezing experiment, respectively"

Figure 3

Profile of the temperature gradient at different times"

Figure 4

Curves of temperature vs. time at different depths"

Figure 5

Distribution of ice content and thermal conductivity during soil freezing"

Figure 6

Distribution of the volumetric apparent heat capacity vs. time at different depths"

Figure 7

Distribution of the heat flux at different times"

Figure 8

Distribution of the pressure head at different times"

Figure 9

(a) Distribution of the water flux at different times and (b) distribution of the water flux vs. time at different depths"

Campbell GS , 1985. Soil Physics with BASIC. Elsevier Elsevier, New York .
Fayer MJ, 2000. UNSAT-H version 3.0: unsaturated soil water and heat flow model theory, user manual, and examples. Richland, WA: Pacific Northwest National Lab. DOI: 10.2 172/15001068.
doi: 10.2 172/15001068.
Hansson K , Šimůnek J , Mizoguchi M , et al ., 2004. Water flow and heat transport in frozen soil. Vadose Zone Journal, 3(2): 693−704. DOI: 10.2136/vzj2004.0693.
doi: 10.2136/vzj2004.0693.
Harlan RL, 1973. Analysis of coupled heat-fluid transport in partially frozen soil. Water Resources Research, 9(5): 1314−1323. DOI: 10.1029/WR009i005p01314.
doi: 10.1029/WR009i005p01314.
Karra S , Painter SL , Lichtner PC , 2014. Three-phase numerical model for subsurface hydrology in permafrost-affected regions. The Cryosphere, (2014−10−23), 8(1):1935−1950.
Koopmans RWR , Miller RD , 1966. Soil freezing and soil water characteristic curves. Soil Science Society of America Journal, 30(6): 680−685. DOI: 10.2136/sssaj1966.0361599500 3000060011x.
doi: 10.2136/sssaj1966.0361599500 3000060011x.
Liu Z , Yu X , 2011. Coupled thermo-hydro-mechanical model for porous materials under frost action: theory and implementation. Acta Geotechnica, 6(2): 51−65. DOI: 10.1007/s11440-011-0135-6.
doi: 10.1007/s11440-011-0135-6.
Lundin LC, 1990. Hydraulic properties in an operational model of frozen soil. Journal of Hydrology, 118(1−4): 289−310. DOI: 10.1016/0022-1694(90)90264-X.
doi: 10.1016/0022-1694(90)90264-X.
Nassar IN , Horton R , 1989. Water Transport in Unsaturated Nonisothermal Salty Soil: II. Theoretical Development. Soil Science Society of America Journal, 53(5): 1330−1337.
Newman GP , Wilson GW , 1997. Heat and mass transfer in unsaturated soils during freezing. Canadian Geotechnical Journal, 34(1): 63−70. DOI: 10.1139/t96-085.
doi: 10.1139/t96-085.
Saito H , Šimůnek J , Mohanty BP , 2006. Numerical analysis of coupled water, vapor, and heat transport in the vadose zone. Vadose Zone Journal, 5(2): 784−800. DOI: 10.2136/vzj20 06.0007.
doi: 10.2136/vzj20 06.0007.
Spaans EJA , Baker JM , 1996. The soil freezing characteristic: its measurement and similarity to the soil moisture characteristic. Soil Science Society of America Journal, 60(1): 13−19. DOI: 10.2136/sssaj1996.03615995006000010005x.
doi: 10.2136/sssaj1996.03615995006000010005x.
Tan XJ , Chen W , Tian H , et al ., 2011. Water flow and heat transport including ice/water phase change in porous media: Numerical simulation and application. Cold Regions Science & Technology, 68(1−2): 74−84.
Van Genuchten MT, 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Science Society of America Journal, 44(5): 892−898. DOI: 10.2136/sssaj1980.03615995004400050002x.
doi: 10.2136/sssaj1980.03615995004400050002x.
Watanabe, 2008. Water and heat flow in a directionally frozen silty soil. In: Saito H , Sakai M , Toride N , et al . (eds.). Proceedings of the Third HYDRUS Workshop. Tokyo, Japan: Tokyo University of Agriculture and Technology.
Watanabe K , Flury M , 2008. Capillary bundle model of hydraulic conductivity for frozen soil. Water Resources Research, 44(12): W12402. DOI: 10.1029/2008WR007012.
doi: 10.1029/2008WR007012.
Watanabe K , Wake T , 2008. Hydraulic conductivity in frozen unsaturated soil. In: Proceedings of the 9th International Conference on Permafrost. Fairbanks: University of Alaska Fairbanks.
Williams PJ , Smith MW , 1989. The Frozen Earth—Fundamentals of Geocryology. Cambridge University Press, Oxford.
Watanabe K , Osada Y , 2017. Simultaneous measurement of unfrozen water content and hydraulic conductivity of partially frozen soil near 0 °C. Cold Regions Science and Technology, 142: 79−84. DOI: 10.1016/j.coldregions.2017.08.002.
doi: 10.1016/j.coldregions.2017.08.002.
Williams PJ, 1964. Unfrozen water content of frozen soils and soil moisture suction. Géotechnique, 14(3): 231−246. DOI: 10.1680/geot.1964.14.3.231.
doi: 10.1680/geot.1964.14.3.231.
Xu XZ , Wang JC , Zhang LX , 2001. Permafrostphysics. Beijing: Science Press, pp. 75−98.
No related articles found!
Full text



[1] Mohan Bahadur Chand,Rijan Bhakta Kayastha. Study of thermal properties of supraglacial debris and degree-day factors on Lirung Glacier, Nepal[J]. Sciences in Cold and Arid Regions, 2018, 10(5): 357 -368 .
[2] AiHong Xie,ShiMeng Wang,YiCheng Wang,ChuanJin Li. Comparison of temperature extremes between Zhongshan Station and Great Wall Station in Antarctica[J]. Sciences in Cold and Arid Regions, 2018, 10(5): 369 -378 .
[3] YanZai Wang,YongQiu Wu,MeiHui Pan,RuiJie Lu. Comparison of two classification methods to identify grain size fractions of aeolian sediment[J]. Sciences in Cold and Arid Regions, 2018, 10(5): 413 -420 .
[4] YinHuan Ao,ShiHua Lyu,ZhaoGuo Li,LiJuan Wen,Lin Zhao. Numerical simulation of the climate effect of high-altitude lakes on the Tibetan Plateau[J]. Sciences in Cold and Arid Regions, 2018, 10(5): 379 -391 .
[5] Zhuo Ga,Za Dui,Duodian Luozhu,Jun Du. Comparison of precipitation products to observations in Tibet during the rainy season[J]. Sciences in Cold and Arid Regions, 2018, 10(5): 392 -403 .
[6] Rong Yang,JunQia Kong,ZeYu Du,YongZhong Su. Altitude pattern of carbon stocks in desert grasslands of an arid land region[J]. Sciences in Cold and Arid Regions, 2018, 10(5): 404 -412 .
[7] Yang Qiu,ZhongKui Xie,XinPing Wang,YaJun Wang,YuBao Zhang,YuHui He,WenMei Li,WenCong Lv. Effect of slow-release iron fertilizer on iron-deficiency chlorosis, yield and quality of Lilium davidii var. unicolor in a two-year field experiment[J]. Sciences in Cold and Arid Regions, 2018, 10(5): 421 -427 .
[8] Ololade A. Oyedapo,Joseph M. Agbedahunsi,H. C Illoh,Akinwumi J. Akinloye. Comparative foliar anatomy of three Khaya species (Meliaceae) used in Nigeria as antisickling agent[J]. Sciences in Cold and Arid Regions, 2018, 10(4): 279 -285 .
[9] YuMing Wei,XiaoFei Ma,PengShan Zhao. Transcriptomic comparison to identify rapidly evolving genes in Braya humilis[J]. Sciences in Cold and Arid Regions, 2018, 10(5): 428 -435 .
[10] FangLei Zhong,AiJun Guo,XiaoJuan Yin,JinFeng Cui,Xiao Yang,YanQiong Zhang. Sociodemographic characteristics, cultural biases, and environmental attitudes: An empirical application of grid-group cultural theory in Northwestern China[J]. Sciences in Cold and Arid Regions, 2018, 10(5): 436 -446 .