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.

The existence of a maximum drag enhancement (MDE) asymptote at high rotation ($Ro$) and Weissenberg ($Wi$) numbers in turbulent viscoelastic spanwise-rotating plane Couette flow has been demonstrated. Specifically, it is shown that above a critical $Wi$, drag enhancement plateaus and the MDE asymptote is realized in a broad range of $Ro$. The mean velocity profiles at MDE appear to closely follow a log-law profile that has a nearly identical slope but different intercepts as a function of $Ro$. Much like the maximum drag reduction (MDR) asymptote, the logarithmic function in MDE is closely followed if the mean velocity is plotted using the traditional inner variable scaling; however, the logarithmic function is not well defined when examined by the indicator function. Hence, in this study, we have used the logarithmic fit as a visual guide for the mean velocity profile. Last and perhaps the most intriguing finding of this study is that MDE occurs in the elasto-inertial turbulence (EIT) flow state; hence, it is mainly sustained by elastic forces much like the MDR flow state. To that end, a universal picture of elastically induced drag modification asymptotes is emerging, namely these asymptotic states are an inherent property of the elastically sustained EIT flow state.

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 present an analytical model of nonlinear evolution of two-dimensional single-mode Rayleigh–Taylor instability (RTI) in cylindrical geometry at arbitrary Atwood number for the first time. Our model covers a full scenario of bubble evolution from the earlier exponential growth to the nonlinear regime with the bubbles growing in time as $\frac {1}{2}a_{b}t^2$ for cylindrical RTI, other than as $V_{b}t$ for planar RTI, where $a_{b}$ and $V_{b}$ are the bubble acceleration and velocity, respectively. It is found that from this model the saturating acceleration $a_{b}$ is formulated as a simplified function of the external acceleration, Atwood number and number of perturbation waves. This model's predictions are in good agreement with data from direct numerical simulations.

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.

Drag reduction (DR) in plane Couette flow (PCF) induced by the addition of flexible polymers has been studied via direct numerical simulation (DNS). The similarities and differences in the drag reduction features of PCF and plane Poiseuille flow (PPF) have been examined in detail, particularly in regard to the polymer-induced modification of large-scale structures (LSSs) in the near-wall turbulence. Specifically, it has been demonstrated that in the near-wall region, drag-reduced PCF has features similar to those of drag-reduced PPF; however, in the core region, intriguing differences are found between these two drag-reduced shear flows. Chief among these differences is the significant polymer stretch that arises from the enhanced exchanges between elastic potential energy and turbulent kinetic energy and the commensurate observation of peak values of the conformation tensor components $\unicode[STIX]{x1D60A}_{yy}$ and $\unicode[STIX]{x1D60A}_{zz}$ in this region. This finding is in stark contrast to that of drag-reduced PPF where the polymer stretch and the exchanges between elastic potential energy and turbulent kinetic energy in the core region are insignificant; to this end, in drag-reduced PPF, peak values of the conformation tensor components appear in the near-wall region. Therefore, this study paves the way for understanding the underlying flow physics in drag-reduced PCF, particularly in the context of elastic theory. Moreover, the longitudinal large-scale streaks at the channel centre of drag-reduced PCF are greatly strengthened due to the increased production/dissipation ratio; the LSS imprint effects on the near-wall flow of drag-reduced PCF monotonically increase as the Weissenberg number is enhanced.

Numerical investigation of the underexpanded sonic coaxial jets is carried out using large eddy simulation for three typical inner nozzle lip-thicknesses. Various fundamental mechanisms dictating the flow phenomena including shock structure, shear layer evolution and sound production are investigated. It is found that the inner nozzle lip induces a recirculation zone between inner and outer jets, which significantly influences the behaviors of shock structures and shear layers. The sound properties of the coaxial jets are further analyzed in detail. As the inner lip-thickness increases, the helical screech mode switches to an axisymmetric one and high-frequency screech also occurs with an oscillation frequency of recirculation zone. Based on the temporal Fourier transform and correlation analysis, the primary sources of low- and high-frequency screeches are associated with the downstream shock cells in the jet column and the secondary shock structures in the outer annular jet, respectively. The proper orthogonal decomposition analysis reveals that the dominant structures constructed by the most energetic modes shift from the downstream shock cells region to the upstream secondary shock region as the lip-thickness increases. The results obtained in this study provide physical insight into the understanding of the mechanisms relevant to the coherent structures and sound properties in sonic coaxial jets.

Direct numerical simulation (DNS) of elastically induced turbulent flows has posed great challenges to researchers engaged in developing first-principle models and simulations that can predict faithfully the complex spatio-temporal dynamics of polymeric flows. To this end, DNS of elastically induced turbulent flow states in the Taylor–Couette (TC) flow are reported here with the aim of paving the way for a mechanistic understanding of this new class of flows. Specifically, the DNS not only faithfully reproduce the key feature of elastically induced turbulent flows, namely, substantial excitation of fluid motion at the smallest temporal and spatial scales, but also for the first time demonstrate the existence of three distinct flow regions in the gap for the inertio-elastic turbulence state: (i) a fluid-inertia (or outflow-) dominated inner-wall region; (ii) a fluid-elasticity (or inflow-) dominated outer-wall region; and (iii) an inflow/outflow core region. Based on this observation, a simple mechanism for the inertio-elastic turbulence in the TC flow has been postulated.