Pipeline leaks pose significant challenges in civil engineering due to water loss, the potential for soil erosion, and the consequent settlement or sinkhole formation. While experimental studies have explored this phenomenon, numerical modeling remains less emphasized, particularly in analyzing the effects of leaks in different directions. This research employs numerical simulations to investigate the impact of pipeline leaks on surrounding soil under various flow directions. The study utilizes a coupled Computational Fluid Dynamics (CFD) and Discrete Element Method (DEM) tool, CFDEM, which integrates OpenFOAM for the Eulerian fluid phase and LIGGGHTS for the Lagrangian particulate material. The model comprises a rectangular container with an orifice on one side to simulate water flow near a downward-facing leak. The container is filled with saturated sand particles, each 4 mm in diameter. Preliminary results under laminar flow conditions indicate that for horizontally directed leaks, the streamline curves upward toward the upper surface at a certain distance from the origin. The flow entrains sand particles and displaces them until they reach the surface, forming a cyclic vertical circulation pattern around the leak origin. Equations from Peng & Fan [1], based on Ergun’s [2], effectively predict the onset of erosion in horizontal leaks. [1] Peng, Y., & Fan, L. (1997). Hydrodynamic characteristics of fluidization in liquid-solid tapered beds. Chemical Engineering Science, 52 (14), 2277–2290. [2] Ergun, S. (1952) Fluid Flow through Packed Columns. Chemical Engineering Progress, 48, 89-94.