Direct numerical simulations were performed to characterize fully developed supersonic turbulent channel flows over isothermal rough walls. The effect of roughness was incorporated using a level-set/volume-of-fluid immersed boundary method. Turbulence statistics of five channel flows are compared, including one reference case with both walls smooth and four cases with smooth top walls and bottom walls with two-dimensional (2-D) and three-dimensional (3-D) sinusoidal roughnesses. Results reveal a strong dependence of the turbulence on the roughness topography and the associated shock patterns. Specifically, the 2-D geometries generate strong oblique shock waves that propagate across the channel and are reflected back to the rough-wall side. These strong shocks are absent in the smooth-wall channel and are significantly weaker in cases with 3-D roughness geometries, replaced by weak shocklets. At the impingement locations of the shocks on the top wall in the 2-D roughness cases, localized augmentations of turbulence shear production are observed. Such regions of augmented production also exist for the 3-D cases, at a much weaker level. The oblique shock waves are thought to be responsible for a more significant entropy generation for cases with 2-D surfaces than those with 3-D ones, leading to a higher irreversible heat generation and consequently higher temperature values in 2-D roughness cases. In the present supersonic channels, the effects of roughness extend beyond the near-wall layer due to the shocks. This suggests that outer layer similarity may not fully apply to a rough-wall supersonic turbulent flow.

We consider the stability of laminar high-Reynolds-number flow through a planar channel formed by a rigid wall and a heavy compliant wall under longitudinal tension with motion resisted by structural damping. Numerical simulations indicate that the baseline state (with Poiseuille flow and a flat wall) exhibits two unstable normal modes: the Tollmien–Schlichting (TS) mode and a surface-based mode which manifests as one of two flow-induced surface instabilities (FISI), known as travelling wave flutter (TWF) and static divergence (SD), respectively. In the absence of wall damping the system is unstable to TWF, where the neutrally stable wavelength becomes shorter as the wall mass increases. With wall damping, TWF is restricted to long wavelengths through interaction with the most unstable centre mode, while for wall damping greater than a critical value the system exhibits an SD mode with a two branch neutral stability curve; the critical conditions along the upper and lower branches are constructed in the limit of large wall damping. We compute the Reynolds–Orr and activation energy descriptions of these neutrally stable FISI by continuing the linear stability analysis to the following order in perturbation amplitude. We find that both FISI are primarily driven by the working of normal stress on the flexible wall, lower-branch SD has negative activation energy, while upper-branch SD approaches zero activation energy in the limit of large wall damping. Finally, we elucidate interaction between TS and TWF modes for large wall mass, resulting in stable islands within unstable regions of parameter space.

We derive a scaling law for the characteristic frequencies of wall pressure fluctuations in swept shock wave/turbulent boundary layer interactions in the presence of cylindrical symmetry, based on analysis of a direct numerical simulations database. Direct numerical simulations in large domains show evidence of spanwise rippling of the separation line, with typical wavelength proportional to separation bubble size. Pressure disturbances around the separation line are shown to be convected at a phase speed proportional to the cross-flow velocity. This information is leveraged to derive a simple model for low-frequency unsteadiness, which extends previous two-dimensional models (Piponniau et al., J. Fluid Mech., vol. 629, 2009, pp. 87–108), and which correctly predicts growth of the typical frequency with the sweep angle. Inferences regarding the typical frequencies in more general swept shock wave/turbulent boundary layer interactions are also discussed.

The dynamics of charge-free dielectric droplets under electric fields in weakly conducting fluids is relevant to a large range of technological and environmental applications. In Sorgentone & Vlahovska (J. Fluid Mech., vol. 951, 2022, R2), the authors investigate by asymptotic theory and boundary integral simulations the motility states of a pair of dissimilar droplets at an arbitrary orientation to the electric field and find a state of steady motility that is carefully balanced between repulsion and contact.

Mixing of fluids in a coaxial jet is studied under four distinct viscosity ratios, $m=1$, $10$, $20$ and $40$, using highly resolved large-eddy simulations (LES), particle image velocimetry and planar laser-induced fluorescence. The accuracy of predictions is tested against data obtained by the simultaneous experimental measurements of velocity and concentration fields. For the highest and lowest viscosity ratios, standard RANS models with unclosed terms pertaining to viscosity variations are employed. We show that the standard Reynolds-averaged Navier–Stokes (RANS) approach with no explicit modelling for variable-viscosity terms is not applicable whereas dynamic LES models provide high-quality agreement with the measurements. To identify the underlying mixing physics and sources of discrepancy in RANS predictions, two distinct mixing modes are defined based on the viscosity ratio. Then, for each mode, the evolution of mixing structures, momentum budget analysis with emphasis on variable-viscosity terms, analysis of the turbulent activity and decay of turbulence are investigated using highly resolved LES data. The mixing dynamics is found to be quite distinct in each mixing mode. Variable viscosity manifests multiple effects that are working against each other. Viscosity gradients induce additional instabilities while increasing overall viscosity decreases the effective Reynolds number leading to laminarization of the turbulent jet, explaining the lack of dispersion and turbulent diffusion. Momentum budget analysis reveals that variable-viscosity terms are significant to be neglected. The scaling of the energy spectrum cascade suggests that in the TLL mode the unsteady laminar shedding is responsible for the eddies observed.

