Intermittency in passive scalar advection
  A Lagrangian method for the numerical simulation of the Kraichnan passive scalar model is introduced. The method is based on Monte--Carlo simulations of tracer trajectories, supplemented by a point-splitting procedure for coinciding points. Clean scaling behavior for scalar structure functions is observed. The scheme is exploited to investigate the dependence of scalar anomalies on the scaling exponent $\xi$ of the advecting velocity field. The three-dimensional fourth-order structure function is specifically considered. 
U. Frisch, A. Mazzino, M. Vergassola