The interaction between deformable surfaces and oscillatory driving is known to produce complex secondary time-averaged flows due to inertial and elastic nonlinearities. Here, we revisit the problem of oscillatory flow in a cylindrical tube with a deformable wall, and analyse it under a long-wave theory for small deformations, but for arbitrary Womersley numbers. We find that the oscillatory pressure does not vary linearly along the length of a deformable channel, but instead decays exponentially with spatial oscillations. We show that this decay occurs over an elasto-visco-inertial length scale that depends on the material properties of the fluid and the elastic walls, the geometry of the system, and the frequency of the oscillatory flow, but is independent of the amplitude of deformation. Inertial and geometric nonlinearities associated with the elastic deformation of the channel drive a time-averaged secondary flow. We quantify the flow using numerical solutions of the perturbation theory, and gain insight by developing analytic approximations. The theory identifies a complex non-monotonic dependence of the time-averaged flux on the elastic compliance and inertia, including a reversal of the flow. Finally, we show that our analytic theory is in excellent quantitative agreement with the three-dimensional direct numerical simulations of Pande et al. (Phys. Rev. Fluids, vol. 8, no. 12, 2023, 124102).

We show how a complete mathematical model of a physical process can be learned directly from data via a hybrid approach combining three simple and general ingredients: physical assumptions of smoothness, locality and symmetry, a weak formulation of differential equations and sparse regression. To illustrate this, we extract a complete system of governing equations of fluid dynamics – the Navier–Stokes equation, the continuity equation and the boundary conditions – as well as the pressure-Poisson and energy equations, from numerical data describing a highly turbulent channel flow in three dimensions. Whether they represent known or unknown physics, relations discovered using this approach take the familiar form of partial differential equations, which are easily interpretable and readily provide information about the relative importance of different physical effects. The proposed approach offers insight into the quality of the data, serving as a useful diagnostic tool. It is also remarkably robust, yielding accurate results for very high noise levels, and should thus be well suited for analysis of experimental data.

The dewetting of thin nanofilms is significantly affected by thermal fluctuations, liquid–solid slip and disjoining pressure, which can be described by lubrication equations augmented by appropriately scaled noise terms, known as stochastic lubrication equations. Here molecular dynamics simulations along with a newly proposed slip-generating method are adopted to study the instability of nanofilms with arbitrary slip. These simulations show that strong-slip dewetting is distinct from weak-slip dewetting by faster growth of perturbations and fewer droplets after dewetting, which cannot be predicted by the existing stochastic lubrication equation. A new stochastic lubrication equation considering the strong-slip boundary condition is thus derived using a long-wave approximation to the equations of fluctuating hydrodynamics. The linear stability analysis of this equation, i.e. surface spectrum, agrees well with molecular simulations. Interestingly, strong slip can break down the usual Stokes limits adopted in weak-slip dewetting and bring the inertia into effect. The evolution of the standard deviation of the film height $W^2(t)={\overline {h^2}-{\bar {h}}^2}$ at the initial stage of the strong-slip dewetting is found to be $W\sim t^{1/4}$ in contrast to $W\sim t^{1/8}$ for the weak-slip dewetting.

The overall system of interest is an infinite half-space in which a compressible ocean is the top layer and an elastic seafloor (together with the crust beneath) forms a semi-infinite bottom layer. Whereas water-column compression waves and seafloor waves individually have received considerable attention, not much is known about their propagation as groups. This work utilizes the group behaviour of these waves to derive energy balance relations for wavenumber spectra for wave groups propagating through a mildly non-uniform water-column–seafloor system. Dispersion relations for the coupled system are derived using known kinematic and kinetic conditions at the interface, and free and forced wave solutions for the wavenumber spectra are obtained, with particular attention to the case when certain frequency–wavenumber combinations in the forcing excite the two-media system into resonance. Wavenumber spectra predicted using the theory for mildly non-uniform media are found to be close to those predicted assuming uniform media, though the effect of non-uniformity becomes more noticeable as the groups propagate farther from the generation area. Here, nonlinear interactions among stationary, random multi-directional surface-wave fields provide the forcing for groups of compression waves in the water and surface waves on the seafloor. The formulation includes the cumulative effect of multiple generation areas along the group propagation direction. Comparisons with observational data from a sensor array in the Atlantic Ocean indicate that the theory can be applied to reconstruct plausible combinations of generation areas and interaction times that are consistent with the measured data, for deriving approximate predictions at down-wave distances along the group propagation directions. Implications of this and other findings are discussed for (i) the potential for energy conversion from water-column compression waves on the seafloor, (ii) tracking of tropical cyclones from the seafloor, and (iii) quantification and comparative assessment of low-frequency mid-ocean ambient noise and microseism activity.

