Computers and Geotechnics最新文献_第2页

The performance of quadratic finite-discrete element method (qFDEM) and its potential advantages 二次有限离散元法（qFDEM）的性能及其潜在优势

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-01-28 DOI: 10.1016/j.compgeo.2026.107925

Zhonghao Li , Xiaofeng Li , Haibo Li , Qi Zhao , Giovanni Grasselli

The combined finite–discrete element method (FDEM) has proven effective for simulating crack initiation, propagation, and coalescence in brittle solids. However, existing FDEM frameworks remain limited to constant-strain elements, leading to restricted capability in representing complex stress fields, pronounced sensitivity to shear and volumetric locking, and a strong tendency toward numerical dispersion in dynamic problems. To overcome these limitations, this study develops a high order element-based framework incorporating a novel quadratic cohesive element to enhance model accuracy and continuity. The proposed quadratic cohesive element ensures uniform traction distribution along edges, avoiding the mid-node stress concentrations that typically lead to mesh incompatibility and artificial strength reduction. Three quasi-static loading tests and one wave propagation test are performed to compare quadratic and linear models. The results show that the quadratic model consistently outperforms the linear counterpart in stress path, crack propagation, and mitigating numerical dispersion. In quasi-static loading, the new quadratic model exhibits a lower error in stress, predicts a more precise crack initiation load, and provides more reliable crack path predictions compared with previous models. In dynamic conditions, it can effectively mitigate the numerical dispersion of high-frequency wave components that low-order elements struggle with and provide more stable wave propagation simulations. Moreover, the quadratic elements FDEM framework offers an economical alternative for enhancing the fidelity of crack simulations: compared to mesh refinement, quadratic elements achieve comparable accuracy in crack initiation load prediction with only 50–60% of the computational cost.

有限-离散元组合方法（FDEM）已被证明是模拟脆性固体裂纹萌生、扩展和合并的有效方法。然而，现有的FDEM框架仍然局限于恒应变单元，导致表征复杂应力场的能力有限，对剪切和体积锁定的敏感性明显，并且在动态问题中有很强的数值离散倾向。为了克服这些限制，本研究开发了一个基于高阶元素的框架，其中包含了一个新的二次内聚元素，以提高模型的准确性和连续性。所提出的二次内聚单元确保沿边缘均匀的牵引力分布，避免了通常导致网格不兼容和人工强度降低的中间节点应力集中。进行了3次准静态加载试验和1次波传播试验，对二次模型和线性模型进行了比较。结果表明，二次型模型在应力路径、裂纹扩展和减轻数值弥散方面均优于线性模型。在准静态加载下，新二次元模型的应力误差更小，预测裂纹起裂载荷更精确，裂纹路径预测更可靠。在动态条件下，它可以有效地缓解低阶元与高频波分量的数值频散，提供更稳定的波传播模拟。此外，二次元FDEM框架为提高裂纹模拟的保真度提供了一种经济的替代方案：与网格细化相比，二次元在裂纹起裂载荷预测方面达到了相当的精度，而计算成本仅为50-60%。

{"title":"The performance of quadratic finite-discrete element method (qFDEM) and its potential advantages","authors":"Zhonghao Li , Xiaofeng Li , Haibo Li , Qi Zhao , Giovanni Grasselli","doi":"10.1016/j.compgeo.2026.107925","DOIUrl":"10.1016/j.compgeo.2026.107925","url":null,"abstract":"<div><div>The combined finite–discrete element method (FDEM) has proven effective for simulating crack initiation, propagation, and coalescence in brittle solids. However, existing FDEM frameworks remain limited to constant-strain elements, leading to restricted capability in representing complex stress fields, pronounced sensitivity to shear and volumetric locking, and a strong tendency toward numerical dispersion in dynamic problems. To overcome these limitations, this study develops a high order element-based framework incorporating a novel quadratic cohesive element to enhance model accuracy and continuity. The proposed quadratic cohesive element ensures uniform traction distribution along edges, avoiding the mid-node stress concentrations that typically lead to mesh incompatibility and artificial strength reduction. Three quasi-static loading tests and one wave propagation test are performed to compare quadratic and linear models. The results show that the quadratic model consistently outperforms the linear counterpart in stress path, crack propagation, and mitigating numerical dispersion. In quasi-static loading, the new quadratic model exhibits a lower error in stress, predicts a more precise crack initiation load, and provides more reliable crack path predictions compared with previous models. In dynamic conditions, it can effectively mitigate the numerical dispersion of high-frequency wave components that low-order elements struggle with and provide more stable wave propagation simulations. Moreover, the quadratic elements FDEM framework offers an economical alternative for enhancing the fidelity of crack simulations: compared to mesh refinement, quadratic elements achieve comparable accuracy in crack initiation load prediction with only 50–60% of the computational cost.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107925"},"PeriodicalIF":6.2,"publicationDate":"2026-01-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080670","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}

引用次数: 0

On the integration of an ACST-based bounding surface model 基于acst的边界面模型的集成

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-01-27 DOI: 10.1016/j.compgeo.2026.107911

Srinivas Vivek Bokkisa , Jorge Macedo , Pedro Arduino

Anisotropic critical state theory (ACST) provides a framework for incorporating fabric effects in constitutive models. However, most previous efforts have focused on constitutive aspects with comparatively limited attention to numerical implementations. This study presents a comprehensive assessment of explicit and implicit implementations of the ACST-based bounding surface model, SANISAND-F. Assessments are conducted in terms of stability, accuracy, computational efficiency, and both local and global performance.

