We consider transverse shear flow over a $2L$-periodic array of air-filled grooves, which are separated by solid slats. We address the scenario where the liquid partially invades the grooves. Following Crowdy's solution of the longitudinal problem (J. Fluid Mech. vol. 925, 2021, R2), we focus upon the simplified geometry of infinitely narrow slats. For flat menisci, the ratio $\lambda$ of the slip length to $L$ depends upon a single geometric parameter, namely the ratio $H$ of the invasion depth to $L$. We analyse the singular limit $H\ll 1$ using matched asymptotic expansions, with an outer region on the scale of the period and an inner region on the scale of a single slat. In the outer region, we employ the complex-variable formulation of Stokes flows, obtaining the complex velocity as a divergent series. This formal representation is transformed to an admissible closed-form expression. In the inner region, the problem is equivalent to that of Stokes flow past a finite line segment. Asymptotic matching yields $\lambda = [\ln (2/{\rm \pi} H)-1]/{\rm \pi}$.

In particle-laden turbulent wall flows, transport of particles towards solid walls is phenomenologically thought to be governed by the wall-normal turbulence intensity supporting the underlying particle–eddy interactions that are usually modelled by a combination of turbophoresis and turbulent diffusion. We estimate the turbophoretic and turbulent diffusive coefficients as a function of wall-normal coordinate directly from a generated direct numerical simulation (DNS) database of low volume fraction point particles in a turbulent pipe flow. These coefficients are then used in an advection–diffusion equation to estimate the particle concentration as a function of wall-normal distance and time, with favourable comparison against DNS for smaller Stokes number ($St^+$) particles suggesting a limitation of the common gradient diffusion hypothesis for larger $St^+$ particles. Using DNS we explore the non-trivial effects of $St^+$, pipe wall condition (particle absorbing or elastic) as well as the influence of drag and lift force on the velocity and particle statistics giving rise to different particle concentrations. We then appraise various Eulerian-based models of turbophoretic and turbulent diffusive coefficients and, finally, use physical insights from Lagrangian correlation times, conditional quadrant analysis and flow topology to shed further light on the particle transport as a function of various parameters and the limits of gradient diffusion hypothesis.

Diffusiophoresis takes place when a particle in solution moves due to the presence of a solute concentration gradient. This phenomenon is often studied under some simplifying assumptions, such as negligible diffusive layer thickness or infinite diffusion coefficient. In this work we simulate diffusiophoresis without these simplifications. The goal of this numerical study is to investigate equilibrium and fully developed states of non-electrolyte phoretic systems. Simulation results show that equilibrium states depend on solute diffusivity and on a reference solute concentration far from the particle. An expression is regressed that gives the (equilibrium) diffusiophoretic velocity as a function of solute concentration gradient, solute diffusion coefficient and the reference solute concentration far from the particle. A different set of results reveals that the state of phoretic systems does not depend on the initial conditions when time goes to infinity. This motivates the definition of fully developed states, designating those systems whose properties no longer depend on initial conditions. Apart from these findings, this work also depicts the effect of solute–interface interactions on diffusiophoresis. Simulation results for two solid particles with different interaction potentials are used to illustrate particle separation via diffusiophoresis. Finally, values of particle mobility are calculated for different solute–interface attraction strengths. These results are compared with another work in the literature, which studies polymer diffusiophoresis via molecular simulations (Ramírez-Hinestrosa et al., J. Chem. Phys., vol. 152, 2020, p. 164901).

We perform two-way coupled direct numerical simulation of particle-laden flow in an open channel at a friction Reynolds number ($Re_{\tau }$) of 5186, which exhibits many characteristics of high-Reynolds-number wall-bounded turbulence, such as the distinct separation of scales in the inner and outer layers. Three representative cases, an unladen case and low- and high-Stokes-number particle-laden cases, are performed to investigate the turbulent modification by particles. To this end, we compare several statistical quantities to understand the particle effect on momentum exchange and interphasial energy transfer. The modulation of large-scale motions (LSMs) and very-large-scale motions (VLSMs) are analysed using spectral information, and we find that the LSMs and VLSMs are generally weakened in the inner and outer layers, which is qualitatively different from similar simulations at lower Reynolds numbers ($Re_{\tau } \approx 500$). The spatial structures are investigated with correlation analysis, and inclined VLSMs are observed in the near-wall region, with decreased inclination angles by particles. The particles tend to widen and shorten the spanwise and streamwise extent of coherent structures, respectively. Furthermore, we find that the vorticity vector displays a preferential alignment with the eigenvector corresponding to the intermediate eigenvalue of the strain-rate tensor, independent of the particle Stokes number.

