Cette thèse fait partie du projet beyond the sea® qui a pour but de développer la traction par cerf-volant à boudins gonflés (kite) comme système de propulsion auxiliaire des navires. Comme le kite est une structure souple, il est nécessaire de mettre en place une boucle d’interaction fluide-structure pour calculer la géométrie du kite en vol et ses performances aérodynamiques. Un modèle de Ligne Portante 3D Non-Linéaire a été développé pour pouvoir gérer ces ailes non planes, avec des angles de dièdre et de flèche qui varient le long de l’envergure, et également pour pouvoir prendre en compte la non-linéarité du coefficient de portance de la section aérodynamique. Le modèle a été vérifié par des simulations RANSE sur différentes géométries et donne des résultats satisfaisants pour des angles d’incidence et de dérapage variant jusqu’à 15°, avec des différences relatives de quelques pour cent pour l’estimation de la portance globale de l’aile. Les résultats locaux sont aussi correctement estimés, le modèle est capable d’estimer la position du minimum et du maximum de chargement local, selon l’envergure de l’aile, et cela même pour une aile en dérapage. En parallèle, un modèle structure a été développé. L’idée principale du modèle Kite as a Beam est de réduire le kite à un ensemble d’éléments poutre, chacun équivalent à une partie du kite composé d’une section du boudin d’attaque, de deux lattes gonflées et de la canopée correspondante. Le modèle Kite as a Beam a été comparé à un modèle éléments finis complet du kite sur des cas de déplacements élémentaires. Les résultats montrent certaines différences de comportement entre les deux modèles, avec notamment une surestimation de la raideur en torsion pour le modèle Kite as a Beam. Finalement, le modèle Kite as a Beam a été couplé avec la Ligne Portante 3D Non-Linéaire, puis comparé au modèle éléments finis, couplé également avec la Ligne Portante. La réduction du temps de calcul est réellement importante mais les résultats de la comparaison montrent la nécessité de calibrer le modèle Kite as a Beam pour pouvoir retrouver correctement les résultats du modèle éléments finis. / The present thesis is part of the beyond the sea® project which aims to develop tethered kite systems as auxiliary devices for ship propulsion. As a kite is a flexible structure, fluid-structure interaction has to be taken into account to calculate the flying shape and aerodynamic performances of the wing. A 3D Non-Linear Lifting Line model has been developed to deal with non-straight kite wings, with dihedral and sweep angles variable along the span and take into account the non-linearity of the section lift coefficient. The model has been checked with 3D RANSE simulations over various geometries and produces satisfactory results for range of incidence and sideslip up to 15°, with typical relative differences of few percent for the overall lift. The local results are also correctly estimated, the model is able to predict the position of the minimum and maximum loading along the span, even for a wing in sideslip. Simultaneously, a structure model has been developed. The core idea of the Kite as a Beam model is to approximate a Leading Edge Inflatable kite by an assembly of beam elements, equivalent each to a part of the kite composed of a portion of the inflatable leading edge, two inflatable battens and the corresponding canopy. The Kite as a Beam model has been compared to a complete kite Finite Element model over elementary comparison cases. The results show the behaviour differences of the two models, for example the torsion stiffness is globally overestimated by the Kite as a Beam model. Eventually, the Kite as a Beam model coupled with the 3D Non-Linear Lifting Line model is compared to the complete finite element model coupled with the 3D Non-Linear Lifting Line model. The gain in computation time is really significant but the results show the necessity of model calibration if the Kite as a Beam model should be used to predict the results of the complete finite element model.

