We consider direct statistical simulation (DSS) of a paradigm system of convection interacting with mean flows. In the Busse annulus model, zonal jets are generated through the interaction of convectively driven turbulence and rotation; non-trivial dynamics including the emergence of multiple jets and bursting ‘predator–prey’ type dynamics can be found. We formulate the DSS by expanding around the mean flow in terms of equal-time cumulants and arrive at a closed set of equations of motion for the cumulants. Here, we present results using an expansion terminated at the second cumulant (CE2); it is fundamentally a quasilinear theory. We focus on particular cases including bursting and bistable multiple jets and demonstrate that CE2 can reproduce the results of direct numerical simulation if particular attention is given to symmetry considerations.

The unsteady dispersion of a solute band by a steady pressure-driven flow in a semicircular microchannel is theoretically studied via the generalized dispersion model. Considering an irreversible first-order reaction at the curved wall while assuming a no-flux boundary condition at the flat wall, analytical solutions are obtained for the exchange, convection and dispersion coefficients as well as the dimensionless forms of solute concentration and mean solute concentration. The solutions are obtained assuming an initial solute band of arbitrary cross-sectional shape and axial distribution and the results are presented for both circular and semicircular shapes with the uniform distribution being a special case of the latter. Besides the general solutions, special solutions are also derived for uniform velocity and no-reaction cases. In the following, the influences of the initial concentration distribution and the Damköhler number, a measure of the reaction rate, on the transport coefficients and the concentration distribution are investigated in depth. It is demonstrated that the combination of the initial concentration distribution and the Damköhler number specifies the variations of the transport coefficients in the short term but the Damköhler number is the only parameter dominating the long-term values: the exchange and convection coefficients are increasing functions of the Damköhler number whereas the opposite is true for the dispersion coefficient. Moreover, we show that, provided the solute injection is appropriately positioned and shaped, liquid-phase transportation with little dispersion is possible in typical microchannels utilizing semicircular geometry.

The flow over a forward-facing step (FFS) at $Ma_\infty =1.7$ and $Re_{\delta _0}=1.3718\times 10^{4}$ is investigated by well-resolved large-eddy simulation. To investigate effects of upstream flow structures and turbulence on the low-frequency dynamics of the shock wave/boundary layer interaction (SWBLI), two cases are considered: one with a laminar inflow and one with a turbulent inflow. The laminar inflow case shows signs of a rapid transition to turbulence upstream of the step, as inferred from the streamwise variation of $\langle C_f \rangle$ and the evolution of the coherent vortical structures. Nevertheless, the separation length is more than twice as large for the laminar inflow case, and the coalescence of compression waves into a separation shock is observed only for the fully turbulent inflow case. The dynamics at low and medium frequencies is characterized by a spectral analysis, where the lower frequency range is related to the unsteady separation region, and the intermediate one is associated with the shedding of shear layer vortices. For the turbulent inflow case, we furthermore use a three-dimensional dynamic mode decomposition to analyse the individual contributions of selected modes to the unsteadiness of the SWBLI. The separation shock and Görtler-like vortices, which are induced by the centrifugal forces in the separation region, are strongly correlated with the low-frequency unsteadiness in the current FFS case. Similarly as observed previously for the backward-facing steps, we observe a slightly higher non-dimensional frequency (based on the separation length) of the low-frequency mode than for SWBLI in flat plate and ramp configurations.

Shock-tube experiments on various two-bubble and two-spike interfaces are performed to examine the dependence of bubble competition and spike competition on the initial spectra and density ratio of the interface. The differences in the influences of bubble competition and spike competition on the Richtmyer–Meshkov instability are highlighted for the first time. The bubble-competition effect is mainly dependent on the initial spectra of the two-bubble configuration. In contrast, the spike-competition effect is determined by both the initial spectra and density ratio. The extended buoyancy–drag model is adopted to explain the variation of the drag force imposed on the long-wavelength and short-wavelength structures as the initial conditions change. Based on the spectrum analysis, it is found that the constituent modes of two-bubble and two-spike interfaces have different contributions to the long-wavelength and short-wavelength perturbation growths. A generalised, nonlinear, analytical model is then established to quantify the bubble-competition effect and spike-competition effect considering arbitrary initial spectra and density ratio. The bubble-competition effect is believed to be stronger than the spike-competition effect at a high density ratio because it suppresses the high-frequency perturbation growth more evidently.

