Small-scale dynamos are expected to operate in all astrophysical fluids that are turbulent and electrically conducting, for example the interstellar medium, stellar interiors, and accretion disks, where they may also be affected by or competing with large-scale dynamos. However, the possibility of small-scale dynamos being excited at small and intermediate ratios of viscosity to magnetic diffusivity (the magnetic Prandtl number) has been debated, and the possibility of them depending on the large-scale forcing wavenumber has been raised. Here we show, using four values of the forcing wavenumber, that the small-scale dynamo does not depend on the scale-separation between the size of the simulation domain and the integral scale of the turbulence, i.e., the forcing scale. Moreover, the spectral bottleneck in turbulence, which has been implied as being responsible for raising the excitation conditions of small-scale dynamos, is found to be invariant under changing the forcing wavenumber. However, when forcing at the lowest few wavenumbers, the effective forcing wavenumber that enters in the definition of the magnetic Reynolds number is found to be about twice the minimum wavenumber of the domain. Our work is relevant to future studies of small-scale dynamos, of which several applications are being discussed.

It is known that unconfined dust explosions consist of a relatively weak primary (turbulent) deflagrations followed by a devastating secondary explosion. The secondary explosion may propagate with a speed of up to 1000 m/s producing overpressures of over 8-10 atm. Since detonation is the only established theory that allows a rapid burning producing a high pressure that can be sustained in open areas, the generally accepted view was that the mechanism explaining the high rate of combustion in dust explosions is deflagration to detonation transition. In the present work we propose a theoretical substantiation of the alternative propagation mechanism explaining origin of the secondary explosion producing the high speeds of combustion and high overpressures in unconfined dust explosions. We show that clustering of dust particles in a turbulent flow gives rise to a significant increase of the thermal radiation absorption length ahead of the advancing flame front. This effect ensures that clusters of dust particles are exposed to and heated by the radiation from hot combustion products of large gaseous explosions sufficiently long time to become multi-point ignition kernels in a large volume ahead of the advancing flame front. The ignition times of fuel-air mixture by the radiatively heated clusters of particles is considerably reduced compared to the ignition time by the isolated particle. The radiation-induced multi-point ignitions of a large volume of fuel-air ahead of the primary flame efficiently increase the total flame area, giving rise to the secondary explosion, which results in high rates of combustion and overpressures required to account for the observed level of overpressures and damages in unconfined dust explosions, such as e.g. the 2005 Buncefield explosion and several vapor cloud explosions of severity similar to that of the Buncefield incident.

It is known that unconfined dust explosions typically starts off with a relatively weak primary flame followed by a severe secondary explosion. We show that clustering of dust particles in a temperature stratified turbulent flow ahead of the primary flame may give rise to a significant increase in the radiation penetration length. These particle clusters, even far ahead of the flame, are sufficiently exposed and heated by the radiation from the flame to become ignition kernels capable to ignite a large volume of fuel-air mixtures. This efficiently increases the total flame surface area and the effective combustion speed, defined as the rate of reactant consumption of a given volume. We show that this mechanism explains the high rate of combustion and overpressures required to account for the observed level of damage in unconfined dust explosions, e.g., at the 2005 Buncefield vapor-cloud explosion. The effect of the strong increase of radiation transparency due to turbulent clustering of particles goes beyond the state-of-the-art of the application to dust explosions and has many implications in atmospheric physics and astrophysics.

The phenomenon of turbulent thermal diffusion in temperature-stratified turbulence causing a non-diffusive turbulent flux of inertial and non-inertial particles in the direction of the turbulent heat flux is found using direct numerical simulations (DNS). In simulations with and without gravity, this phenomenon is found to cause a peak in the particle number density around the minimum of the mean fluid temperature for Stokes numbers less than 1, where the Stokes number is the ratio of particle Stokes time to turbulent Kolmogorov time at the viscous scale. Turbulent thermal diffusion causes the formation of large-scale inhomogeneities in the spatial distribution of inertial particles. The strength of this effect is maximum for Stokes numbers around unity, and decreases again for larger values. The dynamics of inertial particles is studied using Lagrangian modelling in forced temperature-stratified turbulence, whereas non-inertial particles and the fluid are described using DNS in an Eulerian framework.

Direct numerical simulations and mean-field theory are used to model reactive front propagation in a turbulent medium. In the mean-field approach, memory effects of turbulent diffusion are taken into account to estimate the front speed in cases when the Damkohler number is large. This effect is found to saturate the front speed to values comparable with the speed of the turbulent motions. By comparing with direct numerical simulations, it is found that the effective correlation time is much shorter than for non-reacting flows. The nonlinearity of the reaction term is found to make the front speed slightly faster.

