Low-frequency phenomena in an incompressible pressure-induced laminar separation bubble (LSB) on a flat plate is investigated using direct numerical simulation. The LSB configuration of Spalart and Strelets (J. Fluid Mech., vol. 403, 2000, pp. 329–349) is used. Wall pressure spectra indicate low-frequency-flapping $(St \sim 0.08)$ and high-frequency-shedding $(St \sim 1.52)$ regimes. Conditional velocity averages based on the fraction of reversed flow reveal the low frequency as an expansion/contraction of the LSB. While the high frequency only exhibits exponential growth within the LSB up to breakdown of the spanwise rollers, the low frequency and velocity fluctuations exhibit exponential growth upstream of separation. Instantaneous flow fields reveal large streamwise streaky structures forming within the LSB and extending past reattachment, much like high and low speed streaks in turbulent boundary layers. A predominance of sweep-like events ($Q4$) is observed during contraction and of ejection-like events ($Q2$) during expansion. These motions appear as dominant low-frequency modes in three-dimensional proper orthogonal and dynamic mode decompositions, exhibiting spatial amplification from separation to reattachment. The advection of a group of spanwise alternating streaky structures past the LSB results in an overall contraction after which the bubble expands to its ‘unforced’ state in the absence of the streaks. The low frequency then corresponds to the time it takes for streaks to form, amplify and advect past the LSB from separation to reattachment. This behaviour is linked to the mean flow deformation reported by Marxen and Rist (J. Fluid Mech., vol. 660, 2010, pp. 37–54), where the presence of streaks results in reduced mean bubble size. The formation of these streaky structures, in the absence of free stream turbulence, may be attributed to an absolute instability of the LSB due to the development of a secondary bubble within the primary.

The propagation paths of oceanic internal tides are influenced by their interactions with vortices. We examine the scattering effect that an isolated vortex in (cyclo)geostrophic balance has on a rotating shallow-water plane wave. We run a suite of simulations in which we vary the non-dimensional vorticity of the vortex, $Ro$, the relative scale of the vortex size to the Rossby radius of deformation, $Bu$, and the size of the vortex compared with the plane wave wavelength, $K$. We compare the scattered wave flux pattern with ray-tracing predictions. Ray-tracing predictions are relatively insensitive to $K$ in the $1< K<4$ range we investigate; however, they generally underestimate the broad angles of the shallow-water wave scattering patterns, especially for the lower end of the $K$ range. We then measure the ratio of the scattered wave energy flux to the incoming wave energy flux, denoted by $S$, for each simulation. We find that $S$ follows a power law $S \propto (FrK)^2$ when $S < 0.2$, where $Fr = Ro/\sqrt {Bu}$ is the Froude number. When $S>0.2$, it starts plateauing.

For several applications there are advantages in writing turbulent flow equations in a coordinate frame aligned with the streamlines and several two-dimensional examples of this approach have appeared in the literature. In this paper, we extend this approach to general three-dimensional flows. We find that, in any flow that has a component of its vorticity aligned in the streamline direction, congruences of its streamlines do not form integrable manifolds. This limits the development of a streamline coordinate description of such flows, although some useful results can still be obtained. However, in the case of general three-dimensional complex-lamellar flows, where the mean velocity and mean vorticity are everywhere orthogonal, a complete streamline coordinate description can be derived. Furthermore, we show that general complex-lamellar flows are a good approximation to boundary layers and thin free shear layers. We derive the underlying true coordinate system for such flows, where the orthogonal coordinate surfaces are two stream surfaces and a modified potential surface. From this we obtain physical equations, where flow variables have the same dimensions they would have in a Cartesian coordinate frame. Finally, we show that rational approximations to these equations, which describe small-perturbation flows, contain some terms that have been ignored in previous applications and we detail some practical applications of the theory in modelling and analysis.

The dimensional transition in turbulent jets of a shear-thinning fluid is studied via direct numerical simulations. Our findings reveal that under vertical confinement, the flow exhibits a unique mixed-dimensional (or 2.5-dimensional) state, where large-scale two-dimensional and small-scale three-dimensional structures coexist. This transition from three-dimensional turbulence near the inlet to two-dimensional dynamics downstream is dictated by the level of confinement: weak confinement guarantees turbulence to remain three-dimensional, whereas strong confinement forces the transition to two dimensions; the mixed-dimensional state is observed for moderate confinement and it emerges as soon as flow scales are larger than the vertical length. In this scenario, we observed that the mixed-dimensional state is an overall more energetic state, and it shows a multi-cascade process, where the direct cascade of energy at small scales and the direct cascade of enstrophy at large scales coexist. The results provide insights into the complex dynamics of confined turbulent flows, relevant in both natural and industrial settings.