Oblique breakdown in a Mach 2.0 supersonic boundary layer controlled by a local cooling strip with a temperature jump is investigated using direct numerical simulations and linear stability theory. The effect of temperature on the stability of the fundamental oblique waves is first studied by linear stability theory. It is shown that the growth rate of fundamental oblique waves will decrease monotonically as the temperature decreases. However, the results of the direct numerical simulations indicate that transition reversal will occur as the growth rate of the fundamental oblique waves of cooled case becomes faster compared with that of baseline case downstream of the cooling strip. When the cooling strip is in the linear region, the transition is delayed due to the suppression effect of the cooling strip on the fundamental oblique waves. When the cooling strip is located in the early nonlinear region, the fundamental oblique waves will be suppressed by higher spanwise wavenumber steady modes generated by the mutual and self-interaction between the fundamental oblique waves and harmonic modes, which is first called the self-suppression effect (SSE) in the present study. Further research indicated that the meanflow distortion generated by steady modes plays an important role in the SSE. Compared with the stabilization effect of the cooling strip, the SSE is more effective. Moreover, the SSE might provide a new idea on the instability control, as it is observed that the SSE works three times leading to the growth rate of fundamental oblique waves slowing down at three different regions, respectively.

The linear stability of the axisymmetric steady thermocapillary flow in a liquid bridge made from 2 cSt silicone oil (Prandtl number 28) is investigated numerically in the framework of the Boussinesq approximation. The flow and temperature fields in the surrounding gas phase (air) are taken into account for a generic cylindrical container hosting the liquid bridge. The flows in the liquid and in the gas are fully coupled across the hydrostatically deformed liquid–gas interface, neglecting dynamic interface deformations. Originating from a common reference case, the linear stability boundary is computed varying the length of the liquid bridge (aspect ratio), its volume and the gravity level, providing accurate critical data. The qualitative dependence of the critical threshold on these parameters is explained in terms of the characteristics of the critical mode. The heat exchange between the ambient gas and the liquid bridge that is fully resolved has an important influence on the critical conditions.

The unsteady effects of transitioning surface tension, $\gamma (t)$, on the dynamics of capillary imbibition in channels of arbitrary shape are analytically investigated with a focus on rectangular and circular channels. With proper scaling, two unsteady models for $\gamma (t)$ are defined and used to highlight this transient behaviour. The convoluted dynamics at the flow front are correctly captured in the governing equations, which are rigorously analysed using unsteady eigenfunction expansion. Then, the final solution and data are obtained by employing the Runge–Kutta fourth-order scheme elegantly applied simultaneously to two derived nonlinear ordinary differential equations. Ultimately, the results are more accurate compared with previous studies. Dynamics and kinematic similarity between rectangular and circular channels are investigated and discussed and the conditions for equivalence in both channels are highlighted. Using a small parameter ($\epsilon$) that stretches the time scale, we successfully use a robust asymptotic analysis to develop and capture the long-time dynamics. Ultimately, we recover the Lucas–Washburn regime analysed in Washburn (Phys. Rev., vol. 17, 1921, pp. 273–283), Lucas (Kolloidn. Z., vol. 23, 1918, pp. 15–22) for steady surface tension where the variations of depth and rate with time result in $h\thicksim t^{1/2}$ and $v \thicksim t^{-1/2}$. In the end, the three forces, namely the inertia, $F_{v}$, the viscous, $F_{\mu }$, and the surface tension, $F_{\gamma }$, are briefly analysed and used to highlight three main distinct regimes. We show that at early times, $F_{v}/F_{\gamma } \thicksim 1$, whereas at a long time, $F_{\mu }/F_{\gamma } \thicksim -1$.

