This work investigates the compressible turbulence induced by Richtmyer–Meshkov (RM) instability using high-resolution Navier–Stokes simulations. Special attention is paid to the characteristics of RM turbulence including the mixing width growth, the turbulent kinetic energy (TKE) decay, the mixing degree, inhomogeneity and anisotropy. Three distinct initial perturbation spectra are designed at the interface to reveal the initial condition imprint on the RM turbulence. Results show that cases with initial large-scale perturbations present a stronger imprint on statistical characteristics and also a quicker growth of mixing width, whereas cases with small-scale perturbations present a faster TKE decay, greater mixing level, higher isotropy and homogeneity. A thorough analysis on the inter-scale energy transfer in RM turbulence is also presented with the coarse-graining approach that exposes the two subfilter-scale (SFS) energy fluxes (i.e. deformation work and baropycnal work). A strong correlation between the nonlinear model of baropycnal work (Fluids, 4(2), 2019) and the simulation results is confirmed for the first time, demonstrating its potential in modelling the RM turbulence. Two primary mechanisms of baropycnal work (the straining and baroclinic generation processes) are explored with this nonlinear model. The evolutions of two SFS energy fluxes exhibit distinct behaviours at various filter scales, in different flow regions and under various flow motions (strain and rotation). It is found that all three cases share the common inter-scale energy transfer dynamics, which is important for modelling the RM turbulence.

We report an experimental study of the formation and evolution of laminar thermal structures generated by a small heat source, with a focus on their correlation to the thermal boundary layer and effects of heating time $t_{heat}$. The experiments are performed over the flux Rayleigh number ($Ra_f$) range $2.1\times 10^6 \leq Ra_f \leq 3.6\times 10^{7}$ and the Prandtl number ($Pr$) range $28.6 \leq Pr \leq 904.7$. The corresponding Rayleigh number ($Ra= t_{heat}\,Ra_{f}/\tau _0\,Pr$) range is $900 \leq Ra \leq 4\times 10^{4}$, where $\tau _0$ is a diffusion time scale. For thermal structures generated by continuous heating (i.e. starting plumes), their formation process exists three characteristic times that are well reflected by changes in the thermal boundary layer thickness. These characteristic times, denoted as $t_{emit}$, $t_{recover}$ and $t_{static}$, correspond to the moments when the plume emission begins and completes, and when the thermal boundary layer becomes quasi-static, respectively. Their $Ra_f$–$Pr$ dependencies are found to be $t_{emit}/\tau _0\sim Ra_f^{-0.41}\,Pr^{0.41}$, $t_{recover}/\tau _0\sim Ra_f^{-0.48}\,Pr^{0.48}$ and $t_{static}/\tau _0\sim Ra_f^{-0.49}\,Pr^{0.33}$, respectively. Thermal structures generated by finite $t_{heat}$ exhibit similar evolution dynamics once $t_{heat} \ge t_{emit}$, with the accelerating stage behaving like starting plumes and the decay stage like thermals (i.e. a finite amount of buoyant fluids). It is further found that their maximum rising velocity experiences a transition in the $Ra$-dependence from $Ra$ to $(Ra\ln Ra)^{0.5}$ at $Ra \simeq 6000$; and their maximum acceleration reaches the value of starting plumes at $t_{heat}\simeq t_{recover}$, and remains unchanged for larger $t_{heat}$. In particular, the maximum rising velocity for the cases with $t_{heat} = t_{recover}$ follows a scaling relation $Ra_f^{0.37}\,Pr^{-0.37}$, in contrast to the relation $Ra_f^{0.48}\,Pr^{-0.48}$ for starting plumes. This study provides a more comprehensive understanding of laminar thermal structures, which are relevant to a range of processes in nature and laboratory systems such as Rayleigh–Bénard convection.

