Study on rainfall infiltration characteristics of unsaturated fractured soil based on COMSOL Multiphysics

A numerical simulation study on the rainfall infiltration characteristics of unsaturated fractured soil was performed based on COMSOL Multiphysics software. By discretizing the fracture and the matrix into finite elements, a discrete fracture-porous medium model was established to fully simulate the...

Full description

Bibliographic Details
Main Authors: HOU Xiao-ping, FAN Heng-hui
Format: Article
Language:English
Published: SCIENCE PRESS , 16 DONGHUANGCHENGGEN NORTH ST, BEIJING, PEOPLES R CHINA, 100717 2022-02-01
Series:Rock and Soil Mechanics
Subjects:
Online Access:http://rocksoilmech.whrsm.ac.cn/EN/10.16285/j.rsm.2021.5854
_version_ 1818015482310033408
author HOU Xiao-ping
FAN Heng-hui
author_facet HOU Xiao-ping
FAN Heng-hui
author_sort HOU Xiao-ping
collection DOAJ
description A numerical simulation study on the rainfall infiltration characteristics of unsaturated fractured soil was performed based on COMSOL Multiphysics software. By discretizing the fracture and the matrix into finite elements, a discrete fracture-porous medium model was established to fully simulate the fracture flow, matrix flow and fracture–matrix flow exchange in the soil. The upper boundary of the fractured soil was simulated by using the concept of “air element”. This method can describe not only the phenomenon of preferential infiltration for rainwater along the fractures at the beginning of rainfall, but also the phenomenon of rainwater flowing away along the surface when the rainfall is greater than the infiltration of the fractured soil. By simulating the low-permeability fractured soil at a depth of 2 m below the ground surface, the influences of the geometric characteristics of the fracture, the hydraulic properties of the matrix, the previous moisture condition and the rainfall intensity on the rainfall infiltration process of the unsaturated fractured soil were investigated. The results show two main seepage processes in the unsaturated fractured soil: firstly, water flows preferentially along the fracture; secondly, water is continuously imbibed into the matrix from the fracture, and the matrix imbibition inhibits the development of preferential flow in the fracture. Compared with the geometric characteristics of the fracture, the hydraulic properties of the matrix have a greater influence on the seepage flow in the unsaturated fractured soil. Increasing the saturated permeability coefficient of the matrix may change the seepage flow dominated by the fracture flow into that dominated by the matrix flow. The unsaturated properties of the matrix as well as the initial water content of the fractured soil change the soil moisture storage capacity, thereby accelerating or delaying the time for rainwater to infiltrate into a certain depth. The rainfall intensity has an influence on both the infiltration rate and infiltration amount in the soil. When it exceeds the infiltration capacity of the fractured soil, the excess rainwater flows away along the surface, and the infiltration rate across the section tends to stabilize with time.
first_indexed 2024-04-14T06:57:39Z
format Article
id doaj.art-e87c90b6124a4e91a5a789e97a59ffff
institution Directory Open Access Journal
issn 1000-7598
language English
last_indexed 2024-04-14T06:57:39Z
publishDate 2022-02-01
publisher SCIENCE PRESS , 16 DONGHUANGCHENGGEN NORTH ST, BEIJING, PEOPLES R CHINA, 100717
record_format Article
series Rock and Soil Mechanics
spelling doaj.art-e87c90b6124a4e91a5a789e97a59ffff2022-12-22T02:06:50ZengSCIENCE PRESS , 16 DONGHUANGCHENGGEN NORTH ST, BEIJING, PEOPLES R CHINA, 100717Rock and Soil Mechanics1000-75982022-02-0143256357210.16285/j.rsm.2021.5854Study on rainfall infiltration characteristics of unsaturated fractured soil based on COMSOL MultiphysicsHOU Xiao-ping0FAN Heng-hui1Key Laboratory of Agricultural Soil and Water Engineering in Arid and Semiarid Areas of Ministry of Education, Northwest A&F University, Yangling, Shaanxi 712100, ChinaKey Laboratory of Agricultural Soil and Water Engineering in Arid and Semiarid Areas of Ministry of Education, Northwest A&F University, Yangling, Shaanxi 712100, ChinaA numerical simulation study on the rainfall infiltration characteristics of unsaturated fractured soil was performed based on COMSOL Multiphysics software. By discretizing the fracture and the matrix into finite elements, a discrete fracture-porous medium model was established to fully simulate the fracture flow, matrix flow and fracture–matrix flow exchange in the soil. The upper boundary of the fractured soil was simulated by using the concept of “air element”. This method can describe not only the phenomenon of preferential infiltration for rainwater along the fractures at the beginning of rainfall, but also the phenomenon of rainwater flowing away along the surface when the rainfall is greater than the infiltration of the fractured soil. By simulating the low-permeability fractured soil at a depth of 2 m below the ground surface, the influences of the geometric characteristics of the fracture, the hydraulic properties of the matrix, the previous moisture condition and the rainfall intensity on the rainfall infiltration process of the unsaturated fractured soil were investigated. The results show two main seepage processes in the unsaturated fractured soil: firstly, water flows preferentially along the fracture; secondly, water is continuously imbibed into the matrix from the fracture, and the matrix imbibition inhibits the development of preferential flow in the fracture. Compared with the geometric characteristics of the fracture, the hydraulic properties of the matrix have a greater influence on the seepage flow in the unsaturated fractured soil. Increasing the saturated permeability coefficient of the matrix may change the seepage flow dominated by the fracture flow into that dominated by the matrix flow. The unsaturated properties of the matrix as well as the initial water content of the fractured soil change the soil moisture storage capacity, thereby accelerating or delaying the time for rainwater to infiltrate into a certain depth. The rainfall intensity has an influence on both the infiltration rate and infiltration amount in the soil. When it exceeds the infiltration capacity of the fractured soil, the excess rainwater flows away along the surface, and the infiltration rate across the section tends to stabilize with time.http://rocksoilmech.whrsm.ac.cn/EN/10.16285/j.rsm.2021.5854fractured soilunsaturated seepagerainfall infiltrationfracture-matrix flow exchangenumerical simulation
spellingShingle HOU Xiao-ping
FAN Heng-hui
Study on rainfall infiltration characteristics of unsaturated fractured soil based on COMSOL Multiphysics
Rock and Soil Mechanics
fractured soil
unsaturated seepage
rainfall infiltration
fracture-matrix flow exchange
numerical simulation
title Study on rainfall infiltration characteristics of unsaturated fractured soil based on COMSOL Multiphysics
title_full Study on rainfall infiltration characteristics of unsaturated fractured soil based on COMSOL Multiphysics
title_fullStr Study on rainfall infiltration characteristics of unsaturated fractured soil based on COMSOL Multiphysics
title_full_unstemmed Study on rainfall infiltration characteristics of unsaturated fractured soil based on COMSOL Multiphysics
title_short Study on rainfall infiltration characteristics of unsaturated fractured soil based on COMSOL Multiphysics
title_sort study on rainfall infiltration characteristics of unsaturated fractured soil based on comsol multiphysics
topic fractured soil
unsaturated seepage
rainfall infiltration
fracture-matrix flow exchange
numerical simulation
url http://rocksoilmech.whrsm.ac.cn/EN/10.16285/j.rsm.2021.5854
work_keys_str_mv AT houxiaoping studyonrainfallinfiltrationcharacteristicsofunsaturatedfracturedsoilbasedoncomsolmultiphysics
AT fanhenghui studyonrainfallinfiltrationcharacteristicsofunsaturatedfracturedsoilbasedoncomsolmultiphysics