In this study, we develop an analytical model to predict the turbulent boundary layer downstream of a step-change in the surface roughness where upstream flow conditions are given. We first revisit the classical model of Elliott (Trans. Am. Geophys. Union, vol. 39, 1958, pp. 1048–1054), who modelled the velocity distribution within and above the internal layer with a simple piecewise logarithmic profile, and evolved the velocity profile using the streamwise momentum equation. Elliott's model was originally developed for an atmospheric surface layer, and to make the model applicable to a spatially developing turbulent boundary layer with finite thickness, we propose a number of more physical refinements, including adding a wake function to the velocity profile, considering the growth of the entire boundary layer in the streamwise direction, and using a more realistic shear stress profile in the momentum equation. In particular, we implement the blending model (Li et al., J. Fluid Mech., vol. 923, 2021, p. A18) to account for the deviation of the mean flow within the internal layer from a canonical velocity profile based on the local wall condition. These refinements lead to improved agreement between the prediction and the measurement, especially in the vicinity of the rough-to-smooth change.

We numerically investigated the opposing-buoyancy mixed convection through and around square arrays of $10\times 10$ heated circular cylinders with the solid fraction ($\phi$) ranging from 0.0079 to 0.66 and the Richardson number ($Ri$) varying from 0 to 1 at a fixed Reynolds number ($Re$) of 100. Our simulations revealed that the large mean recirculation in the far wake can be detached from or connected with the vortex pair in the near wake for different combinations of $Ri$ and $\phi$. Also, it was found that the array with relatively small $\phi$ can significantly promote flow instability even at moderate $Ri$. The instability, which is closely related to the fluctuating heat flux, develops from the lateral sides to the downstream side of the array and gives rise to the large mean recirculation in the far wake. The power spectra density of the array-scale force coefficients demonstrates that the flow undergoes different bifurcation behaviours under various parameter combinations, which reflects the interaction between the near-wake and far-wake vortexes. Interestingly, the Strouhal–Richardson number curves can be collapsed onto the same curve when $Ri$ is increased by a $\phi$-dependent factor. Also, for $\phi \leqslant 0.22$, both the mean drag coefficient and the mean Nusselt number of the array were found to decrease linearly with $Ri$ since the buoyancy within the array becomes prominent in this range of $\phi$.

We investigate the transient dynamics of a sheared suspension of neutrally buoyant particles in pressure-imposed rheology configuration, subject to a sudden change in shear rate or external pressure. Discrete element method simulations show that, depending on the flow parameters (particle and system size, initial volume fraction), the early stress response of the suspension may differ strongly from the prediction of the suspension balance model based on the steady-state rheology. We show that a two-phase model incorporating the Reynolds-like dilatancy law of Pailha & Pouliquen (J. Fluid Mech., vol. 633, 2009, pp. 115–135), which prescribes the dilation rate of the suspension over a strain scale $\gamma _0$, captures quantitatively the suspension dilation/compaction over the whole range of parameters investigated. Together with the Darcy flow induced by the pore pressure gradient during dilation or compaction, this Reynolds-like dilatancy implies that the early stress response of the suspension is non-local, with a non-local length scale $\ell$ that scales with the particle size and diverges algebraically at jamming. In regions affected by $\ell$, the stress level is fixed, not by the steady-state rheology, but by the Darcy fluid pressure gradient resulting from the dilation/compaction rate. Our results extend the validity of the Reynolds-like dilatancy flow rule, initially proposed for jammed suspensions, to flowing suspension below the critical volume fraction at which the suspension jams, thereby providing a unified framework to describe dilation and shear-induced migration. They pave the way for understanding more complex unsteady flows of dense suspensions, such as impacts, transient avalanches or the impulsive response of shear-thickening suspensions.

