A hyperbolic two-fluid model for gas–particle flow derived using the Boltzmann–Enskog kinetic theory is generalized to include added mass. In place of the virtual-mass force, to guarantee indifference to an accelerating frame of reference, the added mass is included in the mass, momentum and energy balances for the particle phase, augmented to include the portion of the particle wake moving with the particle velocity. The resulting compressible two-fluid model contains seven balance equations (mass, momentum and energy for each phase, plus added mass) and employs a stiffened-gas model for the equation of state for the fluid. Using Sturm's theorem, the model is shown to be globally hyperbolic for arbitrary ratios of the material densities $Z = \rho _f / \rho _p$ (where $\rho _f$ and $\rho _p$ are the fluid and particle material densities, respectively). An eight-equation extension to include the pseudo-turbulent kinetic energy (PTKE) in the fluid phase is also proposed; however, PTKE has no effect on hyperbolicity. In addition to the added mass, the key physics needed to ensure hyperbolicity for arbitrary $Z$ is a fluid-mediated contribution to the particle-phase pressure tensor that is taken to be proportional to the volume fraction of the added mass. A numerical solver for hyperbolic equations is developed for the one-dimensional model, and numerical examples are employed to illustrate the behaviour of solutions to Riemann problems for different material-density ratios. The relation between the proposed two-fluid model and prior work on effective-field models is discussed, as well as possible extensions to include viscous stresses and the formulation of the model in the limit of an incompressible continuous phase.

Predicting and controlling entropy-wave-driven combustion instabilities is challenging, because the production, advection and dispersion of entropy waves in practical systems is difficult to model. The present paper aims to shed new light on this problem by considering a highly turbulent configuration with experiments and large eddy simulations. In this configuration, the decay of entropy waves is not only governed by the shear dispersion of an idealized turbulent pipe flow, as assumed in the recent studies on the topic, but also enhanced by the highly three-dimensional dispersion due to large-scale coherent structures. A novel post-processing approach for background-oriented schlieren (BOS) thermometry is proposed and enables the measurement of entropy waves with high spatio-temporal resolution for wide ranges of entropy wave amplitudes and frequencies. Instantaneous BOS snapshots are arranged in a panoramic coordinate frame using velocity data. This work, therefore, contributes to filling the knowledge gap in experimental data on entropy waves. The new dataset is accompanied with large eddy simulations to further elucidate the mechanisms dominating the amplitude decay of the entropy waves. It is shown that shear dispersion models that are only based on mean profiles of the axial velocity significantly underestimate the decay in the present configuration, where the entropy waves are formed by periodically injecting hot air pockets in the main stream. It is shown that the turbulent nature of the coherent hot pockets plays a key role in the dispersion enhancement.

We consider the dynamics of a vertically stratified, horizontally forced Kolmogorov flow. Motivated by astrophysical systems where the Prandtl number is often asymptotically small, our focus is the little-studied limit of high Reynolds number but low Péclet number (which is defined to be the product of the Reynolds number and the Prandtl number). Through a linear stability analysis, we demonstrate that the stability of two-dimensional modes to infinitesimal perturbations is independent of the stratification, whilst three-dimensional modes are always unstable in the limit of strong stratification and strong thermal diffusion. The subsequent nonlinear evolution and transition to turbulence are studied numerically using direct numerical simulations. For sufficiently large Reynolds numbers, four distinct dynamical regimes naturally emerge, depending upon the strength of the background stratification. By considering dominant balances in the governing equations, we derive scaling laws for each regime which explain the numerical data.

