A computational analysis is performed to study the three-dimensional response of rectangular shear layers to plasma actuator-based control, in the context of sound mitigation of supersonic non-axisymmetric jets. A Mach $1.5$ rectangular jet with an aspect ratio $2:1$ is controlled using experimentally informed actuation patterns, referred to as M0, M1, M2, M3, M${\rm \pi}$ and M+/$-$1. While the first five progressively increase the phase difference between successive actuators thus enhancing three-dimensionality of the shear layer structures, the latter corresponds to the flapping mode of the jet. A preliminary linear analysis identifies that the frequency, $St\sim 1$, has a relatively high overall amplification within the baseline shear layer, and is hence utilized for control in the subsequent nonlinear simulations. Each actuation reveals unique near-field vortical and acoustic responses that have a profound impact on far-field noise levels. The M0 actuation induces circumferentially interconnected strong streamwise vortices, while M1 actuation enhances the circumferential variability in the coherent structures. The M2 actuation encompasses both these effects, and along with a very low tonal impact of forcing, produces the most desirable far-field noise mitigation (${\sim }2.6\,{\rm dB}$), contributed by a broadband reduction around the column-mode peak of the baseline jet. Beyond M2 actuation, effectiveness of control saturates, particularly along the direction of peak noise radiation. Through a near-field analysis of the acoustic component, the efficacy of M2 actuation is attributed to the attenuation of the radiative efficiency of the jet, including reduced energy in the supersonic phase speeds, and redistribution of energy into the higher helical modes. Further, it curtails the nonlinear difference interactions in the plume that energize column-mode frequencies, which often appear as strong intermittent sound-producing events. While the shear layer turbulent kinetic energy decreases with actuation, the controlled jets show minimal variations in mean flow properties, particularly under M2 actuation, suggesting this to be a promising small-perturbation-based noise control strategy.

The present study tests the efficacy of the well-known viscous drag reduction strategy of imposing spanwise wall oscillations to reduce pressure drag contributions in transitional and fully rough turbulent wall flow. This is achieved by conducting a series of direct numerical simulations of a turbulent flow over two-dimensional (spanwise-aligned) semi-cylindrical rods, placed periodically along the streamwise direction with varying streamwise spacing. Surface oscillations, imposed at fixed viscous-scaled actuation parameters optimum for smooth wall drag reduction, are found to yield substantial drag reduction ($\gtrsim$25 %) for all the rough wall cases, maintained at matched roughness Reynolds numbers. While the total drag reduction is due to a drop in both viscous and pressure drag in the case of transitionally rough flow (i.e. with large inter-rod spacing), it is associated solely with pressure drag reduction for the fully rough cases (i.e. with small inter-rod spacing), with the latter being reported for the first time. The study finds that pressure drag reduction in all cases is caused by the attenuation of the vortex shedding activity in the roughness wake, in response to wall oscillation frequencies that are of the same order as the vortex shedding frequencies. Contrary to speculations in the literature, this study confirms that the mechanism behind pressure drag reduction, achieved via imposition of spanwise oscillations, is independent of the viscous drag reduction. This mechanism is responsible for weakening of the Reynolds stresses and increase in base pressure in the roughness wake, explaining the pressure drag reduction observed by past studies, across varying roughness heights and geometries.

