Investigation of a Voronoi-Based Hybrid Seepage Flux Finite Element Method

LUO Chao; ZHANG Rui; GUO Ran; GAI Wenhai

doi:10.21656/1000-0887.460148

Volume 47 Issue 6

Jun. 2026

Turn off MathJax

Article Contents

Article Navigation > Applied Mathematics and Mechanics > 2026 > 47(6): 712-722

LUO Chao, ZHANG Rui, GUO Ran, GAI Wenhai. Investigation of a Voronoi-Based Hybrid Seepage Flux Finite Element Method[J]. Applied Mathematics and Mechanics, 2026, 47(6): 712-722. doi: 10.21656/1000-0887.460148

Citation:

PDF( 4026 KB)

Investigation of a Voronoi-Based Hybrid Seepage Flux Finite Element Method

doi: 10.21656/1000-0887.460148

LUO Chao^1
,,
ZHANG Rui^{1
,
,},
GUO Ran¹,
GAI Wenhai²

1.
Faculty of Civil Engineering and Mechanics, Kunming University of Science and Technology, Kunming 650500, P.R. China
2.
School of Highway and Construction Engineering, Yunnan Communications Vocational and Technical College, Kunming 650500, P.R. China

Received Date: 2025-08-19
Rev Recd Date: 2025-09-12
Publish Date: 2026-06-01

Abstract

Abstract

The soil-rock mixtures and naturally anisotropic dam foundations were treated as composite porous media, and the seepage flow was assumed to follow 2D Darcy's law with continuous hydraulic head distributions within the soil-rock mixture. A hybrid seepage flux finite element method (HS-FEM) model was developed for analyzing seepage fields in soil-rock mixtures with interfaces and anisotropic homogeneous dam foundations. For this model, independently assumed higher-order seepage flux variables was adopted within element domains with hydraulic head values prescribed on element boundaries. Only a limited number of elements were required to effectively simulate confined seepage conditions involving soil-rock interfaces, thereby to overcome the drawback of traditional FEMs necessitating dense mesh refinement at material interfaces. Additionally, the method is capable of solving 2D orthotropic steady-state linear seepage problems. Numerical examples demonstrate that the proposed HS-FEM achieves comparable accuracy to traditional dense-mesh FEMs while maintaining computational efficiency through sparse discretization.
- hybrid seepage flux finite element method (HS-FEM),
- seepage analysis,
- composite material

(Contributed by GUO Ran, Member of the Editorial Board of AMM)

FullText(HTML)

References(31)

References

