Les micronageurs et parmi eux les microangeurs artificiels sont en général, limités à exister dans des écoulements dominés par des forces visqueses. Ces écoulements sont caractérisés par un bas nombre de Reynolds (Re). Cela impacte la stratégie de nage et plus particulièrement les séquences de forme possibles, qui doivent nécessairement être non-réciproques dans l'espace de déformation pour espérer induire un déplacement net non-nul. De plus, due aux forts effets de traînée, les vitesses de nage sont limités à des valeurs faibles.Dans cette thèse, on examine la possibilité d'utiliser un mécanisme de nage basé sur l'instabilité de flambage d'une sphère creuse. Cette instabilité est provoquée en soumettant la sphère à une onde de pression. La particularité de ce mécanisme est qu'il satisfait par construction la condition nécessaire de nage à bas Reynolds exposée précédemment. De plus, la rapidité de la déformation lors de l'instabilité pousse à prévoir l'apparition d'effets inertiels, et ce même à l'échelle microscopique.Une étude expérimentale a été conduite à l'échelle macroscopique dans le but de comprendre la dynamique de l'instabilité et son impact sur le fluide qui entoure la coque creuse. Ces expériences nous permettent de montrer qu'un déplacement net non-nul est produit pour tous les régimes d'écoulements.On met en évidence le rôle de paramètres géométriques, des propriétés du matériau composant la coque creuse et de la rhéologie du fluide sur l'efficacité de la nage.On montre l'existence d'un optimum de déplacement net pour des valeurs intermédiaires du nombre de Reynolds. Pour expliquer cela, on se sert de mesures de PIV résolues temporellement pour mettre en évidence la présence d'effets d'histoire non-triviaux qui augmentent le déplacement net.On dérive un simple modèle en se basant sur les observations expérimentales pour montrer que ce régime optimal de nage est atteignable pour des sphères microscopiques, ceci est possible grâce l'activation rapide de l'instabilité. Cette propriété permet aussi une excitation à haute fréquence en utilisant des ultrasons. Une étude d'échelle nous permet de prédire une vitesse de nage de 1 cm/s pour un micro-robot contrôlé à distance. Cet ordre de grandeur de vitesse est idéal pour des applications biologiques comme la distribution ciblée de médicaments. / Microswimmers, and among them aspirant microrobots, are generally bound to cope with flows where viscous forces are dominant, characterized by a low Reynolds number (Re). This implies constraints on the possible sequences of body motion, which have to be nonreciprocal. Furthermore, the presence of a strong drag limits the range of resulting velocities.Here, we propose a swimming mechanism which uses the buckling instability triggered by pressure waves to propel a spherical hollow shell. The particularity of this mechanism is that it fulfills naturally the necessary condition of swimming at low Re. In addition, the swiftness of the instability might produce inertial effects even at the microscopic scale.With a macroscopic experimental model we show that a net displacement is produced at all Re regimes. We put in evidence the role of geometrical parameters, shell material properties and rheology of the surrounding fluid on the swimming efficiency.An optimal displacement is reached at intermediate Re. Using time-resolved PIV measurements, we explain that non-trivial history effects take place during the instability and enhance net displacement.Using a simple model, derived from the study of shell dynamics, we show that due to the fast activation induced by the instability, this regime is reachable by microscopic shells. The rapid dynamics would also allow high frequency excitation with standard traveling ultrasonic waves. Scale considerations predict a swimming velocity of order 1 cm/s for a remote controlled microrobot, a suitable value for biological applications such as drug delivery.

