The density distribution around a sphere descending in a salt-stratified fluid is measured by the laser-induced fluorescence (LIF) method. The corresponding velocity distribution is measured by particle image velocimetry (PIV), and numerical simulation is also performed to supplement the observations by LIF and PIV. In steady flow, LIF observes a thin and vertically long structure which corresponds to a buoyant jet. The bell-shaped structure, which appears under strong stratification and moderate Reynolds number (Froude number $Fr \lesssim 3$, Reynolds number $50 \lesssim Re \lesssim 500$), is also identified. The measured density distributions in the salinity boundary layer and in the jet agree with the numerical simulations which use the Schmidt number of the fluorescent dye ($Sc \sim 2000$). The initially unsteady process of the jet formation is also investigated. Under weak stratification, the LIF shows an initial development of an axisymmetric rear vortex as observed in homogeneous fluids. However, as time proceeds and the effect of stratification becomes significant, the vortex shrinks and disappears, while the jet extends vertically upward. Under strong stratification, a thin jet develops without generating a rear vortex, since the effect of stratification becomes significant in a short time before the vortex is generated.

Data from direct numerical simulations of disperse bubbly flows in a vertical channel are used to study the effect of the bubbles on the carrier-phase turbulence. We developed a new method, based on an extension of the barycentric map approach, that allows us to quantify and visualize the anisotropy and componentiality of the flow at any scale. Using this we found that the bubbles significantly enhance anisotropy in the flow at all scales compared with the unladen case, and that for some bubble cases, very strong anisotropy persists down to the smallest scales of the flow. The strongest anisotropy observed was for the cases involving small bubbles. Concerning the energy transfer among the scales of the flow, our results indicate that for the bubble-laden cases, the energy transfer is from large to small scales, just as for the unladen case. However, there is evidence of an upscale transfer when considering the transfer of energy associated with particular components of the velocity field. Although the direction of the energy transfer is the same with and without the bubbles, the behaviour of the energy transfer is significantly modified by the bubbles, suggesting that the bubbles play a strong role in altering the activity of the nonlinear term in the flow. The skewness of the velocity increments also reveals a strong effect of the bubbles on the flow, changing both its sign and magnitude compared with the single-phase case. We also consider the normalized forms of the fourth-order structure functions, and the results reveal that the introduction of bubbles into the flow strongly enhances intermittency in the dissipation range, but suppresses it at larger scales. This strong enhancement of the dissipation-scale intermittency has significant implications for understanding how the bubbles might modify the mixing properties of turbulent flows.

We employ a resolvent-based methodology to estimate velocity and pressure fluctuations within turbulent channel flows at friction Reynolds numbers of approximately 180, 550 and 1000 using measurements of shear stress and pressure at the walls, taken from direct numerical simulation (DNS) databases. Martini et al. (J. Fluid Mech., vol. 900, 2021, p. A2) showed that the resolvent-based estimator is optimal when the true space–time forcing statistics are utilised, thus providing an upper bound for the accuracy of any linear estimator. We use this framework to determine the flow structures that can be linearly estimated from wall measurements, and we characterise these structures and the estimation errors in both physical and wavenumber space. We also compare these results to those obtained using approximate forcing models – an eddy-viscosity model and white-noise forcing – and demonstrate the significant benefit of using true forcing statistics. All models lead to accurate results up to the buffer layer, but only using the true forcing statistics allows accurate estimation of large-scale logarithmic-layer structures, with significant correlation between the estimates and DNS results throughout the channel. The eddy-viscosity model displays an intermediate behaviour, which may be related to its ability to partially capture the forcing colour. Our results show that structures that leave a footprint on the channel walls can be accurately estimated using the linear resolvent-based methodology, and the presence of large-scale wall-attached structures enables accurate estimations through the logarithmic layer.