Vapour bubbles produced by long-pulsed laser often have complex non-spherical shapes that reflect some characteristics of the laser beam. The transition between two commonly observed shapes, namely, a rounded pear-like shape and an elongated conical shape, is studied using a new computational model that combines compressible multiphase fluid dynamics with laser radiation and phase transition. Two laboratory experiments are simulated, in which a holmium:YAG or thulium fibre laser is used to generate bubbles of different shapes. In both cases, the predicted bubble nucleation and morphology agree reasonably well with the experimental observation. The full-field results of laser irradiance, temperature, velocity and pressure are analysed to explain bubble dynamics and energy transmission. It is found that due to the lasting energy input, the vapour bubble's dynamics is driven not only by advection, but also by the continued vaporisation at its surface. Vaporisation lasts less than $1~{\rm \mu}$s in the case of the pear-shaped bubble, compared with over $50~{\rm \mu}$s for the elongated bubble. It is thus hypothesised that the bubble's morphology is determined by competition. When the speed of advection is higher than that of vaporisation, the bubble tends to grow spherically. Otherwise, it elongates along the laser beam direction. To test this hypothesis, the two speeds are defined analytically using a model problem, then estimated for the experiments using simulation results. The results support the hypothesis. They also suggest that when the laser's power is fixed, a higher laser absorption coefficient and a narrower beam facilitate bubble elongation.

We investigate a turbulent stratified plane Poiseuille flow using linear models and nonlinear simulations. We propose the first complete explanation for the prolific and coherent backward (BWs)- and forward-propagating waves (FWs), which have been observed in these flows. We demonstrate a significant presence of oblique waves in the channel core, particularly for the FWs. Critically, we show that neglect of spanwise structure leads to a distorted dispersion relation due to its strong dependence on the angle of obliquity. Interestingly, solutions to the Taylor–Goldstein equations show that wave dynamics is strongly dependent on shear, with only a weak dependence on buoyancy for the BWs at low and order-one wavenumbers, when the wavenumber is scaled by the channel half-height. As the wavenumber increases, waves transition from a shear-dominated regime to a buoyancy-dominated regime, with their dispersion relation tending towards that of idealised internal waves subject to a shear-free and constant-buoyancy-gradient flow, with a characteristic velocity and buoyancy frequency corresponding to respective centreline values in the channel. Finally, we show that the dominance of the BWs arises due to the external forcing of the system, whereby turbulent fluid ejected into the core has a lower momentum when compared with the local flow, therefore preferentially generating BWs in the channel. Qualitatively, channel-core dynamics can be reproduced with low-momentum forcing to a velocity profile with a velocity maximum and a corresponding negative second derivative intersecting a region of strong buoyancy gradient. This structure is inherent to a wide variety of jet-like environmental, atmospheric and industrial flows, suggesting that BWs are a critical control on dynamics of such flows.

We study the effect of gas rarefaction on the interaction of small thermodynamic non-uniformities with a finite body. Considering a two-dimensional set-up, the initial system state is modelled as slight perturbations over its uniform density and temperature fields, prescribed in the vicinity of a thin plate. The problem is analysed in the collisionless limit and complemented by direct simulation Monte Carlo computations to cover the entire range of gas rarefaction rates. The high-Knudsen ‘sink-like’ and ‘source-like’ propagation patterns observed in the density- and temperature-driven set-ups, respectively, are discussed, together with the impact of specular (smooth) and diffuse (isothermal) wall reflections. At highly rarefied conditions, the solid body obstructs part of the gas domain, preventing the propagation of acoustic disturbances therein. With decreasing gas rarefaction, the acoustic field penetrates the obscured area via the effect of molecular collisions. Inspecting the near-field description, the propagation of flow disturbances along the plate surface is examined, and the acoustic force on the body is computed. In the thermally excited case, both normal- and shear-force components change sign at late times, attracting the plate towards the initial perturbation location. With reducing gas rarefaction, the shear force diminishes while the normal force sharply increases due to the decrease in signal decay. Finally, we apply the analysis to study the impact of gas rarefaction on acoustic reciprocity. Notably, acoustic reciprocity does not hold at non-continuum conditions over non-specular surfaces, where boundary reflections propagate in the presence of few molecular collisions, insufficient to retain reciprocal symmetry.