We study the decay of compressible magnetohydrodynamic (MHD) turbulence, emphasizing exchanges of energy between compressive and incompressive kinetic energies, magnetic energy, and thermal energy. A recently developed high order finite difference code is employed for compressible runs with a Mach number up to 2. Varying the nature of the initial conditions (magnitudes of velocity and magnetic fluctuations), and initial Mach numbers permits examination of various dynamical regimes characterized here by the changes between different energy reservoirs. Acoustic waves are responsible for the oscillatory exchange between compressive kinetic and thermal energy through the pressure dilatation term. The results indicate that exchange between kinetic and magnetic energy is dominated by interactions involving the solenoidal velocity. Several systematic rapid adjustments are found to be reproducible with simple scalings derived from the empirical data.

The Ginzburg–Landau equation (GLE) can phenomenologically model several key features of non-equilibrium systems including those in fluid mechanics. Its validity in real flows, however, remains questionable. Here, we show that the linear GLE can be formulated such that it has the same Wentzel–Kramers–Brillouin (WKB) approximation as for the linear global stability problem in open shear flows. We use the GLE to model the linear global modes of three different wakes and find that it can accurately capture the linear growth rate and frequency to first order in the WKB approximation. Furthermore, we find that it can also provide the shapes of the direct and adjoint eigenvectors and the regions of maximal structural sensitivity. The proposed model requires only the basic flow as input, but gives robust predictions and is computationally inexpensive. As well as opening up new possibilities for GLE-based control strategies, the proposed model makes accurate stability calculations possible, even for some computationally intractable open shear flows.

The characteristics and dynamics of the uniform-momentum zones (UMZ) and UMZ interfaces in a fully developed turbulent pipe flow are studied using direct numerical simulation at $Re_{\unicode[STIX]{x1D70F}}=500$ . The multiple UMZs detected from the probability density functions of the instantaneous streamwise velocity following de Silva et al. (J. Fluid Mech., vol. 786, 2016, pp. 309–331) showed similarities to both turbulent channel and boundary layer flows (TBL): the hierarchical structural distribution of thinner UMZs with thinner interfaces nearer the wall, accompanied with sharper and larger jumps in the streamwise velocity at the UMZ interface. The conditional average results indicate that channel and pipe are very similar quantitatively whereas pipe and TBL display significant discrepancies. The innermost UMZs in pipe flow exhibit different behaviours to the other UMZs in pipes. The contortion of the UMZ interface representing the meandering of coherent motions with high- and low-momentum streaks is examined three-dimensionally. The meandering of UMZ in both two and three dimensions intensifies away from the wall and is always wavier in the azimuthal direction than the streamwise direction. The UMZs in the near-wall region capture the small-scale velocity fluctuation of the near-wall cycle and show asymmetric modulation of $Q2$ ejections over $Q4$ sweeps. The asymmetric modulation of ejections over sweeps decreases from the wall towards the pipe centre and the opposite trend of elevated $Q4$ sweeps is observed for the innermost UMZs. Near the wall, the ejection regions are very spiky compared to the flat sweep regions whereas, in the pipe centre, the large-scale ejections are relatively flat and the sweep regions are spikier.

The effects of flow topology on the subgrid-scale (SGS) kinetic energy flux in compressible isotropic turbulence is studied. The eight flow topological types based on the three invariants of the filtered velocity gradient tensor are analysed at different scales, along with their roles in the magnitude and direction of kinetic energy transfer. The unstable focus/compressing (UFC), unstable node/saddle/saddle (UN/S/S) and stable focus/stretching (SFS), are the three predominant topological types at all scales; they account for at least 75 % of the flow domain. The UN/S/S and SFS types make major contributions to the average SGS flux of the kinetic energy from large scales to small scales in the inertial range. The unstable focus/stretching (UFS) topology makes a contribution to the reverse SGS flux of kinetic energy from small scales to large scales. In strong compression regions, the average contribution of the stable node/saddle/saddle (SN/S/S) topology to the SGS kinetic energy flux is positive and is predominant over those of other flow topologies. In strong expansion regions, the UFS topology makes a major contribution to the reverse SGS flux of the kinetic energy. As the turbulent Mach number increases, the increase of volume fraction of the UFS topological regions leads to the increase of the SGS backscatter of kinetic energy. The SN/S/S topology makes a dominant contribution to the direct SGS flux of the compressible component of the kinetic energy, while the UFS topology makes a dominant contribution to the reverse SGS flux of the compressible component of the kinetic energy.

