Showing 1–50 of 53 results for author: Gallaire, F

Lagrangian vs Eulerian view on the mean drift and streaming flows in orbital sloshing

Authors: Alessandro Bongarzone, François Gallaire

Abstract: Orbital sloshing is a technique to gently mix a container's liquid content and it is commonly used in fermentation and cell cultivation processes. Besides the rich wave dynamics observed at the interface, Bouvard et al. (2017) [1] unveiled the structure of the Lagrangian mean flow hiding in the fluid bulk. The latter flow shows a global toroidal (azimuthal) rotation co-directed with the wave and n… ▽ More Orbital sloshing is a technique to gently mix a container's liquid content and it is commonly used in fermentation and cell cultivation processes. Besides the rich wave dynamics observed at the interface, Bouvard et al. (2017) [1] unveiled the structure of the Lagrangian mean flow hiding in the fluid bulk. The latter flow shows a global toroidal (azimuthal) rotation co-directed with the wave and nontrivial poloidal vortices near the contact line. Rotating sloshing waves are known to induce a net motion of fluid particles and hence a wave-averaged difference between the Eulerian flow - viscous streaming - and the Lagrangian flow, that is commonly referred to as Stokes drift. Nevertheless, discerning these two components in an experiment is challenging as they scale similarly with the forcing amplitude and frequency. Their relative contributions remain therefore unquantified, particularly for highly viscous fluids, for which prior analysis, based on inviscid arguments, fails. In this work, we construct a truncated asymptotic approximation of the problem, where the solution at each order is computed numerically to describe as accurately as possible the contact line region, otherwise analytically untractable. The results of this weakly nonlinear analysis in terms of first-order wave and second-order mean flow are then thoroughly compared with the experiments by Bouvard et al. (2017) [1], showing a remarkable agreement for off-resonance frequencies. When viscosity matters, our findings suggest that it is incorrect to attribute the poloidal patterns solely to the Eulerian streaming flow and that viscous corrections to the Stokes drift are equally important in the resulting mean Lagrangian flow. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Tumbling elimination induced by permeability: an experimental approach

Authors: Jesús Sánchez-Rodríguez, François Gallaire

Abstract: Archetypal falling behaviors of impervious objects are classified into four modes: fluttering, tumbling, steady descent and chaotic motion. The classical scenario predicts these behaviors to be affected by two dimensionless quantities: dimensionless inertia and Reynolds number. In this article we explore experimentally the effect of permeability and porosity on the falling regimes of porous plates… ▽ More Archetypal falling behaviors of impervious objects are classified into four modes: fluttering, tumbling, steady descent and chaotic motion. The classical scenario predicts these behaviors to be affected by two dimensionless quantities: dimensionless inertia and Reynolds number. In this article we explore experimentally the effect of permeability and porosity on the falling regimes of porous plates. By drilling several hole distributions in rectangular plates, both permeability and porosity are varied systematically. We discover that the introduction of porosity affects the stability of the falling regimes eliminating tumbling. Using a phenomenological model we show that a decrease in circulation induced by the introduction of holes is the primary mechanism for stabilizing the plates' trajectories. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Coating thickness prediction for a viscous film on a rough plate

Authors: Lebo Molefe, Giuseppe A. Zampogna, John M. Kolinski, François Gallaire

Abstract: Surface roughness significantly modifies the liquid film thickness entrained when dip coating a solid surface, particularly at low coating velocity. Using a homogenization approach, we present a predictive model for determining the liquid film thickness coated on a rough plate. A homogenized boundary condition at an equivalent flat surface is used to model the rough boundary, accounting for both f… ▽ More Surface roughness significantly modifies the liquid film thickness entrained when dip coating a solid surface, particularly at low coating velocity. Using a homogenization approach, we present a predictive model for determining the liquid film thickness coated on a rough plate. A homogenized boundary condition at an equivalent flat surface is used to model the rough boundary, accounting for both flow through the rough texture layer, through an interface permeability term, and slip at the equivalent surface. While the slip term accounts for tangential velocity induced by viscous shear stress, accurately predicting the film thickness requires the interface permeability term to account for additional tangential flow driven by pressure gradients along the interface. We find that a greater degree of slip and interface permeability signifies less viscous stress that would promote deposition, thus reducing the amount of free film coated above the textures. The model is found to be in good agreement with experimental measurements and requires no fitting parameters. Furthermore, our model may be applied to arbitrary periodic roughness patterns, opening the door to flexible characterization of surfaces found in natural and industrial coating processes. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Comments: 29 pages, 12 figures

Releasing trapped Taylor bubbles via centrifugation and inclination

Authors: Alice Marcotte, Pier Giuseppe Ledda, Valentin Buriasco, Paul Dené, François Gallaire, Ludovic Keiser

Abstract: In confined systems, the entrapment of a gas volume with an equivalent spherical diameter greater than the dimension of the channel can form extended bubbles that obstruct fluid circuits and compromise performance. Notably, in sealed vertical tubes, buoyant long bubbles -- called Taylor bubbles -- cannot rise if the inner tube radius is below a critical value near the capillary length. This critic… ▽ More In confined systems, the entrapment of a gas volume with an equivalent spherical diameter greater than the dimension of the channel can form extended bubbles that obstruct fluid circuits and compromise performance. Notably, in sealed vertical tubes, buoyant long bubbles -- called Taylor bubbles -- cannot rise if the inner tube radius is below a critical value near the capillary length. This critical threshold for steady ascent is determined by geometric constraints related to matching the upper cap shape with the lubricating film in the elongated part of the bubble. Developing strategies to overcome this threshold and release stuck bubbles is essential for applications involving narrow liquid channels. Effective strategies involve modifying matching conditions with an external force field to facilitate bubble ascent. However, it's unclear how changes in acceleration conditions affect the motion onset of buoyancy-driven long bubbles. This study investigates the mobility of elongated bubbles in sealed tubes with an inner radius near the critical value inhibiting bubble motion in a vertical setting. Two strategies are explored to tune bubble motion, leveraging variations in axial and transversal accelerations: tube rotation around its axis and tube inclination relative to gravity. By revising the geometrical constraints of the simple vertical setting, the study predicts new thresholds based on rotational speed and tilt angle, respectively, providing forecasts for the bubble rising velocity under modified apparent gravity. Experimental measurements of motion threshold and rising velocity compare well with theoretical developments, thus suggesting practical approaches to control and tune bubble motion in confined environments. △ Less

Submitted 27 April, 2024; originally announced April 2024.

Noise-induced transitions past the onset of a steady symmetry-breaking bifurcation: the case of the sudden expansion

Authors: Yves-Marie Ducimetière, Edouard Boujo, François Gallaire

Abstract: We consider fluid flows, governed by the Navier-Stokes equations, subject to a steady symmetry-breaking bifurcation and forced by a weak noise acting on a slow time scale. By generalizing the multiple-scale weakly nonlinear expansion technique employed in the literature for the response of the Duffing oscillator, we rigorously derive a stochastically forced Stuart-Landau equation for the dominant… ▽ More We consider fluid flows, governed by the Navier-Stokes equations, subject to a steady symmetry-breaking bifurcation and forced by a weak noise acting on a slow time scale. By generalizing the multiple-scale weakly nonlinear expansion technique employed in the literature for the response of the Duffing oscillator, we rigorously derive a stochastically forced Stuart-Landau equation for the dominant symmetry-breaking mode. The probability density function of the solution, and of the escape time from one attractor to the other, are then determined by solving the associated Fokker-Planck equation. The validity of this reduced order model is tested on the flow past a sudden expansion, for a given Reynolds number and different noise amplitudes. At a very low numerical cost, the statistics obtained from the amplitude equation accurately reproduce those of long-time direct numerical simulations. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Stick-slip-to-stick transition of liquid oscillation in a U-tube

Authors: Alessandro Bongarzone, François Gallaire

Abstract: The nonlinear decay of oscillations of a liquid column in a U-shaped tube is investigated within the theoretical framework of the projection method formalized by Bongarzone et al. (2021) [1]. Starting from the full hydrodynamic system supplemented by a phenomenological contact line model, this physics-inspired method uses successive linear eigenmode projections to simulate the relaxation dynamics… ▽ More The nonlinear decay of oscillations of a liquid column in a U-shaped tube is investigated within the theoretical framework of the projection method formalized by Bongarzone et al. (2021) [1]. Starting from the full hydrodynamic system supplemented by a phenomenological contact line model, this physics-inspired method uses successive linear eigenmode projections to simulate the relaxation dynamics of liquid oscillations in the presence of sliding triple lines. Each projection is shown to eventually induce a rapid loss of total energy in the liquid motion, thus contributing to its nonlinear damping. A thorough quantitative comparison with experiments by Dollet et al. (2020) [2] demonstrates that, in contradistinction with their simplistic one-degree-of-freedom model, the present approach not only describes well the transient stick-slip dynamics, but it also correctly captures the global stick-slip to stick transition, as well as the secondary bulk motion following the arrest of the contact line, which has been so far overlooked by existing theoretical analyses. This study offers a further contribution to rationalizing the impact of contact angle hysteresis and its associated solidlike friction on the decay of liquid oscillations in the presence of sliding triple lines. △ Less

Submitted 16 November, 2023; originally announced November 2023.

A revised gap-averaged Floquet analysis of Faraday waves in Hele-Shaw cells

Authors: Alessandro Bongarzone, Baptiste Jouron, Francesco Viola, François Gallaire