The entrainment of free-stream unsteady three-dimensional vortical disturbances in the entry region of a channel is studied via matched asymptotic expansions and by numerical means. The interest is in flows at Reynolds numbers where experimental studies have documented the occurrence of intense transient growth, despite the flow being stable according to classical stability analysis. The analytical description of the vortical perturbations at the channel mouth reveals how the oncoming disturbances penetrate into the wall-attached shear layers and amplify downstream. The effects of the channel confinement, the streamwise pressure gradient and the viscous/inviscid interplay between the oncoming disturbances and the boundary-layer perturbations are discussed. The composite perturbation velocity profiles are employed as initial conditions for the unsteady boundary-region perturbation equations. At a short distance from the channel mouth, the disturbance flow is mostly confined within the shear layers and assumes the form of streamwise-elongated streaks, while farther downstream the viscous disturbances permeate the whole channel although the base flow is still mostly inviscid in the core. Symmetrical disturbances exhibit a more significant growth than anti-symmetrical disturbances, the latter maintaining a nearly constant amplitude for several channel heights downstream before growing transiently, a unique feature not reported in open boundary layers. The disturbances are more intense as the frequency decreases or the bulk Reynolds number increases. We compute the spanwise wavelengths that cause the most intense downstream growth and the threshold wall-normal wavelengths below which the perturbations are damped through viscous dissipation.

A two-dimensional model of a hydraulic fracture propagating in a weakly consolidated, highly permeable reservoir rock during a waterflooding operation is described in this paper. The model recognizes the essential differences that exist between this class of fractures and conventional hydraulic fracturing treatments of oil and gas wells, namely: (i) the large-scale perturbations of pore pressure and the associated poroelastic effects caused by extended injection time; (ii) the extremely small volume of fluid stored in the fracture compared with the injected volume; and (iii) the leakage of water from both the borehole and the propagating fracture. The model consists of a set of equations encompassing linear elastic fracture mechanics, porous media flow and lubrication theory. Three asymptotic solutions applicable at different time regimes are found theoretically, and numerical results are obtained from the discretized governing equations. The solution reveals that the injection pressure does not evolve monotonically, as it increases with time in the early time radial-flow regime but decreases in the late time fracture-flow regime. Thus, the peak injection pressure does not correspond to a breakdown of the formation, as usually assumed, but rather to a transition between two regimes of porous media flow. However, this problem exhibits an extreme sensitivity of the time scales on a dimensionless injection rate $\mathcal {I}$. If $\mathcal {I} \lessapprox 1$, the time to reach the peak pressure could become so large that it cannot be observed in field operations, i.e. the fracture remains hydraulically invisible. Finally, it is found that poroelasticity significantly affects the response of the system, by increasing the injection pressure and delaying the time at which the peak pressure takes place.

A mechanical theory is described for a phenomenon in the surgical procedure of resuscitative endovascular balloon occlusion of the aorta (REBOA). In this procedure a balloon is pushed into the aorta by a catheter and then inflated in order to stop haemorrhage. One of the hazards of this procedure is the tendency for the balloon to migrate away from its intended position. This work examines the mechanics of balloon anchoring and migration by analysing the effects of pressure waves, the sheet flow and solid friction in the thin gap between the walls of the aorta and balloon. A viscoelastic model is adopted for the aorta wall for pressure waves between the left ventricle and the balloon. The lubrication approximation is used for blood flow in the thin gap between the walls of the balloon and aorta. Samples of quantitative predictions are discussed on how the inflation pressure and balloon characteristics affect the balloon anchoring and migration. The crucial roles of solid friction and balloon placement are pointed out, which should help in guiding the manufacturing of balloons and their usage in the field.

