Coherence resonance (CR) is a phenomenon in which the response of a stable nonlinear system to external noise exhibits a peak in coherence at an intermediate noise amplitude. We report the first experimental evidence of CR in a hydrodynamic system, a low-density jet capable of undergoing both supercritical and subcritical Hopf bifurcations. By applying noise to the jet in its unconditionally stable regime, we find that, for both types of bifurcation, the coherence factor peaks at an intermediate noise amplitude and increases as the stability boundary is approached. We also find that the autocorrelation function decays differently between the two types of bifurcation, indicating that CR can reveal information about the nonlinearity of a system even before it bifurcates to a limit cycle. We then model the CR dynamics with a stochastically forced van der Pol oscillator calibrated in two different ways: (i) via the conventional method of measuring the amplitude evolution in transient experiments and (ii) via the system-identification method of Lee et al. (J. Fluid Mech., vol. 862, 2019, pp. 200–215) based on the Fokker–Planck equation. We find better experimental agreement with the latter method, demonstrating the deficiency of the former method in identifying the correct form of system nonlinearity. The fact that CR occurs in the unconditionally stable regime, prior to both the Hopf and saddle-node points, implies that it can be used to forecast the onset of global instability. Although demonstrated here on a low-density jet, CR is expected to arise in almost all nonlinear dynamical systems near a Hopf bifurcation, opening up new possibilities for the development of global-instability precursors in a variety of hydrodynamic systems.

Global stability analysis is used to analyse the onset of transonic buffet on infinite swept and unswept wings. This high-Reynolds-number flow is governed by the unsteady Reynolds averaged Navier–Stokes equations. The analysis generalizes earlier studies focused on two-dimensional airfoils. For the unswept wing, results show spanwise-periodic stationary modes in addition to the earlier-observed oscillatory mode. The oscillatory mode is nominally two-dimensional with a spanwise wavelength greater than ten wing chords. The stationary modes of instability exist over two bands of spanwise wavelengths centred around an intermediate wavelength of one wing chord, and around a short wavelength of one tenth of a wing chord. The intermediate-wavelength modes have a flow structure characteristic of airfoil buffeting modes, concentrated at the shock and in the shear layer downstream of the shock. The short-wavelength modes are only concentrated in the shear layer downstream of the shock. These stationary modes can lead to spanwise-periodic flow structures for the unswept wing. For the swept wing, these stationary modes become unsteady travelling modes and contribute to the more complex buffeting-flow structures observed on swept wings as compared with unswept wings. The spanwise-wavelength bands of the travelling modes translate to different frequencies, resulting in a broad-banded unsteady response for the swept wing. For a $30^{\circ }$ swept wing, the frequencies associated with the intermediate-wavelength modes are approximately 10 times higher than the swept-wing generalization of the long-wavelength oscillatory mode, and approximately 6 times higher than the long-wavelength mode for the unswept wing. These instability characteristics are in good agreement with experimental observations.

This paper investigates the large-scale flow reorientations of Rayleigh–Bénard convection in a cubic cell using proper orthogonal decomposition (POD) analysis and modelling. A direct numerical simulation is performed for air at a Rayleigh number of $10^{7}$ and shows that the flow is characterized by four quasi-stable states, corresponding to a large-scale circulation lying in one of the two diagonal planes of the cube with a clockwise or anticlockwise motion, with occasional brief reorientations. Proper orthogonal decomposition is applied to the joint velocity and temperature fields of an enriched database which captures the statistical symmetries of the flow. We found that each quasi-stable state consists of a superposition of four spatial modes representing three types of structures: (i) a mean-flow mode consisting of two stacked counter-rotating torus-like structures; (ii) two large-scale two-dimensional rolls (pair of degenerated modes) which form large-scale diagonal rolls when combined together; and (iii) an eight-roll mode that transports fluid from one corner to the other and strengthens the circulation along the diagonal. In addition, we identified three other modes that play a role in the reorientation process: two boundary-layer modes (pair of degenerated modes) that connect the core region with the horizontal boundary layers and one mode associated with corner rolls. The symmetries of the different POD modes are discussed, as well as their temporal dynamics. A description of the reorientation process in terms of POD modes is provided and compared with other modal approaches available in the literature. Finally, Galerkin projection is used to derive a POD-based reduced-order model. Unresolved modes are accounted for in the model by an extra dissipation term and the addition of noise. A seven-mode model is able to reproduce the low-frequency dynamics of the large-scale reorientations as well as the high-frequency dynamics associated with the large-scale circulation rotation. Linear stability analysis and sensitivity analysis confirm the role of the boundary-layer modes and the corner-rolls mode in the reorientation process.

