Direct numerical simulations in a low-curvature viscoelastic turbulent Taylor vortex flow, with Reynolds numbers ranging from 1500 to 8000 and maximum chain extensibility ($L$) from 50 to 200, reveal a maximum drag reduction (MDR) asymptote. Compared with the classical MDR observed in planar wall-bounded shear flows, that is, drag reduction (DR) is $\sim -80\, \%$, this MDR state achieves only moderate levels of DR ($\sim -60\,\%$). This is due to the existence of large-scale structures (LSSs). A careful examination of the flow structures reveals that the polymer–turbulence interaction suppresses small-scale vortices and stabilizes the LSSs. These structural changes in turn lead to a reduction of Reynolds stress, and consequently to a DR flow state. Although Reynolds stress does not vanish as observed in classical MDR states, the small-scale vortices that heavily populate the near-wall region are also almost completely eliminated in this flow state. Concurrently, significant polymer stresses develop as a consequence of the interaction between polymer chains and LSSs that partially offset the magnitude of DR, leading to MDR asymptotes with moderate levels of DR. Moreover, we demonstrate that polymer deformation, i.e. deviation from the equilibrium state, is directly correlated with the LSSs dynamics, while the polymer deformation fluctuation displays a universal property in the MDR state. Hence, it is not surprising that the extent of DR exhibits a non-monotonic dependence on the maximum chain extensibility. Specifically, the variation in $L$ alters the incoherent and coherent angular momentum transport by small- and large-scale flow structures, respectively. To that end, the most DR flow state occurs at a moderate value $L=100$. Overall, this study further supports the universal property of polymer-induced asymptotic states in wall-bounded turbulence and paves the way for mechanistic understanding of drag modification that arises from the interaction of polymers with small- and large-scale flow structures.

Turbulent flow induced by elastorotational instability in viscoelastic Taylor–Couette flow (TCF) with Keplerian rotation is analogous to a turbulent accretion disk destabilized by magnetorotational instability. We examine this novel viscoelastic Keplerian turbulence via direct numerical simulations (DNS) for the shear Reynolds number ($Re$) ranging from $10^2$ to $10^4$. The observed characteristic flow structure consists of penetrating streamwise vortices with axial length scales much smaller than the gap width, distinct from the classic centrifugally induced Taylor vortices, which have axial lengths of the gap width. These intriguing vortices persist for the wide $Re$ range considered and give rise to intriguing scaling behaviour in key flow quantities. Specifically, the characteristic axial length of the penetrating vortices is shown to scale as $Re^{-0.22}$; the angular momentum transport scales as $Re^{0.42}$; the kinetic and elastic boundary-layer thicknesses based on angular velocity and hoop stress near the inner cylinder wall scale as $Re^{-0.48}$ and $Re^{-0.49}$, respectively. This implies that the viscoelastic Keplerian turbulence belongs to the classical turbulent regime of TCF with the Prandtl–Blasius-type boundary layer. Furthermore, we present an analytical relation between the viscous and elastic dissipation rates of kinetic energy and the angular momentum transport and in turn demonstrate its validity using our DNS data. This study has paved the way for future research to explore astrophysics-related Keplerian turbulence and angular momentum transport via the scaling relations of the analogous TCF of dilute polymeric solutions.

Direct numerical simulation (DNS) is performed to explore turbulent Rayleigh–Bénard convection in spherical shells. Our simulations cover six distinct values of radius ratio, $\eta = r_i/r_o = 0.2$, 0.3, 0.4, 0.5, 0.6 and 0.8, under the assumption of a centrally condensed mass with gravity profile $g \sim 1/r^{2}$; where $r_i$, $r_o$ and $r$ denote the inner shell radius, the outer shell radius and the local radial coordinate, respectively. The Prandtl number is kept constant at unity while the Rayleigh number ($Ra$) is varied from $3 \times 10^{3}$ to $5 \times 10^8$. Our primary aim is to analyze how the radius ratio influences the global transport properties and flow physics. To gain insights into the scaling behaviour of the Nusselt number ($Nu$) and the Reynolds number ($Re$) with respect to $Ra$ and $\eta$, we apply the Grossmann–Lohse (GL) theory (Grossmann & Lohse, J. Fluid Mech., vol. 407, 2000, pp. 27–56) to the system. It is observed that the scaling exponents for $Nu$ and $Re$ in relation to $Ra$ are more significant for smaller $\eta$ values, suggesting that the simulations with smaller $\eta$ reach the classical $Nu\sim Ra^{1/3}$ regime at a relatively lower $Ra$. This observation could also imply the systems with smaller $\eta$ might transition to the ultimate regime earlier at a smaller $Ra$. Based on our extensive DNS data, we establish that the thickness of the inner thermal boundary, $\lambda _{\vartheta }^{i}$, follows a scaling relationship of $\lambda _{\vartheta }^{i} \sim \eta ^{1/2}$. This relationship, in turn, leads to a scaling law for $Nu$ in the form of $Nu \sim f(\eta ) Ra^{\gamma }$, where the function $f(\eta )$ is defined as $f(\eta ) = {\eta ^{1/2}}/{(1+\eta ^{4/3})}$, and the exponent $\gamma$ depends on both $Ra$ and $\eta$. Additionally, we characterize and explain the asymmetry in the velocity field by introducing the separate Reynolds numbers for the inner and outer shells. The asymmetry of the kinetic and thermal energy dissipation rates in the inner and outer boundary layers (BLs) is also quantified.