Le transport de micro-organismes et de fluides biologiques au moyen de cils et flagelles est un phénomène universel que l’on retrouve chez presque tous les êtres vivants. Le but de cette thèse est la modélisation, l’analyse mathématique et la simulation numérique de problèmes d’interaction fluide-structure qui font intervenir des structures actives, capables de se déformer d’elles-mêmes grâce à des contraintes internes, et un fluide à faible nombre de Reynolds, modélisé par les équations de Stokes. Le Chapitre 2 traite de la modélisation de ces structures actives en considérant la loi de Saint Venant-Kirchhoff dans les équations de l’élasticité et en ajoutant un terme d’activité au second tenseur de contraintes de Piola-Kirchhoff. Les équations fluide et structures sont couplées à l’interface fluide-structure et l’étude mathématique d’un problème linéarisé et discrétisé en temps est ensuite réalisée. Une reformulation sous forme d’un problème point-selle est proposée et utilisée pour la simulation numérique du problème. Le Chapitre 3 s’intéresse à l’analyse du problème d’interaction fluide-structure quasi-statique avec une structure active, pour lequel nous montrons l’existence et l’unicité, pour des données petites, d’une solution forte localement en temps. Le Chapitre 4 présente une nouvelle méthode de type domaine fictif (la méthode de prolongement régulier ) pour la résolution numérique de problèmes de transmission. La méthode est d’abord développée pour un problème de transmission de Laplace, puis étendue aux problèmes de transmission de Stokes et d’interaction fluide-structure. / The transport of microorganisms and biological fluids by means of cilia and flagella is an universal phenomenon found in almost all living beings. The aim of this thesis is to model, analyze and simulate mathematical fluid-structure interaction problems involving active structures, capable of deforming themselves through internal stresses, and a low Reynolds number fluid, modeled by Stokes equations. In Chapter 2, these active structures are modeled as elastic materials satisfying Saint Venant-Kirchhoff law for elasticity whose activity comes from the addition of an activity term to the second Piola-Kirchhoff stress tensor. Elasticity and Stokes equations are coupled on the fluid-structure interface and the mathematical study of the linearized problem discretized in time is realized. Then, the problem is formulated as a saddle-point problem which isused for numerical simulations. Chapter 3 focuses on the analysis of a quasi-static fluid-structure with an active structure, for which we show existence and uniqueness, for small data, of a strong solution locally in time. Chapter 4 presents a new fictitious domain method (the smooth extension method) for the numerical resolution of transmission problems. The method is first developed for a Laplace transmission problem and further extended to Stokes transmission and fluid-structure interaction problems.

Ces travaux visent à améliorer la biofidélité des modèles virtuels de l’être Humain. Les statistiques montrent que la tête humaine est fréquemment sujette à des traumatismes cérébraux, des lésions et autres blessures. Une attention particulière sera donc donnée à la modélisation de la tête. Afin de mieux prédire les mécanismes lésionnels de la tête, la biofidélite des modèles doit être améliorée, pour cela les effets du fluide situé à l’intérieur de la tête doivent être pris en compte. Cependant, la modélisation des interactions entre un fluide corporel visqueux et un matériau mou comme le cerveau reste un verrou scientifique. Il est proposé d’étudier en détail la modélisation des interactions fluide-structure entre un fluide et un corps mou. Premièrement, une étude bibliographique détaillée sur les méthodes numériques de modélisation des interactions fluides-structure a permis d’évaluer chacune d’elles et de juger de celle qui est la mieux adaptée pour la résolution de la problématique. Deuxièmement, lors de travaux de thèse précédents, une expérience a été réalisée montrant l’influence du liquide cérébrospinal sur la cinématique du cerveau lors d’un chargement dynamique. Cette expérience est utilisée dans un premier temps pour caractériser numériquement le gel silicone Sylgard 527 utilisé comme substitut de cerveau. Dans un second temps des méthodes de couplage partitionné disponible dans le code commercial LS-Dyna ICFD sont utilisées pour modéliser l’expérience. Bien que les modèles de gel précédemment caractérisés ont été utilisés, la version avec fluide n’a pas pu être modélisée avec succès. Troisièmement, un code de couplage partitionné est donc développé. Il consiste en un middleware écrit en C++ couplant deux codes éprouvés, OpenFOAM et LS-Dyna pour la modélisation du fluide et du solide respectivement. De plus, parce que très peu d’essais expérimentaux utilisables pour la validation de code d’interaction fluide-structure sont disponibles dans la littérature, une expérience permettant cela a été réalisée dans une soufflerie. La comparaison des prédictions numériques avec les résultats expérimentaux est prometteuse et donne des résultats globaux satisfaisants. Les points qui ne peuvent pas être validés nécessitent de plus amples investigations et permettront d’améliorer les techniques de modélisations et le développement du code. / The purpose of this work is to improve the biofidelity of the human body models. The work is focused on the human head as it is one of the most injured part. In order to improve the traumatic brain injury onset and mechanism, the biofidelity of the head models has to be increased, thus, the fluids embedded inside the head has to be taken into account. Nevertheless, the modelling of the interactions occurring between the viscous corporal fluids and the soft matter as the brain remains a challenge. This study intends to investigate the fluid-structure interactions between a soft structure and a fluid. Firstly, in order to found the most relevant methods to solve the problem, a deep literature survey has pointed-out all the numerical methods available nowadays. Secondly, in a previous PhD work an experimental test has been carried-out to demonstrate the influence of the cerebrospinal fluid on the brain kinematics under dynamical load case. On one hand, the Silicon Sylgard 527 gel used as brain substitut has been characterized . Subsequently the partitioned coupling methods available in LS-Dyna ICFD have been assessed to model the experiment. Although, the previous characterized gel model has been used, the experiment has been unsuccessfully completed. Accordingly, it has been decided to develop an in-house coupling code. Thirdly, a partitioned coupling code has been developed. It is a middleware in C++ between two well establishing solvers OpenFOAM and LS-Dyna respectively for the fluid and the solid. Because there is very few experimental tests for the coupling code validation, it has been carried-out in this work a fluid-structure interaction experiment involving a soft plate in a wind channel. This appealing experiment allows the scientific community to validate easily their coupling algorithms. Subsequently, the developed coupling code is used to model the wind channel. The results depict a good overall agreement between the experiment and the simulation. Nonetheless, in order to get validated results further investigation are required mainly about the flow modelling.

