Function crpropa::initTurbulence

Function Documentation

void crpropa::initTurbulence(ref_ptr < Grid3f > grid, double Brms, double lMin, double lMax, double alpha = -11/3., int seed = 0)

Create a random initialization of a turbulent field.

  • lMin: Minimum wavelength of the turbulence
  • lMax: Maximum wavelength of the turbulence
  • alpha: Power law index of <B^2(k)> ~ k^alpha (alpha = -11/3 corresponds to a Kolmogorov spectrum)
  • Brms: RMS field strength
  • seed: Random seed