The flows around a NACA 0018 airfoil at a chord-based Reynolds number of $Re=10\,250$ and angles of attack of $\unicode[STIX]{x1D6FC}=0^{\circ }$ and $\unicode[STIX]{x1D6FC}=10^{\circ }$ are modelled using resolvent analysis and limited experimental measurements obtained from particle image velocimetry. The experimental mean velocity fields are data assimilated so that they are solutions of the incompressible Reynolds-averaged Navier–Stokes equations forced by Reynolds stress terms which are derived from experimental data. Resolvent analysis of the data-assimilated mean velocity fields reveals low-rank behaviour only in the vicinity of the shedding frequency for $\unicode[STIX]{x1D6FC}=0^{\circ }$ and none of its harmonics. The resolvent operator for the $\unicode[STIX]{x1D6FC}=10^{\circ }$ case, on the other hand, identifies two linear mechanisms whose frequencies are a close match with those identified by spectral proper orthogonal decomposition. It is also shown that the second linear mechanism, corresponding to the Kelvin–Helmholtz instability in the shear layer, cannot be identified just by considering the time-averaged experimental measurements as an input for resolvent analysis due to missing data near the leading edge. For both cases, resolvent modes resemble those from spectral proper orthogonal decomposition when the resolvent operator is low rank. The $\unicode[STIX]{x1D6FC}=0^{\circ }$ case is classified as an oscillator and its harmonics, where the resolvent operator is not low rank, are modelled using parasitic modes as opposed to classical resolvent modes which are the most amplified. The $\unicode[STIX]{x1D6FC}=10^{\circ }$ case behaves more like an amplifier and its nonlinear forcing is far less structured. The two cases suggest that resolvent-based modelling can be achieved for more complex flows with limited experimental measurements.

Spray cooling of a hot target is characterized by strong heat flux and fast change of the temperature of the wall interface. The heat flux during spray cooling is determined by the instantaneous substrate temperature, which is illustrated by boiling curves. The variation of the heat flux is especially notable during different thermodynamic regimes: film, transitional and nucleate boiling. In this study transient boiling curves are obtained by measurement of the local and instantaneous heat flux produced by sprays of variable mass flux, drop diameter and impact velocity. These spray parameters are accurately characterized using a phase Doppler instrument and a patternator. The hydrodynamic phenomena of spray impact during various thermodynamic regimes are observed using a high-speed video system. A theoretical model has been developed for heat conduction in the thin expanding thermal boundary layer in the substrate. The theory is able to predict the evolution of the target temperature in time in the film boiling regime. Moreover, a remote asymptotic solution for the heat flux during the fully developed nucleate boiling regime is developed. The theoretical predictions agree very well with the experimental data for a wide range of impact parameters.

We create a highly controlled laboratory environment – accessible to both global and local monitoring – to analyse turbulent boiling flows and in particular their shear stress in a statistically stationary state. By precisely monitoring the drag of strongly turbulent Taylor–Couette flow (the flow in between two coaxially rotating cylinders, Reynolds number $Re\approx 10^{6}$) during its transition from non-boiling to boiling, we show that the intuitive expectation, namely that a few volume per cent of vapour bubbles would correspondingly change the global drag by a few per cent, is wrong. Rather, we find that for these conditions a dramatic global drag reduction of up to 45 % occurs. We connect this global result to our local observations, showing that for major drag reduction the vapour bubble deformability is crucial, corresponding to Weber numbers larger than one. We compare our findings with those for turbulent flows with gas bubbles, which obey very different physics from those of vapour bubbles. Nonetheless, we find remarkable similarities and explain these.

