Les fluctuations thermiques présentes dans les tés de mélange provoquent des contraintes thermiques qui peuvent mener à l'apparition de fissures qui se propagent plus ou moins rapidement dans la structure. Une possibilité pour réduire ces risques est d'installer des mélangeurs statiques (statics mixers) pour accroître le mélange. Une telle technologie a été utilisée par Utveckling AB depuis 1980 dans des installations nucléaires. Toutefois, ces technologies sont très coûteuses. C'est pour cette raison que plusieurs investigations numériques ont été faites pour prédire les fluctuations de température causées par le mélange turbulent dans cette configuration d'écoulement. On effectue la simulation numérique de l'écoulement sur deux types de té de mélange. L'un avec des bords droits et une paroi en Plexiglas, et l'autre avec des bords arrondis et une paroi en inox 304L. Dans le premier cas la condition de paroi est adiabatique et dans le second cas on effectue un couplage entre le code CFD (Computational Fluid Dynamic) Code_Saturne et le code SYRTHES pour l'analyse de la température dans le solide. L'apport principal de la thèse est la détermination des statistiques temporelles d'ordre élevé dans une configuration aussi complexe. En effet, les équations de transport de l'énergie cinétique turbulente, de la variance de la température et des flux thermiques turbulents sont déterminées dans les deux configurations (adiabatique et avec des parois en inox), ce qui montre l'influence de la paroi sur le transfert de chaleur en région proche paroi. / Thermal fluctuations present in mixing tees cause thermal stresses that can lead to the appearation of cracks, which spread more or less quickly in the structure. One possibility to reduce these risks is to set static mixers (statics mixers) to increase the mixture. Such technology has been used by Utveckling AB since 1980 in nuclear installations. However, these technologies are very expensive. It is for this reason that many numerical investigations have been made to predict temperature fluctuations caused by turbulent mixing in this configuration flow. The resolution of the conservation equations is made with a finite volume approach using large eddy simulation or LES . The subgrid models used are Smagorinsky, WALE (Wall Adapted Local Eddy) and dynamic Smagorinsky. The SGDH model (Simple Gradient Di? Usion Hypothesis) is used for modeling greeting thermal subgrid and the turbulent Prandtl number is fixed one. Generation turbulence input field is made using the SEM method (Synthetic Eddy Method). The main contribution of this thesis is the determination of time turbulent statistic in a complex configuration. Indeed, the transport equations of turbulent kinetic energy, temperature variance and turbulent heat flux are determined in both configurations (adiabatic walls and stainless steel), which shows the influence of the wall on heat transfer in near-wall region.

Autour de nous, les systèmes granulaires sont omniprésents et rarement dissociés d'une phase liquide ou gazeuse. L'interaction entre les deux phases (solide et fluide) implique une dynamique d'écoulement complexe de ces systèmes couplés. Cette thèse expérimentale a pour but de caractériser la dynamique d'effondrement gravitaire instationnaire d'une colonne fluide-grains dans différentes configurations de mélange diphasique. Considérant uniquement l'effet des grains au sein du fluide, la dynamique d'écoulement d'une suspension isodense dans l'air est d'abord abordée avec une description macroscopique. En particulier, La rhéologie du fluide apparent est extraite à partir de l'évolution temporelle du front en utilisant des solutions auto-similaires, modèles de la dynamique de propagation aux temps longs. Sur la gamme du nombre de Reynolds étudiée, la dynamique d'écoulement de la suspension est décrite comme celle d'un fluide apparent Newtonien ou non-Newtonien (rhéofluidifiant, rhéoépaississant, viscoplastique) dépendant des paramètres considérés (fraction volumique, viscosité du fluide interstitiel, diamètre des grains, protocole de préparation). Afin de décrire la configuration opposée d'un massif sédimentaire pour lequel les interactions entre particules solides deviennent dominantes, la configuration d'une colonne granulaire dense saturée en fluide, i.e. où les grains sont plus lourds que le fluide porteur, est ensuite étudiée. Un travail préliminaire est consacré à la caractérisation de l'effondrement d'une colonne granulaire totalement immergée dans un fluide. Le rôle dissipatif de ce dernier sur le milieu granulaire est mis en évidence par une analyse de la dynamique d'effondrement et des caractéristiques du dépôt final. Cette caractérisation permet de classifier les régimes d'effondrement en fonction de la viscosité et de la masse volumique du fluide environnant, i.e. du nombre de Stokes et du rapport de masse volumique entre le fluide et les grains. Dans le cas triphasique, i.e. lorsque la colonne fluide-grains s'effondre dans l'air, la dynamique peut être fortement affectée par les effets capillaires à travers le nombre de Bond qui contrôle alors la mise en mouvement de la colonne initiale. Quand ces effets deviennent négligeables à l'échelle de la colonne et à celle du grain, le fluide interstitiel peut jouer un rôle moteur ou dissipatif vis-à-vis dumilieu granulaire conduisant à une longueur d'étalement plus ou moins importante en comparaison au cas sec. Le rôle du fluide interstitiel dépend essentiellement de sa viscosité modifiant, par la même occasion, la dynamique d'effondrement. Enfin, une étude préliminaire sur la dynamique d'écoulement d'une suspension non-isodense, initialement homogène, est réalisée. Cette configuration, à la transition des situations décrites précédemment, permet d'aborder le couplage de la dynamique de sédimentation des grains et celle du courant. En particulier, nous observons que la vitesse de sédimentation décroît avec l'augmentation de la fraction volumiqueinitiale en particules.

