Reactive Transport of Suspended Particles

Monte Carlo simulation of the total probability density function for each constituent easily handles non-linear reactions and has a cost that increases only linearly with the number of reactants. This methodology had been extended to handle suspended particles in which different size classes considered as separate constituents are allowed to settle and react differentially. The coupling between turbulent shear and coagulation is simulated by treating local turbulent dissipation as a randomly distributed quantity. Mixing on a molecular scale is computed using a mixing and coalescing model as improved by S. B. Pope.