Le comportement des matériaux multiphasiques couvre une multitude de phénomènes qui suscitent un grand intérêt dans le domaine scientifique et professionnel. Les propriétés mécaniques de ces types de matériau trouvent leur origine dans les phases dont ils sont composés, leur distribution et interaction. Un nouveau modèle hydrodynamique couplé est présenté dans ce travail de thèse, à appliquer à l'analyse de l'hydrodynamique des milieux granulaires saturés. Le modèle associe la méthode des éléments discrets (DEM) pour la modélisation de la phase solide, avec une formulation en volumes finis, à l'échelle des pores (PFV), du problème de l'écoulement. Une importance particulière est donné à la description de l'interaction entre les phases, avec la détermination des forces fluides à appliquer sur chacune des particule, tout en assurant un coût de calcul abordable, qui permet la modélisation de plusieurs milliers des particules en trois dimensions. Le milieux est considéré saturé par un fluide incompressible. Les pores et leur connectivité est basée sur une triangulation régulière des assemblages. L'analogie de cette formulation avec la théorie classique de Biot est présenté. Le modèle est validé par la comparaison des résultats numériques obtenus pour un problème de consolidation d'un sol granulaire avec la solution analytique de Terzaghi. Une approche pour analyser l'hydrodynamique d'un sédiment granulaire est finalement présenté. La reproduction du phénomène de liquéfaction d'un sol est également présentée. / The behaviour of multiphase materials covers a wide range of phenomena of interest to both scientists and engineers. The mechanical properties of these materials originate from all component phases, their distribution and interaction. A new coupled hydromechanical model is presented in this work, to be applied to the analysis of the hydrodynamics of saturated granular media. The model associates the discrete element method (DEM) for the solid phase, and a pore-scale finite volume (PFV) formulation of the flow problem. The emphasis of this model is, on one hand, the microscopic description of the interaction between phases, with the determination of the forces applied on solid particles by the fluid; on the other hand, the model involve affordable computational costs, that allow the simulation of thousands of particles in three dimensional models. The medium is assumed to be saturated of an incompressible fluid. Pore bodies and their connections are defined locally through a regular triangulation of the packings. The analogy of the DEM-PFV model and the classic Biot's theory of poroelasticity is discussed. The model is validated through comparison of the numerical result of a soil consolidation problem with the Terzaghi's analytical solution. An approach to analyze the hydrodynamic of a granular sediment is finally presented. The reproduction of the phenomenon of soil liquefaction is analysed and discussed.

Afin d’augmenter le rendement et la puissance spécifique des turbines à gaz, la température en entrée de la turbine haute pression a été continuellement augmentée, au point d’atteindre une valeur de l’ordre de la température de fusion des matériaux. La capacité à prédire l’écoulement (température en particulier) dans la turbine est donc un élément essentiel pour la conception des turbines à gaz, particulièrement pour celle qui ne sont pas refroidies. Toutefois, cette tâche est rendue difficile par l’extrême complexité de l’environnement dans lequel la turbine évolue (interaction avec la chambre de combustion, effets technologiques, couplage thermique fluide / solide, etc.).L’objectif de ce travail de thèse est de proposer, d’une part une stratégie de simulation numérique afin de prédire au mieux les champs aérothermiques dans une turbine haute-pression, d’autre part une méthode d’analyse permettant de quantifier l’impact des différents éléments de l’environnement sur les performances. Pour répondre à ces objectifs, des simulations instationnaires d’un étage de turbine ont été réalisées avec le code elsA, prenant en compte les effets technologiques (baignoire, refroidissement externe grâce aux évents du distributeur et aux cavités moyeu et carter du rotor qui sont alimentés par le système d’air secondaire, congés de raccordement) et les transferts thermiques conjugués. Une analyse de la production locale d’entropie a aussi été menée afin de comparer les performances aérodynamiques pour plusieurs niveaux d’approximations dans la définition de la turbine (prise en compte ou non de tel ou tel effet technologique) et de localiser l’origine de ces différences. L’analyse des résultats a montré que la hauteur de jeu et le système de refroidissement externe ont l’impact le plus significatif sur les performances aérodynamiques de la turbine. La température de paroi de la pale de rotor est de son côté fortement influencée par l’écoulement issu de la chambre de combustion, le refroidissement externe et le couplage thermique fluide / solide. Ce travail est un premier pas vers la réalisation de simulations totalement intégrées qui doivent permettre d’améliorer la précision des conceptions. / In order to increase the thermodynamic efficiency of gas turbine engines, the high-pressure turbine inlet temperature has been continually increased up to reach levels of the order of magnitude of the vanes and blades melting temperatures. The ability of predicting the flow through the turbine(especially the temperature) is a key point for the design of gas turbines, especially for the uncooled ones. However, this is challenging because of the complex environment that interacts with the turbine(hot-streak migration, technological details, fluid/solid thermal coupling . . .).The aim of this work is to develop a strategy based on CFD in order to predict aerothermal fields in a high-pressure turbine as well as an analysis allowing to quantify the impact of the environment on the turbine performances. To achieve these goals, the elsA code has been used to perform unsteady simulation of a turbine stage, taking into account technological details (squealer tip, external cooling thanks to vane trailing edge cooling and rotor hub and shroud cooling cavities, fillets) and conjugate heat transfers. An analysis of the local entropy production rate was also used to compare the performances resulting from different modeling of the turbine (taking into account or not one particular technological detail) and to localize the origin of these discrepancies. The results show that the tip height and the external cooling have the greatest impact on the turbine performances. The rotor blade temperature is mainly affected by the flow coming from the combustion chamber, the external cooling and the fluid / solid thermal coupling. This work is the first step towards the realization of integrated simulations allowing to improve the accuracy of design.