Three-dimensional direct numerical simulations of rotating Rayleigh–Bénard convection in the planar geometry with no-slip top and bottom and periodic lateral boundary conditions are performed for a broad parameter range with the Rayleigh number spanning in $5\times 10^{6}\leq Ra \leq 5\times 10^{13}$, Ekman number within $5\times 10^{-9}\leq Ek \leq 5\times 10^{-5}$ and Prandtl number $Pr=1$. The thermal and Ekman boundary layer (BL) statistics, temperature drop within the thermal BL, interior temperature gradient and scaling behaviours of the heat and momentum transports (reflected in the Nusselt $Nu$ and Reynolds numbers $Re$) as well as the convective length scale are investigated across various flow regimes. The global and local momentum transports are examined via the $Re$ scaling derived from the classical theoretical balances of viscous–Archimedean–Coriolis (VAC) and Coriolis–inertial–Archimedean (CIA) forces. The VAC-based $Re$ scaling is shown to agree well with the data in the cellular and columnar regimes, where the characteristic convective length scales as the onset length scale ${\sim } Ek^{1/3}$, while the CIA-based $Re$ scaling and the inertia length scale $\sim (ReEk)^{1/2}$ work well in the geostrophic turbulence regime for $Ek\leq 1.5\times 10^{-8}$. The examinations of $Nu$, global and local $Re$, and convective length scale as well as the temperature drop within the thermal BL and its thickness scaling behaviours, indicate that for extreme parameters of $Ek\leq 1.5\times 10^{-8}$ and $80\lesssim RaEk^{4/3}\lesssim 200$, we have reached the diffusion-free geostrophic turbulence regime.

The geostrophic turbulence in rapidly rotating thermal convection exhibits characteristics shared by many highly turbulent geophysical and astrophysical flows. In this regime, the convective length and velocity scales and heat flux are all diffusion-free, i.e. independent of the viscosity and thermal diffusivity. Our direct numerical simulations (DNS) of rotating Rayleigh–Bénard convection in domains with no-slip top and bottom and periodic lateral boundary conditions for a fluid with the Prandtl number $Pr=1$ and extreme buoyancy and rotation parameters (the Rayleigh number up to $Ra=3\times 10^{13}$ and the Ekman number down to $Ek=5\times 10^{-9}$) indeed demonstrate all these diffusion-free scaling relations, in particular, that the dimensionless convective heat transport scales with the supercriticality parameter $\widetilde {Ra}\equiv Ra\, Ek^{4/3}$ as $Nu-1\propto \widetilde {Ra}^{3/2}$, where $Nu$ is the Nusselt number. We further derive and verify in the DNS that with the decreasing $\widetilde {Ra}$, the geostrophic turbulence regime undergoes a transition into another geostrophic regime, the convective heat transport in this regime is characterized by a very steep $\widetilde {Ra}$-dependence, $Nu-1\propto \widetilde {Ra}^{3}$.

Three-dimensional elastic turbulence in Taylor–Couette flows of dilute polymer solutions has been realized and thoroughly investigated via direct numerical simulations. A novel flow transition pathway from elastically dominated turbulence to solitary vortex pairs (or diwhirls) and eventually to elastic turbulence is observed by decreasing the fluid inertia ($Re$) over seven orders of magnitude, i.e. from $Re=1000$ to $0.0001$. The dominant spatio-temporal flow features in the elastic turbulence regime are those of large-scale unsteady diwhirls and small-scale axial and azimuthal travelling waves in the outer and inner halves of the gap, respectively. Moreover, it is conclusively shown that production of turbulent kinetic energy in purely elastic turbulence solely arises due to the stochastic nature of polymer stretch/relaxation. Overall, based on this comprehensive numerical investigation, the differences in the underlying fluid physics that give rise to turbulent fluctuations in elastically dominated and purely elastic turbulence have been delineated.