In the implicit implementation, the critical importance of accurate gradient calculations is highlighted, introducing a verification procedure that enables quadratic convergence. The explicit and implicit implementations exhibit stability, producing smooth and bounded responses across a wide range of strain increments and numerical tolerances. However, their accuracy differs significantly. The implicit implementation is sensitive to the initial loading state, strain increment, and loading direction, showing minor dependence on the solver tolerance. In contrast, the explicit implementation is influenced by both strain increment and substepping tolerance, and at practical tolerance and strain increment levels, it often outperforms the implicit scheme in accuracy. Regarding efficiency, the explicit implementation proves more efficient at the local integration level. However, at the global level, the implicit implementation with the consistent tangent exhibits a faster rate of convergence in global equilibrium iterations. Nonetheless, the overall computational cost at the global level is not definitive when comparing explicit and implicit schemes; it varies with simulations and loading-specific parameters, as demonstrated through the included boundary-value problems.

各向异性临界状态理论（ACST）为将织物效应纳入本构模型提供了一个框架。然而，大多数以前的努力都集中在本构方面，对数值实现的关注相对有限。本研究对基于acst的边界面模型SANISAND-F的显式和隐式实现进行了全面评估。评估是根据稳定性、准确性、计算效率以及本地和全局性能进行的。在隐式实现中，强调了精确梯度计算的关键重要性，引入了一个验证程序，使二次收敛成为可能。显式和隐式实现表现出稳定性，在广泛的应变增量和数值公差范围内产生平滑和有界的响应。然而，它们的准确性差别很大。隐式实现对初始加载状态、应变增量和加载方向敏感，对求解器公差的依赖较小。相比之下，显式实现受应变增量和步进公差的影响，在实际公差和应变增量水平上，显式实现的精度往往优于隐式方案。关于效率，显式实现证明在本地集成级别更有效。然而，在全局水平上，具有一致切线的隐式实现在全局平衡迭代中表现出更快的收敛速度。尽管如此，当比较显式和隐式方案时，全球一级的总体计算成本并不确定；它随模拟和加载特定参数的变化而变化，正如通过所包含的边值问题所证明的那样。

{"title":"On the integration of an ACST-based bounding surface model","authors":"Srinivas Vivek Bokkisa , Jorge Macedo , Pedro Arduino","doi":"10.1016/j.compgeo.2026.107911","DOIUrl":"10.1016/j.compgeo.2026.107911","url":null,"abstract":"<div><div>Anisotropic critical state theory (ACST) provides a framework for incorporating fabric effects in constitutive models. However, most previous efforts have focused on constitutive aspects with comparatively limited attention to numerical implementations. This study presents a comprehensive assessment of explicit and implicit implementations of the ACST-based bounding surface model, SANISAND-F. Assessments are conducted in terms of stability, accuracy, computational efficiency, and both local and global performance.</div><div>In the implicit implementation, the critical importance of accurate gradient calculations is highlighted, introducing a verification procedure that enables quadratic convergence. The explicit and implicit implementations exhibit stability, producing smooth and bounded responses across a wide range of strain increments and numerical tolerances. However, their accuracy differs significantly. The implicit implementation is sensitive to the initial loading state, strain increment, and loading direction, showing minor dependence on the solver tolerance. In contrast, the explicit implementation is influenced by both strain increment and substepping tolerance, and at practical tolerance and strain increment levels, it often outperforms the implicit scheme in accuracy. Regarding efficiency, the explicit implementation proves more efficient at the local integration level. However, at the global level, the implicit implementation with the consistent tangent exhibits a faster rate of convergence in global equilibrium iterations. Nonetheless, the overall computational cost at the global level is not definitive when comparing explicit and implicit schemes; it varies with simulations and loading-specific parameters, as demonstrated through the included boundary-value problems.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107911"},"PeriodicalIF":6.2,"publicationDate":"2026-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080667","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}

引用次数: 0

Extension of explicit Runge-Kutta substepping stress integration for viscoplastic model of saturated soils 饱和土粘塑性模型显式Runge-Kutta分段应力积分的推广

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-01-27 DOI: 10.1016/j.compgeo.2026.107937

Wei CHENG , Zhen-Yu YIN

Stable integration schemes are critically important for rate-dependent constitutive models, serving as a cornerstone for ensuring accuracy, efficiency, and robustness in finite element implementations. This paper investigates the numerical performance of explicit stress integration schemes with adaptive substepping for integrating a newly proposed fractional consistency two-surface viscoplastic model for saturated clays. The incremental stress–strain-strain rate relation of the model can be linearized following the consistency condition of the rate-dependent loading surface and subsequently integrated using four distinct explicit Runge-Kutta substepping integration algorithms (i.e., RK12, RK23, RK34, RK45) with automatic error control and stress drift correction techniques. The overall numerical performance of the algorithms in terms of accuracy and efficiency is evaluated at both the material point level (i.e., isotropic, oedometric, and triaxial compression tests) and the boundary-value problem level (i.e., piezocone penetration and underground gallery excavation), which demonstrates that the RK23 and RK34 algorithms perform excellently in balancing accuracy and computational cost. The proposed algorithms provide a versatile and adaptive framework for integrating time-dependent constitutive equations, particularly those based on the consistency viscoplastic approaches commonly used in advanced rate-dependent modeling, allowing for a wide range of geotechnical engineering applications.

