How locally injected turbulence spreads in space is investigated with direct numerical simulations. We consider a turbulent flow in a long triply periodic box generated by a forcing that is localized in space. The forcing is such that it does not inject any mean momentum into the flow. We show that at long times a statistically stationary state is reached where the turbulent energy density in space fluctuates around a mean profile that peaks at the forcing location and decreases fast away from it. We measure this profile as a function of the distance from the forcing region for different values of the Reynolds number. It is shown that, as the Reynolds number is increased, it converges to a Reynolds-number-independent profile, implying that turbulence spreads due to self-advection and not due to molecular diffusion. In this limit, therefore, turbulence plays the simultaneous role of cascading the energy to smaller scales and transporting it to larger distances. The two effects are shown to be of the same order of magnitude. Thus a new turbulent state is reached where turbulent transport and turbulent cascade are equally important and control its properties.

We derive expressions relating the entrainment fluxes of momentum and kinetic energy, relative to the mass flux entrained into a turbulent wake exposed to a turbulent background. These expressions contain correlations between the entrainment velocity and the turbulent fluctuations within the background. We perform high-resolution, simultaneous particle image velocimetry and planar laser-induced fluorescence experiments, and observe these correlations to be negligible in the far wake, such that momentum and kinetic energy are entrained into the wake with the same relative efficiency to mass as from an idealised, non-turbulent background. This is a useful result in the context of modelling, since the entrainment hypothesis (Turner, J. Fluid Mech., vol. 173, 1986, pp. 431–471) can still be used to model the entrainment of momentum and kinetic energy. Nevertheless, the entrainment rate of mass is shown to vary spatially, and with the specific nature of the background turbulence, so this in turn drives a spatial/background-turbulence-specific entrainment rate of momentum/kinetic energy. Contrastingly, in the near wake, whilst momentum is entrained from a turbulent background with the same relative efficiency to mass as for an idealised non-turbulent background, this is not the case for kinetic energy. Owing to the sum of multiple positive, small-valued correlations between the fluctuations in the background and the entrainment velocity, kinetic energy is entrained more efficiently than in the idealised case. This includes entrainment from a non-turbulent background, where small correlations are observed between the irrotational background fluctuations and the entrainment velocity. Evidence is also presented that the entrainment velocity scales with the Kolmogorov velocity scale when the background is turbulent.

The present paper concerns the linear fate of transverse perturbations in a gravity-driven, thin-film flow over a soluble substrate. We propose a reduced-order model, based on a boundary-layer treatment of the solute transport and a depth-integration of the Stokes equations, using two extended lubrication methodologies found in the literature. We obtain a closed-form dispersion relation, which we compare to a previous, fully resolved analytical investigation (Bertagni and Camporeale, J. Fluid Mech., vol. 913, 2021, A34). The results allow us to distil the essential physical mechanisms behind the instability.

The scaling and mechanism of the propagation speed of turbulent fronts in pipe flow with the Reynolds number has been a long-standing problem in the past decades. Here, we derive an explicit scaling law for the upstream front speed, which approaches a power-law scaling at high Reynolds numbers, and we explain the underlying mechanism. Our data show that the average wall distance of low-speed streaks at the tip of the upstream front, where transition occurs, appears to be constant in local wall units in the wide bulk-Reynolds-number range investigated, between 5000 and 60 000. By further assuming that the axial propagation of velocity fluctuations at the front tip, resulting from streak instabilities, is dominated by the advection of the local mean flow, the front speed can be derived as an explicit function of the Reynolds number. The derived formula agrees well with the speed measured by front tracking. Our finding reveals a relationship between the structure and speed of a front, which enables a close approximation to be obtained of the front speed based on a single velocity field without having to track the front over time.

We introduce a geometric analysis of turbulent mixing in density-stratified flows based on the alignment of the density gradient in two orthogonal bases that are locally constructed from the velocity gradient tensor. The first basis connects diapycnal mixing to rotation and shearing motions, building on the recent ‘rortex–shear decomposition’ in stratified shear layers (Jiang et al., J. Fluid Mech., vol. 947, 2022, A30), while the second basis connects mixing to the principal axes of the viscous dissipation tensor. Applying this framework to datasets taken in the stratified inclined duct laboratory experiment reveals that density gradients in locations of high shear tend to align preferentially (i) along the direction of minimum dissipation and (ii) normal to the plane spanned by the rortex and shear vectors. The analysis of the local alignment across increasingly turbulent flows offers new insights into the intricate relationship between the density gradient and dissipation, and thus diapycnal mixing.

