Undulating and breaking bores are generated in the laboratory using a programmable long-stroke wavemaker. By changing the stroke length and the speed of the wavemaker, both non-decaying and decaying bores are generated and studied. Bore strength, height and duration are measured and compared with the solutions derived by using the method of characteristics, with excellent agreement. The measurements for inundation depth, runup height and flood duration are checked with the formulas presented in Barranco & Liu (J. Fluid Mech., vol. 915, 2021). The comparisons show that the formulas are also accurate for the non-decaying bores generated by the wavemaker. The maximum inundation depth predicted by the formula for zero bore length at the beach toe agrees with the laboratory observations for decaying bores. Using a high-speed particle image velocimetry system, the ensemble-averaged velocities and fluctuating velocities under undulating bores and breaking bores are measured in constant water depth and in the vicinity of the still water shoreline. Detailed analyses of the velocity fields are presented and discussed. For the undulating bore a long quiescent flood duration is observed, while for the breaking bore the up-rush flow changes into down-rush flow almost linearly.

Bubble–particle collisions in turbulence are central to a variety of processes such as froth flotation. Despite their importance, details of the collision process have not received much attention yet. This is compounded by the sometimes counter-intuitive behaviour of bubbles and particles in turbulence, as exemplified by the fact that they segregate in space. Although bubble–particle relative behaviour is fundamentally different from that of identical particles, the existing theoretical models are nearly all extensions of theories for particle–particle collisions in turbulence. The adequacy of these theories has yet to be assessed as appropriate data remain scarce to date. In this investigation, we study the geometric collision rate by means of direct numerical simulations of bubble–particle collisions in homogeneous isotropic turbulence using the point-particle approach over a range of the relevant parameters, including the Stokes and Reynolds numbers. We analyse the spatial distribution of bubble and particles, and quantify to what extent their segregation reduces the collision rate. This effect is countered by increased approach velocities for bubble–particle compared to monodisperse pairs, which we relate to the difference in how bubbles and particles respond to fluid accelerations. We found that in the investigated parameter range, these collision statistics are not altered significantly by the inclusion of a lift force or different drag parametrisations, or when assuming infinite particle density. Furthermore, we critically examine existing models and discuss inconsistencies therein that contribute to the discrepancy.

Recent studies of pressure-driven flows of dilute polymer solutions in straight channels demonstrated the existence of two-dimensional coherent structures that are disconnected from the laminar state and appear through a subcritical bifurcation from infinity. These travelling-wave solutions were suggested to organise the phase-space dynamics of purely elastic and elasto-inertial chaotic channel flows. Here, we consider a wide range of parameters, covering the purely elastic and elasto-inertial cases, and demonstrate that the two-dimensional travelling-wave solutions are unstable when embedded in sufficiently wide three-dimensional domains. Our work demonstrates that studies of purely elastic and elasto-inertial turbulence in straight channels require three-dimensional simulations, and no reliable conclusions can be drawn from studying strictly two-dimensional channel flows.

The influence of viscosity on the Mach reflection of shock waves in a steady flow of a monatomic gas is studied by solving the Navier–Stokes equations numerically. Based on the nested block grid refinement technique, the flow near the shock wave intersection is simulated, and its behaviour with increasing Reynolds number is studied. The computations are performed for the interaction of both strong (free-stream Mach number $M_\infty = 4$) and weak ($M_\infty = 1.7$) shock waves. In the strong reflection of shock waves at all Reynolds numbers in the examined range, it is found that there exists a small-size zone behind the shock wave intersection where the flow parameters differ from those predicted by the Rankine–Hugoniot relations and hence deviate from the predictions of the inviscid three-shock theory. The structure of this zone is self-similar: in coordinates normalised to the mean free path of molecules in the free stream. The structure is identical at all Reynolds numbers considered in the study. As the Reynolds number increases, the size of this zone in physical coordinates decreases, but the maximum difference between the viscous and inviscid solutions in this zone remains constant, reaching approximately $10\,\%$ for pressure. In the weak reflection of shock waves, the flow structure behind the shock wave intersection is not self-similar, i.e. the flow fields at different Reynolds numbers do not coincide in the normalised coordinates, but converge, as the Reynolds number increases, to the parameters predicted by the inviscid three-shock theory.