稳定的积分方案对于速率相关的本构模型至关重要，是确保有限元实现中准确性、效率和鲁棒性的基石。本文研究了带自适应分步法的显式应力积分方案对新提出的饱和黏土分数浓度双表面粘塑性模型的数值性能。模型的增量应力-应变-应变率关系可根据速率相关加载面一致性条件进行线性化处理，并采用四种不同的显式龙格-库塔分步积分算法（RK12、RK23、RK34、RK45），结合自动误差控制和应力漂移校正技术进行积分。在材料点水平（即各向同性、测径和三轴压缩试验）和边值问题水平（即压电锥穿透和地下巷道开挖）上对算法在精度和效率方面的整体数值性能进行了评估，表明RK23和RK34算法在平衡精度和计算成本方面表现出色。所提出的算法为积分时变本构方程提供了一个通用的自适应框架，特别是那些基于先进速率相关建模中常用的一致性粘塑性方法的框架，允许广泛的岩土工程应用。

{"title":"Extension of explicit Runge-Kutta substepping stress integration for viscoplastic model of saturated soils","authors":"Wei CHENG , Zhen-Yu YIN","doi":"10.1016/j.compgeo.2026.107937","DOIUrl":"10.1016/j.compgeo.2026.107937","url":null,"abstract":"<div><div>Stable integration schemes are critically important for rate-dependent constitutive models, serving as a cornerstone for ensuring accuracy, efficiency, and robustness in finite element implementations. This paper investigates the numerical performance of explicit stress integration schemes with adaptive substepping for integrating a newly proposed fractional consistency two-surface viscoplastic model for saturated clays. The incremental stress–strain-strain rate relation of the model can be linearized following the consistency condition of the rate-dependent loading surface and subsequently integrated using four distinct explicit Runge-Kutta substepping integration algorithms (i.e., RK12, RK23, RK34, RK45) with automatic error control and stress drift correction techniques. The overall numerical performance of the algorithms in terms of accuracy and efficiency is evaluated at both the material point level (i.e., isotropic, oedometric, and triaxial compression tests) and the boundary-value problem level (i.e., piezocone penetration and underground gallery excavation), which demonstrates that the RK23 and RK34 algorithms perform excellently in balancing accuracy and computational cost. The proposed algorithms provide a versatile and adaptive framework for integrating time-dependent constitutive equations, particularly those based on the consistency viscoplastic approaches commonly used in advanced rate-dependent modeling, allowing for a wide range of geotechnical engineering applications.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107937"},"PeriodicalIF":6.2,"publicationDate":"2026-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080669","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}

引用次数: 0

Relating solution tests to pore-scale CaCO3 crystal growth: Numerical simulation based on the phase field method 孔尺度CaCO3晶体生长的相关溶液试验：基于相场法的数值模拟

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-01-27 DOI: 10.1016/j.compgeo.2026.107947

Guoliang Ma , Yang Xiao , Jinxuan Zhang , Zhichao Zhang

Microbially induced carbonate precipitation (MICP) has been widely studied for its broad application potential in geotechnical engineering. Exploring the pore-scale CaCO₃ precipitation processes is one of the key aspects of fully understanding the mechanism of MICP for enhancing mechanical properties and controlling hydraulic properties of soils and rocks. Additionally, solution tests are typically the first step to understanding the reaction processes of MICP before conducting pore-scale and macro-scale tests. However, most explanations remain qualitative and do not establish a direct link between solution tests, pore-scale processes, and macro-scale performance. In this study, a saturation-dependent kinetic model is developed to describe the chemical reactions involved in MICP, and a phase-field model is employed to simulate the growth of CaCO₃ crystals at the pore scale. The two models are coupled through the CaCO₃ precipitation rate, enabling direct translation of solution test parameters into pore-scale simulations. Model predictions are validated against experimental data from microdroplet tests, showing good agreement in crystal growth dynamics under varying bacterial densities. The results demonstrate the capability of the proposed approach to quantitatively link solution chemistry with pore-scale mineralization processes. This work provides a robust framework for further investigation of MICP-induced pore sealing and particle bridging, contributing to the optimization and design of bio-mediated ground improvement strategies.

微生物诱导碳酸盐沉淀（MICP）因其在岩土工程中的广泛应用潜力而受到广泛的研究。探索孔隙尺度CaCO3的沉淀过程是充分认识MICP增强岩土力学性能和控制水力性能机理的关键方面之一。此外，在进行孔隙尺度和宏观尺度测试之前，溶液测试通常是了解MICP反应过程的第一步。然而，大多数解释仍然是定性的，并没有在溶液测试、孔隙尺度过程和宏观尺度性能之间建立直接联系。在本研究中，建立了一个饱和度相关的动力学模型来描述MICP中涉及的化学反应，并采用相场模型来模拟CaCO3晶体在孔隙尺度上的生长。这两个模型通过CaCO3沉淀速率耦合，可以将溶液测试参数直接转换为孔隙尺度模拟。模型预测与微滴测试的实验数据相对照，显示出在不同细菌密度下晶体生长动力学的良好一致性。结果表明，该方法能够定量地将溶液化学与孔隙尺度的矿化过程联系起来。这项工作为进一步研究micp诱导的孔隙密封和颗粒桥接提供了一个强大的框架，有助于优化和设计生物介导的地面改善策略。