We study experimentally the aeroelastic instability boundaries and three-dimensional vortex dynamics of pitching swept wings, with the sweep angle ranging from 0$^\circ$ to 25$^\circ$. The structural dynamics of the wings are simulated using a cyber-physical control system. With a constant flow speed, a prescribed high inertia and a small structural damping, we show that the system undergoes a subcritical Hopf bifurcation to large-amplitude limit-cycle oscillations (LCOs) for all the sweep angles. The onset of LCOs depends largely on the static characteristics of the wing. The saddle-node point is found to change non-monotonically with the sweep angle, which we attribute to the non-monotonic power transfer between the ambient fluid and the elastic mount. An optimal sweep angle is observed to enhance the power extraction performance and thus promote LCOs and destabilize the aeroelastic system. The frequency response of the system reveals a structural-hydrodynamic oscillation mode for wings with relatively high sweep angles. Force, moment and three-dimensional flow structures measured using multi-layer stereoscopic particle image velocimetry are analysed to explain the differences in power extraction for different swept wings. Finally, we employ a physics-based force and moment partitioning method to correlate quantitatively the three-dimensional vortex dynamics with the resultant unsteady aerodynamic moment.

The stability and postcritical behaviour of a horizontal flag undergoing gravity-induced deformation and periodic contact with a nearby horizontal rigid wall are experimentally investigated. The results elucidate the combined effects of gravity and contact on flutter, and reveal design principles for application to triboelectric energy harvesting. By varying the free-stream velocity, flag thickness and distance between the flagpole and the wall, the dynamics of the flag are classified into quasistatic equilibrium, flutter, partial contact and saturated contact modes. Considering the significance of gravitational effects, a new dimensionless flow velocity is proposed to identify the distribution of the dynamic modes, and its definition varies according to whether the wall is placed above or below the flag. The critical conditions for transitions between the dynamic modes are determined from the balance of fluid dynamic and gravitational effects. The distance from the flagpole to the wall is found to be more critical for transitions in the lower-wall configuration than in the upper-wall configuration. The peak contact force as well as the oscillation amplitude and frequency at postequilibrium exhibits remarkably different trends depending on the location of the wall. The peak contact force imposed on the wall by the fluttering flag weakens as the distance to the wall increases in the case of an upper wall, whereas it becomes stronger in the case of a lower wall.

The turbulent external flow around a three-dimensional stepped cylinder is studied by means of direct numerical simulations with the adaptive mesh refinement technique. We give a broad perspective of the flow regimes from laminar to turbulent wake at $Re_D=5000$, which is the highest ever considered for this flow case. In particular, we focus on the intermediate Reynolds number $Re_D=1000$ that reveals a turbulent wake coupled with a stable cylinder shear layer (subcritical regime). This flow shows a junction dynamics similar to the laminar $Re_D=150$, where no hairpin vortex appears around the edges, and just two horseshoe vortices are visible. A new stable vortex in the form of a ring, which coils around the rear area, is also identified. In the turbulent wake, the presence of three wake cells is pointed out: the large and small cylinder cells together with the modulation region. However, the modulation dynamics varies between the subcritical and turbulent regimes. A time-averaged, three-dimensional set of statistics is computed, and spatially coherent structures are extracted via proper orthogonal decomposition (POD). The POD identifies the (long-debated) connection between the N-cell and the downwash behind the junction. Furthermore, as the Reynolds number increases, the downwash phenomenon becomes less prominent. Eventually, a reduced-order reconstruction with the most energetically relevant modes is defined to explain the wake vortex interactions. This also serves as a valuable starting point for simulating the stepped cylinder wake behaviour within complex frameworks, e.g. fluid–structure interaction.

Linear stability theory (LST) is often used to model the large-scale flow structures in the turbulent mixing region and near pressure field of high-speed jets. For perfectly expanded single round jets, these models predict the dominance of azimuthal wavenumbers $m=0$ and $m = 1$ helical modes for the lower frequency range, in agreement with empirical data. When LST is applied to twin-jet systems, four solution families appear following the odd/even behaviour of the pressure field about the symmetry planes. The interaction between the unsteady pressure fields of the two jets also results in their coupling. The individual modes of the different solution families no longer correspond to helical motions, but to flapping oscillations of the jet plumes. In the limit of large jet separations, when the jet coupling vanishes, the eigenvalues corresponding to the $m=1$ mode in each family are identical, and a linear combination of them recovers the helical motion. Conversely, as the jet separation decreases, the eigenvalues for the $m=1$ modes of each family diverge, thus favouring a particular flapping oscillation over the others and preventing the appearance of helical motions. The dominant mode of oscillation for a given jet Mach number $M_j$ and temperature ratio $T_R$ depends on the Strouhal number $St$ and jet separation $s$. Increasing both $M_j$ and $T_R$ independently is found to augment the jet coupling and modify the $(St,s)$ map of the preferred oscillation mode. Present results predict the preference of two modes when the jet interaction is relevant, namely varicose and especially sinuous flapping oscillations on the nozzles’ plane.

