Simulations of Idealised Turbulence

To investigate the correctness of the RMS method, we use idealised simulations of subsonic turbulence. This allows us to compare the turbulent energy, estimated by the RMS method, with an estimate of the same energy in Fourier space.

Initial Conditions

We set up initial conditions in a periodic box, mimicking physical conditions of a galaxy cluster (see table 7.1). Positions are set using a glass file with 100³ particles to achieve a homogeneous Poisson sampling of the volume. To initialise particle velocities we set up a random vector field on a grid with 100³ grid-points, assuming a Kolmogorov-like power spectrum with a spectral index of −11/3. This vector field is then mapped to the particles using the nearest grid-point method (Hockney & Eastwood, 1988). The simulation is evolved for a few Gyr.

The particle velocities of the outputs are then mapped in a mass weighted way to a grid of equal resolution, using the cloud-in-cell technique.

7.2 Local Turbulence in SPH 135

-1.0 -0.5 0.0 0.5 1.0

x [h_sml] 0.0

0.5 1.0 1.5 2.0 2.5 3.0

3D SPH Kernel [1/h3 ]

Figure 7.2: Standard spline SPH kernel, as implemented in GADGET. Scales with strong smoothing are marked gray.

Errors

The procedure described above introduces considerable errors to the final velocity grid, which modify the power spectrum from the initial grid on smallest scales. As we are interested in these scales, we have to carefully treat most of them:

1. The mapping procedure from the initial grid to the particles involves a kernel es-timate, which modifies the power spectrum near the grid scale. We found that we can minimise this effect by using the Nearest-Grid-Point (NGP) method. It is of no importance to the energy estimate, because this spectrum is measured by the RMS method as well as the Fourier method. However, it explains why we are not recovering exactly the seeded Kolmogorov law everywhere.

2. The positions have to be homogeneous, isotropic and Poisson-like to properly sample the smallest scales. A grid-like particle distribution is unsuitable for these purposes.

We use a glass to achieve best results.

3. The finite amount of particles introduces the shot-noise error below the SPH kernel scale. This effect can be accounted for (Jing, 2005) in the case of homogeneous sampling.

4. The particle mapping to the final grid involves another kernel estimate that acts on the power-spectrum as well as the shot-noise. This can be accounted for by

136 7. Reacceleration of Cosmic Ray Electrons in Cluster Mergers

deconvolving with the kernel Fourier transform in k-space, which is given for the cloud-in-cell (CIC) kernel in e.g. Hockney & Eastwood (1988).

Power Spectra

In figure 7.3 and 7.4 we show velocity power spectra of the sub/supersonic simulations at different times. In each panel we plot the seeded power spectrum (black), the uncorrected power spectrum from the grid (blue), the deconvolved power spectrum (red, 7.2.1, point 4), and the deconvolved, shot noise correct spectrum (green). As red, red dashed lines we mark the mean,min and max hsml scale in the simulation, respectively. In addition we over-plot in the first panel (t = 0) the CIC kernel (black, dashed) and the shot noise (black, dashed dotted).

After correcting for CIC kernel and shotnoise (green) the initial power spectrum follows the seeded power spectrum (black) very well. Slight differences arise below the kernel scale.

We attribute this to the NGP mapping of the seeded power spectrum to the particles. At this point due to the homogeneous sampling the effect of kernel and shot noise correction on the power spectrum are rather small. Due to the higher energy content the effect of the shot-noise on the power spectrum is less prominent in the supersonic case.

As the simulations evolve, the power spectrum becomes increasingly shot noise dom-inated on the smallest scales, until after one sound-crossing time, when the whole kernel scale is shot noise dominated. Again the effect is more prominent in the subsonic case.

Above the hsml scale the spectrum steepens with time. Additionally, at intermediate times the spectrum shows wiggles near the maximum kernel scale, i.e. in the most under-dense regions. The supersonic simulation shows a larger density contrast than the subsonic case, which leads to a less wiggly spectrum at intermediate times.

In both cases the increasing importance of the shot noise suggests that on smallest scales correction terms become dominant with time. Additionally the density contrast in the box and inhomogeneity of the sampling increases, as shown by the spread in hsml.

That makes the sampling noise spatially dependent and our shotnoise correction more and more incorrect. Therefore we conclude that any estimation of the power spectrum or the turbulent energy on smallest scales via a grid/FFT based technique becomes more uncertain with higher density gradients. This has to be kept in mind, when comparing a particle based methods with FFT/grid based methods. So far there is no way to correct for spatially inhomogeneous shot noise in Fourier space.

Energy Evolution

In figure 7.5 we show the evolution of the turbulent energy on the mean kernel scale inside the box, the RMS method (green) and the FFT method (purple). We correct the RMS velocity dispersion for the spread in hsml. Using Parsevals theorem and assuming a Kolmogorov power-law, v_rms scales like:

v_rms,0² =v_rms,hsml² k₀

k_hsml −2/3

, (7.2)

7.2 Local Turbulence in SPH 137

Figure 7.3: Power spectra of simulations of ideal subsonic turbulence E_turb = 0.1E_thermal. We plot the seeded power spectrum (black), the uncorrected spectrum (blue), the decon-volved spectrum (red) and the decondecon-volved and shot noise corrected spectrum (green).

We also over-plot the CIC kernel (black dashed) and Poisson shot noise term (black, dot-ted dashed) in the first panel. The red vertical lines represent the kernel scale and its maximum/minimum.

138 7. Reacceleration of Cosmic Ray Electrons in Cluster Mergers

Figure 7.4: Power spectra of simulations of ideal supersonic turbulence (E_turb = 0.5E_thermal). We plot the seeded power spectrum (black), the uncorrected spectrum (blue), the deconvolved spectrum (red) and the deconvolved and shot noise corrected spectrum (green). We also over-plot the CIC kernel (black dashed) and Poisson shot noise term (black, dotted dashed) in the first panel. The red vertical lines represent the kernel scale and its maximum/minimum.

7.2 Local Turbulence in SPH 139

Figure 7.5: Time evolution of the turbulent energy on mean kernel scale in a simulation of ideal subsonic turbulence. In green we plot the measured turbulent energy using the particle based RMS approach. In purple we plot the same, but measured using the standard Fourier transform approach.

where the index 0 is to the reference scale and the index ”hsml” refers to the particle scale.

While this correction can be significant for the velocity dispersion of a single particle, it is only a few percent of the total energy in the simulated volume.

To estimate the turbulent energy in Fourier space, we bin the mass weighted particle velocity to a grid of the same size as the particle number. Similar to the previous section, we Fourier transform the grid and correct for the CIC kernel and the shot noise of a homogeneous particle distribution in k-space. The resulting power grid is then filtered, using an ideal high pass filter, and the total energy in all modes k_hsml is calculated.

Over about ten Gyr, the turbulent energies on the kernel scale agree reasonably well, as measured by both techniques. This is remarkable, considering the increasing uncertainties regarding corrections to the power on the grid in Fourier space. We found the relative error between both methods to be between 4% at early times and 20 % at late times. Given the discussion above, we attribute these errors to the uncertainties in the velocity grid and the inhomogeneous sampling. We conclude that the local turbulent energy in an SPH simulation can be measured reasonably well by the velocity dispersion around the kernel mean.

140 7. Reacceleration of Cosmic Ray Electrons in Cluster Mergers