Relatively strongly stratified turbulent flows tend to self-organise into a ‘layered anisotropic stratified turbulence’ (LAST) regime, characterised by relatively deep and well-mixed density ‘layers’ separated by relatively thin ‘interfaces’ of enhanced density gradient. Understanding the associated mixing dynamics is a central problem in geophysical fluid dynamics. It is challenging to study LAST mixing, as it is associated with Reynolds numbers and Froude numbers ( and being characteristic velocity and length scales, the kinematic viscosity and the buoyancy frequency). Since a sufficiently large dynamic range (largely) unaffected by stratification and viscosity is required, it is also necessary for the buoyancy Reynolds number , where is the (appropriately volume-averaged) turbulent kinetic energy dissipation rate. This requirement is exacerbated for oceanically relevant flows, as the Prandtl number in thermally stratified water (where is the thermal diffusivity), thus leading (potentially) to even finer density field structures. We report here on four forced fully resolved direct numerical simulations of stratified turbulence at various Froude () and Prandtl () numbers forced so that , with resolutions up to . We find that, as increases, emergent ‘interfaces’ become finer and their contribution to bulk mixing characteristics decreases at the expense of the small-scale density structures populating the well-mixed ‘layers’. However, extreme mixing events (as quantified by significantly elevated local destruction rates of buoyancy variance ) are always preferentially found in the (statically stable) interfaces, irrespective of the value of .
We study numerically the microjetting mode obtained when a fluid is injected through a tube submerged in a uniaxial extensional flow. The steady solution to the full nonlinear Navier–Stokes equations is calculated. We obtain the linear global modes determining the linear stability of the steady solution. For sufficiently large outer viscosity, the flow remains stable for infinitely small values of the injected flow rate. This implies that jets with vanishing diameters can be produced regardless of the jet viscosity and outer flow strength. For a sufficiently small inner-to-outer viscosity ratio, the microjetting instability is associated only with the flow near the entrance of the jet. The tapering meniscus stretches and adopts a slender quasiconical shape. Consequently, the cone tip is exposed to an intense outer flow, which stabilizes the flow in the cone–jet transition region. This work presents the first evidence that fluid jets with arbitrarily small diameters can be stably produced via tip streaming. The results are related to those of a droplet in a uniaxial extensional flow with its equator pinned to an infinitely thin ring. The pinning of the equator drastically affects the droplet stability and breakup.
Amplification of velocity gradients, a key feature of turbulent flows, is affected by the non-local character of the incompressible fluid equations expressed by the second derivative (Hessian) of the pressure field. By analysing the structure of the flow in regions where the vorticity is the highest, we propose an approximate expression for the pressure Hessian in terms of the local vorticity, consistent with the existence of intense vortex tubes. Contrary to the often used simplification of an isotropic form for the pressure Hessian, which in effect inhibits vortex stretching, the proposed approximate form of the pressure Hessian enables much stronger vortex stretching. The prediction of the approximation proposed here is validated with results of direct numerical simulations of turbulent flows.
Subcritical pipe flow transition has received a great deal of attention over the past decades, as it constitutes a quintessential bifurcation process between two metastable fluid states: the laminar and turbulent solutions. Coherent lower-branch structures, forming flow states that facilitate between these two attracting equilibria, have been proposed that together form an edge manifold in phase space separating relaminarizing from transitioning perturbations. Typically, direct numerical simulations or low-dimensional model equations have been used to study this edge manifold with bisection methods. In the article by Kaszás & Haller (J. Fluid Mech., vol. 979, 2024, A48), an effective nonlinear invariant-manifold technique has been applied to extract a low-dimensional, global representation of the phase-space dynamics directly from simulation data. It allows the computation of the intersection of the edge manifold with a low-dimensional surface that is strikingly accurate in predicting the long-term dynamics of perturbations about the lower-branch solution and thus provides an accessible parameterization of the edge manifold for subcritical pipe flow transition.
The theory of transient growth describes how linear mechanisms can cause temporary amplification of disturbances even when the linearized system is asymptotically stable as defined by its eigenvalues. This growth is traditionally quantified by finding the initial disturbance that generates the maximum response at the peak time of its evolution. However, this can vastly overstate the growth of a real disturbance. In this paper, we introduce a statistical perspective on transient growth that models statistics of the energy amplification of the disturbances. We derive a formula for the mean energy amplification and spatial correlation of the growing disturbance in terms of the spatial correlation of the initial disturbance. The eigendecomposition of the correlation provides the most prevalent structures, which are the statistical analogue of the standard left singular vectors of the evolution operator. We also derive accurate confidence bounds on the growth by approximating the probability density function of the energy. Applying our analysis to Poiseuille flow yields a number of observations. First, the mean energy amplification is often drastically smaller than the maximum. In these cases, it is exceedingly unlikely to achieve near-optimal growth due to the exponential behaviour observed in the probability density function. Second, the characteristic length scale of the initial disturbances has a significant impact on the expected growth, with large-scale initial disturbances growing orders of magnitude more than small-scale ones. Finally, while the optimal growth scales quadratically with Reynolds number, the mean energy amplification scales only linearly for certain reasonable choices of the initial correlation.
The squirmer is a popular model to analyse the fluid mechanics of a self-propelled object, such as a micro-organism. We demonstrate that some fore–aft symmetric squirmers can spontaneously self-propel above a critical Reynolds number. Specifically, we numerically study the effects of inertia on spherical squirmers characterised by an axially and fore–aft symmetric ‘quadrupolar’ distribution of surface-slip velocity; under creeping-flow conditions, such squirmers generate a pure stresslet flow, the stresslet sign classifying the squirmer as either a ‘pusher’ or ‘puller’. Assuming axial symmetry, and over the examined range of the Reynolds number (defined based upon the magnitude of the quadrupolar squirming), we find that spontaneous symmetry breaking occurs in the puller case above , with steady swimming emerging from that threshold consistently with a supercritical pitchfork bifurcation and with the swimming speed growing monotonically with .
The synchronisation between rotating turbulent flows in periodic boxes is investigated numerically. The flows are coupled via a master–slave coupling, taking the Fourier modes with wavenumber below a given value as the master modes. It is found that synchronisation happens when exceeds a threshold value , and depends strongly on the forcing scheme. In rotating Kolmogorov flows, does not change with rotation in the range of rotation rates considered, being the Kolmogorov length scale. Even though the energy spectrum has a steeper slope, the value of is the same as that found in isotropic turbulence. In flows driven by a forcing term maintaining constant energy injection rate, synchronisation becomes easier when rotation is stronger. Here, decreases with rotation, and it is reduced significantly for strong rotations when the slope of the energy spectrum approaches . It is shown that the conditional Lyapunov exponent for a given is reduced by rotation in the flows driven by the second type of forcing, but it increases mildly with rotation for the Kolmogorov flows. The local conditional Lyapunov exponents fluctuate more strongly as rotation is increased, although synchronisation occurs as long as the average conditional Lyapunov exponents are negative. We also look for the relationship between and the energy spectra of the Lyapunov vectors. We find that the spectra always seem to peak at approximately , and synchronisation fails when the energy spectra of the conditional Lyapunov vectors have a local maximum in the slaved modes.
The axially symmetric deformation of a drop in a viscous fluid, under the influence of an externally imposed flow having simultaneous rotating and compressional or extensional components, is addressed. In the previous studies, two families of stationary drop shapes were constructed by simulating the dynamics of drop deformation: stable singly connected shapes with respect to axisymmetric disturbances, and unstable toroidal shapes. These two branches coexist at the same flow conditions, but were not connected. In this study, we obtain a new family of branches of unstable highly deformed stationary drops connecting with the stable flattened shapes and the toroidal ones. We use a method based on classical control theory. The controller is designed for a two-state dynamic model of the system and is employed on a high-order nonlinear dynamic model of the drop deformation. Through this feedback-control-centred approach, an extended collection of unstable stationary solutions is constructed, which spans the range from the loss of stability to the dimpled shapes almost collapsed at the centre. In the latter region, which was never obtained in previous studies, a multiplicity of solutions is identified.
The development of simple, low-order and accurate unsteady aerodynamic models represents a crucial challenge for the design optimisation and control of fluid dynamical systems. In this work, wind tunnel experiments of a pitching NACA 0018 aerofoil conducted at a Reynolds number and at different free-stream turbulence intensities are used to identify data-driven nonlinear state-space models relating the time-varying angle of attack of the aerofoil to the lift coefficient. The proposed state-space neural network (SS-NN) modelling technique explores an innovative methodology, which brings the flexibility of artificial neural networks into a classical state-space representation and offers new insights into the construction of reduced-order unsteady aerodynamic models. The work demonstrates that this technique provides accurate predictions of the nonlinear unsteady aerodynamic loads of a pitching aerofoil for a wide variety of angle-of-attack ranges and frequencies of oscillation. Results are compared with a modified version of the Goman–Khrabrov dynamic stall model. It is shown that the SS-NN methodology outperforms the classical semi-empirical dynamic stall models in terms of accuracy, while retaining a fast evaluation time. Additionally, the proposed models are robust to noisy measurements and do not require any pre-processing of the data, thus involving only a limited user interaction. Overall, these features make the SS-NN technique an excellent candidate for the construction of accurate data-driven models from experimental fluid dynamics data, and pave the way for their adoption in applications entailing design optimisation and real-time control of systems involving lift.
In this paper, we study wave scattering and radiation by a surface-piercing vertical truncated metamaterial cylinder composed of a closely spaced array of thin vertical barriers, between which fluid can flow. A theoretical model is developed under full depth-dependent linearised water wave theory, where an effective medium equation and effective boundary conditions are employed, respectively, to describe the fluid motion inside the cylinder and match the flow between the fluid regions in and outside the metamaterial cylinder. A damping mechanism is introduced at the surface of the fluid occupied by the metamaterial cylinder to consider the wave power dissipation in narrow gaps between the thin vertical plates. The wave excitation forces acting on the cylinder and the hydrodynamic coefficients can be calculated straightforwardly in terms of the velocity potential inside the cylinder. An alternative way is by using the velocity potential outside the cylinder, the expression of which has the reduction of the integral and an infinite accumulation that are included in the straightforward expression. The results highlight the patterns of the radiated waves induced by the oscillation of the cylinder and the characteristics of the hydrodynamic coefficients. The metamaterial cylinder when fixed in place and with a damping mechanism included is found to capture more wave power than that of a traditional axisymmetric heaving wave energy converter over a wide range of wave frequencies.
The wake of two tandem square cylinders of identical width (d) is experimentally studied, with a view to understanding the dependence of the flow structure, aerodynamics forces and Strouhal number on the centre-to-centre spacing ratio L/d and Reynolds number Re, where L is the distance between the cylinder centres. Extensive measurements are carried out, using hot-wire, particle imaging velocimetry, laser-induced fluorescence flow visualization, surface-oil-flow visualization and surface pressure scanning techniques, for L/d = 1.0 ~ 5.0 and Re ≡ U∞d/ν = 2.8 × (103 ~ 104), where U∞ is the free-stream velocity and ν is the kinematic viscosity of the fluid. The flow is classified into four regimes, i.e. the extended-body (L/d ≤ 1.5–2.0), reattachment (1.5–2.0 < L/d < 2.7–3.2), co-shedding (L/d ≥ 3.0–3.4) and transition (2.7 ≤ L/d ≤ 3.3) where both reattachment and co-shedding phenomena may take place. The mean drag and fluctuating drag and lift exhibit distinct features for different flow regimes, which is fully consistent with the proposed flow classification. Comparison is made between this flow and the wake of two tandem circular cylinders, which provides valuable insight into the profound effect of the flow separation point and the presence of sharp corners on the flow development and classification.
The evaporation process in porous media typically experiences three main periods, among which the first period, named the constant rate period (CRP), performs most efficiently in removing liquid. We aim to prolong the CRP to very low degrees of saturation (S) and increase its evaporation rate by playing with heterogeneity in wettability and pore size. First, we show that a porous medium with a smaller contact angle at the surface and increasing contact angle towards the inside generally dries out faster compared with that with uniform contact angle. Second, a constant contact angle porous medium with smaller/larger pores in the surface/inside part dries out faster than a medium with uniform pore size. The underlying mechanism is the occurrence of a capillary pressure jump at the border between the two layers accompanied by enhanced capillary pumping, increasing/maintaining the interfacial area in the surface pores. Harnessing the potential of this mechanism, we propose an optimized strategy by combining two heterogeneity effects: increasing contact angle and pore size towards the inside. This strategy is found to be robust both for multilayer and larger systems. In this case, a small drying front first penetrates fast towards the inside and then expands, followed by a horizontal drying front moving back layer by layer to the surface. Quantitatively, compared with evaporation from a homogeneously porous medium with uniform contact angle where CRP stops at , our optimized design can extend the CRP down to , and decrease five-fold the drying time needed to reach .
In order to improve the predictive abilities of weather and climate models, it is essential to understand the behaviour of wind stress at the ocean surface. Wind stress is contingent on small-scale interfacial dynamics typically not directly resolved in numerical models. Although skin friction contributes considerably to the total stress up to moderate wind speeds, it is notoriously challenging to measure and predict using physics-based approaches. This work proposes a supervised machine learning (ML) model that estimates the spatial distribution of the skin-friction drag over wind waves using solely wave elevation and wave age, which are relatively easy to acquire. The input–output pairs are high-resolution wave profiles and their corresponding surface viscous stresses collected from laboratory experiments. The ML model is built upon a convolutional neural network architecture that incorporates the Mish nonlinearity as its activation function. Results show that the model can accurately predict the overall distribution of viscous stresses; it captures the peak of viscous stress at/near the crest and its dramatic drop to almost null just past the crest in cases of intermittent airflow separation. The predicted area-aggregate skin friction is also in excellent agreement with the corresponding measurements. The proposed method offers a practical pathway for estimating both local and area-aggregate skin friction and can be easily integrated into existing numerical models for the study of air–sea interactions.
This paper investigates the effects of a light-actuated photosurfactant on the canonical problem of the linear stability of a viscous thread surrounded by a dynamically passive fluid. A model consisting of the Navier–Stokes equations and a set of molar concentration equations is presented that capture light-induced switching between two stable surfactant isomer states, trans and cis. These two states display significantly different interfacial properties, allowing for some external control of the stability behaviour of the thread via incident light. Normal modes are used to generate a generalized eigenvalue problem for the growth rate which is solved with a hybrid analytical and numerical method. The results are validated with appropriate analytical solutions of increasing complexity, beginning with a solution to a clean interface, then analytical solutions for one insoluble surfactant, one soluble surfactant and a special case of two photosurfactants with a spatially uniform undisturbed state. Presenting each of these cases allows for a holistic discussion of the effect of surfactants in general on the stability of a liquid thread. Finally, the numerical solutions in the presence of two photosurfactants that display radially non-uniform undisturbed states are presented, and details of the impact of the illumination on the linear stability of the thread are discussed.
We investigate the linear stability analysis of rotating electro-osmotic flow in confined and unconfined configurations by appealing to the Debye–Hückel approximation. Pertaining to flow in confined and unconfined domains, the stability equations are solved using the Galerkin method to obtain the stability picture. Both qualitative and quantitative aspects of Ekman spirals are examined in stable and unstable scenarios within the unconfined domain. Within the confined domain, the variation of the real growth rate and the transition to instability are analysed using the modified Routh–Hurwitz criteria, employed for the first time in this context. The stability of the underlying flow, characterized by the number of roots with a positive real part, is determined by establishing a Routhian table. The inferences of this analysis show that the velocity plane produces intriguing closed Ekman spirals, which diminish in size with an increase in the rotation speed . The Ekman spirals in the stable region exhibit a distinct discontinuity, indicating the dissipation of disturbances over time. In the confined domain, the flow appears consistently stable for a set of involved parameters pertinent to this analysis, such as electrokinetic parameter and rotational parameter approximately up to . However, the flow instabilities become evident for and .
Flow physics vary in different regimes across the full Mach number range, with our knowledge being particularly poor about the hypersonic regime. An Eulerian realization of the particles on demand method, a kinetic model formulated in the comoving reference frame, is proposed to simulate hypersonic compressible flows. The present model allows for flux evaluation in different reference frames, in this case rescaled and shifted by local macroscopic quantities, i.e. fluid speed and temperature. The resulting system of coupled hyperbolic equations is discretized in physical space with a finite volume scheme ensuring exact conservation properties. Regularization via Grad expansion is introduced to implement distribution function and flux transformation between different reference frames. It is shown that the proposed method possesses Galilean invariance at a Mach number up to . Different benchmarks including both inviscid and viscous flows are reproduced with the Mach number up to and pressure ratio up to . Finally, the new model is demonstrated to be capable of simulating hypersonic reactive flows, including one-dimensional and two-dimensional detonations. The developed methodology opens up possibilities for the simulation of the full range of compressible flows, without or with chemical reactions, from the subsonic to hypersonic regimes, leading to enhanced understanding of flow behaviours across the full Mach number range.
The properties of multiphase flows are challenging to measure, and yet effective properties are fundamental to modelling and predicting flow behaviour. The current study is motivated by rheometric measurements of a gas-fluidized bed using a coaxial rheometer in which the fluidization rate and the rotational speed can be varied independently. The measured torque displays a range of rheological states: quasistatic, dense granular flow behaviour at low fluidization rates and low-to-moderate shear rates; turbulent toroidal-vortex flow at high shear rates and moderate-to-high fluidization rates; and viscous-like behaviour with rate-dependent torque at high shear rates and low fluidization or at low shear rates and high fluidization. To understand the solid-like to fluid-like transitions, additional experiments were performed in the same rheometer using single-phase liquid and liquid–solid suspensions. The fluidized bed experiments are modelled as a Bingham plastic for low fluidization rates, and as a shear-thinning Carreau liquid at high fluidization rates. The suspensions are modelled using the Krieger–Dougherty effective viscosity. The results demonstrate that, by using the effective properties, the inverse Bingham number marks the transition from solid-like to viscous-flow behaviour; a modified gap Reynolds number based on the thickness of the shear layer specifies the transition from solid-like to turbulent vortical flow; and a gap Reynolds number distinguishes viscous behaviour from turbulent vortical flow. The results further demonstrate that these different multiphase flows undergo analogous flow transitions at similar Bingham or Reynolds numbers and the corresponding dimensionless torques show comparable scaling in response to annular shear.
There is currently considerable interest in the guided-jet mode, as a result of recent works demonstrating it being the upstream component of various resonant systems in high-speed flows. For given jet operating conditions, the mode is known to exist over only a finite-frequency range that, for a twin-jet system, has been observed to vary with both jet separation and solution symmetry. Vortex-sheet and finite-thickness linear stability models are here employed to consider the behaviour of the guided-jet mode as the two jets are brought together, for both a planar and round twin-jet system. It is demonstrated that in both cases as the twin-jet system merges it forms a higher-order mode of an equivalent single-jet geometry. This then imposes a constraint on the guided-jet mode as the finite-frequency range must change to meet that of the equivalent geometry the system merges to, explaining the previously observed dependence on jet separation.
A spatially developing flat-plate boundary layer free from and two-way coupled with inertial solid particles is simulated to investigate the interaction between particles and the turbulent/non-turbulent interface. Particle Stokes numbers based on the outer scale are (low), 11 (moderate) and 53 (high). The Eulerian–Lagrangian point-particle approach is deployed for the simulation of particle-laden flow. The outer edge of the turbulent/non-turbulent interface layer is detected as an iso-surface of vorticity magnitude. Results show that the particles tend to accumulate below the interface due to the centrifugal effect of large-scale vortices in the outer region of wall turbulence and the combined barrier effect of potential flow. Consequently, the conditionally averaged fluid velocity and vorticity vary more significantly across the interface through momentum exchange and the feedback of force in the enstrophy transport. The large-scale structures in the outer layer of turbulence become smoother and less inclined in particle-laden flow due to the modulation of turbulence by the inertial particles. As a result, the geometric features of the interface layer are changed, namely, the spatial undulation increases, the fractal dimension decreases and the thickness becomes thinner in particle-laden flow as compared with unladen case. These effects become more pronounced as particle inertia increases.
This study investigates the scenario of flow past a circular cylinder in proximity to a moving wall (or equally a body translating in still fluid parallel to a stationary wall). Fifty high-fidelity three-dimensional direct numerical simulations are performed over a parameter space of turbulent Reynolds numbers (Re) of 300 to 1000 combined with gap-to-diameter ratios (G/D) of 0.2 to 3. The flow, hydrodynamic and turbulence characteristics over the (Re, G/D) parameter space are examined in detail. Small-scale rib-like mode B structures and alternate vortex shedding are observed for all cases. The streamwise location for the vortex shedding (quantified by the wake recirculation length Lr) varies strongly with both Re and G/D. The variation of Lr with Re and G/D can be explained by the spanwise circulation fed into the wake, where the variation trends of Lr and are inversely correlated. The variations of the mean drag and fluctuating lift coefficients are also governed by the variations in Lr and . The total kinetic energy in the wake region reduces drastically as G/D reduces below 0.8, which is contributed collectively by (i) reduction in the strength of the shed vortices, (ii) downstream movement of the location of vortex shedding and (iii) associated delayed generation of streamwise vortices. The present results on a moving wall also help to explain several flow and hydrodynamic characteristics reported in the literature for a stationary wall, because the moving wall eliminates the complex wall boundary layer and retains a ‘clean’ near-wall effect.
An impulsively starting motion of two cylindrical bodies floating on a free liquid surface is considered. The shape of the cross-section of each body and the distance between them are arbitrary. The integral hodograph method is advanced to derive the complex velocity potential defined in a rectangular parameter region in terms of the elliptic quasi-doubly periodic Jacobi theta functions. A system of singular integral equations in the velocity magnitude on the free surface and in the slope of the wetted part of each body is derived using the kinematic boundary condition, which is then solved numerically. The velocity field, the pressure impulse on the bodies and the added mass coefficients of each body immediately after the impact are determined in a wide range of distances between the bodies and for cross-sectional shapes such as the flat plate and half-circle.
Direct numerical simulations are utilised to investigate mass-transfer processes at gas-evolving electrodes that experience successive formation and detachment of bubbles. The gas–liquid interface is modelled employing an immersed boundary method. We simulate the growth phase of the bubbles followed by their departure from the electrode surface in order to study the mixing induced by these processes. We find that the growth of the bubbles switches from a diffusion-limited mode at low to moderate fractional bubble coverages of the electrode to a reaction-limited growth dynamics at high coverages. Furthermore, our results indicate that the net transport within the system is governed by the effective buoyancy driving induced by the rising bubbles and that mechanisms commonly subsumed under the term ‘microconvection’ do not significantly affect the mass transport. Consequently, the resulting gas transport for different bubble sizes, current densities and electrode coverages can be collapsed onto one single curve and only depends on an effective Grashof number. The same holds for the mixing of the electrolyte when additionally taking the effect of surface blockage by attached bubbles into account. For the gas transport to the bubble, we find that the relevant Sherwood numbers also collapse onto a single curve when accounting for the driving force of bubble growth, incorporated in an effective Jakob number. Finally, linking the hydrogen transfer rates at the electrode and the bubble interface, an approximate correlation for the gas-evolution efficiency has been established. Taken together, these findings enable us to deduce parametrisations for all response parameters of the systems.
We present evidence for layered anisotropic stratified turbulence (LAST) and mixing produced in a freely evolving uniformly stratified shear layer where the direction of shear is orthogonal to gravity. As originally reported by Basak & Sarkar (J. Fluid. Mech., vol. 568, 2006, pp. 19–54), such a flow develops a rich three-dimensional structure in the form of interlocking columnar vortices formed by horizontal shear instability that remain coherent at large scales due to the stabilising vertical stratification. Here, we modify the initial velocity field by introducing additional small-amplitude vertical perturbations designed to be representative of pre-existing horizontal layers often observed in strongly stratified ocean environments. This reveals a novel finite amplitude, non-normal mode growth mechanism through which the vertical shear between layers may be rapidly amplified by its interaction with the horizontal shear layer prior to the growth of shear instability, leading to a rapid turbulent transition instigated by the subsequent interaction of the layers with the emerging columnar vortices. Through a consideration of relevant flow statistics and associated dimensionless parameters, we demonstrate that turbulence can enter the LAST regime, thereby indicating a generic mechanism leading to the transient development of regions of strongly stratified turbulence in the ocean. We discuss the properties of mixing and the parameterisation of mixing efficiency in terms of the relationship between turbulent length scales in the flow, in particular highlighting links to models based on the classical vertical shear instability paradigm typically associated with more weakly stratified flows that produce isolated turbulent ‘patches’.
This study presents a noise-robust closed-loop control strategy for wake flows employing model predictive control. The proposed control framework involves the autonomous offline selection of hyperparameters, eliminating the need for user interaction. To this purpose, Bayesian optimization maximizes the control performance, adapting to external disturbances, plant model inaccuracies and actuation constraints. The noise robustness of the control is achieved through sensor data smoothing based on local polynomial regression. The plant model can be identified through either theoretical formulation or using existing data-driven techniques. In this work we leverage the latter approach, which requires minimal user intervention. The self-tuned control strategy is applied to the control of the wake of the fluidic pinball, with the plant model based solely on aerodynamic force measurements. The closed-loop actuation results in two distinct control mechanisms: boat tailing for drag reduction and stagnation point control for lift stabilization. The control strategy proves to be highly effective even in realistic noise scenarios, despite relying on a plant model based on a reduced number of sensors.
The familiar process of bubbles generated via breaking waves in the ocean is foundational to many natural and industrial applications. In this process, large pockets of entrained gas are successively fragmented by the ambient turbulence into smaller and smaller bubbles. The key question is how long it takes for the bubbles to reach terminal sizes for a given system. Despite decades of effort, the reported breakup time from multiple experiments differs significantly. Here, to reconcile those results, rather than focusing on one scale, we measure multiple time scales associated with the process through a unique experiment that resolves bubbles’ local deformation and curvature. The results emphasize that the scale separation among various time scales is controlled by the Weber number, similar to how the Reynolds number determines the scale separation in single-phase turbulence, but shows a distinct transition at a critical Weber number.
Fine-scale flow motions are measured in a Newtonian and polymer drag-reduced turbulent boundary layer (TBL) at a common momentum thickness Reynolds number of 2300. Relative to the Newtonian TBL, the polymer-laden flow has a 33 % lower skin-friction coefficient. Three-dimensional (3-D) particle tracking velocimetry is used to measure the components of the velocity gradient tensor (VGT), rate of deformation tensor (RDT) and rate of rotation tensor (RRT). The invariants in these tensors are then used to distinguish the different types of fine-scale flow motions – a method called the -criterion. Joint probability density functions (j.p.d.f.s) of the VGT invariants, and , for the Newtonian TBL produce the familiar tear-drop pattern, commonly seen in direct numerical simulations of Newtonian turbulence. Relative to the Newtonian TBL, the polymer-laden flow has significantly attenuated values of , implying an overall reduction in fluid stretching. The invariants in the RDT, and , imply that straining motions of the polymeric flow are more two dimensional compared with the Newtonian flow. Moreover, j.p.d.f.s of and the invariant in the RRT , suggest that the flow consists of fewer biaxial extensional events and more shear-dominated flow. Few, if any, experimental investigations have measured the 3-D structure of fine-scale motions in a Newtonian and polymer drag-reduced TBL using the -criterion. We provide the first experimental evidence that supports the notion that an attenuation of fluid stretching, particularly biaxial straining motions, is central to the mechanism of polymer drag reduction.
The impact of a liquid droplet with another droplet or onto a solid surface are important basic processes that occur in many applications such as agricultural sprays and inkjet printing, and in nature such as pathogens transport by raindrops. We investigated the head-on collision of unequal-size droplets of the same liquid on wetting surfaces using the direct numerical simulations technique at different size ratios. The unsteady Navier–Stokes equations are solved and the liquid–gas interface is tracked using the geometric volume-of-fluid method. The numerical model is validated by comparing simulation results of two extreme cases of droplets bouncing with the experimental data from previous studies and the agreement is quite accurate. The validated model is employed to simulate droplets bouncing at several size ratios at different Weber numbers and Ohnesorge number. Two distinct regimes are identified, namely, the inertial regime, where the restitution coefficient is a constant value close to 0.3, the viscous regime, where the restitution coefficient declines. To understand the bouncing behaviour, the velocity field is analysed and an energy budget calculation is performed. The distribution of the sessile droplet energy is found to be important and the sessile droplet surface energy is calculated by its deformation characteristics such as crater depth. Finally, a scaling analysis is performed to rationalize the insensitivity of the coefficient of restitution in the inertial regime, and its decline in the viscous regime, at large size ratios.
Large eddy simulation of flow past a cricket ball with its seam at to the free stream is carried out for . Three regimes of flow are identified on the basis of the time-averaged swing force coefficient () – no swing (NS), conventional swing (CS, ) and reverse swing (RS, ). The effect of seam on the boundary layer is investigated. Contrary to the popular belief, the boundary layer does not transition to a turbulent state in the initial stages of CS. The seam energizes the laminar boundary layer and delays its separation. The delay is significantly larger in a region near the poles, whose extent increases with an increase in causing to increase. Here assumes a near constant value in the later stage of CS. The boundary layer transitions to a turbulent state via formation of a laminar separation bubble (LSB) in the equatorial region and directly, without a LSB, in the polar region. The extent of the LSB shrinks while the region of direct transition near the poles increases with an increase in . A LSB forms on the non-seam side of the ball in the RS regime. A secondary vortex is observed in the wake bubble. While it exists on the non-seam side for the entire range of considered, the mixing in the flow introduced by the seam causes it to disappear beyond a certain on the seam side. The pressure difference between the seam and non-seam sides sets up wing-tip-like vortices. Their polarity reverses with the switch from the CS to RS regime.
The cubic interactions in a discrete system of four weakly nonlinear waves propagating in a conservative dispersive medium are studied. By reducing the problem to a single ordinary differential equation governing the motion of a classical particle in a quartic potential, the complete explicit branches of solutions are presented, either steady, periodic, breather or pump, thus recovering or generalizing some already published results in hydrodynamics, nonlinear optics and plasma physics, and presenting some new ones. Various stability criteria are also formulated for steady equilibria. Theory is applied to deep-water gravity waves for which models of isolated quartets are described, including bidirectional standing waves and quadri-directional travelling waves, steady or not, resonant or not.
This paper examines the rigid body motion of a spheroid sedimenting in a Newtonian fluid with a spatially varying viscosity field. The fluid is at zero Reynolds number, and the viscosity varies linearly in space in an arbitrary direction with respect to the external force. First, we obtain the correction to the spheroid's rigid body motion in the limit of small viscosity gradients, using a perturbation expansion combined with the reciprocal theorem. Next, we determine the general form of the particle's mobility tensor relating its rigid body motion to an external force and torque. The viscosity gradient does not alter the force/translation and torque/rotation relationships, but introduces new force/rotation and torque/translation couplings that are determined for a wide range of particle aspect ratios. Finally, we discuss results for the spheroid's rotation and centre-of-mass trajectory during sedimentation. A steady orientation arises at long time whose value depends on the viscosity gradient direction and particle shape. These results are significantly different than when no viscosity gradient is present, where the particle stays at its initial orientation for all times. We summarize the observations for prolate and oblate spheroids for different viscosity gradient directions and provide plots for the orientation and centre-of-mass trajectory versus time. We also provide guidelines to extend the analysis when the viscosity gradient exhibits a more complicated spatial behaviour.
The aeroacoustic feedback loops in high-speed circular jets that impinge on a large flat plate are investigated via acoustic measurements and schlieren visualizations. In the present experiments, the nozzle pressure ratio ranges from 1.39 to 2.20, the corresponding ideally expanded jet Mach number is from 0.70 to 1.12 and the nozzle-to-plate distance () is from 4.0 to 6.0, where is the nozzle exit diameter. The results of acoustic measurements show that the strongest tones are generated in a limited frequency band. The empirical dispersion relations obtained from the fluctuating greyscales along the jet centreline of time-resolved schlieren images have good agreement with the dispersion relations from the vortex-sheet model. The coherent flow structures at tonal frequencies are extracted by spectral proper orthogonal decomposition and are analysed in detail. For the jets, the upstream-propagating guided jet mode is progressively confined to the potential core of jets with increasing tonal frequency, which provides the first direct experimental support for theoretical results. The evolution in the structures of acoustic resonance loops is studied along a single frequency stage of axisymmetric impinging tones. When the acoustic resonance between the upstream- and downstream-propagating guided jet modes is formed at tonal frequencies, the impinging tones are intenser. Slightly underexpanded impinging jets can simultaneously produce impingement tones and screech tones. Shock-cell structures have modulatory effects on the downstream-propagating Kelvin–Helmholtz wavepacket and the upstream- and downstream-propagating guided jet modes. Due to the interaction between the flow structures at the frequencies of impinging and screech tones, tones of axisymmetric modes can be produced outside the frequency ranges in which the axisymmetric upstream-propagating guided jet modes are supported by jets.