Pipeline leaks present significant challenges in civil engineering, including water loss, soil erosion, and the formation of settlements or sinkholes. While experimental studies have addressed this phenomenon, numerical modeling remains underexplored, particularly regarding the effects of leaks in different directions and burial depths. This research employs a coupled Computational Fluid Dynamics (CFD) and Discrete Element Method (DEM) framework to investigate the impact of pipeline leaks on sandy soil. The model consists of a rectangular container filled with particles, with a defective pipeline buried at varying depths. The analysis evaluates the damage potential as a function of leak area, direction, flow rate, burial depth, and particle size. Results show that water streamlines curve upwards, dragging particles and forming a protruding mound on the surface. Above the leak, particle restitution creates a cyclic vertical circulation pattern, leading to surface settlement at the leak origin. Analytical models predicting pressure loss through the soil are also assessed based on leak configuration. These findings provide valuable insights into the failure mechanism of sandy soils surrounding pipeline leaks, contributing to improved risk assessment and mitigation strategies.