A linear stability analysis of plane Poiseuille flow of an upper-convected Maxwell (UCM) fluid, bounded between rigid plates separated by a distance $2L$, has been carried out to investigate the interplay of elasticity and inertia on flow stability. The stability is governed by the following dimensionless groups: the Reynolds number $Re=\unicode[STIX]{x1D70C}U_{max}L/\unicode[STIX]{x1D702}$ and the elasticity number $E\equiv W/Re=\unicode[STIX]{x1D706}\unicode[STIX]{x1D702}/(\unicode[STIX]{x1D70C}L^{2})$, where $W=\unicode[STIX]{x1D706}U_{max}/L$ is the Weissenberg number. Here, $\unicode[STIX]{x1D70C}$ is the fluid density, $\unicode[STIX]{x1D702}$ is the fluid viscosity, $\unicode[STIX]{x1D706}$ is the micro-structural relaxation time and $U_{max}$ is the maximum base-flow velocity. The stability is analysed for two-dimensional perturbations using both pseudo-spectral and shooting methods. We also analyse the linear stability of plane Couette flow which, along with the results for plane Poiseuille flow, yields insight into the structure of the complete elasto-inertial eigenspectrum. While the general features of the spectrum for both flows remain similar, plane Couette flow is found to be stable over the range of parameters examined ($Re\leqslant 10^{4},E\leqslant 0.01$). On the other hand, plane Poiseuille flow appears to be susceptible to an infinite hierarchy of elasto-inertial instabilities. Over the range of parameters examined, there are up to seven distinct neutral stability curves in the $Re$–$k$ plane (here $k$ is the perturbation wavenumber in the flow direction). Based on the symmetry of the eigenfunctions for the streamwise velocity about the centreline, four of these instabilities are antisymmetric, while the other three are symmetric. The neutral stability curve corresponding to the first antisymmetric mode is shown to be a continuation (to finite $E$) of the Tollmien–Schlichting (TS) instability already present for Newtonian channel flow. As $E$ is increased beyond $0.0016$, a new elastic mode appears at $Re\sim 10^{4}$, which coalesces with the continuation of the TS mode for a range of $Re$, thereby yielding a single unstable mode in this range. This trend persists until $E\sim 0.0021$, beyond which this neutral curve splits into two separate ones in the $Re$–$k$ plane. The new elastic mode which arises out of this splitting has been found to be the most unstable, with the lowest critical Reynolds number $Re_{c}\approx 1210.9$ for $E=0.0066$. The neutral curves for both the continuation of the original TS mode, and the new elastic antisymmetric mode, form closed loops upon further increase in $E$, which eventually vanish at sufficiently high $E$. For $E\ll 1$, the critical Reynolds number and wavenumber scale as $Re_{c}\sim E^{-1}$ and $k_{c}\sim E^{-1/2}$ for the first two of the symmetric modal families, and as $Re_{c}\sim E^{-5/4}$ for first two of the antisymmetric modal families; $k_{c}\sim E^{-1/4}$ for the third antisymmetric family. The critical wave speed for all of these unstable eigenmodes scales as $c_{r,c}\sim E^{1/2}$ for $E\ll 1$, implying that the modes belong to a class of ‘wall modes’ in viscoelastic flows with disturbances being confined in a thin region near the wall. The present study shows that, surprisingly, even in plane shear flows, elasticity acting along with inertia can drive novel instabilities absent in the Newtonian limit.

Recent progress in understanding subcritical transition to turbulence is based on the concept of the edge, the manifold separating the basins of attraction of the laminar and the turbulent state. Originally developed in numerical studies of parallel shear flows with a linearly stable base flow, this concept is adapted here to the case of a spatially developing Blasius boundary layer. Longer time horizons fundamentally change the nature of the problem due to the loss of stability of the base flow due to Tollmien–Schlichting (TS) waves. We demonstrate, using a moving box technique, that efficient long-time tracking of edge trajectories is possible for the parameter range relevant to bypass transition, even if the asymptotic state itself remains out of reach. The flow along the edge trajectory features streak switching observed for the first time in the Blasius boundary layer. At long enough times, TS waves co-exist with the coherent structure characteristic of edge trajectories. In this situation we suggest a reinterpretation of the edge as a manifold dividing the state space between the two main types of boundary layer transition, i.e. bypass transition and classical transition.