{"title":"Relating solution tests to pore-scale CaCO3 crystal growth: Numerical simulation based on the phase field method","authors":"Guoliang Ma , Yang Xiao , Jinxuan Zhang , Zhichao Zhang","doi":"10.1016/j.compgeo.2026.107947","DOIUrl":"10.1016/j.compgeo.2026.107947","url":null,"abstract":"<div><div>Microbially induced carbonate precipitation (MICP) has been widely studied for its broad application potential in geotechnical engineering. Exploring the pore-scale CaCO<sub>3</sub> precipitation processes is one of the key aspects of fully understanding the mechanism of MICP for enhancing mechanical properties and controlling hydraulic properties of soils and rocks. Additionally, solution tests are typically the first step to understanding the reaction processes of MICP before conducting pore-scale and macro-scale tests. However, most explanations remain qualitative and do not establish a direct link between solution tests, pore-scale processes, and macro-scale performance. In this study, a saturation-dependent kinetic model is developed to describe the chemical reactions involved in MICP, and a phase-field model is employed to simulate the growth of CaCO<sub>3</sub> crystals at the pore scale. The two models are coupled through the CaCO<sub>3</sub> precipitation rate, enabling direct translation of solution test parameters into pore-scale simulations. Model predictions are validated against experimental data from microdroplet tests, showing good agreement in crystal growth dynamics under varying bacterial densities. The results demonstrate the capability of the proposed approach to quantitatively link solution chemistry with pore-scale mineralization processes. This work provides a robust framework for further investigation of MICP-induced pore sealing and particle bridging, contributing to the optimization and design of bio-mediated ground improvement strategies.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107947"},"PeriodicalIF":6.2,"publicationDate":"2026-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080673","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}

引用次数: 0

A graph-based algorithm for the Continuous-Projection Embedded Discrete Fracture Model (CpEDFM-U) to simulate two-phase flows in naturally fractured porous media using the MPFA-D method on general unstructured meshes 基于图的连续投影嵌入式离散裂缝模型（CpEDFM-U）在一般非结构化网格上采用MPFA-D方法模拟天然裂缝多孔介质中的两相流动

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-01-27 DOI: 10.1016/j.compgeo.2026.107918

Gabriel Martins Cavalcanti Feitosa , Emanoel Rodrigues dos Santos , Pedro Victor Paixão Albuquerque , Artur Castiel Reis de Souza , Darlan Karlo Elisiário de Carvalho , Paulo Roberto Maciel Lyra

Modeling fluid flow in naturally fractured porous media is crucial for applications such as hydrocarbon production and CO2 sequestration. However, accurately simulating these flows remains challenging due to fractures with complex permeability distributions. The Embedded Discrete Fracture Model (EDFM) has been widely used but has limitations in representing fractures acting as barriers, especially in multiphase flows. To overcome these challenges, the Projection-based Embedded Discrete Fracture Model (pEDFM) was developed, offering better handling of fractures with lower permeability than the matrix. However, it can still exhibit discontinuities in fracture projections. To address these limitations, we propose the Continuous-Projection Embedded Discrete Fracture Model (CpEDFM-U), a graph-based algorithm that guarantees continuous fracture projections in both structured and unstructured 2D meshes. The CpEDFM-U uses Dijkstra’s algorithm to find the shortest path between fracture tips and applies the MPFA-D method for matrix flow and TPFA for fractures. In numerical simulations analyzed, CpEDFM-U outperforms EDFM and pEDFM, demonstrating lower errors and robust performance across different fracture types and mesh resolutions.

模拟天然裂缝多孔介质中的流体流动对于油气生产和二氧化碳封存等应用至关重要。然而，由于具有复杂渗透率分布的裂缝，准确模拟这些流动仍然具有挑战性。嵌入式离散裂缝模型（EDFM）已被广泛应用，但在描述裂缝作为屏障时存在局限性，特别是在多相流中。为了克服这些挑战，开发了基于投影的嵌入式离散裂缝模型（pEDFM），可以更好地处理渗透率较低的裂缝。然而，它仍然可以在裂缝突出处表现出不连续性。为了解决这些限制，我们提出了连续投影嵌入式离散裂缝模型（CpEDFM-U），这是一种基于图的算法，可以保证在结构化和非结构化二维网格中连续的裂缝投影。CpEDFM-U使用Dijkstra算法寻找裂缝尖端之间的最短路径，并对基质流动采用MPFA-D方法，对裂缝采用TPFA方法。在数值模拟分析中，CpEDFM-U优于EDFM和pEDFM，在不同裂缝类型和网格分辨率下显示出更低的误差和更强的性能。

{"title":"A graph-based algorithm for the Continuous-Projection Embedded Discrete Fracture Model (CpEDFM-U) to simulate two-phase flows in naturally fractured porous media using the MPFA-D method on general unstructured meshes","authors":"Gabriel Martins Cavalcanti Feitosa , Emanoel Rodrigues dos Santos , Pedro Victor Paixão Albuquerque , Artur Castiel Reis de Souza , Darlan Karlo Elisiário de Carvalho , Paulo Roberto Maciel Lyra","doi":"10.1016/j.compgeo.2026.107918","DOIUrl":"10.1016/j.compgeo.2026.107918","url":null,"abstract":"<div><div>Modeling fluid flow in naturally fractured porous media is crucial for applications such as hydrocarbon production and CO2 sequestration. However, accurately simulating these flows remains challenging due to fractures with complex permeability distributions. The Embedded Discrete Fracture Model (EDFM) has been widely used but has limitations in representing fractures acting as barriers, especially in multiphase flows. To overcome these challenges, the Projection-based Embedded Discrete Fracture Model (pEDFM) was developed, offering better handling of fractures with lower permeability than the matrix. However, it can still exhibit discontinuities in fracture projections. To address these limitations, we propose the Continuous-Projection Embedded Discrete Fracture Model (CpEDFM-U), a graph-based algorithm that guarantees continuous fracture projections in both structured and unstructured 2D meshes. The CpEDFM-U uses Dijkstra’s algorithm to find the shortest path between fracture tips and applies the MPFA-D method for matrix flow and TPFA for fractures. In numerical simulations analyzed, CpEDFM-U outperforms EDFM and pEDFM, demonstrating lower errors and robust performance across different fracture types and mesh resolutions.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107918"},"PeriodicalIF":6.2,"publicationDate":"2026-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080672","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}

