Linear stability and sensitivity of the flow past a fixed oblate spheroidal bubble

The stability properties of the wake past an oblate spheroidal bubble held fixed in a uniform stream are studied in the framework of a global linear analysis. In line with previous studies, provided the geometric aspect ratio of the bubble, χ, is large enough, the wake is found to be unstable only within a finite range of Reynolds number, Re. The neutral curves corresponding to the occurrence of the first two unstable modes are determined over a wide range of the (χ, Re) domain and the structure of the modes encountered along the two branches of each neutral curve is discussed. Then, using an adjoint-based approach, a series of sensitivity analyses of the flow past the bubble is carried out in the spirit of recent studies devoted to twodimensionaland axisymmetric rigid bodies. The regions of the flow most sensitiveto an external forcing are found to be concentrated in the core or at the periphery of the standing eddy, as already observed with bluff bodies at the surface of which the flow obeys a no-slip condition. However, since the shear-free condition allows the fluid to slip along the bubble surface, the rear half of this surface turns out to be also significantly sensitive to disturbances originating in the shear stress, a finding which may be related to the well-known influence of surfactants on the structure and stability properties of the flow past bubbles rising in water

Bifurcations and symmetry breaking in the wake of axisymmetric bodies

We consider the generic problem of wake instabilities past fixed axisymmetric bodies, and focus on the extreme cases of a sphere and a flat disk. Numerical results reveal that the wakes of these two bodies evolve differently as the Reynolds number is increased. Especially, two new vortex shedding modes are identified behind a disk. To interpret these results, we introduce a model based on the theory of mode interactions in presence of O(2) symmetry. This model, which was initially developed for the Taylor–Couette system, allows us to explain the structural differences observed in the evolution of the two types of wakes and to accurately predict the evolution of the lift force

Falling styles of disks

We numerically investigate the dynamics of thin disks falling under gravity in a viscous fluid medium at rest at infinity. Varying independently the density and thickness of the disk reveals the influence of the disk aspect ratio which, contrary to previous belief, is found to be highly significant as it may completely change the route to non-vertical paths as well as the boundaries between the various path regimes. The transition from the straight vertical path to the planar fluttering regime is found to exhibit complex dynamics: a bistable behaviour of the system is detected within some parameter range and several intermediate regimes are observed in which, although the wake is unstable, the path barely deviates from vertical. By varying independently the body-to-fluid inertia ratio and the relative magnitude of inertial and viscous effects over a significant range, we set up a comprehensive map of the corresponding styles of path followed by an infinitely thin disk. We observe the four types of planar regimes already reported in experiments but also identify two additional fully three-dimensional regimes in which the body experiences a slow horizontal precession superimposed onto zigzagging or tumbling motions

The onset of unsteadiness of two-dimensional bodies falling or rising freely in a viscous fluid: a linear study

We consider the transition between the steady vertical path and the oscillatory path of two-dimensional bodies moving under the effect of buoyancy in a viscous fluid. Linearization of the Navier–Stokes equations governing the flow past the body and of Newton’s equations governing the body dynamics leads to an eigenvalue problem, which is solved numerically. Three different body geometries are then examined in detail, namely a quasi-infinitely thin plate, a plate of rectangular cross-section with an aspect ratio of 8, and a rod with a square cross-section. Two kinds of eigenmodes are observed in the limit of large body-to-fluid mass ratios, namely ‘fluid’ modes identical to those found in the wake of a fixed body, which are responsible for the onset of vortex shedding, and four additional ‘aerodynamic’ modes associated with much longer time scales, which are also predicted using a quasi-static model introduced in a companion paper. The stability thresholds are computed and the nature of the corresponding eigenmodes is investigated throughout the whole possible range of mass ratios. For thin bodies such as a flat plate, the Reynolds number characterizing the threshold of the first instability and the associated Strouhal number are observed to be comparable with those of the corresponding fixed body. Other modes are found to become unstable at larger Reynolds numbers, and complicated branch crossings leading to mode switching are observed. On the other hand, for bluff bodies such as a square rod, two unstable modes are detected in the range of Reynolds number corresponding to wake destabilization. For large enough mass ratios, the leading mode is similar to the vortex shedding mode past a fixed body, while for smaller mass ratios it is of a different nature, with a Strouhal number about half that of the vortex shedding mode and a stronger coupling with the body dynamics

Waves and instabilities in rotating free surface flows

The stability properties of the rotating free surface flow in a cylindrical container is studied using a global stability approach, considering succesively three models. For the case of solid body rotation (Newton’s bucket), all eigenmodes are found to be stable, and are classified into three families : gravity waves, singular inertial modes, and Rossby waves. For the case of a potential flow, an instability is found. The mechanism is explained as a resonance between gravity waves and centrifugal waves, and is thought to be at the origin of the ”rotating polygon instability” observed in experiments where the flow is driven by rotation of the bottom plate (see [9]). Finally, we give some preliminary results concerning a third model : the Rankine vortex

Linear and nonlinear dynamics of axisymmetric waves in the hollow core vortex

