In this study, an improved multi-scale algorithm was developed to analyze the effects of temperature on the bioclogging processes at pore and Representative Elementary Volume (REV) scales. In this algorithm, Immersed Boundary-Lattice Boltzmann Method- Cellular Automata (IB-LBM-CA) model for pore-scale simulation and Discrete Unified Gas-Kinetic Scheme-Cellular Automata (DUGKS-CA) model for REV simulation

Hydrodynamic simulation based on shallow water wave equations (SWE) is one of the most useful methods for flood routing analysis. However, its widespread application in large-scale basin flood disaster prevention is hindered by significant computational efficiency challenges. Addressing this issue, a hydrodynamic algorithm based on multi-GPUs is provided, and the program named as CoSim-SWE is developed

Understanding CO2 capillary trapping and dissolution trapping behaviors in deep saline aquifers is essential for improvement of sequestration efficiency. This study investigated the impact of rock wettability on CO2 capillary and dissolution trapping in a Ketton carbonate rock through direct numerical simulation. Based on experimentally measured CO2-brine-rock physics data, we performed the pore-scale

We present a pioneering experimental study of stratified, rotating exchange flows interacting with a bottom, mobile sediment bed that simulates large estuaries. Two-dimensional velocity fields are coupled with bed scan that allows to reconstruct the bed morphology. The experiments span a large parameter range, notably laminar to turbulent Ekman layer regimes (33

Understanding the drivers of drought variability is crucial for developing effective adaptation and management strategies. This study develops a two-step modelling approach to characterize and predict hydrological droughts in a nonstationary context. First, a multivariate Hidden Markov Model (HMM) is used to classify low-water level time series into Dry, Normal, and Wet years, identifying Dry years

This experimental study systematically investigates the influence of the Peclet number (Pe) and transport length on the transition to Fickian transport in homogeneous sand packs. Through Darcy column experiments with varying lengths and two distinct sediment sizes (d50), we analyzed breakthrough curves (BTCs) to quantify non-Fickian characteristics and transport parameters. The dispersion coefficient

A falaj (plural aflaj) (also known as qanat) is a type of horizontal or nearly horizontal well that extracts water from an aquifer by gravity. Despite their importance, aflaj with arbitrary trajectories has not yet been mathematically modeled. Moreover, the available analytical models do not include the effects of several head losses, including head loss inside the falaj due to friction with the inner

Systems involving advection-dominated transport through heterogeneous porous and fractured media are ubiquitous in subsurface engineering applications. However, upscaling such systems continues to challenge rigorous modeling efforts, particularly when advection is stronger than diffusion at fine spatial scales (i.e., when the Péclet number is greater than one at length scales that characterize a system’s

The insecticide λ-cyhalothrin a type II synthetic parathyroid. It is a hydrophobic and highly effective broad-spectrum insecticide commonly used in pest management. However, its presence in subsurface formations poses environmental toxicity risks and potentially may have adverse effects on humans. Chitosan, a polymer with unique physicochemical and absorption properties, is utilized in numerous various

In surf/swash sediment transport in coastal processes, the influence of pore flow near the surface of the sediment bed cannot be ignored. Understanding how the pore flow affects the turbulent structure of the oscillating boundary layer is crucial for comprehending the mechanisms of bottom sediment transport. However, measurement is difficult and the turbulent structure in the porous sublayer is not

It is believed that channel networks across fluvial landscapes are self-organized into fractal patterns in order to minimize energy expenditure, as evidenced by the similarities between computer-generated optimal channel networks (OCNs) and real networks. However, the specific mechanisms driving such energy minimization remain largely elusive. Here, we propose that the energy minimization tendency

CO2 diffusive leakage, or diffusive transport, through intrinsic or induced caprock fractures poses a significant concern for the security of CO2 sequestration in saline aquifers. Although this issue has garnered considerable interest and has been the subject of many numerical analyses, experimental studies remain limited. We present the first experimental investigation of CO2 diffusive leakage through

The subsurface environment’s complex heterogeneous structure poses challenges for accurately modeling transport phenomena due to limited data and measurement errors, leading to uncertainties in solute transport predictions. This study proposes a computational framework to semi-analytically compute the cumulative distribution function (CDF) of solute concentration in heterogeneous aquifers. We investigate

Solute dispersion in the aqueous phase within unsaturated porous media has critical implications to various natural and engineered systems, such as nutrient and contaminant transport in the vadose zone. This work developed pore-scale lattice Boltzmann (LB) modeling to simulate water-air multiphase flow and solute transport in the aqueous phase, which unraveled the role of the water saturation and Peclet

The macroscopic behavior of multiphase flow systems is governed by interfacial dynamics, which are a strong function of the synergy between viscous forces and wettability. However, the effects of wettability heterogeneity, an inherent feature of natural porous media, remain poorly understood. In this study, we incorporated spatial-related wettability into phase field model and systematically examined

Characterizing velocity and pressure fields in aquatic systems is crucial for understanding fundamental processes such as sediment transport, hyporheic flow, air-water exchange, and habitat quality. While large obstacles like vegetation stalks are known to create significant localized pressure gradients, the role of small-scale bed roughness in generating local pressure gradients remains poorly understood

This study investigates the interplay among water saturation, hydraulic conductivity, and mechanical dispersion in heterogeneous porous media at the continuum scale. Mechanical dispersion of dissolved chemical tracers is influenced directly by water velocity variability, which is governed by the porous structure, the distribution of the water phase within it, and its corresponding conductivity field

The pore water retained in unsaturated soil includes film water attached on the solid surface and capillary water in corners or pores, which are mainly controlled by adsorptive force from the solid surface and capillary force from the water-gas interface, respectively. The soil water retention (SWR) curve represents the fundamental characteristic of unsaturated soil, in which the connected capillary

In many environmental flow situations, a solid body emerges from a sediment bed. This may occur in natural systems, for example when a tree-root emerges from the bed, or around a man-made structure such as a bridge pier or a wind turbine foundation. When this situation occurs, various flow-structure interactions such as the horseshoe vortex or vortex shedding lead to the scour process, namely the erosion

The current study is dedicated to the formal derivation of a hierarchic of asymptotic models that approximate the groundwater waves problem within the Dupuit–Forchheimer regime, over a regular, non-planar substratum. The derivation methodology employed bears resemblance to the techniques utilized in hierarchic of asymptotic models for approximating the water waves problem in the shallow water regime

Geological carbon storage (GCS) is critical for sequestering CO2 deep underground. GCS projects may face environmental challenges, such as leakage risks, adverse pressure buildup, and groundwater contamination. Numerical simulators play a vital role in accurate forecasting but can be computationally expensive. In this work, we leveraged an updated Fourier Neural Operator (FNO) which includes data sparsity

Female salmonids bury their eggs in streambed gravel by digging a pit where they lay their eggs, which they then cover with gravel from a second pit, forming a rough-surfaced dune-like structure called a redd. The interaction between a redd and the stream flow induces surface water to flow into the sediment, through egg pockets, and reemerge downstream of the crest. These downwelling and upwelling

Seawater intrusion can be mitigated by extracting saltwater, creating a negative hydraulic barrier that reduces the extent of saltwater in coastal aquifers. The effects of seawater extraction are analyzed in the current study through a semi-analytical methodology based on sharp-interface, steady-state conditions. The methodology is based on the Dupuit-Forchheimer approximation and applies a power series

We performed two-dimensional (2D) pore-scale simulations of primary CO2 injection using a weakly compressible scheme for geological carbon sequestration (GCS) applications. The aim was to analyze pore-scale relative permeability and saturation of CO2 under wide-ranging injection velocities and wettabilities. The results show that saturation is highest for viscous fingering, lowest for crossover (−

Colloidal-suspension-nano flows with varying ionic strength are widely present in nature and industry. The variation of brine salinity, which highly affects electrostatic particle-rock interaction, triggers fines detachment and consequent rock alteration. The microscale models for fines detachment at the pore-particle and at rock-reservoir scales are widely used to predict core and field behaviour

Understanding effective permeability is crucial for predicting fluid migration and trapping in subsurface reservoirs. The Bunter Sandstone of northwestern Europe hosts major groundwater and geothermal resources and is targeted for CO2 storage projects. Here the effective permeability of fluvial facies within the Bunter Sandstone Formation was assessed using facies-scale models. Twelve lithofacies were

We develop a phase-field model for evaporation from a porous medium by explicitly considering a vapor component together with the liquid and gas phases in the system. The phase-field model consists of the conservation of mass (for phases and vapor component), momentum, and energy. In addition, the evolution of the phase field is described by the Allen–Cahn equation. In the limit of vanishing interface

In this study, based on the minimum operating power principle, the potential control mechanism of nonmonotonic wetting effects in porous media is analyzed. When different wetting conditions are applied to the system, the contribution weights of different energy contribution terms are different during the process when the system approaches the minimum operating power state. For weak drainage or weak

The presence of tides and subsurface dams adds complexity to the migration and mixing processes of land-sourced contaminant in coastal aquifers. While prior studies have explored individual effect of tides and subsurface dams, their combined impact on the transport characteristics of land-sourced contaminant remains unclear. This study conducted laboratory experiments and numerical simulations to thoroughly

Estimations of multi-phase flow properties, mainly permeability, are crucial for several applications, such as CO2 sequestration, efficient oil and gas recovery, and groundwater contaminant treatment. Current methods for estimating the sub-core scale properties rely on numerical simulations, which can be time-consuming. A suitable alternative to numerical simulations is Deep Neural Networks (DNN),

A 3D macroscopic gradient-dynamics model is developed and applied to sandy soil in presence of exopolysaccharides (EPS), to mimic a soil influenced by root exuded mucilages (rhizospheric soil). Depending on water content, amphiphilic soil has a hydrophilic or hydrophobic behavior which impacts water transfer and retention. To model this saturation-dependent wettability, we propose a nonequilibrium

Coastal regions are vulnerable to flood risks due to the combined effects of storm surges, riverine flooding, upstream reservoir releases, and inland rainfall. Traditional models often fail to integrate these critical factors, leading to inaccuracies in flood extent forecasting. This study addresses this gap by developing a comprehensive coastal flood inundation forecasting framework for a region impacted

We present an experimental technique for determining the pH and the total carbon concentration when CO2 diffuses and flows in water. The technique employs three different pH indicators, which, when combined with an image analysis technique, provides a dynamic range in pH from 4.0 to 9.5. In contrast to usual techniques in which a single pH indicator is used, the methodology presented allows not only

Drainage of a liquid by a gas in porous media can be broken down into two main mechanisms: a primary piston-like displacement of the interfaces through the bulk of pore bodies and throats, and a secondary slow flow through corners and films in the wake of the invasion front. In granular porous media, this secondary drainage mechanism unfolds in connected pathways of pendular structures, such as capillary

In contrast to homogeneous soil deposits, stratified layering introduces vertical heterogeneity, resulting in not only greater spatial variability but also more complex structural responses. This complexity is further exacerbated by gravitational compaction, which gives rise to distinct fluid flow and solid deformation mechanics within each variably saturated layer and at the interfaces between layers—markedly

Kornev's (1935, see e.g. p.74, Fig. 34 - right panel) “open system” of capillarity-driven wetting of a fine-textured soil from a buried ditch filled by a coarse porous material is modeled analytically, using the methods of hodograph, and numerically, with the help of HYDRUS2D. Gravity, Darcian resistance of the soil at full saturation but negative pressure, and capillarity are three physical competing

CO2 injection into subsurface reservoirs provides a long-term solution to anthropogenic emissions. A variable injection rate (such as ramping the flow rate up or down) provides flexibility to injection sites, and could influence the amount of residual trapping. Observations made in cm-scale samples showed that starting at a low flow rate established a flow pathway across the core at a low capillary

We present a new approach to the discretisation and solution of the Shallow Water Equations (SWE) based on the finite element (FE) method. The discretisation is expressed as the convolutional layer of a neural network whose weights are determined by integrals of the FE basis functions. The resulting system can be solved with explicit or implicit methods. Expressing and solving discretised systems with

Aquitards, which widely occur throughout sedimentary basins or alluvial plains, play important roles in groundwater storage and contaminant transport. In this study, a mathematical model for one-dimensional contaminant transport which considered the combined effects of diffusion, adsorption and nonlinear consolidation deformation processes in an aquitard (NCD model) was formulated. An analytical solution

When a partially miscible fluid dissolves into a host phase below it, buoyancy-driven fingering develops when the diffusive boundary solution created is denser than the underlying solvent. In many situations, the interface between the two fluids may present level variations introduced by geometrical irregularities. We study here numerically the influence of this interface undulation on the properties

Hydraulic barriers are useful for manipulating groundwater flow to mitigate and contain harmful environmental effects. Injecting air into the aquifer has been suggested as a cost-efficient, sustainable, and reversible hydraulic barrier. In this, the injected air reduces the conductance of the aquifer to water flow. However, this practice is not commonly used despite its potential, probably due to a

Models of colloid transport in porous media that assume constant fractional loss per grain passed fail in the presence of repulsive barriers to attachment, under which condition experiments produce profiles of colloid concentrations with distance from source that are nonexponential. Nonexponential removal is hypothesized to arise from variable likelihood of encountering nanoscale regions of attraction

When performing upscaling of transport phenomena in multiscale systems it is not uncommon that terms of different physical nature than those present at the underlying scale arise in the resulting averaged differential equations. For diffusive species mass transfer with heterogeneous reaction and conductive heat transfer, additional terms result from upscaling using the volume averaging method, which

Dissolution in porous media is widespread in natural and engineered systems, accompanied by the evolution of geometric structure, permeability and surface area of the porous matrix. Although extensive research has examined dissolution dynamics in porous media, there is a lack of quantitative characterization of the relationships among permeability, surface area and porosity, which depend on dissolution

Alluvial fan aquifers are among the most productive aquifers around the world. At present, in mathematical models of these types of aquifers, the unsaturated flow is usually neglected, and the instantaneous or delayed drainage boundary conditions are considered at the water table. The goal of this work is to present an analytical model of the drawdown due to a partially penetrating well in an unconfined

In this paper, we develop a well-balanced, conservative, high-order finite difference weighted essentially non-oscillatory (WENO) method for the shallow water equations. Our approach exactly preserves the moving-water equilibria of the shallow water equations with non-flat bottom topography. The proposed method consists of two key components. First, we reformulate the source term into a flux-gradient

The optimization of well controls over time constitutes an essential step in the design of cost-effective and safe geological carbon sequestration (GCS) projects. However, the computational expense of these optimization problems, due to the extensive number of simulation evaluations, presents significant challenges for real-time decision-making. In this paper, we propose a hybrid CNN-Transformer surrogate

Simulating subsurface contaminant transport at the kilometer-scale often entails modeling reactive flow and transport within and through complex geologic structures. These structures are typically meshed by hand and as a result geologic structure is usually represented by one or a few deterministically generated geological models for uncertainty studies of flow and transport in the subsurface. Uncertainty

In this study, we present a general random walk model for upscaling colloid transport and retention in two-dimensional porous media. The model is based on direct sampling from spatial and temporal jump distributions of single-interceptors, colloids that first enter the near-surface zone within 200 nm of a collector surface, derived from mechanistic pore-assembly trajectory simulations. This framework

Electrokinetic (EK) technologies are promising solutions for the remediation of contaminated sites, particularly in low-permeability porous media. However, their widespread adoption is hindered by the challenge of predicting the complex, coupled physico-chemical processes triggered by the application of external electric fields in the subsurface. Numerical models therefore represent essential tools

Foam flow in porous media is important in various engineering applications, including soil remediation, carbon dioxide sequestration, and enhanced oil recovery. This study explores the relationship between bubble density and permeability in foam flow models, focusing on how different approaches capture foam formation in highly permeable regions. We compare two mechanistic models numerically. The first

Understanding water movement in catchments subsurface is crucial for numerous applications such as pollutant contamination, nutrient loss, water resource management and ecosystem functioning. Among the variables of particular interest, the transit times of water particles and their statistical distribution are a desirable output. Nevertheless, past approaches assume explicitly the form of the transit

The use of geoelectrical monitoring of groundwater quality and contamination is a growing and promising topic. Nowadays, geoelectrical methods are mostly used as qualitative detection tools. This study aims to better use geoelectrical signals as a complementary tool for the quantitative characterization of chemical species transport and reaction in the porous matrix by developing a coupled mechanistic

Subterranean Estuaries (STEs) are important biogeochemical reactors at the land-ocean interface. They transform dissolved species prior to discharge, thereby influencing chemical fluxes from land to sea. The coupling between physical flow and biogeochemical reactions in the STE is complex, and a deeper process understanding demands the application of reactive transport modeling (RTM). Most previous

Porous media is widely distributed in nature, found in environments such as soil, rock formations, and plant tissues, and is crucial in applications like subsurface oil and gas extraction, medical drug delivery, and filtration systems. Understanding the properties of porous media, such as the permeability and formation factor, is crucial for comprehending the physics of fluid flow within them. We present

After sequestration of CO2 into subsurface saline aquifer, CO2 cap forms at the top of a stratum. As overlaying CO2 dissolves into brine, precipitation reactions between CO2 and in-situ ions emerge. The role of reactive precipitation during this process has long been under debate due to the lack of direct observation. Here we conduct visualized experiments on high-pressure CO2 convective dissolution

Biofilms and their formation dynamics are ubiquitous and complex in porous media. The mechanism of biofilm formation on solute transport behavior remains limited, which inhibits potential biofilm applications such as bioremediation. In this study, we present a new numerical solver, BioFOAM, based on the micro-continuum theory, to simulate the coupled pore-scale processes of biofilm formation, fluid