Lagrangian method for multiple correlations in passive scaler advection

A Lagrangian method is introduced for calculating simultaneous n-point correlations of a passive scalar advected by a random velocity field, with random forcing and finite molecular diffusivity kappa, The method, which is here presented in detail, is particularly well suited for studying the kappa-->0 limit when the velocity field is not smooth. Efficient Monte Carlo simulations based on this method are applied to the Kraichnan model of passive scalar and lend to accurate determinations of the anomalous intermittency corrections in the fourth-order structure function as a function of the scaling exponent xi of the velocity field in two and three dimensions. Anomalous corrections are found to vanish in the limits xi-->0 and xi-->2, as predicted by perturbation theory, (C) 1999 American Institute of Physics, [S1070-6631(99)01808-5].