Ekman layers over a rough surface are studied using direct numerical simulation. The roughness takes the form of periodic two-dimensional bumps whose non-dimensional amplitude is fixed at a small value ($h^+=15$) and whose mean slope is gentle. The neutral Ekman layer is subjected to a stabilizing cooling flux for approximately one inertial period ($2{\rm \pi} /f$) to impose the stratification. The Ekman boundary layer is in a transitionally rough regime and, without stratification, the effect of roughness is found to be mild in contrast to the stratified case. Roughness, whose effect increases with the slope of the bumps, changes the boundary layer qualitatively from the very stable (Mahrt, Theor. Comput. Fluid Dyn., vol. 11, issue 3–4, 1998, pp. 263–279) regime, which has a strong thermal inversion and a pronounced low-level jet, in the flat case to the stable regime, which has a weaker thermal inversion and stronger surface-layer turbulence, in the rough cases. The flat case exhibits initial collapse of turbulence which eventually recovers, albeit with inertial oscillations in turbulent kinetic energy. The roughness elements interrupt the initial collapse of turbulence. In the quasi-steady state, the thickness of the turbulent stress profiles and of the near-surface region with subcritical gradient Richardson number increase in the rough cases. Analysis of the turbulent kinetic energy (TKE) budget shows that, in the surface layer, roughness counteracts the stability-induced reduction of TKE production. The flow component, coherent with the surface undulations, is extracted by a triple decomposition, and leads to a dispersive component of near-surface turbulent fluxes. The significance of the dispersive component increases in the stratified cases.

In this work, we study the ascent dynamics of a liquid Taylor drop formed from a lock-exchange configuration in a closed vertical pipe. We focus on the buoyancy-driven motion of an elongated drop surrounded by a denser fluid when viscous forces dominate over inertial and surface tension effects. While gaseous Taylor bubbles have been studied extensively, a liquid Taylor drop moving in a closed pipe is less well understood. We formulate an analytical model for estimating the ascent speed and drop thickness from first principles. First, we use a lubrication approximation to solve for the velocity profiles in the two fluids. Then, we analyse the mechanical energy balance of the whole system, including the effect of viscous dissipation, to understand how the ascent speed and drop thickness scale with the viscosity ratio. We show that a drop with density ratio $\mathcal {R}$ reaches a stationary state with a uniform dimensionless thickness of $\sqrt {2}/{2}$ in the absence of dissipation and $\sqrt {2\mathcal {R}}/{2}$ in the dissipative regime. Through a comparison with existing experimental data, we demonstrate that our model correctly predicts the ascent speed of a Taylor drop if the material properties of the fluids and the geometry of the conduit are known. Our theoretical framework can be generalized to an isolated Taylor drop rising in a vertical pipe.

Previous experimental (Kühnen et al., Flow Turb. Combust., vol. 100, 2018, pp. 919–943) and numerical (Marensi et al., J. Fluid Mech., vol. 863, 2019, pp. 850–875) studies have demonstrated that a streamwise-localised baffle can fully relaminarise pipe flow turbulence at Reynolds numbers of $O(10^4)$. Optimising the design of the baffle involves tackling a complicated variational problem built around time stepping turbulent solutions of the Navier–Stokes equations which is difficult to solve. Here instead, we investigate a much simpler ‘spectral’ approach based upon maximising the energy stability of the baffle-modified laminar flow. The ensuing optimal problem has much in common with the variational procedure to derive an upper bound on the energy dissipation rate in turbulent flows (e.g. Plasting & Kerswell, J. Fluid Mech., vol. 477, 2003, pp. 363–379) so well-honed techniques developed there can be used to solve the problem here. The baffle is modelled by a linear drag force $-F(\boldsymbol {x}) \boldsymbol {u}$ (with $F(\boldsymbol {x}) \ge 0 \ \forall \boldsymbol {x}$) where the extent of the baffle is constrained by an $L_{\alpha }$ norm with various choices explored in the range $1 \leq \alpha \leq 2$. An asymptotic analysis demonstrates that the optimal baffle is always axisymmetric and streamwise independent, retaining just radial dependence. The optimal baffle which emerges in all cases has a similar structure to that found to work in experiments: the baffle retards the flow in the pipe centre causing the flow to become faster near the wall thereby reducing the turbulent shear there. Numerical simulations demonstrate that the designed baffle can relaminarise turbulence efficiently at moderate Reynolds numbers ($Re \le 3500$), and an energy saving regime has been identified. Direct numerical simulation at $Re=2400$ also demonstrates that the drag reduction can be realised by truncating the energy-stability-designed baffle to finite length.

