This paper focuses on the solution of the dispersed phase of Eulerian–Lagrangian one-way coupled particle laden flows. An efficient two-constraint domain partitioning for 2D and 3D unstructured hybrid meshes is proposed and implemented in distributed memory architectures. A preliminary simulation, using a set of representative particles, is performed first to suitably tag the cells with a weight proportional to the probability of being crossed by a particle. In addition, an innovative parallel ray-tracing location algorithm is presented. A global identifier is assigned to each particle resulting in a significant reduction of the overall communication among processes.The proposed approaches are verified against two steady reference cases for ice accretion simulation: a NACA 0012 profile and a NACA 64A008 swept horizontal tail. Furthermore, a cloud droplet impact test case starting from an unsteady flow around a 3D cylinder is performed to evaluate the code performances on unsteady problems.