Pore-scale modeling, aided by imaging advancements and computational power, is now a vital tool for comprehending fluid flow and transport in porous media. It allows detailed exploration of pore structures and analysis of fluid dynamics and mass transport at smallest scales. This work presents the Pore Region Growing (PREGO) Algorithm, a novel approach to address segmentation challenges encountered

Pore-scale modelling is an important tool to improve our understanding of multiphase flow in porous media. Slow fluid invasion is commonly modelled using quasi-static pore network models (PNM). These models simulate the invasion in a network of simplified pores and throats by invading network elements in order of the quasi-static “invasion” capillary pressure needed for the invading fluid to enter

This paper focus on the numerical approximation of two-layer shallow water system. First, a new approximation of the eigenvalues of the system is proposed, which satisfies some interesting properties. From this approximation, we give an accurate estimation of the non-hyperbolic region, which improves significantly the one computed with the classic eigenvalues approximation. In particular, we estimate

Droughts in Africa have a disproportionate detrimental influence on food security, water supply, social equilibrium and the economy compared to other regions in the World. Crop growth is principally rain-fed with highly seasonal rainfall regimes over most of the Continent. Hydraulic infrastructure and harvest storage are scarce. Drought is primarily initiated by a lack of precipitation, and it spreads

Pore-scale modeling has limited applicability at large scales due to its high computational cost. One common approach to upscale pore-scale models is the use of effective medium theories, which homogenize small-scale features in a porous structure and characterize the medium by macroscale properties (e.g., permeability) and equations (e.g., Darcy’s law). However, there are classes of physical processes

We evaluate masked block-structured grids for ocean domains which allow to represent small-scale geometric features without resorting to very small blocks or excessive mesh resolution. The considered approach aims to combine the geometric flexibility of unstructured meshes with the computational efficiency of stencil-based discretizations and is implemented and tested in a quadrature-free discontinuous

In this study, a new method of extraction of pore networks from tomography images of geological porous material has been introduced. Superpixels, a classic method of image segmentation, has been employed as a part of the proposed algorithm to analyze 3D micro-structures and construct a pore network model that is representative of the hydraulic and electrical behavior of the studied porous mediums.

Understanding hydrological connectivity is essential to investigate ecological processes in river catchments and floodplains. Assessing flooding behavior, including flooded areas and connection times, is required to analyze hydrological connectivity in river floodplains. Deep learning, especially Convolutional Neural Networks (CNNs), is an attractive alternative to hydrodynamic modeling, which is more

Geologic carbon dioxide (CO2) storage is a promising strategy for reducing its atmospheric levels, but concerns regarding potential leakage through faults or high-permeability inclusions in caprock highlight the importance of understanding its multiphase flow mechanism. In this study, a comprehensive experimental approach is introduced to characterize the pore structure of heterogeneous and fractured

A series of laboratory experiments was conducted to understand the role of unsteady coherent motions of turbulent flow in the onset of particle motion. Turbulent eddies having coherent structures constantly interact with sediment particles within the bottom boundary layer. Therefore, a 3D turbulent flow field and grain trajectories were obtained by the simultaneous operation of three-dimensional (3D)

Bubbles trapped inside porous materials occur in many applications ranging from geologic CO2 sequestration, underground hydrogen storage (UHS), groundwater remediation, to fuel cells and electrolyzers. If partially miscible in their surrounding wetting phase, bubbles can evolve through a process called Ostwald ripening, in which those with a high interfacial curvature dissolve fastest and feed into

Liquid-filled corners and capillary bridges can establish networks connecting seemingly isolated clusters during drainage in porous media. Coupled with drainage through the bulk of pores and throats, the flow through these networks constitutes a secondary drainage mechanism that can significantly affect fluid configurations and residual saturations. In order to investigate the prevalence of this drainage

Hydrologists have been actively exploring the utility of machine learning (ML) models for predicting streamflow. While ML methods have proven to be as accurate as conventional modeling techniques for streamflows well represented in the training set, they continue to lack satisfactory skills for extreme events. In this study, a novel ‘data reformation’ technique is proposed based on the Relative Strength

Viscous coupling effect on pore-scale flow capacity and its impact on macroscopic flow properties in porous media remain subjects of ongoing debate. This study employs the Stokes equation in conjunction with the zero-velocity interface and continuity interface to investigate viscous coupling effect during two-phase flow through a square capillary tube. Our findings substantiate the significance of