We present a systematic simulation campaign to investigate the pairwise interaction of two mobile, monodisperse particles submerged in a viscous fluid and subjected to monochromatic oscillating flows. To this end, we employ the immersed boundary method to geometrically resolve the flow around the two particles in a non-inertial reference frame. We neglect gravity to focus on fluid–particle interactions associated with particle inertia and consider particles of three different density ratios aligned along the axis of oscillation. We systematically vary the initial particle distance and the frequency based on which the particles show either attractive or repulsive behaviour by approaching or moving away from each other, respectively. This behaviour is consistently confirmed for the three density ratios investigated, although particle inertia dictates the overall magnitude of the particle dynamics. Based on this, threshold conditions for the transition from attraction to repulsion are introduced that obey the same power law for all density ratios investigated. We furthermore analyse the flow patterns by suitable averaging and decomposition of the flow fields and find competing effects of the vorticity induced by the fluid–particle interactions. Based on these flow patterns, we derive a circulation-based criterion that provides a quantitative measure to categorize the different cases. It is shown that such a criterion provides a consistent measure to distinguish the attractive and repulsive arrangements.

Direct numerical simulations of hypersonic turbulent Couette flows are performed for top-wall Mach numbers of 6, 7 and 8, inspired by non-reactive high-enthalpy wind tunnel free-stream conditions, with the goal of analysing the physical processes driving the sub-filter-scale (SFS) stresses to inform development of large-eddy simulation techniques for hypersonic wall-bounded flows. Semi-local scaling laws collapse mean profiles and second-order turbulent statistics well, in spite of the strong wall-normal gradients of temperature and density. On the other hand, the SFS shear stresses exhibit an unexpected profile characterized by a region of pronounced shear stress deficit, which becomes more pronounced for higher Mach numbers. Instantaneous visualizations suggest that the SFS shear stress deficit is induced by the counter-gradient resolved momentum transport driven by residual velocity motions at the interface between high-density low-speed streaks being ejected away from the wall, and low-density high-speed ones replacing the displaced fluid, qualifying this as a compressibility effect. It is shown that the SFS shear stresses are primarily driven by second-order interactions between residual velocities, in spite of their triply nonlinear nature. This, in turn, motivated a statistical quadrant analysis revealing the presence of a SFS shear stress deficit, that is, SFS processes driving momentum transport of the resolved field towards the top wall. Additionally, higher-Reynolds-number simulations reveal that there is an upper limit of spatial filter widths to resolve large-scale structures, and such deficit is observable at any Reynolds number when a reasonable spatial filter is applied.

Exact steady eastward-moving vortex-dipole solutions, the equatorial modons, are constructed in the asymptotic limit of low divergence and small temperature variations in the thermal rotating shallow water (TRSW) model on the equatorial beta-plane. This regime is known to be relevant for the tropical atmosphere. The model itself is a generalization, allowing for horizontal temperature gradients, of the classical rotating shallow water model. The asymptotic modons can carry a temperature anomaly and exist also on the inhomogeneous temperature background. The modon configurations are then used to initialize numerical simulations, in order to check whether such coherent structures can exist in the full TRSW model. The results show that this is, indeed, the case. The parameter regimes and limitations on the structure of the temperature anomaly inside, in order for the modons to persist, are established. It is also shown that the modons keep their coherence even while evolving on the background of meridionally inhomogeneous temperature fields, or while interacting with sharp temperature fronts. A general scenario of disaggregation of the modons, if they are out of the stability domain, is exhibited and analysed.

We study structure formation in two-dimensional turbulence driven by an external force, interpolating between linear instability forcing and random stirring, subject to nonlinear damping. Using extensive direct numerical simulations, we uncover a rich parameter space featuring four distinct branches of stationary solutions: large-scale vortices, hybrid states with embedded shielded vortices (SVs) of either sign, and two states composed of many similar SVs. Of the latter, the first is a dense vortex gas where all SVs have the same sign and diffuse across the domain. The second is a hexagonal vortex crystal forming from this gas when the linear instability is sufficiently weak. These solutions coexist stably over a wide parameter range. The late-time evolution of the system from small-amplitude initial conditions is nearly self-similar, involving three phases: initial inverse cascade, random nucleation of SVs from turbulence and, once a critical number of vortices is reached, a phase of explosive nucleation of SVs, leading to a statistically stationary state. The vortex gas is continued in the forcing parameter, revealing a sharp transition towards the crystal state as the forcing strength decreases. This transition is analysed in terms of the diffusivity of individual vortices using ideas from statistical physics. The crystal can also decay via an inverse cascade resulting from the breakdown of shielding or insufficient nonlinear damping acting on SVs. Our study highlights the importance of the forcing details in two-dimensional turbulence and reveals the presence of non-trivial SV states in this system, specifically the emergence and melting of a vortex crystal.