Vortex-induced vibration of a curved flexible cylinder placed in the test section of a recirculating water tunnel and fixed at both ends is studied experimentally. Both the concave and the convex orientations (with respect to the incoming flow direction) are considered. The cylinder was hung by its own weight with a dimensionless radius of curvature of $R/D=66$, and a low mass ratio of $m^{*} = 3.6$. A high-speed imaging technique was employed to record the oscillations of the cylinder in the cross-flow direction for a reduced velocity range of $U^{*} = 3.7 - 48.4$, corresponding to a Reynolds number range of $Re= 165 - 2146$. Mono- and multi-frequency responses as well as transition from low-mode-number to high-mode-number oscillations were observed. Regardless of the type of curvature, both odd and even mode shapes were excited in the cross-flow directions. However, the response of the system, in terms of the excited modes, amplitudes and frequencies of the oscillations, was observed to be sensitive to the direction of the curvature (i.e. concave vs convex), in particular at higher reduced velocities, where mode transition occurred. Hydrogen bubble flow visualization exhibited highly three-dimensional vortex shedding patterns in the wake of the cylinder, where there existed spatial and temporal evolution of the vortex shedding modes along the length of the cylinder. The time-varying intermittent vortex shedding in the wake of the cylinder was linked to the spanwise travelling wave behaviour of the vortex-induced vibration response. The observed spatially altering wake corresponded to the multi-modal excitation and mode transition along the length of the cylinder.

Many biological microswimmers locomote by periodically beating the densely packed cilia on their cell surface in a wave-like fashion. While the swimming mechanisms of ciliated microswimmers have been extensively studied both from the analytical and the numerical point of view, optimisation of the ciliary motion of microswimmers has received limited attention, especially for non-spherical shapes. In this paper, using an envelope model for the microswimmer, we numerically optimise the ciliary motion of a ciliate with an arbitrary axisymmetric shape. Forward solutions are found using a fast boundary-integral method, and the efficiency sensitivities are derived using an adjoint-based method. Our results show that a prolate microswimmer with a $2\,{:}\,1$ aspect ratio shares similar optimal ciliary motion as the spherical microswimmer, yet the swimming efficiency can increase two-fold. More interestingly, the optimal ciliary motion of a concave microswimmer can be qualitatively different from that of the spherical microswimmer, and adding a constraint to the cilia length is found to improve, on average, the efficiency for such swimmers.

The wake of an oscillating teardrop hydrofoil with combined heaving and pitching motion was studied numerically at Reynolds number of 8000 and Strouhal numbers of $St=0.21{-}0.94$. The lower Strouhal number exhibited high efficiency propulsion with small thrust generation. However, larger thrust generation at high $St$ required more power, which lowered the propulsive efficiency. Quantitative assessment of vortex evolution, along with qualitative investigation of the formation and interaction of primary structures, revealed the association with elliptic instability characteristics for both co-rotating and counter-rotating vortex structures in both wakes. With respect to advection of the leading-edge vortex, the pressure distribution further depicted evidence of spanwise instability with distinct temporal evolution along the suction and pressure surfaces of the oscillating foil. Three-dimensional assessment of wake structures located downstream of the trailing edge depicted the existence of dislocations associated with primary vortex ‘rollers’. At low $St$, these were limited to fine spanwise corrugations (valleys and bulges) on weaker leading edge rollers, which enlarged as the rollers advected downstream. In contrast, at high $St$, the wake exhibited conjoint hairpin-horseshoe vortex structures that led to stronger deformations on the coupled vortex rollers. The statistical characteristics of secondary structures resembled the long wavelength mode and mode A identified previously for purely pitching and heaving foils, respectively. They also mimicked mode B for stationary cylinders. Novel wake models are introduced based on a complete vivid three-dimensional depiction of coherent wake structures.

The paper deals with the steady axially symmetric flow of a viscous liquid layer over a rotating disk. The liquid is fed near the axis of rotation and spreads due to inertia and the centrifugal force. The viscous shallow-water approach gives a system of ordinary differential equations governing the flow. We consider inertia, gravity, centrifugal and Coriolis forces and estimate the effect of surface tension. We found four qualitatively different flow regimes. Transition through these regimes shows the continuous evolution of the flow structure from a hydraulic jump on a static disk to a monotonic thickness decrease on a fast rotating one. We show that, in the absence of surface tension, the intensity of the jump gradually vanishes at a finite distance from the axis of rotation while the angular velocity increases. The surface tension decreases the jump radius and destroys the steady solution for a certain range of parameters.

