It is well known that when a flexible or flexibly-mounted structure is placed perpendicular to the flow of a Newtonian fluid, it can oscillate due to the shedding of separated vortices at high Reynolds numbers. If the same flexible object is placed in non-Newtonian flows, however, the structure's response is still unknown. The main objective of this thesis is to introduce a new field of viscoelastic fluid-structure interactions by showing that the elastic instabilities that occur in the flow of viscoelastic fluids can drive the motion of a flexible structure placed in its path. Unlike Newtonian fluids, the flow of viscoelastic fluids can become unstable at infinitesimal Reynolds numbers due to the onset of a purely elastic flow instability. This instability occurs in the absence of nonlinear effects of fluid inertia and the Reynolds number of the flows studied here are in the order of 10-4. When such an elastic flow instability occurs in the vicinity of a flexible structure, the fluctuating fluid forces exerted on the structure grow large enough to cause a structural instability which in turn feeds back into the fluid resulting in a flow instability. Nonlinear periodic oscillations of the flexible structure are observed which have been found to be coupled to the time-dependent growth and decay of viscoelastic stresses in the wake of the structure. Presented in this thesis are the results of an investigation of the interaction occurring in the flow of a viscoelastic wormlike micelle solution past a flexible rectangular sheet. The structural geometries studied include: flexible sheet inclinations at 20°, 45° and 90° and flexible sheet widths of 5mm and 2.5mm. By varying the flow velocity, the response of the flexible sheet has been characterized in terms of amplitude and frequency of oscillations. Steady and dynamic shear rheology and filament stretching extensional rheology measurements are conducted in order to characterize the viscoelastic wormlike micelle solution. Bright field images show the deformation of the flexible sheet during an unstable oscillation while flow-induced birefringence images highlight the viscoleastic fluid stresses produced in the wake of the flexible sheet.

Une capsule est une goutte de liquide enveloppée par une membrane fine et déformable. Les propriétés mécaniques de la membrane sont essentielles pour le mouvement de la capsule. L’analyse de l’écoulement d’une suspension de capsules dans un canal microfluidique au moyen d’un modèle mécanique est une technique permettant de déterminer les propriétés élastiques de la membrane. Un modèle numérique tridimensionnel a été développé pour résoudre ce problème d’interaction fluide-structure en écoulement confiné. Il couple une méthode des intégrales de frontières pour les écoulements des fluides et une méthode éléments finis pour la déformation de la membrane. Le modèle est utilisé pour étudier l’écoulement d’une capsule initialement sphérique dans des canaux de différentes sections. Dans un canal cylindrique, on montre que l’effet de confinement du canal conduit à la compression de la capsule. Cela engendre la formation de plis sur la membrane autour de l’axe de l’écoulement, phénomène également observé expérimentalement. Dans un canal de section carrée, les effets de la loi constitutive de la membrane, du rapport de taille et du débit d’écoulement sur la déformation de la capsule sont systématiquement étudiés. La comparaison entre les résultats expérimentaux et numériques nous permet de déduire les propriétés mécaniques de la membrane d’une population de capsules artificielles. Ce travail démontre la faisabilité de la mesure de propriétés mécaniques d’une membrane en utilisant une technique microfluidique en canal carré. Il pourrait être étendu par l’étude d’écoulements instationnaires dans un canal de section variable ou avec bifurcations. / A capsule is a liquid droplet enclosed by a thin and deformable membrane. The membrane mechanical properties are critical for the deformation and motion of capsules. The flow of a capsule suspension through a microfluidic channel with dimensions comparable to those of the suspended particles can be used to infer the membrane elastic properties. However a mechanical model of the process is necessary. We present a three-dimensional numerical model to simulate such fluid-structure interaction problem. We use a novel numerical model that couples a boundary integral method for the internal and external fluid flows and a finite element method for the membrane deformation. The model is applied to study the flow of an initially spherical capsule in channels with different cross-sections. In a cylindrical channel with circular cross-section, we show that the confinement effect leads to the compression of the capsule in the hoop direction. The membrane tends to buckle and to fold as observed experimentally. In a microfluidic channel with a square cross-section, the effects of the membrane constitutive law, size ratio and flow strength on the capsule deformation are systematically studied. The comparison between experimental and numerical results allows us to deduce the membrane mechanical properties of a population of artificial capsules. The present work shows that it is possible to measure the membrane mechanical properties by using a microfluidic channel with a square cross-section. It can be extended to unsteady capsule flows in a channel with variable cross-sections or bifurcations.

