Computers and Geotechnics最新文献

{"title":"A highly scalable numerical framework for reservoir simulation on UG4 platform","authors":"Shuai Lu","doi":"10.1016/j.compgeo.2025.107597","DOIUrl":"10.1016/j.compgeo.2025.107597","url":null,"abstract":"<div><div>The modeling and simulation of multiphase fluid flow have received significant attention from the reservoir engineering research community. Many time discretization schemes for multiphase flow equations are either explicit or semi-implicit, relying on the decoupling between the saturation equation and the pressure equation. In this study, we delve into a fully coupled and fully implicit framework for simulating multiphase flow in heterogeneous porous media, considering both the gravity and capillary effects. We utilize the Vertex-Centered Finite Volume Method for spatial discretization and propose an efficient implementation of capillary barrier condition for heterogeneous porous media within the current scheme. Notably, we introduce the Linearly Implicit Extrapolation Method (LIMEX) with an error estimator, adapted for the first time to multiphase flow problems. To solve the resulting linear system, we employ the BiCGSTAB method with the Geometric Multigrid (GMG) preconditioner. The implementations of the models and methods are based on the open-source software: UG4. The results from parallel computations on the supercomputer demonstrate that the scalability of our proposed framework is sufficient, supporting a scale of thousands of processors with Degrees of Freedom (DoF) extending into the billions.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"188 ","pages":"Article 107597"},"PeriodicalIF":6.2,"publicationDate":"2025-09-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144932072","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Multiscale numerical approach to assess impact forces on cylindrical barriers by debris flows","authors":"Woojae Jang ,&nbsp;Hwan-Hui Lim ,&nbsp;Enok Cheon ,&nbsp;Tae-Hyuk Kwon ,&nbsp;Chan-Young Yune ,&nbsp;Beom-Jun Kim ,&nbsp;Shin-Kyu Choi","doi":"10.1016/j.compgeo.2025.107605","DOIUrl":"10.1016/j.compgeo.2025.107605","url":null,"abstract":"<div><div>Debris flow is a flow-like landslide involving a mixture of water and earth materials that rapidly moves downslope. Conducting real-scale experiments is costly and resource-intensive, making numerical analysis an efficient alternative for evaluating the performance of countermeasures. This study investigates impact forces on barriers caused by debris flows using a multiscale numerical approach. First, a 3D depth-averaged model simulates debris flow characteristics at a channel-length scale. Then, a 3D Smoothed Particle Hydrodynamics (SPH) model examines debris-barrier interactions at a channel-width scale. Two real-scale debris flow experiments conducted in Pyeongchang, Korea, were used for validation: one without barriers for flow calibration and another with cylindrical baffle-type countermeasures for impact force validation. The depth-averaged model was calibrated to replicate observed runout distances and velocities, while the SPH model simulated interactions with cylindrical baffles. Due to the curved terrain near the first baffle array, debris primarily impacts outer baffles, resulting in higher dynamic impact coefficients. The parametric study reveals that the initial velocity has a negative correlation with the dynamic impact coefficient but a positive correlation with the Froude number. Increased density reduces both the dynamic impact coefficient and pressure gradient. In addition, previously proposed semi-empirical models effectively reproduce the summed impact forces across all baffles upon calibration against the numerical results, but capturing individual baffle’s responses and estimating localized impact forces in open-type barriers still remain challenging. These findings highlight the critical influence of local terrain and flow dynamics on impact force distribution. demonstrating the values of combining numerical simulations and real-scale experiments in geohazard risk assessment and mitigation strategy development.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"188 ","pages":"Article 107605"},"PeriodicalIF":6.2,"publicationDate":"2025-09-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144925234","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Fracture behavior of heterogeneous layered sandstone under mixed loading after thermal treatment: A thermal–mechanical grain-based FDEM model","authors":"Wei Kuang ,&nbsp;Tengfei Guo ,&nbsp;Xu Chang ,&nbsp;Yichao Rui","doi":"10.1016/j.compgeo.2025.107599","DOIUrl":"10.1016/j.compgeo.2025.107599","url":null,"abstract":"<div><div>In deep rock engineering, temperature and bedding plane properties have a critical influence on the mechanical behavior of layered sandstone. Therefore, this study develops a cracked straight through Brazilian disc (CSTBD) numerical model for layered sandstone using the combined finite-discrete element method (FDEM), incorporating material inhomogeneity and bedding plane characteristics. The fracture behavior of layered sandstone under splitting conditions after high-temperature treatment was examined, with model validation achieved through laboratory experiments. The results indicate that the microcrack density in layered sandstone increases and the crack distribution becomes more disordered as the temperature rises. During the heating phase, initial damage predominantly appears as tensile cracks, whereas in the cooling phase, the damage shifts to shear cracks, with the growth rates of microcracks generally slowing down. On the microscale, high thermal treatment temperatures cause thermal cracking of mineral grains, which damages the bedding planes. On the macroscale, thermal damage results in a reduction of peak load, with a notable transition from brittle to ductile fracture occurring between 600 °C and 800 °C. The crack propagation in layered sandstone is synergistically influenced by a combination of factors, including thermal treatment temperature, loading angle, bedding plane strength, and prefabricated crack length. The fracture toughness diminishes with rising temperatures or prefabricated crack length, and rises with increasing strength of the bedding plane. Furthermore, the effects of bedding plane strength ratio, bedding plane distribution pattern, and prefabricated crack length on the fracture toughness of layered sandstone were examined under diverse thermal treatment temperatures and loading modes. The results indicate that these factors exert significant control over the fracture mode and fracture resistance of layered sandstone.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"188 ","pages":"Article 107599"},"PeriodicalIF":6.2,"publicationDate":"2025-08-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144916771","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Intelligent inversion analysis of seepage parameters for deep overburden dam foundations based on an improved grey wolf optimization algorithm","authors":"Qiaogang Yin ,&nbsp;Yanlong Li ,&nbsp;Wenwei Li ,&nbsp;Lifeng Wen ,&nbsp;Ye Zhang ,&nbsp;Ting Wang ,&nbsp;Tao Yang ,&nbsp;Tao Zhou","doi":"10.1016/j.compgeo.2025.107595","DOIUrl":"10.1016/j.compgeo.2025.107595","url":null,"abstract":"<div><div>The seepage parameters of dam foundations are essential for conducting three-dimensional seepage safety analyses. However, traditional methods often fail to achieve efficient and accurate parameter inversion because of the challenges posed by complex geological conditions and limited in situ measurements. To address this issue, a novel dual-layer collaborative intelligent inversion framework was proposed. This approach integrated an orthogonal experimental design with finite element modeling to construct a training dataset. An Improved Chaotic Grey Wolf Optimization (ICGWO) algorithm was then employed to optimize the hyperparameters of a Light Gradient Boosting Machine (LightGBM), resulting in the development of a high-precision surrogate model. Within a reasonable range of seepage parameter values, ICGWO was further applied to explore the solution space of the surrogate model, which enabled intelligent inversion of permeability parameters. A case study conducted on a deep overburden dam foundation at the Y Hydropower Station demonstrated that the ICGWO-LightGBM model accurately predicted borehole water levels. Moreover, ICGWO exhibited notable efficiency in optimizing stratified permeability parameters. Forward simulations using the permeability coefficients derived from the inversion process yielded a maximum absolute error of 7.92 m and a relative error of only 0.26 % in borehole water level predictions. The simulated seepage field displayed physically consistent distribution patterns that aligned with the typical mountain-type seepage behavior, confirming the accuracy, robustness, and engineering applicability of the method. The proposed approach can provide an efficient and innovative pathway for dam seepage safety assessment under complex geological conditions, offering considerable theoretical value and broad potential for engineering practice.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"188 ","pages":"Article 107595"},"PeriodicalIF":6.2,"publicationDate":"2025-08-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144916772","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A refined fictitious-pile-method to evaluate the vibration mitigation effect of pile arrays in a layered half-space","authors":"Chao He ,&nbsp;Zeyu Zhang ,&nbsp;Yuhao Peng ,&nbsp;Shunhua Zhou","doi":"10.1016/j.compgeo.2025.107578","DOIUrl":"10.1016/j.compgeo.2025.107578","url":null,"abstract":"<div><div>This study focuses on the mitigation of ground-borne vibrations with a frequency range of 1–80 Hz by pile arrays, consistent with typical rail transit systems. A refined fictitious pile method is therefore developed to calculate vibration mitigation performance of pile arrays in horizontally stratified soils. The proposed semi-analytical framework accurately captures wave scattering at the pile–soil interface without relying on periodic or symmetric simplifications. It enables the simulation of arbitrarily distributed pile arrays with varying dimensions and material properties. The thin-layer method is applied to calculate the Green’s function for a layered half-space in the space–frequency domain, eliminating the need for double discrete Fourier transforms of spatial coordinates. Consequently, the proposed method demonstrates superior computational efficiency. The validation of the proposed approach was conducted through comparative analysis with existing numerical methods. Numerical studies reveal the interplay between Bragg scattering and local resonance mechanisms, highlighting the influence of pile dimensions and array characteristics on vibration mitigation performance. Soil stratification exhibits a significant influence on the mitigation performance of pile arrays. Detailed dispersion analysis and circular interface statistical interference mapping elucidate the underlying physical phenomena governing the formation of isolation bands.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"188 ","pages":"Article 107578"},"PeriodicalIF":6.2,"publicationDate":"2025-08-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144916769","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Material characteristic length insensitive nonlocal modelling: A computationally efficient scaled nonlocal integral method","authors":"Xin Zhou ,&nbsp;Dechun Lu ,&nbsp;Jidong Zhao ,&nbsp;Yaning Zhang ,&nbsp;Zhiwei Gao ,&nbsp;Timon Rabczuk ,&nbsp;Xiuli Du","doi":"10.1016/j.compgeo.2025.107587","DOIUrl":"10.1016/j.compgeo.2025.107587","url":null,"abstract":"<div><div>Nonlocal modelling has achieved notable progress in resolving mesh dependence but remains constrained by two persistent challenges: sensitivity to characteristic length parameters and high computational costs. This study presents a scaled nonlocal integral formulation coupled with an optimized computational framework to simultaneously address two limitations. We first analytically demonstrate that variations in characteristic length induce proportional scaling of load–displacement curves, revealing that apparent changes in structure softening rate are artifacts of this scaling. Building on this insight, a dimensionless scaling factor is derived to systematically eliminate characteristic length dependence, enabling consistent predictions across parameter choices. The proposed method is integrated with a Mohr-Coulomb plasticity damage model, employing a return mapping algorithm for plasticity and a novel hybrid local-nonlocal solver accelerated by octree spatial partitioning for damage evolution. Three benchmark boundary value problems, evaluated across diverse element sizes, characteristic lengths, and softening laws, validate the robustness of the method. The results demonstrate that the proposed nonlocal method achieves mesh- and length-invariant load–displacement responses while accommodating arbitrary softening functions. The presented nonlocal computation method also shows a remarkable computational efficiency compared to the traditional nonlocal computation method.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"188 ","pages":"Article 107587"},"PeriodicalIF":6.2,"publicationDate":"2025-08-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144919812","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A refined incremental homogenization approach for heterogeneous nonlinear geomaterials without isotropization","authors":"Yongdong Li ,&nbsp;Yuanming Lai ,&nbsp;Zhaomin Lv ,&nbsp;Xinwen Wang ,&nbsp;Yian Zhou ,&nbsp;Zipeng Qin ,&nbsp;Lingfeng Guo ,&nbsp;Lunyang Zhao","doi":"10.1016/j.compgeo.2025.107589","DOIUrl":"10.1016/j.compgeo.2025.107589","url":null,"abstract":"<div><div>This study presents a refined incremental homogenization (RIH) framework for predicting the macroscopic nonlinear mechanical behavior of heterogeneous geomaterials composed of inelastic constituents. The proposed approach builds upon the classical Hill incremental (CHI) method while incorporating the reference elastic medium (REM) concept to enhance computational efficiency. A key advantage of the RIH approach is its ability to bypass isotropization procedures, thereby providing a more consistent and physically representative homogenization scheme for inelastic geomaterials. The approach is applied to various geomaterials characterized by a continuous matrix (either elastoplastic or elastoplastic porous) embedded with spherical inclusions (elastic or elastoplastic-damage) or pores. The predictive capability of the RIH approach is assessed through comparisons with full-field finite element (FE) simulations under conventional triaxial compression loading conditions, demonstrating its accuracy in capturing the nonlinear mechanical response. Furthermore, the model’s effectiveness is validated by comparing its predictions with experimental data on Callovo-Oxfordian (COx) argillite with varying mineral compositions. These results highlight the robustness and efficiency of the proposed homogenization framework, making it a valuable tool for modeling the macroscopic behavior of complex geomaterials with inelastic heterogeneities.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"188 ","pages":"Article 107589"},"PeriodicalIF":6.2,"publicationDate":"2025-08-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144916768","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"An Improved Fractional Damage Creep Model for Shale and Its Application to the Long-term Fracture Conductivity Prediction","authors":"Jiale Fu ,&nbsp;Shu Jiang ,&nbsp;Linhao Zhang ,&nbsp;Pengfei Xiong ,&nbsp;Kai Zhang","doi":"10.1016/j.compgeo.2025.107601","DOIUrl":"10.1016/j.compgeo.2025.107601","url":null,"abstract":"<div><div>Formation creep not only reduces the effective gas flow channels but may also lead to engineering disasters, such as wellbore narrowing and collapse. Therefore, it is crucial to conduct in-depth research on the viscoelastic-plastic behavior of shale under complex conditions and its impact on the long-term fracture conductivity. This study combines statistical damage mechanics and fractional calculus theory to establish a shale fractional damage creep model (FDC model) that considers the combined effects of stress, temperature, and time. FDC model is applied to evaluate the fracture aperture, and further, a long-term fracture conductivity model (LFC model) is proposed, which accounts for proppant deformation, embedment, and rock creep. The accuracy and rationality of proposed FDC and LFC model are evaluated using experimental data published in literature. The results indicate that higher external stress accelerates the transition into the tertiary creep stage, leading to the pronounced strain growth and significant reductions in fracture conductivity. Compared to the initial condition, the conductivity decreases by two orders of magnitude after 8.50 d under closure pressure. Moreover, fracture conductivity is positively correlated with proppant particle size and negatively correlated with closure pressure. Both FDC model for characterizing viscoelastic-plastic behavior and LFC model for predicting long-term fracture conductivity demonstrate good applicability and accuracy. The study suggests that both excessively high and low proppant elastic moduli are unfavorable for maintaining stable fracture aperture. This is because, although increasing the proppant elastic modulus reduces its deformation, it simultaneously exacerbates its embedment in formation. In this study, when rock elastic modulus is 13,000 MPa, proppant with an elastic modulus ranging between 80,000 MPa and 110,000 MPa are found to be the most suitable. This research provides new insights for optimizing fracturing design and improving shale gas long-term recovery.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"188 ","pages":"Article 107601"},"PeriodicalIF":6.2,"publicationDate":"2025-08-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144912824","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Corrigendum to “Numerical investigation of phase relationship between kinematic and inertial loads of piles behind a quay wall in liquefiable ground” [Comput. Geotech. 170 (2024) 106265]","authors":"Chunhui Liu ,&nbsp;Jianliao Weng ,&nbsp;Lijun Deng ,&nbsp;Yongzhi Wang ,&nbsp;Qinghe Fang","doi":"10.1016/j.compgeo.2025.107590","DOIUrl":"10.1016/j.compgeo.2025.107590","url":null,"abstract":"","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"188 ","pages":"Article 107590"},"PeriodicalIF":6.2,"publicationDate":"2025-08-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145048773","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A novel experimental–numerical approach to quantitatively evaluate open-ended piles in unsaturated ground","authors":"Junnan Ma ,&nbsp;Xi Xiong ,&nbsp;Jianfeng Xue ,&nbsp;Feng Zhang","doi":"10.1016/j.compgeo.2025.107581","DOIUrl":"10.1016/j.compgeo.2025.107581","url":null,"abstract":"<div><div>The groundwater level (GWL) fluctuates significantly with seasons, which may alter the suction distribution within the ground. Studies have shown that the variation of soil suction can affect the mechanical behavior of the piles. However, few pile load tests have been conducted under controlled suction distribution conditions. Thus, there is a lack of quality data to numerically quantify the influence of suction on the bearing capacity of piles in unsaturated ground. In this study, we propose a novel experimental–numerical approach to quantitatively examine the effects of suction on open-ended single piles in unsaturated ground. A novel model test device was first designed in which the suction distribution of the soil was uniformly controlled using hanging water column technique with a microporous membrane filter. Using this device, a series of load tests were conducted on a model open-ended pile in sandy ground under three specific suction conditions. A soil–water–air coupled finite-element-finite-difference analysis method was used to simulate the model test results. The test and simulation results showed that the suction generated by the changes in GWL increases the Bishop-type skeleton stress and dilatancy of the soil, thereby strengthening the bearing capacity of the pile. Additionally, the numerical method accurately captured the bearing characteristics of the pile in unsaturated ground, making it possible for engineers to efficiently estimate the bearing capacity of the pile in real unsaturated ground, which still remains a puzzle problem and is being avoided constantly till now.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"188 ","pages":"Article 107581"},"PeriodicalIF":6.2,"publicationDate":"2025-08-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144907509","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}