We investigate experimentally the effects of micro-grooves on the development of a zero pressure gradient turbulent boundary layer at two different values of the friction Reynolds number. We consider both the well-known streamwise aligned riblets as well as wavy riblets, characterized by a sinusoidal pattern in the mean flow direction. Previous investigations by the authors showed that sinusoidal riblets yield larger values of drag reduction with respect to the streamwise aligned ones. We perform new particle image velocimetry experiments on wall-parallel planes to get insights into the effect of the sinusoidal shape on the near-wall organisation of the boundary layer and the structures responsible for the friction drag reduction and the turbulence generation. Conditional averages, aimed at identifying the topology of the low-speed streaks in the turbulent boundary layer, reveal that the flow is highly susceptible to wall manipulation. This is particularly evident in the cases that are associated with greater values of drag reduction. The results suggest a fragmentation and/or weakening of the streaks in the sinusoidal cases, that is triggered by the larger values of the wall-normal vorticity found at the streaks’ edges. The results are also confirmed by applying the variable interval spatial averaging events eduction technique. The turbulent kinetic energy budget also shows that the sinusoidal geometry significantly attenuates the turbulence production, hence supporting the idea of the manipulation of the turbulence regeneration cycle.

The phenomenon of the ‘coil-stretch’ (C-S) transition, wherein a long-chain polymer initially in a coiled state undergoes a sudden configuration change to become fully stretched under steady elongational flows, has been widely recognized. This transition can display intricate hysteresis behaviours under specific critical conditions, giving rise to unique rheological characteristics in dilute polymer solutions. Historically, microscopic stochastic models and Brownian dynamics simulations have shed light on the underlying mechanisms of the transition by uncovering bistable configurations of polymer chains. Following the initial work by Cerf (J. Chem. Phys., vol. 20, 1952, pp. 395–402), we introduce a continuum model in this study to investigate the C-S transition in a constant uniaxial elongational flow. Our approach involves approximating the unfolding process of the polymer chain as an axisymmetric deformation of an elastic particle. We make the assumption that the particle possesses uniform material properties, which can be represented by a nonlinear, strain-hardening constitutive equation to replicate the finite extensibility of the polymer chain. Subsequently, we analytically solve for the steady-state deformation using a polarization method. By employing this reduced model, we effectively capture the C-S transition and establish its specific correlations with material and geometric properties. The hysteresis phenomena can be comprehended through a force-balance analysis, which involves comparing the externally applied viscous forces with the intrinsic elastic responsive forces. We demonstrate that our model, while simple, unveils rich elastohydrodynamics of the C-S transition.

Soft intertwined channel systems are frequently found in fluid flow networks in nature. The passage geometry of these systems can deform due to fluid flow, which can cause the relationship between flow rate and pressure drop to deviate from the Hagen–Poiseuille linear law. Although fluid–structure interactions in single deformable channels have been extensively studied, such as in Starling's resistor and its variations, the flow transport capacity of an intertwined channel with multiple self-intersections (a ‘hydraulic knot’), is still an open question. We present experiments and theory on soft hydraulic knots formed by interlinked microfluidic devices comprising two intersecting channels separated by a thin elastomeric membrane. Our experiments show flow–pressure relationships similar to flow limitation, where the limiting flow rate depends on the knot configuration. To explain our observations, we develop a mathematical model based on lubrication theory coupled with tension-dominated membrane deflections that compares favourably with our experimental data. Finally, we present two potential hydraulic knot applications for microfluidic flow rectification and attenuation.