Morphodynamic equations governing the behaviour of active nematic fluids on deformable curved surfaces are constructed in the large deformation limit. Emphasis is placed on the formulation of objective rates that account for normal deformations whilst ensuring that tangential flows are Eulerian, and the use of the surface derivative (rather than the covariant derivative) in the nematic free energy, which elastically couples local order to out-of-plane bending of the surface. Focusing on surface geometry and its dynamical interplay with the hydrodynamics, several illustrative instabilities are then characterised. These include cases where the role of the Scriven–Love number and its nematic analogue are non-negligible, and where the active nematic forcing can be characterised by an analogue of the Föppl–von Kármán number. For the former, flows and changes to the nematic texture are coupled to surface geometry by viscous dissipation. This is shown to result in non-trivial relaxation dynamics for a nematic tube. For the latter, the nematic active forcing couples to the surface bending terms of the nematic free energy, resulting in extensile (active ruffling) and contractile (active pearling) instabilities in the tube shape, as well as active bend instabilities in the nematic texture. In comparison to the flat case, such bend instabilities now have a threshold set by the extrinsic curvature of the tube. Finally, we examine a topological defect located on an almost flat surface, and show that there exists a steady state where a combination of defect elasticity, activity and non-negligible spin connection drive a shape change in the surface.

The stability and sensitivity of two- and three-dimensional global modes developing on steady spanwise-homogeneous laminar separated flows around NACA 4412 swept wings are numerically investigated for different Reynolds numbers ${\textit {Re}}$ and angles of attack $\alpha$. The wake dynamics is driven by the two-dimensional von Kármán mode whose emergence threshold in the $\alpha \unicode{x2013}{\textit {Re}}$ plane is computed with that of the three-dimensional centrifugal mode. At the critical Reynolds number, the Strouhal number, the streamwise wavenumber of the von Kármán mode and the spanwise wavenumber of the leading three-dimensional centrifugal mode scale as a power law of $\alpha$. The introduction of a sweep angle attenuates the growth of all unstable modes and entails a Doppler effect in the leading modes’ dynamics and a shift towards non-zero frequencies of the three-dimensional centrifugal modes. These are found to be non-dispersive as opposed to the von Kármán modes. The sensitivity of the leading global modes is investigated in the vicinity of the critical conditions through adjoint-based methods. The growth-rate sensitivity map displays a region on the suction side of the wing, wherein a streamwise-oriented force has a net stabilising effect, comparable to what could have been obtained inside the recirculation bubble. In agreement with the predictions of the sensitivity analysis, a spanwise-homogeneous force suppresses the Hopf bifurcation and stabilises the entire branch of von Kármán modes. In the limit of small amplitudes, passive control via spanwise-wavy forcing produces a stabilising effect similar to that of a spanwise-homogeneous control and is more effective than localised spherical forces.

This study proposes a newly developed deep-learning-based method to generate turbulent inflow conditions for spatially developing turbulent boundary layer (TBL) simulations. A combination of a transformer and a multiscale-enhanced super-resolution generative adversarial network is utilised to predict velocity fields of a spatially developing TBL at various planes normal to the streamwise direction. Datasets of direct numerical simulation (DNS) of flat plate flow spanning a momentum thickness-based Reynolds number, $Re_\theta = 661.5\unicode{x2013}1502.0$, are used to train and test the model. The model shows a remarkable ability to predict the instantaneous velocity fields with detailed fluctuations and reproduce the turbulence statistics as well as spatial and temporal spectra with commendable accuracy as compared with the DNS results. The proposed model also exhibits a reasonable accuracy for predicting velocity fields at Reynolds numbers that are not used in the training process. With the aid of transfer learning, the computational cost of the proposed model is considered to be effectively low. Furthermore, applying the generated turbulent inflow conditions to an inflow–outflow simulation reveals a negligible development distance for the TBL to reach the target statistics. The results demonstrate for the first time that transformer-based models can be efficient in predicting the dynamics of turbulent flows. They also show that combining these models with generative adversarial networks-based models can be useful in tackling various turbulence-related problems, including the development of efficient synthetic-turbulent inflow generators.

