Developing a model to describe the shock-accelerated cylindrical fluid layer with arbitrary Atwood numbers is essential for uncovering the effect of Atwood numbers on the perturbation growth. The recent model (J. Fluid Mech., vol. 969, 2023, p. A6) reveals several contributions to the instability evolution of a shock-accelerated cylindrical fluid layer but its applicability is limited to cases with an absolute value of Atwood numbers close to $1$, due to the employment of the thin-shell correction and interface coupling effect of the fluid layer in vacuum. By employing the linear stability analysis on a cylindrical fluid layer in which two interfaces separate three arbitrary-density fluids, the present work generalizes the thin-shell correction and interface coupling effect, and thus, extends the recent model to cases with arbitrary Atwood numbers. The accuracy of this extended model in describing the instability evolution of the shock-accelerated fluid layer before reshock is confirmed via direct numerical simulations. In the verification simulations, three fluid-layer configurations are considered, where the outer and intermediate fluids remain fixed and the density of the inner fluid is reduced. Moreover, the mechanisms underlying the effect of the Atwood number at the inner interface on the perturbation growth are mainly elucidated by employing the model to analyse each contribution. As the Atwood number decreases, the dominant contribution of the Richtmyer–Meshkov instability is enhanced due to the stronger waves reverberated inside the layer, leading to weakened perturbation growth at initial in-phase interfaces and enhanced perturbation growth at initial anti-phase interfaces.

The efficacy of steady large-amplitude blowing/suction on instability and transition control for a hypersonic flat plate boundary layer with Mach number 5.86 is investigated systematically. The influence of the blowing/suction flux and amplitude on instability is examined through direct numerical simulation and resolvent analysis. When a relatively small flux is used, the two-dimensional instability critical frequency that distinguishes the promotion/suppression mode effect closely aligns with the synchronisation frequency. For the oblique wave, as the spanwise wavenumber increases, the suppression effects would become weaker and the mode suppression bandwidth diminishes/increases in general in the blowing/suction control. Increasing the blowing/suction flux can effectively broaden the frequency bandwidth of disturbance suppression. The influence of amplitude on disturbance suppression is weak in a scenario of constant flux. To gain a deeper insight into disturbance suppression mechanism, momentum potential theory (MPT) and kinetic energy budget analysis are further employed in analysing disturbance evolution with and without control. When the disturbance is suppressed, the blowing induces the transport of certain acoustic components along the compression wave out of the boundary layer, whereas the suction does not. The velocity fluctuations are derived from the momentum fluctuations of the MPT. Compared with the momentum fluctuations, the evolutions indicated by each component's velocity fluctuations greatly facilitate the investigations of the acoustic nature of the second mode. The rapid variation of disturbance amplitude near the blowing is caused by the oscillations of the acoustic component and phase speed differences between vortical and thermal components. Kinetic energy budget analysis is performed to address the non-parallel effect of the boundary layer introduced by blowing/suction, which tends to suppress disturbances near the blowing. Moreover, viscous effects leading to energy dissipation are identified to be stronger in regions where the boundary layer is rapidly thickening. Finally, it is demonstrated that a flat plate boundary layer transition triggered by a random disturbance can be delayed by a blowing/suction combination control. The resolvent analysis further demonstrates that disturbances with frequencies that dominate the early transition stage are dampened in the controlled base flow.

This study presents a comprehensive analysis on the extreme positive and negative events of wall shear stress and heat flux fluctuations in compressible turbulent boundary layers (TBLs) solved by direct numerical simulations. To examine the compressibility effects, we focus on the extreme events in two representative cases, i.e. a supersonic TBL of Mach number $M=2$ and a hypersonic TBL of $M=8$, by scrutinizing the coherent structures and their correlated dynamics based on conditional analysis. As characterized by the spatial distribution of wall shear stress and heat flux, the extreme events are indicated to be closely related to the structural organization of wall streaks, in addition to the occurrence of the alternating positive and negative structures (APNSs) in the hypersonic TBL. These two types of coherent structures are strikingly different, namely the nature of wall streaks and APNSs are shown to be related to the solenoidal and dilatational fluid motions, respectively. Quantitative analysis using a volumetric conditional average is performed to identify and extract the coherent structures that directly account for the extreme events. It is found that in the supersonic TBL, the essential ingredients of the conditional field are hairpin-like vortices, whose combinations can induce wall streaks, whereas in the hypersonic TBL, the essential ingredients become hairpin-like vortices as well as near-wall APNSs. To quantify the momentum and energy transport mechanisms underlying the extreme events, we proposed a novel decomposition method for extreme skin friction and heat flux, based on the integral identities of conditionally averaged governing equations. Taking advantage of this decomposition method, the dominant transport mechanisms of the hairpin-like vortices and APNSs are revealed. Specifically, the momentum and energy transports undertaken by the hairpin-like vortices are attributed to multiple comparable mechanisms, whereas those by the APNSs are convection dominated. In that, the dominant transport mechanisms in extreme events between the supersonic and hypersonic TBLs are indicated to be totally different.