This study explores the application of a wall-attached ferrofluid film to decrease skin-friction drag in turbulent channel flow. We conduct experiments using water as a working fluid in a turbulent channel flow set-up, where one wall is coated with a ferrofluid layer held in place by external permanent magnets. Depending on the flow conditions, the interface between the two fluids is observed to form unstable travelling waves. While ferrofluid coating has been previously employed in laminar and moderately turbulent flows (Reynolds number $Re<4000$) to reduce drag by creating a slip condition at the fluid interface, its effectiveness in fully developed turbulent conditions, particularly when the interface exhibits instability, remains uncertain. Our primary objective is to assess the effectiveness of ferrofluid coating in reducing turbulent drag with particular focus on scenarios when the ferrofluid layer forms unstable waves. To achieve this, we measure flow velocity using two-dimensional particle tracking velocimetry (2D-PTV), and the interface contour between the fluids is determined using an interface tracking algorithm. Our results reveal the significant potential of ferrofluid coating for drag reduction, even in scenarios where the interface between the surrounding fluid and ferrofluid exhibits instability, with observed drag reduction rates up to 95 %. In particular, waves with an amplitude significantly smaller than a viscous length scale positively contribute to drag reduction, while larger waves are detrimental, because of induced turbulent fluctuations. However, for the latter case, slip outcompetes the extra turbulence so that drag is still reduced.

The coupling of Richtmyer–Meshkov instability (RMI) and Kelvin–Helmholtz instability (KHI), referred to as RM-KHI, on a shock-accelerated inclined single-mode air–SF$_6$ interface is studied through shock-tube experiments, focusing on the evolution of the perturbation distributed along the inclined interface. To clearly capture the linear (overall linear to nonlinear) evolution of RM-KHI, a series of experiments with a weak (relatively strong) incident shock is conducted. For each series of experiments, various $\theta _{i}$ (angle between incident shock and equilibrium position of the initial interface) are considered. The nonlinear flow features manifest earlier and develop faster when $\theta _{i}$ is larger and/or shock is stronger. In addition, the interface with $\theta _{i}>0^{\circ }$ evolves obliquely along its equilibrium position under the effect of KHI. RMI dominates the early-time amplitude evolution regardless of $\theta _{i}$ and shock intensity, which arises from the discrepancy in the evolution laws between RMI and KHI. KHI promotes the post-early-stage amplitude growth and its contribution is related positively to $\theta _{i}$. An evident exponential-like amplitude evolution behaviour emerges in RM-KHI with a relatively strong shock and large $\theta _{i}$. The linear model proposed by Mikaelian (Phys. Fluids, vol. 6, 1994, pp. 1943–1945) is valid for RM-KHI within the linear period. In contrast, the adaptive vortex model (Sohn et al., Phys. Rev. E, vol. 82, 2010, p. 046711) can effectively predict both the interface morphology and overall amplitude evolutions from the linear to nonlinear regimes.

Frozen water might appear opaque since gas bubbles can get trapped in the ice during the freezing process. They nucleate and then grow near the advancing solidification front, due to the formation of a gas supersaturation region in its vicinity. A delicate interplay between the rate of mass transfer and the rate of freezing dictates the final shapes and sizes of the entrapped gas bubbles. In this work, we experimentally and numerically investigate the initial growth of such gas bubbles that nucleate and grow near the advancing ice front. We show that the initial growth of these bubbles is governed by diffusion and is enhanced due to a combination of the presence of the background gas concentration gradient and the motion of the approaching front. Additionally, we recast the problem into that of mass transfer to a moving spherical object in a homogeneous concentration field, finding good agreement between our experimental data and the existing scaling relations for that latter problem. Lastly, we address how fluid flow around the bubble might further affect this growth and qualitatively explore this through numerical simulations.