In self-excited combustion systems, the application of open-loop forcing is known to be an effective strategy for controlling periodic thermoacoustic oscillations, but it is not known whether and under what conditions such a strategy would work on thermoacoustic oscillations that are not simply periodic. In this study, we experimentally examine the effect of periodic acoustic forcing on a prototypical thermoacoustic system consisting of a ducted laminar premixed flame oscillating quasiperiodically on an ergodic $\mathbb{T}^{2}$ torus at two incommensurate natural frequencies, $f_{1}$ and  $f_{2}$ . Compared with that of a classical period-1 system, complete synchronization of this $\mathbb{T}_{1,2}^{2}$ system is found to occur via a more intricate route involving three sequential steps: as the forcing amplitude, $\unicode[STIX]{x1D716}_{f}$ , increases at a fixed forcing frequency, $f_{f}$ , the system transitions first (i) to ergodic $\mathbb{T}_{1,2,f}^{3}$ quasiperiodicity; then (ii) to resonant $\mathbb{T}_{1,f}^{2}$ quasiperiodicity as the weaker of the two natural modes, $f_{2}$ , synchronizes first, leading to partial synchronization; and finally (iii) to a $P1_{f}$ limit cycle as the remaining natural mode, $f_{1}$ , also synchronizes, leading to complete synchronization. The minimum $\unicode[STIX]{x1D716}_{f}$ required for partial and complete synchronization decreases as $f_{f}$ approaches either $f_{1}$ or $f_{2}$ , resulting in two primary Arnold tongues. However, when forced at an amplitude above that required for complete synchronization, the system can transition out of $P1_{f}$ and into $\mathbb{T}_{1,2,f}^{3}$ or $\mathbb{T}_{2,f}^{2}$ . The optimal control strategy is to apply off-resonance forcing at a frequency around the weaker natural mode ( $f_{2}$ ) and at an amplitude just sufficient to cause $P1_{f}$ , because this produces the largest reduction in thermoacoustic amplitude via asynchronous quenching. Analysis of the Rayleigh index shows that this reduction is physically caused by a disruption of the positive coupling between the unsteady heat release rate of the flame and the $f_{1}$ and $f_{2}$ acoustic modes. If the forcing is applied near the stronger natural mode ( $f_{1}$ ), however, resonant amplification can occur. We then phenomenologically model this $\mathbb{T}_{1,2}^{2}$ thermoacoustic system as two reactively coupled van der Pol oscillators subjected to external sinusoidal forcing, and find that many of its synchronization features – such as the three-step route to $P1_{f}$ , the double Arnold tongues, asynchronous quenching and resonant amplification – can be qualitatively reproduced. This shows that these features are not limited to our particular system, but are universal features of forced self-excited oscillators. This study extends the applicability of open-loop control from classical period-1 systems with just a single time scale to ergodic $\mathbb{T}^{2}$ quasiperiodic systems with two incommensurate time scales.

