Massively parallel simulation of multiphase flows in porous media is challenging due to the highly nonlinear governing equations with the complexity of various modeling features and the significant heterogeneity of material coefficients. In this paper, we introduce a fully implicit discontinuous Galerkin (DG) reservoir simulator designed for parallel computers to handle large-scale multiphase flow

We present various analytical solutions for a two-dimensional steady air source that demonstrate classical flow physics’ usefulness in air-sparging evaluation. These solutions capture different flow geometry and setup effects while offering accuracy, speed, and a deeper understanding of governing physics. Compared to empirical models, they excel with fewer physical parameters. We validate their accuracy

Extreme rainfall events trigger natural hazards, including floods and debris flows, posing serious threats to society and the economy. Accurately quantifying extreme rainfall return levels in ungauged locations is crucial for improving flood protection infrastructure and mitigating water-related risks. This paper quantifies the added value of combining rainfall observations with Convection Permitting

Hydrological model calibration plays a crucial role in estimating optimal parameters for accurate simulation. Estimation of parameters is inevitable in hydrological modeling due to the challenge of directly measuring them, as most parameters are conceptual descriptions of physical processes. Modelers commonly employ optimization algorithms for calibrating hydrological models. However, these algorithms

Traditional approaches to mathematically describe spontaneous imbibition are usually based on either macro-scale models, such as Richards equation, or simplified pore-scale models, such as the bundle of capillary tubes (BCTM) or pore-network modeling (PNM). It is well known that such models cannot provide full microscopic details of the multiphase flow processes and that many pore-scale mechanisms

With the continuous improvement of mathematical modeling technology, reconstructing the three-dimensional structure of media from two-dimensional reference images has become an important research method for the three-dimensional modeling of multi-porous media. Deep-learning-based methods are currently popular and form the focus of this research field. However, the performance of deep learning in reconstructing

We develop an enhanced reduced model for single-phase flow in fractured porous media capable of incorporating more realistic interface conditions at the fracture terminations. In addition to the traditional dimensional model reduction, where the elements of the discrete fracture network are treated as lower dimensional manifolds embedded in the porous matrix, we explore the microscale behavior of the

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known

CO2 capture and storage is a viable solution in the effort to mitigate global climate change. Deep saline aquifers, in particular, have emerged as promising storage options, owing to their vast capacity and widespread distribution. However, the task of proficiently monitoring and simulating CO2 behavior within these formations poses significant challenges. To address this, we introduce the physics-constraint

Droughts are recurrent phenomena that present a large variety of space and time patterns making rather difficult the assessment of their rarity and the comparison between events.

Land-based pollutants threaten coastal aquifers, highlighting the need to protect groundwater and nearshore marine ecosystems. While aquifer heterogeneity has been recognized as a significant factor affecting solute behavior, the impact of fractures on land-sourced solute transport in coastal aquifers remains unclear. This study attempted to address this issue through laboratory experiments and discrete

The coupled model of water flow and solute transport in unsaturated soils is addressed in this study. Building upon previous research findings by Keita, Beljadid, and Bourgault, we investigate a class of second-order time-stepping techniques where two free parameters are introduced, to identify the most stable and accurate scheme. The spatial discretization of the Richards equation is accomplished

Advancements in satellite technology yield environmental data with ever improving spatial coverage and temporal resolution. This necessitates the development of techniques to discern actionable information from large amounts of such data. We explore the potential of dynamic mode decomposition (DMD) to discover the dynamics of spatially correlated structures present in global-scale data, specifically

History matching is crucial for reliable numerical simulation of geological carbon storage (GCS) in deep subsurface aquifers. This study focuses on inferring highly complex aquifer permeability fields with multi- and intra-facies heterogeneity to improve the characterization of CO2 plume migration. We propose a deep learning (DL)-based parameterization strategy combined with the ensemble smoother with

The sediment transport direction is affected by the bed slope. This effect is of crucial importance for two- and three-dimensional modelling of the interaction between the flow of water and the alluvial bed. It is not uncommon to find applications of numerical morphodynamic models in the literature that exaggerate the effects of transverse bed slopes on sediment transport compared to results from laboratory

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

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

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

Partially miscible bubbles (e.g., CO2) 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

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

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

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 (log10k) plays an important role in the onset, growth

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

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

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

Distributed recharge is commonly predicted from groundwater level data by adopting the water table fluctuation method (WTFM). The simplicity of the technique makes it attractive for groundwater management applications seeking sustainable levels of extraction. While there are variations to the WTFM, the classic approach extends the antecedent recession curve (prior to recharge events) to allow for the

The electrical conductivity of a porous medium is strongly controlled by the structure of the medium at the microscale as the pore configuration governs the distribution of the conductive fluid. The pore structure thus plays a key role since different geometries translate in variations of the fluid distribution, causing different behaviors measurable at the macroscale. In this study, we present a new