Numerical simulations have been conducted to identify the dominant mechanism responsible for driving secondary flow motions in horizontal particle-laden pipe flows, based on an analysis of the forces acting on each phase. A four-way coupling Euler–Lagrangian approach was employed, using direct numerical simulations for the gas phase and Lagrangian particle tracking to account for the drag, gravitational and lift forces, together with the interactions that occur for both particle–wall and inter-particle collisions. The four different flow regimes, which had been identified previously as depending on various combinations of flow parameters and are characterised by the secondary flow structures of both the fluid and particle phases, were identified via varying the mass loading alone from $\varPhi _m=0.4$ to $\varPhi _m=1.8$. The distribution of the divergence of Reynolds stresses was used to help characterise the classes of the secondary fluid flow. This shows that secondary fluid flows of both the first and second kinds can either exist separately or co-exist in such flows. The forces exerted on the fluid phase by the pressure gradient and fluid–particle interactions were examined qualitatively and quantitatively to identify their contribution to the secondary fluid flow motions. A similar study was also applied to the drag, lift and gravitational forces exerted on the particle phase for the secondary particle flow motions. These were found to explain the secondary flows of both the fluid and particle phases with regard to both the flow direction and magnitude, together with the interaction between the two phases.

We perform a three-dimensional direct numerical simulation of flow over the Tacoma Narrows Bridge to understand the vertical and torsional vibrations that occurred before its collapse in 1940. Real-scale structural parameters of the bridge are used for the simulation. The Reynolds number based on the free-stream velocity and height of the deck fence is lower (${Re}=10\ 000$) than the actual one on the day of its collapse (${Re}=3.06 \times 10^{6}$), but the magnitude of a fluid property is modified to provide the real-scale aerodynamic force and moment on the deck. The vertical and torsional vibrations are simulated through two-way coupling of the fluid flow and structural motion. The vertical vibration occurs from the frequency lock-in with the vortex shedding, and its wavelength and frequency agree well with the recorded data in 1940. After saturation of the vertical vibration, a torsional vibration resulting from the aeroelastic fluttering grows exponentially in time, with its wavelength and frequency in excellent agreement with the recorded data of the incident. The critical flutter wind speed for the growth of torsional vibration is obtained as $3.56 < U_c / (f_{nat} B) \le 4$, where $U_{c}$ is the critical flutter wind speed, $f_{nat}$ is the natural frequency of the torsional vibration and $B$ is the deck width. Finally, apart from the actual vibration process in 1940, we perform more numerical simulations to investigate the roles of the free-stream velocity and vertical vibration in the growth of the torsional vibration.

Global stability analysis and direct numerical simulation (DNS) are used to study boundary layer flows with an isolated roughness element. The aspect ratio of the element ($\eta$) is small, while the ratio of element height to displacement boundary layer thickness ($h/\delta ^{*}$) is large. Both steady base flows and time-averaged mean flows are able to capture the frequencies of the primary vortical structures and mode shapes. Global stability results highlight that although the varicose instability is dominant for large $h/\delta ^{*}$, sinuous instability becomes more pronounced as $Re_h$ increases for the thin geometry ($\eta =0.5$), due to increased spanwise shear in the near-wake region. Wavemaker results indicate that $\eta$ affects the convective nature of the shear layer more than the type of instability. DNS results show that different instability mechanisms lead to different development and evolution of vortical structures in the transition process. For $\eta =1$, the varicose instability is associated with the periodic shedding of hairpin vortices, and its stronger spatial transient growth indicated by wavemaker results aids the formation of hairpin vortices farther downstream. In contrast, for $\eta =0.5$, the interplay between varicose and sinuous instabilities results in a broader-banded energy spectrum and leads to the sinuous wiggling of hairpin vortices in the near wake when $Re_h$ is sufficiently high. A sinuous mode with a lower frequency captured by dynamic mode decomposition analysis, and associated with the ‘wiggling’ of streaks, persists far downstream and promotes transition to turbulence. A new regime map is developed to classify and predict instability mechanisms based on $Re_{hh}^{1/2}$ and $d/\delta ^{*}$ using a logistic regression model. Although the mean skin friction demonstrates different evolutions for the two geometries, both of them efficiently trip the flow to turbulence at $Re_h=1100$. An earlier location of a fully-developed turbulent state is established for $\eta =1$ at $x \approx 110h$.