La dynamique de filaments flexibles individuels en écoulement visqueux est une étape essentielle pour comprendre et contrôler la rhéologie de nombreux fluides complexes. Cette dynamique sous-tend également une multitude de processus biophysiques allant de la propulsion des micro-organismes aux écoulements intracellulaires. Cette thèse présente des expériences systématiques permettant d’étudier la dynamique de filaments flexibles browniens dans un écoulement visqueux. Nous avons choisi d’utiliser un biopolymère, l’actine, comme système modèle de filaments. Sa longueur typique varie de 1 à 100 μm, il est flexible à ces échelles avec une longueur de persistance de l’ordre de 20μm, à cause de ses petites dimensions, il est soumis aux forces Browniennes avec des fluctuations en flexion, et enfin il peut être marqué en fluorescence. Nous utilisons des dispositifs microfluidiques associés à des systèmes de contrôle d’écoulements, un microscope optique équipé́ avec une platine motorisée pour réaliser des expériences contrôlées permettant de suivre la dynamique des filaments d’actine dans un écoulement de cisaillement pur et dans un écoulement élongationnel. Pour les expériences en cisaillement pur, des simulations reproduisant les conditions expérimentales ont aussi été́ menées en utilisant la théorie des poutres inextensibles d’Euler-Bernoulli et la théorie non locale des corps élancés en présence de fluctuations Browniennes et sont en accord quantitatif avec les résultats expérimentaux. Nous montrons que la dynamique des filaments dans ce système est principalement régie par le nombre élasto-visqueux, nombre sans dimension comparant les forces de trainée visqueuses aux forces de flexion élastiques, les fluctuations thermiques ne jouant qu’un rôle secondaire. Nous présentons une caractérisation complète des différents modes de déformation subies par le filament pendant une rotation ainsi que des transitions entre les différents modes. Dans la géométrie élongationnelle, nous avons choisi un canal hyperbolique optimisé pour permettre de longs temps de résidence sous taux de déformation constant. Nous avons observé́ directement la suppression des fluctuations transverse dans la partie extensionnelle tandis que nous observons, dans la partie compressive la formation de structures hélocoïdales tridimensionnelles après le flambage du filament. Pour finir, ce manuscrit de thèse décrit des développements expérimentaux permettant de fabriquer des suspensions de filaments d’actine relativement monodisperse en taille ainsi que des résultats préliminaires sur des effets rhéofluidifiants. Au bilan, les résultats présentes dans ce manuscrit pose les premières pierres de travaux futurs en direction de l’étude de la dynamique de ces filaments dans des écoulements plus complexes comme des écoulements de Poiseuille ou oscillants. Ils permettent aussi d’envisager des études sur le lien entre déformations de particules et propriétés des suspensions diluées d’objets flexibles et Brownien, lien encore peu étudié du point de vue expérimental. / The dynamics of individual flexible filament in a viscous flow is the key to deciphering the rheolog- ical behavior of many complex fluids and soft materials. It also underlies a wealth of biophysical processes from flagellar propulsion to intracellular streaming. This thesis presents systematic exper- iments to investigate the dynamics of flexible and Brownian filaments in viscous flows. Biopolymer actin has been chosen to be our experimental model filament: its typical length can be varied from 1 to 100 μm, it is flexible at these dimensions with a persistence length in the order of 20μm, it is Brow- nian due to its small diameter with bending fluctuations and it can be labelled by fluorescent dye. Microfluidic channels and flow control systems are combined to optical microscope with automated stage to carry out well-controlled experiments on the diverse dynamics of actin filaments in shear flow and pure straining flow. In shear flow, simulations matching the experimental conditions have also been performed using inextensible Euler-Bernoulli beam theory and non-local slender body hy- drodynamics in the presence of thermal fluctuations and agree quantitatively with the experimental results. We demonstrate that filament dynamics in this flow geometry is primarily governed by a dimension- less elasto-viscous number comparing viscous forces to elastic forces with thermal fluctuations only playing a secondary role. We present a complete characterization of the different modes of defor- mation undergone by the filament while rotating as well as of the transitions between these different modes. In pure straining flow, we opt to use an optimized hyperbolic channel to allow long resi- dence time at constant strain rate to be applied. We directly observe the suppression of transverse fluctuations in the extensional part of the hyperbolic channel while we observe, in the compressive part of the flow, the formation of three dimensional helical structures subsequent to the initial buck- ling of the filament. Finally, this thesis manuscript also reports on experimental developments to fabricate suspensions of actin filament with a narrow distribution of lengths and on preliminary re- sults on shear-thinning effects. All together the results presented here pave the way of future studies towards the understanding of filament dynamics in more complex flows, as Poiseuille flows or oscil- latory flows, as well as towards establishing the link between filament deformations and rheological response in dilute suspensions of flexible Brownian filaments, which remains nearly unexplored from an experimental point of view.