Instability evolutions of shock-accelerated thin cylindrical SF$_6$ layers surrounded by air with initial perturbations imposed only at the outer interface (i.e. the ‘Outer’ case) or at the inner interface (i.e. the ‘Inner’ case) are numerically and theoretically investigated. It is found that the instability evolution of a thin cylindrical heavy fluid layer not only involves the effects of Richtmyer–Meshkov instability, Rayleigh–Taylor stability/instability and compressibility coupled with the Bell–Plesset effect, which determine the instability evolution of the single cylindrical interface, but also strongly depends on the waves reverberated inside the layer, thin-shell correction and interface coupling effect. Specifically, the rarefaction wave inside the thin fluid layer accelerates the outer interface inward and induces the decompression effect for both the Outer and Inner cases, and the compression wave inside the fluid layer accelerates the inner interface inward and causes the decompression effect for the Outer case and compression effect for the Inner case. It is noted that the compressible Bell model excluding the compression/decompression effect of waves, thin-shell correction and interface coupling effect deviates significantly from the perturbation growth. To this end, an improved compressible Bell model is proposed, including three new terms to quantify the compression/decompression effect of waves, thin-shell correction and interface coupling effect, respectively. This improved model is verified by numerical results and successfully characterizes various effects that contribute to the perturbation growth of a shock-accelerated thin heavy fluid layer.

A novel data-driven modal analysis method, reduced-order variational mode decomposition (RVMD), is proposed, inspired by the Hilbert–Huang transform and variational mode decomposition (VMD), to resolve transient or statistically non-stationary flow dynamics. First, the form of RVMD modes (referred to as an ‘elementary low-order dynamic process’, ELD) is constructed by combining low-order representation and the idea of intrinsic mode function, which enables the computed modes to characterize the non-stationary properties of space–time fluid flows. Then, the RVMD algorithm is designed based on VMD to achieve a low-redundant adaptive extraction of ELDs in flow data, with the modes computed by solving an elaborate optimization problem. Further, a combination of RVMD and Hilbert spectral analysis leads to a modal-based time-frequency analysis framework in the Hilbert view, providing a potentially powerful tool to discover, quantify and analyse the transient and non-stationary dynamics in complex flow problems. To provide a comprehensive evaluation, the computational cost and parameter dependence of RVMD are discussed, as well as the relations between RVMD and some classic modal decomposition methods. Finally, the virtues and utility of RVMD and the modal-based time-frequency analysis framework are well demonstrated via two canonical problems: the transient cylinder wake and the planar supersonic screeching jet.

The highly nonlinear evolution of the single-mode stratified compressible Rayleigh–Taylor instability (RTI) is investigated via direct numerical simulation over a range of Atwood numbers ($A_T=0.1$–$0.9$) and Mach numbers ($Ma=0.1$–$0.7$) for characterising the isothermal background stratification. After the potential stage, it is found that the bubble is accelerated to a velocity which is well above the saturation value predicted in the potential flow model. Unlike the bubble re-acceleration behaviour in quasi-incompressible RTI with uniform background density, the characteristics in the stratified compressible RTI are driven by not only vorticity accumulation inside the bubble but also flow compressibility resulting from the stratification. Specifically, in the case of strong stratification and high $A_T$, the flow compressibility dominates the bubble re-acceleration characters. To model the effect of flow compressibility, we propose a novel model to reliably describe the bubble re-acceleration behaviours in the stratified compressible RTI, via introducing the dilatation into the classical model that takes into account only vorticity accumulation.

Pressure fluctuations play an essential role in the transport of turbulent kinetic energy and vibrational loading. This study focuses on examining the effect of wall cooling on pressure fluctuations in compressible turbulent boundary layers by high-fidelity direct numerical simulations. Pressure fluctuations result from the vorticity mode and the acoustic mode that are both closely dependent on compressibility. To demonstrate the effects of wall cooling at various compressibility intensities, three free-stream Mach numbers are investigated, i.e. $M_\infty =0.5$, 2.0 and 8.0, with real gas effects being absent for $M_\infty =8.0$ due to a low enthalpy inflow. Overall, opposite effects of wall cooling on pressure fluctuations are found between the subsonic/supersonic cases and the hypersonic case. Specifically, the pressure fluctuations normalized by wall shear stress $p^\prime _{rms}/\tau _w$ are suppressed in the subsonic and supersonic cases, while enhanced in the hypersonic case near the wall. Importantly, travelling-wave-like alternating positive and negative structures (APNS), which greatly contribute to pressure fluctuations, are identified within the viscous sublayer and buffer layer in the hypersonic cases. Furthermore, generating mechanisms of pressure fluctuations are explored by extending the decomposition based on the fluctuating pressure equation to compressible turbulent boundary layers. Pressure fluctuations are decomposed into five components, in which rapid pressure, slow pressure and compressible pressure are dominant. The suppression of pressure fluctuations in the subsonic and supersonic cases is due to both rapid pressure and slow pressure being suppressed by wall cooling. In contrast, wall cooling strengthens compressible pressure for all Mach numbers, especially in the hypersonic case, resulting in increased wall pressure fluctuations. Compressible pressure plays a leading role in the hypersonic case, mainly due to the APNS. Essentially, the main effects of wall cooling can be interpreted by the suppression of the vorticity mode and the enhancement of the acoustic mode.

