A coupled Immersed Boundary – Lattice Boltzmann method for incompressible flows through moving porous media A coupled Immersed Boundary -Lattice Boltzmann method for incompressible flows through moving porous media

International audienceIn this work, we propose a numerical framework to simulate fluid flows in interaction with moving porous media of complex geometry. It is based on the Lattice Boltzmann method including porous effects via a Brinkman-Forchheimer-Darcy force model coupled to the Immersed Boundary method to handle complex ge-ometries and moving structures. The coupling algorithm is described in detail and it is validated on well-established literature test cases for both stationary and moving porous configurations. The proposed method is easy to implement and efficient in terms of CPU cost and memory management compared to alternative methods which can be used to deal with moving immersed porous media, e.g. re-meshing at each time step or use of a moving/chimera mesh. An overall good agreement was obtained with reference results, opening the way to the numerical simulation of moving porous media for flow control applications

Transport efficiency of metachronal waves in 3d cilia arrays immersed in a two-phase flow

The present work reports the formation and the characterization of antipleptic and symplectic metachronal waves in 3D cilia arrays immersed in a two-fluid environment, with a viscosity ratio of 20. A coupled lattice-Boltzmann-Immersed-Boundary solver is used. The periciliary layer is confined between the epithelial surface and the mucus. Its thickness is chosen such that the tips of the cilia can penetrate the mucus. A purely hydrodynamical feedback of the fluid is taken into account and a coupling parameter $\alpha$ is introduced allowing the tuning of both the direction of the wave propagation, and the strength of the fluid feedback. A comparative study of both antipleptic and symplectic waves, mapping a cilia inter-spacing ranging from 1.67 up to 5 cilia length, is performed by imposing the metachrony. Antipleptic waves are found to systematically outperform sympletic waves. They are shown to be more efficient for transporting and mixing the fluids, while spending less energy than symplectic, random, or synchronized motions

Explicit and viscosity-independent immersed-boundary scheme for the lattice Boltzmann method

International audienceViscosity independence of lattice-Boltzmann methods is a crucial issue to ensure the physical relevancy of the predicted macroscopic flows over large ranges of physical parameters. The immersed-boundary (IB) method, a powerful tool that allows one to immerse arbitrary-shaped, moving, and deformable bodies in the flow, suffers from a boundary-slip error that increases as a function of the fluid viscosity, substantially limiting its range of application. In addition, low fluid viscosities may result in spurious oscillations of the macroscopic quantities in the vicinity of the immersed boundary. In this work, it is shown mathematically that the standard IB method is indeed not able to reproduce the scaling properties of the macroscopic solution, leading to a viscosity-related error on the computed IB force. The analysis allows us to propose a simple correction of the IB scheme that is local, straightforward and does not involve additional computational time. The derived method is implemented in a two-relaxation-time D2Q9 lattice-Boltzmann solver, applied to several physical configurations, namely, the Poiseuille flow, the flow around a cylinder towed in still fluid, and the flow around a cylinder oscillating in still fluid, and compared to a noncorrected immersed-boundary method. The proposed correction leads to a major improvement of the viscosity independence of the solver over a wide range of relaxation times (from 0.5001 to 50), including the correction of the boundary-slip error and the suppression of the spurious oscillations. This improvement may considerably extend the range of application of the IB lattice-Boltzmann method, in particular providing a robust tool for the numerical analysis of physical problems involving fluids of varying viscosity interacting with solid geometries

Contrôle d'écoulements : approche expérimentale et modélisation de dimension réduite

L'aptitude à déterminer une loi de contrôle optimale pour un objectif donné et construire un actionneur adapté, représente un intérêt industriel et écologique croissant, ainsi qu'un défi scientifique majeur auquel cette thèse apporte une contribution au moyen de deux approches méthodologiquement différentes. La première, basée sur l'expérience, vise à caractériser, par mesures PIV et analyse POD, un actionneur de soufflage capable de contrôler le décollement autour d'un profil d'aile. La seconde approche proposée pour optimiser le contrôle consiste à déterminer un modèle POD de dimension réduite de la dynamique de l'écoulement à contrôler, pour calculer le contrôle optimal à moindres coûts numériques. Des méthodes de calibration de modèles associés à différentes dynamiques sont développées pour améliorer la qualité des approximations. ABSTRACT : The ability to determine the optimal control for a given objective and construct an appropiate actuator, represents growing industrial and economical issues as well as a major scientific challenge. This thesis provides a contribution to that problem using two approaches, based on different methodologies. The first, experimental, aims at characterizing, using PIV measurements and POD analysis, a blowing actuator to control boundary layer separation on an airfoil. The idea of the second approach proposed here is to establish a POD Reduced-Order Model (POD ROM) of the flow dynamics to be controlled, in order to compute the optimal control with low CPU and memory costs. Calibration methods are then developped to improve the accuracy of this approximation based on different dynamics

Why antiplectic metachronal cilia waves are optimal to transport bronchial mucus

International audienceThe coordinated beating of epithelial cilia in human lungs is a fascinating problem from the hydrodynamics perspective. The phase lag between neighboring cilia is able to generate collective cilia motions, known as metachronal waves. Different kinds of waves can occur, antiplectic or symplectic, depending on the direction of the wave with respect to the flow direction. It is shown here, using a coupled lattice Boltzmann-immersed boundary solver, that the key mechanism responsible for their transport efficiency is a blowing-suction effect that displaces the interface between the periciliary liquid and the mucus phase. The contribution of this mechanism on the average flow generated by the cilia is compared to the contribution of the lubrication effect. The results reveal that the interface displacement is the main mechanism responsible for the better efficiency of antiplectic metachronal waves over symplectic ones to transport bronchial mucus. The conclusions drawn here can be extended to any two-layer fluid configuration having different viscosities, and put into motion by cilia-shaped or comb-plate structures, having a back-and-forth motion with phase lags