La suffusion est un cas particulier d'érosion interne qui apparait dans les ouvrages hydrauliques. Elle se caractérise par le détachement et le transport des particules fines à travers les gros grains sous l'action d'un écoulement hydraulique en laissant derrière un squelette granulaire dont les caractéristiques en termes de densité et d'arrangement géométrique granulaire ont été changées. De telles modifications dans la micro-structure du sol peuvent conduire à des déformations à l'échelle macroscopique et peuvent influencer le comportement mécanique du sol. Ce travail a été consacré à l'étude du mécanisme de suffusion et à son impact sur les propriétés mécaniques des sols non-cohésifs. Pour atteindre cet objectif, deux approches ont été suivies dans ce travail: l'une numérique et l'autre expérimentale. Un modèle numérique discret a été défini pour décrire quantitativement le comportement macroscopique du sol et analyser la micro-structure des assemblages granulaires. Il est basé sur l'utilisation de particules sphériques et la résistance au roulement aux contacts, ainsi que sur une nouvelle méthode de compactage des échantillons pour simuler celle utilisée au laboratoire, et atteindre une large gamme de densités initiales. Le modèle a été validé par comparaison de résultats numériques et de données expérimentales. Ce modèle a d'abord été appliqué sur des assemblages granulaires avec des teneur en fines différentes afin d'étudier le rôle des particules fines pour une micro-structure de sol donnée, sans tenir compte d'un processus d'érosion. Il a été montré qu'il existe un pourcentage de fines, en dessous de la teneur seuil, où des particules fines peuvent commencer à participer au transfert de force, et que, si la suffusion s'initie et mobilise ces particules, ceci pourra affecter le comportement macroscopique du sol. Par la suite, une cinétique simplifiée d'extraction du grain a été proposée pour décrire le processus de suffusion. Elle est basée sur une approche de couplage fluide-solide partielle. L'intérêt de cette procédure d'extraction est que, d'une part, elle tient en compte des critères géométriques et hydrauliques complexes; d'autre part, le modèle induit un coût de calcul raisonnable. Cette procédure a été validée sur la base d'essais expérimentaux de suffusion effectués avec un dispositif nouvellement développé. Les conséquences de cette érosion sur les propriétés mécaniques des sols ont été caractérisées par la réalisation en laboratoire et la simulation d'essais de compression triaxiale. Le développement de la suffusion et ses conséquences sur le comportement macroscopique sont discutés en terme de densité et granulométrie initiale du sol, de chargement hydraulique et du type de particules érodées (actives ou inactive vis-à-vis du transfert de force). / Suffusion is a particular case of internal erosion taking place in hydraulic earth structures. It is characterized by the detachment and migration of fine particles by interstitial flow leaving behind the granular skeleton. Such modifications in the soil microstructure may lead to deformations at the macroscopic scale and may influence the mechanical behavior of the soil. This research was devoted to investigate the suffusion mechanism and its impact on the mechanical properties of cohesionless soils. To achieve this objective, two approaches were followed in this work: numerical and experimental approaches. A discrete numerical model was defined to describe quantitatively the soil macroscopic behavior and to analyze the microstructure of granular assemblies. It is based on the use of spherical particles and contact rolling resistance, as well as a new method of compaction to mimic the one used in laboratory and to reach a wide range of initial densities. The model was validated through comparison of numerical results with experimental data. This model was first applied on granular assemblies with different fines contents to study the role of fine particles with respect to a given soil microstructure, without taking into account an erosion process. It was shown that there exists a fines content, below the threshold, where fine particles may start to participate in the force transfer that if suffusion initiates and mobilizes these particles, it may affect the macroscopic behavior of the soil. Thereafter, a simplified kinetics of grain extraction was proposed to describe the suffusion process. It was based on a one-way fluid-solid coupling approach. The importance of this extraction procedure is that, on one hand, it takes into account complex geometric and hydraulic criteria; on the other hand, the model involves an affordable computational cost. This procedure was validated based on experimental tests carried out using a newly developed suffusion apparatus. Consequences of this erosion on the mechanical properties of soils were characterized from experimental and simulated triaxial compression tests. Suffusion development and its consequences on the macroscopic behavior are discussed from the results obtained in terms of the initial soil density and gradation, the hydraulic loading and the type of eroded particles (active or inactive in the force transfer).