引用次数: 0

Coupled mechanisms of load transfer and fracture in bonded granular materials using CT image-based discrete element method 基于CT图像离散元法的粘结颗粒材料载荷传递与断裂耦合机制研究

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-01-23 DOI: 10.1016/j.compgeo.2025.107891

Hu Yang , Yiik Diew Wong , Liyan Shan , Lingwen Li

Bonded granular materials (BGMs) consist of the skeleton of granular aggregate particles and a cementitious agent between the particles. This study presents a damage analysis framework for BGMs using a computed tomography (CT) image-based discrete element method. X-ray CT scanning and image processing technology are employed to construct in situ discrete element models, categorizing particle bonding into interface bonding and mortar paste bonding. The force chains are decomposed into compressive, tensile, and shear chains, and their evolution is analyzed using complex network theory. The study investigates the coupled mechanisms between force chain evolution and bonding failure in BGMs under loading. Key findings include the critical role of fracture resistance differences between interface bonding and mortar paste bonding in determining BGM strength. The bonding-failure rate, particularly at peak load, reliably indicates material strength, with interface bonding failure surpassing mortar paste bonding failure. Force chain evolution shows rapid concentration in the loading zone, forming a stable backbone structure despite post-peak degradation. Shear and tensile force chain networks evolve with crack propagation, and force transmission paths are reorganized in the final phase. The critical points of force chain network parameters align with macroscopic load response and fracture evolution, thereby offering insights into coupled load transfer and fracture mechanisms. The proposed framework not only advances the understanding of force chain dynamics but also supports damage prediction and structural optimization by providing a comprehensive tool to track materials behavior under loading.

粘结颗粒材料（BGMs）由颗粒骨料骨架和颗粒之间的胶凝剂组成。本研究提出了一种基于计算机断层扫描（CT）图像的离散元方法的bgm损伤分析框架。采用x射线CT扫描和图像处理技术构建原位离散单元模型，将颗粒粘结分为界面粘结和砂浆粘结。将力链分解为压缩链、拉伸链和剪切链，并利用复杂网络理论分析了它们的演化过程。研究了应力链演化与粘结破坏的耦合机制。主要发现包括界面粘结和砂浆粘结之间的抗断裂性差异在决定BGM强度中的关键作用。粘结失败率，特别是在峰值荷载下，可靠地反映了材料的强度，界面粘结破坏超过砂浆粘结破坏。力链演化在加载区迅速集中，峰后退化形成稳定的骨干结构。剪切和拉伸力链网络随着裂纹扩展而演化，力传递路径在最后阶段重新组织。力链网络参数的临界点与宏观载荷响应和断裂演化一致，从而深入了解载荷耦合传递和断裂机制。提出的框架不仅促进了对力链动力学的理解，而且通过提供一个全面的工具来跟踪材料在载荷下的行为，支持损伤预测和结构优化。

{"title":"Coupled mechanisms of load transfer and fracture in bonded granular materials using CT image-based discrete element method","authors":"Hu Yang , Yiik Diew Wong , Liyan Shan , Lingwen Li","doi":"10.1016/j.compgeo.2025.107891","DOIUrl":"10.1016/j.compgeo.2025.107891","url":null,"abstract":"<div><div>Bonded granular materials (BGMs) consist of the skeleton of granular aggregate particles and a cementitious agent between the particles. This study presents a damage analysis framework for BGMs using a computed tomography (CT) image-based discrete element method. X-ray CT scanning and image processing technology are employed to construct in situ discrete element models, categorizing particle bonding into interface bonding and mortar paste bonding. The force chains are decomposed into compressive, tensile, and shear chains, and their evolution is analyzed using complex network theory. The study investigates the coupled mechanisms between force chain evolution and bonding failure in BGMs under loading. Key findings include the critical role of fracture resistance differences between interface bonding and mortar paste bonding in determining BGM strength. The bonding-failure rate, particularly at peak load, reliably indicates material strength, with interface bonding failure surpassing mortar paste bonding failure. Force chain evolution shows rapid concentration in the loading zone, forming a stable backbone structure despite post-peak degradation. Shear and tensile force chain networks evolve with crack propagation, and force transmission paths are reorganized in the final phase. The critical points of force chain network parameters align with macroscopic load response and fracture evolution, thereby offering insights into coupled load transfer and fracture mechanisms. The proposed framework not only advances the understanding of force chain dynamics but also supports damage prediction and structural optimization by providing a comprehensive tool to track materials behavior under loading.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107891"},"PeriodicalIF":6.2,"publicationDate":"2026-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146026151","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}

引用次数: 0

Utilizing physics-informed neural network and geotechnical distance field for solving three-dimensional nonlinear consolidation 利用物理信息神经网络和岩土距离场求解三维非线性固结

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-01-23 DOI: 10.1016/j.compgeo.2026.107926

Tran-Gia-Khiem Nguyen, Jongmuk Won