The dynamics of trailing vortices are under constant investigation during last decades since it is of considerable interest to reduce aircraft wakes and associated hazards to forthcoming planes. The isolated axisymmetric vortex is the commonly used simplest elementary model when considering such issue. Although asymptotically stable, recent studies have revealed its sensitiveness to specific perturbations, leading in some cases to considerable gains of energy. Albeit of evident interest, the underlying mechanisms of energy growth are considered in the linear regime. The nonlinear dynamics of such vortices need also to be considered in order to complete the picture. Rather than performing direct numerical simulations3, an interesting way to investigate it is to consider the nonlinear interactions of waves. This approach is motivated by the possible existence of resonance between wave components. For this purpose, the base flow model is simplified by considering the hollow core vortex. Arising naturally when a tank is drained (bath-tube vortex), it presents simpler dynamics than the Lamb-Oseen vortex as it only possesses two families of waves. This point is of crucial importance for the tractability of the problem. In this work, the nonlinear temporal evolution of axisymmetric waves are investigated through numerical integration when the flow is submitted to various initial conditions (travelling or standing wave, pinching of the free surface, wave trains). We focus on wave trains as important energy exchanges between the main component and its sideband waves are observed. This phenomenon is related to the Benjamin-Feir instability4 (triadic resonance) occurring for wave trains on deep water

The steady oblique path of buoyancy-driven disks and spheres

We consider the steady motion of disks of various thicknesses in a weakly viscous flow, in the case where the angle of incidence $\alpha$ (defined as that between the disk axis and its velocity) is small. We derive the structure of the steady flow past the body and the associated hydrodynamic force and torque through a weakly nonlinear expansion of the flow with respect to $\alpha$. When buoyancy drives the body motion, we obtain a solution corresponding to an oblique path with a non-zero incidence by requiring the torque to vanish and the hydrodynamic and net buoyancy forces to balance each other. This oblique solution is shown to arise through a bifurcation at a critical Reynolds number ReSO which does not depend upon the body-to-fluid density ratio and is distinct from the critical Reynolds number ReSS corresponding to the steady bifurcation of the flow past the body held fixed with $\alpha=0$. We then apply the same approach to the related problem of a sphere that weakly rotates about an axis perpendicular to its path and show that an oblique path sets in at a critical Reynolds number ReSO slightly lower than ReSS , in agreement with available numerical studies

Stochastic forcing of the Lamb–Oseen vortex

The aim of the present paper is to analyse the dynamics of the Lamb–Oseen vortex when continuously forced by a random excitation. Stochastic forcing is classically used to mimic external perturbations in realistic configurations, such as variations of atmospheric conditions, weak compressibility effects, wing-generated turbulence injected in aircraft wake, or free-stream turbulence in wind tunnel experiments. The linear response of the Lamb–Oseen vortex to stochastic forcing can be decomposed in relation to the azimuthal symmetry of the perturbation given by the azimuthal wavenumber m. In the axisymmetric case m = 0, we find that the response is characterised by the generation of vortex rings at the outer periphery of the vortex core. This result is consistent with recurrent observations of such dynamics in the study of vortex-turbulence interaction. When considering helical perturbations m = 1, the response at large axial wavelengths consists of a global translation of the vortex, a feature very similar to the phenomenon of vortex meandering (or wandering) observed experimentally, corresponding to an erratic displacement of the vortex core. At smaller wavelengths, we find that stochastic forcing can excite specific oscillating modes of the Lamb–Oseen vortex. More precisely, damped critical-layer modes can emerge via a resonance mechanism. For perturbations with higher azimuthal wavenumber m > 2, we find no structure that clearly dominates the response of the vortex

A quasi-static approach to the stability of the path of heavy bodies falling within a viscous fluid

We consider the gravity-driven motion of a heavy two-dimensional rigid body freely falling in a viscous fluid.We introduce a quasi-static linear model of the forces and torques induced by the possible changes in the body velocity,or by the occurrence of a nonzero incidence angle or a spanwise rotation of the body. The coefficients involved in this model are accurately computed over a full range of Reynolds number by numerically resolving the Navier 13Stokes equations, considering three elementary situations where the motion of the body is prescribed. The falling body is found to exhibit three distinct eigenmodes which are always damped in the case of a thin plate with uniform mass loading or a circular cylinder,but may be amplified for other geometries,such as in the case of a square cylinder

Linear stability of disks falling or rising freely in a viscous fluid

The problem of the instability of a solid body moving in a viscous fluid has been extensively studied in the case where the body is fixed and its wake undergoes instabilities when some control parameter is increased (e:g: 1;2 for two-dimensional configurations and 3;4 for three-dimensional bodies). On the other hand, the deeply different and frequent case of a body freely falling or rising under buoyancy in a fluid otherwise at rest hasreceived much less attention. Yet, results from experiments5;6 and direct numerical simulations7;8 have evidenced a large variety of motion styles ranging from steady oblique to chaotic, though zig-zag and tumbling. To gain more insight into the nature of the instability that drives the departure from a straight vertical path, we recently carried out a linear study of the coupled fluid-body problem in two dimensions9. Here we extend this work to axisymmetric bodies. The linear stability analysis is performed for a disk of finite thickness with an aspect ratio (defined as the ratio of the diameter over thickness) X = 3 and another disk with X = 10. We select these two values because past studies10 led us to consider the corresponding two bodies as prototypes of general thick and thin bodies. This talk mainly focuses on the parametric modal stability analysis of the coupled fluid-disk problem. We present neutral stability curves in the phase space (pb=pf ,Re) and (pb=pf ; St), where pb=pf is the body-to-fluid density ratio, St the non-dimensional frequency of the unstable mode and Re the Reynolds number based on the disk diameter and the relative velocity between the fluid and body. These curves reveal how rich the dynamics of the problem are, including features such as destabilization-restabilization regions and abrupt jumps of the marginal frequency. The types and thresholds of the linear instability show close agreement with existing DNS and experimental results, as well as with theoretical predictions from asymptotic theories11;12. Finally, we consider the structure of the global modes. This allows us to extract information such as the origin of the instability (i.e. we disentangle \body-related modes" from \ fluid-related modes") and the characteristic footprints of both low-frequency and high-frequency modes