We study the flow topology dynamics in terms of the paramount non-linearities of enstrophy and strain production at hard turbulent regimes of Rayleigh-Bénard convection (RBC). To do so, a data-set of direct numerical simulations for an air turbulent RBC at Rayleigh numbers Ra = {10 8 , 10 10 , 10 11 } is analysed. Considering the bulk dynamics therein, the classical 2D mean Lagrangian evolution of Q G and R G invariants of G ≡ ∇u is extended to 3D by decomposing R G into two parts: the strain production, R S , and the enstrophy production, tr(Ω 2 S). In this way, the 3D phase space (Q G , R S , tr(Ω 2 S)) allows to identify separately the non-linear straining and rotational mechanisms in turbulence. The main resultant observations attest that, when the turbulent regime is notably hard, a rising local self-amplification of the velocity gradient takes place in strain-dominated areas. This process is strongly aided by vortex contraction. In concomitant, a pronounced increase in the linear contributions of vortex stretching is also identified, particularly relevant to strain-dominated slots.
I. INTRODUCTIONMany circulations in nature and industry, such as convection in the outer layer of the Sun, coherent structures in the Earth's atmosphere and oceans, mantle convection in the Earth's core, circulations in nuclear reactors and solar thermal power plants, are ruled by Rayleigh-Bénard convection (RBC). Namely, the turbulent dynamics therein is mainly stemmed from buoyancy variations in the dynamo of a thermally driven flow heated from below and cooled from above [1][2][3]. Besides the onset of flow structures, this dynamics becomes of significant complexity when the grade of turbulence and thermal forcing is very high, i.e. Rayleigh number Ra > 10 10 . For instance, the rising self-sustained instabilities of RBC induce augmenting counter-gradient diffusion and energetic nonequilibriums between the buoyant production and viscous dissipation, which are mainly compensated by pressure fluctuations [4,5]. Although important features have been explored using direct numerical simulation (DNS) of RBC at hard turbulent regimes [6,7], such as the stable boundary layers at Ra = 2 × 10 12 [8] and the thermal plumes statistics at Ra = 10 12 [9]; many