The highly optimized performance of nature’s creations and biological assemblies has inspired the development of their engineered counter parts that can potentially outperform conventional systems. In particular, bat wings are populated with air flow hair receptors which feedback the information about airflow over their surfaces for enhanced stability and maneuverability during their flight. The hairs in the bat wing membrane play a role in the maneuverability tasks, especially during low-speed flight. The developments of artificial hair sensors (AHS) are inspired by biological hair cells in aerodynamic feedback control designs. Current mathematical models for hair receptors are limited by strict simplifying assumptions of creeping flow hair Reynolds number on AHS fluid-structure interaction (FSI), which may be violated for hair structures integrated on small-scaled Unmanned Aerial Vehicles (UAVs). This study motivates by an outstanding need to understand the dynamic response of hair receptors in flow regimes relevant to bat-scaled UAVs. The dynamic response of the hair receptor within the creeping flow environment is investigated at distinct freestream velocities to extend the applicability of AHS to a wider range of low Reynolds number platforms. Therefore, a threedimensional FSI model coupled with a finite element model using the computational fluid dynamics (CFD) is developed for a hair-structure and multiple hair-structures in the airflow. The Navier-Stokes equations including continuity equation are solved numerically for the CFD model. The grid independence of the FSI solution is studied from the simulations of the hairstructure mesh and flow mesh around the hair sensor. To describe the dynamic response of the hair receptors, the natural frequencies and mode shapes of the hair receptors, computed from the finite element model, are compared with the excitation frequencies in vacuum. This model is described with both the boundary layer effects and effects of inertial forces due to fluid-structure xiv interaction of the hair receptors. For supporting the FSI model, the dynamic response of the hair receptor is also validated considering the Euler-Bernoulli beam theory including the steady and unsteady airflow.

Cette thèse est motivée par un problème aéronautique: le ballottement du carburantdans des réservoirs d’ailes d’avion très flexibles. Les vibrations induites par le couplagedu fluide avec la structure peuvent conduire à des problèmes tels que l’inconfort des passagers,une manoeuvrabilité réduite, voire même provoquer un comportement instable. Cette thèse apour objectif de développer de nouveaux modèles d’interaction fluide-structure, en mettant enoeuvre la théorie des systèmes Hamiltoniens à ports d’interaction (pHs). Le formalisme pHsfournit d’une part un cadre unifié pour la description des systèmes multi-physiques complexeset d’autre part une approche modulaire pour l’interconnexion des sous-systèmes grâce auxports d’interaction. Cette thèse s’intéresse aussi à la conception de contrôleurs à partir desmodèles pHs. Des modèles pHs sont proposés pour les équations de ballottement du liquide en partantdes équations de Saint Venant en 1D et 2D. L’originalité du travail est de donner des modèlespHs pour le ballottement dans des réservoirs en mouvement. Les ports d’interaction sont utiliséspour coupler la dynamique du ballottement à la dynamique d’une poutre contrôlée par desactionneurs piézo-électriques, celle-ci étant préalablement modélisée sous forme pHs. Aprèsl’écriture des équations aux dérivées partielles dans le formalisme pHs, une approximation endimension finie est obtenue en utilisant une méthode pseudo-spectrale géométrique qui conservela structure pHs du modèle continu au niveau discret. La thèse propose plusieurs extensionsde la méthode pseudo-spectrale géométrique, permettant la discrétisation des systèmesavec des opérateurs différentiels du second ordre d’une part et avec un opérateur d’entrée nonborné d’autre part. Des essais expérimentaux ont été effectués sur une structure constituéed’une poutre liée à un réservoir afin d’assurer la validité du modèle pHs du ballottementdu liquide couplé à la poutre flexible, et de valider la méthode pseudo-spectrale de semi-discrétisation.Le modèle pHs a finalement été utilisé pour concevoir un contrôleur basé surla passivité pour réduire les vibrations du système couplé. / This thesis is motivated by an aeronautical issue: the fuel sloshing in tanksof very flexible wings. The vibrations due to these coupled phenomena can lead to problemslike reduced passenger comfort and maneuverability, and even unstable behavior. Thisthesis aims at developing new models of fluid-structure interaction based on the theory ofport-Hamiltonian systems (pHs). The pHs formalism provides a unified framework for thedescription of complex multi-physics systems and a modular approach for the coupling ofsubsystems thanks to interconnection ports. Furthermore, the design of controllers using pHsmodels is also addressed. PHs models are proposed for the equations of liquid sloshing based on 1D and 2D SaintVenant equations and for the equations of structural dynamics. The originality of the workis to give pHs models of sloshing in moving containers. The interconnection ports are used tocouple the sloshing dynamics to the structural dynamics of a beam controlled by piezoelectricactuators. After writing the partial differential equations of the coupled system using thepHs formalism, a finite-dimensional approximation is obtained by using a geometric pseudospectralmethod that preserves the pHs structure of the infinite-dimensional model at thediscrete level. The thesis proposes several extensions of the geometric pseudo-spectral method,allowing the discretization of systems with second-order differential operators and with anunbounded input operator. Experimental tests on a structure made of a beam connected to atank were carried out to validate both the pHs model of liquid sloshing in moving containersand the pseudo-spectral semi-discretization method. The pHs model was finally used to designa passivity-based controller for reducing the vibrations of the coupled system.