Three examples of non-dissipative yet cumulative interaction between a single wavetrain and a single vortex are analysed, with a focus on effective recoil forces, local and remote. Local recoil occurs when the wavetrain overlaps the vortex core. All three examples comply with the pseudomomentum rule. The first two examples are two-dimensional and non-rotating (shallow water or gas dynamical). The third is rotating, with deep-water gravity waves inducing an Ursell ‘anti-Stokes flow’. The Froude or Mach number, and the Rossby number in the third example, are assumed small. Remote recoil is all or part of the interaction in all three examples, except in one special limiting case. That case is found only within a severely restricted parameter regime and is the only case in which, exceptionally, the effective recoil force can be regarded as purely local and identifiable with the celebrated Craik–Leibovich vortex force – which corresponds, in the quantum fluids literature, to the Iordanskii force due to a phonon current incident on a vortex. Another peculiarity of that exceptional case is that the only significant wave refraction effect is the Aharonov–Bohm topological phase jump.

In Rayleigh–Bénard convection experiments, the thermal coupling between the sidewall and fluid is unavoidable. As a result, the thermal properties of the sidewall can influence the flow structure that develops. To get a better understanding of the influence of the sidewall, we performed a one-to-one comparison between experiments and direct numerical simulations (DNS) in aspect ratio (diameter over height) $\unicode[STIX]{x1D6E4}=1.00$ samples. We focus on the global heat transport, i.e. the Nusselt number $Nu$, and the local vertical temperature gradients near the horizontal mid-plane on the cylinder axis and close to the sidewall. The data cover the range $10^{5}\lesssim Ra\lesssim 10^{10}$ where $Ra$ is the Rayleigh number. The $Nu$ number obtained from experimental measurements and DNS, in which we use an adiabatic sidewall, agree well. The experiments are performed with several gases, which have widely varying thermal conductivities, but all have a Prandtl number $Pr\approx 0.7$. For $Ra\gtrsim 10^{7}$, both experiments and DNS reveal a stabilizing (positive) temperature gradient at the cylinder axis. This phenomenon was known for high $Pr$, but had not been observed for small $Pr\approx 0.7$ before. The experiments reveal that the temperature gradient decreases with decreasing $Ra$ and eventually becomes destabilizing (negative). The decrease appears at a higher $Ra$ when the sidewall admittance, which measures how easily the heat transfers from the fluid to the wall, is smaller. However, the simulations with an adiabatic sidewall do not reproduce the destabilizing temperature gradient at the cylinder axis in the low $Ra$ number regime. Instead, these simulations show that the temperature gradient increases with decreasing $Ra$. We find that the simulations can reproduce the experimental findings on the temperature gradient at the cylinder axis qualitatively when we consider the physical properties of the sidewall and the thermal shields. However, the temperature gradients obtained from experiments and simulations do not agree quantitatively. The reason is that it is incredibly complicated to reproduce all experimental details accurately due to which it is impossible to reproduce all experimental measurement details. The simulations show, in agreement with the models of Ahlers (Phys. Rev. E, vol. 63 (1), 2000, 015303) and Roche et al. (Eur. Phys. J. B, vol. 24 (3), 2001, pp. 405–408), that the sidewall can act as an extra heat conductor, which absorbs heat from the fluid near the bottom plate and releases it into the fluid near the top plate. The importance of this effect decreases with increasing $Ra$. A crucial finding of the simulations is that the thermal coupling between the sidewall and fluid can strongly influence the flow structure, which can result in significant changes in heat transport. Since this effect goes beyond a simple short circuit of the heat transfer through the sidewall, it is impossible to correct experimental measurements for this effect. Therefore, careful design of experimental set-ups is required to minimize the thermal interaction between the fluid and sidewall.