Input–output analysis of transitional channel flows has proven to be a valuable analytical tool for identifying important flow structures and energetic motions. The traditional approach abstracts the nonlinear terms as forcing that is unstructured, in the sense that this forcing is not directly tied to the underlying nonlinearity in the dynamics. This paper instead employs a structured-singular-value-based approach that preserves certain input–output properties of the nonlinear forcing function in an effort to recover the larger range of key flow features identified through nonlinear analysis, experiments and direct numerical simulation (DNS) of transitional channel flows. Application of this method to transitional plane Couette and plane Poiseuille flows leads to not only the identification of the streamwise coherent structures predicted through traditional input–output approaches, but also the characterization of the oblique flow structures as those requiring the least energy to induce transition, in agreement with DNS studies, and nonlinear optimal perturbation analysis. The proposed approach also captures the recently observed oblique turbulent bands that have been linked to transition in experiments and DNS with very large channel size. The ability to identify the larger amplification of the streamwise varying structures predicted from DNS and nonlinear analysis in both flow regimes suggests that the structured approach allows one to maintain the nonlinear effects associated with weakening of the lift-up mechanism, which is known to dominate the linear operator. Capturing this key nonlinear effect enables the prediction of a wider range of known transitional flow structures within the analytical input–output modelling paradigm.

What is the most robust way to communicate flow trajectories? To answer this question, we employ two neural networks to respectively deconstruct (the encoder) and reconstruct (the decoder) trajectories, where information is passed between the two networks through a low-dimensional latent space in a set-up known as an autoencoder. To ensure that their communications are robust, we add noise to the coded information passed through this latent space. In the low-noise limit the latent space structures are non-spatial in nature, resembling modes of a principle component analysis (PCA). However, as the signal-to-noise ratio is decreased, we uncover Lagrangian coherent structures (LCS) as the most compact representations which still allow the decoder to accurately reconstruct trajectories. This relationship offers increased interpretability to both PCA and LCS analysis, and helps to bridge the gap between two methods of flow analysis.

Implicit large eddy simulation is performed to investigate large-scale characteristics of a temporally evolving, stably stratified turbulent shear layer arising from the Kelvin–Helmholtz instability. The shear layer at late time has two energy-containing length scales: the scale of the shear layer thickness, which characterizes large-scale motions (LSM) of the shear layer; and the larger streamwise scale of elongated large-scale structures (ELSS), which increases with time. The ELSS forms in the middle of the shear layer when the Richardson number is sufficiently large. The contribution of the ELSS to velocity and density variances becomes relatively important with time although the LSM dominate the momentum and density transport. The ELSS have a highly anisotropic Reynolds stress, to a degree similar to the near-wall region of turbulent boundary layers, while the Reynolds stress of the LSM is as anisotropic as in the outer region. Peaks in the spectral energy density associated with the ELSS emerge because of the slow decay of turbulence at very large scales. A forward interscale energy transfer from large to small scales occurs even at a small buoyancy Reynolds number. However, an inverse transfer also occurs for the energy of spanwise velocity. Negative production of streamwise velocity and density spectra, i.e. counter-gradient transport of momentum and density, is found at small scales. These behaviours are consistent with channel flows, indicating similar flow dynamics in the stratified shear layer and wall-bounded shear flows. The structure function exhibits a logarithmic law at large scales, implying a $k^{-1}$ scaling of energy spectra.