Melting icebergs provide nearly half of the total freshwater flux from ice shelves to the ocean, but the availability of accurate, data-constrained melting rate parametrisations limits the correct representation of this process in ocean models. Here, we investigate the melting of a vertical ice face in a warm, salt-stratified environment in a laboratory setting. Observations of the depth-dependent melting rates ${m}$ and boundary layer flow speed $U$ are reported for a range of initially uniform far-field ambient temperatures $T_a$ above ${10}\,^{\circ }{\rm C}$. Ice scallops are characteristic features observed in all experiments, with the width of the scallops consistent with the theory of double-diffusive layers. The morphology of the scallops changes from symmetric about the scallop centre in the colder experiments to asymmetric in the warmer experiments. Observed melting rates are consistent with a melting rate scaling of the form ${m}\propto U\,\Delta T_a$ proposed by previous work in less extreme parameter regimes, where $\Delta T_a$ is the magnitude of thermal driving between the ambient and ice–fluid interface. Our results indicate that ice scalloping is closely linked to the naturally convecting flow of the ambient fluid. Depth-averaged melting rates depend on the buoyancy frequency in the ambient fluid, and double-diffusive convection promotes a turbulent-flux regime distinct from that explained previously in an unstratified regime. These findings have implications for parametrising melting rates of icebergs and glaciers in numerical models or potential freshwater harvesting operations, and provide insights into the interplay between stratification and ice melting.

This study presents the interplay of flow and acoustics within tandem deep cavities, focusing on the resonance mechanism occurring between turbulent shear layers and acoustic eigenmodes. The arrangement inside the tandem deep cavities includes both close and remote configurations. A combined fully coupled and decoupled aeroacoustic simulation strategy was devised. Employing an advanced high-order spectral/hp element method in conjunction with implicit large eddy simulation, the nonlinear compressible Navier–Stokes equations were solved to acquire internal flow–acoustic resonant field. In parallel, the linearized Navier–Stokes equations were tackled to determine coherent shear layer perturbations with external acoustic forcing. Based on acoustic measurements, the mainstream Reynolds number approaches approximately $R{e_{in}} = {O}({10^5})$, where we identified the presence of frequency lock-in and a resonance range. Aeroacoustic noise sources were examined by implementing spectral proper orthogonal decomposition to decompose the pressure fields into hydrodynamic and acoustic components. As feedback intensified, the flow characteristics by the acoustic forcing effect and the flow-interactive effect were categorized according to the development of concurrent turbulent shear layers. Subsequently, the alternating and synchronous behaviours of concurrent shear layers resonated with the out-of-phase and in-phase acoustic eigenmodes were identified, and the corresponding large-scale counter-rotating vortex pairs and co-rotating vortex structures at the cavity entrances were extracted. The acoustic power generated by the Coriolis force was calculated using Howe's vortex-sound analogy, and the aeroacoustic energy transfer mechanism between large-scale shear layer vortices with acoustic eigenmodes was further explored. Finally, a linear response of coherent perturbations of the concurrent shear layers by external acoustic forcing was established. The amplification of flow in the streamwise direction toward the main duct led to the formation of coherent vortex structures, accompanied by separation bubbles into the main duct.

A variational principle is proposed to derive the governing equations for the problem of ocean wave interactions with a floating ice shelf, where the ice shelf is modelled by the full linear equations of elasticity and has an Archimedean draught. The variational principle is used to form a thin-plate approximation for the ice shelf, which includes water–ice coupling at the shelf front and extensional waves in the shelf, in contrast to the benchmark thin-plate approximation for ocean wave interactions with an ice shelf. The thin-plate approximation is combined with a single-mode approximation in the water, where the vertical motion is constrained to the eigenfunction that supports propagating waves. The new terms in the approximation are shown to have a major impact on predictions of ice shelf strains for wave periods in the swell regime.