In this study, the behaviours of subgrid-scale (SGS) turbulence are investigated with direct numerical simulations when an isotropic turbulence is brought to interact with imposed rapid waves. A partition of the velocity field is used to decompose the SGS stress into three parts, namely, the turbulent part $\unicode[STIX]{x1D749}^{T}$, the wave-induced part $\unicode[STIX]{x1D749}^{W}$ and the cross-interaction part $\unicode[STIX]{x1D749}^{C}$. Under strong wave straining, $\unicode[STIX]{x1D749}^{T}$ is found to follow the Kolmogorov scaling $\unicode[STIX]{x1D6E5}_{c}^{2/3}$, where $\unicode[STIX]{x1D6E5}_{c}$ is the filter width. Based on the linear Airy wave theory, $\unicode[STIX]{x1D749}^{W}$ and the filtered strain-rate tensor due to the wave motion, $\tilde{\unicode[STIX]{x1D64E}}^{W}$, are found to have different phases, posing a difficulty in applying the usual eddy-viscosity model. On the other hand, $\unicode[STIX]{x1D749}^{T}$ and the filtered strain-rate tensor due to the turbulent motion, $\tilde{\unicode[STIX]{x1D64E}}^{T}$, are only weakly wave-phase-dependent and could be well related by an eddy-viscosity model. The linear wave theory is also used to describe the vertical distributions of SGS statistics driven by the wave-induced motion. The predictions are in good agreement with the direct numerical simulation results. The budget equation for the turbulent SGS kinetic energy shows that the transport terms related to turbulence are important near the free surface and they compensate the imbalance between the energy flux and the SGS energy dissipation.

Based on recent numerical simulations and field experiments, the mechanism behind wake meandering is increasingly accepted to be through the amplification of upstream disturbances owing to the convectively unstable nature of the flow. In this paper, we deduce a low-order phenomenological model for the far-wake region, which is based on a modified form of the complex Ginzburg–Landau (CGL) equation for flows that are in the amplifier regime, i.e. are only convectively unstable. The model reproduces the main qualitative features of wake meandering: (i) its origin via amplification of upstream structures, (ii) dependence of oscillation frequency on the upstream disturbance amplitude (higher amplitudes lead to lower frequencies), (iii) shift towards lower frequencies as the wake flow evolves in the streamwise direction and, to an extent, (iv) the transfer of energy from very low frequencies towards relatively higher frequencies. Additionally, the model also predicts the increase in the meandering amplitude and an advancement in its onset with increasing thrust coefficient. To our knowledge, this is the first low-order dynamical system in the literature that models wake meandering. The model coefficients are obtained from the mean flow local stability results that we show correctly account for the changing operating conditions and thus pave way for the prediction of wake meandering features. Its low order makes it suitable to use inside an energy farm design model, where it can help to mitigate the adverse effects of wake meandering.

In an earlier paper (Wan et al., J. Fluid Mech., vol. 697, 2012, pp. 296–315), the authors showed that a similarity solution for anisotropic incompressible three-dimensional magnetohydrodynamic (MHD) turbulence, in the presence of a uniform mean magnetic field $\boldsymbol{B}_{0}$ , exists if the ratio of parallel to perpendicular (with respect to $\boldsymbol{B}_{0}$ ) similarity length scales remains constant in time. This conjecture appears to be a rather stringent constraint on the dynamics of decay of the energy-containing eddies in MHD turbulence. However, we show here, using direct numerical simulations, that this hypothesis is indeed satisfied in incompressible MHD turbulence. After an initial transient period, the ratio of parallel to perpendicular length scales fluctuates around a steady value during the decay of the eddies. We show further that a Taylor–Kármán-like similarity decay holds for MHD turbulence in the presence of a mean magnetic field. The effect of different parameters, including Reynolds number, mean field strength, and cross-helicity, on the nature of similarity decay is discussed.

Cascades of temperature and entropy fluctuations are studied by numerical simulations of stationary three-dimensional compressible turbulence with a heat source. The fluctuation spectra of velocity, compressible velocity component, density and pressure exhibit the $-5/3$ scaling in an inertial range. The strong acoustic equilibrium relation between spectra of the compressible velocity component and pressure is observed. The $-5/3$ scaling behaviour is also identified for the fluctuation spectra of temperature and entropy, with the Obukhov–Corrsin constants close to that of a passive scalar spectrum. It is shown by Kovasznay decomposition that the dynamics of the temperature field is dominated by the entropic mode. The average subgrid-scale (SGS) fluxes of temperature and entropy normalized by the total dissipation rates are close to 1 in the inertial range. The cascade of temperature is dominated by the compressible mode of the velocity field, indicating that the theory of a passive scalar in incompressible turbulence is not suitable to describe the inter-scale transfer of temperature in compressible turbulence. In contrast, the cascade of entropy is dominated by the solenoidal mode of the velocity field. The different behaviours of cascades of temperature and entropy are partly explained by the geometrical properties of SGS fluxes. Moreover, the different effects of local compressibility on the SGS fluxes of temperature and entropy are investigated by conditional averaging with respect to the filtered dilatation, demonstrating that the effect of compressibility on the cascade of temperature is much stronger than on the cascade of entropy.