Solving three-dimensional (3D) nonlinear consolidation is complex and computationally expensive. This study proposes a framework for solving 3D nonlinear consolidation by utilizing an improved physics-informed neural network with hard constraints coupling with machine learning based geotechnical distance functions for three-dimensional spatial interpolation. The performance of the developed framework was assessed by comparing pore water pressure data between the developed framework and those obtained from COMSOL Multiphysics. In addition, the impact of vertical hydraulic conductibility heterogeneity, compression index, and void ratio on long-term settlement was also evaluated and discussed. It was found that the proposed framework showed a reliable estimation of the 3D distribution of pore water pressure across the 3D domain, achieving results that are comparable to data obtained from COMSOL. In addition, the heterogeneity of hydraulic conductivity can be successfully considered using the developed framework, which enables assessing the long-term settlement of a clay deposit with high uncertainty of hydraulic conductivity. Overall, the developed framework shown in this study can be applied to complex consolidation problems with low computational costs and high accuracy.

求解三维（3D）非线性固结是复杂和计算昂贵的。本研究提出了一个解决三维非线性固结的框架，该框架利用改进的物理信息神经网络与硬约束耦合，基于机器学习的岩土距离函数用于三维空间插值。通过比较开发框架与COMSOL Multiphysics获得的孔隙水压力数据，对开发框架的性能进行了评估。此外，还评价和讨论了竖向导电性非均质性、压缩指数和孔隙比对长期沉降的影响。研究发现，所提出的框架可以可靠地估计孔隙水压力在三维域中的三维分布，其结果与COMSOL获得的数据相当。此外，利用所开发的框架可以成功地考虑水力导电性的非均质性，从而能够评估具有高水力导电性不确定性的粘土沉积物的长期沉降。总的来说，本研究中显示的开发框架可以应用于复杂的固结问题，计算成本低，精度高。

{"title":"Utilizing physics-informed neural network and geotechnical distance field for solving three-dimensional nonlinear consolidation","authors":"Tran-Gia-Khiem Nguyen, Jongmuk Won","doi":"10.1016/j.compgeo.2026.107926","DOIUrl":"10.1016/j.compgeo.2026.107926","url":null,"abstract":"<div><div>Solving three-dimensional (3D) nonlinear consolidation is complex and computationally expensive. This study proposes a framework for solving 3D nonlinear consolidation by utilizing an improved physics-informed neural network with hard constraints coupling with machine learning based geotechnical distance functions for three-dimensional spatial interpolation. The performance of the developed framework was assessed by comparing pore water pressure data between the developed framework and those obtained from COMSOL Multiphysics. In addition, the impact of vertical hydraulic conductibility heterogeneity, compression index, and void ratio on long-term settlement was also evaluated and discussed. It was found that the proposed framework showed a reliable estimation of the 3D distribution of pore water pressure across the 3D domain, achieving results that are comparable to data obtained from COMSOL. In addition, the heterogeneity of hydraulic conductivity can be successfully considered using the developed framework, which enables assessing the long-term settlement of a clay deposit with high uncertainty of hydraulic conductivity. Overall, the developed framework shown in this study can be applied to complex consolidation problems with low computational costs and high accuracy.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107926"},"PeriodicalIF":6.2,"publicationDate":"2026-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146026134","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}

引用次数: 0

High-fidelity digital modeling and comparative morphological analysis of Chang’E-5 and Chang’E-6 lunar grains “嫦娥五号”和“嫦娥六号”月球颗粒高保真数字建模及形态对比分析

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-01-23 DOI: 10.1016/j.compgeo.2026.107933

Ziwei Tian , Songzheng Yu , Hanyu Wang , Guang Zhang , Yiwei Liu , Yizhuo He , Quan Zheng , Guodong Liu , Xin Liu , Ronghua Pang , Guanghui Liu , Shijing He , Yang Li , Peng Zhang

Using high-resolution X-ray microscopy (XRM), the 3D structures of lunar regolith from the Chang’E-5 and Chang’E-6 missions were reconstructed, and a four-stage CNN-based volumetric segmentation framework was established to enable accurate extraction of tens of thousands of particles from both samples. Quantitative analysis shows that Chang’E-5 particles exhibit more elongated and irregular shapes (Aspect Ratio = 0.612, Elongation Index = 0.727, Flatness Index = 0.845) and higher energy contained in the high-degree spherical harmonic spectrum (0.0346) than Chang’E-6 (Aspect Ratio = 0.623, Elongation Index = 0.732, Flatness Index = 0.854; high-degree spectral energy = 0.0299), indicating rougher surfaces and lower maturity. XRM-derived porosity measurements further reveal a pronounced structural contrast, with porosities of 38.37 % for Chang’E-5 and 54.44 % for Chang’E-6, which leads to a substantial difference in their estimated bearing capacities. These results establish a direct quantitative link between micro-scale grain morphology, macro-scale regolith structure, and surface mechanical behavior, providing critical constraints for landing-site evaluation and lunar surface infrastructure design.

利用高分辨率x射线显微镜（XRM），重建了嫦娥五号和嫦娥六号任务月球风化层的三维结构，并建立了一个基于cnn的四级体积分割框架，以便从两个样本中精确提取数万个颗粒。定量分析表明，与嫦娥6号（长径比= 0.623，延伸指数= 0.732，平整度指数= 0.854，高度光谱能量= 0.0299）相比，嫦娥5号颗粒的形状更拉长、不规则（宽高比= 0.612，延伸指数= 0.727，平整度指数= 0.845），高度球谐波谱能量（0.0346）更高，表面粗糙，成熟度较低。xrm孔隙度测量进一步揭示了明显的结构差异，嫦娥五号和嫦娥六号的孔隙度分别为38.37%和54.44%，这导致了两者估算的承载能力存在很大差异。这些结果在微观尺度的颗粒形态、宏观尺度的风化层结构和表面力学行为之间建立了直接的定量联系，为着陆点评估和月球表面基础设施设计提供了关键约束。