The release of anthropogenic chemicals to streams, stemming from contaminated sites, agriculture, urban sources, or accidental input, represents a significant threat to water resources and thus the health of humans and aquatic ecosystems. Predicting the transport and fate of chemicals is critical to quantifying contaminant concentrations and developing environmental quality standards (EQS). Tracer

In this paper, a numerical model is presented to simulate the density-driven flow in heterogeneous porous media with micro- and macro-fractures. The micro-fractures (like fissures) with relatively small fracture lengths compared to the reservoir are modeled using an equivalent continuum model. The macro-fractures (like faults) with large fracture lengths are modeled using the extended finite element

A set of analytical solutions have been presented for non-Darcian flow to a well in a leaky confined aquifer system considering leakage from the storage of semipervious layers/adjacent aquifer lying above and below the main aquifer. Flow in the main aquifer has been assumed to be non-Darcian and horizontal, whereas flow in the semipervious layers is assumed to be Darcian and vertical. The Izbash equation

This paper presents formulations for modeling the hydromechanical behavior in unsaturated soils using different coupling techniques. For this purpose, hydraulic constitutive relationships based on classical empirical models and mechanical constitutive relationships based on generalized effective stress are adopted. The coupling terms are defined using fixed stress rates. They are included in finite-element

In dilated fractures in the Earth’s crust fluid flow in combination with precipitation or dissolution processes can occur, which in turn influences the mechanical and transport properties of the rock system. We use a phase-field modeling framework to investigate these processes of epitaxial crystal growth and dissolution on the fracture walls of crystalline rock systems on microscale. Fluid flow simulations

Current implementations of Physics Informed Neural Networks (PINNs) can experience convergence problems in simulating fluid flow in porous media with highly heterogeneous domains. The objective of this work is to develop an accurate implementation of PINNs to model fluid flow in heterogeneous porous media that avoids alternative approaches in the literature that tend to impose unphysical continuity

In the present paper, we study the mathematical and numerical modelling of rapid transients at partially lifted sluice gates in the context of one-dimensional Shallow water Equations. First, improved exact solutions of the dam-break problem with partially lifted gate are presented, assuming (i) the dependence of the gate contraction coefficient on the upstream flow depth, and (ii) a viable definition

This paper presents a numerical investigation of multiple identical swash events to study the swash–swash interaction processes and their impacts on beachface evolution. The numerical model, based on the Nonlinear Shallow Water Equations, is first calibrated/validated against two different single-event-based data-sets. Multiple swash events are generated by identical solitary waves separated by different

Treated wastewater contains vital crop nutrients and offers the possibility for not only irrigation but also fertilization from a renewable source. The feasibility of replacing conventional chemical fertilizers with wastewater reuse is investigated, focusing on the implications of coupling irrigation and fertilization. We consider a multi-objective dynamic optimization problem of crop irrigation and

In this paper, we propose a novel full finite volume method to solve the advection–dispersion transport equation, where the pressure and concentration equations are solved implicitly and explicitly, respectively. The advective term is discretized using a high-order MUSCL-type finite volume method (Monotonic Upstream-centered Scheme for Conservation Laws) and to avoid numerical oscillations, we use

Recent studies argue that hydrological models must be configured for the specific objectives they need to address. Model configuration involves multiple model decision steps, such as choice of model structure, spatial discretization, spatial representation of forcing and performance metrics to calibrate them. We investigated the influence of each of the decisions in turn in a standardized framework

Data assimilation (DA) techniques such as the Ensemble Kalman filter (EnKF) and its extensions allow for real-time corrections of state-space models and model parameters based on an assumption of Gaussian error. The hydrological DA literature primarily documents applications of the EnKF to solve sequential state estimation problems. Recent advances in the DA literature demonstrate the potential of

Contaminant source and aquifer characterization (CSAC) is critical in groundwater pollution evaluation and remediation. The ensemble smoother with multiple data assimilation (ES-MDA) is utilized to jointly identify contaminant source information and hydraulic conductivities by assimilating time-lapse electrical resistivity tomography (ERT) data. In a synthetic profile with a time-varying release history