L'avance technologique a encouragé, au cours de la dernière décennie, la fabrication et l'utilisation de nouveaux matériaux, ayant des propriétés innovantes et un comportement complexe. Leur développement a été conditionné par la nécessité d'augmenter l'efficacité de certains procédés technologiques, dans des domaines importants tels que la médecine, la fabrication de médicaments, de produits cosmétiques, de produits alimentaires, la production d'énergie, le recyclage des matériaux usagés, par exemple par l'industrie des transports dont la construction automobile. Dans le cadre de cette étude, nous nous sommes intéressés au comportement de fluides issus de ces matériaux complexes. En effet, leur développement, leur production et leur utilisation sont conditionnés par la compréhension et la modélisation de leur rhéologie, tant au niveau macroscopique que microscopique, voire même nanométrique. La rhéologie d'un fluide est implicitement déterminée par ses propriétés de viscosité, d'élasticité et de plasticité. Un grand intérêt a toujours été porté au comportement viscoélastique des matériaux, la viscoélasticité d'une matière étant la propriété à manifester simultanément un comportement rhéologique visqueux et élastique, quand elle est soumise à une déformation, en raison de l'application d'une tension. La caractérisation rhéologique d'un liquide suppose la détermination du tenseur de tension supplémentaire, présent dans l'équation de mouvement, provenant de mesures rhéométriques. L'expression de ce tenseur, qui dépend de la déformation et de la vitesse de déformation, s'appelle la relation ou loi constitutive. La détermination du comportement rhéologique des fluides, / This paper was dedicated exclusively to the rheological characterization of simple and complex fluids and a detailed analyze of available rheometric testing procedures. Throughout this thesis the rheological behavior of both simple and complex fluids has been studied and modeled in both shear and the complex motions (shear tests and squeezing tests). By linking conventional experimental test methods with numerical simulations of real flows, this thesis introduces a new concept in rheology: Computational Rheometry. All experimental investigations carried out in this study for the squeezing flows were accompanied by numerical simulations. For the oscillatory squeezing flow the influence of initial film thikness, oscillatory amplitude and frequency, computational time step was investigated by comparison with the theoretical predictions of squeeze force and a Genrealized Reynolds Equation inclued in a finite element code in Fortran. A validity domain was established for the analitical formulation of squeezing force. The constant velocity squeeze flow was investigated also using a quasi-steady approximation of the motion, which brings a significant reduction of the computational time, and a very good correlation with the transient (deformable mesh) approximation and the analytical predictions. The investigation of free surface influence on the distribution of normal force in both constant velocity and oscillatory squeeze flow was analyzed. In the case of constant velocity squeeze flow, the numerical simulations coupled with the free surface evolution and measured normal force during experimental investigations are suggesting the presence of a partial slip during the experimen

