We develop a generalized Rotne–Prager–Yamakawa approximation for the dipolar components of the inverse friction matrix and use it for calculating the intrinsic viscosity of rigidly connected bead conglomerates. Such bead models are commonly used in the calculation of hydrodynamic properties of macromolecules. We consider both the case of non-overlapping constituent beads as well as overlapping beads of different sizes. We demonstrate the accuracy of the approximation in two test cases and show that it performs well even if the distances between the beads are small or if the beads overlap. Robust performance of this approximation in the case of overlapping beads stems from its correct limiting behaviour at a complete overlap, with one sphere fully immersed in the other. The generalized Rotne–Prager–Yamakawa approximation is thus well suited for evaluation of intrinsic viscosity, which is a key quantity in characterizing molecular conformations of biological macromolecules.

We demonstrate that the critical magnetic Reynolds number $Rm_{c}$ for a turbulent non-helical dynamo in the limit of low magnetic Prandtl number $Pm$ (i.e. $Pm=Rm/Re\ll 1$ ) can be significantly reduced if the flow is subjected to global rotation. Even for moderate rotation rates the required energy injection rate can be reduced by a factor of more than $10^{3}$ . This strong decrease in the onset is attributed to the transfer of energy to the large scales, forming a large-scale condensate, and the reduction in the turbulent fluctuations that cause the flow to have a much larger cutoff length scale than in a non-rotating flow of the same Reynolds number. The dynamo thus behaves as if it is driven just by the large scales that act as a laminar flow (i.e. it behaves as a high $Pm$ dynamo) even though the actual Reynolds number is much higher than the magnetic Reynolds number (i.e. low $Pm$ ). Our finding thus points to a new paradigm for the design of new experiments on liquid metal dynamos.

Using the energy method we investigate the stability of pure conduction in Pearson’s model for Bénard–Marangoni convection in a layer of fluid at infinite Prandtl number. Upon extending the space of admissible perturbations to the conductive state, we find an exact solution to the energy stability variational problem for a range of thermal boundary conditions describing perfectly conducting, imperfectly conducting, and insulating boundaries. Our analysis extends and improves previous results, and shows that with the energy method global stability can be proven up to the linear instability threshold only when the top and bottom boundaries of the fluid layer are insulating. Contrary to the well-known Rayleigh–Bénard convection set-up, therefore, energy stability theory does not exclude the possibility of subcritical instabilities against finite-amplitude perturbations.

The main objective of this work is to study experimentally the primary instability of non-Newtonian film flows down an inclined plane. We focus on low-concentration shear-thinning aqueous solutions obeying the Carreau law. The experimental study essentially consists of measuring wavelengths in marginal conditions, which yields the primary stability threshold for a given slope. The experimental results for neutral curves presented in the $(Re,f_{c})$ and $(Re,k)$ planes (where $f_{c}$ is the driving frequency, $k$ is the wavenumber and $Re$ is the Reynolds number) are in good agreement with the numerical results obtained by a resolution of the generalized Orr–Sommerfeld equation. The long-wave asymptotic extension of our results is consistent with former theoretical predictions of the critical Reynolds number. This is the first experimental evidence of the destabilizing effect of the shear-thinning behaviour in comparison with the Newtonian case: the critical Reynolds number is smaller, and the ratio between the critical wave celerity and the flow velocity at the free surface is larger.

The problem of subsonic flow past micro-electro-mechanical-system-type (MEMS-type) heating elements placed on a flat surface, where the MEMS devices have hump-shaped surfaces, is investigated using triple-deck theory. The compressible Navier–Stokes equations supplemented by the energy equation are considered in the limit that the Reynolds number is large. The triple-deck problem is formulated, and the linear and nonlinear analysis and results are presented. The current work is a generalisation of the problem discussed by Koroteev & Lipatov (J. Fluid Mech., vol. 707, 2012, pp. 595–605; Z. Angew. Math. Mech., vol. 77, 2013, pp. 486–493), where the MEMS devices have flat-shaped surfaces. The results show that the hump-shaped heating elements enhance large drops in pressure, and peaks and troughs in the skin friction over the centre of the hump compared with the flat-shaped devices, which may be useful for controlling the flow.