Direct numerical simulation (DNS) of rotating pipe flows up to $Re_\tau \approx 3000$ is carried out to investigate drag reduction effects associated with axial rotation, extending previous studies carried out at a modest Reynolds number (Orlandi & Fatica, J. Fluid Mech., vol. 343, 1997, pp. 43–72; Orlandi & Ebstein, Intl J. Heat Fluid Flow, vol. 21, 2000, pp. 499–505). The results show that the drag reduction, which we theoretically show to be equivalent to net power saving assuming no mechanical losses, monotonically increases as either the Reynolds number or the rotation number increases, proportionally to the inner-scaled rotational speed. Net drag reduction up to approximately $70\,\%$ is observed, while being far from flow relaminarisation. Scaling laws for the mean axial and azimuthal velocity are proposed, from which a predictive formula for the friction factor is derived. The formula can correctly represent the dependency of the friction factor on the Reynolds and rotation numbers, maintaining good accuracy for low-to-moderate rotation numbers. Examination of the turbulent structures highlights the role of rotation in widening and elongating the small-scale streaks, with subsequent suppression of sweeps and ejections. In the core part of the flow, clear weakening of large-scale turbulent motions is observed at high Reynolds numbers, with subsequent suppression of the outer-layer peak in the pre-multiplied spectra. The Fukagata–Iwamoto–Kasagi decomposition indicates that, consistent with a theoretically derived formula, the outer layer yields the largest contribution to drag reduction at increasingly high Reynolds numbers. In contrast, both the inner and the outer layers contribute to drag reduction as the rotation number increases.

The development of a bubble plume from a vertical gas-evolving electrode is driven by buoyancy and hydrodynamic bubble dispersion. This canonical fluid mechanics problem is relevant for both thermal and electrochemical processes. We adopt a mixture model formulation for the two-phase flow, considering variable density (beyond Boussinesq), viscosity and hydrodynamic bubble dispersion. Introducing a new change of coordinates, inspired by the Lees–Dorodnitsyn transformation, we obtain a new self-similar solution for the laminar boundary layer equations. The results predict a wall gas fraction and gas plume thickness that increase with height to the power of 1/5 before asymptotically reaching unity and scaling with height to the power 2/5, respectively. The vertical velocity scales with height to the power of 3/5. Our analysis shows that self-similarity is only possible if gas conservation is entirely formulated in terms of the gas specific volume instead of the gas fraction.

Direct numerical simulations of the droplet impact on a flat solid surface with an annular part are conducted. We investigate droplet impact on a superhydrophobic substrate with a superhydrophilic annulus to understand the formation conditions of droplets in different states. The location and size of superhydrophilic annulus are carried out through the phase diagram. We describe the formation process of droplets in three different states and the spreading radius with time to catch the rupture time of the film. Two different ruptures occur in the spreading stage or the retraction stage, respectively. The rupture times from these two mechanisms observed numerically are found to be a key factor resulting in partial rebound and lens-shaped/ring-shaped droplets. Finally, the influence of non-dimensional numbers on the formation of the ring-shaped droplet is demonstrated. The Weber number can alter the amplitude of the up and down oscillation on the droplet's upper surface, while the Froude number affects primarily the time to form the central penetrating hole. This gives the guidance and method to control the ring-shaped droplets formation time.