By examining a systematic set of direct numerical simulations, we develop a model which captures the effect of roughness density on global and local heat transfer in forced convection. The surfaces considered are zero-skewed three-dimensional sinusoidal rough walls with solidities, $\varLambda$ (defined as the frontal area divided by the total plan area), ranging from low $\varLambda = 0.09$, medium $\varLambda = 0.18$ to high $\varLambda = 0.36$. For each solidity, we vary the roughness height characterised by the roughness Reynolds number, $k^+$, from transitionally rough to fully rough conditions. The findings indicate that, as the fully rough regime is approached, there is a pronounced breakdown in the analogy between heat and momentum transfer, whereby the velocity roughness function $\Delta U^+$ continues to increase and the temperature roughness function $\Delta \varTheta ^+$ attains a peak with increasing $k^+$. This breakdown occurs at higher sand-grain roughness Reynolds numbers ($k_s^+$) with increasing solidity. Locally, we find that the heat transfer can be meaningfully partitioned into two categories: exposed, high-shear regions experiencing higher heat transfer obeying a local Reynolds analogy and sheltered, reversed-flow regions experiencing lower and spatially uniform heat transfer. The relative contribution of these distinct mechanisms to the global heat transfer depends on the fraction of the total surface area covered by these regions, which ultimately depends on $\varLambda$. These insights enable us to develop a model for the rough-wall heat-transfer coefficient, ${C_{h,k}(k^+, \varLambda, Pr)}$, where $Pr$ is the molecular Prandtl number, that assumes different heat-transfer laws in exposed and sheltered regions. We show that the exposed–sheltered surface-area fractions can be modelled through simple ray tracing that is solely dependent on the surface topography and a prescribed sheltering angle. Model predictions compare well when applied to heat-transfer data of traverse ribs from the literature.

The statistical properties of uniform momentum zones (UMZs) are extracted from laboratory and field measurements in rough wall turbulent boundary layers to formulate a set of stochastic models for the simulation of instantaneous velocity profiles. A spatiotemporally resolved velocity dataset, covering a field of view of $8 \times 9\,{\rm m}^2$, was obtained in the atmospheric surface layer using super-large-scale particle image velocimetry (SLPIV), as part of the Grand-scale Atmospheric Imaging Apparatus (GAIA). Wind tunnel data from a previous study are included for comparison (Heisel et al., J. Fluid Mech., vol. 887, 2020, R1). The probability density function of UMZ attributes such as their thickness, modal velocity and averaged vertical velocity are built at varying elevations and modelled using log-normal and Gaussian distributions. Inverse transform sampling of the distributions is used to generate synthetic step-like velocity profiles that are spatially and temporally uncorrelated. Results show that in the wide range of wall-normal distances and $Re_\tau$ up to $\sim O(10^6)$ investigated here, shear velocity scaling is manifested in the velocity jump across shear interfaces between adjacent UMZs, and attached eddy behaviour is observed in the linear proportionality between UMZ thickness and their wall normal location. These very same characteristics are recovered in the generated instantaneous profiles, using both fully stochastic and data-driven hybrid stochastic (DHS) models, which address, in different ways, the coupling between modal velocities and UMZ thickness. Our method provides a stochastic approach for generating an ensemble of instantaneous velocity profiles, consistent with the structural organisation of UMZs, where the ensemble reproduces the logarithmic mean velocity profile and recovers significant portions of the Reynolds stresses and, thus, of the streamwise and vertical velocity variability.

The finitely extensible nonlinear elastic with Peterlin closure (FENE-P) and simplified Phan-Thien–Tanner (sPTT) viscoelastic models are used widely for modelling of complex fluids. Although they are derived from distinct micro-structural theories, these models can become mathematically identical in steady and homogeneous flows with a particular choice of the values of the model parameters. However, even with this choice of parameter values, the model responses are known to differ from each other in transient flows. In this work, we investigate the responses of the FENE-P and sPTT constitutive models in large-amplitude oscillatory shear (LAOS). In steady shear, the shear stress scales with the non-dimensional group $Wi/(aL)$ ($Wi\,\sqrt {\epsilon }$) for the FENE-P (sPTT) model, where $Wi$ is the Weissenberg number, $L^2$ is the limit of extensibility in the FENE-P model ($a$ being $L^2/(L^2-3)$), and $\epsilon$ is the extensibility parameter in the sPTT model. Our numerical and analytical results show that in LAOS, the FENE-P model shows this universality only for large values of $L^2$, whereas the sPTT model shows it for all values of $\epsilon$. In the strongly nonlinear region, there is a drastic difference between the responses of the two models, with the FENE-P model exhibiting strong shear stress overshoots that manifest as self-intersecting secondary loops in the viscous Lissajous curves. We quantify the nonlinearity exhibited by each constitutive model using the sequence of physical processes framework. Despite the high degree of nonlinearity exhibited by the FENE-P model, we also show using fully nonlinear one-dimensional simulations that it does not shear band in LAOS within the range of conditions studied.