A l'échelle de plusieurs dizaines de millions d'années, un système de subduction implique de grandes déformations de la plaque plongeante assimilée un solide viscoélastique, et du manteau supérieur assimilé à un fluide newtonien. L'objectif de ce travail est de développer une stratégie de couplage solide-fluide appliquée à l'étude de l'interaction lithosphère-asthénosphère. Cette stratégie est basée sur l'utilisation de maillages non-conformes aux interfaces et d'une méthode de domaines fictifs (MDF) pour la résolution du problème fluide. Pour l'efficience des modèles 3D, nous employons une formulation simplifiée de la méthode de domaines fictifs par multiplicateurs de Lagrange. La MDF développée est validée par des comparaisons avec des solutions analytiques qui montrent que la méthode est d'ordre 1. La stratégie de couplage est également validée par la comparaison avec d'autres méthodes de couplage solide-fluide. Une première étude est ensuite menée pour analyser l'influence de certains paramètres rhéologiques et cinématiques sur la dynamique d'une subduction contrôlée par les vitesses des plaques. Cette étude, en 2D, concerne plus spécifiquement le mécanisme de plissement périodique du slab lorsque celui-ci est ancré à 660 km de profondeur. Ce mécanisme induit des variations de pendage du slab générant des variations de l'état de contrainte de la plaque chevauchante. Un intérêt particulier est porté sur l'influence de la viscosité du manteau sur les plissements. Dans ce cadre, nous réalisons une application à la subduction andine. / Over the time scale of tens of millions of years, a subduction system involves large deformations of tectonics plates, as one plate sinks into the Earth's mantle. The aim of this work was to develop a soli-fluid coupling method applied to the lithosphere-asthenosphere interaction in the context of subduction zones. Plates were assumed to behave as viscoelastic bodies, while the upper mantle was assimilated to a newtonian fluid. The method developped here is based on the use of non-matching interface meshes and a fictitious domain method (FDM) for the fluid problem. To optimize the computational efficiency of 3D model, we used a simplified version of the Lagrange multipliers fictitious domain method. The developped FDM has been benchmarked with analytical solutions and we showed that this FDM is a first-order method. The coupling method has also been compared to other fluid-solid coupling methods using matching interfaces meshes. A first two-dimensional study was performed in order to evaluate the influence of some rheological and kinematic parameters on the dynamics of a subduction controlled by the velocity of the plates. This study aimed at investigating cyclic slab folding over a rigid 660 km depth transition zone. This folding mechanism induces variations in slab dip that generate variations in the stress state of the overriding plate. We focussed on the influence of the upper mantle viscosity on slab folding. We also applied this model to the Andean subduction zone. Several studies have determined a cyclic variation of the South-American tectonic regime (period of 30-40~Myrs) which may have been related to the slab dip evolution.

Cette thèse est consacrée à trois différentes équations d’évolution non-linéaires dans le cadre de mécanique des fluides : le système fluide-solide, le système de Boussinesq et un modèle de films liquides. Pour le système fluide-solide, nous étudions l’évolution d’un petit solide en mouvement dans un fluide newtonien incompressible dans le cas où l’obstacle se contracte vers un point. En supposant que la densité du solide tend vers l’infini, nous montrons la convergence des solutions du système fluide-solide vers une solution des équations de Navier-Stokes dans $\mathbb{R}^d$ , avec $d^2$ et 3. Pour le problème de convection, nous travaillons sur l’unicité des solutions ‘mild’ du système de Boussinesq et généralise de plusieurs manières différentes des résultats classiques d’unicité pour les équations de Navier-Stokes. Dans la dernière partie, nous exposons nos contributions à l’étude des interface 2D de films liquides en dimension trois. Nous montrons qu’une variante 2D, non-local, de l’équation de Kuramoto-Sivashinsky admet un attracteur globale compact et obtenons enfin une majoration du nombre d’oscillations spatiales des solutions / This thesis is devoted to three different non-linear evolution equations in fluid mechanics : the fluid-solid system, the Boussinesq system and a falling films model. For the fluid-solid system, we study the evolution of a small moving solid in incompressible viscous fluid in the case the obstacle converges to a point. Assuming that the density of the solid tends to infinity, we prove that the rigid body has no influence on the limit equation by showing the convergence of solutions of the fluid-solid system towards to a solution of the Navier-Stokes equations in the full $\mathbb{R}^d$ , avec $d^2$ et 3. For the convection problem, we provide several uniqueness classes on the velocity and the temperature and generalize some classical uniqueness result for ‘mild’ solutions of the Navier-Stokes equations. We then work on a falling films model in three dimensions (2D interface). We show that a non-local variant of the Kuramoto-Sivashinsky equation admits a compact global attractor and we study the number of spatial oscillations of the solutions