It was previously observed by Krug et al. (J. Fluid Mech., vol. 765, 2015, pp. 303–324) that the surface area $A_{\unicode[STIX]{x1D702}}$ of the turbulent/non-turbulent interface (TNTI) in gravity currents decreases with increasing stratification, significantly reducing the entrainment rate. Here, we consider the multiscale properties of this effect using direct numerical simulations of temporal gravity currents with different gradient Richardson numbers $Ri_{g}$. Our results indicate that the reduction of $A_{\unicode[STIX]{x1D702}}$ is caused by a decrease of the fractal scaling exponent $\unicode[STIX]{x1D6FD}$, while the scaling range remains largely unaffected. We further find that convolutions of the TNTI are characterized by different length scales in the streamwise and wall-normal directions, namely the integral scale $h$ and the shear scale $l_{Sk}=k^{1/2}/S$ (formed using the mean shear $S$ and the turbulent kinetic energy $k$) respectively. By recognizing that the anisotropy implied by the different scaling relations increases with increasing $Ri_{g}$, we are able to model the $Ri_{g}$ dependence of $\unicode[STIX]{x1D6FD}$ in good agreement with the data.

This work presents the analysis of a transverse combustion instability in a reduced-scale rocket engine. The study is conducted on a time-resolved database of three-dimensional fields obtained via large-eddy simulation. The physical mechanisms involved in the response of the coaxial hydrogen/oxygen flames are discussed through the analysis of the Rayleigh term in the disturbance-energy equation. The interaction between acoustics and vorticity, also explicit in the disturbance-energy balance, is shown to be the main damping mechanism for this instability. The relative contributions of Rayleigh and damping terms, depending on the position of the flame with respect to the acoustic field, are discussed. The results give new insight into the phenomenology of transverse combustion instabilities. Finally, the applicability of spectral analysis on the nonlinear Rayleigh and dissipation terms is discussed.

The flow around different prolate (needle-like) and oblate (disc-like) spheroids is studied using a multi-relaxation-time lattice Boltzmann method. We compute the mean drag coefficient $C_{D,\unicode[STIX]{x1D719}}$ at different incident angles $\unicode[STIX]{x1D719}$ for a wide range of Reynolds numbers ($\mathit{Re}$). We show that the sine-squared drag law $C_{D,\unicode[STIX]{x1D719}}=C_{D,\unicode[STIX]{x1D719}=0^{\circ }}+(C_{D,\unicode[STIX]{x1D719}=90^{\circ }}-C_{D,\unicode[STIX]{x1D719}=0^{\circ }})\sin ^{2}\unicode[STIX]{x1D719}$ holds up to large Reynolds numbers, $\mathit{Re}=2000$. Further, we explore the physical origin behind the sine-squared law, and reveal that, surprisingly, this does not occur due to linearity of flow fields. Instead, it occurs due to an interesting pattern of pressure distribution contributing to the drag at higher $\mathit{Re}$ for different incident angles. The present results demonstrate that it is possible to perform just two simulations at $\unicode[STIX]{x1D719}=0^{\circ }$ and $\unicode[STIX]{x1D719}=90^{\circ }$ for a given $\mathit{Re}$ and obtain particle-shape-specific $C_{D}$ at arbitrary incident angles. However, the model has limited applicability to flatter oblate spheroids, which do not exhibit the sine-squared interpolation, even for $\mathit{Re}=100$, due to stronger wake-induced drag. Regarding lift coefficients, we find that the equivalent theoretical equation can provide a reasonable approximation, even at high $\mathit{Re}$, for prolate spheroids.

We show that an electric field parallel to an electrically neutral surface can generate flow of electrolytic mixtures in small channels. We term this solvo-osmotic flow, since the flow is induced by the asymmetric preferential solvation of ions at the liquid–solid interface. The generated flow is comparable in magnitude to the ubiquitous electro-osmotic flow at charged surfaces, but for a fixed surface charge density, it differs qualitatively in its dependence on ionic strength. Solvo-osmotic flow can also be sensitively controlled with temperature. We derive a modified Helmholtz–Smoluchowski equation that accounts for these effects.

