  We present an analytical scheme, easily implemented numerically, to generate
synthetic Gaussian turbulent flows by using a linear Langevin equation, where
the noise term acts as a stochastic stirring force. The characteristic
parameters of the velocity field are well introduced, in particular the
kinematic viscosity and the spectrum of energy. As an application, the
diffusion of a passive scalar is studied for two different energy spectra.
Numerical results are compared favorably with analytical calculations.