[1]	NEUMAN S P, WITHERSPOON P A. Finite element method of analyzing steady seepage with a free surface[J]. Water Resources Research, 1970, 6(3): 889-897. doi: 10.1029/WR006i003p00889
[2]	ZHANG W, DAI B, LIU Z, et al. Unconfined seepage analysis using moving Kriging mesh-free method with Monte Carlo integration[J]. Transport in Porous Media, 2017, 116(1): 163-180. doi: 10.1007/s11242-016-0769-9
[3]	DAI Q, LEI Y, ZHANG B, et al. A practical adaptive moving-mesh algorithm for solving unconfined seepage problem with Galerkin finite element method[J]. Scientific Reports, 2019, 9: 6988. doi: 10.1038/s41598-019-43391-4
[4]	SHARMA V, FUJISAWA K, MURAKAMI A. Space-time finite element method for transient and unconfined seepage flow analysis[J]. Finite Elements in Analysis and Design, 2021, 197: 103632. doi: 10.1016/j.finel.2021.103632
[5]	CHUGH A K, FALVEY H T. Seepage analysis in a zoned anisotropic medium by the boundary element method[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1984, 8(4): 399-407. doi: 10.1002/nag.1610080407
[6]	柴军瑞, 仵彦卿. 采用边界元法确定渗流自由面的一种改进方法[J]. 陕西水力发电, 2000(1): 11-12. CHAI Junrui, WU Yanqing. A modified iterative approach for solving seepage free surface by BEM[J]. Journal of Shaanxi Water Power, 2000(1): 11-12. (in Chinese)
[7]	葛锦宏, 李广信, 介玉新. 无单元法在有自由面渗流计算中的应用[J]. 计算力学学报, 2003, 20(2): 241-243. GE Jinhong, LI Guangxin, JIE Yuxin. Application of element-free method to the seepage calculation with free surface[J]. Chinese Journal of Computational Mechanics, 2003, 20(2): 241-243. (in Chinese)
[8]	张石峰, 李茜, 高佩玲. 天然差分方法在地下水渗流问题中的应用[J]. 计算物理, 2007, 24(3): 307-312. ZHANG Shifeng, LI Qian, GAO Peiling. Natural difference method for groundwater seepage[J]. Chinese Journal of Computational Physics, 2007, 24(3): 307-312. (in Chinese)
[9]	BARDET J P, TOBITA T. A practical method for solving free-surface seepage problems[J]. Computers and Geotechnics, 2002, 29(6): 451-475. doi: 10.1016/S0266-352X(02)00003-4
[10]	李凤志. 渗流自由面分析的比例边界有限元法[J]. 计算物理, 2009, 26(5): 665-670. LI Fengzhi. Scaled boundary finite-element method for seepage free surfaces analysis[J]. Chinese Journal of Computational Physics, 2009, 26(5): 665-670. (in Chinese)
[11]	李珏, 金吾根, 薛馨, 等. Trefftz间接法解自由面渗流问题[J]. 计算力学学报, 2005, 22(3): 295-298. LI Jue, JIN Wugen, XUE Xin, et al. Solution of free surface seepage problems by Trefftz indirect method[J]. Chinese Journal of Computational Mechanics, 2005, 22(3): 295-298. (in Chinese)
[12]	薛馨, 金吾根, 李珏, 等. Trefftz直接法解渗流自由面问题[J]. 岩石力学与工程学报, 2005, 24(13): 2322-2326. XUE Xin, JIN Wugen, LI Jue, et al. Solution of free surface seepage problems by Trefftz direct method[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(13): 2322-2326. (in Chinese)
[13]	徐文杰, 王永刚. 土石混合体细观结构渗流数值试验研究[J]. 岩土工程学报, 2010, 32(4): 542-550. XU Wenjie, WANG Yonggang. Meso-structural permeability of S-RM based on numerical tests[J]. Chinese Journal of Geotechnical Engineering, 2010, 32(4): 542-550. (in Chinese)
[14]	YAN L, MENG Q X, XU W Y, et al. A numerical method for analyzing the permeability of heterogeneous geomaterials based on digital image processing[J]. Journal of Zhejiang University (Science A): Applied Physics & Eingineering, 2017, 18(2): 124-137.
[15]	WEI W, JIANG Q, YE Z, et al. Equivalent fracture network model for steady seepage problems with free surfaces[J]. Journal of Hydrology, 2021, 603: 127156. doi: 10.1016/j.jhydrol.2021.127156
[16]	陈远强, 郑宏, 屈新. 基于数值流形法的降雨入渗与坡面径流耦合算法研究[J]. 应用数学和力学, 2023, 44(12): 1499-1511. CHEN Yuanqiang, ZHENG Hong, QU Xin. A coupling analysis of rainfall infiltration and slope surface runoff based on the numerical manifold method[J]. Applied Mathematics and Mechanics, 2023, 44(12): 1499-1511. (in Chinese)
[17]	JIA Z, ZHENG H. A new procedure for locating free surfaces of complex unconfined seepage problems using fixed meshes[J]. Computers and Geotechnics, 2024, 166: 106032. doi: 10.1016/j.compgeo.2023.106032
[18]	ZHAO L, ZHANG H, MAO J, et al. Numerical simulation of the bubbling phenomenon ahead of the dam using a two-phase seepage free surface model[J]. Computers and Geotechnics, 2023, 155: 105239. doi: 10.1016/j.compgeo.2022.105239
[19]	ZHAO L, ZHANG H, MAO J, et al. An ICLS-based method for solving two-phase seepage free surface considering compressible gas in porous media[J]. Computers and Geotechnics, 2022, 141: 104528. doi: 10.1016/j.compgeo.2021.104528
[20]	李树光, 曲凯. 多孔介质中单相气体局部流动的均质化建模[J]. 应用数学和力学, 2024, 45(2): 175-183. doi: 10.21656/1000-0887.440246 LI Shuguang, QU Kai. Homogenization modeling of single-phase gas local flow in porous media[J]. Applied Mathematics and Mechanics, 2024, 45(2): 175-183. (in Chinese) doi: 10.21656/1000-0887.440246
[21]	PIAN T H H. Derivation of element stiffness matrices by assumed stress distributions[J]. AIAA Journal, 1964, 2(7): 1333-1336. doi: 10.2514/3.2546
[22]	PIN T, PIAN T H H. A variational principle and the convergence of a finite-element method based on assumed stress distribution[J]. International Journal of Solids and Structures, 1969, 5(5): 463-472. doi: 10.1016/0020-7683(69)90036-5
[23]	PIAN T H H, TONG P. Basis of finite element methods for solid continua[J]. International Journal for Numerical Methods in Engineering, 1969, 1(1): 3-28. doi: 10.1002/nme.1620010103
[24]	GHOSH S, MUKHOPADHYAY S N. A material based finite element analysis of heterogeneous media involving Dirichlet tessellations[J]. Computer Methods in Applied Mechanics and Engineering, 1993, 104(2): 211-247. doi: 10.1016/0045-7825(93)90198-7
[25]	GHOSH S, MALLETT R L. Voronoi cell finite elements[J]. Computers & Structures, 1994, 50(1): 33-46.
[26]	GUO R, SHI H J, YAO Z H. Modeling of interfacial debonding crack in particle reinforced composites using Voronoi cell finite element method[J]. Computational Mechanics, 2003, 32(1): 52-59.
[27]	ZHANG R, GUO R. Voronoi cell finite element method for fluid-filled materials[J]. Transport in Porous Media, 2017, 120(1): 23-35. doi: 10.1007/s11242-017-0898-9
[28]	ZHANG R, GUO R. Determination of crack tip stress intensity factors by singular Voronoi cell finite element model[J]. Engineering Fracture Mechanics, 2018, 197: 206-216. doi: 10.1016/j.engfracmech.2018.04.042
[29]	KARCHEWSKI B. Multi-scale modelling of geomechanical behaviour using the Voronoi cell finite element method (VCFEM) and finite-discrete element method (VCFEM-DEM)[D]. Hamilton, Ontario: McMaster University, 2015.
[30]	FELIPPA C A. Parametrized multifield variational principles in elasticity, Ⅰ: mixed functionals[J]. Communications in Applied Numerical Methods, 1989, 5(2): 79-88. doi: 10.1002/cnm.1630050204
[31]	FELIPPA C A. A survey of parametrized variational principles and applications to computational mechanics[J]. Computer Methods in Applied Mechanics and Engineering, 1994, 113(1/2): 109-139.