Exact solutions of the time-dependent three-dimensional nonlinear vorticity equation for Euler flows with spherical geometry are provided. The velocity solution is the sum of a multipolar oscillatory function and a rigid cylindrical motion with swirl. The multipolar oscillation is a velocity mode whose radial and angular dependencies are given by the spherical Bessel functions and vector spherical harmonics, respectively. The local frequency of the velocity oscillations equals the angular speed of the rigid flow times the angular azimuthal wavenumber of the oscillating flow. The unsteady motion corresponds to inertial oscillations in multipolar flows with spatial azimuthal waves (non-vanishing azimuthal wavenumber) in the presence of a background flow with constant axial vorticity. In these nonlinear solutions, the curl of the Lamb vector has a linear dependence with the oscillation velocity, a property that makes it possible for the oscillating motion to satisfy different linear wave equations. Based on these inviscid time-dependent velocity modes, new exact solutions to the time-dependent Navier–Stokes equation are also provided.

Surface distortions to an otherwise planar channel flow introduce vorticity perturbations. In Newtonian fluids, the vorticity induced by small surface undulations on the lower wall is advected by the background flow and diffuses into the fluid. When the fluid is viscoelastic, we identify new mechanisms by which significant vorticity perturbations can be generated in both inertialess and elasto-inertial channel flows. We focus on the case where the lengthscale of the surface distortion is much longer than the channel depth, where we find significant departure from plane shear (Page & Zaki, J. Fluid Mech., vol. 901, 2016, pp. 392–429) due to the non-monotonic base-flow streamwise-normal elastic stress. In inertialess flows, a purely elastic response results in streamlines deforming to match the bottom topography in the lower half the channel. However, the vanishing stress at the centreline introduces a blocking effect, and the associated $O(1)$ jump in normal velocity is balanced by a large-amplitude streamwise-oscillating ‘jet’ in a boundary layer, resulting in a localised, chevron-shaped vorticity perturbation field. In elasto-inertial flows, resonance between the frequency of elasto-inertial ‘Alfvén’ waves and the frequency apparent to an observer moving with the fluid results in vorticity amplification in a pair of critical layers on either side of the channel. The vorticity in both layers is equal in magnitude, to leading order in Weissenberg number, and as such the perturbation vorticity field penetrates the full channel depth even when inertia is dominant. The results demonstrate that long-wave distortions, which are relatively innocuous in Newtonian fluids, can drive a significant flow distortion in viscoelastic fluids for a wide range of parameter values.

The width effect on contact angle hysteresis in a microchannel with patterned heterogeneous surfaces is systematically investigated. In the model, identical defects periodically appear on the background surface. To this end, a droplet's evaporation and condensation processes inside the microchannel are studied by theoretical analysis and numerical simulation based on a diffuse-interface lattice Boltzmann method. The microchannel width effect on the system's equilibrium properties is studied. The results demonstrate that the number of equilibrium configurations increases linearly with the microchannel width ($b$), and has a quadratic relationship with the cosine of the reference contact angle and the heterogeneity strength ($\varepsilon$). The average most stable contact angle is independent of $b$ and is always equal to the contact angle predicted by the Cassie–Baxter equation. For contact angle hysteresis ($H$), when the microchannels are narrow and wide, there are individual-effect-dominated hysteresis (IDH) and collective-effect-dominated hysteresis (CDH), respectively. The IDH and CDH are hysteresis modes corresponding to the jumping behaviour of contact lines affected by individual defects and two neighbouring defects, respectively. Based on the graphical force balance approach, we establish a scaling law to quantify the connection between $H$, $b$ and $\varepsilon$. Specifically, in the IDH mode, $H\sim b \varepsilon ^2$, while in the CDH mode, $H$ increases linearly with $\varepsilon$ but nonlinearly with $b$.

