Chuhao Huang , Jun Liu , Zhiqiang Hu , Lei Xiong , Zhi Liu , Lei Gan , Wenbin Ye , Gao Lin
{"title":"Predicting steady-state seepage flow in complex half-space using a modified scaled boundary finite element method","authors":"Chuhao Huang , Jun Liu , Zhiqiang Hu , Lei Xiong , Zhi Liu , Lei Gan , Wenbin Ye , Gao Lin","doi":"10.1016/j.compgeo.2025.107228","DOIUrl":null,"url":null,"abstract":"<div><div>Seepage problems in half-space domains are crucial in hydrology, environmental, and civil engineering, involving groundwater flow, pollutant transport, and structural stability. Typical examples include seepage through dam foundations, coastal aquifers, and levees under seepage forces, requiring accurate numerical modeling. However, existing methods face challenges in handling complex geometries, heterogeneous media, and anisotropic properties, particularly in multi-domain half-spaces. This study addresses these challenges by extending the modified scaled boundary finite element method (SBFEM) and using this method to explore steady seepage problems in complex half-space domain. In the modified SBFEM framework, segmented straight lines or curves, parallel to the far-field infinite boundary, are introduced as scaling lines, with a one-dimensional discretization applied to them, thereby reducing computational costs.Then the weighted residual method is applied to obtain the modified SBFEM governing equations and boundary conditions of steady-state seepage problem according to the Laplace diffusion equation and Darcy’s law. Furthermore, the steady seepage matrix at infinity is obtained by solving the eigenvalue problem of Schur decomposition and then the 4th-order Runge-Kutta algorithm is used to iteratively solve until the seepage matrix at the boundary lines is reached. Comparisons between the present numerical results and solutions available in the published work have been conducted to demonstrate the efficiency and accuracy of this method. At the same time, the influences of the geometric parameters and complex half-space domain on the seepage flow characteristics in complex half-space domain are investigated in detail.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"183 ","pages":"Article 107228"},"PeriodicalIF":5.3000,"publicationDate":"2025-03-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computers and Geotechnics","FirstCategoryId":"5","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0266352X25001776","RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
Seepage problems in half-space domains are crucial in hydrology, environmental, and civil engineering, involving groundwater flow, pollutant transport, and structural stability. Typical examples include seepage through dam foundations, coastal aquifers, and levees under seepage forces, requiring accurate numerical modeling. However, existing methods face challenges in handling complex geometries, heterogeneous media, and anisotropic properties, particularly in multi-domain half-spaces. This study addresses these challenges by extending the modified scaled boundary finite element method (SBFEM) and using this method to explore steady seepage problems in complex half-space domain. In the modified SBFEM framework, segmented straight lines or curves, parallel to the far-field infinite boundary, are introduced as scaling lines, with a one-dimensional discretization applied to them, thereby reducing computational costs.Then the weighted residual method is applied to obtain the modified SBFEM governing equations and boundary conditions of steady-state seepage problem according to the Laplace diffusion equation and Darcy’s law. Furthermore, the steady seepage matrix at infinity is obtained by solving the eigenvalue problem of Schur decomposition and then the 4th-order Runge-Kutta algorithm is used to iteratively solve until the seepage matrix at the boundary lines is reached. Comparisons between the present numerical results and solutions available in the published work have been conducted to demonstrate the efficiency and accuracy of this method. At the same time, the influences of the geometric parameters and complex half-space domain on the seepage flow characteristics in complex half-space domain are investigated in detail.
期刊介绍:
The use of computers is firmly established in geotechnical engineering and continues to grow rapidly in both engineering practice and academe. The development of advanced numerical techniques and constitutive modeling, in conjunction with rapid developments in computer hardware, enables problems to be tackled that were unthinkable even a few years ago. Computers and Geotechnics provides an up-to-date reference for engineers and researchers engaged in computer aided analysis and research in geotechnical engineering. The journal is intended for an expeditious dissemination of advanced computer applications across a broad range of geotechnical topics. Contributions on advances in numerical algorithms, computer implementation of new constitutive models and probabilistic methods are especially encouraged.