Relative dispersion experiments are often analysed using theoretical predictions from two- and three-dimensional turbulence. These apply to infinite inertial ranges, assuming the same dispersive behaviour over all scales. With finite inertial ranges, the metrics are less conclusive. We examine this using pair separation probability density functions (PDFs), obtained by integrating a Fokker–Planck equation with different diffusivity profiles. We consider time-based metrics, such as the relative dispersion, and separation-based metrics, such as the finite scale Lyapunov exponent (FSLE). As the latter cannot be calculated from a PDF, we introduce a new measure, the cumulative inverse separation time (CIST), which can. This behaves like the FSLE, but advantageously has analytical solutions in the inertial ranges. This allows the establishment of consistency between the time- and space-based metrics, something which has been lacking previously. We focus on three dispersion regimes: non-local spreading (as in a two-dimensional enstrophy inertial range), Richardson dispersion (as in an energy inertial range) and diffusion (for uncorrelated pair motion). The time-based metrics are more successful with non-local dispersion, as the corresponding PDF applies from the initial time. Richardson dispersion is barely observed, because the self-similar PDF applies only asymptotically in time. In contrast, the separation-based CIST correctly captures the dependencies, even with a short (one decade) inertial range, and is superior to the traditional FSLE at large scales. Nevertheless, it is advantageous to use all measures together, to seek consistent indications of the dispersion.

The well-known energy dissipation anomaly in the inviscid limit, related to velocity singularities according to Onsager, still needs to be demonstrated by numerical experiments. The present work contributes to this topic through high-resolution numerical simulations of the inviscid three-dimensional Taylor–Green vortex problem using a novel high-order discontinuous Galerkin discretisation approach for the incompressible Euler equations. The main methodological ingredient is the use of a discretisation scheme with inbuilt dissipation mechanisms, as opposed to discretely energy-conserving schemes, which – by construction – rule out the occurrence of anomalous dissipation. We investigate effective spatial resolution up to $8192^3$ (defined based on the $2{\rm \pi}$-periodic box) and make the interesting phenomenological observation that the kinetic energy evolution does not tend towards exact energy conservation for increasing spatial resolution of the numerical scheme, but that the sequence of discrete solutions seemingly converges to a solution with non-zero kinetic energy dissipation rate. Taking the fine-resolution simulation as a reference, we measure grid-convergence with a relative $L^2$-error of $0.27\,\%$ for the temporal evolution of the kinetic energy and $3.52\,\%$ for the kinetic energy dissipation rate against the dissipative fine-resolution simulation. The present work raises the question of whether such results can be seen as a numerical confirmation of the famous energy dissipation anomaly. Due to the relation between anomalous energy dissipation and the occurrence of singularities for the incompressible Euler equations according to Onsager's conjecture, we elaborate on an indirect approach for the identification of finite-time singularities that relies on energy arguments.

The behaviour of turbulent, buoyant, planar plumes is fundamentally coupled to the environment within which they develop. The effect of a background stratification directly influences a plumes buoyancy and has been the subject of numerous studies. Conversely, the effect of an ambient co-flow, which directly influences the vertical momentum of a plume, has not previously been the subject of theoretical investigation. The governing conservation equations for the case of a uniform co-flow are derived and the local dynamical behaviour of the plume is shown to be characterised by the scaled source Richardson number and the relative magnitude of the co-flow and plume source velocities. For forced, pure and lazy plume release conditions the co-flow acts to narrow the plume and reduce both the dilution and the asymptotic Richardson number relative to the classic zero co-flow case. Analytical solutions are developed for pure plumes from line sources, and for highly forced and highly lazy releases from sources of finite width in a weak co-flow. Contrary to releases in quiescent surroundings, our solutions show that all classes of release can exhibit plume contraction and the associated necking. For entraining plumes, a dynamical invariance spatially only occurs for pure and forced releases and we derive the co-flow strengths that lead to this invariance.