The high-Reynolds-number stratified wake of a slender body is studied using a high-resolution hybrid simulation. The wake generator is a 6 : 1 prolate spheroid with a tripped boundary layer, the diameter-based body Reynolds number is ${Re}= U_\infty D/\nu = 10^5$, and the body Froude numbers are ${Fr}=U_\infty /ND=\{2,10,\infty \}$. The wake defect velocity decays following three stages with different wake decay rates (Spedding, J. Fluid Mech., vol. 337, 1997, pp. 283–301) as for a bluff body. However, the transition points among stages do not follow the expected $Nt = Nx/U_\infty$ values. Comparison with the wake of a circular disk in similar conditions (Chongsiripinyo & Sarkar, J. Fluid Mech., vol. 885, 2020) quantifies the influence of the wake generator – bluff versus slender – in stratified flow. The strongly stratified ${Fr}=2$ wake is in a resonant state. The steady lee waves strongly modulate the mean flow, and relative to the disk, the 6 : 1 spheroid (a high-aspect-ratio shape) wake at ${Fr}=2$ shows an earlier transition from the non-equilibrium (NEQ) stage to the quasi-two-dimensional (Q2D) stage. The NEQ–Q2D transition is followed by a sharp increase in the turbulent kinetic energy and horizontal wake meanders. At ${Fr}=10$, the start of the NEQ stage is delayed for the spheroid. Transfers between kinetic energy and potential energy reservoirs (both mean and turbulence) are analysed, and the flows are compared in phase space (with local Froude and Reynolds numbers as coordinates). Overall, the results of this study point to the difficulty of finding a universal framework for stratified wake evolution, independent of the features of the body, and provide insights into how buoyancy effects depend on the wake generator.

In this work, the stability of hypersonic flow over a curved compression ramp is studied using several stability analysis tools and direct numerical simulations (DNS). The free-stream Mach number and the unit Reynolds number are 7.7 and $4.2 \times 10^6$ m$^{-1}$, respectively. Corner rounding is considered to alter the separation bubble flow so as to suppress the intrinsic instability of the compression-ramp flow. The variation of intrinsic instability is confirmed by global stability analysis. Subsequently, resolvent analysis is employed to examine the response of intrinsically stable flows to external disturbances. It is shown that the considered flows strongly amplify low-frequency streamwise streaks with a preferential spanwise wavelength. This result is verified using DNS by introducing a random forcing upstream of the separation point. Furthermore, both resolvent analysis and DNS demonstrate that the separation bubble contributes little to the selection of the spanwise wavelength of streamwise streaks. The combined effects of convective and intrinsic instabilities are also explored using DNS. A better agreement with experimental data is achieved after introducing upstream disturbances in an inherently unstable flow.

When a viscous fluid spreads underneath a deformable surface skin or crust, the peeling dynamics at the fluid front can control the rate of advance rather than bulk self-similar flow. For an elastic skin, this control results in a quasi-static interior blister held at constant pressure that is matched to a narrow peeling region behind the fluid front. In this paper, the analogous problem is considered for a skin that deforms either viscously or plastically, or both. In particular, the deformable surface is assumed to be a thin plate of material governed by the Herschel–Bulkley constitutive law. We examine how such a skin controls viscous flow underneath, fed at constant flux and spreading as either a planar or axisymmetric current. As for an elastic skin, the peeling dynamics at the viscous fluid front again controls the rate of spreading. However, contrary to that situation, the mathematical matching problem for viscoplastic peeling is simplified considerably as a result of an integral constraint. Despite this, the structure of the peeling region is complicated significantly by any plasticity in the skin, which can create a convoluted peeling wave ahead of the main blister that features interwoven yielded and plugged sections of the plate.