In this work, we compute numerically breakage/coalescence rates and size distribution of surfactant-laden droplets in turbulent flow. We use direct numerical simulation of turbulence coupled with a two-order-parameter phase-field method to describe droplets and surfactant dynamics. We consider two different values of the surface tension (i.e. two values for the Weber number, $We$, the ratio between inertial and surface tension forces) and four types of surfactant (i.e. four values of the elasticity number, $\unicode[STIX]{x1D6FD}_{s}$, which defines the strength of the surfactant). Stretching, breakage and merging of droplet interfaces are controlled by the complex interplay among shear stresses, surface tension and surfactant distribution, which are deeply intertwined. Shear stresses deform the interface, changing the local curvature and thus surface tension forces, but also advect surfactant over the interface. In turn, local increases of surfactant concentration reduce surface tension, changing the interface deformability and producing tangential (Marangoni) stresses. Finally, the interface feeds back to the local shear stresses via the capillary stresses, and changes the local surfactant distribution as it deforms, breaks and merges. We find that Marangoni stresses have a major role in restoring a uniform surfactant distribution over the interface, contrasting, in particular, the action of shear stresses: this restoring effect is proportional to the elasticity number and is stronger for smaller droplets. We also find that lower surface tension (higher $We$ or higher $\unicode[STIX]{x1D6FD}_{s}$) increases the number of breakage events, as expected, but also the number of coalescence events, more unexpected. The increase of the number of coalescence events can be traced back to two main factors: the higher probability of inter-droplet collisions, favoured by the larger number of available droplets, and the decreased deformability of smaller droplets. Finally, we show that, for all investigated cases, the steady-state droplet size distribution is in good agreement with the $-10/3$ power-law scaling (Garrett et al., J. Phys. Oceanogr., vol. 30 (9), 2000, pp. 2163–2171), conforming to previous experimental observations (Deane & Stokes, Nature, vol. 418 (6900), 2002, p. 839) and numerical simulations (Skartlien et al., J. Chem. Phys., vol. 139 (17), 2013).

This paper uses a tight mathematical bound on the degree of the non-normality of the turbulent velocity gradient tensor to classify flow behaviour within vortical regions (where the eigenvalues of the tensor contain a conjugate pair). Structures attaining this bound are preferentially generated where enstrophy exceeds total strain and there is a positive balance between strain production and enstrophy production. Lagrangian analysis of homogeneous, isotropic turbulence shows that attainment of this bound is associated with relatively short durations and an upper limit to the spatial extent of the flow structures that is similar to the Taylor scale. An analysis of the dynamically relevant terms using a recently developed formulation (Keylock, J. Fluid Mech., vol. 848, 2018, pp. 876–904), highlights the controls on this dynamics. In particular, in high enstrophy regions it is shown that the bound is attained when normal strain decreases rather than when non-normality increases. The near absence of normal total strain results in a source of intermittency in the dynamics of dissipation that is hidden in standard analyses. It is shown that of the two terms that contribute to the non-normal production dynamics, it is the one that is typically smallest in magnitude that is of greatest importance within these $\ell =1$ filaments. The typical distance between filament centroids is just less than a Taylor scale, implying a connection to the manner in which flow topology at the Taylor scale explains dissipation at smaller scales.

Compliant lift-generating surfaces have widespread applications as marine propellers, hydrofoils and control surfaces, and the fluid–structure interactions (FSI) of such systems have important effects upon their performance and stability. Multi-phase flows like cavitation and ventilation alter the hydrodynamic and hydroelastic behaviours of lifting surfaces in ways that are not yet completely understood. This paper describes experiments on one rigid and two flexible variants of a vertical surface-piercing hydrofoil in wetted, ventilating and cavitating conditions. Tests were conducted in a towing tank and a free-surface cavitation channel. This work, which is Part 1 of a two-part series, examines the passive, or flow-induced, fluid–structure interactions of the hydrofoils. Four characteristic flow regimes are described: fully wetted, partially ventilated, partially cavitating and fully ventilated. Hydroelastic coupling is shown to increase the hydrodynamic lift and yawing moments across all four flow regimes by augmenting the effective angle of attack. The effective angle of attack, which was derived using a beam model to account for the effect of spanwise twisting deflections, effectively collapses the hydrodynamic load coefficients for the three hydrofoils. A generalized cavitation parameter, using the effective angle of attack, is used to collapse the lift and moment coefficients for all trials at a single immersed aspect ratio, smoothly bridging the four distinct flow regimes. None of the hydrofoils approached the static divergence condition, which occurs when the hydrodynamic stiffness negates the structural stiffness, but theory and experiments both show that ventilation increases the divergence speed by reducing the hydrodynamic twisting moment about the elastic axis. Coherent vortex shedding from the blunt trailing edge of the hydrofoil causes vortex-induced vibration at an approximately constant Strouhal number of 0.275 (based on the trailing edge thickness), and leads to amplified response at lock-in, when the vortex-shedding frequency approaches one of the resonant modal frequencies of the coupled fluid–structure system.