Hydraulic conductivity exhibits a high spatial variability due to the heterogeneity and discontinuity of the geologic environments and their constituent materials. Representing such variability is problematic when implementing groundwater flow models, especially in geological media such as fractured rocks, fractured porous media, and karstic media, where the scale of observation is important when defining

Simulating the transfer of mass between particles is not straightforwardly parallelized because it involves the calculation of the influence of many particles on each other. Engdahl et al. (2019) intuited that the number of matrix operations used for mass transfer grows quadratically with the number of particles, so that dividing the domain geometrically into sub-domains will give speed and memory

Coreflooding experiments are one of the primary methods for reservoir rock characterization and have been developing in recent years, providing increasing detail. An advanced analysis of coreflooding experiments consists of constructing three dimensional permeability (k) maps of the core sample with sub-core resolution. Such detailed characterizations provide important information on the core heterogeneity

A conceptual mathematical model was developed to simulate the transport of migrating nanoparticles in homogeneous, water saturated, 1-dimensional porous media. The model assumes that nanoparticles can collide with each other and aggregate. Nanoparticles can be found attached reversibly and/or irreversibly onto the solid matrix of the aquifer or suspended in aqueous phase. Attached particles may either

Many traditional models that predict plant–groundwater use based on groundwater level variations, such as the White method, make various simplifying assumptions. For example, these models often neglect the role of plant hydraulic redistribution, a process that can contribute up to 80% of transpiration. Thus, this work aims to avoid such assumptions and subsequently explore the dynamic interactions

Carbon capture and storage (CCS) in deep geological aquifers has shown to be the most viable option for mitigating the greenhouse gas effect of carbon dioxide (CO2) at a large scale. However, the underground formations often possess discontinuities in the caprocks, leaking the stored CO2. Potential leakage paths, such as abandoned wells, have been growing due to excessively unplugged oil and gas exploration

This study explores the suitability of quasi-static pore-network modeling for simulating the transport of hydrogen in networks with box-shaped pores and square cylinder throats. The dynamic pore-network modeling results are compared with quasi-static pore-network modeling, and a good agreement is observed when the simulations reach steady-state, for a capillary number of Nc≤10−7. This finding suggests

Digital rock analysis has become increasingly crucial in earth sciences and geological engineering. However, the multiscale characteristics of rock pores often exceed the capabilities of single-resolution imaging, which is inadequate for a comprehensive description of their characteristics. To address this issue, we introduce a novel multiscale rock image fusion framework based on a generative adversarial

We study the abilities of pore network models of different complexities to determine the flow statistics and structural controls for single-phase flow in partially saturated porous media. The medium permeability and hydraulic tortuosity are the basic parameters for upscaling flow problems from the pore to the Darcy scale. They represent average flow properties. However, upscaling and predicting dispersion

Dam-break waves are a major concern for communities and infrastructures in flood-prone areas. The impact of dam-break waves against rigid obstacles after propagation on a mobile bed is lacking both in experimental datasets and in numerical investigations aimed at assessing the capabilities and limitations of available morphodynamic models. To fill these gaps, a novel data set from experiments of dam-break

Simulating non-wetting fluid invasion in volumetric images of porous materials is of broad interest in applications as diverse as electrochemical devices and CO2 sequestration. Among available methods, image-based algorithms offer much lower computational cost compared to direct numerical simulations. Recent work has extended image-based methods to incorporate more physics such as gravity and volume-controlled

The dynamics and aggregation of microplastics in marine environments are investigated through high-fidelity direct numerical simulations with Lagrangian point-particle tracking. The properties of microplastics and biogenic particles, including size, density, and concentration, align with scenarios typical of seawater systems. The stickiness nature of microplastics, induced by biofilm formation (biofouling)

Physics-informed neural networks (PINNs) are receiving increased attention in modeling flow in porous media because they can surpass purely data-driven approaches. However, in heterogeneous domains, PINNs often face convergence challenges due to discontinuities in rock properties. A promising alternative is the mixed formulation of PINNs, which utilizes pressure head and velocity fields as primary

A sudden surge of data has created new challenges in water management, spanning quality control, assimilation, and analysis. Few approaches are available to integrate growing volumes of data into interpretable results. Process-based hydrologic models have not been designed to consume large amounts of data. Alternatively, new machine learning tools can automate data analysis and forecasting, but their

Freshwater lenses in otherwise saline environments contain an important source of fresh water for natural vegetation and agricultural crops. Such lenses are regularly found in areas where both upward seeping saline groundwater and downward infiltrating fresh recharge water occur simultaneously during part of the year, resulting in shallow freshwater lenses which are highly susceptible to changes in