Les écoulements autour de structures en mouvement font l'objet de plusieurs travaux numériques et expérimentaux. L'objectif de ce travail de thèse consiste à montrer la pertinence de la combinaison de la pénalisation volumique avec la méthode de lattice Boltzmann (LBM), dans l'étude du mouvement d'obstacles mobiles dans un écoulement, et de leur interaction avec celui-ci. La LBM,qui est simple et précise à mettre en œuvre, a prouvé ces dernières années son efficacité dans le domaine de la mécanique des fluides. Par ailleurs, la méthode de pénalisation volumique consiste à introduire un terme de pénalisation dans l'équation que l'on souhaite résoudre, afin de prendre en compte l'influence de l'obstacle sur le domaine fluide. Comme cette équation est résolue sur l'ensemble du domaine composé du fluide et du solide, les conditions aux limites à l'interface fluide-solide sont appliquées de façon naturelle. Il semble donc aisé de combiner cette technique avec la méthode de lattice Boltzmann. Nous avons dans un premier temps rappelé les notions de base et les principales caractéristiques de la méthode de lattice Boltzmann. On a présenté quelques exemples d'applications sur des cas tests, que nous avons programmés. Ensuite, une étude bibliographique faisant état des différentes approches qui utilisent la LBM dans l'étude des problèmes d'interaction fluide structure (IFS) a été réalisée. Puis, la combinaison de la pénalisation volumique avec la LBM a été testée avec succès sur l'équation de Burgers monodimensionnelle. La validation s'est portée en premier lieu, sur un écoulement autour d'un solide fixe, puis sur un écoulement autour d'une structure dont le mouvement est imposé, et finalement sur un problème d'IFS de type masse-ressort. La méthode développée a été ensuite testée sur les équations de Navier-Stokes, en considérant un fluide incompressible et une structure rigide. La validation s'est portée tout d'abord sur un écoulement autour d'obstacles immobiles (carré et cylindre), puis autour d'un cylindre mobile en oscillations forcées et libres. Enfin, une dernière application a été portée sur un écoulement entre deux plaques mobiles dans un canal. Nous avons montré que pour tous les cas étudiés, l'approche développée donne de bons résultats, elle reproduit avec précision les résultats de référence. / Flows around moving bodies are the subject of several numerical and experimental studies. The work presented in this document deals with the implementation of a volume penalization technique in a lattice Boltzmann model (LBM), in order to compute flows around moving obstacles. The LBM, which is accurate and easy to implement, has been successfully applied in fluid mechanics during the last decades. It was thus chosen in the present work, for flow computation. Furthermore, the volume penalization technique consists in introducing a volume penalization term into the equation that needs to be solved, in order to take into account the influence of the obstacle on the fluid domain. Since this equation is solved on both fluid and solid domains, the boundary conditions at the fluid-solid interface are naturally applied. Hence this technique seems easy to implement in a lattice Boltzmann framework. In the first chapter, the foundations and the main features of the lattice Boltzmann method are recalled, and several test cases that we simulated are presented. The second chapter deals with a literature review of the techniques developed for the simulation of fluid structure interaction problems in combination with the LBM. In the third chapter, the volume penalization method combined with the LBM was first applied to the one dimensional Burgers equation, considering motionless and moving obstacles (forced motion, and coupling between the fluid force calculated with the penalized Burgers equation and the motion of the obstacle). The combination of the volume penalization approach and the LBM was then employed to solve the incompressible NavierStokes equations, for cases of flows past motionless obstacles (flows over a square obstacle, and past a circular cylinder), and past an oscillating cylinder (where forced and free oscillations of the cylinder were simulated). Finally, this method was also applied to a symmetric Couette flow. For all these simulated cases, a good agreement with numerical results obtained with other techniques, and with results found in literature, was obtained.