Les travaux de cette thèse s’inscrivent dans le cadre d’une collaboration CIFRE entre Naval Group, le Laboratoire d’Hydrodynamique, d’Energétique et d’Environnement Atmosphérique (LHEEA) de l’Ecole Centrale de Nantes et de l’Institut de Génie des Matériaux (GeM) de l’ICAM de Nantes, sur la problématique de développement d’hélices marines plus efficaces, plus discrètes, et plus respectueuses de l’environnement. Une des solutions passe par le développement de structures composites plus légères et plus flexibles, capables de se déformer passivement pour s’auto-adapter à l’écoulement incident grâce à leurs propriétés spécifiques de couplage flexion-torsion. En lien direct avec cette problématique, nous avons réalisé les travaux de thèse en deux temps. Nous avons dans un premier temps mis en place, montré la faisabilité et validé une méthode de couplage numérique fluidestructure implicite fort entre les codes commerciaux de CFD Starccm+ et de CSD Abaqus pour un cas d'application issu de résultats expérimentaux disponibles dans la littérature sur deux hydrofoils flexibles déformables. Puis dans un second temps nous avons développé, réalisé et testé, dans le tunnel hydrodynamique de l’Ecole Navale, un profil portant composite spécifiquement conçu pour s’approcher du comportement d’une hélice. L’étude expérimentale et numérique de cette configuration nous a permis i) d’éprouver l'utilisation en milieu académique et industriel de nouvelles méthodes expérimentales d'instrumentation d'une pièce composite par insertion directe de fibres optiques dans les plis de composite, et d'une méthode mixte numérique - expérimentale de calibration fine d'un modèle numérique structure ; ii) d’apporter un éclairage sur la physique de l'interaction fluidestructure se produisant sur une surface portante composite ; et iii) de préciser les limitations actuelles concernant la diffusion en milieu industriel de cette méthode de couplage numérique fluide-structure. / This Ph.D is sponsored by the French company Naval Group in collaboration with LHEEA Laboratory from Ecole Centrale de Nantes and GeM Institute from ICAM de Nantes, and deals with the development of new composite marine propellers with improved efficiency, improved acoustic discretion and more environment-friendly. One of the key solutions lies in the application of composite materials to marine structures, in order to benefit from their reduced weight, increased flexibility and bend-twist coupling capacity. Indeed, the latter enables the shape-adaptability of the structure to passively adapt to the incoming flow. To meet this challenge, we first set-up a tightly coupled numerical fluid-structure method using two commercial CFD (Starccm+) and CSD (Abaqus) solvers on two flexible hydrofoils and we validated this method against experimental results available in the literature. Second, we specifically developed a composite hydrofoil to behave closely like a marine propeller and tested it in the hydrodynamic tunnel of the Ecole Navale. Thanks to the combined experimental and numerical analysis of this composite hydrofoil we reached the following conclusions: i) we helped demonstrate the industrial application of a state-of-the-art strain measurement technique using optical fibers directly embedded within the composite plies, ii) we provided some insights into the physics of the fluid-structure interaction occurring on composite hydrofoils and iii) we presented the current limitations of this coupled numerical fluid-structure method relatively to its industrial application.