The presence of very-large-scale motions in wall-bounded turbulent flows is commonly associated with their footprint in the form of the superposition of the large scales at the wall and the additional amplitude modulation of small-scale near-wall turbulence. These two phenomena are currently understood to be interlinked, with the superposed large-scale velocity gradient causing the modulation of small-scale activity in the proximity of the wall. To challenge this idea, we devise a numerical strategy that selectively suppresses either superposition or amplitude modulation, in an effort to isolate and study the remaining phenomenon. Results from our direct numerical simulations indicate that a positive correlation between the amplitude of the small scales in the near-wall region and the large-scale signal in the outer flow persists even when near-wall large-scale motions are suppressed – i.e. in absence of superposition. Clearly, this kind of correlation cannot be caused by the near-wall large-scale velocity or its gradients, as both are absent. Conversely, when modulation is blocked, the near-wall footprints of the large scales seem to disappear. This study has been carried out on channel flows at friction Reynolds number $Re_\tau =1000$ in both standard simulation domains and minimal streamwise units (MSUs), where the streamwise fluctuation energy is enhanced. The consistency of the results obtained by the two approaches suggests that MSUs can capture correctly this kind of scale interaction at a much reduced cost.

The transition to dripping in the gravity-driven flow of a liquid film under an inclined plate is investigated at zero Reynolds number. Computations are carried out on a periodic domain assuming either a fixed fluid volume or a fixed flow rate for a hierarchy of models: two lubrication models with either linearised curvature or full curvature (the LCM and FCM, respectively), and the full equations of Stokes flow. Of particular interest is the breakdown of travelling-wave solutions as the plate inclination angle is increased. For any fixed volume, the LCM reaches the horizontal state where it attains a cosine-shaped profile. For sufficiently small volume, the FCM and Stokes solutions attain a weak Young–Laplace equilibrium profile, the approach to which is described by an asymptotic analysis generalising that of Kalliadasis & Chang (J. Fluid Mech., vol. 261, 1994, pp. 135–168) for the LCM. For large volumes, the bifurcation curves for the FCM and Stokes model have a turning point so that the fully inverted state is never reached. For fixed flow rate, the LCM blows up at a critical angle that is well predicted by asymptotic analysis. The bifurcation curve for the FCM either has a turning point or else reaches a point at which the surface profile has an infinite slope singularity, indicating the onset of multi-valuedness. The latter is confirmed by the Stokes model, which can be continued to obtain overturning surface profiles. Overall, the thin-film models either provide an accurate prediction for dripping onset or else supply an upper bound on the critical inclination angle.

Machine learning frameworks such as genetic programming and reinforcement learning (RL) are gaining popularity in flow control. This work presents a comparative analysis of the two, benchmarking some of their most representative algorithms against global optimization techniques such as Bayesian optimization and Lipschitz global optimization. First, we review the general framework of the model-free control problem, bringing together all methods as black-box optimization problems. Then, we test the control algorithms on three test cases. These are (1) the stabilization of a nonlinear dynamical system featuring frequency cross-talk, (2) the wave cancellation from a Burgers’ flow and (3) the drag reduction in a cylinder wake flow. We present a comprehensive comparison to illustrate their differences in exploration versus exploitation and their balance between ‘model capacity’ in the control law definition versus ‘required complexity’. Indeed, we discovered that previous RL control attempts of controlling the cylinder wake were performing linear control and that the wide observation space was limiting their performances. We believe that such a comparison paves the way towards the hybridization of the various methods, and we offer some perspective on their future development in the literature of flow control problems.