Various systems of mechanical devices and natural organisms use the repeated expansion and contraction of deformable structures to draw in and blow out fluid. Instead of such deformable continuous structures, a system that consists of multiple discrete bodies can also generate a directional flow through the cooperative movement of the individual bodies. In this study, we numerically investigate the collective effects of a multi-body system composed of eight circular cylinders, each of which oscillates separately in the radial direction to generate thrust. The cylinder array performs cooperative motion regulated by three motion parameters: phase difference, oscillation amplitude and frequency at a low Reynolds number ($Re = 10$). The phase difference between the cylinders is critical in determining the extent of the directional flow and the time-averaged thrust. The optimal phase difference that yields the maximum time-averaged thrust is consistent, regardless of the oscillation amplitude and frequency. However, the thrust generation performance becomes significantly weaker at a higher Reynolds number ($Re = 100$). This highlights that the hydrodynamic blockage in gaps between cylinders, which is induced by strong viscous diffusion at the low Reynolds number, is essential for the cooperative force generation of multiple closely spaced bodies. A new dimensionless geometric parameter based on the motion of the array is proposed to characterize the degree of bias in the generated flow and it successfully predicts the trend in the time-averaged thrust at the low Reynolds number with strong hydrodynamic blockage.

Jetting of collapsing bubbles is a key aspect in cavitation-driven fluid–solid interactions as it shapes the bubble dynamics and additionally due to its direct interaction with the wall. We study experimentally and numerically the near-wall collapse and jetting of a single bubble seeded into the stagnation flow of a wall jet, i.e. a jet that impinges perpendicular onto a solid wall. High-speed imaging shows rich and rather distinct bubble dynamics for different wall jet flow velocities and bubble-to-wall stand-off distances. The simulations use a volume-of-fluid method that allows us to numerically determine the microscopic and transient pressures and shear stresses on the wall. It is shown that a wall jet at moderate flow velocities of a few metres per second already shapes the bubble ellipsoidally inducing a planar and convergent jet flow. The distinct bubble dynamics allow us to tailor the wall interaction. In particular, the shear stresses can be increased by orders of magnitude without increasing impact pressures the same way. Interestingly, at small seeding stand-offs, the bubble during the final collapse stage can lift off the wall and migrate against the flow direction of the wall jet such that the violent collapse occurs away from the wall.

Viscous streaming is an efficient rectification mechanism to exploit flow inertia at small scales for fluid and particle manipulation. It typically entails a fluid vibrating around an immersed solid feature that, by concentrating stresses, modulates the emergence of steady flows of useful topology. Motivated by its relevance in biological and artificial settings characterized by soft materials, recent studies have theoretically elucidated, in two dimensions, the impact of body elasticity on streaming flows. Here, we generalize those findings to three dimensions, via the minimal case of an immersed soft sphere. We first improve existing solutions for the rigid-sphere limit, by considering previously unaccounted terms. We then enable body compliance, exposing a three-dimensional, elastic streaming process available even in Stokes flows. Such effect, consistent with two-dimensional analyses but analytically distinct, is validated against direct numerical simulations and shown to translate to bodies of complex geometry and topology, paving the way for advanced forms of flow control.

In this two-part study, we investigate the motion of rigid, active objects in shear Stokes flow, focusing on bodies that induce rapid rotation as part of their activity. In Part 2, we derive and analyse governing equations for rapidly spinning complex-shaped particles – general helicoidal objects with chirality. Using the multiscale framework that we develop in Part 1 (Dalwadi et al., J. Fluid Mech., vol. 979, 2024, A1), we systematically derive emergent equations of motion for the angular and translational dynamics of these chiral spinning objects. We show that the emergent dynamics due to rapid rotation can be described by effective generalised Jeffery's equations, which differ from the classic versions via the inclusion of additional terms that account for chirality and other asymmetries. Furthermore, we use our analytic results to characterise and quantify the explicit effect of rotation on the effective hydrodynamic shape of the chiral objects, expanding significantly the scope of Jeffery's seminal study.