There is an overabundance of microswimmers in nature, including bacteria, algae, mammalian cells and so on. They use flagellum, cilia or global shape changes (amoeboid motion) to move forward. In the presence of confining channels, these swimmers exhibit often non-trivial behaviours, such as accumulation at the wall, navigation and so on, and their swimming speed may be strongly influenced by the geometric confinement. Several numerical studies have reported that the presence of walls either enhances or reduces the swimming speed depending on the nature of the swimmer, and also on the confinement. The purpose of this paper is to provide an analytical explanation of several previously obtained numerical results. We treat the case of amoeboid swimmers and the case of squirmers having either a tangential (the classical situation) or normal velocity prescribed at the swimmer surface (pumper). For amoeboid motion we consider a quasi-circular swimmer which allows us to tackle the problem analytically and to extract the equations of the motion of the swimmer, with several explicit analytical or semi-analytical solutions. It is found that the deformation of the amoeboid swimmer as well as a high enough order effect due to confinement are necessary in order to account for previous numerical results. The analytical theory accounts for several features obtained numerically also for non-deformable swimmers.

A simplified model is introduced to study finite-amplitude thermo-acoustic oscillations in $N$-periodic annular combustion devices. Such oscillations yield undesirable effects and can be triggered by a positive feedback between heat-release and pressure fluctuations. The proposed model, comprising the governing equations linearized in the acoustic limit, and with each burner modelled as a one-dimensional system with acoustic damping and a compact heat source, is used to study the instability caused by cross-sector coupling. The coupling between the sectors is included by solving the one-dimensional acoustic jump conditions at the locations where the burners are coupled to the annular chambers of the combustion device. The analysis takes advantage of the block-circulant structure of the underlying stability equations to develop an efficient methodology to describe the onset of azimuthally synchronized motion. A modal analysis reveals the dominance of global instabilities (encompassing the large-scale dynamics of the entire system), while a non-modal analysis reveals a strong response to harmonic excitation at forcing frequencies far from the eigenfrequencies, when the overall system is linearly stable. In all presented cases, large-scale, azimuthally synchronized (coupled) motion is observed. The relevance of the non-modal response is further emphasized by demonstrating the subcritical nature of the system’s Hopf point via an asymptotic expansion of a nonlinear model representing the compact heat source within each burner.

This research presents a new theory that explains analytically the behaviour of any fully developed incompressible turbulent pipe flow, steady or unsteady. We propose the name of theory of underlying laminar flow (TULF), because its main consequence is the description of any turbulent pipe flow as the sum of two components: the underlying laminar flow (ULF) and the purely turbulent component (PTC). We use the framework of the TULF to explain analytically most of the important and interesting phenomena reported in He & Jackson (J. Fluid Mech., vol. 408, 2000, pp. 1–38). To do so, we develop a simple model for the pressure gradient and Reynolds shear stress that could be applied to the linearly accelerated pipe flow described by He & Jackson (2000). The following features of the unsteady flow are explained: the deformation undergone by the mean velocity profiles during the transient, the velocity overshoot observed in the more rapid excursions, the dual deformation of mean velocity profiles when overshoots are present, the laminarisation effects described during acceleration, the rapid decrease of the Reynolds shear stress upon approaching the wall that brings forth the laminar sublayer, and some other minor effects. A new field is defined to characterise the degree of turbulence within the flow, directly calculable from the theory itself. Arguably, this new field would describe the degree of turbulence in a pipe flow more accurately than the familiar turbulence intensity parameter. Finally, a paradox is found in the deformation of the unsteady mean velocity profiles with respect to equal-Reynolds-number steady profiles, which is duly explained. The research also predicts the occurrence of mean velocity undershoots if the flow is decreased rapidly enough.