Predicting the performance of large turbine arrays requires the understanding of many physical factors, such as array geometry, turbine operation, inflow conditions and turbulent wake mixing. Due to the large parameter space that an array may be optimised over, low-order models with low computational cost are often employed. This paper extends one of these models, the inviscid–viscous coupled model, for multi-row turbine modelling. Firstly, an extension to the inviscid actuator disc theory is presented by removing the limit on the number of discrete streamtubes computed. The extended model allows for the quantification of the impact of shear in the bypass and core flows separately. In particular, it is shown that averaging a sheared bypass flow profile can result in a substantial over-prediction of the power of a turbine in a laterally bounded flow as the effective blockage of the flow increases. The model is also used to confirm that an approximation using a limited number of streamtubes in some previous applications of the inviscid–viscous approach has a negligible impact on the results. Secondly, we explore the performance of a multi-row array with either uniform or varying turbine resistance across different rows. Results suggest that by varying resistance across rows, the array may outperform the uniform resistance case. The performance gain is dependent, however, on the arrangement and inter-turbine spacing both in the spanwise and streamwise directions.

The present study investigates the linear spatio-temporal and weakly nonlinear stability of a pressure-driven two-layer channel flow subjected to a wall-normal temperature gradient commonly encountered in industrial applications. The liquid–liquid interface tension is assumed to be a linearly decreasing function of temperature. The study employs both numerical (pseudo-spectral method) and long-wave approaches. The general linear stability analysis (GLSA) predicts shear-flow and thermocapillary modes that arise due to the imposed pressure and temperature gradients, respectively. The previous stability analyses of the same problem predicted a negligible effect of the pressure-driven flow on the linear stability of the system. However, the GLSA reveals stabilising and destabilising effects of the pressure-driven flow depending on the viscosity ratio ($\mu _r$), thermal conductivity ratio ($\kappa _r$), interface position ($H$) and the sign of the imposed temperature gradient ($\beta _1$). The analysis predicts a range of $H$ for given $\mu _r$ and $\kappa _r$, which can not be stabilised by the thermocapillarity. The numerically predicted long-wave instability is then captured using the long-wave asymptotic approach. The arguments based on the physical mechanism further successfully explain the role of $\mu _r$, $\kappa _r$, $H$, the sign of $\beta _1$ and the interaction between the velocity and temperature perturbations in stabilising/destabilising the flow. The spatio-temporal analysis reveals the dominance of the spanwise mode in causing the absolutely unstable flow. The weakly nonlinear analysis reveals a subcritical pitchfork bifurcation without shear flow. However, with the shear flow, the streamwise mode undergoes a supercritical Hopf bifurcation.

We propose an analytical model to estimate the interface temperature $\varTheta _{\varGamma }$ and the Nusselt number $Nu$ for an evaporating two-layer Rayleigh–Bénard configuration in statistically stationary conditions. The model is based on three assumptions: (i) the Oberbeck–Boussinesq approximation can be applied to the liquid phase, while the gas thermophysical properties are generic functions of thermodynamic pressure, local temperature and vapour composition, (ii) the Grossmann–Lohse theory for thermal convection can be applied to the liquid and gas layers separately and (iii) the vapour content in the gas can be taken as the mean value at the gas–liquid interface. We validate this setting using direct numerical simulations in a parameter space composed of the Rayleigh number ($10^6\leq Ra\leq 10^8$) and the temperature differential ($0.05\leq \varepsilon \leq 0.20$), which modulates the variation of state variables in the gas layer. To better disentangle the variable property effects on $\varTheta _\varGamma$ and $Nu$, simulations are performed in two conditions. First, we consider the case of uniform gas properties except for the gas density and gas–liquid diffusion coefficient. Second, we include the variation of specific heat capacity, dynamic viscosity and thermal conductivity using realistic equations of state. Irrespective of the employed setting, the proposed model agrees very well with the numerical simulations over the entire range of $Ra$–$\varepsilon$ investigated.