We study the effects of fluid–particle and particle–particle interactions in a three-dimensional monodispersed reactor with unstable fluidization. Simulations were conducted using the immersed boundary method for particle Reynolds numbers of 20–70 with an Archimedes number of 23 600. Two different flow regimes were identified as a function of the particle Reynolds number. For low particle Reynolds numbers ($20 < Re_p < 40$), the porosity is relatively low and the particle dynamics are dominated by interparticle collisions that produce anisotropic particle velocity fluctuations. The relative importance of hydrodynamic effects increases with increasing particle Reynolds number, leading to a minimized anisotropy in the particle velocity fluctuations at an intermediate particle Reynolds number. For high particle Reynolds numbers ($Re_p > 40$), the particle dynamics are dominated by hydrodynamic effects, leading to decreasing and more anisotropic particle velocity fluctuations. A sharp increase in the anisotropy occurs when the particle Reynolds number increases from 40 to 50, corresponding to a transition from a regime in which collision and hydrodynamic effects are equally important (regime 1) to a hydrodynamic-dominated regime (regime 2). The results imply an optimum particle Reynolds number of roughly 40 for the investigated Archimedes number of 23 600 at which mixing in the reactor is expected to peak, which is consistent with reactor studies showing peak performance at a similar particle Reynolds number and with a similar Archimedes number. Results also show that maximum effective collisions are attained at intermediate particle Reynolds number. Future work is required to relate optimum particle Reynolds number to Archimedes number.

Supercritical pressure fluids are widely used in heat transfer and energy systems. The benefit of high heat transfer performance and the successful avoidance of phase change from the use of supercritical pressure fluids are well-known, but the complex behaviours of such fluids owing to dramatic thermal property variations pose strong challenges to the design of heat transfer applications. In this paper, the turbulent flow and heat transfer of supercritical pressure $\textrm {CO}_2$ in a small vertical tube influenced by coupled effects of buoyancy and thermal acceleration are numerically investigated using direct numerical simulation. Both upward and downward flows with an inlet Reynolds number of 3540 and pressure of 7.75 MPa have been simulated and the results are compared with corresponding experimental data. The flow and heat transfer results reveal that under buoyancy and thermal acceleration, the turbulent flow and heat transfer exhibit four developing periods in which buoyancy and thermal acceleration alternately dominate. The results suggest a way to distinguish the dominant factor of heat transfer in different periods and a criterion for heat transfer degradation under the complex coupling of buoyancy and thermal acceleration. An analysis of the orthogonal decomposition and the generative mechanism of turbulent structures indicates that the flow acceleration induces a stretch-to-disrupt mechanism of coherent turbulent structures. The significant flow acceleration can destroy the three-dimensional flow structure and stretch the vortices resulting in dissipation.

Unsteady bluff-body flows exhibit dominant oscillatory behaviour owing to periodic vortex shedding. The ability to manipulate this vortex shedding is critical to improving the aerodynamic performance of bodies in a flow. This goal requires a precise understanding of how the perturbations affect the asymptotic behaviour of the oscillatory flow and of the ability to control transient dynamics. In this work, we develop an energy-efficient flow-control strategy to alter the oscillation phase of time-periodic fluid flows rapidly. First, we perform a phase-sensitivity analysis to construct a reduced-order model for the response of the flow oscillation to impulsive control inputs at various phases. Next, we introduce a real-time optimal phase-control strategy based on the phase-sensitivity function obtained by solving the associated Euler–Lagrange equations as a two-point boundary-value problem. Our approach is demonstrated for the incompressible laminar flow past a circular cylinder and an airfoil. We show the effectiveness of phase control with different actuation inputs, including blowing and rotary control. Moreover, our control approach is a sensor-based approach without the need for access to high-dimensional measurements of the entire flow field.

A formal framework to characterize and control/optimize the flow past permeable membranes by means of a homogenization approach is proposed and applied to the wake flow past a permeable cylindrical shell. From a macroscopic viewpoint, a Navier-like effective stress jump condition is employed to model the presence of the membrane, in which the normal and tangential velocities at the membrane are respectively proportional to the so-called filtrability and slip numbers multiplied by the stresses. Regarding the particular geometry considered here, a characterization of the steady flow for several combinations of constant filtrability and slip numbers shows that the flow morphology is dominantly influenced by the filtrability and exhibits a recirculation region that moves downstream of the body and eventually disappears as this number increases. A linear stability analysis further shows the suppression of vortex shedding as long as large values of the filtrability number are employed. In the control/optimization phase, specific objectives for the macroscopic flow are formulated by adjoint methods. A homogenization-based inverse procedure is proposed to obtain the optimal constrained microscopic geometry from macroscopic objectives, which accounts for fast variations of the filtrability and slip profiles along the membrane. As a test case for the proposed design methodology, a cylindrical membrane is designed to maximize the resulting drag coefficient.