Avec la hausse des niveaux de dimensionnement sismique il est devenu nécessaire de limiter les chargements internes dans les structures, notamment en utilisant des dispositifs glissants. Ces dispositifs plafonnent les efforts internes en déclenchant un glissement de la structure. Il devient cependant nécessaire d'estimer l'amplitude des déplacements de corps rigide, notamment pour les structures stockées dans des réservoirs. Dans ce cas, il est nécessaire de prévenir les impacts entre la structure glissante et les bords du réservoir pour contrôler les risques de fuite. Parmi les structures glissantes immergées, on citera les ponts, les structures côtières en maçonnerie, les râteliers de stockage de combustible nucléaire, etc...Les équations de dynamique associées au comportement de ces structures sont non-linéaires et nécessitent l'utilisation de simulations numériques coûteuses en temps de calcul et ne permettant pas de faire des études de sensibilité rapides. On propose donc une méthode de résolution quasi-analytique de ces équations en traitant dans un premier temps, l'évaluation analytique des matrices de masses ajoutées du couplage fluide-structure, dans un second temps, une méthode de résolution quasi-analytique du glissement d'une structure quelconque immergée dans un fluide avec une actualisation de la géométrie de lames d'eau. Les résultats obtenus présentent une bonne adéquation avec des simulations numériques et offrent un temps de calcul quasiment instantané compatible avec une étude paramétrique ou stochastique de ces structures / As the seismic loadings are increasing in accordance to the recent regulations regarding Earthquake design, the use of sliding devices in structures is becoming more common. These devices limitate the internal forces by creating a rigid body sliding. It is then necessary to estimate the global displacement of the structure, especially concerning structures that are immersed in a reservoir. In this case, the displacement must be well estimated in order to prevent impacts between the sliding structure and the boundaries of the reservoir. We can find such structures in : bridges, costal structures in brick and masonry, or in the nuclear industry with the underwater fuel storage racks, ...The governing equations for the behaviour of these structures are non linear and must be solved using time-consuming computer simulations which are not fit for a stochastic study. Our method consists in, firstly, evaluating analytically the added masses of the fluid-structure interaction, secondly, a semi-analytical solving of the governing equations including the updating of the dimensions of the fluid layers surrounding the sliding structure. The results of this new method are in accordance with the numerical simulations and can be obtained in a short time (1 or 2 seconds) which offers the possibility to make a stochastic analysis of the non linear behaviour