Nous étudions de manière expérimentale un cas d’interaction fluide/structure entre une fibre élastique et un milieu granulaire dense proche du jamming. Dans notre dispositif expérimental horinzontal, un milieu granulaire bidimensionnel, aléatoire et dense est en interaction avec une fibre élastique encastrée. L'extrémité encastrée de cette fibre est en translation par rapport aux grains tandis que son extrémité libre fait initialement face au flux de grains dans une géométrie de pénétration, avant de se défléchir par interaction avec les grains. Les paramètres de contrôle de notre système sont la compacité du milieu granulaire, la longueur et la rigidité de la fibre. Notre dispositif nous permet de mesurer la déformée de la fibre, les efforts que les grains exercent sur la fibre, notamment la force de traînée et de suivre les réorganisation du milieu granulaire. Sous certaines conditions expérimentales, nous avons observé que ce système présente plusieurs régimes d'évolution. A faible et à forte déflexions de la fibre, le système connaît une évolution stochastique, liée à la nature aléatoire du milieu granulaire : respectivement de petites fluctuations de la fibre autour de sa position d'équilibre et des événements brusques de balayage de la fibre en flexion associés à des avalanches granulaires. La transition entre ces deux régimes est associée à une brisure de symétrie qui affecte tant la fibre que le milieu granulaire. Dans une première partie, nous étudions comment la déflexion de la fibre est influencée par les paramètres de contrôle, en particulier par l'écart au jamming. Dans une seconde partie, nous montrons qu'une asymétrie de densité se développe dans le milieu granulaire lorsque la fibre est suffisamment défléchie. Dans une troisième partie, nous tentons de comprendre l'évolution des forces mesurées à l'aide de simulations numériques de l'Elastica sur une poutre sous chargement. Enfin, nous proposons un modèle de compaction des grains par la fibre qui est en bon accord avec l'évolution expérimentale du système. Une application de cette thèse est la compréhension des mécanismes de pénétration d'une racine ou d'une tige dans un sol granulaire. / We study experimentally a fluid/structure interaction between an elastic fiber and a dense granular medium near the jamming transition. In our horizontal setup, a bidimensionnal and dense random packing is interacting with a clamped elastic fiber. The clamped end of the fiber is in translation relative to the the grains while the free end faces initially the grain flux in a penetration geometry before being deflected by grain/fiber interactions. The control parameters of our system are the packing fraction of the medium, the length and the rigidity of the fiber. Our setup allows us to measure the deflected shape of the fiber, the forces exerted by the grains on the fiber (including the drag force) and to follow the reorganizations of the granular medium. Under some experimental conditions, we observed that this system exhibits several regimes of evolution. For low and large fiber deflections, the system evolves stochastically as a consequence of the random nature of the granular medium : respectively small fluctuations of the fiber around its equilibrium position and sudden events of unlozding of the fiber in flexion associated to granular avalanches. The transition between those two regimes is associated with a symmetry breaking that affects the granular medium as well as the fiber. In a first part, we study how the fiber deflection is influenced by the control parameters, particularly by the distance to jamming. In a second part, we show that a packing fraction asymmetry appears in the granular medium when the fiber is deflected enough. In a third part, we try to understand the evolution of the measured forces using numerical simulations of the Elastica on a constrained beam. Finally, we propose a model of compaction of the grains by the fiber that is in a good aggreement with the experimental evolution of the system. An application of this thesis is the understanding of the mechanisms involved in the penetration of a root or a beam in a granular soil.