A high-order scheme for direct numerical simulations of turbulent combustion is discussed. Its implementation in the massively parallel and publicly available Pencil Code is validated with the focus on hydrogen combustion. Ignition delay times (0D) and laminar flame velocities (1D) are calculated and compared with results from the commercially available Chemkin code. The scheme is verified to be fifth order in space. Upon doubling the resolution, a 32-fold increase in the accuracy of the flame front is demonstrated. Finally, also turbulent and spherical flame front velocities are calculated and the implementation of the non-reflecting so-called Navier-Stokes Characteristic Boundary Condition is validated in all three directions.

Direct and large eddy simulations of hydrodynamic and hydromagnetic turbulence have been performed in an attempt to isolate artifacts from real and possibly asymptotic features in the energy spectra. It is shown that in a hydrodynamic turbulence simulation with a Smagorinsky subgrid scale model using 512^3 meshpoints two important features of the 4096^3 simulation on the Earth simulator (Kaneda et al. 2003, Phys. Fluids 15, L21) are reproduced: a k^{-0.1} correction to the inertial range with a k^{-5/3} Kolmogorov slope and the form of the bottleneck just before the dissipative subrange. Furthermore, it is shown that, while a Smagorinsky-type model for the induction equation causes an artificial and unacceptable reduction in the dynamo efficiency, hyper-resistivity yields good agreement with direct simulations. In the large-scale part of the inertial range, an excess of the spectral magnetic energy over the spectral kinetic energy is confirmed. However, a trend towards spectral equipartition at smaller scales in the inertial range can be identified. With magnetic fields, no explicit bottleneck effect is seen.

Random motions can occur in the intergalactic gas of galaxy clusters at all stages of their evolution. Depending on the poorly known value of the Reynolds number, these motions can or cannot become turbulent, but in any case they can generate random magnetic fields via dynamo action. We argue that magnetic fields inferred observationally for the intracluster medium require dynamo action, and then estimate parameters of random flows and magnetic fields at various stages of the cluster evolution. Polarization in cluster radio halos predicted by the model would be detectable with the SKA.

We discuss, using simple analytical models and MHD simulations, the origin and parameters of turbulence and magnetic fields in galaxy clusters. Three physically distinct regimes can be identified in the evolution of cluster turbulence and magnetic fields. Firstly, the fluctuation dynamo will produce microgauss-strong, random magnetic fields during cluster formation and major mergers. Turbulent velocity of about 300 km/s can be maintained at scales 100-200 kpc. The magnetic field is intermittent, has a smaller scale of 20-30 kpc and average strength of 2 microgauss. Secondly, when major mergers end, turbulent speed and magnetic field undergo a power-law decay, decreasing in strength but increasing in scale by a factor of about two. Thirdly, smaller-mass subclusters and cluster galaxies produce turbulent wakes, with turbulent speeds and magnetic field strengths similar to those quoted above. The velocity scales are about 200 kpc and 10 kpc respectively, and the magnetic field scale is about 6 times smaller. Although these wakes may fill only a small fraction of the cluster volume, their area covering factor can be close to unity. So one can potentially reconcile observations that indicate the coexistence of turbulence with ordered filamentary gas structures, as in the Perseus cluster. Random Faraday rotation measure is estimated to be typically 100-200 rad/m^2, in agreement with observations. We predict detectable synchrotron polarization from cluster radio halos at wavelengths 3-6 cm, if observed at sufficiently high resolution (abridged).

We study numerically the dependence of the critical magnetic Reynolds number Rmc for the turbulent small-scale dynamo on the hydrodynamic Reynolds number Re. The turbulence is statistically homogeneous, isotropic, and mirror--symmetric. We are interested in the regime of low magnetic Prandtl number Pm=Rm/Re<1, which is relevant for stellar convective zones, protostellar disks, and laboratory liquid-metal experiments. The two asymptotic possibilities are Rmc->const as Re->infinity (a small-scale dynamo exists at low Pm) or Rmc/Re=Pmc->const as Re->infinity (no small-scale dynamo exists at low Pm). Results obtained in two independent sets of simulations of MHD turbulence using grid and spectral codes are brought together and found to be in quantitative agreement. We find that at currently accessible resolutions, Rmc grows with Re with no sign of approaching a constant limit. We reach the maximum values of Rmc~500 for Re~3000. By comparing simulations with Laplacian viscosity, fourth-, sixth-, and eighth-order hyperviscosity and Smagorinsky large-eddy viscosity, we find that Rmc is not sensitive to the particular form of the viscous cutoff. This work represents a significant extension of the studies previously published by Schekochihin et al. 2004, PRL 92, 054502 and Haugen et al. 2004, PRE, 70, 016308 and the first detailed scan of the numerically accessible part of the stability curve Rmc(Re).