Les écoulements chargés en particules sont présents dans de nombreuses applications industrielles telles que le transport de boues ou l’industrie chimique en général. Dans des mélanges constitués de particules solides immergées dans un fluide visqueux, les interactions entre particules jouent un rôle essentiel dans la viscosité globale du mélange.Le phénomène de suspension est causé par des interactions hydrodynamiques à courte distance, connues sous le nom de lubrification. Les forces de lubrification sont généralement sous-estimées en raison de leur nature et de la discrétisation spatiale du problème.Dans cette thèse, nous proposons un modèle de lubrification qui estime les forces et couples hydrodynamiques non résolues par un solveur couplant la résolution des équations de Navier-Stokes incompressible par une méthode de volumes pénalisés, à la résolution de la dynamique des particules par une méthode aux éléments discrets. Les corrections des contraintes hydrodynamiques sont faites localement sur la surface des particules en interaction sans aucune hypothèse sur la forme générale des particules. La version finale du modèle de lubrification proposée peut être utilisée pour des suspensions de particules convexes sans aucune tabulation. La méthode numérique a été validée avec des particules sphériques et des ellipsoïdes, en comparant des simulations à des données expérimentales.Dans le cas de particules sphériques, le modèle de lubrification est aussi précis que les modèles de lubrification existants qui sont limités à ce type de géométrie. La compatibilité du modèle avec des particules convexes a été validée en comparant des simulations,utilisant des ellipsoïdes, à des mesures expérimentales que nous avons réalisées. / Particle-laden flows can be found in many industrial applications such as slurry transport or the chemical industry in general. In mixtures made of solid particles emerged in a viscous fluid, particle interactions play an essential role in the overall mixture viscosity. The suspension phenomenon is caused by short-range hydrodynamic interactions, known as lubrication. Lubrication forces are usually underestimated due to their singularities and the spatial discretization of the numerical schemes. In this thesis, we propose a lubrication model for a coupled volume penalization method and discrete element method solver that estimates the unresolved hydrodynamic forces and torques in incompressible Navier-Stokes flows. Corrections are made locally on the surfaces of the interacting particles without any assumption on the global particle shapes. The final version of the local lubrication model can be used for suspension of convex particles without any tabulations. The numerical method has been validated against experimental data with spherical and ellipsoidal particles. With spherical particles, the lubrication model performs as well as existing numerical models that are limited to this specific particle shape. The model compatibility with convex particles has been validated by comparing simulations using ellipsoids to experimental measurements we made.

Recent technological advances are based on effectively using complex multiphysics concepts. Therefore, there is an ever increasing need for accurate numerical al-gorithms of reduced complexity for solving multiphysics problems. Traditional mesh-based simulation methods depend on a neighbor connectivity information for formulation of operators like derivatives. In large deformation problems, de-pendence on a mesh could prove a limitation in terms of accuracy and cost of preprocessing. Meshless methods obviate the need to construct meshes thus al-lowing simulations involving severe geometric deformations such as breakup of a contiguous domain into multiple fragments. Smoothed Particle Hydrodynamics (SPH) is a meshless particle based Lagrangian numerical method that has the longest continuous history of development ever since it was introduced in 1977. Commensurate with the significant growth in computational power, SPH has been increasingly applied to solve problems of greater complexity in fluid mechanics, solid mechanics, interfacial flows and astrophysics to name a few. The SPH approximation of the continuity and momentum equations govern-ing fluid flow traditionally involves a stiff equation of state relating pressure and density, when applied to incompressible flow problems. Incompressible Smoothed Particle Hydrodynamics (ISPH) is a variant of SPH that replaces this weak com-pressibility approach with a pressure equation that gives a hydrostatic pressure field which ensures a divergence-free (or density invariant) velocity field. The present study explains the development of an ISPH algorithm and its implementa-tion with focus on application to free surface flows, interaction of fluid with rigid bodies and coupling of incompressible fluids with a compressible second phase. Several improvements to the exiting ISPH algorithm are proposed in this study. A semi-analytic free surface model which is more accurate and robust compared to existing algorithms used in ISPH methods is introduced, validated against experi-ments and grid based CFD results. A surface tension model with specific applica-bility to free surfaces is presented and tested using 2D and 3D simulations. Using theoretical arguments, a volume conservation error in existing particle methods in general is demonstrated. A deformation gradient based approach is used to derive a new pressure equation which reduces these errors. The method is ap-plied to both free surface and internal flow problems and is shown to have better volume conservation and therefore reduced density fluctuations. Also, comments on instabilities arising from particle distributions are made and the role of the smoothing functions in such instabilities is discussed. The challenges in imple-menting the ISPH algorithm in a computer code are discussed and the experience of developing an in-house ISPH code is described. A parametric study on water entry of cylinders of different shapes, angular velocity and density is performed and aspects such as surface profiles, impact pressures and penetration velocities are compared. An analysis on the energy transfer between the solid and the fluid is also performed. Low Froude number water entry of a sphere is studied and the impact pressure is compared with the theoretical estimates. The Incompressible SPH formulation, employing the proposed improvements from the study is then coupled with a compressible SPH formulation to perform two phase flow simulations interacting compressible and incompressible fluids. To gain confidence in its applicability, the simulations are compared against the theoretical predication given by the Rayleigh-Plesset equation for the problem of compressible drop in an incompressible fluid.