Les écoulements pyroclastiques sont des mélanges à haute température de gaz et de particules volcaniques qui peuvent se propager sur de très grandes distances. Cette forte « mobilité » est souvent attribuée à leur capacité à se fluidiser, c’est-à-dire à générer et conserver une forte pression interstitielle de gaz qui réduit les forces de friction interne. L’objectif principal de cette thèse est de comprendre comment les irrégularités des terrains sur lesquelles se propagent les écoulements pyroclastiques peuvent favoriser leur fluidisation. Une première série d’expériences de laboratoire a consisté à générer des écoulements de particules fines (diamètre de 45-90 μm) sur des substrats de différentes rugosités. Les résultats montrent que la distance de parcours des écoulements augmente avec la rugosité, allant jusqu’à doubler par rapport à la distance de parcours sur fond lisse. Des analyses de vidéos haute vitesse et des mesures de pression interstitielle d’air à la base des écoulements montrent que la tête (partie antérieure) des écoulements qui se propagent sur un substrat rugueux s’auto-fluidisent en conséquence de la sédimentation des particules dans les interstices du substrat, chassant l’air qui remonte et percole dans l’écoulement. Ce mécanisme d’auto-fluidisation est efficace pour toutes les inclinaisons étudiées (0-30°), suggérant qu’il est susceptible de se produire tout au long de la mise en place d’un écoulement pyroclastique. Une seconde étude a consisté à faire chuter des lits de particules dans une colonne statique. Les résultats montrent que même pour une hauteur de relâchement relativement faible (20 cm), le mélange peut entièrement s’auto-fluidiser durant sa chute. Quand les particules sont suffisamment fines (<100 μm) la pression interstitielle dans le dépôt diffuse pendant plusieurs secondes, la durée de cette diffusion augmentant avec l’augmentation de l’épaisseur du lit et la diminution de taille des particules. Les temps de diffusions les plus longs sont observés avec un matériau provenant d’un dépôt d’écoulement pyroclastique (~30 s pour des lits de 28.5 cm d’épaisseur). Ces résultats suggèrent que les écoulements pyroclastiques qui se propagent sur des terrains accidentés peuvent s’auto-fluidiser et conserver une faible friction au cours de leur mise en place. / Pyroclastic flows are hot mixtures of gas and particles that can propagate over large distances. This high “mobility” is often attributed to their ability to be fluidized, that is, to generate and retain high gas pore pressure that reduces internal friction forces. The main objective of this thesis is to understand how irregularities of substrates on which pyroclastic flows propagate can enhance their fluidization. A first set of laboratory experiments consisted of the generation of fine-grained flows (diameter of 45-90 μm) on substrate of various roughness. Results show that the flow runout distance increases with the substrate roughness, and is up to twice the runout on a smooth substrate. High speed video analyses and air pore pressure measurements at the flow base show that the flow head propagating over a rough substrate can auto-fluidize because of particles sedimentation into the substrate interstices, which forces the air to escape upward and percolate through the flow. This auto-fluidization mechanism is efficient at all inclinations investigated (0-30°), suggesting that it could occur during the whole emplacement of a pyroclastic flow. A second study consisted of the vertical release of beds of particles in a static column. Results show that the granular mixture can be fully fluidized, even when collapsing from a relatively low height (20 cm). When particles are fine enough (<100 μm), pore pressure in the deposit diffuses for several seconds, the diffusion duration increasing with increasing bed thickness and decreasing particle size. The longest diffusion durations are observed for pyroclastic flow deposit materials (~30 s for 28.5 cm thick beds). These results suggest that pyroclastic flows propagating on irregular terrains can auto-fluidize and preserve low internal friction during their emplacement.

Un modèle mathématique d'interfaces diffuses pour l'interaction de N solides élasto-plastiques a été construit. C'est une extension du modèle développé par Favrie & Gavrilyuk (2012) pour l'interaction d'un fluide et d'un solide. En dépit du grand nombre d'équations présentes dans ce modèle, deux propriétés remarquables ont été démontrées : ce modèle est hyperbolique (quelles que soient les déformations admissibles) et il vérifie le second principe de la thermodynamique. En dépit du grand nombre d'équations présentes dans ce modèle, deux propriétés remarquables ont été démontrées: ce modèle est hyperbolique (quelles que soient les déformations admissibles) et il vérifie le second principe de la thermodynamique. L'énergie interne de chaque solide est prise sous forme séparable: c'est la somme d'une énergie hydrodynamique qui ne dépend que de la densité et de l'entropie, et d'une énergie de cisaillement. L'équation d'état de chaque solide est telle que si nous prenons le module de cisaillement du solide égale à zéro, on retrouve les équations de la mécanique des fluides. Ce modèle permet, en particulier, de:- prédire les déformations de solides élasto-plastiques en petites déformations et en très grandes déformations.- prédire l'interaction d'un nombre arbitraire de solides élasto-plastiqueset de fluides. L'aptitude de ce modèle à résoudre des problèmes complexes a été démontrée. Sans être exhaustif, on peut citer:-le phénomène d'écaillage dans les solides.- La fracturation et la fragmentation dynamique dans les solides. / A mathematical model of diffuse interface for the interaction of N elasto-plastic solidS was built. It is an extension of the model developed by Favrie & Gavrilyuk (2012) for a fluid-solid interaction. Despite the large number of equations present in this model, two remarkable properties have been demonstrated: it is hyperbolic for any admissible deformations and satisfies the second principle of thermodynamics. In this model, the internal energy of each solid is taken in separable form: it is the sum of a hydrodynamic energy (which depends only on the density and entropy) and shear energy. The equation of state of each solid is such that if we take the shear modulus of the solid vanishes, we find the equations of fluid mechanics. This model allows, in particular:- predict the deformation of elastic-plastic solids in small and very large deformations.- predict the interaction of an arbitrary number of elasto-plastic solids and fluids.The ability of this model to solve complex problems has been demonstrated. Without being exhaustive, one can mention:- the spall phenomenon in solids.- fracturing and fragmentation in solids.