This paper presents some results coming from a linear stability analysis of turbulent depth-averaged open-channel flows (OCFs) with secondary currents. The aim was to identify plausible mechanisms underpinning the formation of large-scale turbulence structures, which are commonly referred to as large-scale motions (LSMs) and very-large-scale motions (VLSMs). Results indicate that the investigated flows are subjected to a sinuous instability whose longitudinal wavelength compares very well with that pertaining to LSMs. In contrast, no unstable modes at wavelengths comparable to those associated with VLSMs could be found. This suggests that VLSMs in OCFs are triggered by nonlinear mechanisms to which the present analysis is obviously blind. We demonstrate that the existence of the sinuous instability requires two necessary conditions: (i) the circulation of the secondary currents $\omega$ must be greater than a critical value $\omega _c$; (ii) the presence of a dynamically responding free surface (i.e. when the free surface is modelled as a frictionless flat surface, no instabilities are detected). The present paper draws some ideas from the work by Cossu, Hwang and co-workers on other wall flows (i.e. turbulent boundary layers, pipe, channel and Couette flows) and somewhat supports their idea that LSMs and VLSMs might be governed by an outer-layer cycle also in OCFs. However, the presence of steady secondary flows makes the procedure adopted herein much simpler than that used by these authors.

The growth of attached bubbles during the electrochemical evolution of hydrogen at a horizontal cathode at the base of a quiescent, dilute aqueous solution is analysed using a simple model of the process that includes the Butler–Volmer reaction model, the diffusion and migration of electroactive species and a symmetry condition that approximately accounts for the presence of periodically spaced bubbles on the electrode surface. The diffusion controlled growth of a bubble approximately follows a $t^{1/2}$ law when the spacing of the bubbles on the electrode is large, departing slightly from it due to the non-uniformity of the concentration of dissolved hydrogen in the supersaturated solution into which the bubble grows, and approaches a $t^{1/3}$ law when the spacing decreases. The space- and time-averaged current density increases exponentially with the applied voltage for an alkaline solution when the consumption of water in the reaction is not taken into account. For an acidic solution, the average current density saturates to a transport limited value that depends on bubble spacing. For a given voltage, the presence of attached bubbles increases the average current density due to the decrease of the concentration overpotential caused by the bubbles. The spacing of the bubbles on the electrode surface decreases when the voltage increases if the maximum supersaturation at the electrode is imposed to be constant. The result suggests that coalescence of attached bubbles will occur above a certain voltage.

Numerical experiments on the turbulent entrainment and mixing of scalars in a incompressible flow have been performed. These simulations are based on a scale decomposition of the velocity field, thus allowing the establishment from a dynamic point of view of the evolution of scalar fields under the separate action of large-scale coherent motions and small-scale fluctuations. The turbulent spectrum can be split into active and inactive flow structures. The large-scale engulfment phenomena actively prescribe the mixing velocity by amplifying inertial fluxes and by setting the area and the fluctuating geometry of the scalar interface. On the contrary, small-scale isotropic nibbling phenomena are essentially inactive in the mixing process. It is found that the inertial mechanisms initiate the process of entrainment at large scales to be finally processed by scalar diffusion at the molecular level. This last stage does not prescribe the amount of mixing but adapts itself to the conditions imposed by the coherent anisotropic motion at large scales. The present results may have strong repercussions for the theoretical approach to scalar mixing, as anticipated here by simple heuristic arguments which are shown able to reveal the rich dynamics of the process. Interesting repercussions are also envisaged for turbulence closures, in particular for large-eddy simulation approaches where only the large scales of the velocity field are resolved.