Small scale turbulent dynamo action in compressible transonic turbulence is discussed. It is shown that the critical value of the magnetic Reynolds number displays a bimodal behavior and changes from a typical value of 35 for small Mach numbers to about 80 for larger Mach numbers. The transition between the two regimes is relatively sharp. The direct simulations are then compared with simulations where shocks are captured using a shock viscosity that becomes large at locations where there are shocks. In the presence of shear it is shown that large scale dynamo action is possible.

Two different computational approaches to the magnetorotational instability (MRI) are pursued: the shearing box approach which is suited for local simulations and the embedding box approach whereby a Taylor Couette flow is embedded in a box so that numerical problems with the coordinate singularity are avoided. New shearing box simulations are presented and differences between regular and hyperviscosity are discussed. Preliminary simulations of spherical nonlinear Taylor Couette flow in an embedding box are presented and the effects of an axial field on the background flow are studied.

Three closely related stumbling blocks of solar mean field dynamo theory are discussed: how dominant are the small scale fields, how is the alpha effect quenched, and whether magnetic and current helicity fluxes alleviate the quenching? It is shown that even at the largest currently available resolution there is no clear evidence of power law scaling of the magnetic and kinetic energy spectra in turbulence. However, using subgrid scale modeling, some indications of asymptotic equipartition can be found. The frequently used first order smoothing approach to calculate the alpha effect and other transport coefficients is contrasted with the superior minimal tau approximation. The possibility of catastrophic alpha quenching is discussed as a result of magnetic helicity conservation. Magnetic and current helicity fluxes are shown to alleviate catastrophic quenching in the presence of shear. Evidence for strong large scale dynamo action, even in the absence of helicity in the forcing, is presented.

Numerical turbulence with hyperviscosity is studied and compared with direct simulations using ordinary viscosity and data from wind tunnel experiments. It is shown that the inertial range scaling is similar in all three cases. Furthermore, the bottleneck effect is approximately equally broad (about one order of magnitude) in these cases and only its height is increased in the hyperviscous case--presumably as a consequence of the steeper decent of the spectrum in the hyperviscous subrange. The mean normalized dissipation rate is found to be in agreement with both wind tunnel experiments and direct simulations. The structure function exponents agree with the She-Leveque model. Decaying turbulence with hyperviscosity still gives the usual t^{-1.25} decay law for the kinetic energy, and also the bottleneck effect is still present and about equally strong.

Non-helical hydromagnetic turbulence with an externally imposed magnetic field is investigated using direct numerical simulations. It is shown that the imposed magnetic field lowers the spectral magnetic energy in the inertial range. This is explained by a suppression of the small scale dynamo. At large scales, however, the spectral magnetic energy increases with increasing imposed field strength for moderately strong fields, and decreases only slightly for even stronger fields. The presence of Alfven waves is explicitly confirmed by monitoring the evolution of magnetic field and velocity at one point. The frequency omega agrees with vA k1, where vA is the Alfven speed and k1 is the smallest wavenumber in the box.

The effect of compressibility on the onset of nonhelical turbulent dynamo action is investigated using both direct simulations as well as simulations with shock-capturing viscosities, keeping however the regular magnetic diffusivity. It is found that the critical magnetic Reynolds number increases from about 35 in the subsonic regime to about 70 in the supersonic regime. In the high resolution direct simulations the shock structures are much sharper than in the low resolution shock-capturing simulations. Nevertheless, the magnetic field looks roughly similar in both cases and does not show shock structures. Similarly, the onset of dynamo action is not significantly affected by the shock-capturing viscosity.

Nonhelical hydromagnetic forced turbulence is investigated using large scale simulations on up to 256 processors and $1024^3$ meshpoints. The magnetic Prandtl number is varied between 1/8 and 30, although in most cases it is unity. When the magnetic Reynolds number is based on the inverse forcing wavenumber, the critical value for dynamo action is shown to be around 35 for magnetic Prandtl number of unity. For small magnetic Prandtl numbers we find the critical magnetic Reynolds number to increase with decreasing magnetic Prandtl number. The Kazantsev $k^{3/2}$ spectrum for magnetic energy is confirmed for the kinematic regime, i.e. when nonlinear effects are still unimportant and when the magnetic Prandtl number is unity. In the nonlinear regime, the energy budget converges for large Reynolds numbers (around 1000) such that for our parameters about 70% is in kinetic energy and about 30% is in magnetic energy. The energy dissipation rates are converged to 30% viscous dissipation and 70% resistive dissipation. Second order structure functions of the Elsasser variables give evidence for a $k^{-5/3}$ spectrum. Nevertheless, the three-dimensional spectrum is close to $k^{-3/2}$, but we argue that this is due to the bottleneck effect. The bottleneck effect is shown to be equally strong both for magnetic and nonmagnetic turbulence, but it is far weaker in one-dimensional spectra that are normally studied in laboratory turbulence. Structure function exponents for other orders are well described by the She-Leveque formula, but the velocity field is significantly less intermittent and the magnetic field is more intermittent than the Elsasser variables.