Abstract: Existing theoretical analyses of Faraday waves in Hele-Shaw cells rely on the Darcy approximation and assume a parabolic flow profile in the narrow direction. However, Darcy's model is known to be inaccurate when convective or unsteady inertial effects are important. In this work, we propose a gap-averaged Floquet theory accounting for inertial effects induced by the unsteady terms in the Navier-S… ▽ More Existing theoretical analyses of Faraday waves in Hele-Shaw cells rely on the Darcy approximation and assume a parabolic flow profile in the narrow direction. However, Darcy's model is known to be inaccurate when convective or unsteady inertial effects are important. In this work, we propose a gap-averaged Floquet theory accounting for inertial effects induced by the unsteady terms in the Navier-Stokes equations, a scenario that corresponds to a pulsatile flow where the fluid motion reduces to a two-dimensional oscillating Poiseuille flow, similarly to the Womersley flow in arteries. When gap-averaging the linearized Navier-Stokes equation, this results in a modified damping coefficient, which is a function of the ratio between the Stokes boundary layer thickness and the cell's gap, and whose complex value depends on the frequency of the wave response specific to each unstable parametric region. We first revisit the standard case of horizontally infinite rectangular Hele-Shaw cells by also accounting for a dynamic contact angle model. A comparison with existing experiments shows the predictive improvement brought by the present theory and points out how the standard gap-averaged model often underestimates the Faraday threshold. The analysis is then extended to the less conventional case of thin annuli. A series of dedicated experiments for this configuration highlights how Darcy's thin-gap approximation overlooks a frequency detuning that is essential to correctly predict the locations of the Faraday tongues in the frequency-amplitude parameter plane. These findings are well rationalized and captured by the present model. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Stability of liquid film coating a horizontal cylinder: interplay of capillary and gravity forces

Authors: Shahab Eghbali, Simeon Djambov, François Gallaire

Abstract: We study the drainage of a viscous liquid film coating outside a horizontal cylinder. We first study the evolution of the axially invariant draining flow, initiated at rest with uniform film thickness $δ$. Non-linear simulations indicate that for each $δ$, there is a threshold in the Bond number ($Bo$), which compares the gravitational effects with surface tension, above which the draining liquid… ▽ More We study the drainage of a viscous liquid film coating outside a horizontal cylinder. We first study the evolution of the axially invariant draining flow, initiated at rest with uniform film thickness $δ$. Non-linear simulations indicate that for each $δ$, there is a threshold in the Bond number ($Bo$), which compares the gravitational effects with surface tension, above which the draining liquid bulk ruptures. This critical $Bo$ is found to scale inversely with $δ$, defining the existence of a quasi-stationary pendant curtain sustained below the cylinder by surface tension. The interface of the pendant curtain is unconditionally linearly unstable and is prone to Rayleigh-Plateau-like, capillarity-driven, and Rayleigh-Taylor, gravity-driven, instabilities. The linear stability of the quasi-static state along with an energy analysis of the unstable mode illustrates that while the Rayleigh-Taylor instability is always present, capillary effects dominate the instability at small $Bo$, which promotes the formation of pearls enveloping the cylinder. In contrast, at large $Bo$, capillarity acts in a stabilising way and the instability is purely gravity-driven, forming underside modulations. We present the asymptotic energy repartition representing the different physical mechanisms at play in the instability of the saturated curtains for a wide range of $\{Bo,δ\}$. The results of the linear analysis agree with the pre-existing experiments of~\citeauthor{Bruyn1997} and non-linear simulations of~\citeauthor{Weidner1997} on thin film and extend the results for thick films. Additionally, based on the volume made available for droplet growth by the development of the most linearly amplified wavelength, we build a tentative regime diagram predicting the final patterns emerging from the pendant curtain: an array of pearls or pendant drops or three-dimensional droplet pinch-off. △ Less

Submitted 15 May, 2023; originally announced May 2023.

Swirling against the forcing: evidence of stable counter-directed sloshing waves in orbital-shaken reservoirs

Authors: Alice Marcotte, François Gallaire, Alessandro Bongarzone

Abstract: We study the free surface response in a cylindrical container undergoing an elliptic periodic orbit. For small forcing amplitudes and deep liquid layers, we quantify the effect of orbit's aspect ratio onto the surface dynamics in the vicinity of the fluid system's lowest natural frequency. We provide experimental evidences of the existence of a frequency range where stable swirling can be either c… ▽ More We study the free surface response in a cylindrical container undergoing an elliptic periodic orbit. For small forcing amplitudes and deep liquid layers, we quantify the effect of orbit's aspect ratio onto the surface dynamics in the vicinity of the fluid system's lowest natural frequency. We provide experimental evidences of the existence of a frequency range where stable swirling can be either co- or counter-directed with respect to the container's direction of motion. Our findings are successfully predicted by an inviscid asymptotic model, amended with a heuristic damping. △ Less

Submitted 3 March, 2023; v1 submitted 28 February, 2023; originally announced February 2023.

Super-harmonically resonant swirling waves in longitudinally forced circular cylinders

Authors: Alice Marcotte, François Gallaire, Alessandro Bongarzone

Abstract: Resonant sloshing in circular cylinders was studied by Faltinsen et al. (2016), whose theory was used to describe steady-state resonant waves due to time-harmonic container's elliptic orbits. In the limit of longitudinal container motions, a symmetry-breaking of the planar wave solution occurs, with clockwise and anti-clockwise swirling equally likely. In addition to this primary harmonic dynamics… ▽ More Resonant sloshing in circular cylinders was studied by Faltinsen et al. (2016), whose theory was used to describe steady-state resonant waves due to time-harmonic container's elliptic orbits. In the limit of longitudinal container motions, a symmetry-breaking of the planar wave solution occurs, with clockwise and anti-clockwise swirling equally likely. In addition to this primary harmonic dynamics, previous experiments have unveiled that diverse super-harmonic dynamics are observable far from primary resonances. Among these, the so-called double-crest (DC) dynamics, first observed by Reclari et al. (2014) for rotary sloshing, is particularly relevant, as its manifestation is the most favored by the spatial structure of the external driving. Following Bongarzone et al. (2022a), in this work we develop a weakly nonlinear (WNL) analysis to describe the system response to super-harmonic longitudinal forcing. The resulting system of amplitude equations predicts that a planar wave symmetry-breaking via stable swirling may also occur under super-harmonic excitation. This finding is confirmed by our experimental observations, which identify three possible super-harmonic regimes, i.e. (i) stable planar DC waves, (ii) irregular motion and (iii) stable swirling DC waves, whose corresponding stability boundaries in the forcing frequency-amplitude plane quantitatively match the present theoretical estimates. △ Less

Submitted 5 February, 2023; originally announced February 2023.

On the liquid film instability of an internally coated horizontal tube

Authors: Shahab Eghbali, Yves-Marie Ducimetiere, Edouard Boujo, Francois Gallaire

Abstract: We study numerically and theoretically the gravity-driven flow of a viscous liquid film coating the inner side of a horizontal cylindrical tube and surrounding a shear-free dynamically inert gaseous core. The liquid-gas interface is prone to the Rayleigh-Plateau and Rayleigh-Taylor instabilities. Here, we focus on the limit of low and intermediate Bond numbers, Bo, where the capillary and gravitat… ▽ More We study numerically and theoretically the gravity-driven flow of a viscous liquid film coating the inner side of a horizontal cylindrical tube and surrounding a shear-free dynamically inert gaseous core. The liquid-gas interface is prone to the Rayleigh-Plateau and Rayleigh-Taylor instabilities. Here, we focus on the limit of low and intermediate Bond numbers, Bo, where the capillary and gravitational forces are comparable and the Rayleigh-Taylor instability is known to be suppressed. We first study the evolution of the axially invariant draining flow, initiating from a uniform film thickness until reaching a quasi-static regime as the bubble approaches the upper tube wall. We then investigate the flow linear stability within two frameworks: frozen time-frame (quasi-steady) stability analysis and transient growth analysis. We explore the effect of the surface tension (Bo) and inertia (measured by the Ohnesorge number, Oh) on the flow and its stability. The linear stability analysis suggests that the interface deformation at large Bo results in the suppression of the Rayleigh-Plateau instability in the asymptotic long-time limit. Furthermore, the transient growth analysis suggests that the initial flow evolution does not lead to any considerable additional amplification of initial interface perturbations, a posteriori rationalising the quasi-steady assumption. The present study yields a satisfactory prediction of the stabilisation threshold found experimentally by Duclaux et al. [15]. △ Less

Submitted 21 November, 2022; originally announced November 2022.

Numerical estimate of the viscous damping of capillary-gravity waves: A macroscopic depth-dependent slip-length model

Authors: Alessandro Bongarzone, Francois Gallaire

Abstract: We propose a numerical approach to regularize the contact line singularity appearing in the computation of viscous capillary-gravity waves with moving contact line in cylindrical containers. The linearized Navier-Stokes equations are complemented by a macroscopic Navier-like slip condition on the container side wall, with a depth-varying slip-length bridging a free-slip condition at the contact li… ▽ More We propose a numerical approach to regularize the contact line singularity appearing in the computation of viscous capillary-gravity waves with moving contact line in cylindrical containers. The linearized Navier-Stokes equations are complemented by a macroscopic Navier-like slip condition on the container side wall, with a depth-varying slip-length bridging a free-slip condition at the contact line to a no-slip condition further away from the meniscus. In accordance with suggestions from the literature, this characteristic penetration depth is chosen as the typical Stokes layer thickness pertaining to the eigenfrequency. Since the latter is unknown, the resulting nonlinear eigenvalue problem is first simplified using the inviscid eigenfrequency to estimate the Stokes layer thickness. The solution is then shown to provide consistent results when compared to the asymptotic approaches of the literature using the inviscid eigenmodes to estimate the different sources of dissipation with the exception of the neglected corner regions. This demonstrates that such numerical regularization scheme is suitable to retrieve physically meaningful results. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Breaking one into three: surface-tension-driven droplet breakup in T-junctions

Authors: Jiande Zhou, Yves-Marie Ducimetière, Daniel Migliozzi, Ludovic Keiser, Arnaud Bertsch, François Gallaire, Philippe Renaud

Abstract: Droplet breakup is an important phenomenon in the field of microfluidics to generate daughter droplets. In this work, a novel breakup regime in the widely studied T-junction geometry is reported, where the pinch-off occurs laterally in the two outlet channels, leading to the formation of three daughter droplets, rather than at the center of the junction for conventional T-junctions which leads to… ▽ More Droplet breakup is an important phenomenon in the field of microfluidics to generate daughter droplets. In this work, a novel breakup regime in the widely studied T-junction geometry is reported, where the pinch-off occurs laterally in the two outlet channels, leading to the formation of three daughter droplets, rather than at the center of the junction for conventional T-junctions which leads to two daughter droplets. It is demonstrated that this new mechanism is driven by surface tension, and a design rule for the T-junction geometry is proposed. A model for low values of the capillary number $Ca$ is developed to predict the formation and growth of an underlying carrier fluid pocket that accounts for this lateral breakup mechanism. At higher values of $Ca$, the conventional regime of central breakup becomes dominant again. The competition between the new and the conventional regime is explored. Altogether, this novel droplet formation method at T-junction provides the functionality of alternating droplet size and composition, which can be important for the design of new microfluidic tools. △ Less

