Direct numerical simulations of temporally developing mixing layers have been performed to investigate the effects of compressibility on statistics and structures near the interfaces of high- and low-speed large-scale structures (LSSs), covering a range of convective Mach numbers from $M_c=0.2$ to $1.8$ and Taylor Reynolds numbers up to 290. The interfaces of LSSs are directly defined by the isosurface of zero fluctuating streamwise velocity. The characteristic velocity jump at the interfaces grows rapidly in the transition stage and then decreases until reaching a gradual self-similar stage. The evolution of velocity jump is evidently delayed as the convective Mach number increases. The interface layer is formed by the shearing motion of neighbouring LSSs. A small vortical motion characterized by the Kolmogorov scale is induced in the interface layer by shear-dominated outer regions. The strengths of outer shearing motion and central vortical motion are greater at the leading edge, but smaller at the trailing edge, which is also reflected in a larger velocity jump at the leading edge. As the convective Mach number increases, the small-scale vortical structure is obviously suppressed by compressibility. At high convective Mach number $M_c=1.8$, the compressive shear-dominated outer regions are linked with a sheet-like structure passing through the centre of the expansion region corresponding to a small-scale vortical structure. The compressibility and shearing strength near the interface are highly dependent on the interface orientation.

The tangential strain rate in premixed flames impacts significantly the flame surface area generation and thus the combustion process. Studies on incompressible isotropic turbulence have revealed that the mean tangential strain rate at material and iso-scalar surfaces is positive and exhibits a universal value when normalized by the Kolmogorov time. This is associated with the preferential alignment of the surface normal with the most compressive principal strain rate. The present study investigates such effects in premixed hydrogen and iso-octane flame kernels using direct numerical simulations. It is shown that the normalized mean tangential strain rate of the investigated flames has a very similar value compared with the incompressible flows. However, in the reaction zone, the flame surface normal aligns preferentially with the most extensive principal strain rate. Furthermore, this alignment depends on the reaction progress variable and the Lewis number, while the tangential strain rate remains independent of these parameters. Such counter-intuitive behaviour is systematically investigated by decomposing the effects of dilatation and residual solenoidal turbulence. It is found that the solenoidal turbulence influences significantly the tangential strain rate. A general effect of turbulence on the tangential strain rate is identified, which is consistent with incompressible flows and independent of the Lewis number and the reaction progress variable. This is a remarkable finding indicating that models of the tangential strain rate developed based on incompressible flows apply also to premixed flames with different Lewis numbers, and, for the modelling, only the solenoidal turbulence should be considered.

Splashing of impacting drops produces a myriad of secondary spray droplets, which generate aerosols during rain on the ocean and can cause health hazards during the spraying of pesticides or enhance the droplet transmission of disease. Determining the size and number of the finest splashed droplets is therefore of practical interest. Herein, we use a novel experimental facility with a 26 m tall vacuum tube, to study well-controlled drop impacts at velocities as high as 22 m s$^{-1}$, where we reach parameter regimes not studied before using freely falling drops. Using extreme video frame rates, we pinpoint the primary source of the finest spray, coming from the catastrophic bending and rupture of the sub-micron-thick ejecta sheet, which emerges at a high speed from the neck connecting the drop and pool. The axisymmetric bending and convoluted ejecta shapes are driven primarily by resistance from the surrounding air, but also depend on the viscosity difference between drop and pool, which influences the initial ejection angle of the sheet. These extreme impact conditions provide new insights into general spray formation, through a sequence of bucklings of the rising ejecta, which dances next to the drop surface and can also form an enclosed air torus.

Industrial applications of flow through fractures such as geothermal energy or hydraulic stimulation involve forcing large flow rates through small fractures, thereby inducing inertial fluid behaviours and turbulence. The most common fracture flow model, Poiseuille flow (the cubic law), is incapable of capturing these phenomena and thus the impact of inertial and turbulent forces in fracture flow has remained relatively unexplored. The GG22 flow model is a newly derived fracture flow model that is capable of capturing inertial, transient and turbulent forces. In this article, we apply the GG22 flow model to hydraulic stimulation of radial fractures for the first time to determine how these phenomena manifest. We show that inertia and turbulence only manifest near the wellbore (within 30 radii) and lead to changes in fracture shape and injection pressure but have little impact on tip behaviour. Turbulence increases wellbore pressure and aperture while inertia decreases wellbore pressure and aperture. The majority of the pressure loss along the fracture occurs near the wellbore and is captured by turbulence where entrance correction factors would otherwise be needed. Using water, turbulence is the dominant mechanism that causes departures from Poiseuille flow at high $Re$. The solution departs immediately upon the manifestation of turbulence ($Re\geq 2\times 10^3$), while inertial effects manifest at higher flow rates ($Re\geq 2\times 10^4$). Using slickwater, the opposite trend is observed: inertial effects manifest first ($Re\geq 5\times 10^3$), while turbulent effects are delayed ($Re\geq 10^4$). In both cases, the threshold for departures from the Poiseuille flow solution are low and the differences are large.

