Shuang Cindy Cao , Mengzhen Cao , Pengpeng Zhang , Jongwon Jung , Xiaoshuang Li
{"title":"Impact of fines migration and clogging on gas production in hydrate-bearing sediments: A 2D microfluidic pore model and 3D DEM-CFD study","authors":"Shuang Cindy Cao , Mengzhen Cao , Pengpeng Zhang , Jongwon Jung , Xiaoshuang Li","doi":"10.1016/j.compgeo.2025.107685","DOIUrl":"10.1016/j.compgeo.2025.107685","url":null,"abstract":"<div><div>Natural gas hydrate, as a promising unconventional energy source, often experiences the migration and accumulation of fine particles during the gas extraction process, leading to pore clogging and permeability decline. These effects reduce gas production efficiency and compromise reservoir stability. However, a systematic understanding of the microscopic mechanisms governing fine particles migration remains lacking. This study investigates the influence of fine type, pore-throat size, pore fluid chemistry, as well as hydrate saturation, and depressurization conditions on fines migration, clogging, and permeability evolution using both two-dimensional microfluidic pore model experiments and three-dimensional DEM-CFD coupled numerical simulations. The experimental results demonstrate that the electrical sensitivity of fines and pore-throat geometry jointly determine the critical pore-clogging concentration. Cohesive representative fines—montmorillonite—tends to form flocculation structures in saltwater, while non-cohesive fines—silica particles—exhibit high mobility. Numerical simulations further reveal that a higher depressurization gradient and lower saturation increase the fines migration rate and permeability while slowing the clogging evolution process during hydrate dissociation. The findings indicate that fines migration behavior has strong multi-physical field coupling characteristics. The proposed microfluidic pore model experiments and DEM-CFD modeling effectively elucidate microscopic mechanisms and quantify macroscopic effects, providing theoretical support for optimizing hydrate extraction strategies and ensuring reservoir stability.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"190 ","pages":"Article 107685"},"PeriodicalIF":6.2,"publicationDate":"2025-10-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145247975","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":"Thermal fracture modelling of borehole-heated rocks: A focus on inhomogeneities, physical properties, and heating rates","authors":"Zhengkuo Ma , Chunshun Zhang , Congying Li","doi":"10.1016/j.compgeo.2025.107670","DOIUrl":"10.1016/j.compgeo.2025.107670","url":null,"abstract":"<div><div>The thermal cracking characteristics of borehole-heated rocks are significantly affected by inhomogeneities, physical properties, and heating rates. While most current numerical methods using Weibull-based models account for material inhomogeneity, they often overlook the simulation of real cracks. Although the finite discrete element method (FDEM) is effective for fracture simulation, this study enhances the process by integrating Voronoi diagrams and Weibull distributions to model rock inhomogeneity. Additionally, the finite discrete element method coupled with thermo-mechanical damage (FDEM-TMD) is employed to simulate realistic thermal cracking processes. This comprehensive model considers the distribution of temperature and thermal stress fields throughout the entire process from crack initiation to propagation. The model’s validity is confirmed through analytical solutions, experimental data, and other numerical simulations, which collectively demonstrate the influence of inhomogeneities, physical properties, and heating rates on rock fracturing. The results show that varying degrees of inhomogeneity can affect rock rupture behavior, while physical properties can impact energy consumption and heat transfer, ultimately influencing fracture efficiency. Furthermore, rapid heating may prolong stress accumulation and increase the probability of cracking. These insights enhance our understanding of heating-induced fracturing mechanisms and offer valuable theoretical guidance for applications in nuclear waste disposal, geothermal energy extraction, and laser-assisted rock fragmentation.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107670"},"PeriodicalIF":6.2,"publicationDate":"2025-10-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145220999","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}
Liting Zhang , Hang Zhou , Jianxin Wang , Chunyong Jiang
{"title":"Numerical investigation of advancement ratio and helix plate number on stress response and soil disturbance during helical pile installation in sand","authors":"Liting Zhang , Hang Zhou , Jianxin Wang , Chunyong Jiang","doi":"10.1016/j.compgeo.2025.107676","DOIUrl":"10.1016/j.compgeo.2025.107676","url":null,"abstract":"<div><div>The installation process of helical piles induces substantial changes in the surrounding soil, including stress redistribution, void ratio variation, and strain localization, all of which significantly influence pile performance. This study investigates the impact of advancement ratio (AR) and helix plate number on stress responses and soil disturbance during pile installation in dense sand. A density- and stress-dependent hypoplastic constitutive model is implemented within a Coupled Eulerian-Lagrangian (CEL) framework to capture the nonlinear and non-axisymmetric soil behavior. Parametric simulations involving single, double, and triple plate piles with varying AR values reveal that AR strongly governs the extent of stress concentration, dilation-induced loosening, and effective energy transfer into deep soil layers. Higher AR values intensify vertical compression, while multiple helix plates promote uniform compaction and reduce stress anisotropy. These asymmetric stress patterns shift stress centroids and generate directional stress concentration zones, which are critical for understanding helical pile-soil interaction.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107676"},"PeriodicalIF":6.2,"publicationDate":"2025-10-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145221000","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}
Wenjun Lu , Dongting Cai , Jinhui Li , Aiting Wang
{"title":"The role of shared suction anchors for mitigating cascading failure in floating offshore wind farms","authors":"Wenjun Lu , Dongting Cai , Jinhui Li , Aiting Wang","doi":"10.1016/j.compgeo.2025.107655","DOIUrl":"10.1016/j.compgeo.2025.107655","url":null,"abstract":"<div><div>Shared anchors offer significant potential for reducing mooring costs and enhancing floating wind competitiveness, yet industry adoption remains hindered by cascading failure risks in multi-anchor configurations. Critical research gaps persist regarding the initiation mechanisms, propagation dynamics, and mitigation strategies for such failures in shared-anchor floating wind farms. This study addresses these gaps through a novel systematic methodology that derives shared suction anchor loading conditions using inverted catenary theory. Contrary to conventional assumptions, the analysis reveals that multi-line suction anchors exhibit no mechanical superiority over conventional foundations. Instead, they demonstrate less rational force distribution and greater failure susceptibility under environmental loading. Following initial anchor failure, platform displacement in shared-anchor wind farm decreases by 80% relative to conventional floating wind farms using strengthened chains, but the peak tension rises to 1.25 times owing to the dynamic amplification effect. Although variations in soil capacity affect mooring chain inclination angles, 3-line anchors remain most prone to initial failure across all soil conditions. After the initial failure of adjacent anchors, the susceptibility to secondary failure of 3-line anchors markedly decreases in stiffer soils. To prevent cascading failures, suction anchor designs must satisfy dual safety criteria: withstand intact-system loads and accommodate post-initial-failure amplified secondary anchor loads.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107655"},"PeriodicalIF":6.2,"publicationDate":"2025-10-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145221001","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}
Haoting Li , Yuting Zhao , Jipeng Li , Haoyang Yin , Cen Song , Min Hao
{"title":"CFD-DEM investigation on continuous invasion behaviors of particles and evolution of pore clogging in porous media","authors":"Haoting Li , Yuting Zhao , Jipeng Li , Haoyang Yin , Cen Song , Min Hao","doi":"10.1016/j.compgeo.2025.107675","DOIUrl":"10.1016/j.compgeo.2025.107675","url":null,"abstract":"<div><div>The phenomenon of particle migration and clogging in porous media is prevalent in geotechnical and petroleum engineering, particularly when involving continuous particle invasion and agglomeration processes that necessitate in-depth investigation. In this work, the two-way coupling Computational Fluid Dynamics-Discrete Element Method (CFD-DEM) is employed to investigate the microscopic dynamic behaviors of particle continuous invasion and clogging by multiple invasion cycles at different fluid viscosities and friction coefficients of particles. The evolution of agglomerate formation and pore clogging are predicted through metrics including aggregation growth rate and development curves of penetrating particles. Results reveal that the particle clogging mechanisms transitions from point clogging within individual pores to surface clogging across multiple pores with the increase of invasion cycle. The retention of fine particles showed a continuous increment with an accumulation of over 80 % to 95 %. The severe damage after the second invasion cycle is predicted according to the decrease of normalized permeability to 0.33. The average translational velocity of particles increases together with fluid viscosity, thus reducing the particle penetration time and particle retention rate. At low fluid viscosity, particles exhibit a multi-point dispersed retention pattern, which promotes the formation of particle aggregates. Besides, the high-friction-coefficient particles experience stronger static friction force and sliding resistance at the pore throat, inhibiting their continued migration with fluid, which also decreases the average invasion depth and increases the particle retention rate. Decreasing the friction coefficients will improve particle aggregate growth rates and promote the formation of agglomeration.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107675"},"PeriodicalIF":6.2,"publicationDate":"2025-10-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145220997","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}
Xin-Dong Wei , Muhammad Shoaib , Tao Zhou , Yuhang Wu , Zhe Li , Yang Zou , Gao-Feng Zhao
{"title":"Numerical simulation of rock creep and failure using a Norton-Bailey-based four-dimensional lattice spring model","authors":"Xin-Dong Wei , Muhammad Shoaib , Tao Zhou , Yuhang Wu , Zhe Li , Yang Zou , Gao-Feng Zhao","doi":"10.1016/j.compgeo.2025.107669","DOIUrl":"10.1016/j.compgeo.2025.107669","url":null,"abstract":"<div><div>Creep deformation significantly influences the long-term stability of rock structures subjected to sustained loading. This study presents a novel computational framework that integrates the Zienkiewicz-Pande (ZP) plasticity model and the Norton-Bailey (NB) creep law within a four-dimensional lattice spring model (4D-LSM), enabling simulation of the coupled elastic, plastic, and creep responses of geomaterials. To capture tertiary creep characterized by an accelerated deformation rate, the NB model is extended with a strain-based elasto-brittle failure criterion, leveraging the inherent fracture-handling capabilities of the 4D-LSM. The proposed framework is validated against experimental and analytical benchmarks through single-stage uniaxial and multi-stage triaxial compression simulations, demonstrating good agreement and predictive accuracy. The model’s applicability is further demonstrated in tunnel excavation analyses, showcasing its effectiveness in predicting long-term time-dependent deformation. A three-point bending test is also simulated to assess the model’s capability in reproducing creep-induced fracture. Overall, the ZP-NB-4D-LSM framework offers a robust and efficient tool for simulating the full creep-fracture process in rock engineering applications.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107669"},"PeriodicalIF":6.2,"publicationDate":"2025-10-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145221002","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}
Boris Kratz , Pierre Jehel , Maxime Tatin , Emmanuel Vazquez
{"title":"Parameter influence analysis in a 3D TBM model via sensitivity analysis and GP metamodels","authors":"Boris Kratz , Pierre Jehel , Maxime Tatin , Emmanuel Vazquez","doi":"10.1016/j.compgeo.2025.107650","DOIUrl":"10.1016/j.compgeo.2025.107650","url":null,"abstract":"<div><div>Urban tunnel excavation with tunnel boring machines induces ground movements that can affect nearby structures. Three-dimensional finite element models (FEM) are widely used to predict these settlements, but their high computational cost limits direct exploration of parameter influence. This work presents a 3D FEM simulator of mechanized tunneling and a methodology to quantify the impact of both numerical and physical inputs on settlement predictions. First, an <em>accuracy-cost model reduction</em> study evaluates the effect of domain dimensions and mesh densities on a small number of scalar quantities of interest extracted from simulated settlement fields. Empirical error models are fitted and used to select a reduced configuration that balances accuracy and runtime. Second, Gaussian process models are trained on simulation data from the reduced configuration and validated using exact leave-one-out cross-validation. These metamodels enable the computation of Sobol’ sensitivity indices with quantified uncertainty, identifying the most influential geological, operational, and loading parameters. The proposed framework reduces the cost of sensitivity analysis for computationally intensive 3D tunneling simulations, supporting input screening and dimensionality reduction for design and calibration.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107650"},"PeriodicalIF":6.2,"publicationDate":"2025-10-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145221004","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":"Collapse mechanism of stratum induced by subsurface excavation for a shallow buried tunnel in water-rich region","authors":"Fu Huang, Yongtao Wang, Min Zhang, Lu Chen","doi":"10.1016/j.compgeo.2025.107672","DOIUrl":"10.1016/j.compgeo.2025.107672","url":null,"abstract":"<div><div>The collapse of stratum induced by subsurface excavation for a shallow tunnel in water-rich region is common in urban tunnel engineering. Although scholars have conducted numerous research on the stability of the soil mass above shallow-buried tunnel roof, the collapse characteristics of the soil mass induced by shallow-buried tunnel construction in water-rich strata have not been well investigated. In this work, a new three-dimensional collapse mechanism of the stratum above shallow-buried tunnel roof is generated “point by point” based on the spatial discretization technique for the first time. A scaled model test is designed to verify the validity of the constructed collapse mechanism. By introducing the pore water pressure which is obtained from numerical simulation as an external force into the virtual work equation, the safety factor of the stratum is derived from the upper bound theorem of limit analysis. Finally, a comparison between the theoretical solutions and numerical solutions have been conducted to prove the effectiveness of the proposed method</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107672"},"PeriodicalIF":6.2,"publicationDate":"2025-10-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145220998","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}
He-mei Sun , Xin Kang , Xiong-ying Ma , Shu-ying Wang , Ren-peng Chen
{"title":"A Mesoscale Perspective on Micron-sized Clay Aggregates Compression: Insight from Coarse-Grained Molecular Dynamics","authors":"He-mei Sun , Xin Kang , Xiong-ying Ma , Shu-ying Wang , Ren-peng Chen","doi":"10.1016/j.compgeo.2025.107668","DOIUrl":"10.1016/j.compgeo.2025.107668","url":null,"abstract":"<div><div>In this paper, coarse-grained molecular dynamics (CGMD) employing the Gay-Berne (GB) potential was applied to investigate the mesoscale compression behavior of micron-sized clay aggregates fulfilled with explict pore water fluids. A systematic set of CGMD-compatible methods was proposed to characterize the mesoscale structure of clay. Simulation results revealed a topology transformation from loose “honeycomb” structure, featuring homogeneous clay platelets orientations, diverse pore morphologies, and normally distributed pore size, into dense “band type” structure, characterized by uniform clay platelet orientation, flattened pore morphology, and micron-sized pores. The pore water migration showed a three-stage progression. Initially, pore water fluids fully or partially filled the clay pores of various shapes. Subsequently, pores flattend gradually while pore water molecules adhered onto clay platelets. Finally, clay platelets “lay down” to form a flat micropore full of water. Meanwhile, the clay skeleton transitioned from stability to collapse and ultimately to failure with the increase in pore pressure heterogeneity. With axial strain increasing, the clay aggregates initially experienced elastic deformation, followed by plastic deformation and lattice fracture. These findings suggest that the microscopic failure mechanism of clayey soils under external stress involves the development of heterogeneous microstructural alterations and non-uniform pore water pressure gradients. In conclusion, this research offers novel insights into the micro-mechanisms controlling the macroscale compression behaviors of clayey soils and highlights CGMD technique as a powerful tool for investigating microscopic response of clay under different conditions, including high geostress, high hydrostatic pressure, freeze–thaw cycles, and cyclic loading–unloading.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107668"},"PeriodicalIF":6.2,"publicationDate":"2025-09-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145220995","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":"Influence of gradation on shear behavior and pore water response of sand","authors":"Mohammed Elnur, Khalid A. Alshibli","doi":"10.1016/j.compgeo.2025.107674","DOIUrl":"10.1016/j.compgeo.2025.107674","url":null,"abstract":"<div><div>This paper examines the influence of particle size distribution (PSD) on the shearing behavior of sand using axisymmetric triaxial compression (ATC) testing in conjunction with 3D imaging using synchrotron micro-computed tomography (SMT). Angular sand was sieved to obtain a narrow gradation (uniform) and a wide gradation (non-uniform) with uniformity coefficients (<span><math><mrow><msub><mi>C</mi><mi>u</mi></msub></mrow></math></span>) of 1.2 and 3.0, respectively. Experiments were conducted under drained conditions and a constant effective cell pressure (<span><math><mrow><msubsup><mi>σ</mi><mrow><mn>3</mn></mrow><mo>′</mo></msubsup></mrow></math></span>) of 50 kPa at low and high back pressures (BP). Results of paired ATC-SMT were analyzed in terms of global constitutive behavior, strain localization, and localized flow characteristics. A new particle tracking algorithm is developed to expose strain localization. Relative particle translation gradient (RPTG) results revealed that gradation affects how early the shear band develops and a tendency for narrow gradation to exhibit varying degrees of diffused shear band at high axial strain (<span><math><mrow><msub><mi>ε</mi><mn>1</mn></msub></mrow></math></span>) even though a clear single shear band was observed at the peak strength. Geometric tortuosity (<span><math><mrow><mi>τ</mi></mrow></math></span>) is used as a quantitative parameter to assess flow changes and the efficiency of the flow path. Representative elementary volumes (REVs) inside and outside the shear band were analyzed and the results showed a variation in flow path characteristics where wide gradation specimens had preferential flow paths.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107674"},"PeriodicalIF":6.2,"publicationDate":"2025-09-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145220883","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}