The observed behaviour of passive objects in simple flows can be surprisingly intricate, and is complicated further by object activity. Inspired by the motility of bacterial swimmers, in this two-part study we examine the three-dimensional motion of rigid active particles in shear Stokes flow, focusing on bodies that induce rapid rotation as part of their activity. In Part 1 we develop a multiscale framework to investigate these emergent dynamics and apply it to simple spheroidal objects. In Part 2 (Dalwadi et al., J. Fluid Mech., vol. 979, 2024, A2) we apply our framework to understand the emergent dynamics of more complex shapes; helicoidal objects with chirality. Via a multiple scales asymptotic analysis for nonlinear systems, we systematically derive emergent equations of motion for long-term trajectories that explicitly account for the strong (leading-order) effects of fast spinning. Supported by numerical examples, we constructively link these effective dynamics to the well-known Jeffery's orbits for passive spheroids, deriving an explicit closed-form expression for the effective shape of the active particle, broadening the scope of Jeffery's seminal study to spinning spheroids.

The growth of perturbations subject to the Crow instability along two vortex rings of equal and opposite circulation undergoing a head-on collision is examined. Unlike the planar case for semi-infinite line vortices, the zero-order geometry of the flow (i.e. the ring radius, core thickness and separation distance) and by extension the growth rates of perturbations vary in time. The governing equations are therefore temporally integrated to characterize the perturbation spectrum. The analysis, which considers the effects of ring curvature and the distribution of vorticity within the vortex cores, explains several key flow features observed in experiments. First, the zero-order motion of the rings is accurately reproduced. Next, the predicted emergent wavenumber, which sets the number of secondary vortex structures emerging after the cores come into contact, agrees with experiments, including the observed increase in the number of secondary structures with increasing Reynolds number. Finally, the analysis predicts an abrupt transition at a critical Reynolds number to a regime dominated by a higher-frequency, faster-growing instability mode that may be consistent with the experimentally observed rapid generation of a turbulent puff following the collision of rings at high Reynolds numbers.

The linear stability of the axisymmetric steady flow in a thermocapillary liquid bridge made from 2-cSt silicone oil $({Pr}=28)$ is investigated numerically. The liquid bridge is heated either from above or below and exposed to an axial air flow which is confined to a concentric tube surrounding the bridge. At the annular inlet, the air flow is fully developed and has the same temperature as the adjacent support rod. Using an extended Oberbeck–Boussinesq approximation in which the density of both fluids depends linearly on the temperature in all equations, critical thermocapillary Reynolds numbers are obtained depending on the strength of the imposed axial air flow. The critical conditions are sensitive with respect to the direction of a weak air flow, because the air flow changes the plateau value of the interfacial temperature midway between the hot and cold ends. For stronger air flow the critical thermocapillary Reynolds number almost saturates at moderate values. Throughout, the instability arises as a hydrothermal wave with the gas phase being passive. The dynamic interface deformations for axisymmetric flow caused by the thermocapillarity flow in the liquid and by the stresses from the air flow are considered separately. Apart from turning points of the critical curve, the impact of dynamic surface deformations on the critical thermocapillary Reynolds number is moderate.

The progressive destabilisation of the incompressible flow in a cubical cavity driven by a constant shear stress is investigated numerically. To that end, one of the square faces of the cube is subjected to a constant shear stress parallel to two opposite edges of that face. The three-dimensional steady basic flow loses its mirror symmetry through a supercritical pitchfork bifurcation leading to a pair of steady stable non-symmetric flow states that are antisymmetric to each other. Upon increase of the strength of the driving, these non-symmetric equilibria become unstable via a Hopf bifurcation resulting in two limit cycles. The bifurcations are investigated using classical linear stability analyses as well as nonlinear simulations. Upon a further increase of the driving shear stress, the limit cycles destabilise through bursts triggering a complex interplay between the unstable equilibria. The transition to turbulence resembles the Pomeau–Manneville scenario.