L’éolien offshore est l’énergie marine la plus avancée et utilisée dans le monde. Afin d’accroître l’énergie extraite du vent, les dimensions des éoliennes deviennent plus importantes et les parcs éoliens sont installées de plus en plus loin des côtes, où les mers sont plus agitées et les vents plus forts. De fait, les opérations marines sont plus complexes et plus chères et les fenêtres météo sont écourtées et se raréfient. Dans le cadre de cette thèse, un logiciel de simulation numérique des opérations marines est développé, en particulier pour des applications de descentes et de remontées de colis lourds. L’Algorithme aux Corps Rigides Composites, implémenté dans le logiciel InWave, est utilisé pour modéliser le système multicorps. Un modèle de câble et de treuil est développé, suivant la théorie multicorps utilisée, et comparé à la théorie câble classique dite « lumped mass ». Les efforts hydrodynamiques ainsi que les interactions hydrodynamiques sont modélisés par une théorie potentiel instationnaire satisfaisant l’hypothèse de faible perturbation, dite « weak-scatterer ». L’approche « weak-scatterer » du logiciel WS_CN est étendue aux simulations multi-flotteurs et validée par comparaison avec des données expérimentales. InWave et WS_CN sont couplés afin de résoudre l’interaction houle-structure pour des systèmes multicorps articulés en mer. Un couplage fort est adopté pour sa robustesse. L’équation de couplage est établie et validée via des comparaisons avec WS_CN. Le logiciel ainsi crée se nomme InWaveS_CN et utilise un code d’intégration en Python. Une nouvelle stratégie de maillage, basée sur un algorithme de découpe de maillages et une méthode par avance de front, est développée dans WS_CN. Enfin, des essais en bassin d’une opération de redressement ont été menés à l’ECN. La comparaison entre les simulations numériques et les données expérimentales offre une première et prometteuse validation d’InWaveS_CN. / Offshore wind represents the most advanced and used marine energy in the world. To increase the wind power extraction, turbines grow in size and wind farms are installed further offshore in presence of rough seas and strong winds. Marine operations become more challenging and expensive, weather windows are shorter and less frequent. This PhD work focuses on the development of a numerical tool to simulate marine operations with consistency, in particular lowering and lifting operations. The Composite-Rigid-Body Algorithm, implemented in the numerical tool InWave, is used to model multibody systems. A cable model and a winch model are developed following this multibody approach and compared to the classical low-order lumped mass theory. Hydrodynamic loads and hydrodynamic interactions are simulated using an unsteady potential flow theory based on the weakscatterer hypothesis, implemented in the numerical tool WS_CN. This approach is extended to multibody simulations and validated with comparisons to experimental data. InWave and WS_CN are coupled to solve wavestructure interaction for articulated multibody systems with large relative motions in waves. A tight coupling is selected for its robustness. The coupling equation is derived and validated from comparisons with WS_CN. This leads to the creation of a new numerical tool, InWaveS_CN, using Python as glue code language. A new mesh strategy, based on the coupling between a panel cutting method and an advance front method, is developed in WS_CN. Experiments of an upending operation were conducted at Ecole Centrale de Nantes. The comparison between the numerical simulations and the experimental data leads to a first and promising validation of InWaveS_CN.