Solute–surface interactions have garnered considerable interest in recent years as a novel control mechanism for driving unique fluid dynamics and particle transport with potential applications in fields such as biomedicine, the development of microfluidic devices and enhanced oil recovery. In this study, we will discuss dispersion induced by the diffusioosmotic motion near a charged wall in the presence of a solute concentration gradient. Here, we introduce a plug of salt with a Gaussian distribution at the centre of a channel with no background flow. As the solute diffuses, the concentration gradient drives a diffusioosmotic slip flow at the walls, which results in a recirculating flow in the channel; this, in turn, drives an advective flux of the solute concentration. This effect leads to cross-stream diffusion of the solute, altering the effective diffusivity of the solute as it diffuses along the channel. We derive theoretical predictions for the solute dynamics using a multiple-time-scale analysis to quantify the dispersion driven by the solute–surface interactions. Furthermore, we derive a cross-sectionally averaged concentration equation with an effective diffusivity analogous to that from Taylor dispersion. In addition, we use numerical simulations to validate our theoretical predictions.

Existing experimental results show that swirling flames in annular combustors respond with a different gain to acoustic azimuthal modes rotating in either the clockwise or anti-clockwise direction. The ratio $R$ of these two gains is introduced, with $R=1$ being the conventional case of flames responding the same to the two forcing directions. To allow a difference in response to the different directions ($R\neq 1$), a multiple-input single-output azimuthal flame describing function is successfully implemented in a quaternion valued low-order model of an annular combustion chamber in the current work. Theoretical studies have explored this kind of symmetry breaking between the two acoustic wave directions in the past, but it has not been backed by experimental data. One of the main features of the new model proposed in this work is the potential difference in mode shapes between the acoustic and the heat release rate modes, which has recently been observed experimentally. This results in a gain-dependent equation for the nature of the mode, which has a significant influence on the fixed points of the system. For example, one of the spinning solutions and the standing solution can disappear through a saddle node bifurcation as the parameters are varied. The presence of only a single direction for the spinning solution matches experimental observations better than the conventional models, and the proposed model is shown to qualitatively describe experimental measurements well.

Invariant solutions of the Navier–Stokes equations play an important role in the spatiotemporally chaotic dynamics of turbulent shear flows. Despite the significance of these solutions, their identification remains a computational challenge, rendering many solutions inaccessible and thus hindering progress towards a dynamical description of turbulence in terms of invariant solutions. We compute equilibria of three-dimensional wall-bounded shear flows using an adjoint-based matrix-free variational approach. To address the challenge of computing pressure in the presence of solid walls, we develop a formulation that circumvents the explicit construction of pressure and instead employs the influence matrix method. Together with a data-driven convergence acceleration technique based on dynamic mode decomposition, this yields a practically feasible alternative to state-of-the-art Newton methods for converging equilibrium solutions. We compute multiple equilibria of plane Couette flow starting from inaccurate guesses extracted from a turbulent time series. The variational method outperforms Newton(-hookstep) iterations in converging successfully from poor initial guesses, suggesting a larger convergence radius.

Considerable effort has been directed towards the characterization of chiral mesoscale structures, as shown in chiral protein assemblies and carbon nanotubes. Here, we establish a thermally driven hydrodynamic description for the actuation and separation of mesoscale chiral structures in a fluid medium. Cross-flow of a Newtonian liquid with a thermal gradient gives rise to an effective torque that propels each particle of a chiral suspension according to its handedness. In turn, the chiral suspension alters the liquid flow, which thus acquires a transverse (chiral) velocity component. Since observation of the predicted effects requires a low degree of sophistication, our work provides an efficient and inexpensive approach to test and calibrate chiral particle propulsion and separation strategies.