In this paper we discuss the dynamics of vorticity at partial-slip boundaries. We consider the total vector circulation, which includes both the total vorticity of the fluid and the slip velocity at the boundary (the interface vortex sheet). The generation of vector circulation is an inviscid process, which does not depend on either viscosity or the slip length at the boundary. Vector circulation is generated by the inviscid relative acceleration between the fluid and the solid, due to either tangential pressure gradients or tangential acceleration of the partial-slip wall. While the slip length does not affect the creation of vector circulation, it governs how vector circulation is distributed between the total vorticity of the fluid and the interface vortex sheet. Specifically, the partial-slip boundary condition prescribes the ratio between boundary vorticity and the strength of the interface vortex sheet, and the viscous boundary flux transfers vector circulation between the interface vortex sheet and the fluid interior to maintain this condition. The interaction between a vortex ring and a partial-slip wall is examined to highlight various aspects of this formulation. For the head-on collision, the quantity of vector circulation diffused into the fluid as secondary vorticity increases as the slip length is decreased, resulting in a stronger secondary vortex and increased rebound of the vortex ring. For the oblique interaction, the extent to which the vortex ring connects to the boundary increases as the slip length is increased.

Direct numerical simulations of turbulent flows over highly permeable porous walls were performed at various Reynolds numbers to examine the effects of the Reynolds number on permeable wall turbulence. The porous medium consisted of Kelvin cell arrays with porosity $0.95$, and the permeability Reynolds number $Re_K$ ranged from approximately 7 to 50. Simulations with thin and thick porous walls were performed to investigate the effects of spanwise roller vortices associated with the Kelvin–Helmholtz instability. The results show that the effect of the Kelvin–Helmholtz instability becomes more significant with increasing the permeability Reynolds number, and spanwise rollers, for which length scale is an order of channel height, dominate turbulence when $Re_K \gtrsim 30$. Spanwise rollers reinforce the negative correlation between the wall-normal and streamwise velocity fluctuations close to the porous/fluid interface, and intensify the turbulent velocity fluctuations away from the porous walls, leading to increased frictional resistance. An investigation of the Reynolds number dependence of the modified logarithmic law indicates that the zero-plane displacement and equivalent roughness height are proportional to the square root of permeability, whereas the von Kármán constant increases with the permeability Reynolds number because of the increased mixing length resulting from the relatively large-scale velocity fluctuations induced by spanwise rollers. We developed a model for the modified log law for permeable wall turbulence based on permeability, and confirmed that the skin friction coefficient obtained from the model reasonably predicts the skin friction coefficient for several types of high-porosity porous media. Hence, permeability is a key parameter that characterizes the logarithmic mean velocity profiles over a variety of porous media with high porosity.

Beneitez et al. (Phys. Rev. Fluids, vol. 8, 2023, L101901) have recently discovered a new linear ‘polymer diffusive instability’ (PDI) in inertialess rectilinear viscoelastic shear flow using the finitely extensible nonlinear elastic constitutive model of Peterlin (FENE-P) when polymer stress diffusion is present. Here, we examine the impact of inertia on the PDI for both plane Couette and plane Poiseuille flows under varying Weissenberg number ${W}$, polymer stress diffusivity $\varepsilon$, solvent-to-total viscosity ratio $\beta$ and Reynolds number ${Re}$, considering the FENE-P and simpler Oldroyd-B constitutive relations. Both the prevalence of the instability in parameter space and the associated growth rates are found to significantly increase with ${Re}$. For instance, as $Re$ increases with $\beta$ fixed, the instability emerges at progressively lower values of $W$ and $\varepsilon$ than in the inertialess limit, and the associated growth rates increase linearly with $Re$ when all other parameters are fixed. For finite $Re$, it is also demonstrated that the Schmidt number $Sc=1/(\varepsilon Re)$ collapses curves of neutral stability obtained across various $Re$ and $\varepsilon$. The observed strengthening of PDI with inertia and the fact that stress diffusion is always present in time-stepping algorithms, either implicitly as part of the scheme or explicitly as a stabilizer, implies that the instability is likely operative in computational work using the popular Oldroyd-B and FENE-P constitutive models. The fundamental question now is whether PDI is physical and observable in experiments, or is instead an artifact of the constitutive models that must be suppressed.