Submitted 30 June, 2022; originally announced July 2022.

An amplitude equation modeling the single-to-double crest wave transition in orbital shaken cylindrical containers

Authors: Alessandro Bongarzone, Margherita Guido, François Gallaire

Abstract: The container motion along a planar circular trajectory at a constant angular velocity, i.e. orbital shaking, is of interest in several industrial applications, e.g. for fermentation processes or in cultivation of stem cells, where good mixing and efficient gas exchange are the main targets. Under these external forcing conditions, the free surface typically exhibits a primary steady state motion… ▽ More The container motion along a planar circular trajectory at a constant angular velocity, i.e. orbital shaking, is of interest in several industrial applications, e.g. for fermentation processes or in cultivation of stem cells, where good mixing and efficient gas exchange are the main targets. Under these external forcing conditions, the free surface typically exhibits a primary steady state motion through a single-crest dynamics, whose wave amplitude, as a function of the external forcing parameters, shows a Duffing-like behaviour. However, previous experiments in lab-scale cylindrical containers have unveiled that, owing to the excitation of super-harmonics, diverse dynamics are observable in certain driving-frequency ranges. Among these super-harmonics, the double-crest dynamics is particularly relevant, as it displays a notably large amplitude response, that is strongly favored by the spatial structure of the external forcing. In the inviscid limit and with regards to circular cylindrical containers, we formalize here a weakly nonlinear analysis via multiple timescale method of the full hydrodynamic sloshing system, leading to an amplitude equation suitable to describe such a super-harmonic dynamics and the resulting single-to-double crest wave transition. The weakly nonlinear prediction is shown to be in fairly good agreement with previous experiments described in the literature. Lastly, we discuss how an analogous amplitude equation can be derived by solving asymptotically for the first super-harmonic of the forced Helmholtz-Duffing equation with small nonlinearities. △ Less

Submitted 21 March, 2022; originally announced March 2022.

Sub-harmonic parametric instability in nearly-brimful circular-cylinders: a weakly nonlinear analysis

Authors: Alessandro Bongarzone, Francesco Viola, Simone Camarri, François Gallaire

Abstract: In lab-scale Faraday experiments, meniscus waves respond harmonically to small-amplitude forcing without threshold, hence potentially cloaking the instability onset of parametric waves. Their suppression can be achieved by resorting to a contact line pinned at the container brim with static contact angle $θ_s=90^{\circ}$ (brimful condition). However, tunable meniscus waves are desired in some appl… ▽ More In lab-scale Faraday experiments, meniscus waves respond harmonically to small-amplitude forcing without threshold, hence potentially cloaking the instability onset of parametric waves. Their suppression can be achieved by resorting to a contact line pinned at the container brim with static contact angle $θ_s=90^{\circ}$ (brimful condition). However, tunable meniscus waves are desired in some applications as those of liquid-based biosensors, where they can be controlled adjusting the shape of the static meniscus by slightly under/over-filling the vessel ($θ_s\ne90^{\circ}$) while keeping the contact line fixed at the brim. Here, we refer to this wetting condition as nearly-brimful. Although classic inviscid theories based on Floquet analysis have been reformulated for the case of a pinned contact line (Kidambi 2013), accounting for (i) viscous dissipation and (ii) static contact angle effects, including meniscus waves, makes such analyses practically intractable and a comprehensive theoretical framework is still lacking. Aiming at filling this gap, in this work we formalize a weakly nonlinear analysis via multiple timescale method capable to predict the impact of (i) and (ii) on the instability onset of viscous sub-harmonic standing waves in both brimful and nearly-brimful circular-cylinders. Notwithstanding that the form of the resulting amplitude equation is in fact analogous to that obtained by symmetry arguments (Douady 1990), the normal form coefficients are here computed numerically from first principles, thus allowing us to rationalize and systematically quantify the modifications on the Faraday tongues and on the associated bifurcation diagrams induced by the interaction of meniscus and sub-harmonic parametric waves. △ Less

Submitted 28 February, 2022; originally announced February 2022.

Gravity-driven coatings on curved substrates: a differential geometry approach

Authors: Pier Giuseppe Ledda, Matteo Pezzulla, Etienne Jambon-Puillet, Pierre-Thomas Brun, François Gallaire

Abstract: Although the drainage and spreading processes of thin liquid films on substrates have received growing attention during the last decades, the study of three-dimensional cases is limited to a few studies on flat and axisymmetric substrates. In this work, we exploit differential geometry to study the drainage and spreading of thin films on generic curved substrates. We initially investigate the drai… ▽ More Although the drainage and spreading processes of thin liquid films on substrates have received growing attention during the last decades, the study of three-dimensional cases is limited to a few studies on flat and axisymmetric substrates. In this work, we exploit differential geometry to study the drainage and spreading of thin films on generic curved substrates. We initially investigate the drainage and spreading processes on spheroidal and paraboloidal substrates by employing an asymptotic expansion in the vicinity of the pole and a self-similar approach, finding that the thickness distribution is set by the substrate metric and tangential gravity force components. Spheroids with a large ratio between height and equatorial radius are characterized by a growing thickness moving away from the pole, and vice versa. The non-symmetric coating on a toroidal substrate shows that larger thicknesses and a faster spreading are attained on the inner region than on the outer region of the torus. An ellipsoid with three different axes is chosen as a testing ground for three-dimensional drainage and spreading. Modulations in the drainage solution are observed, with a different variation of the thickness along the two axes. By imposing the conservation of mass, an analytical solution for the average spreading front is obtained. The analytical and numerical results are in good agreement. The resulting drainage solutions show also a good agreement with experimental measurements obtained from the coating of a curing polymer on diverse substrates. △ Less

Submitted 9 July, 2022; v1 submitted 23 February, 2022; originally announced February 2022.

Comments: 49 pages, 28 figures (including Supporting Information)

Fiber buckling in confined viscous flows: an absolute instability described by the Ginzburg-Landau equation

Authors: Jean Cappello, Olivia du Roure, François Gallaire, Camille Duprat, Anke Lindner

Abstract: We explore the dynamics of a flexible fiber transported by a viscous flow in a Hele-Shaw cell of height comparable to the fiber height. We show that long fibers aligned with the flow experience a buckling instability. Competition between viscous and elastic forces leads to the deformation of the fiber into a wavy shape convolved by a Bell-shaped envelope. We characterize the wavelength and phase v… ▽ More We explore the dynamics of a flexible fiber transported by a viscous flow in a Hele-Shaw cell of height comparable to the fiber height. We show that long fibers aligned with the flow experience a buckling instability. Competition between viscous and elastic forces leads to the deformation of the fiber into a wavy shape convolved by a Bell-shaped envelope. We characterize the wavelength and phase velocity of the deformation as well as the growth and spreading of the envelope. Our study of the spatio-temporal evolution of the deformation reveals a linear and absolute instability arising from a local mechanism well described by the Ginzburg-Landau equation. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Weak nonlinearity for strong nonnormality

Authors: Yves-Marie Ducimetière, Edouard Boujo, François Gallaire

Abstract: We propose a theoretical approach to derive amplitude equations governing the weakly nonlinear evolution of nonnormal dynamical systems when they experience transient growth or respond to harmonic forcing. This approach reconciles the nonmodal nature of these growth mechanisms and the need for a center manifold to project the leading-order dynamics. Under the hypothesis of strong nonnormality, we… ▽ More We propose a theoretical approach to derive amplitude equations governing the weakly nonlinear evolution of nonnormal dynamical systems when they experience transient growth or respond to harmonic forcing. This approach reconciles the nonmodal nature of these growth mechanisms and the need for a center manifold to project the leading-order dynamics. Under the hypothesis of strong nonnormality, we take advantage of the fact that small operator perturbations suffice to make the inverse resolvent and the inverse propagator singular, which we encompass in a multiple-scale asymptotic expansion. The methodology is outlined for a generic nonlinear dynamical system, and four application cases highlight common nonnormal mechanisms in hydrodynamics: the streamwise convective nonnormal amplification in the flow past a backward-facing step, and the Orr and lift-up mechanisms in the plane Poiseuille flow. △ Less

Submitted 19 April, 2022; v1 submitted 15 October, 2021; originally announced October 2021.

Comments: 37 pages, 14 figures

Prediction of Self-Assembled Dewetted Nanostructures for Photonics Applications via a Continuum Mechanics Framework

Authors: Louis Martin-Monier, Pier Giuseppe Ledda, Pierre-Luc Piveteau, Francois Gallaire, Fabien Sorin

Abstract: When a liquid film lies on a non-wettable substrate, the configuration is unstable and the film then retracts from a solid substrate to form droplets. This phenomenon, known as dewetting, commonly leads to undesirable morphological changes. Nevertheless, recent works have demonstrated the possibility to harness dewetting by employing templated substrates with a degree of precision on par with adva… ▽ More When a liquid film lies on a non-wettable substrate, the configuration is unstable and the film then retracts from a solid substrate to form droplets. This phenomenon, known as dewetting, commonly leads to undesirable morphological changes. Nevertheless, recent works have demonstrated the possibility to harness dewetting by employing templated substrates with a degree of precision on par with advanced lithographic processes for high-performance nanophotonic applications. Since resonant behavior is highly sensitive to geometrical changes, predicting quantitatively dewetting dynamics is of high interest. In this work, we develop a continuum model that predicts the evolution of a thin film on a patterned substrate, from the initial reflow to the nucleation and growth of holes. We provide an operative framework based on macroscopic measurements to model the intermolecular interactions at the origin of the dewetting process, involving length scales that span from sub-nanometer to micron range. A comparison of experimental and simulated results shows that the model can accurately predict the final distributions, thereby offering novel predictive tools to tailor the optical response of dewetted nanostructures. △ Less

Submitted 4 May, 2021; originally announced May 2021.

Comments: 19 pages, 7 figures

Journal ref: Phys. Rev. Applied 16, 034025 (2021)

The effects of spanwise confinement on stratified shear instabilities

Authors: Yves-Marie Ducimetière, François Gallaire, Adrien Lefauve, Colm-cille P. Caulfield