We experimentally explored the effect of single-sidewall cooling on Rayleigh–Bénard (RB) convection. Canonical RB was also studied to aid insight. The scenarios shared tank dimensions and bottom and top wall temperatures; the single sidewall cooling had the top wall temperature. Turbulence was explored at two canonical Rayleigh numbers, $Ra=1.6\times 10^{10}$ and $Ra=2\times 10^9$ under Prandtl number $Pr=5.4$. Particle image velocimetry described vertical planes parallel and perpendicular to the sidewall cooling. The two $Ra$ scenarios reveal pronounced changes in the flow structure and large-scale circulation (LSC) due to the sidewall cooling. The density gradient induced by the sidewall cooling led to asymmetric descending and ascending flows and irregular LSC. Flow statistics departed from the canonical case, exhibiting lower buoyancy effects, represented by an effective Rayleigh number with effective height dependent on the distance from the lateral cooling. Velocity spectra show two scalings, $\varPhi \propto f^{-5/3}$ Kolmogorov (KO41) and $\varPhi \propto f^{-11/5}$ Bolgiano (BO59) in the larger $Ra$; the latter was not present in the smaller set-up. The BO59 scaling with sidewall cooling appears at higher frequencies than its canonical counterpart, suggesting weaker buoyancy effects. The LSC core motions allowed us to identify a characteristic time scale of the order of vortex turnover time associated with distinct vortex modes. The velocity spectra of the vortex core oscillation along its principal axis showed a scaling of $\varPhi _c \propto f^{-5/3}$ for the single sidewall cooling, which was dominant closer there. It did not occur in the canonical case, evidencing the modulation of LSC oscillation on the flow.

Particle effects on the amplitude modulation are investigated in this study based on observational data with various mass loading acquired from long-term measurements of aeolian sandstorms in high-Reynolds-number ($Re_{\tau }\sim O(10^6)$) near-neutral atmospheric surface layers. In both particle-laden and unladen flows, in addition to the positive top–down modulation behaviour in the logarithmic region, a significant modulation effect that exists for some specific motions is also found for the single-point amplitude modulation. The most energetic turbulent motions exhibit the strongest modulation effect, and the modulating signals do not change with the small-scale motions being modulated. In particle-laden flows, the length of the most energetic structure is almost constant, thus the scales of the modulating signal and carrier signal are hardly affected by particles. However, the addition of particles changes the distribution of energy between multi-scale turbulent motions. The kinetic energy of the large-scale component is less enhanced than the total kinetic energy by particles. This leads to a reduced energy proportion of the large-scale component and an augmented one of the small-scale component. Moreover, the particles produce a large damping in the degree of the amplitude modulation and move down the positions of the modulating signals and carrier signals corresponding to the strongest inter-layer modulation, but the damping is weakened with the wall-normal distance due to the decreased mass loading. This study may provide a more general insight into the modulation mechanism between multi-scale turbulent motions and the effect of particles on turbulence.

A closed-form solution is presented for internal compressible flows with an arbitrary combination of area change, heat addition, friction and non-uniform flow. The approach assumes steady flow of a calorically perfect ideal gas using the integral form of the governing equations. It is shown that the equations reduce to a simple, bi-quadratic equation in exit Mach number. This solution is compared to previously known solutions for Rayleigh flow, Fanno flow, isentropic flow and normal shock relations. Theoretical solutions for sudden expansion, sudden contraction, and supersonic–supersonic two-stream constant area mixing are also presented. These solutions were found to be in agreement with data available in the literature.

Experiments are reported that explore the onset of motion of bubbles in a model yield stress fluid, Carbopol gel. Starting from a trapped spherical bubble in a gel, the yielding limit for the bubble motion is obtained by gradually expanding the bubble via a stepwise decrease in pressure. Our results show that at the yielding limit bubbles are longer and thinner when they are in a higher concentrated gel. This is suggestive of a link between the shape and size of the bubbles at the onset of motion and the rheology of the material, in particular elastic behaviour below the yielding point. Particular attention has been paid to investigating the dynamic response of gel during the bubble growth. Subjecting the bubble to a periodic change in the pressure confirms the irreversibility of the gel deformation and its hysteresis, which are hallmarks of nonlinear viscoelastic behaviour of the gel before yielding. In this context, the periodic expansion and contraction of the bubbles leave residual deformation (stresses) in the gel which facilitates the liberation of bubbles.