Fast and reliable estimation of engineered fracture geometries is a key factor in controlling undesirable fractures and enhancing stimulation design. Measuring the surface deformation gradient (tilt) for engineered fractures in shallow depths (<1000 m) has been proven a reliable source of data to infer fracture geometry, thanks to the impressive resolution of tiltmeter units (in the order of nano-radians)

Understanding the flow dynamics of hydrogen-water in aquifers is critical to maximizing hydrogen storage and recovery. By using direct computational fluid dynamics at elevated pressure, this paper aims at detailed pore-scale investigation of the effect of flow regime, compressibility, and hysteresis on flow pattern, trapping mechanisms and the efficiency of hydrogen storage and recovery. The results

Skeletonization, a crucial step in pore network modeling, traditionally involves the extraction of skeleton pixels from binarized, segmented X-ray images of porous materials. However, this conventional approach often suffers from user bias during segmentation, potentially leading to the loss of essential image details. This study addresses this limitation by developing deep learning model, called PoreSkel

Salt precipitation during brine evaporation in porous media is an important phenomenon in a variety of natural and engineering scenarios. This work establishes a multiphase multicomponent lattice Boltzmann (LB) method with phase change for simulating salt precipitation during brine evaporation. In the proposed LB models, the gas–brine multiphase flow, brine evaporation, salt concentration evolution

The hydrological cycle comprises several components, including precipitation, infiltration, evaporation, transpiration, and runoff, with the latter being a vital part of the cycle that connects precipitation to bodies of water. modelling is necessary due to limitations in standard hydrological measurement techniques, providing a better understanding of hydrological systems that were previously inaccessible

Previous studies of transverse hyporheic exchange in compound rivers mainly focused on processes driven by riverbed morphology, and neglected the effects of regional groundwater discharge (RGD). This study numerically investigated the transverse hyporheic exchange subjected to the joint influences of bedform and RGD in a conceptual compound channel. The results show the existence of an interfacial

Modeling dissolution processes in discrete fracture networks (DFNs) is a challenging task. Challenges are related to the highly nonlinear coupling between flow, mass transport, and reactive processes associated with fracture aperture evolution by dissolution. Further, advection-dominated transport due to fast fluid flow in fractures renders the problem more complex from a computational point of view

Ocean tides cause water table overheight near the seaward boundary of coastal aquifers which can have a large impact on groundwater flows and saltwater intrusion. Despite this, regional-scale coastal groundwater flow models often neglect the effects of ocean tides or alternatively assume that the influences of ocean tides and sea level are constant over time. These simplifications are often required

This paper expands the Consistent Particle Method (CPM) to simulation of turbidity currents with focus on turbulent water-sediment flow. A mixture model is adopted in solving the sediment transport equation to simulate the sediment mass transfer due to settling convection and turbulent diffusion. To compute the gravity-driven convection term, the proposed CPM-upwind scheme is relatively easy to implement

Surrogate modelling has been used successfully to alleviate the computational burden that results from high-fidelity numerical models of seawater intrusion in simulation-optimization routines. Nevertheless, little attention has been given to multi-fidelity modelling methods to address cases where only limited runs with computationally expensive seawater intrusion models are considered affordable imposing

The geomorphic relations of hydraulic geometry at-a-station and downstream, expressed as power laws, are derived for a straight, single thread, gravel river channel. The cross-section has a plane central bed, narrow or wide, flanked by curved sides or banks. The analysis employs mechanistic formulae for fixed and mobile bed flow resistance, downchannel bedslope, bedload transport and the novel introduction

In the context of a T-shaped open-channel confluence, where the tributary flow separates at the downstream junction edge and creates a separation zone in the post-confluence channel, this paper provides a critical reflection upon the terminology and the delineation methods for the separation zone and its proxies: the recirculation zone and the reverse flow zone. Following the fluid-mechanical literature

As a gridless computational method, physics-informed neural networks (PINNs) have been widely applied and are a research hotspot in the field of numerical simulation due to their potential to provide approximations that are infinitely close to the analytical solutions. To address the complex hydrodynamic simulation scenarios with a large solution domain, where the conventional PINNs are not fully applicable

The connectivity of the more conductive hydrofacies strongly determines flow and transport in heterogeneous media. Here we study solute transport in 3D binary isotropic samples with a proportion p of a high hydraulic conductivity facies (k+), and (1−p) of a low (k−) one. The k+ facies is characterized by two connectivity parameters: a connectivity structure type (no, low, intermediate and high), that