Direct numerical simulation of spanwise-rotation-driven flow transitions in viscoelastic plane Couette flow from a drag-reduced inertial to a drag-enhanced elasto-inertial turbulent flow state followed by full relaminarization is reported for the first time. Specifically, this novel flow transition begins with a drag-reduced inertial turbulent flow state at a low rotation number $0\leqslant Ro \leqslant 0.1$, and then transitions to a rotation/polymer-additive-driven drag-enhanced inertial turbulent regime, $0.1\leqslant Ro \leqslant 0.3$. In turn, the flow transitions to a drag-enhanced elasto-inertial turbulent state, $0.3\leqslant Ro \leqslant 0.9$, and eventually relaminarizes at $Ro=1$. In addition, two novel rotation-dependent drag enhancement mechanisms are proposed and substantiated. (1) The formation of large-scale roll cells results in enhanced convective momentum transport along with significant polymer elongation and stress generated in the extensionally dominated flow between adjacent roll cells at $Ro\leqslant 0.2$. (2) Coriolis-force-generated turbulent vortices cause strong incoherent transport and homogenization of significant polymer stress in the bulk via their vortical circulations at $Ro=0.5 - 0.9$.

A high-order transition route from inertial to elasticity-dominated turbulence (EDT) in Taylor–Couette flows of polymeric solutions has been discovered via direct numerical simulations. This novel two-step transition route is realized by enhancing the extensional viscosity and hoop stresses of the polymeric solution via increasing the maximum chain extension at a fixed polymer concentration. Specifically, in the first step inertial turbulence is stabilized to a laminar flow much like the modulated wavy vortex flow. The second step destabilizes this laminar flow state to EDT, i.e. a spatially smooth and temporally random flow with a $-3.5$ scaling law of the energy spectrum reminiscent of elastic turbulence. The flow states involved are distinctly different to those observed in the reverse transition route from inertial turbulence via a relaminarization of the flow to elasto-inertial turbulence in parallel shear flows, underscoring the importance of polymer-induced hoop stresses in realizing EDT that are absent in parallel shear flows.

The flow physics of inertio-elastic turbulent Taylor–Couette flow for a radius ratio of $0.5$ in the Reynolds number ($Re$) range of $500$ to $8000$ is investigated via direct numerical simulation. It is shown that as $Re$ is increased the turbulence dynamics can be subdivided into two distinct regimes: (i) a low $Re \leqslant 1000$ regime where the flow physics is essentially dominated by nonlinear elastic forces and the main contribution to transport and mixing of momentum, stress and energy comes from large-scale flow structures in the bulk region and (ii) a high $Re \geqslant 5000$ regime where inertial forces govern the flow physics and the flow dynamics is mainly governed by small-scale flow structures in the near-wall region. Flow–microstructure coupling analysis reveals that the elastic Görtler instability in the near-wall region is triggered via significant polymer extension and commensurately high hoop stresses. This instability gives rise to small-scale elastic vortical structures identified as elastic Görtler vortices which are present at all $Re$ considered. In fact, these vortices develop herringbone streaks near the inner wall that have a longer average life span than their Newtonian counterparts due to their elastic origin. Examination of the budgets of mean streamwise enstrophy, mean kinetic energy, turbulent kinetic energy and Reynolds shear stress demonstrates that increasing fluid inertia hinders the generation of elastic stresses, leading to a monotonic reduction of the elastic-related effects on the flow physics.

Direct numerical simulation of polymer-induced flow relaminarization of turbulent spanwise-rotating plane Couette flow (RPCF) is reported for the first time. Specifically, the reverse transition pathway from a Newtonian turbulent RPCF to a fully relaminarized drag enhanced viscoelastic flow has been elucidated. Evidently, this transition occurs gradually by weakening and eventual elimination of small-scale vortices as the Weissenberg number ($Wi$) is enhanced, paving the way for a two-dimensional laminar flow consisting of large-scale and highly organized roll cells. The influence of polymer additives on convective momentum exchange by large-scale roll cells and small-scale turbulent vortices, namely, the drag reduction (DR) realized by elimination of turbulent vortices and the significant drag enhancement (DE) that results from polymer roll cell interactions has been identified as the mechanism of DE. The observed vortical changes point to a universal mechanism for the coupling of polymer chains and turbulent vortices in wall-bounded viscoelastic DE and DR flows.

We report direct numerical simulation results that clearly elucidate the mechanism that leads to curvature dependence of drag enhancement (DE) in viscoelastic turbulent Taylor–Couette flow. Change in the angular momentum transport and its inherent link to transitions in vortical flow structures have been explored to depict the influence of the curvature of the flow geometry on DE. Specifically, it has been demonstrated that a transition in vortical structures with increasing radius ratio leads to weakening and elimination of the small-scale Görtler vortices and development and better organization (occupying the entire gap) of large-scale Taylor vortices as also evinced by the patterns of angular momentum current. The commensurate change in DE and its underlying mechanism are examined by contributions of convective flux and polymeric stress to the angular momentum current. The present finding paves the way for capturing highly localized elastic turbulence structures in direct numerical simulation by increasing geometry curvature in traditional turbulent curvilinear flows.