We investigate the spatial organization and temporal dynamics of large-scale, coherent structures in turbulent Rayleigh–Bénard convection via direct numerical simulation of a 6.3 aspect-ratio cylinder with Rayleigh and Prandtl numbers of $9.6\times 10^7$ and $6.7$, respectively. Fourier modal decomposition is performed to investigate the structural organization of the coherent turbulent motions by analysing the length scales, time scales and the underlying dynamical processes that are ultimately responsible for the large-scale structure formation and evolution. We observe a high level of rotational symmetry in the large-scale structure in this study and that the structure is well described by the first four azimuthal Fourier modes. Two different large-scale organizations are observed during the duration of the simulation and these patterns are dominated spatially and energetically by azimuthal Fourier modes with frequencies of 2 and 3. Studies of the transition between these two large-scale patterns, radial and vertical variations in the azimuthal energy spectra, as well as the spatial and modal variations in the system's correlation time are conducted. Rotational dynamics are observed for individual Fourier modes and the global structure with strong similarities to the dynamics that have been reported for unit aspect-ratio domains in prior works. It is shown that the large-scale structures have very long correlation time scales, on the order of hundreds to thousands of free-fall time units, and that they are the primary source for a horizontal inhomogeneity within the system that can be observed during a finite, but a very long-time simulation or experiment.

The present paper aims at establishing the synthetic jet technology capabilities in controlling the von Kármán street behind a circular cylinder. The circular cylinder, placed in an open-circuit wind tunnel, presents a slot in its rear position, through which the synthetic jet is issued. The Reynolds number, based on the circular cylinder diameter and the free-stream velocity, is equal to 4600 and the von Kármán street is characterized, in the baseline configuration (i.e. without synthetic jet), by a shedding frequency of 16.2 Hz. Several synthetic jet operating conditions are tested. Therefore, the effects of the momentum coefficient ($C_{\mu } = 5.4$%, 10.8% and 21.6%) and the dimensionless frequency ($f^{+} = 0.49$, 0.98 and 1.96) on the von Kármán street behaviour can be analysed. Instantaneous two-dimensional in-plane velocity fields are measured in a plane containing the synthetic jet slot axis using multigrid/multipass cross-correlation digital particle image velocimetry. These measurements have been used to investigate the mean flow quantities and turbulent statistics of the phenomenon. In addition, the wake extent and behaviour (i.e. symmetric or asymmetric) are analysed as well as the drag coefficient, for each configuration. The extent of the wake region decreases as the momentum coefficient and/or the dimensionless frequency increase, while the symmetric/asymmetric wake behaviour is found to be governed by a different control parameter: the synthetic jet Reynolds number based on its impulse. As regards the drag coefficient, a maximum reduction, of approximately 35%, is found for the configuration at $C_{\mu }=10.8\,\%$ and $f^{+}=0.98$.

A numerical model that allows one to study numerically the evolution of waves along the test section of a wind-wave tank is offered. The simulations are directly related to wind-wave tank experiments carried out for a range of steady wind velocities. At each wind forcing condition, the evolving wind-wave field is strongly non-homogeneous, with wave energy growth along the test section accompanied by frequency downshifting. The wave parameters measured at a short fetch serve as a basis for generating numerous realizations of the initial conditions in the Monte Carlo numerical simulations. The computations are based on a modified unidirectional spatial version of the Zakharov equation that accounts for wind input and dissipation and is applicable for the whole range of wind velocities employed. The model contains two empirical parameters that are selected by comparison of the experimental and numerical results; the same values of those parameters are applied for all wind forcing conditions. The availability of an experimentally verified numerical model allows one to study the contributions of nonlinear wave–wave interactions, dissipation and wind input separately. Special attention is given to accounting for the three-dimensional and random nature of wind waves as observed in experiments. The suggested model combines approaches adopted in the wind-wave growth theories by Miles and Phillips.