An experimental investigation is conducted to study the flow patterns, spectral properties and energy fluxes in thin-layer turbulence with varying system sizes and damping rates. It is found that although a system-size vortex (an indicator of spectral condensation) occurs for small system sizes and does not for large ones, the spectra for different system sizes consistently exhibit a scaling close to $k^{-3}$ in inverse cascade (another indicator of spectral condensation). On the other hand, under a fixed system size larger than the friction-dominated length scale, the energy spectrum in the inverse cascade range changes from $k^{-3}$ to $k^{-5/3}$ as the damping rate increases, suggesting that the friction-dominated length scale may not be a suitable parameter for predicting spectral transition. At lower damping rates and large system sizes, turbulent structures grow larger via inverse cascade, manifesting as long streamers, and the small-scale vortices are suppressed. This suppression leads to a reduction of energy flux at intermediate scales and a change in the spectral shape. The dimensionless Taylor microscale is found to exhibit a monotonic dependence on the damping rate. With the reduction in the damping rate, the Taylor microscale increases to become comparable with the forcing scale, and the spectrum in inverse cascade transits to a steeper scaling, $k^{-3}$, indicating that the dimensionless Taylor microscale may be used as a diagnostic parameter for spectral transition.

Cross-stream migration of a deformable fluid particle is investigated computationally in a pressure-driven channel flow of a viscoelastic fluid via interface-resolved simulations. Flow equations are solved fully coupled with the Giesekus model equations using an Eulerian–Lagrangian method and extensive simulations are performed for a wide range of flow parameters to reveal the effects of particle deformability, fluid elasticity, shear thinning and fluid inertia on the particle migration dynamics. Migration rate of a deformable particle is found to be much higher than that of a solid particle under similar flow conditions mainly due to the free-slip condition on its surface. It is observed that the direction of particle migration can be altered by varying shear thinning of the ambient fluid. With a strong shear thinning, the particle migrates towards the wall while it migrates towards the channel centre in a purely elastic fluid without shear thinning. An onset of elastic flow instability is observed beyond a critical Weissenberg number, which in turn causes a path instability even for a nearly spherical particle. An inertial path instability is also observed once particle deformation exceeds a critical value. Shear thinning is found to be suppressing the path instability in a viscoelastic fluid with a high polymer concentration whereas it reverses its role and promotes path instability in a dilute polymer solution. It is found that migration of a deformable particle towards the wall induces a secondary flow with a velocity that is approximately an order of magnitude higher than the one induced by a solid particle under similar flow conditions.

This study presents direct numerical simulation results of two-layer Rayleigh–Bénard convection, investigating the previously unexplored Rayleigh–Weber parameter space $10^6\leq Ra\leq 10^8$ and $10^2\leq We\leq 10^3$. Global properties, such as the Nusselt and Reynolds numbers, are compared against the extended Grossmann–Lohse theory for two fluid layers, confirming a weak Weber number dependence for all global quantities and considerably larger Reynolds numbers in the lighter fluid. Statistics of the flow reveal that the interface fluctuates more intensely for larger Weber and smaller Rayleigh numbers, something also reflected in the increased temperature root mean square values next to the interface. The dynamics of the deformed two-fluid interface is further investigated using spectral analysis. Temporal and spatial spectrum distributions reveal a capillary wave range at small Weber and large Rayleigh numbers, and a secondary energy peak at smaller Rayleigh numbers. Furthermore, the maxima of the space–time spectra lie in an intermediate dispersion regime, between the theoretical predictions for capillary and gravity-capillary waves, showing that the gravitational energy of the interfacial waves is strongly altered by temperature gradients.

Linear stability analysis currently fails to predict turbulence transition in canonical viscous flows. We show that two alternative models of the boundary condition for incipient perturbations at solid walls produce linear instabilities that could be sufficient to explain turbulence transition. In many cases, the near-wall behaviour of the discovered instabilities is empirically indistinguishable from the classical no-slip condition. The ability of these alternative boundary conditions to predict linear instabilities that are consistent with turbulence transition suggests that the no-slip condition may be an overly simplified model of fluid–solid interface physics, particularly as a description of the flow perturbations that lead to turbulence transition in wall-bounded flows.

