Robustness of vortex populations in the two-dimensional inverse energy cascade

We study how the properties of forcing and dissipation affect the scaling behaviour of the vortex population in the two-dimensional turbulent inverse energy cascade. When the flow is forced at scales intermediate between the domain and dissipation scales, the growth rates of the largest vortex area and the spectral peak length scale are robust to all simulation parameters. For white-in-time forcing the number density distribution of vortex areas follows the scaling theory predictions of Burgess & Scott (J. Fluid Mech., vol. 811, 2017, pp. 742–756) and shows little sensitivity either to the forcing bandwidth or to the nature of the small-scale dissipation: both narrowband and broadband forcing generate nearly identical vortex populations, as do Laplacian diffusion and hyperdiffusion. The greatest differences arise in comparing simulations with correlated forcing to those with white-in-time forcing: in flows with correlated forcing the intermediate range in the vortex number density steepens significantly past the predicted scale-invariant A^-1 scaling. We also study the impact of the forcing Reynolds number Re_j, a measure of the relative importance of nonlinear terms and dissipation at the forcing scale, on vortex formation and the scaling of the number density. As Re_j  decreases, the flow changes from one dominated by intense circular vortices surrounded by filaments to a less structured flow in which vortex formation becomes progressively more suppressed and the filamentary nature of the surrounding vorticity field is lost. However, even at very small Re_j, and in the absence of intense coherent vortex formation, regions of anomalously high vorticity merge and grow in area as predicted by the scaling theory, generating a three-part number density similar to that found at higher Re_j. At late enough stages the aggregation process results in the formation of long-lived circular vortices, demonstrating a strong tendency to vortex formation, and via a route distinct from the axisymmetrization of forcing extrema seen at higher Re_j. Our results establish coherent vortices as a robust feature of the two-dimensional inverse energy cascade, and provide clues as to the dynamical mechanisms shaping their statistics.