Particle-laden jet flows can be observed in many industrial applications. In this investigation, the horizontal motion of particle laden jets is simulated using the Eulerian–Lagrangian framework. The two-way coupling is applied to the model to simulate the interaction between discrete and continuum phase. In order to track the continuum phase, a passive scalar equation is added to the solver. Eddy Life Time (ELT) is employed as a dispersion model. The influences of different non-dimensional parameters, such as Stokes number, Jet Reynolds number and mass loading ratio on the flow characteristics, are studied. The results of the simulations are verified with the available experimental data. It is revealed that more gravitational force is exerted on the jet as a result of the increase in mass loading, which deflects it more. Moreover, with an increase in the Reynolds number, the speed of the jet rises, and consequently, the gravitational force becomes less capable of deviating the jet. In addition, it is observed that by increasing the Stokes number, the particles leave the jet at higher speed, which causes a lower deviation of the jet.