In this study, we analyse ‘magneto-Stokes’ flow, a fundamental magnetohydrodynamic (MHD) flow that shares the cylindrical-annular geometry of the Taylor–Couette cell but uses applied electromagnetic forces to circulate a free-surface layer of electrolyte at low Reynolds numbers. The first complete, analytical solution for time-dependent magneto-Stokes flow is presented and validated with coupled laboratory and numerical experiments. Three regimes are distinguished (shallow-layer, transitional and deep-layer flow regimes), and their influence on the efficiency of microscale mixing is clarified. The solution in the shallow-layer limit belongs to a newly identified class of MHD potential flows, and thus induces mixing without the aid of axial vorticity. We show that these shallow-layer magneto-Stokes flows can still augment mixing in distinct Taylor dispersion and advection-dominated mixing regimes. The existence of enhanced mixing across all three distinguished flow regimes is predicted by asymptotic scaling laws and supported by three-dimensional numerical simulations. Mixing enhancement is initiated with the least electromagnetic forcing in channels with order-unity depth-to-gap-width ratios. If the strength of the electromagnetic forcing is not a constraint, then shallow-layer flows can still yield the shortest mixing times in the advection-dominated limit. Our robust description of momentum evolution and mixing of passive tracers makes the annular magneto-Stokes system fit for use as an MHD reference flow.

Lozano-Duran et al. (J. Fluid Mech., vol. 914, 2021, p. A8) have recently identified the ability of streamwise-averaged turbulent streak fields ${\mathcal {U}}(y,z,t)\hat {\boldsymbol {x}}$ in minimal channels to produce short-term transient growth as the key linear mechanism needed to sustain turbulence at $Re_{\tau }=180$. Here, in an attempt to extend this result to larger domains and higher $Re_{\tau }$, we model this streak transient growth as a two-stage linear process by first selecting the dominant streak structure expected to emerge over the eddy turnover time on the turbulent mean profile $U(y)\hat {\boldsymbol {x}}$, and then examining the secondary growth on this (frozen) streak field ${\mathcal {U}}(y,z)\hat {\boldsymbol {x}}$. Choosing the mean streak amplitude and eddy turnover time consistent with simulations captures the growth thresholds found by Lozano-Duran et al. (2021) for sustained turbulence. In a larger domain at $Re_{\tau }=180$, the most energetic near-wall streaks observed in simulations are close to the predicted optimal streaks. This most energetic streak spacing, approaches the optimal streak at $Re_{\tau }=550$ where the secondary growth possible on each also comes together. A key prediction from the model is that the threshold transient growth required to sustain turbulence decreases with increasing $Re_{\tau }$. More fundamentally, the work of Lozano-Duran et al. (2021) and our results suggest a subtle but significant revision of Malkus's (J. Fluid Mech., vol. 1, 1956, pp. 521–539) classic hypothesis concerning realisable turbulent mean profiles. The key property for a realisable turbulent mean profile could be the ability to generate sufficient short-term transient growth rather than dependence on its (long-term) linear stability characteristics, which was Malkus's original idea.

We explore the drawing of a shear-thinning or shear-thickening thread with an axisymmetric hole that evolves due to axial drawing, inertia and surface tension effects. The stress is assumed to be proportional to the shear rate raised to the $n$th power. The presence of non-Newtonian rheology and surface tension forces acting on the hole introduces radial pressure gradients that make the derivation of long-wavelength equations significantly more challenging than either a Newtonian thread with a hole or shear-thinning and shear-thickening threads without a hole. In the case of weak surface tension, we determine the steady-state profiles. Our results show that for negligible inertia the hole size at the exit becomes smaller as $n$ is decreased (i.e. strong shear-thinning effects) above a critical draw ratio, but surprisingly the opposite is true below this critical draw ratio. We determine an accurate estimate of the critical draw ratio and also discuss how inertia affects this process. We further show that the dynamics of hole closure is dominated by a different limit, and we determine the asymptotic forms of the hole closure process for shear-thinning and shear-thickening fluids with inertia. A linear instability analysis is conducted to predict the onset of draw resonance. We show that increased shear thinning, surface tension and inlet hole size all act to destabilise the flow. We also show that increasing shear-thinning effects reduce the critical Reynolds number required for unconditional stability. Our study provides valuable insights into the drawing process and its dependence on the physical effects.