This study extends our previous work (McCloughan & Suslov, J. Fluid Mech., vol. 887, 2020, A23), where the existence of a saddle-node bifurcation of steady axisymmetric electrolyte flows driven by the Lorentz force in a shallow annular domain was first reported. Here we perform further weakly nonlinear analysis over a wider range of the governing parameters to demonstrate that the previously reported saddle-node bifurcation is a local feature of a global fold catastrophe, which, in turn, is a part of cusp catastrophe occurring as the thickness of the fluid layer increases. The amplitude equation characterising multiple flow solutions in the finite vicinity of catastrophe points is derived. The sensitivity of its coefficients and solutions to the distance from the catastrophe points is assessed demonstrating the robustness of the used analytical procedure. The asymptotic flow solution past the catastrophe point is subsequently obtained and its topology is explored confirming the existence of the secondary circulation in the bulk of flow (two-tori background flow structure). The latter is argued to lead to the appearance of experimentally observable vortices on the fluid surface. The rigorous justification of this conjecture is to be given in Part 2 of the study.

The interaction between ingested turbulence and rotating blades is a key source of broadband noise in engineering applications. In this paper, a far-field noise model accounting for source correlation across the span of the blade and between blades is developed and applied to the study of homogeneous isotropic turbulence ingestion by a model cooling fan, wind turbine and aircraft propeller. The widely used theory of Amiet is revisited and it is shown that previous works produce conflicting results when attempting to account for blade-to-blade correlation. Central to the model is the calculation of the time between blade chops of the same turbulent eddy as heard by the observer. In this paper it is shown that, when derived correctly, Amiet's theory accounts for correlated sources between blades and, thus, can predict haystacking tones. Comparisons with the new rotational formulation and with experimental data enable us to show that Amiet's theory can be used to accurately predict turbulence ingestion noise from open rotors. In particular, it is found that the infinite-span assumption in strip theory and the neglect of correlation effects across the span do not undermine the accuracy of this theory. This is of great importance because, unlike Amiet's theory, models which treat rotational effects and source correlation exactly are expensive to evaluate routinely at high frequencies due to the slow convergence of infinite series with Bessel functions.

Cross-flow transition over a delta wing is systematically studied in a Mach 6.5 hypersonic wind tunnel, employing the Rayleigh scattering flow visualisation, high-speed schlieren and fast-response pressure sensors. Direct numerical simulations and analysis based on linear stability theory under the same flow conditions are applied to analyse the transition mechanism. Three unstable modes are identified: the travelling cross-flow instabilities, the second mode and the low-frequency waves. It is shown that the travelling cross-flow vortices first appear in the cross-flow region near the leading edge of the model. These vortices can modulate the mean profile of the flow, which benefits the growth of second mode. A phase-locked interaction mechanism transfers energy from the cross-flow instabilities to the high-frequency second mode, leading to amplification at the expense of the cross-flow instability. As the second mode grows to a critical amplitude, it triggers a $z$-type secondary instability within a similar frequency range, which introduces secondary finger-like structures connecting to the cross-flow vortex. It is further found that the generation of these finger-like structures is related to the expansion and compression of the second mode. These finger vortices further evolve along the streamwise direction into low-frequency waves and corresponding hairpin-like structures that finally trigger turbulence. An interaction mechanism likely exists between the secondary instability and the low-frequency waves, since their phase speeds are approaching each other. These observations of the interaction mechanism are consistent with those of previous studies on hypersonic boundary layers (Zhang et al., Phys. Fluids, vol. 32 (7), 2020, 071702; Li et al., Phys. Fluids, vol. 32 (5), 2020, 051701).

Artificial intelligence (AI) has always drawn inspiration from the brain, from its most basic forms like nodes and layers to more recent advances that mimic individual neurons and various aspects of visual and sensory processing.