Les tubes des générateurs de vapeur des centrales nucléaires vibrent sous l'effet d'écoulement eau/vapeur. Pour appréhender ce phénomène et le comprendre, des expériences à échelles réduites sont réalisées. La simulation numérique a montré son habilité à reproduire l'interaction fluide-structure sur ce type de géométrie pour des écoulements monophasiques. L'objectif est désormais de faire de même en écoulement diphasique et de caractériser les propriétés physiques du mélange liquide/gaz influant sur la vibration.Pour se faire, un code CFD avec une approche bi-fluide est utilisé. Une méthode dite de "Discrete forcing" est implémentée pour permettre le mouvement imposé de corps solides au sein d'un écoulement à plusieurs phases. Celle-ci est alros validée sur des cas simples et intégraux avec une comparaison systématique à des résultats expérimentaux ou théoriques.En se basant sur un algorithme implicite existant dans la littérature, un couplage fluide-structure utilisant cette méthode de suivi d'interface est implémenté. Validé sur des cas monophasiques et diphasiques, ce couplage offre désormais la possibilité de déplacer un solide en fonction des forces fluides diphasiques qui lui sont appliquées.Les différentes méthodes numériques présentes dans NEPTUNE_CFD sont ensuite évaluées pour un écoulement fréon/fréon au travers d'un faisceau de tubes inclinés. La nécessité d'utiliser des modèles dit "multi-régime" est mis en avant.Afin de déterminer l'influence sur l'écoulement des différentes propriétés physiques d'un mélange diphasique, plusieurs cas simples sont réalisés.Finalement, l'application industrielle cible, un écoulement eau/fréon dans un faisceau de tubes à pas carré, est simulée et comparée à un écoulement en conditions réelles (eau/vapeur à 70 bar). Les vibrations induites par écoulement monophasique puis diphasique sont correctement reproduites sur des cas dit de "faisabilité". / In nuclear power plants, steam generator tubes vibrate because of steam/water cross-flows. In order to understant this phenomenon, reduced-scale experiments are performed. Numerical simulations have shown their ability to accurately reproduce the vibration induced by a single phase flow in a tube bundle. The aim of the present work is to do the same with two-phase flow and to characterize the effect of the mixture physical properties on vibration.To do so, a CFD code based on a two-fluid approach is used. A "discrete forcing" method is implemented in order to allow solid body motion in a two-phase flow. The validation is performed with simple and industrial cases using experimental and theoretical results.Using an existing implicit algorithm, a fluid-structure coupling based on the developed interface tracking method is implemented. Validated for single and two-phase flows, it is now possible to have solid motion induced by fluid forces.The different numerical models dedicated to two-phase flows are then evaluated on a freon/freon flow across an inclined tube bundle. The use of a multi-regime model is required. In order to investigate the role of the different physical properties on the vibration, three simple studies are performed.Finally, the industrial application, a freon/water flow across a square pitch tube bundle, is performed. First, it is compared to a steam/water flow in order to characterize the discrepancies when we are using a modeling mixture. Then, the vibration induced by single- and two-phase flows is reproduced by the developed method on feasibility test cases.

La modélisation des problèmes d'interactions sol-structure et fluide-structure couvre plusieurs domaines de recherche très actifs qui traitent une multitude d'aspects tels que la géométrie non bornée du sol et dans certains cas du fluide stocké, les effets dissipatifs visqueux et radiatifs, l'application du chargement sismique, le choix des variables de base, les propriétés algébriques des systèmes d'équations résultant du couplage,... etc. Dans le présent travail, différents modèles numériques de couplage sol-structure et fluide-structure ont été examinés. Les limites de troncature géométrique du sol et du fluide on été traitées avec des éléments infinis dont les performances ont été comparées à celles des conditions de radiations. Le problème de vibrations libres couplées des systèmes fluide-structure a été résolu en introduisant de nouvelles techniques de symétrisation efficaces. De plus, une nouvelle formulation symétrique en éléments de frontière a été proposée. Cette formulation permet de produire une matrice symétrique définie positive et aboutit ainsi à un système algébrique similaire à celui qui découle de la discrétisation en éléments finis. La matrice bâtie dite "raideur équivalente" peut facilement être assemblée ou couplée avec les matrices de la formulation en éléments finis. Toutes les applications qui ont servi soit à comparer des modèles soit à valider les programmes développés, ont été effectué es dans le cas des barrages poids en béton. Ce cas constitue un problème de couplage fluide-sol-structure typique / Modeling fluid-structure and soil-structure interaction problems covers several research fields dealing with multiple aspects such as : unbounded geometry of soil media and in some cases of retained fluid, viscous and radiation dissipative effects, application of seismic loading, choice of the basic variables, algebraic properties of the resulting coupled system, ... etc. In this work, different numerical models of soil-structure and fluid-structure coupling have been studied. The truncation boundaries of the soil and of the fluid domains have been considered by using infinite elements whose performances were compared to those of radiation conditions. The problem of the free coupled vibrations of the fluid-structure systems has been solved by introducing efficient symmetrization techniques. In addition, a new symmetric boundary element formulation is proposed. It allows to produce a positive definite and symmetric matrix and therefore to conduct to an algebraic system similar to the one obtained from finite element discretization. The produced matrix called "equivalent stiffness matrix" can easily be assembled or coupled to finite element matrices. All of the applications which have been done in order to compare models or to validate the developed programs were built in the case of concrete gravity dams, which constitute typical coupled fluid- soil-structure problems

