The exhaust from combustion engines contains particulate matter (PM), which poses potential health risks to human lungs. Current emission laws place increasingly strict limitations on both PM and particle number, leading to the necessity of using wall-flow filters to separate out a significant amount of the introduced PM. As this leads to an increase in the filter's loading, it is regenerated continuously or periodically, leading to the rearrangement of individual particulate structures inside the filter channels. Such rearrangement events cause the formation of specific deposition patterns, which affect the filter's pressure drop, its loading capacity and the separation efficiency. In order to derive predictions on the formation of specific deposition patterns, the transient behaviour of individual particle structures needs to be examined. The present work investigates the detachment and transport of particle structures during filter regeneration with three-dimensional surface-resolved simulations using a lattice Boltzmann method. The goal of this work is the determination of relevant key quantities and their interpretation with respect to predictions regarding the resulting deposition patterns. In this context, it is shown that lift forces are not the predominant detachment forces for non-spherical particle structures, and that the stopping distance of such structures is too long to avoid back-end deposition.

We numerically study the melting process of a solid layer heated from below such that a liquid melt layer develops underneath. The objective is to quantitatively describe and understand the emerging topography of the structures (characterized by the amplitude and wavelength), which evolve out of an initially smooth surface. By performing both two-dimensional (achieving Rayleigh number up to $Ra=10^{11}$) and three-dimensional (achieving Rayleigh number up to $Ra=10^9$) direct numerical simulations with an advanced finite difference solver coupled to the phase-field method, we show how the interface roughness is spontaneously generated by thermal convection. With increasing height of the melt the convective flow intensifies and eventually even becomes turbulent. As a consequence, the interface becomes rougher but is still coupled to the flow topology. The emerging structure of the interface coincides with the regions of rising hot plumes and descending cold plumes. We find that the roughness amplitude scales with the mean height of the liquid layer. We derive this scaling relation from the Stefan boundary condition and relate it to the non-uniform distribution of heat flux at the interface. For two-dimensional cases, we further quantify the horizontal length scale of the morphology, based on the theoretical upper and lower bounds given for the size of convective cells known from Wang et al. (Phys. Rev. Lett., vol. 125, 2020, 074501). These bounds agree with our simulation results. Our findings highlight the key connection between the morphology of the melting solid and the convective flow structures in the melt beneath.

