  We demonstrate the utility of the equation free methodology developed by one
of the authors (I.G.K) for the study of scalar conservation laws with
disordered initial conditions. The numerical scheme is benchmarked on exact
solutions in Burgers turbulence corresponding to Levy process initial data. For
these initial data, the kinetics of shock clustering is described by
Smoluchowski's coagulation equation with additive kernel. The equation free
methodology is used to develop a particle scheme that computes self-similar
solutions to the coagulation equation, including those with fat tails.