Decaying turbulence is studied numerically using as initial condition a random flow whose shell-integrated energy spectrum increases with wavenumber k like k^q. Alternatively, initial conditions are generated from a driven turbulence simulation by simply stopping the driving. It is known that the dependence of the decaying energy spectrum on wavenumber, time, and viscosity can be collapsed onto a unique scaling function that depends only on two parameters. This is confirmed using three-dimensional simulations and the dependence of the scaling function on its two arguments is determined.

The dimensionless kinetic energy dissipation rate C_epsilon is estimated from numerical simulations of statistically stationary isotropic box turbulence that is slightly compressible. The Taylor microscale Reynolds number Re_lambda range is 20 < Re_lambda < 220 and the statistical stationarity is achieved with a random phase forcing method. The strong Re_lambda dependence of C_epsilon abates when Re_lambda approx. 100 after which C_epsilon slowly approaches approx 0.5 a value slightly different to previously reported simulations but in good agreement with experimental results. If C_epsilon is estimated at a specific time step from the time series of the quantities involved it is necessary to account for the time lag between energy injection and energy dissipation. Also, the resulting value can differ from the ensemble averaged value by up to +-30%. This may explain the spread in results from previously published estimates of C_epsilon.

Nonhelical hydromagnetic turbulence without an imposed magnetic field is considered in the case where the magnetic Prandtl number is unity. The magnetic field is entirely due to dynamo action. The magnetic energy spectrum peaks at a wavenumber of about 5 times the minimum wavenumber in the domain, and not at the resistive scale, as has previously been argued. Throughout the inertial range the spectral magnetic energy exceeds the kinetic energy by a factor of about 2.5, and both spectra are approximately parallel. At first glance, the total energy spectrum seems to be close to k^{-3/2}, but there is a strong bottleneck effect and it is suggested that the asymptotic spectrum is k^{-5/3}. This is supported by the value of the second order structure function exponent that is found to be \zeta_2=0.70, suggesting a k^{-1.70} spectrum.

At numerical resolutions around $512^3$ and above, three-dimensional energy spectra from turbulence simulations begin to show noticeably shallower spectra than $k^{-5/3}$ near the dissipation wavenumber (`bottleneck effect'). This effect is shown to be significantly weaker in one-dimensional spectra such as those obtained in wind tunnel turbulence. The difference can be understood in terms of the transformation between one-dimensional and three-dimensional energy spectra under the assumption that the turbulent velocity field is isotropic. Transversal and longitudinal energy spectra are similar and can both accurately be computed from the full three-dimensional spectra. Second-order structure functions are less susceptible to the bottleneck effect and may be better suited for inferring the scaling exponent from numerical simulation data.

According to the kinematic theory of nonhelical dynamo action the magnetic energy spectrum increases with wavenumber and peaks at the resistive cutoff wavenumber. It has previously been argued that even in the dynamical case the magnetic energy peaks at the resistive scale. Using high resolution simulations (up to 1024^3 meshpoints) with no large scale imposed field we show that the magnetic energy peaks at a wavenumber that is independent of the magnetic Reynolds number and about 5 times larger than the forcing wavenumber. Throughout the inertial range the spectral magnetic energy exceeds the kinetic energy by a factor of 2 to 3. Both spectra are approximately parallel. The total energy spectrum seems to be close to k^{-3/2}, but there is a strong bottleneck effect and it is suggested that the asymptotic spectrum is instead k^{-5/3}. This is supported by the value of the second order structure function exponent that is found to be zeta_2=0.70, suggesting a k^{-1.70} spectrum. The third order structure function scaling exponent is very close to unity, in agreement with Goldreich-Sridhar theory. Adding an imposed field tends to suppress the small scale magnetic field. We find that at large scales the magnetic power spectrum follows then a k^{-1} slope. When the strength of the imposed field is of the same order as the dynamo generated field, we find almost equipartition between the magnetic and kinetic power spectra.

Isotropic homogeneous hydromagnetic turbulence is studied using numerical simulations at resolutions of up to 1024^3 meshpoints. It is argued that, in contrast to the kinematic regime, the nonlinear regime is characterized by a spectral magnetic power that is decreasing with increasing wavenumber, regardless of whether or not the turbulence has helicity. This means that the root-mean-square field strength converges to a limiting value at large magnetic Reynolds numbers. The total (magnetic and kinetic) energy spectrum tends to be somewhat shallower than k^{-5/3}, in agreement with the findings of other groups. In the presence of helicity, an inverse cascade develops, provided the scale separation between the size of the computational box and the scale of the energy carrying eddies exceeds a ratio of at least two. Finally, the constraints imposed by magnetic helicity conservation on mean-field theory are reviewed and new results of simulations are presented.