Consecutive oscillatory eruptions of a mixture of gas and liquid in urban stormwater systems, commonly referred to as sewer geysers, are investigated using transient three-dimensional (3D) computational fluid dynamics (CFD) models. This study provides a detailed mechanistic understanding of geyser formation under partially filled dropshaft conditions, an area not previously explored in depth. The maximum geyser eruption velocities were observed to reach 14.58 m/s under fully filled initial conditions (hw/hd = 1) and reduced to 5.17 m/s and 3.02 m/s for partially filled conditions (hw/hd = 0.5 and 0.23, respectively). The pressure gradients along the horizontal pipe drove slug formation and correlated directly with the air ingress rates and dropshaft configurations. The influence of the dropshaft diameter was also assessed, showing a 116% increase in eruption velocity when the dropshaft to horizontal pipe diameter ratio (Dd/Dt) was reduced from 1.0 to 0.5. It was found that the strength of the geyser (as represented by the eruption velocity from the top of the dropshaft) increased with an increase in the initial water depth in the dropshaft and a reduction in the dropshaft diameter. Additionally, the Kelvin–Helmholtz instability criteria were satisfied during transitions from stratified to slug flow, and they were responsible for the jump and transition of the flow during the initial rise and fallback of the water in the dropshaft. The present study shows that, under an initially lower water depth in the dropshaft, immediate spillage is not guaranteed. However, the subsequent mixing of air from the horizontal pipe generated a less dense mixture, causing a change in pressure distribution along the tunnel, which drove the entire geyser mechanism. This study underscores the critical role of the initial conditions and geometric parameters in influencing geyser dynamics, offering practical guidelines for urban drainage infrastructure.
