A. Yuksel-Ozan, K. Chang, G. Gunes, G. Constantinescu
{"title":"Flow Hydrodynamics, Sediment Erosion Potential and Drag Forces In An Open Channel With an Array of Vegetation Patches at Its Banks","authors":"A. Yuksel-Ozan, K. Chang, G. Gunes, G. Constantinescu","doi":"10.1016/j.advwatres.2026.105481","DOIUrl":"https://doi.org/10.1016/j.advwatres.2026.105481","url":null,"abstract":"Detached eddy simulations are conducted using a viscous flow solver to understand flow hydrodynamics and turbulence structure in a horizontal bed, straight open channel of width W containing an array of equally spaced, semi-circular porous patches at its two banks placed in staggered way. In the numerical experiments, the streamwise distance between successive patches on opposite banks, L, is equal to the diameter of the patch, D<ce:inf loc=\"post\">p</ce:inf>. The patches containing uniformly distributed emerged, solid cylinders mimic vegetation canopies in river channels that can be used in river restoration projects to induce meandering of straight, degraded channels or to stabilize alternate bars. The paper investigates how the relative patch size, D<ce:inf loc=\"post\">p</ce:inf>/W, and the solid volume fraction of the patches, SVF, affect the patch shear layers and, for sufficiently high SVF, the continuous mixing layers forming on both sides of the meandering corridor of high velocities generated inside the central part of the channel, the sediment entrainment capacity of the flow beneath the corridor of high velocities, the flow discharges through the patches, and the additional channel drag due to the cylinders/plant stems. Results show that the second patch in the array is subject to the largest streamwise drag force and this force can be up to 3 to 4 times larger than the force acting on the higher-rank patches. The paper also discusses analogies and differences between the mixing layers generated by the two arrays of porous patches and mixing layers generated at the interface between a groyne field and the open water region of a channel.","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"60 1","pages":""},"PeriodicalIF":4.7,"publicationDate":"2026-08-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148884816","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A Novel Two-Scale Continuum Simulation of Thermo-Hydro-Chemical Phenomena of Matrix Acidizing Operation using Element Free Galerkin (EFG) Method","authors":"Seyed Raef Sheikhahmadi, Ali Pak","doi":"10.1016/j.advwatres.2026.105478","DOIUrl":"https://doi.org/10.1016/j.advwatres.2026.105478","url":null,"abstract":"Numerical modeling of reservoir stimulation using matrix acidizing requires a thorough understanding of chemical reactions, fluid flow, heat transfer, and their interactions in porous rock. This research introduces a robust in-house mesh-free code utilizing Element Free Galerkin (EFG) method to address Chemical, Hydro-Chemical (HC), and Thermo-Hydro-Chemical (THC) problems, particularly used for simulation of matrix acidizing operations with a two-scale continuum approach. Incorporating thermal interactions improves the realism of simulations by considering the heat from exothermic reactions. This study also introduces a combination of Equivalent Continuum Method (ECM) with EFG method to examine the impact of the rock microcracks on the matrix acidizing process. The verification of the developed EFG code is carried out first by solving an analytical chemical transport problem (C). Then two HC, and one THC matrix acidizing problems are solved and the obtained results confirm its ability to model acid flow, chemical and thermal effects, highlighting its efficacy for simulating various matrix acidizing problems. In the EFG-ECM simulation of the effects of microcracks in a coupled THC study, 0.15%, 0.71% and 12.28% reductions in Pore Volume to Break Through (PVBT) were observed for 0<ce:sup loc=\"post\">o</ce:sup>, 45<ce:sup loc=\"post\">o</ce:sup> and 90° crack orientations respectively, in comparison with similar cases without thermal effects. The results highlight a minimal temperature impact on the PVBT for 0<ce:sup loc=\"post\">o</ce:sup> and 45° crack inclination cases, due to unchanged dissolution pattern; however, temperature variation influences acid kinetics and dissolution pattern in the vertical cracks case, resulting in a noticeably lower PVBT.","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"225 1","pages":""},"PeriodicalIF":4.7,"publicationDate":"2026-08-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148884817","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Anfal Al Zarafi, Ahmed AlZaabi, Rukuan Chai, Branko Bijeljic, Martin J. Blunt
{"title":"Low CO₂ relative permeability in Bentheimer sandstone revealed by steady-state measurements and pore-scale imaging","authors":"Anfal Al Zarafi, Ahmed AlZaabi, Rukuan Chai, Branko Bijeljic, Martin J. Blunt","doi":"10.1016/j.advwatres.2026.105479","DOIUrl":"https://doi.org/10.1016/j.advwatres.2026.105479","url":null,"abstract":"Quantifying pore-scale fluid displacement mechanisms in CO<ce:inf loc=\"post\">2</ce:inf>/brine systems is essential for predicting multiphase flow and trapping during CO<ce:inf loc=\"post\">2</ce:inf> storage. In this study, we imaged the steady-state flow of brine and CO<ce:inf loc=\"post\">2</ce:inf> in a water-wet Bentheimer sandstone sample. We devised an experimental method based on differential imaging to examine the pore-scale behaviour of CO<ce:inf loc=\"post\">2</ce:inf> within the system at reservoir conditions during drainage. This allowed us to measure relative permeability and capillary pressure while observing the evolution of the gas ganglia in the pore space at different fractional flows under capillary-dominated conditions.","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"29 1","pages":""},"PeriodicalIF":4.7,"publicationDate":"2026-08-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148884819","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Foojan Kazemzadeh Haghighi, Achyut Mishra, Jay R. Black, Edward M. Hinton, Ralf R. Haese
{"title":"Heterogeneity of pore geometry and directional fluid flow behavior in vesicular basalts","authors":"Foojan Kazemzadeh Haghighi, Achyut Mishra, Jay R. Black, Edward M. Hinton, Ralf R. Haese","doi":"10.1016/j.advwatres.2026.105473","DOIUrl":"https://doi.org/10.1016/j.advwatres.2026.105473","url":null,"abstract":"Characterizing pore volume, connectivity, and multi-phase flow properties in subsurface rocks is fundamental to understanding fluid transport mechanisms. The latter is crucially important for the effective development of strategies for the utilization of subsurface resources such as groundwater aquifers, hydrocarbon and CO<ce:inf loc=\"post\">2</ce:inf> storage reservoirs. This study investigates morphological and flow attributes in vesicular basalts through an analysis of pore-scale architecture and directional fluid flow utilizing processed and segmented micro-CT images. Incorporating pore shape clustering, shape factor analysis and direct numerical simulations, we characterize the internal anisotropy and heterogeneity of these basalts and compare their properties with the commonly studied sedimentary lithologies. The investigation of morphological heterogeneity revealed a size-dependent transition in pore geometries. Following an increase in the longest Feret diameter, a convergence toward cubic and cuboidal voids occurs, suggesting a geometric limit of complexity within the volcanic matrix, defining a unique structural signature of these rocks. This plateau limit value varies over different samples, highlighting the existing heterogeneity. Permeability anisotropy is assessed through center-based fluid injection simulations in micro-CT images, designed to replicate radial flow conditions representative of subsurface reservoir scenarios, which are validated against pore network-derived directional metrics. Analysis of azimuthal and polar distributions of both velocity vectors and pore-center-to-pore-center vectors shows well-defined principal flow orientations across all investigated sub-volumes. Furthermore, the full permeability tensors resolved for each sub-volume demonstrate extreme transverse hydrodynamic coupling, with off-diagonal components exceeding principal diagonal entries by several orders of magnitude in select sub-volumes, confirming the existence of oriented hydraulic conduits in vesicular basalts. In contrast, such behavior was not observed in sandstone and limestone. The substantial variability in permeability and principal flow axes across and within samples represents extreme hydraulic heterogeneity and anisotropy in vesicular basalts with significant implications for fluid flow.","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"38 1","pages":""},"PeriodicalIF":4.7,"publicationDate":"2026-08-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148853059","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Yi-fan Xia , Zhong-kai Feng , Yang Xiao , Tao-tao Zhang , Ling-zhong Kong , Wen-jing Niu , Hui-ming Zhang
{"title":"Mixed Copula for streamflow simulation based on intelligent knowledge set and parameter calibration using cooperation search algorithm","authors":"Yi-fan Xia , Zhong-kai Feng , Yang Xiao , Tao-tao Zhang , Ling-zhong Kong , Wen-jing Niu , Hui-ming Zhang","doi":"10.1016/j.advwatres.2026.105237","DOIUrl":"10.1016/j.advwatres.2026.105237","url":null,"abstract":"<div><div>Streamflow simulation plays a crucial role in the scientific management of reservoir operations and water resources. This study proposes a novel streamflow simulation method using a mixed Copula approach, coupled with an intelligent knowledge set and a cooperation search algorithm for parameter optimization. This method addresses the limitations of traditional Copula-based models, such as limited expressive capability, an excessive number of model parameters, and long computation times. Initially, an intelligent knowledge set is established using a time series model. Next, the cooperation search algorithm is employed to estimate model parameters. Finally, simulated streamflow values are generated via conditional Copula. Theoretical analysis shows significant reduction in computational complexity compared to traditional methods. Engineering applications at multiple hydrological stations validate that the proposed method reduces parameters by over 99%, shortens estimation time to 3–4% of traditional methods, and improves accuracy by 5–15%. This novel method integrates an intelligent knowledge set with parameter optimization to improve the precision of streamflow simulation, providing a valuable technical tool for watershed management.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"210 ","pages":"Article 105237"},"PeriodicalIF":4.2,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146138401","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Haotian Li , Saideep Pavuluri , Harris Sajjad Rabbani , Bicheng Yan
{"title":"Multi-resolution transfer learning for rapid prediction of pore-scale multiphase flow","authors":"Haotian Li , Saideep Pavuluri , Harris Sajjad Rabbani , Bicheng Yan","doi":"10.1016/j.advwatres.2026.105236","DOIUrl":"10.1016/j.advwatres.2026.105236","url":null,"abstract":"<div><div>Deciphering multiphase flow patterns at the pore-scale is fundamentally essential for upscaling to determine macro-scale flow parameters. Direct numerical simulations provide detailed insights related to pore-scale flow physics but are computationally expensive. To reduce computational costs, coarser meshes may be used at the expense of accuracy. This study presents a deep learning framework that leverages multi-resolution data for two-dimensional pore-scale two-phase flow at fixed capillary number: low resolution simulations generate large training datasets, while high resolution simulations offer the required supervision to capture the pore-scale flow physics. The model effectively transfers the flow physics learned from low resolution dataset to the high resolution cases, requiring only limited high-fidelity data for adaptation. By combining computational efficiency with predictive accuracy, the proposed framework facilitates rapid and accurate pore-scale flow analysis, addressing a critical need in multiphase flow research.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"210 ","pages":"Article 105236"},"PeriodicalIF":4.2,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146134973","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Xianmeng Meng , Lintao Shen , Xiaoxuan Liu , Qu Wang , Maosheng Yin , Dengfeng Liu
{"title":"One-dimensional non-Darcian flow model incorporating the impact of nonlinear clay consolidation on the threshold hydraulic gradient","authors":"Xianmeng Meng , Lintao Shen , Xiaoxuan Liu , Qu Wang , Maosheng Yin , Dengfeng Liu","doi":"10.1016/j.advwatres.2026.105239","DOIUrl":"10.1016/j.advwatres.2026.105239","url":null,"abstract":"<div><div>The process of clay consolidation can alter hydraulic conductivity, which subsequently impacts the threshold hydraulic gradient, thereby affecting seepage flow. Existing seepage consolidation models have not considered the impact of changing hydraulic conductivity on the threshold hydraulic gradient. To address this issue, this paper establishes a one-dimensional non-Darcian flow model that accounts for the changes in the threshold hydraulic gradient due to the nonlinear consolidation characteristics of clay. The model is solved using the finite difference method, and the results are compared with those from a model that neglects the changes in the threshold hydraulic gradient. The results indicate that when the changes in the threshold hydraulic gradient are taken into account, both the rate of movement of the seepage moving boundary and the seepage flow velocity are reduced. The hydraulic head calculated with consideration of the threshold hydraulic gradient changes is higher than that calculated without considering such changes. The discrepancies in the position of the seepage moving boundary and the threshold hydraulic gradient collectively dictate the variations in the hydraulic head difference and the seepage flow velocity difference. When the initial hydraulic conductivity is small, the initial void ratio is large, the compression index is large, and the permeability index is small, the differences in hydraulic head between the model accounting for changes in the threshold hydraulic gradient and the one that does not are more significant. Ultimately, a laboratory experiment is used to validate the developed model. Experimental simulation results indicate that ignoring the variation in the threshold hydraulic gradient in long-term seepage simulations leads to a flow prediction error of approximately 15%.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"210 ","pages":"Article 105239"},"PeriodicalIF":4.2,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146160585","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A reinterpreted discrete fracture model for Darcy–Brinkman flow in fractured porous media and its extension on nonconforming meshes","authors":"Jingyao Liu, Hui Guo , Zhaoqin Huang , Yang Yang","doi":"10.1016/j.advwatres.2026.105240","DOIUrl":"10.1016/j.advwatres.2026.105240","url":null,"abstract":"<div><div>A novel hybrid-dimensional model for Darcy–Brinkman flow in fractured porous media is proposed, inherently compatible with non-conforming grids. Building upon the reinterpreted discrete fracture model (RDFM) for Darcy flow, which introduces a Dirac–<span><math><mi>δ</mi></math></span> function to unify matrix-fracture flow, we develop the hybrid-dimensional RDFM for Darcy–Brinkman flow. We also rigorously establish its mathematical equivalence with the classical interface model. For numerical discretization, a hybrid scheme combining Local Discontinuous Galerkin (LDG) and standard Galerkin finite element methods (FEM) is employed. This approach overcomes key limitations of the LDG method in modeling one-dimensional fractures, such as difficulties with numerical flux selection and auxiliary variable specification, while maintaining computational efficiency. To solve the resulting coupled system, we introduce a pseudo-time and advance the solution in time toward a stationary state. Validation through coupled tracer transport simulations confirms the model’s robustness and applicability on non-matching grids.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"210 ","pages":"Article 105240"},"PeriodicalIF":4.2,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147386189","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Shuyao Niu , Zhike Zou , Longcang Shu , Giovanni Michele Porta
{"title":"Single collector capture model analysis of particle deposition: Effects of hydrodynamics, particle size, and collector geometry","authors":"Shuyao Niu , Zhike Zou , Longcang Shu , Giovanni Michele Porta","doi":"10.1016/j.advwatres.2026.105233","DOIUrl":"10.1016/j.advwatres.2026.105233","url":null,"abstract":"<div><div>Particle deposition in porous media is critical in natural and engineered systems and has significant implications for groundwater recharge, colloid filtration, and contaminant transport. Therefore, a comprehensive understanding of pore-scale dynamics is essential. Current research mainly focuses on deposition outcomes, whereas dynamic deposition processes require further research, particularly into the coupled effects of flow velocity, particle size, and collector geometry. This study establishes a single-collector model under favorable conditions, using Lagrangian simulations to analyze pore-scale force interactions and particle trajectories. This approach facilitates a systematic investigation of how flow velocity, particle size, and collector geometry influence deposition time, transport trajectory length, and spatial distribution of deposition particles, as well as force analysis of particles. The results showed that particle size and flow velocity primarily control η while collector geometry considerably affects deposition kinetics and spatial patterns. A distinct size-dependent transition in deposition mechanisms was observed, with sub-micron particles (≤1 μm) exhibiting stochastic, diffusion-dominated behavior with dispersed deposition times and trajectories and larger particles (≥2 μm) undergoing rapid, deterministic deposition driven by gravity. Irregular collector shapes enhanced the dispersion and delay of small particle deposition due to complex flow fields with vortices and tortuous streamlines, causing localized deposition hotspots in surface concavities. Force balance analysis shows that in the concave regions of irregular collectors, strong competition occurred Brownian motion and drag force, resulting in more random and tortuous deposition paths. As particle size increased, deposition behavior became shifts to gravity-dominated, with irregular flow fields exhibiting less influence. This study reveals the key role of local flow features in shaping particle deposition behavior, thus providing a deeper understanding of the particle deposition process.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"210 ","pages":"Article 105233"},"PeriodicalIF":4.2,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146110457","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Coupled two-phase flow and surfactant/PFAS transport in porous media with angular pores: From pore-scale physics to Darcy-scale modeling","authors":"Sidian Chen , Bo Guo , Tianyuan Zheng","doi":"10.1016/j.advwatres.2026.105222","DOIUrl":"10.1016/j.advwatres.2026.105222","url":null,"abstract":"<div><div>Two-phase surfactant-laden fluid flow and surfactant transport in porous media are critical to many natural and engineering applications. Surfactants modify two-phase flow by altering interfacial tension and wettability, while two-phase flow controls surfactant transport pathways and adsorption sites. These coupled processes are commonly modeled by combining Darcy-type two-phase flow equations with advection–dispersion–adsorption transport equations, with capillary pressure–saturation relationships scaled using the Leverett <span><math><mi>J</mi></math></span>-function. However, the Leverett <span><math><mi>J</mi></math></span>-function simplifies the porous medium as bundles of cylindrical tubes and decouples interfacial tension and wettability, limiting representation of angular pore geometries and coupled interfacial tension and wettability effects. We present a modeling framework that incorporates pore angularity and interfacial tension–wettability coupling effect into Darcy-scale surfactant-laden fluid flow and surfactant transport models. Within this framework, we derive two-phase flow properties for angular pores, upscale them across pore size distributions, and obtain explicit and closed-form expressions for the upscaled properties. These expressions are incorporated into a coupled flow–transport model for simulating transient two-phase flow and surfactant transport processes. Modeling results suggest a nonmonotonic and nonlinear dependence of two-phase flow properties on pore structure (angularity and size distribution) and interfacial tension (controlled by surfactant type and concentration). Example simulations of water flow and PFAS (surfactant-like contaminants) migration in unsaturated soils indicate that surfactant-induced flow effects on PFAS leaching are generally minor under typical site conditions, whereas pore angularity exerts dominant control on water flow, interfacial area, and consequently PFAS retention. Overall, the upscaling framework offers a more physically grounded approach for modeling two-phase surfactant-laden fluid flow and surfactant transport in porous media.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"209 ","pages":"Article 105222"},"PeriodicalIF":4.2,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146072684","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}