The present study investigates the quality of four three-dimensional (3D) printing technologies to accurately reproduce the complex pore structure of a real undisturbed soil sample for laboratory experiments of transport in porous media at a 1:1 scale. Four state-of-the-art 3D printing technologies were evaluated (digital light synthesis, PolyJet with gel support material, low-force stereolithography

Thin water lenses floating on top of the main groundwater body are important for many natural and agricultural systems, owing to their different properties in terms of chemical composition or density compared to the surrounding groundwater. In settings with upward seeping groundwater, lenses may form that have thicknesses ranging from tens of centimeters to a few meters, making them prone to changing

Deep limestone aquifers are potential CO2 storage sites, but CO2-saturated brine reacts with the carbonate rock, changing its transport and storage properties. This study provided a preliminary investigation of the optimal injection rate of CO2-saturated brine in carbonate rocks. Indiana limestone cores were subjected to CO2-saturated brine injection at varied rates using an HPHT coreflooding setup

Mixing-limited reactions are central to a wide range of processes in natural and engineered porous media. Recent advances have shown that concentration gradients sustained by flow at the pore-scale influence macroscopic reaction rates over a large range of reactive transport regimes. Yet, resolving concentration gradients driven by fluid mixing at the pore-scale is challenging with current simulation

Non-Darcy flow through porous media is of great significance in hydraulics, oil and gas engineering, biomedical science, chemical and civil engineering etc. However, it is difficult to fully grasp the nature of fluid flow through porous media from macroscopic scale alone. Based on the statistically fractal scaling laws of pore structures, a new fractal pore-throat chain model (FPTCM) for non-Darcy

Increasing flood risk due to urbanization and climate change poses a significant challenge to societies at global scale. Hydrologic information that is required for understanding flood processes and for developing effective warning procedures is currently lacking in most parts of the world. Procedures that can combine global climate dataset from satellite and reanalysis with fast and low computational

We unravel the nonmodal stability of a three-dimensional nonstratified Poiseuille flow in a saturated hyperporous medium constrained by impermeable rigid parallel plates. The primary objective is to broaden the scope of previous studies that conducted modal stability analysis for two-dimensional disturbances. Here, we explore both temporal and spatial transient disturbance energy growths for three-dimensional

Pore networks play a key role in understanding and quantifying flow and transport processes in complex porous media. Realistic pore-spaces may be characterized by singular regions, that is, isolated subnetworks that do not connect inlet and outlet, resulting from unconnected porosity or multiphase configurations. The robust identification of these features is critical for the characterization of network

Hydraulic tomography (HT) has been proven effective for characterizing aquifer hydraulic heterogeneity for decades. Many have also proposed using electrical resistivity tomography (ERT) to gain prior information about large-scale layer structures to improve the HT estimates, when the number of pumping tests and drawdown measurements is limited. This study investigates the merits of such proposals via

The overall benefits of geological carbon storage (GCS) depend primarily on CO storability and injectability, expressed as saturation and relative permeability, respectively. The effects of GCS schemes on these two properties, the macroscopic response indicators of a two-phase seepage system, are closely related to pore-scale two-phase behaviors. However, the comprehensive effects of capillary number

Multispecies transport analytical models that solve advection-dispersion equations (ADEs) are efficient tools for evaluating the transport of decaying contaminants and their sequential products. This study develops a novel semi-analytical model to simulate the multispecies transport of decaying contaminants, considering nonequilibrium sorption and decay in both dissolved and sorbed phases. First-order

Numerical simulations explore for the first time the role of mega columns and aperture variability on particle transport through mature volumetric contraction networks as informed by a unique synthesis of network propagation and maturity. Columnar fracture patterns are generated by updating a series of Voronoi centers to the midpoint of a generated polygon over many iterations, creating 250 network

Modelling non-Darcy flow behaviour in shale rocks, composed of nanometer-sized pores and multi-scale fracture networks, is crucial for various subsurface energy applications. However, incorporating multiple physical mechanisms across numerous scales is not trivial. This work proposes an improved and practical upscaling workflow for coupling an Upscaled Discrete Fracture Matrix (UDFM) model and a pressure-dependent

Transport phenomena in unsaturated porous media still present an important research topic. In particular, in the context of recent environmental concerns, further understanding of contaminant transport in the partially saturated vadose zone is necessary. However, there is currently a lack of understanding of the relationship between water saturation, in particular the two-phase distribution, and dispersion

This paper focuses on the modeling of hydrogen (H2) storage in subsurface formations, particularly focusing on the equilibrium between H2 and brine and its implications for hydrogen transport properties in black-oil reservoir simulations. Initially, we evaluate and calibrate various equations of state (EoS) for H2-water and H2-brine mixtures. Our analysis ranges from the molecular-level Perturbed-Chain

Solute transport in unsaturated porous media is of interest in many engineering and environmental applications. The interplay between small-scale, local forces and the porous microstructure exerts a strong control on the transport of fluids and solutes at the larger, macroscopic scales. Heterogeneity in pore geometry is intrinsic to natural materials across a large range of scales. This multiscale

There are still many unclear mechanisms in the multiphase reactive flow with solid dissolution processes. In this study, the reactive transport processes coupled with solid dissolution and self-induced multiphase flow in three-dimensional (3D) structures with increasing complexity is studied by developing a 3D computational microfluidic method, which considers multiphase flow, interfacial mass transport