No abstract

No abstract

Dynamic pore-network model (PNM) has been widely used to model pore-scale two-phase flow. Numerical algorithms commonly used for dynamic PNM including IMPES (implicit pressure explicit saturation) and IMP-SIMS (implicit pressure semi-implicit saturation) can be numerically unstable or inaccurate for challenging flow regimes such as low capillary number (Ca) flow and unfavorable displacements. We perform comprehensive analyses of IMPES and IMP-SIMS for a wide range of flow regimes under drainage conditions and develop a novel fully implicit (FI) algorithm to address their limitations. Our simulations show the following: (1) While IMPES was reported to be numerically unstable for low Ca flow, using a smoothed local pore-body capillary pressure curve appears to produce stable simulations. (2) Due to an approximation for the capillary driving force, IMP-SIMS can deviate from quasi-static solutions at equilibrium states especially in heterogeneous networks. (3) Both IMPES and IMP-SIMS introduce mass conservation errors. The errors are small for networks with cubic pore bodies (less than 1.4% for IMPES and 1.2% for IMP-SIMS). They become much greater for networks with square-tube pore bodies (up to 45% for IMPES and 46% for IMP-SIMS). Conversely, the new FI algorithm is numerically stable and mass conservative regardless of the flow regimes and pore geometries. It also precisely recovers the quasi-static solutions at equilibrium states. The FI framework has been extended to include compressible two-phase flow, multicomponent transport, and phase change dynamics. Example simulations of two-phase displacements accounting for phase change show that evaporation and condensation can suppress fingering patterns generated during invasion.

Air–water interfacial adsorption complicates PFAS transport in vadose zones. Air–water interfaces can arise from pendular rings between soil grains and thin water films on grain surfaces, the latter of which account for over 90% of the total air–water interfaces for most field‐relevant conditions. However, whether all thin‐water‐film air–water interfaces are accessible by PFAS and how mass‐transfer limitations in thin water films control PFAS transport in soils remain unknown. We develop a pore‐scale model that represents both PFAS adsorption at bulk capillary and thin‐water‐film air–water interfaces and the mass‐transfer processes between bulk capillary water and thin water films (including advection, aqueous diffusion, and surface diffusion along air–water interfaces). We apply the pore‐scale model to a series of numerical experiments—constrained by experimentally‐determined hydraulic parameters and air–water interfacial area data sets—to examine the impact of thin‐water‐film mass‐transfer limitations in a sand medium. Our analyses suggest: (1) The mass‐transfer limitations between bulk capillary water and thin water films inside a pore are negligible due to surface diffusion. (2) However, strong mass‐transfer limitations arise in thin water films of pore clusters where pendular rings disconnect. The mass‐transfer limitations lead to early arrival and long tailing behaviors even if surface diffusion is present. (3) Despite the mass‐transfer limitations, all air–water interfaces in the thin water films were accessed by PFAS under the simulated conditions. These findings highlight the importance of incorporating the thin‐water‐film mass‐transfer limitations and surface diffusion for modeling PFAS transport in vadose zones.This article is protected by copyright. All rights reserved.

Increasing impervious land cover has great impacts on groundwater regimes in urbanized areas. Low impact development (LID) is generally regarded as a sustainable solution for groundwater conservation. However, the effects of LID on the spatial-temporal distribution of groundwater are not yet fully understood. In this case study, a coupled Storm Water Management Model (SWMM) and Finite Element Subsurface FLOW system (FEFLOW) model was used to simulate surface and groundwater flow in an urbanized area in Shenzhen, China. After verification, the model was used to analyze the spatial-seasonal variations of groundwater head and hydrological processes under different LID scenarios. The results indicate that if the runoff from 7.5% and 15% of impervious area is treated by LID facilities, the annual surface runoff decreases by 5% and 9%, respectively, and the spatial average groundwater head relative to sea level pressure increases by 0.9 m and 1.7 m in the study area, respectively. The rise in groundwater head generally decreases from the recharge zones to the discharge zones surrounded by the streams and coastal waters. However, the groundwater head change is determined not only by the location in the catchment, but also by the hydraulic conductivity of underlying aquifer and LID infiltration intensity. Moreover, LID significantly enhances groundwater recharge and aquifer storage in the wet seasons; in turn it increases aquifer release and groundwater discharge in the dry seasons. However, LID has the potential to increase the risk of groundwater flooding during wet seasons in areas with poor aquifer drainage capacity and shallow groundwater depth. The findings from this study provide the basis for further assessing the benefit and risk of LID infiltration for groundwater supplementation in the urbanized areas.