Capillary folding is the process of folding planar objects into three-dimensional (3-D) structures using capillary force. It has many important industrial applications, e.g. the fabrication of microelectromechanical systems. In this work, we propose a 3-D model for the capillary folding of thin elastic sheets with pinned contact lines. The energy of the system consists of interfacial energies between the different phases and the elastic energy of the sheet. The latter is described by the nonlinear Koiter's model, which can accommodate large deformations of the sheet. From the energy, we derive the governing equations for the static system using a variational approach. We then develop a numerical method to find equilibrium solutions via a relaxation dynamics. At each time step, we evolve the sheet by using the subdivision element method, and update the droplet surface by minimizing a squared area functional using linear finite elements. Qualitatively, numerical solutions for the equilibrium configurations of the sheet–droplet system agree well with those observed in experiments. Furthermore, we identify the critical bendabilities and present bifurcation diagrams for the folding of a triangular sheet. The results exhibit rich and fully 3-D behaviours that were not captured by previous two-dimensional models. Our results provide new insights into the nonlinear process of capillary folding, and may contribute to the advancement of microfabrication techniques.

The derivation of the Navier–Stokes equation in continuum mechanics leads to a number of consequences which are discussed in depth. In spite of its very high representativity of real flows, this equation presents some artefacts due to the whole notion of the continuous medium. An alternative to the Navier–Stokes equation is proposed, based on the conservation of energy per unit mass instead of momentum. The classical inertial frame of reference is replaced by a set of local frames of reference where interactions are treated as cause and effect. Invoking the principle of equivalence between energy and mass, the latter is eliminated from the quantities used in this new formalism. All quantities, variables and physical properties are thus expressed in units of mass. The law of motion is established in the form of the conservation of acceleration, an energy per unit of mass and length. The acceleration is thus written in the form of a Helmholtz–Hodge decomposition, in two terms, the first curl-free and the second divergence-free as a function of two potentials, scalar and vector. Maxwell's idea of federating the laws of electrodynamics and magnetism to establish electromagnetism is taken up here to establish the new law of motion as a nonlinear wave equation. This approach makes it possible to demonstrate that this law is relativistic from the start. The form of the equation of motion in two Lagrangians gives access to symmetries related to the conservation of certain quantities according to Noether's theorem.

Analytical expressions for the velocity field and the effective slip length of pressure-driven Stokes flow through slippery pipes and annuli with rotationally symmetrical longitudinal slits are derived. Specifically, the developed models incorporate a finite local slip length and constant shear stress along the slits, and thus go beyond the assumption of perfect slip employed commonly for superhydrophobic surfaces. Thereby, they provide the possibility to assess the influence of both the viscosity of the air or other fluid that is modelled to fill the slits as well as the influence of the micro-geometry of these slits. First, expressions for tubes and annular pipes with superhydrophobic or slippery walls are provided. Second, these solutions are combined to a tube-within-a circular-pipe scenario, where one fluid domain provides a slip to the other. This scenario is interesting as an application to achieve stable fluid–fluid interfaces. With respect to modelling, it illustrates the specification of the local slip length depending on a linked flow field. The comparison of the analytically calculated solutions with numerical simulations shows excellent agreement. The results of this paper thus represent an important instrument for the design and optimization of slippage along surfaces in circular geometries.