Understanding and predicting the dynamics of dispersed micro-objects in microfluidics is crucial in numerous natural, industrial and technological situations. In this paper, we experimentally characterized the equilibrium velocity $V$ and lateral position $\varepsilon$ of various dispersed micro-objects, such as beads, bubbles and drops, in a cylindrical microchannel over an unprecedentedly wide range of parameters. By varying the dimensionless object size ($d \in [0.1; 1]$), the viscosity ratio ($\lambda \in [10^{-2}; \infty [$), the density ratio ($\varphi \in [10^{-3}; 2]$), the Reynolds number ($Re \in [10^{-2}; 10^2]$) and the capillary number ($Ca \in [10^{-3}; 0.3]$), we offer an exhaustive parametric study exploring various dynamics from the non-deformable viscous regime to the deformable inertial regime, thus enabling us to highlight the sole and combined roles of inertia and capillary effects on lateral migration. Experiments are compared and agree well with a steady three-dimensional Navier–Stokes model for incompressible two-phase fluids, including the effects of inertia and possible interfacial deformations. This model enables us to propose a correlation for the object velocity $V$ as functions of $d$, $\varepsilon$ and $\lambda$, obtained in the ${Re}={Ca}=0$ limit, but valid for a larger range of values of ${Re}$ and ${Ca}$ delimited by the validity of the linear regime. Next, we present stability maps for the centred position showing that non-deformable objects dominated by inertial effects are only stable if large enough, typically for $d \gtrsim 0.7$, whereas deformable objects dominated by capillary effects can be stable for much smaller sizes, provided the viscosity ratio is outside the range $0.7 \lesssim \lambda \lesssim 10$, in which deformability also plays a destabilizing effect, as for inertia.

We experimentally and numerically study the dynamics of a liquid jet issued from a rotating orifice, whose breakup is regulated by a vibrating piezo element. The helical trajectory of the spiralling jet yields fictitious forces varying along the jet whose longitudinal projections stretch and thin the jet, affecting the growth of perturbations. We show that by quantifying these fictitious forces, one can estimate the jet intact length and size distribution of drops formed at jet breakup. The presence of the locally varying fictitious forces may render high-frequency perturbations, that would otherwise be stable in the abscence of stretching, unstable, as observed similarly in the case of straight jets stretching under gravity. The perturbation amplitude then dictates how strong the perturbation is coupled to the jet compared with random noise that is inherently present in any experimental set-up. In the present study we exploit the slenderness of the jet to separate the calculation of the base flow and the growth of perturbations. The fictitious forces calculated from the base flow trajectory are then used in a nonlinear slender-jet model, which treats the spiralling jet as a quasi-straight jet with locally varying body forces. We show both experimentally and numerically that jet breakup characteristics (e.g. intact length and drop size distribution) can be controlled by finite-amplitude perturbations created by mechanically induced pressure modulations. Finally, we revisit the integrated net gain approach developed for straight jets under gravity and we provide simple analogous relations for spiralling jets.

We experimentally investigate the evaporation of very volatile liquid droplets (Novec 7000 Engineered Fluid, chemical name hydrofluoroethers HFE-7000) in a turbulent spray. Droplets with diameters of the order of a few micrometres are produced by a spray nozzle and then injected into a purpose-built enclosed dodecahedral chamber filled with air containing various amounts of water vapour. The ambient temperature and relative humidity in the chamber are carefully controlled. We observe water condensation on the rapidly evaporating droplet, both for the spray and for a single acoustically levitated millimetric Novec 7000 droplet. We further examine the effect of humidity, and reveal that a more humid environment leads to faster evaporation of the volatile liquid, as well as more water condensation. This is explained by the much larger latent heat of water as compared with that of Novec 7000. We extend an analytical model based on Fick's law to quantitatively account for the data.

Well-resolved direct numerical simulations (DNS) have been performed of the flow in a smooth circular pipe of radius $R$ and axial length $10{\rm \pi} R$ at friction Reynolds numbers up to $Re_\tau =5200$ using the pseudo-spectral code OPENPIPEFLOW. Various turbulence statistics are documented and compared with other DNS and experimental data in pipes as well as channels. Small but distinct differences between various datasets are identified. The friction factor $\lambda$ overshoots by $2\,\%$ and undershoots by $0.6\,\%$ the Prandtl friction law at low and high $Re$ ranges, respectively. In addition, $\lambda$ in our results is slightly higher than in Pirozzoli et al. (J. Fluid Mech., vol. 926, 2021, A28), but matches well the experiments in Furuichi et al. (Phys. Fluids, vol. 27, issue 9, 2015, 095108). The log-law indicator function, which is nearly indistinguishable between pipe and channel up to $y^+=250$, has not yet developed a plateau farther away from the wall in the pipes even for the $Re_\tau =5200$ cases. The wall shear stress fluctuations and the inner peak of the axial turbulence intensity – which grow monotonically with $Re_\tau$ – are lower in the pipe than in the channel, but the difference decreases with increasing $Re_\tau$. While the wall value is slightly lower in the channel than in the pipe at the same $Re_\tau$, the inner peak of the pressure fluctuation shows negligible differences between them. The Reynolds number scaling of all these quantities agrees with both the logarithmic and defect-power laws if the coefficients are properly chosen. The one-dimensional spectrum of the axial velocity fluctuation exhibits a $k^{-1}$ dependence at an intermediate distance from the wall – also seen in the channel. In summary, these high-fidelity data enable us to provide better insights into the flow physics in the pipes as well as the similarity/difference among different types of wall turbulence.

The aeroacoustic characteristics of flying vehicles with pitch-fixed rotors differ from traditional helicopters with pitch-controlled rotor blades. Accurate predictions of rotor noise are still challenging because many uncertainty factors and unsteadinesses exist. This work investigates the aeroacoustic effects of rotational speed deviation, rotation speed fluctuation, blade vibration and blade geometric asymmetry. The analysis is based on the efficient computation of rotor noise under different working conditions. The mean aerodynamic variables are computed using the blade element moment theory, while small-amplitude fluctuations are introduced to account for the unsteadiness and uncertainty factors. It is shown that periodic rotation speed fluctuations and blade vibrations can produce significant extra tones. By contrast, if the fluctuations and vibrations are random, the noise level in a wide frequency range is increased. The intriguing result reminds us of the need to revisit the rotor broadband noise sources commonly attributed to turbulent flows. The influences are observer angle dependent, and the extra noise production is more significant in the upstream and downstream directions. The asymmetric blade geometry can cause extra tonal noise at the harmonics of the blade shaft frequency. The noise features of dual rotors are also investigated. Usually, the noise is sensitive to the initial phase difference and rotation directions due to the interference effect. However, the noise features are vastly altered if there are slight differences in the rotation speeds. Although the influences of some factors on rotor noise were already known, the present study provides a more comprehensive analysis of the problem. The results also highlight the need to consider these practical factors for accurate noise prediction of multi-rotor flying vehicles.

Immiscible and incompressible liquid–liquid flows are considered in a Taylor–Couette geometry and analysed by direct numerical simulations coupled with the volume-of-fluid method and a continuum surface force model. The system Reynolds number $Re \equiv r_i \omega _i d / \nu$ is fixed to $960$, where the single-phase flow is in the steady Taylor vortex regime, whereas the secondary-phase volume fraction $\varphi$ and the system Weber number $We \equiv \rho r_i^2 \omega _i^2 d / \sigma$ are varied to study the interactions between the interface and the Taylor vortices. We show that different Weber numbers lead to two distinctive flow regimes, namely an advection-dominated regime and an interface-dominated regime. When $We$ is high, the interface is easily deformed because of its low surface tension. The flow patterns are then similar to the single-phase flow, and the system is dominated mainly by advection (advection-dominated regime). However, when $We$ is low, the surface tension is so large that stable interfacial structures with sizes comparable to the cylinder gap can exist. The background velocity field is modulated largely by these persistent structures, thus the overall flow dynamics is governed by the interface (interface-dominated regime). The effect of the interface on the global system response is assessed by evaluating the Nusselt number $Nu_{\omega }$ based on the non-dimensional angular velocity transport. It shows non-monotonic trends as functions of the volume fraction $\varphi$ for both low and high $We$. We explain how these dependencies are closely linked to the velocity and interfacial structures.

The dynamics of an evaporating droplet in an unsteady flow is of practical interest in many industrial applications and natural processes. To investigate the transport and evaporation dynamics of such droplets, we present a numerical study of an isolated droplet in an oscillating gas-phase flow. The study uses a one-way coupled two-phase flow model to assess the effect of the amplitude and the frequency of a sinusoidal external flow field on the lifetime of a multicomponent droplet containing a non-volatile solute dissolved in a volatile solvent. The results show that the evaporation process becomes faster with an increase in the amplitude or the frequency of the gas-phase oscillation. The liquid-phase transport inside the droplet also is influenced by the unsteadiness of the external gas-phase flow. A scaling analysis based on the response of the droplet under the oscillating drag force is subsequently carried out to unify the observed evaporation dynamics in the simulations under various conditions. The analysis quantifies the enhancement in the droplet velocity and Reynolds number as a function of the gas-phase oscillation parameters and predicts the effects on the evaporation rate.

A mathematical model is developed to investigate seabed heat transfer processes under long-crested ocean waves. The unsteady convection–diffusion equation for water temperature includes terms depending on the velocity field in the laminar boundary layer, analogous to mass transfer near the seabed. Here we consider regular progressive waves and standing waves reflected from a vertical structure, which complicate the convective term in the governing equation. Rectangular and Gaussian distributions of seabed temperature and heat flux are considered. Approximate analytical solutions are derived for uniform and trapezoidal currents, and compared against predictions from a numerical solver of the full equations. The effects of heat source profile, location and strength on heat transfer dynamics in the thermal boundary layer are explained, providing insights into seabed temperature forced convection mechanisms enhanced by free-surface waves.

After decades of research efforts, wind–wave interaction mechanisms have been recognized as extremely elusive. The reason is the complex nature of the problem, which combines complex coupling mechanisms between turbulent wind and water waves with the presence of multiple governing parameters, such as the friction Reynolds number of the wind, the water depth and the wind fetch. As shown unequivocally here, the use of suitable flow settings allows us to reduce the complex problem of wind–wave interaction to its essential features, mainly as a function of the sole friction Reynolds number of the wind. The resulting numerical solution allows us to study the interactions between water and air layers with their own fluid properties, and to unveil very interesting features, such as an oblique wave pattern travelling upstream and a wave-induced Stokes sublayer. The latter is responsible for a drag reduction mechanism in the turbulent wind. Despite the simulated flow conditions being far from the intense events occurring at the ocean–atmosphere interface, the basic flow phenomena unveiled here may explain some experimental evidence in wind–wave problems. Among other things, the wave-induced Stokes sublayer may shed light on the large scatter of the drag coefficient data in field measurements where swell waves of arbitrary directions are often present. Hence the present results and the developed approach pave the way for the understanding and modelling of the surface fluxes at the ocean–atmosphere interface, which are of overwhelming importance for climate science.

The propagation of waves from a vertical uplift of a slender rectangular fault in a sea of constant depth is discussed, accounting for water compressibility, gravity and seabed elasticity. The compressed water column results in the generation of acoustic–gravity waves that travel at the speed of sound in water. Acoustic–gravity waves are found to terminate after a finite time, with the decay time most influenced by seabed rigidity, which is in contrast to the rigid stationary-phase model where signals persist indefinitely. At certain frequencies acoustic–gravity waves couple with the elastic seabed and travel at the shear velocity (speed of sound in an elastic solid). Improved estimates of the critical frequencies are derived. Moreover, besides the usual tsunami, a second – very small amplitude – surface wave mode travelling at the speed of sound arises under certain frequencies. We derive the cut-off frequency for this mode. The acoustic modes possess a frequency spectrum which depends on the time evolution and spatial properties of the rupture. We find that appropriate filtering of the acoustic–gravity wave signal can reveal characteristic peaks that encode information on the fault's geometry and dynamics.

This paper presents an analytic model for the analysis of co-planar turbine fences that partially span the width of a channel in which the flow is driven by a sinusoidally oscillating driving head. The thrust presented by the turbines reduces the flow rate through the channel leading to a solution for overall power that is dependent upon turbine resistance and flow blockage as well as on channel characteristics. We introduce a return parameter, in terms of power per turbine area, to assess optimum turbine fence deployment for a given channel. We find that the optimal deployment rests on a universal curve independent of the channel characteristics, and that these characteristics – namely the integrated channel bed friction and a modified channel Froude number – move the optimum along this curve. We find that blockage considerations play a large role in the performance of a tidal farm – its achievable power, optimal return, channel flow rate reduction and device thrust – and that the scales of blockage must be considered even when designing relatively unblocked farms. The impact of the channel characteristics on the optimal arrangement, alongside environmental constraints that may limit permissible flow blockage, are quantified and discussed.

Nonlinearity of density stratification modulates buoyancy effects. We report results from a body-inclusive large eddy simulation of a wake in nonlinear stratification, specifically for a circular disk at diameter-based Reynolds number (${\it Re} $) of $5000$. Five density profiles are considered; the benchmark has linear stratification and the other four have hyperbolic tangent profiles of the same thickness to model a pycnocline. The disk moves inside the central core of the pycnocline in two of those four cases and, in the other two cases with a shifted density profile, the disk moves partially/completely outside the pycnocline. The maximum buoyancy frequency ($N_{max}$) for all the profiles is the same. The first part of the study investigates the centred cases. Non-uniform stratification results in increasing wake turbulence relative to the benchmark owing to reduced suppression of turbulence production as well as wave trapping in the pycnocline. Steady lee waves are also quantified to understand the limitations of linear theory. The second part pays attention to the effect of a relative shift between the pycnocline and the disk. The wake defect velocity decays substantially faster in the cases with a shift and the wake has higher turbulence level. The effect of disk location on the Kelvin wake waves (a family of steady waves within the pycnocline) and its modal form is obtained and explained by solving the Taylor–Goldstein equation. The family of unsteady internal gravity waves that are generated by the wake is also studied and the effect of disk shift is quantified.

Drop impact onto a thin liquid film of another liquid is observed and characterized using a high-speed video system. A new mode of splash – a complete, simultaneous corona detachment – has been observed, which is the result of the lamella breakup near the wall film. The abrupt outward and upward displacement of the lamella leads to an extreme stretching of the corona wall, resulting in its rapid thinning and a rupture. This rupture triggers propagating Taylor–Culick rims, which rapidly spread, meet and, thus, undercut simultaneously the entire corona, resulting in its detachment. Special experiments with the spreading corona impingement onto a fixed needle, supplement the physical evidence of the above-mentioned mechanism. A self-consistent theory of the observed phenomena is proposed and compared with experiments, exhibiting good agreement.

The present study investigates the profiles of statistically axisymmetric turbulent jets with arbitrary buoyancy. Analytical expressions for the shape of the radial velocity, Reynolds stress and radial scalar flux profiles are derived from the governing equations by assuming self-similar Gaussian mean velocity and scalar profiles. Previously these have only been derived for the special cases of pure jets and plumes, whereas the present study generalises them to arbitrary buoyancies. These are then used to derive analytical expressions for the turbulent Schmidt/Prandtl numbers, which, along with the mean profiles, are shown to give predictions in agreement with existing literature.

At high surface Péclet numbers, it is common to associate the presence of surfactants with surface immobilization, where a free surface becomes indistinguishable from a no-slip surface. A different mechanism has recently been proposed for longitudinal shear flow along a unidirectional trench (Baier & Hardt, J. Fluid Mech., vol. 949, 2022, A34) wherein, at high Marangoni numbers, the meniscus spanning the finite-length trench becomes a constant-shear-stress surface due to contamination by incompressible surfactant. That model predicts recirculating interfacial flows on the meniscus, a phenomenon that has been observed experimentally (Song et al., Phys. Rev. Fluids, vol. 3, issue 3, 2018, 033303). By finding an explicit solution to the constant-shear-stress model at all protrusion angles and calculating the effective slip length for a dilute mattress of such surfactant-laden trenches, we show that those effective slip lengths are almost indistinguishable from those for a surface whose menisci have the same deflection but have been completely immobilized (i.e. they are no-slip surfaces). This means that, despite the presence of non-trivial recirculating vortical flows on the menisci, the aggregate slip characteristics of such surfaces are that they have been effectively immobilized. This surprising result underscores the need for caution in comparing theory with experiments based on effective slip properties alone.

Direct numerical simulations were performed to characterize fully developed supersonic turbulent channel flows over isothermal rough walls. The effect of roughness was incorporated using a level-set/volume-of-fluid immersed boundary method. Turbulence statistics of five channel flows are compared, including one reference case with both walls smooth and four cases with smooth top walls and bottom walls with two-dimensional (2-D) and three-dimensional (3-D) sinusoidal roughnesses. Results reveal a strong dependence of the turbulence on the roughness topography and the associated shock patterns. Specifically, the 2-D geometries generate strong oblique shock waves that propagate across the channel and are reflected back to the rough-wall side. These strong shocks are absent in the smooth-wall channel and are significantly weaker in cases with 3-D roughness geometries, replaced by weak shocklets. At the impingement locations of the shocks on the top wall in the 2-D roughness cases, localized augmentations of turbulence shear production are observed. Such regions of augmented production also exist for the 3-D cases, at a much weaker level. The oblique shock waves are thought to be responsible for a more significant entropy generation for cases with 2-D surfaces than those with 3-D ones, leading to a higher irreversible heat generation and consequently higher temperature values in 2-D roughness cases. In the present supersonic channels, the effects of roughness extend beyond the near-wall layer due to the shocks. This suggests that outer layer similarity may not fully apply to a rough-wall supersonic turbulent flow.

We present the analytical solution for the two-dimensional velocity and density fields within an approximation for laminar stratified inclined duct (SID) flows where diffusion dominates over inertia in the along-channel momentum equation but is negligible in the density transport equation. We refer to this approximation as the hydrostatic/gravitational/viscous in momentum and advective in density (HGV-A) approximation due to the leading balances in the governing equations. The analytical solution is valid for laminar flows in a two-layer configuration in the limit of long ducts. The non-dimensional volume flux within the HGV-A approximation is given by $Fr^* ={{Re}}_g/(AK)$, which is a control parameter with ${{Re}}_g$ the gravitational Reynolds number, $A$ the aspect ratio of the duct and $K$ a geometrical parameter that depends on the tilt of the duct and is obtained from the analytical solution. This analytical solution was validated against results from laboratory experiments, and allows us to gain new insight into the dynamics and properties of SID flows. Most importantly, constant values of $Fr^*$ describe, in both horizontal and inclined ducts, the transitions between increasingly turbulent flow regimes: from laminar flow, to interfacial waves, to intermittent turbulence and sustained turbulence.