In the present study, a compact and fast MATLAB program coupled with the Taylor-series expansion method of moments (TEMOM) was developed to simulate the effect of coherent structures on the particle Brownian coagulation in the temporal mixing layer. The distributions of number concentration, mass concentration and particle average volume driven by coagulation, advection and diffusion are obtained. The developed TEMOM method has no prior requirement for the particle size distribution (PSD), and it is a promising method to approximate the aerosol general dynamics equation (GDE). The fluid and particle fields are coupled together and are presented with three non-dimensional parameters (i.e.. Reynolds number, Re; Schmidt number based on particle moment, SCM, and Damkohler number, Da) in the governing equations. The temporal evolutions of the first three moments are discussed for different Damkohler numbers (Da=0.5, 1.0, 2.0). As the fluid flow evolves in time, the number concentration of nanoparticles decreases gradually, while the particle average volume increase. The distribution of number concentration, mass concentration, and average volume of nanoparticles are spatially inhomogeneous due to the mixing of coherent vortex structures. Far away from the eddy structure, the effect of the fluid advection on particle coagulation is small: however, the particle coagulation within the eddy core has an obvious wave-like distribution because of the large-scale eddy. The results reveal that the coherent structures play a significant role in the particle Brownian coagulation in the mixing layer. The particle coagulation affects quantitatively the distribution of particle number concentration, volume concentration and average diameter, but the qualitative characteristics of these distributions remain unchanged. (C) 2012 Elsevier Ltd. All rights reserved.