This paper presents a methodology to combine stochastic Lagrangian approach and continuum model to simulate the dispersed phase in gas-particle turbulent flows using that both approaches are based on the same Boltzmann-like kinetic equation governing the joint fluid-particle probability density function (pdf). The proposed hybrid method is based on the separate application of each approach in two adjacent domains and their coupling at the interface via flux boundary conditions. Validation of the method is carried out for non-colliding solid particles suspended in homogeneous turbulent shear flow without two-way coupling. (c) 2007 Elsevier Ltd. All rights reserved.