Numerical simulations of dispersion experiments within the planetary boundary layer are actually feasible making use of Large Eddy Simulations (LES). In Eulerian framework, a conservation equation for a passive scalar may be
superimposed on LES wind/turbulence fields to get a realistic description of timevarying concentration field. Aim of this work is to present a numerical technique to solve the Eulerian conservation equation. The technique is based on Fractional Step/Locally One-Dimensional (LOD) methods. Advection terms are calculated with a
semi-Lagrangian cubic-spline technique, while diffusive terms are calculated with Crank-Nicholson implicit scheme. To test the grid model, the dispersion of contaminants
emitted from an elevated continuous point source in a convective boundary layer is simulated. Results show that the calculated concentration distributions agree quite well with numerical and experimental data found in the literature