Nucleate boiling, a ubiquitous heat transfer mode, involves multiple vapour bubble nucleations on the heater surface and offers high heat transfer coefficients. The bubble growth process on a heating substrate involves the formation of microlayer, a thin liquid film trapped between the growing bubble and the heating substrate, and contributes to the bubble growth phenomenon through evaporation. Microlayer dynamics for a single bubble have been widely investigated in the pool and flow boiling conditions. However, the literature on multiple bubbles interactions and their associated microlayer information is scarce. Notably, in the case of flow boiling, where the microlayer dynamics are not symmetric due to the bubble's movement, the bubbles’ interaction and its influence on associated microlayer dynamics have never been reported. Therefore, microlayer and bubble dynamics in multiple interactions have been investigated experimentally using simultaneous application of thin-film interferometry and high-speed videography techniques in flow boiling with water as the working fluid. Our experimental investigation revealed that the secondary nucleation could cause a reduction in lift-off time and may assist or hinder the movement of the first bubble. The experimental results also demonstrated that the secondary nucleation could deplete the microlayer of the first bubble hydrodynamically even when the bubbles are far apart. Furthermore, it has been found that the microlayer depletion rate depends on the growth rate and the location of the secondary nucleation. Hence, this experimental study emphasises the need to consider the interaction of bubbles while modelling boiling flows to avoid overestimating the contribution of microlayer evaporation.

We study the linear stability of bubbles in a capillary tube under external flow. Yu et al. (J. Fluid Mech., vol. 911, 2021, pp. 1–19) showed that a rich variety of bubble dynamics occurs when a downward external flow is applied, opposing the buoyancy-driven ascent of the bubble. They found experimentally and numerically the existence of two branches of solutions that overlap over a finite range of the capillary number of the downward external flow in cases where the Reynolds number is small and the Bond number is larger than the critical value for which the bubble can rise spontaneously (Bretherton, J. Fluid Mech., vol. 10, issue 2, 1961, pp. 166–188). Furthermore, inertialess, symmetry-breaking steady-state shapes were found as the bubble transits near the tipping points of the solution branches. In this work, using steady axisymmetric simulations, we show that the reported multiplicity of solutions can be described using bifurcation diagrams with three branches of steady axisymmetric solutions and two limit points. The linear global stability analysis of the different branches of the stationary axisymmetric solutions demonstrates that the symmetry breaking is due to the development of three-dimensional instabilities with azimuthal wavenumber $|m|=1$.

This paper presents a numerical investigation of the effects of the moving ground and rotating wheels on the turbulent flow around a 1/10 scaled square-back van model. A comprehensive comparison among the partially averaged Navier–Stokes (PANS), large eddy simulation (LES) and particle image velocimetry (PIV) involving the aerodynamic drag, the wake topology, the velocity and the Reynolds stress profiles in the wake region is conducted. The proper orthogonal decomposition (POD) and fast Fourier transform (FFT) are applied to the shear layers shedding from the trailing edges to comment on the coherent structures and their frequency content. The Reynolds number for both simulations and experiments is set to Re = 2.5 × 105 based on the inlet velocity ${U_{inf}} = 9\;\textrm{m}\;{\textrm{s}^{ - 1}}$ and the width of the model W = 0.17 m. The results show that PANS accurately predicts the flow field measured in experiments and predicted by a resolved LES, even with a low-resolution grid. The superiority of the PANS approach could provide good guidance for industrial research in predicting the turbulent flow around the square-back van model with affordable computational grids. The ground and wheel motion mechanism on the aerodynamic forces has been revealed by analysing the surface pressure distribution, the wheels’ surrounding flow, the underbody flow characteristics and the turbulent wake structures. The effects of the ground and wheel motion on the frequency, evolution and development characteristics of the wake shear layers are analysed, thus providing relevant insights for future experimental investigations of square-back van models.

The dynamics of hydrogen bubbles produced by water electrolysis in an acidic electrolyte is studied using electrochemical and optical methods. A defined cyclic modulation of the electric potential is applied at a microelectrode to produce pairs of interacting H$_2$ bubbles in a controlled manner. Three scenarios of interactions are identified and studied systematically. The most prominent one consists of a sudden reversal in the motion of the first detached bubble, its return to the electrode, and finally its coalescence with the second bubble. Attested by Toepler's schlieren technique, an explanation of contactless motion reversal is provided by the competition between buoyancy and thermocapillary effects.