This paper provides an experimental investigation on the internal shear layers and the edges of the uniform momentum zones (UMZs) in a turbulent pipe flow. The time-resolved stereoscopic particle image velocimetry data are acquired in the cross-section of the pipe, and span the range of Reynolds number $\textit {Re}_\tau =340\text {--}1259$. In the first part of the study, internal shear layers are detected using a three-dimensional detection method, and both their geometry as well as their fingerprint in the flow statistics are examined. Three-dimensional conditional mean flow analysis revealed a strong low-speed region beneath the average shear layers. This low-speed region is associated with positive wall-normal fluctuations, and it is accompanied by two swirling motions having opposite signs on either side in the azimuthal direction. Moreover, the shear layers are stretched by the two opposite azimuthal motions. In the second part of the study, the shear layers are treated as the continuous edges of the UMZs, which are detected using the histogram method following Adrian et al. (J. Fluid Mech., vol. 422, 2000, pp. 1–54) and de Silva et al. (J. Fluid Mech., vol. 786, 2016, pp. 309–331). For this part, two different orientation of the planes are used, i.e. the wall-normal–streamwise plane and the wall-normal–spanwise plane (cross-section of the pipe). Comparison of the detected structures shows that the shear layers mostly overlap with a UMZ edge (in either plane).

We consider the stability of an elastic membrane on the bottom of a uniform horizontal flow of an inviscid and incompressible fluid of finite depth with free surface. The membrane is simply supported at the leading and the trailing edges which attach it to the two parts of the horizontal rigid floor. The membrane has an infinite span in the direction perpendicular to the direction of the flow and a finite length in the direction of the flow. For the membrane of infinite length we derive a full dispersion relation that is valid for arbitrary depth of the fluid layer and find conditions for the flutter of the membrane due to emission of surface gravity waves. We describe this radiation-induced instability by means of the perturbation theory of the roots of the dispersion relation and the concept of negative energy waves and discuss its relation to the anomalous Doppler effect.

We perform large-eddy simulations of turbulent flow in a channel constricted by streamwise periodically distributed hill-shaped protrusions. Two Reynolds number cases, i.e. $Re_h=10\,595$ (Fröhlich et al., J. Fluid Mech., vol. 526, 2005, pp. 19–66) and $Re_h=33\,000$ (Kähler et al., J. Fluid Mech., vol. 796, 2016, pp. 257–284), are repeated and utilized to verify and validate our numerical results, including the pressure and skin friction coefficients on bottom and top walls of the channel, mean velocity profiles and Reynolds stresses. All comparisons show reasonable agreement, providing a measure of validity that enables us to further probe simulation results at higher Reynolds number ($Re_h=10^5$) into aspects of flow physics that are not available from experiments. Effects of variation of Reynolds number are studied, with emphasis on the mean skin friction coefficients, separation bubble size and pressure fluctuations that are related to separation and reattachment. In addition, the main large-scale features of the separation behind the hill, including the scaling of the mean velocity profiles, are discussed. Furthermore, the instantaneous near-wall flow field is analysed in terms of skin friction portraits, and we confirm the existence of the local very small separation bubble on the hill crest as observed in experimental and numerical investigations. The flow field at the top wall, which is generally not given sufficient attention, is evaluated with the empirical friction law and universal logarithmic law as in planar channel flows. It is found that these empirical laws compare well with the large-eddy simulation results, although the hill constrictions behave as a perturbation source and the developed shear layer has some effects on the flow field near the top wall.

We present results on the effect of dispersed droplets in vertical natural convection (VC) using direct numerical simulations based on a two-way fully coupled Euler–Lagrange approach with a liquid phase and a dispersed droplets phase. For increasing thermal driving, characterised by the Rayleigh number, Ra, of the two analysed droplet volume fractions, $\alpha = 5\times 10^{-3}$ and $\alpha = 2\times 10^{-2}$, we find non-monotonic responses to the overall heat fluxes, characterised by the Nusselt number, Nu. The Nu number is larger when the droplets are thermally coupled to the liquid. However, Nu values remain close to the 1/4-laminar VC scaling, suggesting that the heat transport is still modulated by thermal boundary layers. Local analyses reveal the non-monotonic trends of local heat fluxes and wall-shear stresses: whilst regions of high heat fluxes are correlated to increased wall-shear stresses, the spatio-temporal distribution and magnitude of the increase are non-monotonic, implying that the overall heat transport is obscured by competing mechanisms. Most crucially, we find that the transport mechanisms inherently depend on the dominance of droplet driving to thermal driving that can quantified by (i) the bubblance parameter $b$, which measures the ratio of energy produced by the dispersed phase and the energy of the background turbulence, and (ii) ${\textit {Ra}}_d/{\textit {Ra}}$, where ${\textit {Ra}}_d$ is the droplet Rayleigh number, which we introduce in this paper. When $b \lesssim O(10^{-1})$ and ${\textit {Ra}}_d/{\textit {Ra}} \lesssim O(100)$, the Nu scaling is expected to recover to the VC scaling without droplets, and comparison with $b$ and ${\textit {Ra}}_d/{\textit {Ra}}$ from our data supports this notion.

