Abstract : Particle-laden ﬂows occur in industrial applications ranging from droplets in gas turbines to ﬂuidized bed in chemical industry. Prediction of the dispersed phase properties such as concentration and dynamics are crucial for the design of more efﬁcient devices that meet the new pollutant regulations of the European community. The objective of this thesis is to develop an Euler-Lagrange formulation on a parallel and unstructured solver for large-eddy simulation. This work is motivated by the rapid increase in computing power which opens a new way for simulations that were prohibitive one decade ago. Special attention is taken to keep data structure simplicity and code portability. Developments are validated in two conﬁgurations: an academic test of a decaying homogeneous isotropic turbulence and a polydisperse two-phase ﬂow of a conﬁned bluff body. The use of load-balancing capabilities is highlighted as a promising solution in Lagrangian two-phase ﬂow simulations to improve performance when strong imbalance of the dispersed phase is present.