Jinbao Guo , Yixin Zhao , Linfeng He , Xiaodong Guo , Hua Shen , Sen Gao
{"title":"Spontaneous imbibition of unsaturated sandstone under different vertical temperature gradients: neutron radiography experiments and dynamic models","authors":"Jinbao Guo , Yixin Zhao , Linfeng He , Xiaodong Guo , Hua Shen , Sen Gao","doi":"10.1016/j.advwatres.2024.104832","DOIUrl":"10.1016/j.advwatres.2024.104832","url":null,"abstract":"<div><div>To elucidate the imbibition behavior of water in complex temperature and stress environments, spontaneous imbibition experiments were conducted on unsaturated matrix sandstones at different vertical temperature gradients by neutron radiography technology. Additionally, corresponding dynamic models of water imbibition in porous media were established. The research results reveal the phased characteristics of sandstone spontaneous imbibition and the influence of vertical temperature gradient on the evolution of wetting front. Specifically, the initial development speed of the wetting front increases with an increase in the vertical temperature gradient, indicating a direct relationship. However, the growth rate of the wetting front gradually slows down with increasing time, eventually reaching a saturated state. Model validation demonstrates that the traditional Washburn's law is still valid in isothermal conditions without temperature gradient (G=0). Further analysis indicates that the imbibition rate has a direct correlation with linear thermal expansion coefficient (α) and viscosity temperature coefficient (β) across various vertical temperature gradients, and an inverse correlation with surface tension temperature coefficient (γ). Furthermore, when the values of α, β, and γ fall below 0.001, their impact on the imbibition rate becomes negligible. The sensitivity of the imbibition rate to parameters γ, β, and α decreases in that order, with γ being the most sensitive, followed by β, and α being the least sensitive. Moreover, the relative importance of α, β, and γ dictates their specific influence on the imbibition rate, and a synergistic effect exists among these parameters, which collectively influence the water absorption behavior of sandstone.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104832"},"PeriodicalIF":4.0,"publicationDate":"2024-10-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142432342","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":"Upscaling transport in heterogeneous media featuring local-scale dispersion: Flow channeling, macro-retardation and parameter prediction","authors":"Lian Zhou, Scott K. Hansen","doi":"10.1016/j.advwatres.2024.104830","DOIUrl":"10.1016/j.advwatres.2024.104830","url":null,"abstract":"<div><div>Many theoretical treatments of transport in heterogeneous Darcy flows consider advection only. When local-scale dispersion is neglected, flux weighting persists over time; mean Lagrangian and Eulerian flow velocity distributions relate simply to each other and to the variance of the underlying hydraulic conductivity field. Local-scale dispersion complicates this relationship, potentially causing initially flux-weighted solute to experience lower-velocity regions as well as Taylor-type macrodispersion due to transverse solute movement between adjacent streamlines. To investigate the interplay of local-scale dispersion with conductivity log-variance, correlation length, and anisotropy, we perform a Monte Carlo study of flow and advective-dispersive transport in spatially-periodic 2D Darcy flows in large-scale, high-resolution multivariate Gaussian random conductivity fields. We observe flow channeling at all heterogeneity levels and quantify its extent. We find evidence for substantial effective retardation in the upscaled system, associated with increased flow channeling, and observe limited Taylor-type macrodispersion, which we physically explain. A quasi-constant Lagrangian velocity is achieved within a short distance of release, allowing usage of a simplified continuous-time random walk (CTRW) model we previously proposed in which the transition time distribution is understood as a temporal mapping of unit time in an equivalent system with no flow heterogeneity. The numerical data set is modeled with such a CTRW; we show how dimensionless parameters defining the CTRW transition time distribution are predicted by dimensionless heterogeneity statistics and provide empirical equations for this purpose.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104830"},"PeriodicalIF":4.0,"publicationDate":"2024-10-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142441967","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}
Mamta Jotkar , Ilan Ben-Noah , Juan J. Hidalgo , Marco Dentz
{"title":"Diffusiophoresis of colloids in partially-saturated porous media","authors":"Mamta Jotkar , Ilan Ben-Noah , Juan J. Hidalgo , Marco Dentz","doi":"10.1016/j.advwatres.2024.104828","DOIUrl":"10.1016/j.advwatres.2024.104828","url":null,"abstract":"<div><div>It is known that structural heterogeneity induced by the distribution of the water and air phases creates complex flow patterns with a broad distribution of flow velocities, which in turn control key aspects of transport including arrival and residence times, dispersion and spatial distributions of dissolved salts and suspended colloidal particles. Stagnation zones serve as hot spots where colloidal particles can get entrapped, providing a challenging task in controlling their transport. Recent investigations in simple confined geometries suggest diffusiophoresis, the colloid migration driven by local salt gradients, to be an efficient mechanism to control colloidal migration. However, despite its potential, diffusiophoresis in complex porous media remains poorly understood. We use detailed numerical simulations to unravel the effects of diffusiophoresis occurring at pore-scale on the macroscopic dispersion of colloids in partially-saturated porous media with different water-saturation degrees. Diffusiophoresis can promote particle retention or removal, depending on the diffusiophoretic mobility. For fully-saturated media, the pore-scale dynamics due to diffusiophoresis are manifested in the long-time tailing of the breakthrough curves. For partially-saturated media as the degree of water-saturation decreases and flow heterogeneity increases, we observe accumulation and depletion effects in the colloid breakthrough curves which can be traced back to trapping and release in dead-end zones. Finally, our results suggest that colloid mobilisation and retention due to diffusiophoresis can be controlled by the flow rate of the injected salt solution.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104828"},"PeriodicalIF":4.0,"publicationDate":"2024-09-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142419711","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 generalized kinetic theory of Ostwald ripening in porous media","authors":"Nicolas Bueno, Luis Ayala, Yashar Mehmani","doi":"10.1016/j.advwatres.2024.104826","DOIUrl":"10.1016/j.advwatres.2024.104826","url":null,"abstract":"<div><div>Partially miscible bubbles (e.g., CO<span><math><msub><mrow></mrow><mrow><mn>2</mn></mrow></msub></math></span>) trapped inside a porous medium and surrounded by a wetting phase (e.g., water) occur in a number of applications including underground hydrogen storage, geologic carbon sequestration, and the operation of electrochemcial devices such as fuel cells and electrolyzers. Such bubbles evolve due to a process called Ostwald ripening that is driven by differences in their interfacial curvature. For spherical bubbles, small bubbles shrink and vanish while feeding into larger ones, resulting in one large bubble at equilibrium. Within the confinement of a porous medium, however, bubbles can attain a distribution of sizes at equilibrium that have identical curvature. This work concerns itself with the formulation of a kinetic theory that predicts the statistical evolution of bubble <em>states</em>, defined as the sizes of the pores within which bubbles are trapped and the extent to which those pores are saturated with bubbles. The theory consists of a population balance equation and appropriate closure approximations. Systematic comparisons against a previously published pore network model (PNM) are conducted to validate the theory. Our theory generalizes existing variants in the literature limited to spherical bubbles trapped in homogeneous media to non-spherical (deformed) bubbles inside microstructures with arbitrary heterogeneity and spatial correlation in pore/throat sizes. We discuss the applicability, limitations, and implications of the theory towards future extensions.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104826"},"PeriodicalIF":4.0,"publicationDate":"2024-09-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142329763","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}
Subhojit Kadia , I.A. Sofia Larsson , Mats Billstein , Nils Rüther , Leif Lia , Elena Pummer
{"title":"Investigating supercritical flow characteristics and movement of sediment particles in a narrow channel bend using PTV and video footage","authors":"Subhojit Kadia , I.A. Sofia Larsson , Mats Billstein , Nils Rüther , Leif Lia , Elena Pummer","doi":"10.1016/j.advwatres.2024.104827","DOIUrl":"10.1016/j.advwatres.2024.104827","url":null,"abstract":"<div><div>This experimental study investigates the cause of nonuniform invert abrasion observed at sediment bypass tunnel (SBT) bends by examining the variations in velocity distributions, turbulence properties, bed shear stress, and bulk sediment movements under three supercritical bend flow conditions, detailed investigation of such flow is scarce. Using a laboratory-scaled model (1:22) of the downstream bend at Solis SBT, Switzerland, the research utilized particle tracking velocimetry and high-speed cameras with spherical sandstones and glass spheres representing sediments. The results indicate that as the secondary currents develop in the flow direction, the flow properties and sediments redistribute across the channel: the high-momentum fluids are directed toward the outer wall, the bed shear stress increases toward the outer wall, and the sediments are pushed toward the inner wall, which then follow this path downstream, even in straight sections, despite lower bed shear stress. This distribution of sediments, driven by secondary currents, leads to deeper invert abrasions toward the inner wall at SBT bends and downstream sections. Thus, these abrasions are primarily influenced by sediment movement rather than the bed shear stress alone. The study's findings are also valuable for validating future numerical simulations.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104827"},"PeriodicalIF":4.0,"publicationDate":"2024-09-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142419709","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Dynamics of bedload transport under run-up wave by gravel resolved scheme based on 3D DEM-MPS coupling","authors":"Takumi Tazaki, Eiji Harada, Hitoshi Gotoh","doi":"10.1016/j.advwatres.2024.104824","DOIUrl":"10.1016/j.advwatres.2024.104824","url":null,"abstract":"<div><div>Accurate predictions of morphological changes in swash zones require a detailed understanding of sediment transport mechanisms, which are strongly related to bore-induced vortices and turbulence, surface-subsurface interactions, namely, infiltrate/exfiltrate flow, and swash-swash interactions. However, obtaining experimental or field measurements of instantaneous velocity and sediment flux is challenging owing to the suspended sediment, turbulence, and shallow depth characteristics of these regions. The present study simulates the gravel bedload transport under a dam-break bore at a grain-resolved spatial scale. The simulation uses a 3D Lagrangian–Lagrangian solid–fluid coupled model comprising the moving particle semi-implicit (MPS) method for a violent swash flow and the discrete element method (DEM) for gravels. The simulated water depth, velocity, and sediment flux agree with existing experimental results during a run-up. The gravel transport mechanisms for the lower, mid, and upper swash zones were discussed. Discussions on bedload mechanisms reveal that bore-generated horizontal vortices can reduce the onshore velocity near the beach surface, reducing sediment flux in the lower swash zone. Modified Shields numbers investigate the seepage effects: the frequently used standard Shields number value is insufficient to estimate bedload flux under the intense infiltration in the mid-swash zone. The simulation result also elucidates the turbulence characteristics in the upper swash zone.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104824"},"PeriodicalIF":4.0,"publicationDate":"2024-09-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142322774","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":"Density-driven free convection in heterogeneous aquifers with connectivity features","authors":"Hongkai Li , Lei Zhang , Yu Ye , Chunhui Lu","doi":"10.1016/j.advwatres.2024.104825","DOIUrl":"10.1016/j.advwatres.2024.104825","url":null,"abstract":"<div><div>Free convection usually happens in variable-density groundwater flow systems, and it favors contaminant transport by enlarging length scales and shortening timescales compared to advection and diffusion/dispersion alone. Previous studies have demonstrated that heterogeneity with multi-variate Gaussian distribution for logarithmic permeabilities (log<sub>10</sub><em>k</em>) plays an important role in the onset, growth, and/or decay of instability during the density-driven convective process. Nevertheless, the connectivity features (i.e., connected structures of extremely high/low-<em>k</em> values), which are common in natural aquifers, have received little attention. In this study, the classical problem of transient free convection has been modified and numerically simulated by Monte Carlo approach to investigate the effects of connectivity features in heterogeneity on the unstable convective processes. Results show that free convection is promoted by the connected high-<em>k</em> structures and retarded by the connected low-<em>k</em> structures during mainly the early-stage mass loading. The impacts of connectivity features tend to be amplified by higher variation in log<sub>10</sub><em>k</em> distributions, and can be secondarily influenced by correlation lengths and anisotropy. Under the multi-variate Gaussian assumption, the existence of connected high-<em>k</em> structures leads to underestimation of density-driven instability, in which the risk differs based on the statistics of permeability fields, metrics of interest and timeframe. This study highlights the importance of understanding connectivity features in heterogeneous geological media when assessing density-dependent solute transport in groundwater systems.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104825"},"PeriodicalIF":4.0,"publicationDate":"2024-09-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142327841","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}
Youssra Rahham, Stephen Dauphinais, Jeff T. Gostick, Marios A. Ioannidis
{"title":"Nanoparticle transport in partially saturated porous media: Attachment at fluid interfaces","authors":"Youssra Rahham, Stephen Dauphinais, Jeff T. Gostick, Marios A. Ioannidis","doi":"10.1016/j.advwatres.2024.104816","DOIUrl":"10.1016/j.advwatres.2024.104816","url":null,"abstract":"<div><div>Like the solid-water interface (SWI), air-water and oil-water interfaces (AWI and OWI) also act as collectors for nano-sized particles in porous media. The attachment of hydrophobic nanoparticles, which is often favorable and irreversible, is of particular interest because the transport and retention of such particles is closely linked to the fate of nanoplastics in unsaturated subsurface environments and the success of nanoremediation practices. Here, we show how a pore-network model (PNM) can be used to upscale the kinetics and extent of irreversible nanoparticle attachment at a single fluid-fluid interface under conditions of advection and dispersion in a sphere packing. By focusing on a trapped (immobile) non-wetting phase, we highlight a fundamental difference between the single-collector contact efficiency of AWI/OWI and SWI. Namely, AWI/OWI collectors, which are largely by-passed by the flowing aqueous phase, are exposed to a hydrodynamic environment dominated by diffusion. This difference has profound implications for the modelling of nanoparticle transport in porous media at the continuum (Darcy) scale. This study reveals the potential of pore network modelling as an essential complement to continuum models for upscaling the behavior of nanocolloids in porous media.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104816"},"PeriodicalIF":4.0,"publicationDate":"2024-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142327840","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Numerical modeling of transient water table in shallow unconfined aquifers: A hyperbolic theory and well-balanced finite volume scheme","authors":"Ying-Hsin Wu, Eiichi Nakakita","doi":"10.1016/j.advwatres.2024.104820","DOIUrl":"10.1016/j.advwatres.2024.104820","url":null,"abstract":"<div><div>We present a new methodology capable of modeling transient motion of shallow phreatic surface of groundwater in unconfined aquifers. This methodology is founded on a new and comprehensive theory for water table motion and a corresponding efficient numerical scheme. In the theoretical aspect, we derived a new set of governing equations constituted by a depth-averaged continuity equation and momentum equations based on unsteady Darcy’s law. The derived governing equations are of the hyperbolic type and possess stiff terms in the momentum equations due to the inertia motion in a characteristic time scale that is relatively shorter than the time scale of seepage motion. To effectively solve the derived hyperbolic system with stiff terms, in the numerical aspect, we utilize <em>f</em>-wave propagation algorithm, an explicit finite volume method, that can ensure numerical convergence and well-balancing solutions when momentum is rapidly relaxing to an equilibrium of steady state. Verification is successfully performed by comparing the results with analytic solutions to the classic problem of multidimensional spreading of a groundwater mound. This study demonstrates that the proposed methodology can accurately and satisfactorily simulate the spatiotemporal distribution of shallow water table and its wetting front in unconfined aquifers.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104820"},"PeriodicalIF":4.0,"publicationDate":"2024-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0309170824002070/pdfft?md5=a039eb9377fc4092d60f707501aeabf4&pid=1-s2.0-S0309170824002070-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142275669","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Nadav Peleg , Daniel B. Wright , Hayley J. Fowler , João P. Leitão , Ashish Sharma , Francesco Marra
{"title":"A simple and robust approach for adapting design storms to assess climate-induced changes in flash flood hazard","authors":"Nadav Peleg , Daniel B. Wright , Hayley J. Fowler , João P. Leitão , Ashish Sharma , Francesco Marra","doi":"10.1016/j.advwatres.2024.104823","DOIUrl":"10.1016/j.advwatres.2024.104823","url":null,"abstract":"<div><div>Hydrologists and civil engineers often use design storms to assess flash flood hazards in urban, rural, and mountainous catchments. These synthetic storms are not representations of real extreme rainfall events, but rather simplified versions parameterized to mimic extreme precipitation statistics often obtained from intensity–duration–frequency (IDF) curves. To construct design storms for the future climate, it is thus necessary first to recalculate IDF curves to represent rainfall under warmer conditions. We propose a framework for adjusting IDF curves and design storms to future climate conditions using the TENAX model, a novel statistical approach that can provide future short-duration precipitation return levels based on projected temperature changes. For most applications, information from climate models at the daily scale can be used to construct design storms at the sub-hourly scale without any downscaling or bias adjustment. Our approach is illustrated through a re-parameterization of the Chicago Design Storm (CDS) in the context of climate change. As a case study demonstration, we apply the TENAX model to data from the city of Zurich to calculate changes in the historical IDF curve for durations ranging from 10 min to 3 h. We then construct synthetic 100-year return period design storms based on the CDS for present and future climates and use the CAFlood model to produce flood inundation maps to assess changes in flood hazard. The codes for adapting design storms to climate change are simple to implement, easily applicable by practitioners, and made freely available.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104823"},"PeriodicalIF":4.0,"publicationDate":"2024-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0309170824002100/pdfft?md5=bcb162ebb1fe34867d9140461e25ffc6&pid=1-s2.0-S0309170824002100-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142275670","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}