Published online by Cambridge University Press: 23 November 2011
We present high-resolution numerical simulations of the Euler and Navier–Stokes equations for a pair of colliding dipoles. We study the possible approach to a finite-time singularity for the Euler equations, and contrast it with the formation of developed turbulence for the Navier–Stokes equations. We present numerical evidence that seems to suggest the existence of a blow-up of the inviscid velocity field at a finite time () with scaling
,
. This blow-up is associated with the formation of a
spectral range, at least for the finite range of wavenumbers that are resolved by our computation. In the evolution toward
, the total enstrophy is observed to increase at a slower rate,
, than would naively be expected given the behaviour of the maximum vorticity,
. This indicates that the blow-up would be concentrated in narrow regions of the flow field. We show that these regions have sheet-like structure. Viscous simulations, performed at various
, support the conclusion that any non-zero viscosity prevents blow-up in finite time and results in the formation of a dissipative exponential range in a time interval around the estimated inviscid
. In this case the total enstrophy saturates, and the energy spectrum becomes less steep, approaching
. The simulations show that the peak value of the enstrophy scales as
, which is in accord with Kolmogorov phenomenology. During the short time interval leading to the formation of an inertial range, the total energy dissipation rate shows a clear tendency to become independent of
, supporting the validity of Kolmogorov’s law of finite energy dissipation. At later times the kinetic energy shows a
decay for all
, in agreement with experimental results for grid turbulence. Visualization of the vortical structures associated with the stages of vorticity amplification and saturation show that, prior to
, large-scale and the small-scale vortical structures are well separated. This suggests that, during this stage, the energy transfer mechanism is non-local both in wavenumber and in physical space. On the other hand, as the spectrum becomes shallower and a
range appears, the energy-containing eddies and the small-scale vortices tend to be concentrated in the same regions, and structures with a wide range of sizes are observed, suggesting that the formation of an inertial range is accompanied by transfer of energy that is local in both physical and spectral space.