We study the effect of insoluble surfactants on the spatio-temporal evolution of turbulent jets. We use three-dimensional numerical simulations and employ an interface-tracking/level-set method that accounts for surfactant-induced Marangoni stresses. The present study builds on our previous work (Constante-Amores et al., J. Fluid Mech., vol. 922, 2021, A6) in which we examined in detail the vortex–surface interaction in the absence of surfactants. Numerical solutions are obtained for a wide range of Weber and elasticity numbers in which vorticity production is generated by surface deformation and surfactant-induced Marangoni stresses. The present work demonstrates, for the first time, the crucial role of Marangoni stresses, brought about by surfactant concentration gradients, in the formation of coherent, hairpin-like vortex structures. These structures have a profound influence on the development of the three-dimensional interfacial dynamics. We also present theoretical expressions for the mechanisms that influence the rate of production of circulation in the presence of surfactants for a general, three-dimensional, two-phase flow, and highlight the dominant contribution of surfactant-induced Marangoni stresses.

Active diffusion of substances in binary immiscible and incompressible fluids with different densities occurs universally in nature and industry, but relevant mathematical models and numerical simulation have been studied scarcely so far. In this paper, a thermodynamically consistent phase-field model is established to describe the activated solute transport in binary fluids with different densities. A mixed free-energy function for multiple solutes is proposed, which leads to different solute chemical potentials in binary solvent fluids, thus it has the ability to characterize the solubility difference of the solutes in two solvents. The two-phase flow is governed by a general hydrodynamic phase-field model that can account for general average velocity and different densities. The proposed model is derived rigorously using the second law of thermodynamics. Moreover, a general multi-component solute diffusion model is established using the Maxwell–Stefan approach, which involves the crossing influences between different solutes. To solve the model effectively, an efficient, linearized and decoupled numerical method is proposed for the model as well. The proposed numerical method can preserve the thermodynamical consistency, i.e. obeying an energy dissipation law at the discrete level, as well as guarantee the mass conservation law for the solutes and solvent fluids. Numerical experiments are carried out to show that the proposed model and numerical method can simulate various processes of the solute active diffusion in two-phase solvent fluids.

When a cylinder is mounted on an elastic support within a current, vortex-induced vibrations (VIV) may occur down to a Reynolds number (Re) close to $20$, based on the body diameter ($D$) and inflow velocity ($U$), i.e. below the critical value of $47$ reported for the onset of flow unsteadiness when the body is fixed. The impact of a forced rotation of the elastically mounted cylinder on the system behaviour is explored numerically for $Re \leqslant 30$, over wide ranges of values of the rotation rate (ratio between body surface velocity and $U$, $\alpha \in [0,5]$) and reduced velocity (inverse of the oscillator natural frequency non-dimensionalized by $D$ and $U$, $U^\star \in [2,30]$). The influence of the rotation is not monotonic, but the most prominent effect uncovered in this work is a substantial enhancement of the subcritical-Re, flow-induced vibrations beyond $\alpha =2$. This enhancement is twofold. First, the rotation results in a considerable expansion of the vibration/flow unsteadiness region in the $({Re},U^\star )$ domain, down to $Re=4$. Second, the elliptical orbits described by the rotating body are subjected to a major amplification, with a transition from VIV to responses whose magnitude tends to increase unboundedly with $U^\star$, even though still synchronized with flow unsteadiness. The emergence of such galloping-like oscillations close to the onset of vibrations disrupts the scenario of gradual vibration growth with Re, as amplitudes larger than $10$ body diameters may be observed at $Re=10$.

