Abstract A mathematical model for slug (finite liquid volume) motion in not-fully-wettable capillary tubes with sinusoidally varying cross-sectional areas was developed. The model, based on the Navier-Stokes equation, accounts for the full viscous terms due to nonuniform geometry, the inertial term, the slug's front and rear meniscus hysteresis effect, and dependence of contact angle on flow velocity (dynamic contact angle). The model includes a velocity-dependent film that is left behind the advancing slug, reducing its mass. The model was successfully verified experimentally by recording slug movement in uniform and sinusoidal capillary tubes with a gray-scale high-speed camera. Simulation showed that tube nonuniformity has a substantial effect on slug flow pattern: in a uniform tube it is monotonic and depends mainly on the slug's momentary mass/length; an undulating tube radius results in nonmonotonic flow characteristics. The static nonzero contact angle varies locally in nonuniform tubes owing to the additional effect of wall slope. Moreover, the nonuniform cross-sectional area induces slug acceleration, deceleration, blockage, and metastable-equilibrium locations. Increasing contact angle further amplifies the geometry effect on slug propagation. The developed model provides a modified means of emulating slug flow in differently wettable porous media for intermittent inlet water supply (e.g., raindrops on the soil surface).