Kinetic energy transfer in compressible isotropic turbulence is studied using numerical simulations with solenoidal forcing at turbulent Mach numbers ranging from 0.4 to 1.0 and at a Taylor Reynolds number of approximately 250. The pressure dilatation plays an important role in the local conversion between kinetic energy and internal energy, but its net contribution to the average kinetic energy transfer is negligibly small, due to the cancellation between compression and expansion work. The right tail of probability density function (PDF) of the subgrid-scale (SGS) flux of kinetic energy is found to be longer at higher turbulent Mach numbers. With an increase of the turbulent Mach number, compression motions enhance the positive SGS flux, and expansion motions enhance the negative SGS flux. Average of SGS flux conditioned on the filtered velocity divergence is studied by numerical analysis and a heuristic model. The conditional average of SGS flux is shown to be proportional to the square of filtered velocity divergence in strong compression regions for turbulent Mach numbers from 0.6 to 1.0. Moreover, the antiparallel alignment between the large-scale strain and the SGS stress is observed in strong compression regions. The inter-scale transfer of solenoidal and compressible components of kinetic energy is investigated by Helmholtz decomposition. The SGS flux of solenoidal kinetic energy is insensitive to the change of turbulent Mach number, while the SGS flux of compressible kinetic energy increases drastically as the turbulent Mach number becomes larger. The compressible mode persistently absorbs energy from the solenoidal mode through nonlinear advection. The kinetic energy of the compressible mode is transferred from large scales to small scales through the compressible SGS flux, and is dissipated by viscosity at small scales.

Turbulence is ubiquitous in nature and engineering applications. Although Kolmogorov’s (C. R. Acad. Sci. URSS, vol. 30, 1941a, pp. 301–305; Dokl. Akad. Nauk URSS, vol. 30, 1941b, pp. 538–540) theory suggested a unique turbulent state for high Reynolds numbers, multiple states were reported for several flow problems, such as Rayleigh–Bénard convection and Taylor–Couette flows. In this paper, we report that multiple states also exist for turbulent plane Couette flow with spanwise rotation through direct numerical simulations at rotation number $Ro=0.2$ and Reynolds number $Re_{w}=1300$ based on the angular velocity in the spanwise direction and half of the wall velocity difference. With two different initial flow fields, our results show that the flow statistics, including the mean streamwise velocity and Reynolds stresses, show different profiles. These different flow statistics are closely related to the flow structures in the domain, where one state corresponds to two pairs of roll cells, and the other shows three pairs. The present result enriches the studies on multiple states in turbulence.

We argue that the hypothesis of preservation of shape of dimensionless second- and third-order correlations during decay of incompressible homogeneous magnetohydrodynamic (MHD) turbulence requires, in general, at least two independent similarity length scales. These are associated with the two Elsässer energies. The existence of similarity solutions for the decay of turbulence with varying cross-helicity implies that these length scales cannot remain in proportion, opening the possibility for a wide variety of decay behaviour, in contrast to the simpler classic hydrodynamics case. Although the evolution equations for the second-order correlations lack explicit dependence on either the mean magnetic field or the magnetic helicity, there is inherent implicit dependence on these (and other) quantities through the third-order correlations. The self-similar inertial range, a subclass of the general similarity case, inherits this complexity so that a single universal energy spectral law cannot be anticipated, even though the same pair of third-order laws holds for arbitrary cross-helicity and magnetic helicity. The straightforward notion of universality associated with Kolmogorov theory in hydrodynamics therefore requires careful generalization and reformulation in MHD.