Le présent travail porte sur le problème géologique de formation de roches à configuration complexe, appelée skarns. Un des intérêts de ces roches provient de ce qu'elles contiennent des composants ayant une valeur économique certaine. Après avoir mis en place le langage à utiliser et donné une idée des échelles des phénomènes, nous fournirons un modèle qui espère représenter l'évolution des skarns en prenant en considération plusieurs composants chimiques, plusieurs dimensions d'espace et des conditions limites variables en "lieu" et "temps". On aboutit à un système d'équations aux dérivées partielles relativement complet dont on diminue la taille en évoquant des hypothèses géologiques simplificatrices. Dans la mesure où l'on attendra des solutions faibles (avec "chocs") de ces systèmes, on essaiera de préciser la forme exacte des équations à traiter (liée aux quantités "effectivement conservées"), ce qui n'a pas toujours été évident par le passé. Afin de nous donner les moyens d'un traitement numérique de notre problème, nous fournirons une méthode de transport-écroulement en dimension d'''espace" quelconque finie, un prolongement heuristique de cette méthode à des cas plus complexes, y compris aux multiples problèmes de conditions limites rencontrés. Seront aussi présents divers paragraphes techniques utilisés systématiquement dans les autres parties. On fera enfin une programmation des diverses méthodes avancées avec, entre autres, des sorties graphiques couleur. Nous voudrions préciser avant de nous engager dans cette étude, que le travail que nous avons pu effectuer ici n'était pas évident à mener à terme : Il s'agit en effet là d'un travail de pionnier, avec au départ, du point de vue mathématique, très peu de formulations faites.

Application a un cristal liquide nématique dans lequel la déformation du directeur est induite soit par un écoulement de poiseuille plan du cristal liquide, soit par application d'un champ électrique. Détermination de l'énergie d'ancrage du nématique en contact avec son support.

Nous nous intéressons à la modélisation et à la simulation numérique d'écoulements gravitaires transitoires à surface libre, pour des fluides visqueux et incompressibles. La loi de comportement est de type viscoplastique avec transition fluide/solide. Plus précisément, nous considérons la loi rhéologique de Drucker--Prager. Nous nous intéressons tout d'abord le cas unidimensionnel d'un écoulement longitudinal cisaillé. Nous étudions un modèle simplifié, avec terme source empirique, pour lequel nous concevons une méthode numérique pour le suivi de la position de l'interface entre la couche solide et la couche fluide. Nous présentons des résultats numériques, avec divers termes sources, et nous comparons ces résultats, lorsque la viscosité est petite, à la solution analytique non visqueuse. Dans le cas visqueux, nous étudions les phases de démarrage et d'arrêt de l'écoulement. Dans un second temps, nous étudions le cas bidimensionnel d'écoulement de Drucker--Prager avec surface libre. La loi de comportement du fluide est traitée par régularisation, et nous utilisons la méthode ALE pour traiter le mouvement du domaine. Nous présentons des résultats numériques pour l'étude de la mise en mouvement d'un talus.

La détermination de la forme d'un obstacle élastique immergé dans un milieu fluide à partir de mesures du champ d'onde diffracté est un problème d'un vif intérêt dans de nombreux domaines tels que le sonar, l'exploration géophysique et l'imagerie médicale. A cause de son caractère non-linéaire et mal posé, ce problème inverse de l'obstacle (IOP) est très difficile à résoudre, particulièrement d'un point de vue numérique. De plus, son étude requiert la compréhension de la théorie du problème de diffraction direct (DP) associé, et la maîtrise des méthodes de résolution correspondantes. Le travail accompli ici se rapporte à l'analyse mathématique et numérique du DP élasto-acoustique et de l'IOP. En particulier, nous avons développé un code de simulation numérique performant pour la propagation des ondes associée à ce type de milieux, basé sur une méthode de type DG qui emploie des éléments finis d'ordre supérieur et des éléments courbes à l'interface afin de mieux représenter l'interaction fluide-structure, et nous l'appliquons à la reconstruction d'objets par la mise en oeuvre d'une méthode de Newton régularisée.