This work reviews the transient two-phase flow in porous media with engineering applications in Geotechnics. It initially overviews constitutive relationships, conventional theories, and experiments. Then, corresponding limitations are discussed according to conflicting observations and multiphase interfacial dynamics. Based on those findings, the dynamic nonequilibrium effects were so defined, which could also be abbreviated as dynamic/transient effects. Four advanced theories have already been developed to resolve these effects. This review collects them and discusses their pros and cons. In addition, this work further reviews the state-of-art in terms of experimental methods, influential factors in dynamic/transient effects, and modelling performance, as well as micromodel and numerical methods at pore-scale. Last, the corresponding geotechnical applications are reviewed, discussing their applicability in effective stress, shear strength, and deformation. Finally, the entire review is briefed to identify research gaps in Geotechnics.

Tidal wetland ecosystems are sensitive to porewater salinity dynamics. However, it is unclear how salts move and distribute in these wetlands, particularly how the salts accumulated by evaporation get removed from the wetland soil, so that the salinity levels may stabilize to accommodate vegetation. We conducted a combined field and modeling study to identify the porewater flow and salinity patterns in a subtropical wetland subjected to tidal inundation and evaporation. Measured and simulated salinity contours indicated the formation of hypersaline porewater plumes in the upper intertidal zone and the fresher porewater zones with salinity close to that of seawater near the creek and in the supratidal zone. Simulations indicated the discharge of the hypersaline upper intertidal porewater to the creek with a discharge pathway developed under the fresher near‐creek porewater zone, and this was further confirmed by the field‐observed significant salinity gradient under the creek. Our model suggested that both water and salt discharge from the wetland soil occurred predominantly at the creek bank and creek bed. The porewater discharge is more intensive through the creek bank than the creek bed, while the salt discharge across both the creek bank and creek bed was comparable due to the much higher salinity level under the creek bed. Salt discharge driven by density gradients and tidal‐induced porewater circulation provides a mechanism for removing salts accumulated in the upper intertidal zone due to evaporation and could prevent salt flat formation and marsh plants dieback.

The dynamic two-phase flow in porous media was theoretically developed based on mass, momentum conservation, and fundamental constitutive relationships for simulating immiscible fluid-fluid retention behavior and seepage in the natural geomaterial. The simulation of transient two-phase flow seepage is, therefore, dependent on both the hydraulic boundaries applied and the immiscible fluid-fluid retention behavior experimentally measured. Many previous studies manifested the velocity-dependent capillary pressure–saturation relationship (P_c-S) and relative permeability (K_r-S). However, those works were experimentally conducted on a continuum scale. To discover the dynamic effects from the microscale, the Computational Fluid Dynamic (CFD) is usually adopted as a novel method. Compared to the conventional CFD methods solving Naiver–Stokes (NS) equations incorporated with the fluid phase separation schemes, the two-phase Lattice Boltzmann Method (LBM) can generate the immiscible fluid-fluid interface using the fluid-fluid/solid interactions at a microscale. Therefore, the Shan–Chen multiphase multicomponent LBM was conducted in this study to simulate the transient two-phase flow in porous media. The simulation outputs demonstrate a preferential flow path in porous media after the non-wetting phase fluid is injected until, finally, the void space is fully occupied by the non-wetting phase fluid. In addition, the inter-relationships for each pair of continuum state variables for a Representative Elementary Volume (REV) of porous media were analyzed for further exploring the dynamic nonequilibrium effects. On one hand, the simulating outcomes reconfirmed previous findings that the dynamic effects are dependent on both the transient seepage velocity and interfacial area dynamics. Nevertheless, in comparison to many previous experimental studies showing the various distances between the parallelly dynamic and static P_c-S relationships by applying various constant flux boundary conditions, this study is the first contribution showing the P_c-S striking into the nonequilibrium condition to yield dynamic nonequilibrium effects and finally returning to the equilibrium static P_c-S by applying various pressure boundary conditions. On the other hand, the flow regimes and relative permeability were discussed with this simulating results in regards to the appropriateness of neglecting inertial effects (both accelerating and convective) in multiphase hydrodynamics for a highly pervious porous media. Based on those research findings, the two-phase LBM can be demonstrated to be a powerful tool for investigating dynamic nonequilibrium effects for transient multiphase flow in porous media from the microscale to the REV scale. Finally, future investigations were proposed with discussions on the limitations of this numerical modeling method.

The strength of unsaturated soil is defined by the soil water retention behavior and soil suction acting inside the soil matrix. In order to obtain the suction and moisture profile in the vadose zone, specific measuring techniques are needed. Time domain reflectometry (TDR) conventionally measures moisture at individual points only. Therefore, spatial time domain reflectometry (spatial TDR) was developed for characterizing the moisture content profile along the unsaturated soil strata. This paper introduces an experimental set-up used for measuring dynamic moisture profiles with high spatial and temporal resolution. The moisture measurement method is based on inverse modeling the telegraph equation with a capacitance model of soil/sensor environment using an optimization technique. With the addition of point-wise soil suction measurement using tensiometers, the soil water retention curve (SWRC) can be derived in the transient flow condition instead of the static or steady-state condition usually applied for conventional testing methodologies. The experiment was successfully set up and conducted with thorough validations to demonstrate the functionalities in terms of detecting dynamic moisture profiles, dynamic soil suction, and outflow seepage flux under transient flow condition. Furthermore, some TDR measurements are presented with a discussion referring to the inverse analysis of TDR traces for extracting the dielectric properties of soil. The detected static SWRC is finally compared to the static SWRC measured by the conventional method. The preliminary outcomes underpin the success of applying the spatial TDR technique and also demonstrate several advantages of this platform for investigating the unsaturated soil seepage issue under transient flow conditions.