Experimental Investigations on Fluidic Control Over an Airfoil

International audienceThis study presents the development of two fluidic actuators − namely, microjets and tangential blowing actuator (TBA), designed for flow separation control. The developed actuators are compact enough to fit inside an ONERA D profiled wing with a chord of 0.35 m. Test bench experiments showed that the microjets (resp. TBA) were able to produce exit velocities up to 330 m/s (resp. 60 m/s). These actuators were placed in the model and were tested in wind tunnels for various blowing rates. The investigations included the use of force balance measurements, on-surface flow visualization with pigmented oil, off-surface flow visualizations with smoke, surface pressure distribution measurements, and Particle Image Velocimetry (PIV). Most of the tests were performed at free-stream velocities between 20 m/s (for PIV) and 40 m/s, corresponding to Reynolds numbers in the range 0.47 × 10^6−0.93 × 10^6 . The angle of attack varied from −2 to 20 degrees. Experiments were conducted using the naturally occurring laminar boundary layer as well as for a turbulent boundary layer. In such a case, rough strips were used in the vicinity of the leading edge. The present tests show the efficiency of these devices to delay separation and improve aerodynamic performances of the wing: for example, a maximum of 30% gain in CL has been reached using the microjets. Both actuators tend to increase the lift coefficient CL after stall and areas of separated flow have been eliminated by applying control, as suggested by flow visualizations and PIV velocity fields

An immersed boundary-lattice Boltzmann method for single- and multi-component fluid flows

International audienceThe paper presents a numerical method to simulate single-and multi-component fluid flows around moving/deformable solid boundaries, based on the coupling of Immersed Boundary (IB) and Lattice Boltzmann (LB) methods. The fluid domain is simulated with LB method using the single relaxation time BGK model, in which an interparticle potential model is applied for multi-component fluid flows. The IB-related force is directly calculated with the interpolated definition of the fluid macroscopic velocity on the Lagrangian points that define the immersed solid boundary. The present IB-LB method can better ensure the no-slip solid boundary condition, thanks to an improved spreading operator. The proposed method is validated through several 2D/3D single-and multi-component fluid test cases with a particular emphasis on wetting conditions on solid wall. Finally, a 3D two-fluid application case is given to show the feasibility of modeling the fluid transport via a cluster of beating cilia

From a vortex gas to a vortex crystal in instability-driven two-dimensional turbulence

We study structure formation in two-dimensional turbulence driven by an external force, interpolating between linear instability forcing and random stirring, subject to nonlinear damping. Using extensive direct numerical simulations, we uncover a rich parameter space featuring four distinct branches of stationary solutions: large-scale vortices, hybrid states with embedded shielded vortices (SVs) of either sign, and two states composed of many similar SVs. Of the latter, the first is a dense vortex gas where all SVs have the same sign and diffuse across the domain. The second is a hexagonal vortex crystal forming from this gas when the instability is sufficiently weak. These solutions coexist stably over a wide parameter range. The late-time evolution of the system from small-amplitude initial conditions is nearly self-similar, involving three phases: initial inverse cascade, random nucleation of SVs from turbulence and, once a critical number of vortices is reached, a phase of explosive nucleation of SVs, leading to a statistically stationary state. The vortex gas is continued in the forcing parameter, revealing a sharp transition towards the crystal state as the forcing strength decreases. This transition is analysed in terms of the diffusion of individual vortices and tools from statistical physics. The crystal can also decay via an inverse cascade resulting from the breakdown of shielding or insufficient nonlinear damping acting on SVs. Our study highlights the importance of the forcing details in two-dimensional turbulence and reveals the presence of nontrivial SV states in this system, specifically the emergence and melting of a vortex crystal

Spontaneous suppression of inverse energy cascade in instability-driven 2D turbulence

Instabilities of fluid flows often generate turbulence. Using extensive direct numerical simulations, we study two-dimensional turbulence driven by a wavenumber-localised instability superposed on stochastic forcing, in contrast to previous studies of state-independent forcing. As the contribution of the instability forcing, measured by a parameter $\gamma$, increases, the system undergoes two transitions. For $\gamma$ below a first threshold, a regular large-scale vortex condensate forms. Above this threshold, shielded vortices (SVs) emerge within the condensate. At a second, larger value of $\gamma$, the condensate breaks down, and a gas of weakly interacting vortices with broken symmetry spontaneously emerges, characterised by preponderance of vortices of one sign only and suppressed inverse energy cascade. The latter transition is shown to depend on the damping mechanism. The number density of SVs in the broken symmetry state slowly increases via a random nucleation process. Bistability is observed between the condensate and mixed SV-condensate states. Our findings provide new evidence for a strong dependence of two-dimensional turbulence phenomenology on the forcing

Sulla scia di Icaro

La progettazione di velivoli dal minimo attrito tiene impegnati i ricercatori del DICAT di Genova che lavorano al progetto europeo FLUBIO. La realizzazione del "perfetto automa volante" trae ispirazione dalla natura e in particolare dal volo dell'airone, un ottimo modello per velivoli con ali ricoperte da piume che battono al vent