On the small-scale structure of turbulence and its impact on the pressure field

Understanding the small-scale structure of incompressible turbulence and its implications for the non-local pressure field is one of the fundamental challenges in fluid mechanics. Intense velocity gradient structures tend to cluster on a range of scales which affects the pressure through a Poisson equation. Here we present a quantitative investigation of the spatial distribution of these structures conditional on their intensity for Taylor-based Reynolds numbers in the range [160, 380]. We find that the correlation length, the second invariant of the velocity gradient, is proportional to the Kolmogorov scale. It also is a good indicator for the spatial localization of intense enstrophy and strain-dominated regions, as well as the separation between them. We describe and quantify the differences in the two-point statistics of these regions and the impact they have on the non-locality of the pressure field as a function of the intensity of the regions. Specifically, across the examined range of Reynolds numbers, the pressure in strong rotation-dominated regions is governed by a dissipation-scale neighbourhood. In strong strain-dominated regions, on the other hand, it is determined primarily by a larger neighbourhood reaching inertial scales.Comment: Accepted for publication by the Journal of Fluid Mechanic

Resolved energy budget of superstructures in Rayleigh-B\'{e}nard convection

Turbulent superstructures, i.e. large-scale flow structures in turbulent flows, play a crucial role in many geo- and astrophysical settings. In turbulent Rayleigh-B\'{e}nard convection, for example, horizontally extended coherent large-scale convection rolls emerge. Currently, a detailed understanding of the interplay of small-scale turbulent fluctuations and large-scale coherent structures is missing. Here, we investigate the resolved kinetic energy and temperature variance budgets by applying a filtering approach to direct numerical simulations of Rayleigh-B\'{e}nard convection at high aspect ratio. In particular, we focus on the energy transfer rate between large-scale flow structures and small-scale fluctuations. We show that the small scales primarily act as a dissipation for the superstructures. However, we find that the height-dependent energy transfer rate has a complex structure with distinct bulk and boundary layer features. Additionally, we observe that the heat transfer between scales mainly occurs close to the thermal boundary layer. Our results clarify the interplay of superstructures and turbulent fluctuations and may help to guide the development of an effective description of large-scale flow features in terms of reduced-order models

Nonlinear closures for scale separation in supersonic magnetohydrodynamic turbulence

Turbulence in compressible plasma plays a key role in many areas of astrophysics and engineering. The extreme plasma parameters in these environments, e.g. high Reynolds numbers, supersonic and super-Alfvenic flows, however, make direct numerical simulations computationally intractable even for the simplest treatment -- magnetohydrodynamics (MHD). To overcome this problem one can use subgrid-scale (SGS) closures -- models for the influence of unresolved, subgrid-scales on the resolved ones. In this work we propose and validate a set of constant coefficient closures for the resolved, compressible, ideal MHD equations. The subgrid-scale energies are modeled by Smagorinsky-like equilibrium closures. The turbulent stresses and the electromotive force (EMF) are described by expressions that are nonlinear in terms of large scale velocity and magnetic field gradients. To verify the closures we conduct a priori tests over 137 simulation snapshots from two different codes with varying ratios of thermal to magnetic pressure ($\beta_\mathrm{p} = 0.25, 1, 2.5, 5, 25$) and sonic Mach numbers ($M_s = 2, 2.5, 4$). Furthermore, we make a comparison to traditional, phenomenological eddy-viscosity and $\alpha-\beta-\gamma$ closures. We find only mediocre performance of the kinetic eddy-viscosity and $\alpha-\beta-\gamma$ closures, and that the magnetic eddy-viscosity closure is poorly correlated with the simulation data. Moreover, three of five coefficients of the traditional closures exhibit a significant spread in values. In contrast, our new closures demonstrate consistently high correlation and constant coefficient values over time and and over the wide range of parameters tested. Important aspects in compressible MHD turbulence such as the bi-directional energy cascade, turbulent magnetic pressure and proper alignment of the EMF are well described by our new closures.Comment: 15 pages, 6 figures; to be published in New Journal of Physic