L'objectif de ce travail est de contribuer à une meilleure prédiction de la rupture des Anévrismes de l'Aorte Abdominale. Pour ce faire, des simulations par éléments finis ont été mises en oeuvre sur des anévrismes modèles dans des conditions proches de la réalité physiopathologique, i.e. en tenant compte de l'anisotropie de la paroi anévrismale, du caractère poreux du thrombus et des Interactions Fluide-Structure. Dans la première partie, une étude statique en l'absence du thrombus a permis de mettre en évidence l'influence de la géométrie de l'anévrisme et du comportement mécanique (isotrope ou anisotrope) de la paroi sur la distribution des contraintes, i.e. la rupture potentielle de l'anévrisme, ainsi que sur l'évolution du module de Peterson. Dans la seconde partie, une modélisation poro-hyperélastique du thrombus est proposée, en s'appuyant sur des données expérimentales de la littérature. La prise en compte de ce comportement et des Interactions Fluide-Structure montre que la pression intra-thrombus demeure du même ordre de grandeur que la pression intra-luminale, conformément à des mesures in vivo réalisées par ailleurs. Enfin, nous montrons que ceci n'est pas en contradiction avec une réduction du risque de rupture potentielle de l'anévrisme. / The aim of this work is to contribute to a better prediction of the Abdominal Aortic Aneurysm rupture (AAA). For that purpose, finite elements simulations have been performed on idealized AAA models under physiopathological like conditions, by taking into account the aneurysmal wall anisotropy, the intra-luminal thrombus porosity and the Fluid-Structure Interactions. In the first part, the influence of the aneurysm geometry and its wall properties (isotropic or anisotropic hyperelasticity) on the wall stress distribution and the Peterson's modulus has been studied in a static analysis and without taking into account the thrombus. In the second part, based on the experimental results from the litterature, a porohyperelastic model has been proposed for the thrombus. By considering such behavior for the thrombus and the Fluid-Structure Interactions, we observe that the intra-thrombus pressure is the same order as the intra-luminal pressure, which is consistent with in vivo measurements. Our results show that despite this unchanged pressure, the maximum wall stress decreases leading to a decrease of the aneurysm potential rupture.