Abstract: We consider the influence of transverse confinement on the instability properties of velocity and density distributions reminiscent of those pertaining to exchange flows in stratified inclined ducts, such as the recent experiment of Lefauve et al. (J. Fluid Mech. 848, 508-544, 2018). Using a normal mode streamwise and temporal expansion for flows in ducts with various aspect ratios $B$ and non-tri… ▽ More We consider the influence of transverse confinement on the instability properties of velocity and density distributions reminiscent of those pertaining to exchange flows in stratified inclined ducts, such as the recent experiment of Lefauve et al. (J. Fluid Mech. 848, 508-544, 2018). Using a normal mode streamwise and temporal expansion for flows in ducts with various aspect ratios $B$ and non-trivial transverse velocity profiles, we calculate two-dimensional (2D) dispersion relations with associated eigenfunctions varying in the 'crosswise' direction, in which the density varies, and the spanwise direction, both normal to the duct walls and to the flow direction. We also compare these 2D dispersion relations to the so-called one-dimensional (1D) dispersion relation obtained for spanwise invariant perturbations, for different aspect ratios $B$ and bulk Richardson numbers $Ri_b$. In this limited parameter space, the presence of lateral walls has a stabilizing effect. Furthermore, accounting for spanwise-varying perturbations results in a plethora of unstable modes, the number of which increases as the aspect ratio is increased. These modes present an odd-even regularity in their spatial structures, which is rationalized by comparison to the so-called one-dimensional oblique (1D-O) dispersion relation obtained for oblique waves. Finally, we show that in most cases, the most unstable 2D mode is the one that oscillates the least in the spanwise direction, as a consequence of viscous damping. However, in a limited region of the parameter space and in the absence of stratification, we show that a secondary mode with a more complex `twisted' structure dominated by crosswise vorticity becomes more unstable than the least oscillating Kelvin-Helmholtz mode associated with spanwise vorticity. △ Less

Submitted 12 October, 2021; v1 submitted 29 April, 2021; originally announced April 2021.

Comments: 35 pages, 12 figures

Journal ref: Phys. Rev. Fluids 6, 103901 (2021)

Leidenfrost flows: instabilities and symmetry breakings

Authors: Eunok Yim, Ambre Bouillant, David Quéré, François Gallaire

Abstract: Leidenfrost drops were recently found to host strong dynamics. In the present study, we investigate both experimentally and theoretically the {flows structures and stability} inside a Leidenfrost water drop as it evaporates, starting with a large puddle. As revealed by infrared mapping, the drop base is warmer than its apex by typically 10$^{\circ}$C, which is likely to trigger bulk thermobuoyant… ▽ More Leidenfrost drops were recently found to host strong dynamics. In the present study, we investigate both experimentally and theoretically the {flows structures and stability} inside a Leidenfrost water drop as it evaporates, starting with a large puddle. As revealed by infrared mapping, the drop base is warmer than its apex by typically 10$^{\circ}$C, which is likely to trigger bulk thermobuoyant flows and Marangoni surface flows. Tracer particles unveil complex and strong flows that undergo successive symmetry breakings as the drop evaporates. We investigate the linear stability of the baseflows in a non-deformable, quasi-static, levitating drop induced by thermobuoyancy and effective thermocapillary surface stress, using only one adjustable parameter. The stability analysis of nominally axisymmetric thermoconvective flows, parametrized by the drop radius $R$, yields the most unstable, {thus, dominant, azimuthal modes (of wavenumber $m$). Our theory predicts well the radii $R$ for the mode transitions and cascade with decreasing wavenumbers from $m=3$, $m=2$, down to $m=1$ (the eventual rolling mode that entails propulsion) as the drop shrinks in size}. The effect of the escaping vapor is not taken into account here, which may further destabilize the inner flow and couple to the liquid/vapor interface to give rise to motion Bouillant et al. (2018) [8] and Brandao et al. (2020) [9]. △ Less

Submitted 30 November, 2022; v1 submitted 18 December, 2020; originally announced December 2020.

Comments: 13 pages, 8 figures

Too fast to grow: Dynamics of pendant drops sliding on a thin film

Authors: Etienne Jambon-Puillet, Pier Giuseppe Ledda, François Gallaire, P. -T. Brun

Abstract: Pendant drops suspended on the underside of a wet substrate are known to accumulate fluid from the surrounding thin liquid film, a process that often results in dripping. The growth of such drops is hastened by their ability to translate over an otherwise uniform horizontal film. Here we show that this scenario is surprisingly reversed when the substrate is slightly tilted ($\approx 2$ deg) ; drop… ▽ More Pendant drops suspended on the underside of a wet substrate are known to accumulate fluid from the surrounding thin liquid film, a process that often results in dripping. The growth of such drops is hastened by their ability to translate over an otherwise uniform horizontal film. Here we show that this scenario is surprisingly reversed when the substrate is slightly tilted ($\approx 2$ deg) ; drops become too fast to grow and shrink over the course of their motion. Combining experiments and numerical simulations, we rationalize the transition between the conventional growth regime and the previously unknown decay regime we report. Using an analytical treatment of the Landau-Levich meniscus that connects the drop to the film, we quantitatively predict the drop dynamics in the two flow regimes and the value of the critical inclination angle where the transition between them occurs. △ Less

Submitted 30 July, 2021; v1 submitted 10 December, 2020; originally announced December 2020.

Journal ref: Phys. Rev. Lett. 127, 044503 (2021)

Nonlinear evolution of the centrifugal instability using a semi-linear model

Authors: Eunok Yim, Paul Billant, Francois Gallaire

Abstract: We study the nonlinear evolution of the centrifugal instability developing on a columnar anticyclone with a Gaussian angular velocity using a semi-linear approach. The model consists in two coupled equations: one for the linear evolution of the most unstable perturbation on the axially averaged mean flow and another for the evolution of the mean flow under the effect of the axially averaged Reynol… ▽ More We study the nonlinear evolution of the centrifugal instability developing on a columnar anticyclone with a Gaussian angular velocity using a semi-linear approach. The model consists in two coupled equations: one for the linear evolution of the most unstable perturbation on the axially averaged mean flow and another for the evolution of the mean flow under the effect of the axially averaged Reynolds stresses due to the perturbation. Such model is similar to the self-consistent model of \cite{Vlado14} except that the time averaging is replaced by a spatial averaging. The non-linear evolutions of the mean flow and the perturbations predicted by this semi-linear model are in very good agreement with DNS for the Rossby number $Ro=-4$ and both values of the Reynolds numbers investigated: $Re=800$ and $2000$ (based on the initial maximum angular velocity and radius of the vortex). An improved model taking into account the second harmonic perturbations is also considered. The results show that the angular momentum of the mean flow is homogenized towards a centrifugally stable profile via the action of the Reynolds stresses of the fluctuations. The final velocity profile predicted by \cite{Kloosterziel07} in the inviscid limit is extended to finite high Reynolds numbers. It is in good agreement with the numerical simulations. △ Less

Submitted 20 December, 2019; originally announced December 2019.

Frequency selection in a gravitationally stretched capillary jet in the jetting regime

Authors: Isha Shukla, François Gallaire

Abstract: A capillary jet falling under the effect of gravity continuously stretches while thinning downstream. We report here the effect of external periodic forcing on such a spatially varying jet in the jetting regime. Surprisingly, the optimal forcing frequency producing the most unstable jet is found to be highly dependent on the forcing amplitude. Taking benefit of the one-dimensional Eggers & Dupont… ▽ More A capillary jet falling under the effect of gravity continuously stretches while thinning downstream. We report here the effect of external periodic forcing on such a spatially varying jet in the jetting regime. Surprisingly, the optimal forcing frequency producing the most unstable jet is found to be highly dependent on the forcing amplitude. Taking benefit of the one-dimensional Eggers & Dupont (J. Fluid Mech., vol. 262, 1994, 205-221) equations, we investigate the case through nonlinear simulations and linear stability analysis. In the local framework the WKBJ formalism, established for weakly non-parallel flows, fails to capture the nonlinear simulation results quantitatively. However in the global framework, the resolvent analysis supplemented by a simple approximation of the required response norm inducing breakup, is shown to correctly predict the optimal forcing frequency at a given forcing amplitude and the resulting jet breakup length. The results of the resolvent analysis are found to be in good agreement with those of the nonlinear simulations. △ Less

Submitted 10 August, 2019; originally announced August 2019.

Optimal spanwise-periodic control for recirculation length in a backward-facing step flow

Authors: E. Yim, I. Shukla, F. Gallaire, E. Boujo

Abstract: Three-dimensional control is considered in the flow past a backward-facing step (BFS). The BFS flow at Reynolds number $Re=500$ (defined with the step height and the maximum inlet velocity) is two-dimensional and linearly stable but increasingly receptive to disturbances, with a potential for amplification as the recirculation length increases. We compute optimal spanwise-periodic control (steady… ▽ More Three-dimensional control is considered in the flow past a backward-facing step (BFS). The BFS flow at Reynolds number $Re=500$ (defined with the step height and the maximum inlet velocity) is two-dimensional and linearly stable but increasingly receptive to disturbances, with a potential for amplification as the recirculation length increases. We compute optimal spanwise-periodic control (steady wall blowing/suction or wall deformation) for decreasing the recirculation length, based on a second-order sensitivity analysis. Results show that wall-normal velocity control is always more efficient than wall-tangential control. The most efficient spanwise wavelength for the optimal control depends on the location: $β=0.6$ on the upper wall and $β=1$ on the upstream part of the lower wall. The linear amplification of the optimal control resembles the maximum linear gain, which confirms the link between recirculation length and amplification potential in this flow. Sensitivity predictions for blowing/suction amplitudes up to $O(10^{-3})$ and wall deformation amplitudes up to $O(10^{-2})$ are in good agreement with three-dimensional direct numerical simulations. For larger wall deformation amplitudes, the flow becomes unsteady. This study illustrates how the concept of second-order sensitivity and the associated optimization method allow for a systematic exploration of the best candidates for spanwise-periodic control. △ Less

Submitted 11 August, 2020; v1 submitted 3 July, 2019; originally announced July 2019.

Journal ref: Phys. Rev. Fluids 5, 083901 (2020)

Film thickness distribution in gravity-driven pancake-shaped droplets rising in a Hele-Shaw cell

Authors: Isha Shukla, Nicolas Kofman, Gioele Balestra, Lailai Zhu, François Gallaire

Abstract: We study here experimentally, numerically and using a lubrication approach; the shape, velocity and lubrication film thickness distribution of a droplet rising in a vertical Hele-Shaw cell. The droplet is surrounded by a stationary immiscible fluid and moves purely due to buoyancy. A low density difference between the two mediums helps to operate in a regime with capillary number $Ca$ lying betwee… ▽ More We study here experimentally, numerically and using a lubrication approach; the shape, velocity and lubrication film thickness distribution of a droplet rising in a vertical Hele-Shaw cell. The droplet is surrounded by a stationary immiscible fluid and moves purely due to buoyancy. A low density difference between the two mediums helps to operate in a regime with capillary number $Ca$ lying between $0.03-0.35$, where $Ca=μ_o U_d /γ$ is built with the surrounding oil viscosity $μ_o$, the droplet velocity $U_d$ and surface tension $γ$. The experimental data shows that in this regime the droplet velocity is not influenced by the thickness of the thin lubricating film and the dynamic meniscus. For iso-viscous cases, experimental and three-dimensional numerical results of the film thickness distribution agree well with each other. The mean film thickness is well captured by the Aussillous & Quéré (2000) model with fitting parameters. The droplet also exhibits the ''catamaran'' shape that has been identified experimentally for a pressure-driven counterpart (Huerre $\textit{et al}$. 2015). This pattern has been rationalized using a two-dimensional lubrication equation. In particular, we show that this peculiar film thickness distribution is intrinsically related to the anisotropy of the fluxes induced by the droplet's motion. △ Less

Submitted 17 June, 2019; originally announced June 2019.

Bifurcation dynamics of a particle-encapsulating droplet in shear flow

Authors: Lailai Zhu, François Gallaire

Abstract: To understand the behavior of composite fluid particles such as nucleated cells and double-emulsions in flow, we study a finite-size particle encapsulated in a deforming droplet under shear flow as a model system. In addition to its concentric particle-droplet configuration, we numerically explore other eccentric and time-periodic equilibrium solutions, which emerge spontaneously via supercritical… ▽ More To understand the behavior of composite fluid particles such as nucleated cells and double-emulsions in flow, we study a finite-size particle encapsulated in a deforming droplet under shear flow as a model system. In addition to its concentric particle-droplet configuration, we numerically explore other eccentric and time-periodic equilibrium solutions, which emerge spontaneously via supercritical pitchfork and Hopf bifurcations. We present the loci of these solutions around the codimenstion-two point. We adopt a dynamical system approach to model and characterize the coupled behavior of the two bifurcations. By exploring the flow fields and hydrodynamic forces in detail, we identify the role of hydrodynamic particle-droplet interaction which gives rise to these bifurcations. △ Less

Submitted 4 June, 2019; originally announced June 2019.

Journal ref: Phys. Rev. Lett., 119, 064502 (2017)

Unraveling radial dependency effects in fiber thermal drawing

Authors: Alexis G. Page, Mathias Bechert, François Gallaire, Fabien Sorin

Abstract: Fiber-based devices with advanced functionalities are emerging as promising solutions for various applications in flexible electronics and bioengineering. Multimaterial thermal drawing, in particular, has attracted strong interest for its ability to generate fibers with complex architectures. Thus far, however, the understanding of its fluid dynamics has only been applied to single material prefor… ▽ More Fiber-based devices with advanced functionalities are emerging as promising solutions for various applications in flexible electronics and bioengineering. Multimaterial thermal drawing, in particular, has attracted strong interest for its ability to generate fibers with complex architectures. Thus far, however, the understanding of its fluid dynamics has only been applied to single material preforms for which higher order effects, such as the radial dependency of the axial velocity, could be neglected. With complex multimaterial preforms, such effects must be taken into account, as they can affect the architecture and the functional properties of the resulting fiber device. Here, we propose a versatile model of the thermal drawing of fibers, which takes into account a radially varying axial velocity. Unlike the commonly used cross section averaged approach, our model is capable of predicting radial variations of functional properties caused by the deformation during drawing. This is demonstrated for two effects observed, namely, by unraveling the deformation of initially straight, transversal lines in the preform and the dependence on the draw ratio and radial position of the in-fiber electrical conductivity of polymer nanocomposites, an important class of materials for emerging fiber devices. This work sets a thus far missing theoretical and practical understanding of multimaterial fiber processing to better engineer advanced fibers and textiles for sensing, health care, robotics, or bioengineering applications. △ Less

Submitted 31 July, 2019; v1 submitted 4 February, 2019; originally announced March 2019.

Journal ref: Appl. Phys. Lett. 115, 044102 (2019)

Controlling transport dynamics of confined asymmetric fibers

Authors: M. Bechert, J. Cappello, M. Daïeff, F. Gallaire, A. Lindner, C. Duprat

Abstract: Transport properties of particles in confining geometries show very specific characteristics as lateral drift, oscillatory movement between lateral walls or the deformation of flexible fibers. These dynamics result from viscous friction with transversal and lateral channel walls inducing strong flow perturbations around the particles that act like moving obstacles. In this paper, we modify the fib… ▽ More Transport properties of particles in confining geometries show very specific characteristics as lateral drift, oscillatory movement between lateral walls or the deformation of flexible fibers. These dynamics result from viscous friction with transversal and lateral channel walls inducing strong flow perturbations around the particles that act like moving obstacles. In this paper, we modify the fiber shape by adding an additional, small fiber arm, which leads to T and L shaped fibers with only one or, respectively, zero symmetry axes and investigate the transport properties. For this purpose, we combine precise microfluidic experiments and numerical simulations based on modified Brinkman equations. Even for small shape perturbations, the transport dynamics change fundamentally and formerly stable configurations become unstable, leading to non-monotonous fiber rotation and lateral drift. Our results show that the fundamental transport dynamics change with respect to the level of fiber symmetry, which thus enables a precise control of particle trajectories and which can further be used for targeted delivery, particle sorting or capture inside microchannels. △ Less

Submitted 2 April, 2020; v1 submitted 2 March, 2019; originally announced March 2019.

Journal ref: 2019, EPL (Europhysics Letters), Volume 126, Number 4

Particle size selection in capillary instability of locally heated co-axial fiber

Authors: Saviz Mowlavi, Isha Shukla, PT Brun, François Gallaire

Abstract: Harnessing fluidic instabilities to produce structures with robust and regular properties has recently emerged as a new fabrication paradigm. This is exemplified in the work of Gumennik et al. [Nat. Comm. 4:2216, DOI: 10.1038/ncomms3216, (2013)], in which the authors fabricate silicon spheres by feeding a silicon-in-silica co-axial fiber into a flame. Following the localized melting of the silicon… ▽ More Harnessing fluidic instabilities to produce structures with robust and regular properties has recently emerged as a new fabrication paradigm. This is exemplified in the work of Gumennik et al. [Nat. Comm. 4:2216, DOI: 10.1038/ncomms3216, (2013)], in which the authors fabricate silicon spheres by feeding a silicon-in-silica co-axial fiber into a flame. Following the localized melting of the silicon, a capillary instability of the silicon-silica interface induces the formation of uniform silicon spheres. Here, we try to unravel the physical mechanisms at play in selecting the size of these particles, which was notably observed by Gumennik et al. to vary monotonically with the speed at which the fiber is fed into the flame. Using a simplified model derived from standard long-wavelength approximations, we show that linear stability analysis strikingly fails at predicting the selected particle size. Nonetheless, nonlinear simulations of the simplified model do recover the particle size observed in experiments, without any adjustable parameters. This shows that the formation of the silicon spheres in this system is an intrinsically nonlinear process that has little in common with the loss of stability of the underlying base flow solution. △ Less

Submitted 1 March, 2019; originally announced March 2019.

Comments: 23 pages, 15 figures

Viscous growth and rebound of a bubble near a rigid surface

Authors: Sebastien Michelin, Giacomo Gallino, Francois Gallaire, Eric Lauga

Abstract: Motivated by the dynamics of microbubbles near catalytic surfaces in bubble-powered microrockets, we consider theoretically the growth of a free spherical bubble near a flat no-slip surface in a Stokes flow. The flow at the bubble surface is characterised by a constant slip length allowing to tune the hydrodynamic mobility of its surface and tackle in one formulation both clean and contaminated bu… ▽ More Motivated by the dynamics of microbubbles near catalytic surfaces in bubble-powered microrockets, we consider theoretically the growth of a free spherical bubble near a flat no-slip surface in a Stokes flow. The flow at the bubble surface is characterised by a constant slip length allowing to tune the hydrodynamic mobility of its surface and tackle in one formulation both clean and contaminated bubbles as well as rigid shells. Starting with a bubble of infinitesimal size, the fluid flow and hydrodynamic forces on the growing bubble are obtained analytically. We demonstrate that, depending on the value of the bubble slip length relative to the initial distance to the wall, the bubble will either monotonically drain the fluid separating it from the wall, which will exponentially thin, or it will bounce off the surface once before eventually draining the thin film. Clean bubbles are shown to be a singular limit who always monotonically get repelled from the surface. The bouncing events for bubbles with finite slip lengths are further analysed in detail in the lubrication limit. In particular, we identify the origin of the reversal of the hydrodynamic force direction as due to the change in the flow pattern in the film between the bubble and the surface and to the associated lubrication pressure. Finally, the final drainage dynamics of the film is observed to follow a universal algebraic scaling for all finite slip lengths. △ Less

Submitted 20 November, 2018; v1 submitted 23 October, 2018; originally announced October 2018.

Comments: 22 pages; 7 figures; to appear in J. Fluid Mech

Journal ref: J. Fluid Mech., 2019, 860, 172-199

Second-order sensitivity in the cylinder wake: optimal spanwise-periodic wall actuation and wall deformation

Authors: E. Boujo, A. Fani, F. Gallaire

Abstract: Two-dimensional (2D) flows are efficiently controlled with spanwise waviness, i.e. spanwise-periodic (SP) wall blowing/suction/deformation. We tackle the global linear stability of 2D flows subject to small-amplitude 3D SP control. Building on previous work for parallel flows (Boujo et al. 2015), an adjoint method is proposed for computing the 2nd-order eigenvalue sensitivity. Such control has a z… ▽ More Two-dimensional (2D) flows are efficiently controlled with spanwise waviness, i.e. spanwise-periodic (SP) wall blowing/suction/deformation. We tackle the global linear stability of 2D flows subject to small-amplitude 3D SP control. Building on previous work for parallel flows (Boujo et al. 2015), an adjoint method is proposed for computing the 2nd-order eigenvalue sensitivity. Such control has a zero 1st-order linear effect, so the 2nd-order quadratic effect prevails. The sensitivity operator allows one (i) to predict the effect of any control without computing the controlled flow, (ii) to compute the optimal control for growth rate/frequency modification. The method takes advantage of spanwise periodicity to reduce computational complexity (from a 3D problem to a 2D one). The approach is applied to the leading eigenvalue of the laminar flow around a circular cylinder. Two SP controls are explored: wall actuation (blowing/suction) and wall deformation. Decomposing the eigenvalue variation, we find that the 3D contribution (from the SP 1st-order flow modification) is generally larger than the 2D one from the mean flow correction (spanwise-invariant 2nd-order flow modification). The optimal control for stabilization is top-down symmetric, leading to varicose streaks in the wake. Analyzing the competition between amplification and stabilization shows that optimal varicose streaks are not significantly more amplified than sinuous ones but have a stronger stabilizing effect. The optimal wall deformation induces a flow modification similar to that of the optimal wall actuation. In general, spanwise and tangential actuation have a small contribution to the optimal control, so normal-only actuation is a good trade-off between simplicity and effectiveness. Our method opens the way to the systematic design of optimal SP control for a variety of control objectives other than linear stability. △ Less

Submitted 3 October, 2019; v1 submitted 8 August, 2018; originally announced August 2018.

Comments: 27 pages, 23 figures, Physical Review Fluids

Journal ref: Physical Review Fluids 4(5), 053901 (2019)

Transport of flexible fibers in confined micro-channels

Authors: Jean Cappello, Mathias Bechert, Camille Duprat, Olivia du Roure, François Gallaire, Anke Lindner

Abstract: When transported in confined geometries rigid fibers show interesting transport dynamics induced by friction with the top and bottom walls. Fiber flexibility causes an additional coupling between fiber deformation and transport and is expected to lead to more complex dynamics. A first crucial step for their understanding is the characterization of the deformed fiber shape. Here we characterize thi… ▽ More When transported in confined geometries rigid fibers show interesting transport dynamics induced by friction with the top and bottom walls. Fiber flexibility causes an additional coupling between fiber deformation and transport and is expected to lead to more complex dynamics. A first crucial step for their understanding is the characterization of the deformed fiber shape. Here we characterize this shape for a fiber transported in a confined plug flow perpendicular to the flow direction using a combination of microfluidic experiments and numerical simulations. In the experiments, size, initial orientation, and mechanical properties of the fibers are controlled using micro-fabrication techniques and in-situ characterization methods.The numerical simulations use modified Brinkman equations as well as full 3D simulations. We show that the bending of a perpendicular fiber results from the force distribution acting on the elongated object and is proportional to the elasto-viscous number, comparing viscous to elastic forces. We quantitatively characterize the influence of the confinement on the fiber deformation. The precise understanding of the deformation of a flexible fiber in a confined geometry can also be used in future to understand the deformation and transport of more complex deformable particles in confined flows, as for example vesicles or red blood cells. △ Less

Submitted 13 July, 2018; originally announced July 2018.

Journal ref: Phys. Rev. Fluids 4, 034202 (2019)

Edge states control droplet break-up in sub-critical extensional flows

Authors: Giacomo Gallino, Tobias M. Schneider, François Gallaire

Abstract: A fluid droplet suspended in an extensional flow of moderate intensity may break into pieces, depending on the amplitude of the initial droplet deformation. In subcritical uniaxial extensional flow the non-breaking base state is linearly stable, implying that only a finite amplitude perturbation can trigger break-up. Consequently, the stable base solution is surrounded by its finite basin of attra… ▽ More A fluid droplet suspended in an extensional flow of moderate intensity may break into pieces, depending on the amplitude of the initial droplet deformation. In subcritical uniaxial extensional flow the non-breaking base state is linearly stable, implying that only a finite amplitude perturbation can trigger break-up. Consequently, the stable base solution is surrounded by its finite basin of attraction. The basin boundary, which separates initial droplet shapes returning to the non-breaking base state from those becoming unstable and breaking up, is characterized using edge tracking techniques. We numerically construct the edge state, a dynamically unstable equilibrium whose stable manifold forms the basin boundary. The edge state equilibrium controls if the droplet breaks and selects a unique path towards break-up. This path physically corresponds to the well-known end-pinching mechanism. Our results thereby rationalize the dynamics observed experimentally [Stone & Leal, J. Fluid Mech. 206, 223 (1989)] △ Less

Submitted 14 June, 2018; originally announced June 2018.

Journal ref: Phys. Rev. Fluids 3, 073603 (2018)

Suppression of von Kármán vortex streets past porous rectangular cylinders

Authors: Pier Giuseppe Ledda, Lorenzo Siconolfi, Francesco Viola, François Gallaire, Simone Camarri

Abstract: Although the stability properties of the wake past impervious bluff bodies have been widely examined in the literature, similar analyses regarding the flow around and through porous ones are still lacking. In this work, the effect of the porosity and permeability on the wake patterns of porous rectangular cylinders is numerically investigated at low to moderate Reynolds numbers in the framework of… ▽ More Although the stability properties of the wake past impervious bluff bodies have been widely examined in the literature, similar analyses regarding the flow around and through porous ones are still lacking. In this work, the effect of the porosity and permeability on the wake patterns of porous rectangular cylinders is numerically investigated at low to moderate Reynolds numbers in the framework of direct numerical simulation combined with local and global stability analyses. A modified Darcy-Brinkman formulation is employed here so as to describe the flow behavior inside the porous media, where also the convective terms are retained to correctly account for the inertial effects at high values of permeability. Different aspect ratios of the cylinder are considered, varying the thickness-to-height ratios, t/d, from 0.01 (flat plate) to 1.0 (square cylinder). The results show that the permeability of the bodies has a strong effect in modifying the characteristics of the wakes and of the associated flow instabilities, while the porosity weakly affects the resulting flow patterns. In particular, the fluid flows through the porous bodies and, thus, as the permeability is progressively increased, the recirculation regions, initially attached to the rear part of the bodies, at first detach from the body and, eventually, disappear even in the near wakes. Global stability analyses lead to the identification of critical values of the permeability above which any linear instability is prevented. Moreover, a different scaling of the non-dimensional permeability allows to identify a general threshold for all the configurations here studied that ensures the suppression of vortex shedding, at least in the considered parameter space. △ Less

Submitted 30 May, 2018; originally announced May 2018.

Comments: 31 pages and 17 figures

Journal ref: Phys. Rev. Fluids 3, 103901 (2018)

Oscillations of confined fibers transported in microchannels

Authors: Mathias Nagel, Pierre-Thomas Brun, Helene Berthet, Anke Lindner, François Gallaire, Camille Duprat

Abstract: We investigate the trajectories of rigid fibers as they are transported in a pressure-driven flow, at low Reynolds number, in shallow Hele Shaw cells. The transverse confinement and the resulting viscous friction on these elongated objects, as well as the lateral confinement (i.e. the presence of lateral walls), lead to complex fibers trajectories that we characterize with a combination of microfl… ▽ More We investigate the trajectories of rigid fibers as they are transported in a pressure-driven flow, at low Reynolds number, in shallow Hele Shaw cells. The transverse confinement and the resulting viscous friction on these elongated objects, as well as the lateral confinement (i.e. the presence of lateral walls), lead to complex fibers trajectories that we characterize with a combination of microfluidic experiments and simulations using modified Brinkman equations. We show that the transported fiber behaves as an oscillator for which we obtain and analyze a complete state diagram. △ Less

Submitted 1 May, 2018; originally announced May 2018.

Journal ref: Nagel, M., Brun, P., Berthet, H., Lindner, A., Gallaire, F., & Duprat, C. (2018). Oscillations of confined fibres transported in microchannels. Journal of Fluid Mechanics, 835, 444-470. doi:10.1017/jfm.2017.662

Physics of Bubble-Propelled Microrockets

Authors: Giacomo Gallino, Francois Gallaire, Eric Lauga, Sebastien Michelin

Abstract: A popular method to induce synthetic propulsion at the microscale is to use the forces created by surface-produced gas bubbles inside the asymmetric body of a catalytic swimmer (referred to in the literature as microrocket). Gas bubbles nucleate and grow within the catalytic swimmer and migrate toward one of its opening under the effect of asymmetric geometric confinement, thus generating a net hy… ▽ More A popular method to induce synthetic propulsion at the microscale is to use the forces created by surface-produced gas bubbles inside the asymmetric body of a catalytic swimmer (referred to in the literature as microrocket). Gas bubbles nucleate and grow within the catalytic swimmer and migrate toward one of its opening under the effect of asymmetric geometric confinement, thus generating a net hydrodynamic force which propels the device. In this paper we use numerical simulations to develop a joint chemical (diffusive) and hydrodynamic (Stokes) analysis of the bubble growth within a conical catalytic microrocket and of the associated bubble and microrocket motion. Our computational model allows us to solve for the bubble dynamics over one full bubble cycle ranging from its nucleation to its exiting the conical rocket and therefore to identify the propulsion characteristics as function of all design parameters, including geometry and chemical activity of the motor, surface tension phenomena, and all physicochemical constants. Our results suggest that hydrodynamics and chemistry partially decouple in the motion of the bubbles, with hydrodynamics determining the distance travelled by the microrocket over each cycle while chemistry setting the bubble ejection frequency. Our numerical model finally allows us to identify an optimal microrocket shape and size for which the swimming velocity (distance travelled per cycle duration) is maximized. △ Less

Submitted 28 March, 2018; originally announced March 2018.

Viscous Taylor droplets in axisymmetric and planar tubes: from Bretherton's theory to empirical models

Authors: Gioele Balestra, Lailai Zhu, François Gallaire

Abstract: The aim of this study is to derive accurate models for quantities characterizing the dynamics of droplets of non-vanishing viscosity in capillaries. In particular, we propose models for the uniform-film thickness separating the droplet from the tube walls, for the droplet front and rear curvatures and pressure jumps, and for the droplet velocity in a range of capillary numbers, $Ca$, from… ▽ More The aim of this study is to derive accurate models for quantities characterizing the dynamics of droplets of non-vanishing viscosity in capillaries. In particular, we propose models for the uniform-film thickness separating the droplet from the tube walls, for the droplet front and rear curvatures and pressure jumps, and for the droplet velocity in a range of capillary numbers, $Ca$, from $10^{-4}$ to $1$ and inner-to-outer viscosity ratios, $λ$, from $0$, i.e. a bubble, to high viscosity droplets. Theoretical asymptotic results obtained in the limit of small capillary number are combined with accurate numerical simulations at larger $Ca$. With these models at hand, we can compute the pressure drop induced by the droplet. The film thickness at low capillary numbers ($Ca<10^{-3}$) agrees well with Bretherton's scaling for bubbles as long as $λ<1$. For larger viscosity ratios, the film thickness increases monotonically, before saturating for $λ>10^3$ to a value $2^{2/3}$ times larger than the film thickness of a bubble. At larger capillary numbers, the film thickness follows the rational function proposed by Aussillous \& Quéré (2000) for bubbles, with a fitting coefficient which is viscosity-ratio dependent. This coefficient modifies the value to which the film thickness saturates at large capillary numbers. The velocity of the droplet is found to be strongly dependent on the capillary number and viscosity ratio. We also show that the normal viscous stresses at the front and rear caps of the droplets cannot be neglected when calculating the pressure drop for $Ca>10^{-3}$. △ Less

Submitted 2 July, 2018; v1 submitted 28 November, 2017; originally announced November 2017.

Absolute/convective secondary instabilities and the role of confinement in free shear layers

Authors: Cristobal Arratia, Saviz Mowlavi, Francois Gallaire

Abstract: Mixing layers can grow in time or space by vortex pairings that succeed each other in a nearly self-similar way. We use a point vortex model to study how confinement eventually limits this growth process, leading us to propose a wavelength selection mechanism for free shear layers with counterflow. Mixing layers can grow in time or space by vortex pairings that succeed each other in a nearly self-similar way. We use a point vortex model to study how confinement eventually limits this growth process, leading us to propose a wavelength selection mechanism for free shear layers with counterflow. △ Less

Submitted 7 May, 2018; v1 submitted 24 July, 2017; originally announced July 2017.

Comments: 20 pages, 9 figures

Journal ref: Phys. Rev. Fluids 3, 053901 (2018)

Swimming with a cage: Low-Reynolds-number locomotion inside a droplet

Authors: Shang Yik Reigh, Lailai Zhu, François Gallaire, Eric Lauga

Abstract: Inspired by recent experiments using synthetic microswimmers to manipulate droplets, we investigate the low-Reynolds-number locomotion of a model swimmer (a spherical squirmer) encapsulated inside a droplet of comparable size in another viscous fluid. Meditated solely by hydrodynamic interactions, the encaged swimmer is seen to be able to propel the droplet, and in some situations both remain in a… ▽ More Inspired by recent experiments using synthetic microswimmers to manipulate droplets, we investigate the low-Reynolds-number locomotion of a model swimmer (a spherical squirmer) encapsulated inside a droplet of comparable size in another viscous fluid. Meditated solely by hydrodynamic interactions, the encaged swimmer is seen to be able to propel the droplet, and in some situations both remain in a stable co-swimming state. The problem is tackled using both an exact analytical theory and a numerical implementation based on boundary element method, with a particular focus on the kinematics of the co-moving swimmer and droplet in a concentric configuration, and we obtain excellent quantitative agreement between the two. The droplet always moves slower than a swimmer which uses purely tangential surface actuation but when it uses a particular combination of tangential and normal actuations, the squirmer and droplet are able to attain a same velocity and stay concentric for all times. We next employ numerical simulations to examine the stability of their concentric co-movement, and highlight several stability scenarios depending on the particular gait adopted by the swimmer. Furthermore, we show that the droplet reverses the nature of the far-field flow induced by the swimmer: a droplet cage turns a pusher swimmer into a puller, and vice versa. Our work sheds light on the potential development of droplets as self-contained carriers of both chemical content and self-propelled devices for controllable and precise drug deliveries. △ Less

Submitted 11 March, 2017; originally announced March 2017.

Comments: 26pages, 10 figures, Soft Matter 2017

A pancake droplet translating in a Hele-Shaw cell: lubrication film and flow field

Authors: Lailai Zhu, François Gallaire

Abstract: We adopt a boundary integral method to study the dynamics of a translating droplet confined in a Hele-Shaw cell in the Stokes regime. The droplet is driven by the motion of the ambient fluid with the same viscosity. We characterize the three-dimensional (3D) nature of the droplet interface and of the flow field. The interface develops an arc-shaped ridge near the rear-half rim with a protrusion in… ▽ More We adopt a boundary integral method to study the dynamics of a translating droplet confined in a Hele-Shaw cell in the Stokes regime. The droplet is driven by the motion of the ambient fluid with the same viscosity. We characterize the three-dimensional (3D) nature of the droplet interface and of the flow field. The interface develops an arc-shaped ridge near the rear-half rim with a protrusion in the rear and a laterally symmetric pair of higher peaks; this pair of protrusions has been identified by recent experiments (Huerre et al., Phys. Rev. Lett., vol. 115 (6), 2015, 064501) and predicted asymptotically (Burgess and Foster, Phys. Fluids A, vol. 2 (7), 1990, pp. 1105-1117). The mean film thickness is well predicted by the extended Bretherton model (Klaseboer et al., Phys. Fluids, vol. 26 (3), 2014, 032107) with fitting parameters. The flow in the streamwise wall-normal middle plane is featured with recirculating zones, which are partitioned by stagnation points closely resembling those of a two-dimensional droplet in a channel. Recirculation is absent in the wall-parallel, unconfined planes, in sharp contrast to the interior flow inside a moving droplet in free space. The preferred orientation of the recirculation results from the anisotropic confinement of the Hele-Shaw cell. On these planes, we identify a dipolar disturbance flow field induced by the travelling droplet and its $1/r^2$ spatial decay is confirmed numerically. We pinpoint counter-rotating streamwise vortex structures near the lateral interface of the droplet, further highlighting the complex 3D flow pattern. △ Less

Submitted 27 June, 2016; v1 submitted 29 January, 2016; originally announced January 2016.

Journal ref: Journal of Fluid Mechanics, Volume 798, July 2016, pp 955-969

The stability of a rising droplet: an inertialess nonmodal growth mechanism

Authors: Giacomo Gallino, Lailai Zhu, Francois Gallaire

Abstract: Prior modal stability analysis (Kojima et al., Phys. Fluids, vol. 27, 1984) predicted that a rising or sedimenting droplet in a viscous fluid is stable in the presence of surface tension no matter how small, in contrast to experimental and numerical results. By performing a non-modal stability analysis, we demonstrate the potential for transient growth of the interfacial energy of a rising droplet… ▽ More Prior modal stability analysis (Kojima et al., Phys. Fluids, vol. 27, 1984) predicted that a rising or sedimenting droplet in a viscous fluid is stable in the presence of surface tension no matter how small, in contrast to experimental and numerical results. By performing a non-modal stability analysis, we demonstrate the potential for transient growth of the interfacial energy of a rising droplet in the limit of inertialess Stokes equations. The predicted critical capillary numbers for transient growth agree well with those for unstable shape evolution of droplets found in the direct numerical simulations of Koh & Leal (Phys. Fluids, vol. 1, 1989). Boundary integral simulations are used to delineate the critical amplitude of the most destabilizing perturbations. The critical amplitude is negatively correlated with the linear optimal energy growth, implying that the transient growth is responsible for reducing the necessary perturbation amplitude required to escape the basin of attraction of the spherical solution. △ Less

Submitted 30 November, 2015; v1 submitted 31 July, 2015; originally announced August 2015.

Comments: 11pages, 7 figures

Journal ref: Journal of Fluid Mechanics, Rapids., 786, 2016

Breaking anchored droplets in a microfluidic Hele-Shaw cell

Authors: Gabriel Amselem, P. T. Brun, François Gallaire, Charles N. Baroud

Abstract: We study microfluidic self digitization in Hele-Shaw cells using pancake droplets anchored to surface tension traps. We show that above a critical flow rate, large anchored droplets break up to form two daughter droplets, one of which remains in the anchor. Below the critical flow velocity for breakup the shape of the anchored drop is given by an elastica equation that depends on the capillary num… ▽ More We study microfluidic self digitization in Hele-Shaw cells using pancake droplets anchored to surface tension traps. We show that above a critical flow rate, large anchored droplets break up to form two daughter droplets, one of which remains in the anchor. Below the critical flow velocity for breakup the shape of the anchored drop is given by an elastica equation that depends on the capillary number of the outer fluid. As the velocity crosses the critical value, the equation stops admitting a solution that satisfies the boundary conditions; the drop breaks up in spite of the neck still having finite width. A similar breaking event also takes place between the holes of an array of anchors, which we use to produce a 2D array of stationary drops in situ. △ Less

Submitted 6 April, 2015; originally announced April 2015.

Comments: 5 pages, 4 figures, to appear in Phys. Rev. Applied

Rayleigh-Taylor instability under an inclined plane

Authors: P. -T Brun, Adam Damiano, Pierre Rieu, Gioele Balestra, François Gallaire

Abstract: We revisit the canonical Rayleigh-Taylor instability and investigate the case of a thin film of fluid upon the underside of an inclined plane. The presence of a natural flow along the plane competes with the conventional droplet forming instability. In particular, experiments reveal that no drops form for inclinations greater than a critical value. These features are rationalized in the context of… ▽ More We revisit the canonical Rayleigh-Taylor instability and investigate the case of a thin film of fluid upon the underside of an inclined plane. The presence of a natural flow along the plane competes with the conventional droplet forming instability. In particular, experiments reveal that no drops form for inclinations greater than a critical value. These features are rationalized in the context of the absolute/convective analysis conducted in this article. △ Less

Submitted 8 March, 2015; originally announced March 2015.

Boundary elements method for microfluidic two-phase flows in shallow channels

Authors: Mathias Nagel, François Gallaire

Abstract: In the following work we apply the boundary element method to two-phase flows in shallow microchannels, where one phase is dispersed and does not wet the channel walls. These kinds of flows are often encountered in microfluidic Lab-on-a-Chip devices and characterized by low Reynolds and low capillary numbers. Assuming that these channels are homogeneous in height and have a large aspect ratio, w… ▽ More In the following work we apply the boundary element method to two-phase flows in shallow microchannels, where one phase is dispersed and does not wet the channel walls. These kinds of flows are often encountered in microfluidic Lab-on-a-Chip devices and characterized by low Reynolds and low capillary numbers. Assuming that these channels are homogeneous in height and have a large aspect ratio, we use depth-averaged equations to describe these two-phase flows using the Brinkman equation, which constitutes a refinement of Darcy's law. These partial differential equations are discretized and solved numerically using the boundary element method, where a stabilization scheme is applied to the surface tension terms, allowing for a less restrictive time step at low capillary numbers. The convergence of the numerical algorithm is checked against a static analytical solution and on a dynamic test case. Finally the algorithm is applied to the non-linear development of the Saffman-Taylor instability and compared to experimental studies of droplet deformation in expanding flows. △ Less

Submitted 11 November, 2014; originally announced November 2014.

Comments: accepted for publication, Computers and Fluids 2014

Manipulating flow separation: sensitivity of stagnation points, separatrix angles and recirculation area to steady actuation

Authors: Edouard Boujo, François Gallaire

Abstract: A variational technique is used to derive analytical expressions for the sensitivity of several geometric indicators of flow separation to steady actuation. Considering the boundary layer flow above a wall-mounted bump, the six following representative quantities are considered: the locations of the separation point and reattachment point connected by the separatrix, the separation angles at these… ▽ More A variational technique is used to derive analytical expressions for the sensitivity of several geometric indicators of flow separation to steady actuation. Considering the boundary layer flow above a wall-mounted bump, the six following representative quantities are considered: the locations of the separation point and reattachment point connected by the separatrix, the separation angles at these stagnation points, the backflow area and the recirculation area. For each geometric quantity, linear sensitivity analysis allows us to identify regions which are the most sensitive to volume forcing and wall blowing/suction. Validations against full non-linear Navier-Stokes calculations show excellent agreement for small-amplitude control for all considered indicators. With very resemblant sensitivity maps, the reattachment point, the backflow and recirculation areas are seen to be easily manipulated. In contrast, the upstream separation point and the separatrix angles are seen to remain extremely robust with respect to external steady actuation. △ Less

Submitted 10 November, 2014; originally announced November 2014.

Comments: 20 pages, 12 figures

Journal ref: Proc. R. Soc. A 2014, vol. 470, no. 2170, 20140365

Open-loop control of noise amplification in a separated boundary layer flow

Authors: Edouard Boujo, Uwe Ehrenstein, François Gallaire

Abstract: Linear optimal gains are computed for the subcritical two-dimensional separated boundary-layer flow past a bump. Very large optimal gain values are found, making it possible for small-amplitude noise to be strongly amplified and to destabilize the flow. The optimal forcing is located close to the summit of the bump, while the optimal response is the largest in the shear layer. The largest amplific… ▽ More Linear optimal gains are computed for the subcritical two-dimensional separated boundary-layer flow past a bump. Very large optimal gain values are found, making it possible for small-amplitude noise to be strongly amplified and to destabilize the flow. The optimal forcing is located close to the summit of the bump, while the optimal response is the largest in the shear layer. The largest amplification occurs at frequencies corresponding to eigenvalues which first become unstable at higher Reynolds number. Nonlinear direct numerical simulations show that a low level of noise is indeed sufficient to trigger random flow unsteadiness, characterized here by large-scale vortex shedding. Next, a variational technique is used to compute efficiently the sensitivity of optimal gains to steady control (through source of momentum in the flow, or blowing/suction at the wall). A systematic analysis at several frequencies identifies the bump summit as the most sensitive region for control with wall actuation. Based on these results, a simple open-loop control strategy is designed, with steady wall suction at the bump summit. Linear calculations on controlled base flows confirm that optimal gains can be drastically reduced at all frequencies. Nonlinear direct numerical simulations also show that this control allows the flow to withstand a higher level of stochastic noise without becoming nonlinearly unstable, thereby postponing bypass transition. In the supercritical regime, sensitivity analysis of eigenvalues supports the choice of this control design. Full restabilization of the flow is obtained, as evidenced by direct numerical simulations and linear stability analysis. △ Less

Submitted 8 November, 2014; originally announced November 2014.

Comments: 33 pages, 20 figures

Journal ref: Physics of Fluids 25, 124106 (2013)

Sensitivity and open-loop control of stochastic response in a noise amplifier flow: the backward-facing step

Authors: Edouard Boujo, François Gallaire

Abstract: The two-dimensional backward-facing step flow is a canonical example of noise amplifier flow: global linear stability analysis predicts that it is stable, but perturbations can undergo large amplification in space and time as a result of non-normal effects. This amplification potential is best captured by optimal transient growth analysis, optimal harmonic forcing, or the response to sustained noi… ▽ More The two-dimensional backward-facing step flow is a canonical example of noise amplifier flow: global linear stability analysis predicts that it is stable, but perturbations can undergo large amplification in space and time as a result of non-normal effects. This amplification potential is best captured by optimal transient growth analysis, optimal harmonic forcing, or the response to sustained noise. In view of reducing disturbance amplification in these globally stable open flows, a variational technique is proposed to evaluate the sensitivity of stochastic amplification to steady control. Existing sensitivity methods are extended in two ways to achieve a realistic representation of incoming noise: (i) perturbations are time-stochastic rather than time-harmonic, (ii) perturbations are localised at the inlet rather than distributed in space. This allows for the identification of regions where small-amplitude control is the most effective, without actually computing any controlled flows. In particular, passive control by means of a small cylinder and active control by means of wall blowing/suction are analysed for Reynolds number $Re=500$ and step-to-outlet expansion ratio $Γ=0.5$. Sensitivity maps for noise amplification appear largely similar to sensitivity maps for optimal harmonic amplification at the most amplified frequency. This is observed at other values of $Re$ and $Γ$ too, and suggests that the design of steady control in this noise amplifier flow can be simplified by focusing on the most dangerous perturbation at the most dangerous frequency. △ Less

Submitted 7 November, 2014; originally announced November 2014.

Comments: 31 pages, 24 figures

Journal ref: Journal of Fluid Mechanics, 762, 2015

Controlled reattachment in separated flows: a variational approach to recirculation length reduction

Authors: E. Boujo, F. Gallaire

Abstract: A variational technique is used to derive analytical expressions for the sensitivity of recirculation length to steady forcing in separated flows. Linear sensitivity analysis is applied to the two-dimensional steady flow past a circular cylinder for Reynolds numbers $40 \leq Re \leq 120$, both in the subcritical and supercritical regimes. Regions which are the most sensitive to volume forcing and… ▽ More A variational technique is used to derive analytical expressions for the sensitivity of recirculation length to steady forcing in separated flows. Linear sensitivity analysis is applied to the two-dimensional steady flow past a circular cylinder for Reynolds numbers $40 \leq Re \leq 120$, both in the subcritical and supercritical regimes. Regions which are the most sensitive to volume forcing and wall blowing/suction are identified. Control configurations which reduce the recirculation length are designed based on the sensitivity information, in particular small cylinders used as control devices in the wake of the main cylinder, and fluid suction at the cylinder wall. Validation against full non-linear Navier-Stokes calculations shows excellent agreement for small-amplitude control. The linear stability properties of the controlled flow are systematically investigated. At moderate Reynolds numbers, we observe that regions where control reduces the recirculation length correspond to regions where it has a stabilising effect on the most unstable global mode associated to vortex shedding, while this property does not hold any more at larger Reynolds numbers. △ Less

Submitted 2 October, 2014; originally announced October 2014.

Comments: 17 pages, 11 figures

Journal ref: Journal of Fluid Mechanics, 742, 2014

Second-order sensitivity of parallel shear flows and optimal spanwise-periodic flow modifications

Authors: E. Boujo, A. Fani, F. Gallaire

Abstract: The question of optimal spanwise-periodic modification for the stabilisation of spanwise-invariant flows is addressed. A 2nd-order sensitivity analysis is conducted for the linear temporal stability of parallel flows U0 subject to small-amplitude spanwise-periodic modification e*U1, e<<1. Spanwise-periodic modifications have a quadratic effect on stability, i.e. the 1st-order eigenvalue variation… ▽ More The question of optimal spanwise-periodic modification for the stabilisation of spanwise-invariant flows is addressed. A 2nd-order sensitivity analysis is conducted for the linear temporal stability of parallel flows U0 subject to small-amplitude spanwise-periodic modification e*U1, e<<1. Spanwise-periodic modifications have a quadratic effect on stability, i.e. the 1st-order eigenvalue variation is zero. A 2nd-order sensitivity operator is computed from a 1D calculation, allowing one to predict how eigenvalues are affected by any U1, without actually solving for modified eigenvalues/eigenmodes. Comparisons with full 2D stability calculations in a plane channel flow and in a mixing layer show excellent agreement. Next, optimisation is performed on the 2nd-order sensitivity operator: for each eigenmode streamwise wavenumber and base flow modification spanwise wavenumber b, the most stabilising profiles U1 are computed, together with lower bounds for the variation in leading eigenvalue. These bounds increase like b^-2 as b goes to 0, yielding a large stabilising potential. However, 3D modes with wavenumbers |b0|=b and b/2 are destabilised, thus larger control wavenumbers should be preferred. The modification U1 optimised for the most unstable streamwise wavenumber has a stabilising effect on other streamwise wavenumbers too. Finally, the potential of transient growth to amplify perturbations and stabilise the flow is assessed. Combined optimal perturbations that achieve the best balance between transient linear amplification and flow stabilisation are determined. In the mixing layer with b<1.5, these combined optimal perturbations appear similar to transient growth-only optimal perturbations, and achieve a more efficient overall stabilisation than optimal 1D and 2D modifications computed for stabilisation only. This is consistent with the efficiency of streak-based control strategies. △ Less

Submitted 27 October, 2015; v1 submitted 2 October, 2014; originally announced October 2014.

Comments: 23 pages, 15 figures

Journal ref: Journal of Fluid Mechanics, 782, 2015