L'Interaction Fluide-Structure (IFS) décrit une classe très générale de problème physique, ce qui explique la nécessité de développer une méthode numérique capable de simuler le problème FSI. Pour cette raison, un solveur IFS est développé qui peut traiter un écoulement de fluide incompressible en interaction avec des structures différente: élastique ou rigide. Dans cet aspect, le solveur peut couvrir une large gamme d'applications.La méthode proposée est développée dans le cadre d'une formulation monolithique dans un contexte Eulérien. Cette méthode consiste à considérer un seul maillage et résoudre un seul système d'équations avec des propriétés matérielles différentes. La fonction distance permet de définir la position et l'interface de tous les objets à l'intérieur du domaine et de fournir les propriétés physiques pour chaque sous-domaine. L'adaptation de maillage anisotrope basé sur la variation de la fonction distance est ensuite appliquée pour assurer une capture précise des discontinuités à l'interface fluide-solide.La formulation monolithique est assurée par l'ajout d'un tenseur supplémentaire dans les équations de Navier-Stokes. Ce tenseur provient de la présence de la structure dans le fluide. Le système est résolu en utilisant une méthode élément fini et stabilisé suivant la formulation variationnelle multiéchelle. Cette formulation consiste à décomposer les champs de vitesse et pression en grande et petite échelles. La particularité de l'approche proposée réside dans l'enrichissement du tenseur de l'extra contraint.La première application est la simulation IFS avec un corps rigide. Le corps rigide est décrit en imposant une valeur nul du tenseur des déformations, et le mouvement est obtenu par la résolution du mouvement de corps rigide. Nous évaluons le comportement et la précision de la formulation proposée dans la simulation des exemples 2D et 3D. Les résultats sont comparés avec la littérature et montrent que la méthode développée est stable et précise.La seconde application est la simulation IFS avec un corps élastique. Dans ce cas, une équation supplémentaire est ajoutée au système précédent qui permet de résoudre le champ de déplacement. Et la contrainte de rigidité est remplacée par la loi de comportement du corps élastique. La déformation et le mouvement du corps élastique sont réalisés en résolvant l'équation de convection de la Level-Set. Nous illustrons la flexibilité de la formulation proposée par des exemples 2D. / Numerical simulations of fluid-structure interaction (FSI) are of first interest in numerous industrial problems: aeronautics, heat treatments, aerodynamic, bioengineering... Because of the high complexity of such problems, analytical study is in general not sufficient to understand and solve them. FSI simulations are then nowadays the focus of numerous investigations, and various approaches are proposed to treat them. We propose in this thesis a novel monolithic approach to deal with the interaction between an incompressible fluid flow and rigid/ elastic material. This method consists in considering a single grid and solving one set of equations with different material properties. A distance function enables to define the position and the interface of any objects with complex shapes inside the volume and to provide heterogeneous physical properties for each subdomain. Different anisotropic mesh adaptation algorithms based on the variations of the distance function or on using error estimators are used to ensure an accurate capture of the discontinuities at the fluid-solid interface. The monolithic formulation is insured by adding an extra-stress tensor in the Navier-Stokes equations coming from the presence of the structure in the fluid. The system is then solved using a finite element Variational MultiScale (VMS) method, which consists of decomposition, for both the velocity and the pressure fields, into coarse/resolved scales and fine/unresolved scales. The distinctive feature of the proposed approach resides in the efficient enrichment of the extra constraint. In the first part of the thesis, we use the proposed approach to assess its accuracy and ability to deal with fluid-rigid interaction. The rigid body is prescribed under the constraint of imposing the nullity of the strain tensor, and its movement is achieved by solving the rigid body motion. Several test case, in 2D and 3D with simple and complex geometries are presented. Results are compared with existing ones in the literature showing good stability and accuracy on unstructured and adapted meshes. In the second, we present different routes and an extension of the approach to deal with elastic body. In this case, an additional equation is added to the previous system to solve the displacement field. And the rigidity constraint is replaced with a corresponding behaviour law of the material. The elastic deformation and motion are captured using a convected level-set method. We present several 2D numerical tests, which is considered as classical benchmarks in the literature, and discuss their results.