Even though the original paper (Townley, 1995) contained no significant mathematical errors, this Addendum has been prepared (i) to give readers access to 3D graphics that are of much higher quality than are available online, (ii) to provide access to Mathematica notebook files in both .nb and .pdf formats, that demonstrate the use of Mathematica to obtain the six published solutions and to produce

The direct acquisition of the permeability of porous media by digital images helps to enhance our understanding of and facilitate research into the problem of subsurface flow. A complex pore space makes the numerical simulation methods used to calculate the permeability quite time-consuming. Deep learning models represented by three-dimensional convolutional neural networks (3D CNNs), as a promising

Steady, phreatic flows in wedge-shaped aquifers subtended by a tilted bedrock and discharged through a hillslope seepage face are studied by 2-D and 1-D analytical models for fully saturated and a FEM numerical model for saturated-unsaturated flows. The curvilinear triangle in the hodograph domain is conformally mapped onto a reference half-plane where the Hilbert problem is solved for two holomorphic

Sediments with poor drainage conditions are prone to build-up of excess pore water pressure when exposed to dynamic wave loading. This results in accumulated seepage flow that may have a significant impact on migration of solutes in marine sediments. In this study, the complex interplay between solute transport and wave-induced seepage process is investigated through a pioneering approach that combines

We investigated the applicability of Silica encapsulated, superparamagnetic DNA particles (SiDNAmag) in determining aquifer hydraulic parameters at different ionic strengths (1 mM, 5 mM, and 20 mM phosphate buffer) of injection suspension. Thereto, in a homogeneous, unconsolidated sand tank we pulse - injected two uniquely sequenced SiDNAmag at two injection points. At 0.5 m and 0.8 m downstream from

This study develops two new analytical models for radially divergent tracer tests at a fully-penetrating well in two-zone confined aquifers of skin and formation zones. Mixing length in the well is below the well length for inducing vertical dispersion. One model treats skins as an adsorptive Robin boundary condition for reflecting the adsorption effect. The other model applies skin's governing equation

In this study, a novel algorithm coupling the flow at the pore and representative elementary volume (REV) scales was developed to analyze the effect of microbial growth at the pore scale on the permeability decay of porous media at the REV scale. This provides a new method for investigating microbial clogging in porous media. At the pore scale, an immersed boundary-lattice Boltzmann model was used

Shallow groundwater resources overlaying deep saline formations used in carbon storage applications are subjected to a potential contamination threat by CO2/brine leakage via natural or anthropogenically-induced conductive pathways in the confining caprock. Identifying the leakage source location and rate is critical for developing remediation plans and designing corrective actions. Owing to limited

We propose a novel hybrid-dimensional model for the Darcy–Forchheimer flow in fractured rigid porous media, with a natural applicability to non-conforming meshes. Motivated by the previous work on the reinterpreted discrete fracture model (RDFM) for Darcy flows, we extend its key idea to non-Darcy flows. Coupling the Darcy’s law in the matrix and the Forchheimer’s law in the fractures through the introduction

Multi-species diffusion/heterogeneous reaction coupled problem in porous media involving a boundary layer problem is studied. The spectral approach developed in our previous work allows to derive a macroscopic model in the high Damköhler number regime combining homogenization technique and spectral approach. Such a homogenized model cannot describe the perturbation at the external boundary where the

Decision support systems (DSSs) in water resources management often use optimization models to make decisions. Therefore, the quality of these models is crucial for the success of DSSs. However, optimization models can struggle with reaching global optima, resulting in inconsistent and untrustworthy results for decision-makers. Additionally, current solution methods can be slow and inefficient, hindering

A decrease in reservoir pressure can lead to remobilization of residually trapped CO2. In this study, the pore-scale processes related to trapped CO2 remobilization under pressure depletion were investigated with the use of high-resolution 3D X-ray microtomography. The distribution of CO2 in the pore space of Bentheimer sandstone was measured after waterflooding at a fluid pressure of 10 MPa, and then

In this paper we use a statistical framework to analyze the relation between storm properties and the statistics of extreme precipitation. We identify storm events using a 24-hour dry hiatus separation. We investigate the statistics of the hourly maximum intensity for the heaviest storm events at durations of 1 and 24 h. A two-parameter Weibull distribution is used to represent precipitation frequencies