We present measurements of the entrainment of ambient fluid into high-Reynolds-number gravity currents produced by a steady flux of buoyancy. The currents propagate along a horizontal channel and the mixing is measured using a light attenuation technique to obtain the cross-channel average of the density throughout the current. The total volume of the current increases linearly with time, at a rate in the range $(1.8{-}2.1)Q_{o}$ for source Froude numbers, $Fr_{o}$, in the range $0.1{-}3.7$, where $Q_{o}$ is the source volume flux per unit width. Most mixing occurs either immediately downstream of the inflow or near the head of the flow, with an increasing proportion of the entrainment occurring in a mixing zone near the inflow as $Fr_{o}$ increases. A vertical gradient in the density and horizontal velocity develops in this mixing zone. This enables relatively dense fluid at the base of the current to catch up with the head, where it rises and mixes with the ambient fluid which is displaced over the head. The mixed fluid continues forward more slowly than the head, forming the relatively dilute fluid in the upper part of the current. Our data show that the depth and the depth-averaged buoyancy are primarily dependent on the position relative to the front, with the speed of the front being $\unicode[STIX]{x1D706}(Fr_{o})B_{o}^{1/3}$, where $B_{o}$ is the source buoyancy flux per unit width. Here, $\unicode[STIX]{x1D706}(Fr_{o})$ increases from 0.9 to 1.1 as $Fr_{o}$ increases from 0.1 to 3.7, while the Froude number at the head of the flow has a value of $1.1\pm 0.05$.

The present study focuses on the transition between steady convective patterns in fluid-saturated porous media. We conduct experiments to identify the transition point from the single- to double-cell pattern in a two-dimensional porous medium. We then perform a basin stability analysis to assess the relative stability of different convective modes. The resulting basin stability diagram not only provides the domains of coexistence of different modes, but it also shows that the likelihood of finding convective patterns depends strongly on the Rayleigh number. The experimentally observed transition point from single- to double-cell mode agrees well with the stochastically preferred mode inferred from the basin stability diagram.

A cellulose sponge is a mundane porous medium composed of numerous microporous cellulose sheets surrounding macroscale voids. Here, we quantify the capillary rise dynamics of non-aqueous liquids in a sponge using a combination of experiment and theory. Although the classical law of Washburn is obeyed in the early stages, the wet front propagation is no longer diffusive in the late stages and follows a power law, $h\sim t^{1/4}$, with $h$ and $t$ being the capillary rise height and time respectively. The transition of the power law is a consequence of the peculiar heterogeneous pore structure of cellulose sponges. The permeability and driving pressure change at the rise height above which the macro voids can no longer be filled completely due to significant effects of gravity. We rationalize the $t^{1/4}$ law by considering liquid flows along the corners of macro voids driven by capillary pressure of microscale wall pores.

Understanding of the roughness-induced laminar–turbulent transition of supersonic and hypersonic flows is partly challenged by the intricate sensitivities presented by different correlation criteria. We investigate experimentally the effect of height for an isolated roughness element of quadrilateral planform. Heat transfer measurements document the enhancement of roughness-induced disturbances – here the associated heat flux perturbation – along a downstream axisymmetric laminar separation. With increasing element height $k$, a gradual intensification in wake disturbance levels is found for subcritical elements ($k/\unicode[STIX]{x1D6FF}_{k}<0.15$, where $\unicode[STIX]{x1D6FF}_{k}$ is the undisturbed boundary layer thickness) while elements taller than the effective condition ($k/\unicode[STIX]{x1D6FF}_{k}\geqslant 0.32$) bypass the more moderate transition mechanisms to produce a fully turbulent element wake. Results exhibit high sensitivity to flow properties at roughness height between critical and effective conditions. A reduction in wake disturbance levels with increasing height is documented within $0.23\leqslant k/\unicode[STIX]{x1D6FF}_{k}\leqslant 0.32$. This effect coincides with a decrease in kinematic viscosity at roughness height $\unicode[STIX]{x1D708}_{k}$ (as Mach number at height $M_{k}$ increases from 1.52 to 1.96) and is restricted to elements with strong local separation, whereby the influence of local shear effects is enhanced.