Fluid deformable surfaces show a solid–fluid duality which establishes a tight interplay between tangential flow and surface deformation. We derive the governing equations as a thin film limit and provide a general numerical approach for their solution. The simulation results demonstrate the rich dynamics resulting from this interplay, where, in the presence of curvature, any shape change is accompanied by a tangential flow and, vice versa, the surface deforms due to tangential flow. However, they also show that the only possible stable stationary state in the considered setting is a sphere with zero velocity.

In the preparation of café latte, spectacular layer formation can occur between the espresso shot in a glass of milk and the milk itself. Xue et al. (Nat. Commun., vol. 8, 2017, pp. 1–6) showed that the injection velocity of espresso determines the depth of coffee–milk mixture. After a while, when a stable stratification forms in the mixture, the layering process can be modelled as a double diffusive convection system with a stably stratified coffee–milk mixture cooled from the side. More specifically, we perform (two-dimensional) direct numerical simulations of laterally cooled double diffusive convection for a wide parameter range, where the convective flow is driven by a lateral temperature gradient while stabilized by a vertical concentration gradient. Depending on the strength of stabilization as compared to the thermal driving, the system exhibits different flow regimes. When the thermal driving force dominates over the stabilizing force, the flow behaves like vertical convection in which a large-scale circulation develops. However, with increasing strength of the stabilizing force, a meta-stable layered regime emerges. Initially, several vertically-stacked convection rolls develop, and these well-mixed layers are separated by sharp interfaces with large concentration gradients. The initial thickness of these emerging layers can be estimated by balancing the work exerted by thermal driving and the required potential energy to bring fluid out of its equilibrium position in the stably stratified fluid. In the layered regime, we further observe successive layer merging, and eventually only a single convection roll remains. We elucidate the following merging mechanism: as weakened circulation leads to accumulation of hot fluid adjacent to the hot sidewall, larger buoyancy forces associated with hotter fluid eventually break the layer interface. Then two layers merge into a larger layer, and circulation establishes again within the merged structure.

We use experiments and direct numerical simulations to probe the phase space of low-curvature Taylor–Couette flow in the vicinity of the ultimate regime. The cylinder radius ratio is fixed at $\eta =r_i/r_o=0.91$, where $r_i \, (r_o)$ is the inner (outer) cylinder radius. Non-dimensional shear drivings (Taylor numbers ${\textit {Ta}}$) in the range $10^7\leq {\textit {Ta}}\leq 10^{11}$ are explored for both co- and counter-rotating configurations. In the ${\textit {Ta}}$ range $10^8\leq {\textit {Ta}}\leq 10^{10}$, we observe two local maxima of the angular momentum transport as a function of the cylinder rotation ratio, which can be described as either ‘co-’ or ‘counter-rotating’ due to their location or as ‘broad’ or ‘narrow’ due to their shape. We confirm that the broad peak is accompanied by the strengthening of the large-scale structures, and that the narrow peak appears once the driving (Ta) is strong enough. As first evidenced in numerical simulations by Brauckmann et al. (J. Fluid Mech., vol. 790, 2016, pp. 419–452), the broad peak is produced by centrifugal instabilities and that the narrow peak is a consequence of shear instabilities. We describe how the peaks change with ${\textit {Ta}}$ as the flow becomes more turbulent. Close to the transition to the ultimate regime when the boundary layers (BLs) become turbulent, the usual structure of counter-rotating Taylor vortex pairs breaks down and stable unpaired rolls appear locally. We attribute this state to changes in the underlying roll characteristics during the transition to the ultimate regime. Further changes in the flow structure around ${\textit {Ta}}\approx 10^{10}$ cause the broad peak to disappear completely and the narrow peak to move. This second transition is caused when the regions inside the BLs which are locally smooth regions disappear and the whole boundary layer becomes active.