Shock-induced instability developments of two successive interfaces have attracted much attention, but remain a difficult problem to solve. The feedthrough and reverberating waves between two successive interfaces significantly influence the hydrodynamic instabilities of the two interfaces. The evolutions of two successive slow/fast interfaces driven by a weak shock wave are examined experimentally and numerically. First, a general one-dimensional theory is established to describe the movements of the two interfaces by studying the rarefaction waves reflected between the two interfaces. Second, an analytical, linear model is established by considering the arbitrary wavenumber and phase combinations and compressibility to quantify the feedthrough effect on the Richtmyer–Meshkov instability (RMI) of two successive slow/fast interfaces. The feedthrough significantly influences the RMI of the two interfaces, and even leads to abnormal RMI (i.e. phase reversal of a shocked slow/fast interface is inhibited) which is the first observational evidence of the abnormal RMI provided by the present study. Moreover, the stretching effect and short-time Rayleigh–Taylor instability or Rayleigh–Taylor stabilisation imposed by the rarefaction waves on the two interfaces are quantified considering the two interfaces’ phase reversal. The conditions and outcomes of the freeze-out and abnormal RMI caused by the feedthrough are summarised based on the theoretical model and numerical simulation. A specific requirement for the simultaneously freeze-out of the instability of the two interfaces is proposed, which can potentially be used in the applications to suppress the hydrodynamic instabilities.

Recent works on wall-bounded flows have corroborated the coexistence of wall-attached eddies, whose statistical features are predicted through Townsend's attached-eddy hypothesis (AEH), and very-large-scale motions (VLSMs). Furthermore, it has been shown that the presence of wall-attached eddies within the logarithmic layer is linked to the appearance of an inverse-power-law region in the streamwise velocity energy spectra, upon significant separation between outer and viscous scales. In this work, a near-neutral atmospheric surface layer is probed with wind light detection and ranging to investigate the contributions to the streamwise velocity energy associated with wall-attached eddies and VLSMs for a very-high-Reynolds-number boundary layer. Energy and linear coherence spectra (LCS) of the streamwise velocity are interrogated to identify the spectral boundaries associated with eddies of different typologies. Inspired by the AEH, an analytical model for the LCS associated with wall-attached eddies is formulated. The experimental results show that the identification of the wall-attached-eddy energy contribution through the analysis of the energy spectra leads to an underestimate of the associated spectral range, maximum height attained and turbulence intensity. This feature is due to the overlap of the energy associated with VLSMs obscuring the inverse-power-law region. The LCS analysis estimates wall-attached eddies with a streamwise/wall-normal ratio of about 14.3 attaining a height of about 30 % of the outer scale of turbulence.

Using the theory of wave turbulence for rapidly rotating incompressible fluids derived by Galtier (Phys. Rev. E, vol. 68, 2003, 015301), we find the locality conditions that the solutions of the kinetic equation must satisfy. We show that the exact anisotropic Kolmogorov–Zakharov spectrum satisfies these conditions, which justifies the existence of this constant (positive) energy flux solution. Although a direct cascade is predicted in the transverse ($\perp$) and parallel ($\parallel$) directions to the rotation axis, we show numerically that in the latter case some triadic interactions can have a negative contribution to the energy flux, while in the former case all interactions contribute to a positive flux. Neglecting the parallel energy flux, we estimate the Kolmogorov constant at $C_K \simeq 0.749$. These results provide theoretical support for recent numerical and experimental studies.

Geophysical mass flows such as debris flows, dense pyroclastic flows and snow avalanches can self-channelize on shallow slopes. The confinement afforded by formed levees helps to maintain the flow depth, and hence mobility, allowing self-channelized flows to run out significantly farther than unconfined, spreading flows. Levee formation and self-channelization are strongly associated with particle-size segregation, but can also occur in monodisperse flows. This paper uses the monodisperse depth-averaged theory of Rocha et al. (J. Fluid Mech., vol. 876, 2019, pp. 591–641), which incorporates a hysteretic friction law and second-order depth-averaged viscous terms. Both of these are vital for the formation of a travelling wave that progressively deposits a pair of levees just behind the front. The three-dimensional velocity field is reconstructed in a frame moving with the front assuming Bagnold flow. This enables a bidisperse particle-size segregation theory to be used to solve for the large and small particle concentrations and particle paths in three-dimensions, for the first time. The model shows that the large particles tend to segregate to the surface of the flow, forming a carapace that extends over the centre of the channel, as well as along the external sides and base of the levee walls. The small particles segregate downwards, and are concentrated in the main channel and in the inner levee walls. This supports the contention that a low-friction channel lining provides a secondary mechanism for run-out enhancement. It is also shown that the entire theory scales with particle diameter, so experiments with millimetre-sized particles provide important insights into geophysical-scale flows with boulders and smaller rock fragments. The model shows that self-channelization does not need particle-size segregation to occur, but supports the hypothesis that particle-size segregation and the associated frictional feedback can significantly enhance both the flow mobility and the levee strength.