Roll waves produced by a local disturbance comprise a group of shock waves with steep fronts. We used a robust and accurate numerical scheme to capture the steep fronts in a shallow-water hydraulic model of the waves. Our simulations of the waves in clear water revealed the existence of a front runner with an exceedingly large amplitude – much greater than those of all other shock waves in the wave group. The trailing waves at the back remained periodic. Waves were produced continuously within the group due to nonlinear instability. The celerity depended on the wave amplitude. Over time, the instability produced an increasing number of shock waves in an ever-expanding wave group. We conducted simulations for three types of local disturbances of very different duration over a range of amplitudes. We interpreted the simulation results for the front runner and the trailing waves, guided by an analytical solution and the laboratory data available for the smaller waves in the trailing end of the wave group.

The present work studies the nonlinear dynamics of a shear layer, driven by a body force and confined between parallel walls, a simplified setting to study transitional and turbulent shear layers. It was introduced by Nogueira & Cavalieri (J. Fluid Mech., vol. 907, 2021, A32), and is here studied using a reduced-order model based on a Galerkin projection of the Navier–Stokes system. By considering a confined shear layer with free-slip boundary conditions on the walls, periodic boundary conditions in streamwise and spanwise directions may be used, simplifying the system and enabling the use of methods of dynamical systems theory. A basis of eight modes is used in the Galerkin projection, representing the mean flow, Kelvin–Helmholtz vortices, rolls, streaks and oblique waves, structures observed in the cited work, and also present in shear layers and jets. A dynamical system is obtained, and its transition to chaos is studied. Increasing Reynolds number $Re$ leads to pitchfork and Hopf bifurcations, and the latter leads to a limit cycle with amplitude modulation of vortices, as in the direct numerical simulations by Nogueira & Cavalieri. Further increase of $Re$ leads to the appearance of a chaotic saddle, followed by the emergence of quasi-periodic and chaotic attractors. The chaotic attractors suffer a merging crisis for higher $Re$, leading to a chaotic dynamics with amplitude modulation and phase jumps of vortices. This is reminiscent of observations of coherent structures in turbulent jets, suggesting that the model represents a dynamics consistent with features of shear layers and jets.

This work studies the application of a reinforcement learning (RL)-based flow control strategy to the flow past a cylinder confined between two walls to suppress vortex shedding. The control action is blowing and suction of two synthetic jets on the cylinder. The theme of this study is to investigate how to use and embed physical information of the flow in the RL-based control. First, global linear stability and sensitivity analyses based on the time-mean flow and the steady flow (which is a solution to the Navier–Stokes equations) are conducted in a range of blockage ratios and Reynolds numbers. It is found that the most sensitive region in the wake extends itself when either parameter increases in the parameter range we investigated here. Then, we use these physical results to help design RL-based control policies. We find that the controlled wake converges to the unstable steady base flow, where the vortex shedding can be successfully suppressed. A persistent oscillating control seems necessary to maintain this unstable state. The RL algorithm is able to outperform a gradient-based optimisation method (optimised in a certain period of time) in the long run. Furthermore, when the flow stability information is embedded in the reward function to penalise the instability, the controlled flow may become more stable. Finally, according to the sensitivity analyses, the control is most efficient when the probes are placed in the most sensitive region. The control can be successful even when few probes are properly placed in this manner.

We present new experiments of particle-driven turbulent plumes issuing from a constant source of dense particle-laden fluid, with buoyancy flux, $B$, in a uniform horizontal current, $u$. Experiments show that a turbulent, well-mixed plume develops, in which the downward vertical speed $w$ decreases with depth $z$ according to $w = 0.76 (B/uz)^{1/2}$ while the horizontal speed rapidly asymptotes to the current speed $u$, provided that the Stokes settling speed of the particles $v<0.92 w$. For $v > 0.92 w$, the particles separate from the plume fluid, and their depth $z$ increases according to the simple sedimentation trajectory $\textrm {d}z/{\textrm {d}\kern0.7pt x} = v/u$. As the particles sediment, they form clusters of particles, which lead to fluctuations in the particle load with position, but do not appear to change the time-average sedimentation speed. We explore the impact of these results for deep-sea mining, in which the fate of the plume water as well as the particles is key for assessing potential environmental impacts.