Ce travail de thèse s’appuie sur une étude expérimentale en tube à choc, plus précisément, c’est une approche expérimentale de l'étude de l'interaction fluide-structure. Considérons un solide indéformable auquel on laisse un degré de liberté en rotation autour d'un axe. Cette structure ferme un espace clos. Si le contenu de l'espace clos subit le passage d'une onde de choc, ce solide va être mis mouvement et tourner autour de son axe. Concrètement, l'onde de choc va augmenter les caractéristiques physiques, en particulier sa pression, du fluide en contact avec la face impactée de cette porte. La face opposée de la porte ne subissant pas ou que très peu l'influence de l'onde de choc, une seule de ses faces est soumise à la surpression. Au moment de l'impact, le déséquilibre ainsi créé impose une action mécanique sur la porte qui va la faire accélérer et tourner autour de son axe de rotation. Jusqu'à ce stade tout est relativement simple. La difficulté intervient à l'instant où la porte commence à s'ouvrir, car les frontières du volume dans lequel le fluide évolue sont modifiées. Des fuites apparaissent et le gaz qui était maintenu dans un volume clos peut maintenant s'écouler vers un milieu libre. Une communication entre les gaz agissant de chaque coté de la porte est créée modifiant leurs propriétés et par conséquent la pression agissant sur chaque côté de la porte. Les actions mécaniques qui s'appliquent sur la porte ne sont plus les mêmes, et par conséquent l'accélération que la porte subit aussi. Au fur et à mesure que la porte change de position, le problème fluide continue d'être modifié et change en retour son action sur la porte. Cette interaction perdure soit jusqu'à ce que les limites du problème cessent d'être modifiées, la porte ne peut plus bouger, ou bien lorsque les actions mécaniques agissant sur la porte s'équilibrent, les fluides de chaque côté de la porte étant dans le même état physique. Le travail présenté ici est une étude des paramètres du fluide ou du solide en mouvement qui sont les acteurs de la loi comportementale gérant ce système complexe. Pour ce faire, nous avons réalisé une maquette expérimentale mettant en action la physique que nous venons de décrire et nous l'avons adaptée à un tube à choc. En éprouvant de nombreuses configurations expérimentales, nous avons pu déterminer comment l'écoulement interne d'un tube à choc évolue lorsqu'il est plus ou moins ouvert à son extrémité. Comment une porte fermée réagit-elle à l'impact d'une onde de choc et quelles en sont les conséquences sur l'évolution des fluides mis en jeu? Quelles sont les conséquences d'une position différente de la porte au moment de l'impact avec l'onde de choc? Ou encore, quel rôle joue l'intensité de l'onde de choc incidente ou l'inertie de la porte sur toute cette dynamique? / This thesis is based on an experimental study carried out in shock tube; in particular, this is an experimental approach to the study of fluid-structure interaction. Consider a rigid body which is allowed to rotate only around an axis and which closes a confined space. If a shock wave crosses the content of the confined space, the body will accelerate and rotate around its axis. Specifically, the shock wave will increase the physical characteristics, especially its pressure, of the fluid acting on the impacted face of the door. The opposite side of the door is not influenced by the incident shock wave, only one of its faces is subjected to overpressure. Following the first impact, the resulting imbalance imposes a mechanical action on the door that will increase its speed and make it turn around its rotation axis. The difficulty comes when the door begins to open: the volume boundaries in which the fluid is contained are modified. Leaks occur and the gas kept in this closed volume can now flow to the atmosphere. Communication between the gas acting on each side of the door is created modifying their properties and consequently the pressure acting on each side of the door.The mechanical actions that apply to the door are no more the same with time, and therefore the acceleration of the door is changing. As the door moves, the fluid problem continues to be changed and in turn it changes its action on the door. This interaction process continues until either the limits of the problem ceases to be changed, the door cannot move, or when the mechanical actions acting on the door are in equilibrium, fluids on each side of the door are in the same physical state. The presented work is a study of the parameters of the fluid or the solid motion which are main actors in the behavioral law managing this complex system. In this aim, we designed an experimental device involving the physics that we have described and we have adapted it to a shock tube. Testing many experimental configurations, we could determine how the internal flow of a shock tube evolves when the end of this shock tube is more or less open.How a closed door reacts to the impact of a shock wave and what are the implications for the evolution of the involved fluids? What are the consequences of a different position of the door at the instant of the impact with the incident shock wave? What role plays the intensity of the incident shock wave or the inertia of the door on this dynamic?