By extracting unstable invariant solutions directly from body-forced three-dimensional turbulence, we study the dynamical processes at play when the forcing is large scale and unidirectional in either the momentum or the vorticity equations. In the former case, the dynamical processes familiar from recent work on linearly stable shear flows – variously called the self-sustaining process (Waleffe, Phys. Fluids, vol. 9 (4), 1997, pp. 883–900) or vortex–wave interaction (Hall & Smith, J. Fluid Mech., vol. 227, 1991, pp. 641–666; Hall & Sherwin, J. Fluid Mech., vol. 661, 2010, pp. 178–205) – are important even when the base flow is linearly unstable. In the latter case, where the forcing drives Taylor–Green vortices, a number of mechanisms are observed from the various types of periodic orbits isolated. In particular, two different transient growth mechanisms are discussed to explain the more complex states found.

The present work considers the low-Reynolds-number wake flow behind a squareback Ahmed body, in close proximity to a ground. At low Reynolds numbers such wakes are known to undergo a series of bifurcations to a state that breaks reflectional symmetry. The symmetry breaking of the wake also persists at turbulent high Reynolds numbers, where it manifests as bi-modal behaviour with random switching between the asymmetric states. Thus far, it has only been possible to study the low-Reynolds-number sequence of bifurcations experimentally and mathematically. The present work presents the first numerical simulations capturing the sequence of symmetry breaking bifurcations that occur. A study of how the wake topology changes throughout suggests that interaction between the closer top/bottom pair of parallel shear layers can only dominate once there is sufficient underbody flow. When this occurs, the two main vortex structures in the wake switch from being horizontally to vertically aligned. A linear feedback control strategy, designed to attenuate base pressure force fluctuations, is then implemented. This causes an accompanying reduction in drag and re-symmetrisation of the wake. Analysis using the dynamic mode decomposition confirms that the wake shedding mode is re-symmetrised. This work motivates future attempts to capture wake symmetry breaking and bi-modality in numerical simulations, and application of a promising feedback control strategy at higher, turbulent Reynolds numbers.

The effects of insulating lids on the convection beneath were investigated experimentally using rectangular convection cells in the flux Rayleigh number range $2.3\times 10^{9}\leqslant Ra_{F}\leqslant 1.8\times 10^{11}$ and cylindrical cells in the range $1.4\times 10^{10}\leqslant Ra_{F}\leqslant 1.2\times 10^{12}$ with the Prandtl number $Pr$ fixed at 4.3. It is found that the presence of the insulating lids leads to reduction of the global heat transfer efficiency as expected, which primarily depends on the insulating area but is insensitive to the detailed insulating patterns. At the leading-order level, the magnitude of temperature fluctuations in the bulk fluid is, again, found to be insensitive to the insulating pattern and mainly depends on the insulating area; while the temperature probability density function in the bulk is essentially invariant with respect to both the insulating area and the spatial pattern of the lids. The flow dynamics, on the other hand, is sensitive to both the covering area and the spatial distribution of the lids. At fixed $Ra_{F}$ , the flow strength is found to increase with increasing insulating area so as to transfer the same amount of heat through a smaller cooling area. Moreover, for a constant insulating area, a symmetric insulating pattern results in a symmetric flow pattern, i.e. double-roll structure; whereas an asymmetric insulating pattern leads to asymmetric flow, i.e. single-roll structure. It is further found that the symmetry breaking of the insulating pattern leads to a stronger flow that enhances the horizontal velocity more than the vertical velocity.