Nonlinear Hall-magnetohydrodynamic dynamos associated with coherent structures in subcritical shear flows are investigated by using unstable invariant solutions. The dynamo solution found has a relatively simple structure, but it captures the features of the typical nonlinear structures seen in simulations, such as current sheets. As is well known, the Hall effect destroys the symmetry of the magnetohydrodynamic equations and thus modifies the structure of the current sheet and mean field of the solution. Depending on the strength of the Hall effect, the generation of the magnetic field changes in a complex manner. However, a too strong Hall effect always acts to suppress the magnetic field generation. The hydrodynamic/magnetic Reynolds number dependence of the critical ion skin depth at which the dynamos start to feel the Hall effect is of interest from an astrophysical point of view. An important consequence of the matched asymptotic expansion analysis of the solution is that the higher the Reynolds number, the smaller the Hall current affects the flow. We also briefly discuss how the above results for a relatively simple shear flow can be extended to more general flows such as infinite homogeneous shear flows and boundary layer flows. The analysis of the latter flows suggests that interestingly a strong induction of the generated magnetic field might occur when there is a background shear layer.

Direct numerical simulation is used to investigate the integral behaviour of buoyant plumes subjected to a uniform crossflow that are infinitely lazy at the source. Neither a plume trajectory defined by the centre of mass of the plume $z_c$ nor a trajectory defined by the central streamline $z_{U}$ is aligned with the average streamlines inside the plume. Both $z_c$ and $z_{U}$ are shown to correlate with field lines of the total buoyancy flux, which implies that a model for the vertical turbulent buoyancy flux is required to faithfully predict the plume angle. A study of the volume conservation equation shows that entrainment due to incorporation of ambient fluid with non-zero velocity due to the increase in the surface area (the Leibniz term) is the dominant entrainment mechanism in strong crossflows. The data indicate that pressure differences between the top and bottom of the plume play a leading role in the evolution of the horizontal and vertical momentum balances and are crucial for appropriately modelling plume rise. By direct parameterisation of the vertical buoyancy flux, the entrainment and the pressure, an integral plume model is developed which is in good agreement with the simulations for sufficiently strong crossflow. A perturbation expansion shows that the current model is an intermediate-range model valid for downstream distances up to $100\ell _b$–$1000 \ell _b$, where $\ell _b$ is the buoyancy length scale based on the flow speed and plume buoyancy flux.

In this research, highly disturbed turbulent flow of distinct three-dimensional characteristics in a square duct with inclined or V-shaped ribs mounted on one wall is investigated using direct numerical simulation. The turbulence field is highly sensitive to not only the rib geometry but also the boundary layers developed over the side and top walls. In a cross-stream plane secondary flows appear as large longitudinal vortices in both inclined and V-shaped rib cases due to the confinement of four sidewalls of the square duct. However, owing to the difference in the pattern of cross-stream secondary flow motions, the flow physics is significantly different in these two ribbed duct cases. It is observed that the mean flow structures in the cross-stream directions are asymmetrical in the inclined rib case but symmetrical in the V-shaped rib case, causing substantial differences in the momentum transfer across the spanwise direction. The impacts of rib geometry on near-wall turbulence structures are investigated using vortex identifiers, joint probability density functions between the streamwise and vertical velocity fluctuations, statistical moments of different orders, spatial two-point autocorrelations and velocity spectra. It is found that near the leeward and windward rib faces, the mean inclination angle of turbulence structures in the V-shaped rib case is greater than that of the inclined rib case, which subsequently enhances momentum transport between the ribbed bottom wall and the smooth top wall.