Superhydrophobic surfaces are capable of trapping gas pockets within the micro-roughnesses on their surfaces when submerged in a liquid, with the overall effect of lubricating the flow on top of them. These bio-inspired surfaces have proven to be capable of dramatically reducing skin friction of the overlying flow in both laminar and turbulent regimes. However, their effect in transitional conditions, in which the flow evolution strongly depends on the initial conditions, has still not been deeply investigated. In this work the influence of superhydrophobic surfaces on several scenarios of laminar–turbulent transition in channel flow is studied by means of direct numerical simulations. A single phase incompressible flow has been considered and the effect of the micro-structured superhydrophobic surfaces has been modelled imposing a slip condition with given slip length at both walls. The evolution from laminar, to transitional, to fully developed turbulent flow has been followed starting from several different initial conditions. When modal disturbances issued from linear stability analyses are used for perturbing the laminar flow, as in supercritical conditions or in the classical K-type transition scenario, superhydrophobic surfaces are able to delay or even avoid the onset of turbulence, leading to a considerable drag reduction. Whereas, when transition is triggered by non-modal mechanisms, as in the optimal or uncontrolled transition scenarios, which are currently observed in noisy environments, these surfaces are totally ineffective for controlling transition. Superhydrophobic surfaces can thus be considered effective for delaying transition only in low-noise environments, where transition is triggered mostly by modal mechanisms.

The electrohydrodynamics of a weakly conducting buoyant drop under the combined influence of gravity and a uniform electric field is studied computationally, focusing on the inertia-dominated regime. Numerical simulations are performed for both perfectly dielectric and leaky dielectric drops over a wide range of dimensionless parameters to explore the interplay of fluid inertia and electrical stress to govern the drop shape and charge convection. For perfectly dielectric drops, the fluid inertia alters the drop shape and the deformation behaviour of the drop follows a non-monotonic path. The drop shape at steady state exhibits the transition from oblate to prolate shape on increasing the electric field strength, in sharp contrast to the cases concerning the Stokes flow regime. Similar behaviour is also obtained for leaky dielectric drops for certain fluid properties. For leaky dielectric drops, the fluid inertia also affects the convective transport of charges at the drop surface and thereby alters the drop dynamics. Unlike the Stokes flow regime, where surface charge convection has little effect on the settling speed, the same modifies the drop settling speed quite significantly in the finite inertial regime depending on the combination of electrical conductivity ratio and permittivity ratio. For oblate drops at low capillary number, charge convection alters drop shape, while keeping the nature of deformation unaltered. However, for relatively large capillary number, the oblate drop transforms into a dimpled shape due to charge convection. For all cases, an interesting fact is noticed that under the combined action of electric and inertial forces, the resultant deformation is less than the summation of the deformations caused by individual effects, when inertial effects are strong. These results are likely to provide deep insights into the interplay of various nonlinearities towards altering electrohydrodynamic settling of drops and bubbles.

We investigate the effects of surfactant coating on a deformable viscous drop under the combined action of shear flow and a uniform electric field. Employing a comprehensive three-dimensional approach, we analyse the non-Newtonian shearing response of the bulk emulsion in the dilute suspension regime. Our results reveal that the location of the peak surfactant accumulation on the drop surface may get shifted from the plane of shear to a plane orthogonal to it, depending on the tilt angle of the applied electric field and strength of the electrical stresses relative to their hydrodynamic counterparts. The surfactant non-uniformity creates significant alterations in the flow perturbation around the drop, triggering modulations in the bulk shear viscosity. Overall, the shear-thinning or shear-thickening behaviour of the emulsion appears to be greatly influenced by the interplay of surface charge convection and Marangoni stresses. We show that the balance between electrical and hydrodynamic stresses renders a vanishing surface tension gradient on the drop surface for some specific shear rates, rendering negligible alterations in the bulk viscosity. This critical condition largely depends on the electrical permittivity and conductivity ratios of the two fluids and orientation of the applied electric field. Also, the physical mechanisms of charge convection and surface deformation play their roles in determining this critical shear rate. As a consequence, we obtain new discriminating factors, involving electrical property ratios and the electric field configuration, which govern the same. Consequently, the surfactant-induced enhancement or attenuation of the bulk emulsion viscosity depends on the electrical conductivity and permittivity ratios. The concerned description of the drop-level flow physics and its connection to the bulk rheology of a dilute emulsion may provide a fundamental understanding of a more complex emulsion system encountered in industrial practice.