Because of the associated far-field logarithmic divergence, the transport problem governing two-dimensional phoretic self-propulsion lacks a steady solution when the Péclet number $\mathit{Pe}$ vanishes. This indeterminacy, which has no counterpart in three dimensions, is remedied by introducing a non-zero value of $\mathit{Pe}$, however small. We consider that problem employing a first-order kinetic model of solute absorption, where the ratio of the characteristic magnitudes of reaction and diffusion is quantified by the Damköhler number $\mathit{Da}$. As $\mathit{Pe}\rightarrow 0$ the dominance of diffusion breaks down at distances that scale inversely with $\mathit{Pe}$; at these distances, the leading-order transport represents a two-dimensional point source in a uniform stream. Asymptotic matching between the latter region and the diffusion-dominated near-particle region provides the leading-order particle velocity as an implicit function of $\log \mathit{Pe}$. Another scenario involving weak advection takes place under strong reactions, where $\mathit{Pe}$ and $\mathit{Da}$ are large and comparable. In that limit, the breakdown of diffusive dominance occurs at distances that scale as $\mathit{Da}^{2}/\mathit{Pe}$.

The fate of small particles in turbulent flows depends strongly on the velocity gradient properties of the surrounding fluid, such as rotation and strain rates. For non-inertial (fluid) particles, the restricted Euler model provides a simple low-dimensional dynamical system representation of Lagrangian evolution of velocity gradients in fluid turbulence, at least for short times. Here, we derive a new restricted Euler dynamical system for the velocity gradient evolution of inertial particles, such as solid particles in a gas, or droplets and bubbles in turbulent liquid flows. The model is derived in the limit of small (sub-Kolmogorov-scale) particles and low Stokes number. The system exhibits interesting fixed points, stability and invariant properties. Comparisons with data from direct numerical simulations show that the model predicts realistic trends such as the tendency of increased straining over rotation along heavy particle trajectories and, for light particles such as bubbles, the tendency of reduced self-stretching of the strain rate.

Experimental and numerical studies have shown that, with sufficient nonlinearity, the theoretical capillary-wave power-law spectrum derived from the kinetic equation (KE) of weak turbulence theory can be realized. This is despite the fact that the KE is derived assuming an infinite domain with continuous wavenumber, while experiments and numerical simulations are conducted in realistic finite domains with discrete wavenumbers for which the KE theoretically allows no energy transfer. To understand this, we first analyse results from direct simulations of the primitive Euler equations to elucidate the role of nonlinear resonance broadening (NRB) in discrete turbulence. We define a quantitative measure of the NRB, explaining its dependence on the nonlinearity level and its effect on the properties of the obtained stationary power-law spectra. This inspires us to develop a new quasi-resonant kinetic equation (QKE) for discrete turbulence, which incorporates the mechanism of NRB, governed by a single parameter $\unicode[STIX]{x1D705}$ expressing the ratio of NRB and wavenumber discreteness. At $\unicode[STIX]{x1D705}=\unicode[STIX]{x1D705}_{0}\approx 0.02$, the QKE recovers simultaneously the spectral slope $\unicode[STIX]{x1D6FC}_{0}=-17/4$ and the Kolmogorov constant $C_{0}=6.97$ (corrected from the original derivation) of the theoretical continuous spectrum, which physically represents the upper bound of energy cascade capacity for the discrete turbulence. For $\unicode[STIX]{x1D705}<\unicode[STIX]{x1D705}_{0}$, the obtained spectra represent those corresponding to a finite domain with insufficient nonlinearity, resulting in a steeper spectral slope $\unicode[STIX]{x1D6FC}<\unicode[STIX]{x1D6FC}_{0}$ and reduced capacity of energy cascade $C>C_{0}$. The physical insights from the QKE are corroborated by direct simulation results of the Euler equations.

For the first time, a physical mechanism is identified to explain the phase lag term in Powell’s impinging feedback loop equation (Powell, J. Acoust. Soc. Am., vol. 83 (2), 1988, pp. 515–533). Ultra-high-speed schlieren reveals a previously unseen periodic transient shock in the wall jet region of underexpanded impinging flows. The motion of this shock appears to be responsible for the production of the acoustic waves corresponding to the impingement tone. It is suggested that the delay between the inception of the shock and the formation of the acoustic wave explains the phase lag in the aeroacoustic feedback process. This suggestion is quantitatively supported through an assessment of Powell’s feedback equation, using high-resolution particle image velocimetry and acoustic measurements.