It is known that the dripping of a liquid film on the underside of a plate can be suppressed by tilting the plate so as to cause a sufficiently strong flow. This paper uses two-dimensional numerical simulations in a closed-flow framework to study several aspects of this phenomenon. It is shown that, in quasi-equilibrium conditions, the onset of dripping is closely associated with the curvature of the wave crests approaching a well-defined maximum value. When dynamic effects become significant, this connection between curvature and dripping weakens, although the critical curvature remains a useful reference point as it is intimately related to the short length scales promoted by the Rayleigh–Taylor instability. In the absence of flow, when the film is on the underside of a horizontal plate, the concept of a limit curvature is relevant only for small liquid volumes close to a critical value. Otherwise, the drops that form have a smaller curvature and a large volume. The paper also illustrates the peculiarly strong dependence of the dripping transition on the initial conditions of the simulations. This feature prevents the development of phase maps dependent only on the governing parameters (Reynolds number, Bond number, etc.) similar to those available for film flow on the upper side of an inclined plate.

Breaking nonlinear internal waves (NLIWs) of depression on boundary slopes drives mixing in the coastal ocean. Of the different breaker types, fission is most commonly observed on mild slopes of continental margins. However, fission on mild slopes has rarely been investigated in the laboratory owing to limitations on flume length. In the present work, a train of NLIWs of depression is generated in an 18.2 m wave flume and shoaled upon a mild uniform slope. During fission, each NLIW of depression scatters into one or two NLIWs of elevation, which transforms into a bolus at the bolus birth point, where shear instability occurs through the pycnocline. The bolus propagates upslope, decreasing in size until it degenerates by shear and lobe-cleft instability, while losing volume to a return flow along the bed. The location of the bolus birth point, bolus propagation length scale, initial size and the number of boluses from each incident wave are parameterized from the wave half-width and the wave Froude number associated with the incident NLIW. These are compared with the characteristics of boluses generated by other breaking mechanisms on steeper slopes. Some bolus characteristics (height to length ratio, change in size and velocity field) are similar for boluses generated by fission, collapsing sinusoidal waves and internal solitary waves of elevation; however, the number of boluses, their birth point and initial height differ. The boluses formed by fission have more initial energy and no reflection. Further research is required to better quantify bolus-driven mixing on continental margins.

Non-isothermal reactive transport in complicated porous media is diverse in nature and industrial applications. There are challenges in the modelling of multiple physicochemical processes in multiscale pore structures with various length scales ranging from nanometres to micrometres. This study focuses on coke combustion during in situ crude oil combustion techniques. A micro-continuum model was developed to perform an image-based simulation of coke combustion through a multiscale porous medium. The simulation coupled weakly compressible gas flow, species transport, conjugate heat transfer, heterogeneous coke oxidation kinetics and structural evolution. The unresolved nanoporous coke region was treated as a continuum, for which the random pore model, permeability model and species diffusivity model were integrated as sub-grid models to account for the sub-resolution reactive surface area, Darcy flow and Knudsen diffusion, respectively. A Pe–Da diagram was provided to present five characteristic combustion regimes covering the ignition temperature and air flux in realistic field operations and laboratory measurements. The present model proved to achieve more accurate predictions of the feasible ignition temperature than previous models. Compared with the air flux of $\phi \sim O\textrm{(1) s}{\textrm{m}^\textrm{3}}(\textrm{air})\;{({\textrm{m}^\textrm{2}}\ \textrm{h})^{ - 1}}$ in the field, the increasing air flux in the laboratory transformed the combustion regime from diffusion-limited to convection-limited, which led to an overpredicted burning temperature. Reactive fingering combustion was analysed to understand the potential risks in some experimental measurements. The findings provide a better understanding of coke combustion and can help engineers design sustainable combustion methods. The developed image-based model allows other types of multiscale and nonlinear reactive transport to be simulated.