This paper introduces a viscous vortex model for predicting the optimal drag reduction of riblet surfaces, eliminating the need for expensive direct numerical simulations (DNSs) or experiments. The footprint of a typical quasi-streamwise vortex, in terms of the spanwise and wall-normal velocities, is extracted from smooth-wall DNS flow fields in close proximity to the surface. The extracted velocities are then averaged and used as boundary conditions in a Stokes-flow problem, wherein riblets with various cross-sectional shapes are embedded. Here, the same smooth-wall-based boundary conditions can be used for riblets, as we observe from the DNSs that the quasi-streamwise vortices remain unmodified apart from an offset. In particular, the position of these vortices remain unpinned above small riblets. The present approach is compared with the protrusion-height model of Luchini et al. (J. Fluid Mech., vol. 228, 1991, pp. 87–109), which is also based on a Stokes calculation, but represents the vortex with only a uniform spanwise velocity boundary condition. The key novelty of the present model is the introduction of a wall-normal velocity component into the boundary condition, thus inducing transpiration at the riblet crests, which becomes relevant as the riblet size increases. Consequently, the present model allows for the drag-reduction prediction of riblets up to the optimal size. The present approach does not rely on the scale separation formally required by homogenisation techniques, which are only applicable for vanishingly small riblets.

We study the capillary rise of viscous liquids into sharp corners formed by two surfaces whose geometry is described by power laws $h_i(x) = c_i x^n$, $i = 1,2$, where $c_2 > c_1$ for $n \geq 1$. Prior investigations of capillary rise in sharp corners have shown that the meniscus altitude increases with time as $t^{1/3}$, a result that is universal, i.e. applies to all corner geometries. The universality of the phenomenon of capillary rise in sharp corners is revisited in this work through the analysis of a partial differential equation for the evolution of a liquid column rising into power-law-shaped corners, which is derived using lubrication theory. Despite the lack of geometric similarity of the liquid column cross-section for $n>1$, there exist a scaling and a similarity transformation that are independent of $c_i$ and $n$, which gives rise to the universal $t^{1/3}$ power law for capillary rise. However, the prefactor, which corresponds to the tip altitude of the self-similar solution, is a function of $n$, and it is shown to be bounded and monotonically decreasing as $n\to \infty$. Accordingly, the profile of the interface radius as a function of altitude is also independent of $c_i$ and exhibits slight variations with $n$. Theoretical results are compared against experimental measurements of the time evolution of the tip altitude and of profiles of the interface radius as a function of altitude.

When a liquid jet plunges into a pool, it can generate a bubble-laden jet flow underneath the surface. This common and simple phenomenon is investigated experimentally for circular jets to illustrate and quantify the role played by the net gas/liquid void fraction on the maximum bubble penetration depth. It is first shown that an increase in either the impact diameter or the jet fall height to diameter ratio at constant impact momentum leads to a reduction in the bubble cloud size. By measuring systematically the local void fraction using optical probes in the biphasic jet, it is then demonstrated that this effect is a direct consequence of the increase in air content within the cloud. A simple momentum balance model, including only inertia and the buoyancy force, is shown to predict the bubble cloud depth without any fitting parameters. Finally, a Froude number based on the bubble terminal velocity, the cloud depth and also the net void fraction is introduced to propose a simple criterion for the threshold between the inertia-dominated and buoyancy-dominated regimes.

The paper studies the statistics of total and static temperature (total $h_t$ and static $h$ enthalpy) in compressible turbulent plane channel flow using direct numerical simulations (DNS) data covering the range of centreline Mach numbers $0.3\lessapprox \bar {M}_{{{CL}}_x} \lessapprox 2.5$ and Huang–Coleman–Bradshaw friction Reynolds numbers $100\lessapprox Re_{\tau ^\star }\lessapprox 1000$. For this class of very-cold-wall flows, the DNS data for correlation coefficients and joint probability density functions (p.d.f.s) show that $h_t'$ is invariably very strongly correlated with the streamwise velocity fluctuation $u'$, in contrast to static temperature (static enthalpy $h'$) whose correlation with $u'$ weakens rapidly with increasing wall distance. We study various correlations and joint p.d.f.s of $h_t'$ and $h'$ in relation to the fluctuating velocity field, including the turbulent Prandtl number $Pr_{{T}}$, and discuss the predictions of Reynolds analogy. The scaling of the mean enthalpy and the fluctuating enthalpy variance and fluxes with respect to inner and outer velocity scales is investigated. The complex behaviour and scaling of different terms in the transport equations for the enthalpy variance and fluxes are discussed.