Abstract Accurately simulating multiphase flows at low capillary numbers remains challenging due to spurious currents and interface discretization errors. We present a systematic comparison of four volume‐of‐fluid (VOF) solvers within the OpenFOAM framework for pressure‐driven Poiseuille–Bingham flow over grooved surfaces: interFoam with classical continuum, smoothed continuum, sharp surface force formulations, and the multiphaseInterFoam solver. Solver performance is evaluated under capillarity‐dominated conditions over a wide range of Bingham and Reynolds numbers and for both air‐ and oil‐filled grooves. Among all cases, multiphaseInterFoam consistently produces the lowest spurious currents. Its predictions are validated against experimental interface profiles and phase‐field simulations, demonstrating excellent agreement in both interface morphology and axial velocity. The validated solver is then used to examine the influence of geometric and viscoplastic parameters on flow over air‐filled textures. The results show that the velocity field within the groove is strongly affected by the Bingham number, particularly at high Reynolds numbers. The local slip length attains a maximum at the groove centre and decreases smoothly toward the edges, with its distribution governed by the groove aspect ratio, while the effective slip length increases monotonically with Bingham number, slip area fraction, and groove aspect ratio.
Joulaei et al. (Mon,) studied this question.