{"title":"High-fidelity digital modeling and comparative morphological analysis of Chang’E-5 and Chang’E-6 lunar grains","authors":"Ziwei Tian , Songzheng Yu , Hanyu Wang , Guang Zhang , Yiwei Liu , Yizhuo He , Quan Zheng , Guodong Liu , Xin Liu , Ronghua Pang , Guanghui Liu , Shijing He , Yang Li , Peng Zhang","doi":"10.1016/j.compgeo.2026.107933","DOIUrl":"10.1016/j.compgeo.2026.107933","url":null,"abstract":"<div><div>Using high-resolution X-ray microscopy (XRM), the 3D structures of lunar regolith from the Chang’E-5 and Chang’E-6 missions were reconstructed, and a four-stage CNN-based volumetric segmentation framework was established to enable accurate extraction of tens of thousands of particles from both samples. Quantitative analysis shows that Chang’E-5 particles exhibit more elongated and irregular shapes (Aspect Ratio = 0.612, Elongation Index = 0.727, Flatness Index = 0.845) and higher energy contained in the high-degree spherical harmonic spectrum (0.0346) than Chang’E-6 (Aspect Ratio = 0.623, Elongation Index = 0.732, Flatness Index = 0.854; high-degree spectral energy = 0.0299), indicating rougher surfaces and lower maturity. XRM-derived porosity measurements further reveal a pronounced structural contrast, with porosities of 38.37 % for Chang’E-5 and 54.44 % for Chang’E-6, which leads to a substantial difference in their estimated bearing capacities. These results establish a direct quantitative link between micro-scale grain morphology, macro-scale regolith structure, and surface mechanical behavior, providing critical constraints for landing-site evaluation and lunar surface infrastructure design.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107933"},"PeriodicalIF":6.2,"publicationDate":"2026-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146026152","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}

引用次数: 0

A random simplified analysis method to evaluate braced excavation-induced wall deflection considering spatial variability of soil 考虑土体空间变异性的支挖墙挠度随机简化分析方法

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-01-23 DOI: 10.1016/j.compgeo.2026.107917

Xiaolu Gan , Nianwu Liu , Jin Guo , Wenbo Li , Jianlin Yu , Xiaonan Gong

The wall deflection in braced excavations is inevitably and significantly influenced by the inherent spatial variability of soil parameters, and the applicability of current beam-spring model-based simplified methods becomes severely constricted when considering soil spatial variability. Therefore, a random simplified analysis method for the evaluation of the braced excavation-induced wall deflection in spatially variable soil is proposed in this paper. The deterministic model is established based on an improved beam-spring model that incorporates spatial variations in soil parameters, and its finite difference solution is derived. The spatial variability of soil parameters is characterized by the lognormal random field, which is further discretized using the Karhunen–Loève expansion method. The sparse polynomial chaos expansion is employed to construct an efficient metamodel for rapid random analysis. The effectiveness of the proposed method is demonstrated through a case study of a well-documented braced excavation project. Random analysis results reveal that the variability in the normalized maximum wall deflection initially decreases and then increases, while that in its position increases abruptly before progressively decreasing, due to the coupled effect of the increased excavation depth and the addition of internal supports. The ignorance of the soil spatial variability may result in overestimation of the variability in the normalized deflection and underestimation of that in its position. The presented methodology provides a useful reference for the uncertainty analysis of excavation-induced wall deflection during preliminary design stages.

土体参数固有的空间变异性对支撑式基坑墙体挠度的影响是不可避免的和显著的，在考虑土体空间变异性时，现有基于梁-弹簧模型的简化方法的适用性受到严重限制。因此，本文提出了一种评价空间变土中支撑开挖引起的墙体挠度的随机简化分析方法。在考虑土壤参数空间变化的改进梁-弹簧模型的基础上，建立了确定性模型，并推导了其有限差分解。土壤参数的空间变异性表现为对数正态随机场，并采用karhunen - lo展开法对其进行离散化。利用稀疏多项式混沌展开构造了一个快速随机分析的有效元模型。所提出的方法的有效性通过一个有充分记录的支撑挖掘项目的案例研究来证明。随机分析结果表明，由于开挖深度增加和内支护的耦合作用，归一化最大墙体挠度的变异性先减小后增大，而其位置的变异性先增大后逐渐减小。对土壤空间变异性的忽视可能导致对归一化挠度变异性的高估和对其位置变异性的低估。该方法为初步设计阶段开挖引起的墙体挠度的不确定性分析提供了有益的参考。