Ce travail vise à obtenir un matériau possédant des propriétés de complexation grâce aux cyclodextrines (CyD) et à démontrer la faisabilité de la synthèse en milieu CO2 supercritique (scCO2). L'adaptation en scCO2 de la polyaddition directe entre une α-CyD et un diisocyanate n’a montré qu’une faible conversion après 7 heures de réaction (presque totale en 4 heure en milieu liquide), en raison de la complexation possible du diisocyanate par la CyD, de l’insolubilité de la CyD dans le scCO2 ou d’une diminution de la vitesse de réaction nécessitant un temps supérieure pour atteindre une conversion plus élevée. Des études supplémentaires sont nécessaires pour confirmer ses hypothèses. La conception d’un monomère original pouvant être scCO2-soluble, par acétylation des hydroxyles de la CyD, fonctionnalisé par 3 isocyanates a montré, par polyaddition avec une alkyldiamine, en milieu solvant, une forte dépendance de la masse molaire avec la concentration (jusqu’à 16,89 kg/mol à 50 %pds en monomères). En milieu scCO2, la conversion du monomère original a été totale, avec des masses molaires atteintes supérieures à celles en milieu liquide à concentration identique, démontrant l’intérêt d’un tel milieu en remplacement de solvants organiques dans le cas d’une conception bien adaptée. L’extension de la gamme de concentrations et des conditions opératoires (température, pression) reste à réaliser pour conclure à l’intérêt du scCO2 par rapport à ce monomère original à base de CyD / This work aims at obtaining a polymeric material possessing complexation abilities thanks to cyclodextrins (CyD) within its framework and demonstrating the feasibility of the synthesis in supercritical CO2 medium (scCO2). Contrary to liquid medium, the polymerization of α-CyD and a diisocyanate in scCO2 shows only low conversion after 7 hours due to the possible complexation of the diisocyanate by the CyD, the insolubility of the CyD in scCO2 medium or insufficient time to achieve high conversion due to diminution of the reaction rate. New studies have to be done in order to confirm these hypothesis. The design of an original monomer which could be soluble in scCO2 thanks to acetyl groups fixed on the CyD and functionalized by 3 isocyanates, leads to a strong dependency of the molar mass with the concentration of monomers (up to 16.89 kg/mol at 50 %w of monomers) when reacted with a diaminoalkyle chain extender in solvent medium. In scCO2 medium, full conversion of this monomer is achieved, with higher molar masses obtained than in liquid medium at the same concentration. This demonstrates the potential of such a synthesis medium instead of organic solvents when design is well adapted. However, more experiments are needed to study the influence of concentration, temperature and pressure on larger scales to conclude in the interest of scCO2