Cette thèse s’inscrit dans le contexte de la simulation et du contrôle des vibrations de structures montées sur ressort qui peuvent apparaître sous l’effet de l’interaction avec l’écoulement de sillage instationnaire. Le contrôle de ce phénomène, appelé vibrations induites par vortex (VIV), est un enjeu critique dans l’optimisation de nombreux systèmes. Une méthode de frontières immergées (IBM) a été intégrée dans l’algorithme PISO du code OpenFOAM, dédié à la simulation d’écoulements fluides incompressibles. La méthode IBM permet une représentation précise de corps fixes ou en mouvement, tout en conservant des maillages structurés conduisant à des algorithmes plus précis et efficaces en termes de performances numériques. Pour calculer la divergence de l’équation de quantité de mouvement dans une boucle PISO et l’interpolation des flux, un calcul hybride orignal a été proposé avec une résolution analytique utilisant l’équation de la fonction noyau des quantités impliquant le terme force de l’IBM (quantités singulières). La méthode a été étendu au formalisme d’écoulements en régimes turbulents. Une loi de paroi a été intégrée permettant de modéliser la couche limite à grand nombre de Reynolds. Le travail de validation a été réalisé au regard des données expérimentales et numériques disponibles dans la littérature pour l’étude d’écoulements autour de cylindres et de sphères, sur une large gamme de nombres de Reynolds. Avec l’objectif de développer des lois de contrôle optimal pour le VIV, basées sur les mécanismes d’instabilité linéaire du système couplé dans le cadre de la théorie du contrôle, un solveur adjoint a été développé et validé. / This thesis is related to the simulation and the control of the vortex induced vibrations phenomenon (VIV), which can result from the fluid structure interactions between an unsteady wake and the body, when the shedding frequency in the wake is close to the natural frequency of the body. The control of VIV is a critical issue when optimizing many systems. An Immersed Boundaries Method (IBM) was implemented into the PISO algorithm as a new library of OpenFOAM, in order to perform reliable simulations of incompressible flows around bluff bodies.To compute the divergence of the momentum equation and the interpolation of the fluxes, an hybrid calculation with an analytical resolution of the quantities involving the force term (singular quantities) has been proposed. The mesh convergence of several errors was shown by means of a manufactured solution, allowing to analyze both the errors irelated to the discretization and to the IBM. The new algorithm was subsequently extended to the RANS and DDES formalism proposed in OpenFOAM for the simulation of turbulent flows. A wall law was integrated into theIBM method to model the boundary layers that develop around the bodies at large Reynolds numbers. Various 2D and 3D well-documented test cases of academic flows around fixed or moving solid bodies (cylinderand sphere) have been simulated and carefully validated against existing data from the literature in a large range of Reynolds numbers. With the objective of developing optimal control laws for VIV, based on the linear instability mechanisms of the coupled system within the framework of the control theory, a new adjoint solver was also developed and validated in OpenFOAM.

[pt] As interações fluido-estrutura são muito comuns na engenharia mecânica e civil porque muitas estruturas, como pontes, plataformas de petróleo, linhas de transmissão ou turbinas eólicas, estão diretamente em contato com um fluido, que pode ser o ar, no caso de vento, ou água, que irá perturbar a estrutura através de ondas. Um papel importante do engenheiro é prevenir a falha da estrutura devido às instabilidades criadas pelas interações fluidoestrutura. Este trabalho apresentará em primeiro lugar todos os conceitos básicos necessários para o estudo de problemas de interação fluido-estrutura. Assim, é realizada uma análise dimensional visando classificar os problemas de fluido-estrutura. A classificação é baseada na velocidade reduzida, e algumas conclusões sobre as conseqüências das interações fluido-estrutura podem ser feitas em termos de estabilidade ou, o que é mais interessante, de instabilidade. De fato, usando modelos simplificados, pode-se mostrar instabilidades estáticas e dinâmicas, induzidas por fluxo, que podem ser críticas para a estrutura. As partes finais do trabalho apresentarão uma estrutura não-linear específica, uma ponte suspensa. Primeiro, a formulação de um modelo simplificado unidimensional é explicada e, em seguida, através de uma discretização por elementos finitos, é realizado um estudo dinâmico. Além disso, algumas conclusões são apresentadas sobre a dinâmica das pontes suspensas. A última parte deste trabalho apresenta um método que foi uma importante fonte de publicação para nós, o método de decomposição regular. / [en] Fluid-structure interactions are very common in mechanical and civil engineering because many structures, as bridges, offshore risers, transmission lines or wind turbines are directly in contact with a fluid, which can be air, which will be source of wind, or water, which will perturb the structure through waves. An important role of the engineer is to prevent structure failure due to instabilities created by the fluid-structure interactions. This work will first present all the basic concepts needed for the study of fluid-structure interaction problems. Thus, a dimensional analysis of those problems is performed and also all the equations governing such cases are presented. Then, thanks to the dimensional analysis made, a classification of problems, namely based on the reduced velocity, can be done and some conclusions concerning the consequences of the fluid-structure interactions can be drawn in terms of stability or, which is more interesting, instability. Indeed, using simplified models one can show static and dynamic flow-induced instabilities that may be critical for the structure. The final parts of the work will present a specific non-linear structure, a suspension bridge. First the formulation of a simplified one-dimensional model is explained and then, through a finite element discretization, a dynamical study is performed. Also, some conclusions are made concerning the dynamic of suspension bridges. The last part of this work presents a method that was an important source of publication for us, the Smooth Decomposition method.