We examine the effect of Dean number on the inertial focusing of spherical particles suspended in flow through curved microfluidic ducts. Previous modelling of particle migration in curved ducts assumed the flow rate was small enough that a leading-order approximation of the background flow with respect to the Dean number produces a reasonable model. Herein, we extend our model to situations involving a moderate Dean number (in the microfluidics context) while the particle Reynolds number remains small. Variations in the Dean number cause a change in the axial velocity profile of the background flow which influences the inertial lift force on a particle. Simultaneously, changes in the cross-sectional velocity components of the background flow directly affect the secondary flow induced drag. In keeping the particle Reynolds number small, we continue to approximate the inertial lift force using a regular perturbation while capturing the subtle effects from the modified background flow. This approach pushes the limits at which a regular perturbation is applicable to provide some insights into how variations in the Dean number influence particle focusing. Our results illustrate that, as the extrema in the background flow move towards the outside of edge of the cross-section with increasing Dean number, we observe a similar shift in the stable equilibria of some, but not all, particle sizes. This might be exploited to enhance the lateral separation of particles by size in a number of practical scenarios.

We carry out a weakly nonlinear analysis of the centrifugal instability for a columnar vortex in a rotating fluid, and compare the results to those of the semi-linear model derived empirically by Yim et al. (J. Fluid Mech., vol. 897, 2020, A34). The asymptotic analysis assumes that the Reynolds number is close to the instability threshold so that the perturbation is only marginally unstable. This leads to two coupled equations that govern the evolutions of the amplitude of the perturbation and of the mean flow under the effect of the Reynolds stresses due to the perturbation. These equations differ from the Stuart–Landau amplitude equation or coupled amplitude equations involving a mean field that have been derived previously. In particular, the amplitude does not saturate to a constant as in the supercritical Stuart–Landau equation, but decays afterwards reflecting the instability disappearance when the mean flow tends toward a neutrally stable profile in the direct numerical simulations (DNS). These equations resemble those of the semi-linear model except that the perturbation in the weakly nonlinear model keeps at leading order the structure of the eigenmode of the unperturbed base flow. The predictions of the weakly nonlinear equations are compared to those of the semi-linear model and to DNS for the Rossby number $Ro=-4$ and various Reynolds numbers and wavenumbers. They are in good agreement with the DNS when the growth rate is sufficiently small. However, the agreement deteriorates and becomes only qualitative for parameters away from the marginal values, whereas the semi-linear model continues to be in better agreement with the DNS.

Responses of a geometry-induced separation bubble (GISB) and a pressure-induced separation bubble (PISB) at enhanced levels of free-stream turbulence (FST) have experimentally been investigated for a comparative study using the particle image velocimetry (PIV) technique. The outlines of separation bubbles based on the dividing streamlines are self-similar for different levels of FST and Reynolds numbers. The spectral analyses of the time-resolved PIV data show that the vortex shedding frequency of a separated shear layer remains unchanged for the GISB cases even with an enhanced level of FST. In contrast, it is different for the PISB cases. We propose a criterion that determines whether the frequency will remain the same even for the cases with FST. Linear stability analyses reveal that the inviscid-inflectional instability dominates the transition process, and the linear stages of transition are not completely bypassed even at an enhanced level of FST. The most amplified frequencies, while scaling with the displacement thickness and the boundary layer edge velocity, collapse in a single curve for all the cases. Furthermore, measurements in the spanwise plane show that the streamwise velocity streak/Klebanof mode at an enhanced level of FST is not a general flow feature for all types of separation bubbles. However, at an enhanced level of FST for the PISB case, the boundary layer streaks are found to distort the two-dimensional vortex structure associated with the Kelvin–Helmholtz instability, eventually leading to a three-dimensional $\varLambda$-like structure in the spanwise plane.