The interaction between boundary layer turbulence and a porous layer is the cornerstone of interface engineering. In this study, the spatial and spectral-resolved transfer entropy is used to assess the asymmetry of the causal interaction next to the permeable wall. The analysis is based on pore-resolved direct numerical simulation of turbulent channel flow over cylinder arrays. The spatial map of transfer entropy reveals the information flux between the porous medium and arbitrary nearby positions, and paths connecting locations with maximum information transfer are identified. The paths in the ‘top-down’ and ‘bottom-up’ directions, respectively, lean upstream and downstream, demonstrating that the coupling process is directionally dependent. The scale dependence of transfer entropy is inspected with a surrogate data strategy. As wall permeability increases, the active scale range in causal interaction shifts from near-wall vortices to Kelvin–Helmholtz type eddies. In addition, linear stochastic estimation is used to determine the statistical velocity field for a local informative event. In an average sense, the interaction between a convecting sweep or ejection event and the up/down-welling motions at the pore unit is the core mechanism that contributes to the causal coupling. The statistical findings derived from the transfer entropy are then validated using a neural network-based remote sensing model.

Diffusiophoresis refers to the motion of a colloidal particle in a solute concentration gradient, animated by particle–solute interactions. We present a theoretical analysis of the diffusiophoretic motion of a viscous drop in a gradient of neutral solute at zero Reynolds number. In a spatially uniform gradient, the translational velocity of a spherical drop was found by Anderson et al. (J. Fluid Mech., vol. 117, 1982, pp. 107–121). Here, we show additionally that the drop experiences no tendency to deform, regardless of the magnitude of the interfacial tension at the interface of the drop and suspending fluid. Next, we consider a non-uniform gradient, where the ambient solute concentration takes the form of a quadrupole around the drop centroid. This gradient does not induce drop translation, due to symmetry, but does induce a deformation in the drop shape, which is spheroidal to first order in the capillary number $Ca=\beta k_B T R^2 K/\gamma$, where $\beta$ is the magnitude of the quadrupolar variation in solute concentration, $k_B T$ is the thermal energy, $R$ is the drop radius, $K$ is the Gibbs adsorption length, and $\gamma$ is the interfacial tension. Whether the drop becomes prolate or oblate depends on whether the solute–drop interaction is attractive or repulsive. Therefore, our work shows that in principle, a drop could undergo deformation during diffusiophoresis in a non-uniform solute gradient.

The effects of insoluble surfactants, satisfying a Langmuir equation of state, on transverse Stokes flow over a superhydrophobic surface of unidirectional grooves with flat menisci are examined. The phenomenon of surface immobilisation, whereby surfactants cause part or all of the Cassie-state menisci to become effectively no-slip zones thereby degrading the slip properties of the surface, is of primary interest. The study employs a combined analytical and numerical approach allowing an exploration of surfactant effects over the full range of surface Péclet numbers, Marangoni numbers and surfactant loads. For small surface Péclet and Marangoni numbers, perturbation theory is used to gain basic insights into the physical mechanisms at work. That analysis also provides checks on a robust numerical scheme, built around a complex variable formulation of the problem, used to compute solutions across a wide range of non-dimensional parameter values. Two distinct mechanisms are found to be responsible for surface immobilisation. In the first, most commonly seen at higher surface Péclet numbers, a stagnant cap forms because swept surfactant immobilises a section of the meniscus before any portion of the meniscus reaches maximum packing. Such a cap exists even for a linear equation of state and grows in length with increasing Marangoni number. A second immobilisation mechanism is associated with the nonlinear equation of state: an immobilised region forms because the surfactant concentration reaches its maximum value near the downstream edge of the meniscus. Immobilisation of the latter form can occur at much lower surface Péclet numbers, even if the Marangoni number is small, as long as there is sufficient surfactant in the system. The study enhances understanding of how insoluble surfactants can degrade slip over superhydrophobic surfaces.

In this paper, we determine the instability effects of a phase twist superposed on a quantum vortex defect governed by the Gross–Pitaevskii equation. For this, we consider the modified form of the equation in two cases: when a uniform phase twist is present everywhere in the condensate, and when the defect is subject to a localized phase twist confined to the defect healing region. In the first case, we show that a secondary, new defect is produced as manifestation of an Aharonov–Bohm type effect. In the second case, we prove that due to energy minimization, the defect changes its configurational energy by converting localized twist to writhe. This mechanism, typical of classical elastic systems, is shown to occur also in quantum defects, and it may find useful applications in science and technology.