{"title":"A random simplified analysis method to evaluate braced excavation-induced wall deflection considering spatial variability of soil","authors":"Xiaolu Gan , Nianwu Liu , Jin Guo , Wenbo Li , Jianlin Yu , Xiaonan Gong","doi":"10.1016/j.compgeo.2026.107917","DOIUrl":"10.1016/j.compgeo.2026.107917","url":null,"abstract":"<div><div>The wall deflection in braced excavations is inevitably and significantly influenced by the inherent spatial variability of soil parameters, and the applicability of current beam-spring model-based simplified methods becomes severely constricted when considering soil spatial variability. Therefore, a random simplified analysis method for the evaluation of the braced excavation-induced wall deflection in spatially variable soil is proposed in this paper. The deterministic model is established based on an improved beam-spring model that incorporates spatial variations in soil parameters, and its finite difference solution is derived. The spatial variability of soil parameters is characterized by the lognormal random field, which is further discretized using the Karhunen–Loève expansion method. The sparse polynomial chaos expansion is employed to construct an efficient metamodel for rapid random analysis. The effectiveness of the proposed method is demonstrated through a case study of a well-documented braced excavation project. Random analysis results reveal that the variability in the normalized maximum wall deflection initially decreases and then increases, while that in its position increases abruptly before progressively decreasing, due to the coupled effect of the increased excavation depth and the addition of internal supports. The ignorance of the soil spatial variability may result in overestimation of the variability in the normalized deflection and underestimation of that in its position. The presented methodology provides a useful reference for the uncertainty analysis of excavation-induced wall deflection during preliminary design stages.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107917"},"PeriodicalIF":6.2,"publicationDate":"2026-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146026135","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}

引用次数: 0

Stochastic simulation of three-dimensional unsaturated flow in water repellent heterogeneous soil 拒水非均质土中三维非饱和流动的随机模拟

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-01-21 DOI: 10.1016/j.compgeo.2026.107938

Evan John Ricketts , Peter John Cleall , Anthony Jefferson , Pierre Kerfriden , Paul Lyons

It has been established that spatial variability of material properties can lead to distinct unstable hydraulic behaviour, and that this is prominent in soils due to their large spatial heterogeneity. This characteristic can lead to non-uniform wetting behaviour and is particularly prominent when the wettability of the medium is also non-uniform. In water repellent soil, its wettability is often spatially varying, such that a network of flow paths is created where fluid can move preferentially, leading to fingered flow patterns. In this study, the development of a model to represent moisture transport in hydrophobic soil is presented. Local spatial variations in material properties are represented by Gaussian random fields as part of a stochastic finite element based model. Key components of the model include an approach to represent the transition region between wettable and non-wettable layers, and the adoption of a suitable saturation–capillary pressure relationship. For wettable soil, this can be achieved with the standard van Genuchten relation. For hydrophobic soil, this is not applicable; thus, an alternative is employed. The model is then validated against field-scale experimental observations by Lipsius and Mooney (2006), which examined the impact of soil heterogeneity on infiltration profiles. The results demonstrate the model’s ability to capture complex flow dynamics in hydrophobic soils, extending the understanding of moisture transport in heterogeneous soils by explicitly modelling the spatial variability of wettability and its impact on soil hydraulic response.

已经确定，材料特性的空间变异性会导致不同的不稳定水力行为，并且由于其巨大的空间异质性，这在土壤中是突出的。这一特性可能导致不均匀的润湿行为，当介质的润湿性也不均匀时，这一特性尤为突出。在疏水性土壤中，其润湿性通常是空间变化的，因此流体可以优先移动的流动路径网络会形成指状流动模式。在这项研究中，提出了一个模型来表示水分在疏水土壤中的输送。材料性能的局部空间变化用高斯随机场表示，作为随机有限元模型的一部分。该模型的关键组成部分包括一种表示可湿层和不可湿层之间过渡区域的方法，以及采用合适的饱和度-毛细管压力关系。对于可湿性土壤，这可以用标准的van Genuchten关系来实现。对于疏水性土壤，这是不适用的；因此，采用了另一种方法。然后，Lipsius和Mooney（2006）根据现场尺度的实验观察验证了该模型，该实验研究了土壤异质性对入渗剖面的影响。结果表明，该模型能够捕捉疏水性土壤中复杂的流动动力学，通过明确模拟润湿性的空间变异性及其对土壤水力响应的影响，扩展了对非均质土壤中水分输送的理解。

{"title":"Stochastic simulation of three-dimensional unsaturated flow in water repellent heterogeneous soil","authors":"Evan John Ricketts , Peter John Cleall , Anthony Jefferson , Pierre Kerfriden , Paul Lyons","doi":"10.1016/j.compgeo.2026.107938","DOIUrl":"10.1016/j.compgeo.2026.107938","url":null,"abstract":"<div><div>It has been established that spatial variability of material properties can lead to distinct unstable hydraulic behaviour, and that this is prominent in soils due to their large spatial heterogeneity. This characteristic can lead to non-uniform wetting behaviour and is particularly prominent when the wettability of the medium is also non-uniform. In water repellent soil, its wettability is often spatially varying, such that a network of flow paths is created where fluid can move preferentially, leading to fingered flow patterns. In this study, the development of a model to represent moisture transport in hydrophobic soil is presented. Local spatial variations in material properties are represented by Gaussian random fields as part of a stochastic finite element based model. Key components of the model include an approach to represent the transition region between wettable and non-wettable layers, and the adoption of a suitable saturation–capillary pressure relationship. For wettable soil, this can be achieved with the standard van Genuchten relation. For hydrophobic soil, this is not applicable; thus, an alternative is employed. The model is then validated against field-scale experimental observations by <span><span>Lipsius and Mooney (2006)</span></span>, which examined the impact of soil heterogeneity on infiltration profiles. The results demonstrate the model’s ability to capture complex flow dynamics in hydrophobic soils, extending the understanding of moisture transport in heterogeneous soils by explicitly modelling the spatial variability of wettability and its impact on soil hydraulic response.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107938"},"PeriodicalIF":6.2,"publicationDate":"2026-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146022685","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}

引用次数: 0