A direct numerical simulation database of a weakly compressible turbulent channel flow with bulk Mach number 1.56 is studied in detail, including the geometrical relationships between the pressure-Hessian tensor and the vorticity/strain-rate tensor, as well as the mechanism of the pressure-Hessian tensor contributing to the evolution of invariants of the velocity gradient tensor. The results show that the geometrical relationships between the pressure-Hessian tensor and the vorticity/strain-rate tensor in the central region of the channel are consistent with that of isotropic turbulence. However, in the buffer layer with relatively stronger inhomogeneity and anisotropy, the vorticity tends to be aligned with the first or second eigenvector of the pressure-Hessian tensor in the unstable focus/compressing topological region, and tends to be aligned with the first eigenvector of the pressure-Hessian tensor in the stable focus/stretching topological region. In the unstable node/saddle/saddle and stable node/saddle/saddle topological regions, the vorticity prefers to lie in the plane of the first and second eigenvectors of the pressure-Hessian tensor. The strain-rate and the pressure-Hessian tensors tend to share their second principal direction. Moreover, for the coupling between the pressure-Hessian tensor and the principal strain rates, we clarify the influence on dissipation, the nonlinear generation of dissipation and the enstrophy generation. The decomposition of the pressure-Hessian tensor further shows that the slow pressure-related term dominates the pressure-Hessian tensor's contribution, and the influence of inhomogeneity and anisotropy mainly originates from the inhomogeneity and anisotropy of the fluctuating velocity. These statistical properties would be instructive in formulating dynamical models of the velocity gradient tensor for wall turbulence.

Kinetic energy and enstrophy transfer in compressible Rayleigh–Taylor (RT) turbulence were investigated by means of direct numerical simulation. It is revealed that compressibility plays an important role in the kinetic energy and enstrophy transfer based on analyses of transport and large-scale equations. For the generation and transfer of kinetic energy, some findings have been obtained as follows. The pressure-dilatation work dominates the generation of kinetic energy in the early stage of flow evolution. The baropycnal work and deformation work handle the kinetic energy transfer from large to small scales on average for RT turbulence. The baropycnal work is mainly responsible for the kinetic energy transfer on large scales, and the deformation work for the kinetic energy transfer on small scales. The baropycnal work is also disclosed to be related to the compressibility from the finding that the expansion motion enhances the positive baropycnal work and the compression motion strengthens the negative baropycnal work. For the generation and transfer of enstrophy, the horizontal enstrophy is generated by the baroclinic effect and the vertical enstrophy by vortex stretching and tilting. Then the enstrophy is strengthened by the vortex stretching and tilting during the evolution of RT turbulence and the vorticity tends to be isotropic in the turbulent mixing region. The large-scale enstrophy equation in compressible flow has also been derived to deal with the enstrophy transfer. It is identified that the enstrophy is transferred from large to small scales on average and tends to stabilize for RT turbulence.

Locomotion of a passively flapping flat plate has been studied numerically by means of a multiblock lattice Boltzmann method. A flexible plate is modelled by a rigid plate with a torsion spring acting about the pivot at the leading edge of the plate. A dynamic model of this kind is called a lumped-torsional-flexibility model. When the leading edge is forced to heave sinusoidally, the plate pitches passively and propels itself in the horizontal direction as a result of the fluid–plate interaction. We have investigated various aspects of the mechanics behind the behaviour of the flapping plate, including the periodic- and non-periodic-flow states, the spontaneous motion of the plate, vortical structure and how they compare to similar propulsion systems in animals. In the periodic-flow regime, two dynamical responses of the passively pitching plate (forward and backward movements) are observed. Which movement will occur depends only on the frequency ratio F of the natural frequency of the system and the heaving frequency associated with the lumped torsional flexibility. It is found that the plate will select the forward movement when F > 1 and the backward movement when F ≤ 1. In the forward-movement regime, analysis of the dynamical behaviours and propulsive properties of the passively pitching plate indicates that the torsional flexibility can remarkably improve the propulsive performance. In addition, four kinds of vortex structures in the near wake are identified, which mainly depend on the forward speed of the plate. Finally the forward movement is compared to the flapping-based locomotion of swimming and flying animals. The results obtained in this study are consistent with the observations and measurements of swimming and flying animals; thus, they may provide physical insights into understanding of the propulsive mechanisms of the flapping wings and fins of animals.