Recent studies reveal the dramatic impact of seafloor roughness on the dynamics and stability of broad oceanic flows. These findings motivate the development of parameterizations that concisely represent the effects of small-scale bathymetric patterns in theoretical and coarse-resolution numerical circulation models. The previously reported quasi-geostrophic ‘sandpaper’ theory of flow–topography interaction a priori assumes gentle topographic slopes and weak flows with low Rossby numbers. Since such conditions are often violated in the ocean, we now proceed to formulate a more general model based on shallow-water equations. The new version of the sandpaper model is validated by comparing roughness-resolving and parametric simulations of the flow over a corrugated seamount.

Compound bubbles with a liquid coating in another continuous immiscible bulk phase are ubiquitous in a wide range of natural and industrial processes. Their formation, rise and ultimate bursting at the air–liquid interface play crucial roles in the transport and fate of natural organic matter and contaminants. However, the dynamics of compound bubbles has not received considerable attention until recently. Here, inspired by our previous work (Yang et al., Nat. Phys., vol. 19, 2023, pp. 884–890), we investigate the entrainment of daughter oil droplets in bulk water produced by a bursting oil-coated bubble. We document that the size of the entrained daughter oil droplet is affected by the oil coating fraction and the bulk liquid properties. We rationalize this observation by balancing the viscous force exerted by the extensional flow produced by bubble bursting with the capillary force resisting the deformation of the oil coating, and considering the subsequent end-pinching process which finally entrains the daughter oil droplets. We propose a scaling analysis for the daughter oil droplet size that well captures the experimental results for a wide range of oil coating fractions and Ohnesorge numbers of the bulk liquid. In addition, we discuss the non-monotonic variation of daughter droplet size with the Ohnesorge number, and show the eventual absence of daughter droplets because of the strong viscous effect in the high-Ohnesorge-number regime. Our findings may advance the fundamental understanding of compound bubble bursting and provide guidance and modelling constraints for bubble-mediated contaminant transport in liquids.

Recently, a non-local eddy diffusivity model for the turbulent scalar flux was proposed to improve the local model and was validated using direct numerical simulation (DNS) of homogeneous isotropic turbulence with an inhomogeneous mean scalar (Hamba, J. Fluid Mech., vol. 950, 2022, A38). The non-local eddy diffusivity was assumed to be proportional to the two-point velocity correlation that was expressed in terms of the energy spectrum. Because the Fourier transform of velocity in the homogeneous directions was used to define the energy spectrum, it is not yet understood whether the proposed model can be applied to inhomogeneous turbulence. Thus, this study aimed to improve the non-local model using the scale-space energy density instead of the energy spectrum. First, the scale-space energy density based on filtered velocities was examined using the DNS data of homogeneous isotropic turbulence to obtain its simple form corresponding to the Kolmogorov energy spectrum. Subsequently, the two-point velocity correlation was expressed in terms of the scale-space energy density. Using these expressions, a new non-local eddy diffusivity model was proposed and validated using the DNS data. The one-dimensional non-local eddy diffusivity obtained from the new model agrees with the DNS value. The temporal behaviour of the three-dimensional non-local eddy diffusivity was improved compared with the previous model. Because the scale-space energy density was already examined in turbulent channel flow, it is expected that the new non-local model can also be applied to inhomogeneous turbulence and is useful for gaining insight into turbulent scalar transport.

Kinks commonly appear on the convergent shock surface when an internal conical flow deviates from the axisymmetric state. In this paper, the formation mechanisms of kinks on internal conical shocks (ICSs) generated by elliptical ring wedges with typical entry aspect ratios ($AR{\rm s}$) in a Mach 6 flow are revealed using a theoretical method, in which the spatial evolution of the three-dimensional elliptical ICS is converted into a temporal evolution of a two-dimensional elliptical moving shock (EMS) using the hypersonic equivalence principle. To simultaneously track the shock front of the EMS and the disturbances propagating along it, a front-disturbance tracking method (FDTM) based on geometrical shock dynamics is proposed. It is found that the shock–compression disturbances from the same family initially near the major axis catch up with the disturbance initially emitted from the major axis to form kinks on the EMS. The equivalent kink formation positions predicted by the FDTM always lag behind the real kink formation positions on the elliptical ICS because the applicability of the hypersonic equivalence principle decays as the shock strengthens along the incoming flow direction. The accuracy of the equivalent kink formation positions predicted by the FDTM gradually declines with the reduction in $AR$, but it can be significantly improved for all $AR{\rm s}$ after a modification of the equivalent relationship using the shock angle in the major plane of the elliptical ICS, which provides a new way to solve the kinks on the elliptical ICS.