La physique de nage d’une nageoire caudale flexible est étudiée expérimentalement grâce à une plaque flexible immergée dans un écoulement uniforme. Le bord d’attaque est forcé par un mouvement harmonique, tandis que le bord de fuite répond passivement au forçage. Une résonance en amplitude dans la réponse est mise en évidence et révèle une forte interaction entre les fréquences naturelles de la plaque et celles du forçage. Les résultats expérimentaux suggèrent un rôle non trivial de l’amplitude de forçage, qui souligne le rôle des non linéarités dans ce problème. Pour avoir une meilleure compréhension de l'origine de ces non linéarités, un modèle analytique faiblement non linéaire a été développé. Nous supposons une plaque d'épaisseur nulle immergée dans un écoulement potentiel, sujette à une force de traînée résistive. La déflection de la plaque a ensuite été décomposée en un mode rigide, mimant le forçage en pilonnement et en modes propres de flexion d’une poutre dans le vide. Les prédictions de la réponse en amplitude et en fréquence du système forcé sont alors calculées. Les fréquences de résonances, ainsi que l’enveloppe de la déflection, sont en bon accord avec les résultats expérimentaux. Les performances du système, mesurées à travers la poussée générée, est également correctement prédite par la modélisation. L’accord entre les expériences et le modèle est étendu à une étude trouvée dans la littérature. Une optimisation analytique a été conduite et étendue à l’application de la bio-robotique. / The physics of the swimming with a flexible caudal fin is studied experimentally by the means of an elastic plate immersed in a uniform water flow. The leading edge of the plate is forced into a harmonic motion, while its trailing edge responds passively to this actuation. A resonance response in amplitude is demonstrated, revealing a strong interaction between the natural frequencies of the plate and the forcing frequencies. Experimental results suggest a non-trivial role of the forcing amplitude, emphasizing the role of non linearities in this problem. To gain better insight into the origin of these non linearities, a weakly non linear model is developed. We model a quasi two-dimensional plate of zero thickness immersed in a potential flow and subject to a resistive drag-like force. The plate deflection is then decomposed into a forcing heaving mode and natural flexural modes. Predictions of the response in amplitude and frequency to a heave forcing system are then calculated. The frequencies of the resonances, as well as the shapes of the deflection, match the experimental results. The performance of the system measured through the generated thrust, is well predicted by the modelling. The experimental and modelling results presented here show (very) good agreement with the literature. Finally, an analytical optimization is undertaken and potential applications to bio-robotic are suggested.

L'objectif de ce travail est de prédire le gonflement de fûts d'enrobés bitumineux contenant des déchets radioactifs en condition d'entreposage, de transport et de stockage. L'incorporation de radionucléides dans du bitume provoque en effet sa radiolyse et du dihydrogène est ainsi généré dans le matériau. Ce gaz est solubilisé dans la matrice jusqu'à ce que sa concentration atteigne une valeur limite, la saturation. Au-delà de celle-ci, un mécanisme de germination a lieu : des bulles sont formées. Elles vont par la suite croître grâce à l'apport de gaz par radiolyse et migrer sous l'effet de la poussée d'Archimède. Le gonflement du matériau résulte alors de la compétition entre création et évacuation de gaz. Un modèle a été construit pour décrire l'évolution d'une population de bulles. Sa résolution analytique n'étant pas possible, un code numérique a été développé. Les résultats montrent qu'une augmentation de la viscosité du fluide, du terme source de gaz apporté par radiolyse ou de la hauteur du récipient induisent une augmentation du gonflement et qu'une augmentation du coefficient de diffusion du gaz dans la matrice diminue le gonflement. Lorsque le fluide présente un seuil d'écoulement, le comportement de la population de bulles est modifié et l'évacuation de gaz par migration des bulles se fait par à-coups, contrairement au cas d'un fluide newtonien où une évacuation stationnaire est atteinte. Cette dernière a enfin servi à l'optimisation des conditions opératoires. / The aim of this PhD thesis is to predict the swelling of bitumen barrels in which radioactive salts are mixed. The bitumen exposed to radioactivity undergoes a chemical reaction: the radiolys. This implies a generation of dihydrogen. The created is solubilized until the concentration reaches a limit value which is called saturation. Over this limit nucleation of bubbles is observed. Then they will grow thank to the contribution of the gas generated by radiolys and they will be submitted to Archimede's principle so that they will rise in the fluid. The swelling is the result of the competition between generation and evacuation of gas. A model has been built to describe theevolution of a bubble population. Because of it is not possible to solve it analitically, a numerical program was developed. The results show that an increase of the fluid viscosity, the gas generation or the container height lead to an increase of the swelling and that an increase of the diffusioncoefficient contributes to a decrease of the swelling. In the particular case of a yield stress fluid, the behavior of the bubble population is modified and the evacuation of gas is done with several shots, at the opposite of the case of a Newtonian fluid for which a stationary evacuation is reached.