Spray, a 3D Lagrangian particle model for the simulation of complex flow dispersion, is presented. Its performance is tested against the Environmental Protection Agency wind tunnel concentration distribution of passive tracer released from elevated point sources, located in the lee region of a two-dimensional schematic hill, in a neutrally stratified boundary layer. Based on the measured values of the first two moments of the turbulent flow velocity, the mean fields are computed over a regular grid using a mass-consistent model, whereas the turbulence structure is simply interpolated. From these fields, trajectories of tracer particles are computed using a linear formulation of the Langevin equation, with a correlated, skewed forcing. The self-consistence test (well-mixed condition), aimed at maintaining an initially well-mixed particle distribution uniform in time, has shown satisfactory results in the region under study. The computed concentration field turns out to be in good agreement with th... Abstract Spray, a 3D Lagrangian particle model for the simulation of complex flow dispersion, is presented. Its performance is tested against the Environmental Protection Agency wind tunnel concentration distribution of passive tracer released from elevated point sources, located in the lee region of a two-dimensional schematic hill, in a neutrally stratified boundary layer. Based on the measured values of the first two moments of the turbulent flow velocity, the mean fields are computed over a regular grid using a mass-consistent model, whereas the turbulence structure is simply interpolated. From these fields, trajectories of tracer particles are computed using a linear formulation of the Langevin equation, with a correlated, skewed forcing. The self-consistence test (well-mixed condition), aimed at maintaining an initially well-mixed particle distribution uniform in time, has shown satisfactory results in the region under study. The computed concentration field turns out to be in good agreement with th...