A focusing shock wave can be generated during the high-speed impact of a droplet on a $180^\circ$ constrained wall, which can be used to realise energy convergence on a small scale. In this study, to realise high energy convergence and peak pressure amplification, a configuration of droplets embedded with cavities is proposed for high-speed impingement on a $180^\circ$ constrained wall. A multicomponent two-phase compressible flow model considering the phase transition is used to simulate the high-speed droplet impingement process. The properties of the embedded cavities can influence the collapse pressure peak. The collapse of an embedded single air cavity or vapour cavity, as well as the cavities in a tandem array, is simulated in this study. The physical evolution mechanisms of the impinging droplet and the embedded cavities are investigated qualitatively and quantitatively by characterising the focusing shock wave generated inside the droplet and its interaction with different cavity configurations. The interaction dynamics between the cavities is analysed and a theoretical prediction model for the intensity of each cavity collapse in the tandem array is established. With the help of this theoretical model, the influencing factors for the collapse intensities of the tandem cavities are identified. The results reveal that the properties of the initial shock wave and the interval between the cavities are two predominant factors for the amplification of the collapse intensity. This study enhances the understanding of the physical process of shock-induced tandem-cavity collapse.

In this study direct numerical simulations are employed to investigate the effects of various parameters on the laminar-flow-control capabilities of narrowly spaced streaks in a supersonic boundary layer at Mach $2.0$. Previous work by Sharma et al. (J. Fluid Mech., vol. 873, 2019, pp. 1072–1089) has found these streak modes, excited by a spanwise blowing-and-suction strip, to be highly effective at delaying pure oblique-type breakdown. In the present work it is shown that spectrum-enriching subharmonic modes, relevant with increasing running-length Reynolds number, do not destroy the controlling mechanism, and also a complex breakdown scenario, triggered by a multi-frequency point source, is found to be effectively controlled. Moreover, the control-streak excitation by roughness elements is compared in detail with the blowing-and-suction method, revealing relevant differing features.

A strongly nonlinear theory describing the effect of small amplitude boundary forcing in the form of waves on high Reynolds number shear flows is given. The interaction leads to an $O(1)$ change in the unperturbed flow and is relevant to a number of forcing mechanisms. The cases of the shear flow being bounded or unbounded are both considered and the results for the unbounded case apply to quite arbitrary flows. The instability criterion for unbounded flows is expressed in terms of the wall forcing and the friction Reynolds number. As particular examples we investigate wall transpiration or surface undulations as sources of the forcing and both propagating and stationary waves are considered. Results are given for propagating waves with crests perpendicular to the flow direction and for stationary waves with crests no longer perpendicular to the flow direction. In the first of those situations we find the instability induced by transpiration waves is independent of the propagation speed. For wavy walls downstream propagation completely stabilises the flow at a critical speed whereas upstream propagation greatly destabilises the flow. For stationary oblique waves we find that the instability is enhanced and a much wider range of unstable wavenumbers exists. For the bounded case with a wall of fixed wavelength we identify a critical wavelength where the most dangerous mode switches from the aligned to the oblique configuration. For the transpiration problem in the oblique configuration a strong resonance occurs when the vortex wavelength coincides with the spanwise wavelength of the forcing.

The effects of along-shelf barotropic geostrophic currents on internal wave generation by the $K_1$ tide interacting with a shelf at near-critical latitudes are investigated. The horizontal shear of the background current results in a spatially varying effective Coriolis frequency which modifies the slope criticality and potentially creates blocking regions where freely propagating internal tides cannot exist. This paper is focused on the barotropic to baroclinic energy conversion rate, which is affected by a combination of three factors: slope criticality, size and location of the blocking region where the conversion rate is extremely small and the internal tide (IT) beam patterns. All of these are sensitive to the current parameters. In our parameter space, the current can increase the conversion rate up to 10 times.