  We report on a comparison of high-resolution numerical simulations of
Lagrangian particles advected by incompressible turbulent hydro- and
magnetohydrodynamic (MHD) flows. Numerical simulations were performed with up
to $1024^3$ collocation points and 10 million particles in the Navier-Stokes
case and $512^3$ collocation points and 1 million particles in the MHD case. In
the hydrodynamics case our findings compare with recent experiments from
Mordant et al. [1] and Xu et al. [2]. They differ from the simulations of
Biferale et al. [3] due to differences of the ranges choosen for evaluating the
structure functions. In Navier-Stokes turbulence intermittency is stronger than
predicted by a multifractal approach of [3] whereas in MHD turbulence the
predictions from the multifractal approach are more intermittent than observed
in our simulations. In addition, our simulations reveal that Lagrangian
Navier-Stokes turbulence is more intermittent than MHD turbulence, whereas the
situation is reversed in the Eulerian case. Those findings can not consistently
be described by the multifractal modeling. The crucial point is that the
geometry of the dissipative structures have different implications for
Lagrangian and Eulerian intermittency. Application of the multifractal approach
for the modeling of the acceleration PDFs works well for the Navier-Stokes case
but in the MHD case just the tails are well described.
