Coupled two-phase flow and surfactant/PFAS transport in porous media with angular pores: From pore-scale physics to Darcy-scale modeling

Two-phase surfactant-laden flow and transport in porous media are central to many natural and engineering applications. Surfactants alter two-phase flow by modifying interfacial tension and wettability, while two-phase flow controls surfactant transport pathways and interfacial adsorption. These coupled processes are commonly modeled using Darcy-type two-phase flow equations combined with advection--dispersion--adsorption transport equations, with capillary pressure--saturation relationships scaled by the Leverett $J$-function. However, the Leverett $J$-function idealizes porous media as bundles of cylindrical tubes and decouples interfacial tension and wettability, limiting its ability to represent angular pore geometries and interfacial tension--wettability coupling effects. We present a modeling framework that explicitly incorporates pore angularity and interfacial tension--wettability coupling into Darcy-scale surfactant-laden flow and transport models. Two-phase flow properties are derived for angular pores, upscaled across pore size distributions, and formulated as explicit and closed-form expressions. These upscaled relationships are integrated into a coupled flow--transport model to simulate transient two-phase flow and surfactant transport. Results reveal a nonlinear and nonmonotonic dependence of two-phase flow properties on pore angularity, pore size distribution, and interfacial tension. Example simulations of water flow and PFAS migration in unsaturated soils indicate that surfactant-induced flow effects on PFAS leaching are generally minor under typical conditions, whereas pore angularity strongly controls water flow, interfacial area, and PFAS retention. Overall, the proposed framework provides a more physically grounded approach for modeling two-phase surfactant-laden flow and transport in porous media.