Local instantaneous exchanges of volume, momentum and buoyancy across turbulent/non-turbulent interfaces (TNTIs) and turbulent/turbulent interfaces (TTIs) are studied using data from direct numerical simulations of a turbulent forced fountain. We apply a novel algorithm that enables independent calculation of the instantaneous local entrainment and detrainment fluxes, and therefore, for the first time, the entrainment and detrainment coefficients according to the fountain model (Bloomfield & Kerr, J. Fluid Mech., vol. 424, 2000, pp. 197–216) are determined explicitly. Across the interface between the fountain and the ambient fluid, which is a TNTI, only volume entrainment occurs, and it is well predicted by the fountain model. Across the interface between the rising upflow and falling downflow within the fountain, which is a TTI, both entrainment and detrainment of volume, momentum and buoyancy occur – with the magnitude of both entrainment and detrainment typically being large compared with the net for all exchanges. However, the model seems to be unable to capture the momentum exchanges due to its ignorance of the pressure. We find that each conditional entrainment and detrainment rate, of volume, momentum and buoyancy, can be described accurately by Gaussian profiles, while the net exchange that is the superposition of the entrainment and detrainment cannot. Moreover, the entrainment exchange rate has its maximum closer to the fountain centreline than that of detrainment, explaining the tendency for net entrainment closer to the fountain centreline and net detrainment further away.

The breakup of coaxial cone-jet interfaces to compound droplets in axisymmetric co-flow focusing (CFF) upon actuation is studied through numerical simulations. Due to the coupling effect of double interfaces, the response behaviours of coaxial cone-jet flow to actuation are more complex than those of a single-layered interface structure. Particularly, the coaxial jet presents totally different response modes between weak and strong interface coupling situations. In this work, the phase diagrams of response modes for coaxial jet breakup are depicted, considering the effect of perturbation frequency, amplitude and liquid flow rates. In particular, the breakup of a coaxial jet can be synchronized with actuation within a frequency range containing the natural breakup frequency, resulting in uniform compound droplets with a single core inside the shell, and the size of droplets can be adjusted by frequency. As the perturbation frequency exceeds the upper critical value, the external perturbation is unable to dominate the jet breakup, while below the lower critical frequency, the jet breaks up with multiple droplets generated in one period. The perturbation amplitude mainly affects the jet breakup length and also leads to the transition between different response modes. The coaxial cone upstream of the orifice can act as a buffer layer, regulating the perturbation amplitude of the coaxial jet downstream. The degree of buffering effect is affected by the perturbation frequency and amplitude. As the perturbation amplitude approaches unity, the decrease of perturbation frequency leads to the intermittent jet behaviour from the cone tip with a vibrating manner of the coaxial cone. Based on the linear instability analysis on the simplified single jet models for weak-coupled and strong-coupled jets, scaling analyses are carried out, which predict the jet breakup length and the natural frequency and critical frequency for the synchronized breakup. Finally, a strong pulse is added on the perturbation to produce compound droplets with a controllable number of cores. The present work provides valuable guidance for the practical application of on-demand compound droplet generation through active CFF.

The effect of a uniform mean scalar gradient on the small scales of a passive scalar field in statistically stationary homogeneous isotropic turbulence is investigated through the transport equation for the scalar fluctuation. After some manipulation of the equation, it is shown that the effect can be recast in the form $S_\theta ^* {{Pe^{-1}_{\lambda _\theta }}}$ ($S_\theta ^*$ is the non-dimensional scalar gradient, ${{Pe_{\lambda _\theta }}}$ is the turbulent Péclet number). This effect gradually disappears as ${{Pe_{\lambda _\theta }}}$ becomes sufficiently large, implying a gradual approach towards local isotropy of the passive scalar. It is further argued that, for a given $S_\theta ^*$, the normalized odd moments of the scalar derivative tend towards isotropy as ${{Pe^{-1}_{\lambda _\theta }}}$. This is supported by direct numerical simulations data for the normalized odd moments of the scalar derivative at large Péclet numbers. Further, the present derivation leads to the same prediction (${\sim }Sc^{-0.45}$ where Sc is the Schmidt number) as Buaria et al. (Phys. Rev. Lett., vol. 126, no. 3, 2021a, p. 034504) and complements the derivation by the latter authors, which is based on dimensional arguments and the introduction of a new diffusive length scale.