The asymmetries that arise when a mixing layer involves two miscible fluids of differing densities are investigated using incompressible (low-speed) direct numerical simulations. The simulations are performed in the temporal configuration with very large domain sizes, to allow the mixing layers to reach prolonged states of fully turbulent self-similar growth. Imposing a mean density variation breaks the mean symmetry relative to the classical single-fluid temporal mixing layer problem. Unlike prior variable-density mixing layer simulations in which the streams are composed of the same fluids with dissimilar thermodynamic properties, the density variations are presently due to compositional differences between the fluid streams, leading to different mixing dynamics. Variable-density (non-Boussinesq) effects introduce strong asymmetries in the flow statistics that can be explained by the strongest turbulence increasingly migrating to the lighter fluid side as free-stream density difference increases. Interface thickness growth rates also reduce, with some thickness definitions particularly sensitive to the corresponding changes in alignment between density and streamwise velocity profiles. Additional asymmetries in the sense of statistical distributions of densities at a given position within the mixing layer reveal that fine scales of turbulence are preferentially sustained in lighter fluid, which also is where fastest mixing occurs. These effects influence statistics involving density fluctuations, which have important implications for mixing and more complicated phenomena that are sensitive to the mixing dynamics, such as combustion.

We explore the near-resonant interaction of inertial waves with geostrophic modes in rotating fluids via numerical and theoretical analysis. When a single inertial wave is imposed, we find that some geostrophic modes are unstable above a threshold value of the Rossby number $kRo$ based on the wavenumber and wave amplitude. We show this instability to be caused by triadic interaction involving two inertial waves and a geostrophic mode such that the sum of their eigenfrequencies is non-zero. We derive theoretical scalings for the growth rate of this near-resonant instability. The growth rate scaled by the global rotation rate is proportional to $(kRo)^2$ at low $kRo$ and transitions to a $kRo$ scaling for larger $kRo$. These scalings are in excellent agreement with direct numerical simulations. This instability could explain recent experimental observations of geostrophic instability driven by waves.

We investigate the formation of subaqueous transverse bedforms in turbulent open channel flow by means of direct numerical simulations with fully resolved particles. The main goal of the present analysis is to address the question whether the initial pattern wavelength scales with the particle diameter or with the mean fluid height. A previous study (Kidanemariam & Uhlmann, J. Fluid Mech., vol. 818, 2017, pp. 716–743) has observed a lower bound for the most unstable pattern wavelength in the range 75–100 times the particle diameter, which was equivalent to 3–4 times the mean fluid height. In the current paper, we vary the streamwise box length in terms of the particle diameter and of the mean fluid height independently in order to distinguish between the two possible scaling relations. For the chosen parameter range, the obtained results clearly exhibit a scaling of the initial pattern wavelength with the particle diameter, with a lower bound around a streamwise extent of approximately $80$ particle diameters. In longer domains, on the other hand, patterns are observed at initial wavelengths in the range 150–180 times the particle diameter, which is in good agreement with experimental measurements. Variations of the mean fluid height, on the other hand, seem to have no significant influence on the most unstable initial pattern wavelength. We argue that the observed scaling with the particle diameter is due to the wake effect induced by the seeds which are formed by initially dislodged lumps of particles, in accordance with the ideas of Coleman & Nikora (Water Resour. Res., vol. 45, 2009, W04402). Finally, for the cases with the largest relative submergence, we observe spanwise and streamwise sediment waves of similar amplitude to evolve and superimpose, leading to three-dimensional sediment patterns.

Most of our understanding of moving contact lines relies on the limit of small capillary number ${Ca}$. This means the contact line speed is small compared to the capillary speed $\gamma /\eta$, where $\gamma$ is the surface tension and $\eta$ the viscosity, so that the interface is only weakly curved. The majority of recent analytical work has assumed in addition that the angle between the free surface and the substrate is also small, so that lubrication theory can be used. Here, we calculate the shape of the interface near a slip surface for arbitrary angles, and for two phases of arbitrary viscosities, thereby removing a key restriction in being able to apply small capillary number theory. Comparing with full numerical simulations of the viscous flow equation, we show that the resulting theory provides an accurate description up to $Ca \approx 0.1$ in the dip coating geometry, and a major improvement over theories proposed previously.