This work is focused on the numerical modelling of fluid-structure interaction in three dimensions. Both internal and external laminar flow around flexible bodies are considered. The fluid flow simulated is based on the incompressible Navier-Stokes equations and the general focus is on laminar Newtonian flow. The streamline upwind/ pressure stabilising Petrov-Galerkin (SUPG/PSPG) method is employed to achieve a stable low order finite element discretisation of the fluid, while the solid is discretised spatially by a standard Galerkin finite element approach. The behavior of the solid is governed by Neo-Hooke elasticity. For temporal discretisation the discrete implicit generalised-alpha method is employed for both the fluid and the solid domains. The motion of the fluid mesh is solved using an arbitrary Lagrangian-Eulerian (ALE) scheme employing a nonlinear pseudo-elastic mesh update method. The fluid-solid interface is modelled using a finite element interpolation method that allows for non-matching meshes and satisfies the required conservation laws. The resulting sets of fully implicit strongly coupled nonlinear equations are then decomposed into a general framework consisting of fluid, interface and solid domains. These equations are then solved using different solution techniques consisting of strongly coupled monolithic Newton and block Gauss-Seidel methods as well as a weakly coupled novel staggered scheme. These solvers are employed to solve a number of three dimensional numerical examples consisting of: External flow: o a soft elastic beam fixed at both ends o a thin cantilever plate.

Doctor of Philosophy / Department of Mechanical and Nuclear Engineering / Mingjun Wei / In this study, a global Proper Orthogonal Decomposition (POD)-Galerkin based Reduced Order model (ROM) is proposed. It is extended from usual fixed-domain problems to more general fluid-solid systems with moving boundaries/interfaces. The idea of the extension is similar to the immersed boundary method in numerical simulations which uses embedded forcing terms to represent boundary motions and domain changes. This immersed boundary method allows a globally defined fixed domain including both fluid and solid, where POD-Galerkin projection can be directly applied. However, such a modified approach cannot get away with the unsteadiness of boundary terms which appear as time-dependent coefficients in the new Galerkin model. These coefficients need to be pre-computed for prescribed periodic motion, or worse, to be computed at each time step for non-prescribed (e.g. with fluid-structure interaction) or non-periodic situations. Though computational time for each unsteady coefficient is smaller than the coefficients in a typical Galerkin model, because the associated integration is only in the close neighborhood of moving boundaries. The time cost is still much higher than a typical Galerkin model with constant coefficients. This extra expense for moving-boundary treatment eventually undermines the value of using ROMs. An aggressive approach is to decompose the moving boundary/domain to orthogonal modes and derive another low-order model with fixed coefficients for boundary motion. With this domain decomposition, an approach including two coupled low-order models both with fixed coefficients is proposed. Therefore, the new global ROM with decomposed approach is more efficient. Though the model with the domain decomposition is less accurate at the boundary, it is a fair trade-off for the benefit on saving computational cost. The study further shows, however, that the most time-consuming integration in both approaches, which come from the unsteady motion, has almost negligible impact on the overall dynamics. Dropping these time-consuming terms reduces the computation cost by at least one order while having no obvious effect on model accuracy. Based on this global POD-Galerkin based ROM with forcing term, an improved ROM which can handle the parametric variation of body motions in a certain range is also presented. This study shows that these forcing terms not only represent the moving of the boundary, but also decouple the moving parameters from the computation of model coefficients. The decoupling of control parameters provides the convenience to adapt the model for the prediction on states under variation of control parameters. An improved ROM including a shit mode seems promising in model adaptation for typical problems in a fixed domain. However, the benefit from adding a shit mode to model diminishes when the method is applied to moving-boundary problems. Instead, a combined model, which integrates data from a different set of parameters to generate the POD modes, provides a stable and accurate ROM in a certain range of parametric space for moving-boundary problems. By introducing more data from a different set of parameters, the error of the new model can be further reduced. This shows that the combined model can be trained by introducing more and more information. With the idea of the combined model, the improved global ROM with forcing terms shows impressive capability to predict problems with different unknown moving parameters, and can be used in future parametric control and optimization problems.

Submitted by Gabriel Biancolin Moimas (gabrielmoimas18@gmail.com) on 2018-06-18T11:47:19Z No. of bitstreams: 1 Dissertação.pdf: 13607292 bytes, checksum: 4574d9cd3135ad2562a9648fd1641465 (MD5) / Approved for entry into archive by Cristina Alexandra de Godoy null (cristina@adm.feis.unesp.br) on 2018-06-18T12:12:00Z (GMT) No. of bitstreams: 1 moimás_gb_me_ilha.pdf: 13607292 bytes, checksum: 4574d9cd3135ad2562a9648fd1641465 (MD5) / Made available in DSpace on 2018-06-18T12:12:00Z (GMT). No. of bitstreams: 1 moimás_gb_me_ilha.pdf: 13607292 bytes, checksum: 4574d9cd3135ad2562a9648fd1641465 (MD5) Previous issue date: 2018-05-28 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Compressores herméticos alternativos são amplamente utilizados em sistemas de refrigeração por compressão de vapor de pequeno e médio porte. O sistema de válvulas automáticas utilizadas para controlar os processos de sucção e descarga é um componente importante deste tipo de compressor. O estudo experimental ou analítico/numérico destes processos é bastante complexo devido, principalmente, à interação fluido-estrutura entre as válvulas e o escoamento de fluido refrigerante. Uma revisão de literatura mostra que este tema ainda merece muita atenção dos pesquisadores, em particular daqueles que necessitam de dados experimentais acurados para validação de metodologias numéricas. Este é o contexto deste trabalho, cujo objetivo é descrever a dinâmica do movimento de um modelo de válvula usada neste tipo de compressor. Sensores óticos foram utilizados para medir o deslocamento instantâneo de modelos de válvulas fabricadas em aço mola SAE 1070 com 0,4 e 0,5 mm de espessura. Os testes foram realizados para escoamentos com números de Reynolds variando de 1.500 a 19.000. A descrição da dinâmica da válvula foi apresentada pela sua amplitude, frequência, pressão máxima à montante do escoamento da válvula e pelo coeficiente de restituição no impacto entre a válvula e o assento. De uma maneira geral, identificamos que a dinâmica da válvula depende significativamente de sua rigidez, do valor de força de pré-carga que está sujeita e do número de Reynolds do escoamento que incide sobre sua superfície. O impacto entre válvula e assento também altera significativamente a dinâmica do movimento. / Hermetic reciprocating compressors are widely used in small and medium size vapor compression refrigeration systems. An important component of this type of compressor is the automatic valve system used for controlling the suction and discharge processes. The experimental or analytical/numerical study of this system is complex due to the fluid-structure interaction between the valves and the refrigerant flow. A literature review suggests that this theme still deserves attention from the researchers, in particular for those who need accurate experimental data to validate numerical methodologies. We aim to describe the dynamics of a suction valve model used in this type of compressor. Optical sensors were used to measure the instantaneous displacement of a spring-steel SAE 1070 valve model with 0.4 and 0.5 mm thickness. The tests were performed for Reynolds number of the flow ranging from 1,500 to 19,000. We present the dynamic behavior of the valve by its amplitude, frequency, maximum upstream pressure of the valve, and coefficient of restitution of the impact between valve and seat. Generally, we identify that the valve dynamics depends significantly on its stiffness, precharging force and Reynolds numbers of the flow. The impact between valve and seat also significantly change the dynamics of the movement.

Neste trabalho um método de partículas para simular a dinâmica de sólidos elásticos e interação fluido estrutura e implementado. O método e baseado no Moving Particle Semi-implicit (MPS), originalmente desenvolvido para escoamentos incompressíveis com superfície livre. A estratégia principal do MPS e substituir os operadores diferenciais das equações governantes por operadores diferenciais discretos em uma distribuição de nos irregulares, derivados de um modelo de interação entre partículas. Inicialmente são apresentados os detalhes da formulação do método e modelos constitutivos utilizados. Uma condição simplificada de fragmentação e proposta, assim como um algoritmo de detecção de contato, permitindo a fragmentação entre vários sólidos. No caso da interação fluido-estrutura, as partículas de superfície do solido são tratadas como partículas de fluido e as pressões destas partículas são calculadas pela resolução da equação de Poisson para a pressão, tal como as partículas de fluido. Desta forma, o acoplamento entre solido e fluido e realizado utilizando o deslocamento e velocidade do solido elástico, como condições de contorno do fluido, e a pressão na interface, obtida pela resolução do movimento do fluido, e aplicada ao movimento do solido elástico. São apresentados e detalhados os algoritmos de solido elástico, fragmentação, colisão e acoplamento fluido-estrutura. Validações qualitativas e quantitativas do método são realizadas para casos estáticos e dinâmicos sujeitos a diferentes condições de contorno, comparando os resultados numéricos obtidos pelo MPS, outros métodos numéricos, soluções analíticas e medições experimentais presentes na literatura. / In this work a particle method to simulate the dynamics of elastic solids and fluid-structure interaction is implemented. It is based on the Moving Particle Semi-implicit Method (MPS), which was originally developed for incompressible flows with free surface. The main strategy of the MPS is to replace the differential operators of the governing equations by discrete differential operators on irregular nodes, which are derived from a model of interaction between particles. Initially details of the method and constitutive equations are shown. A simplified condition of fragmentation and collision between solids are proposed to allow the investigation of fragmentation amount multiple solids. In case of fluid-structure interaction, the solid\'s surface particles are treated as a fluid particle and the pressures of the surface particles are computed by solving Poisson equation for the pressure, just as the fluid particles. Therefore, the coupling between solid and fluid is done by using the displacement and velocity of elastic solid as the boundary conditions of the fluid, and the pressure at the interface, which is obtained when solving the fluid motion, is used to calculate the motion of the elastic solid. The algorithms for elastic solid, fragmentation, collision and fluid-structure interaction are presented and detailed. The qualitative and quantitative validations of the method are carried out herein considering static and dynamic cases subjected to deferent boundary conditions by comparing the numerical results from MPS with other numerical, analytical and experimental results available in the literature.

Um dos grandes desafios enfrentados pelos fabricantes de turbinas hidráulicas é prevenir o aparecimento de vibrações induzidas pelo escoamento nas travessas do pré-distribuidor e pás do rotor. Considerando apenas as travessas, e atribuídos a tais vibrações, foram relatados 28 casos de trincas ou ruídos anormais nas últimas décadas, que acarretaram enormes prejuízos associados a reparos, atrasos e perda de geração. O estado da arte na prevenção destes problemas baseia-se na utilização de sofisticados, e caros, programas comerciais de dinâmica dos fluidos computacional para o cálculo transiente do fenômeno. Este trabalho faz uma ampla revisão bibliográfica e levantamento de eventos de trincas ou ruídos ocorridos em travessas nos últimos 50 anos. Propõe, então, um enfoque alternativo, baseado exclusivamente em ferramentas de código aberto. A partir de hipóteses simplificadoras devidamente justificadas, o problema é formulado matematicamente de forma bidimensional, no plano da seção transversal da travessa, levando em conta a interação fluido-estrutura. Nesta estratégia, as equações de Navier-Stokes são resolvidas pelo método dos elementos finitos por meio da biblioteca gratuita oomph-lib. Um código especial em C++ é desenvolvido para o problema de interação fluido-estrutura, no qual o fenômeno de turbulência é levado em consideração por meio de um algoritmo baseado no modelo de Baldwin-Lomax. O método proposto é validado por meio da comparação dos resultados obtidos com referências e medições disponíveis na literatura, que tratam de problemas de barras retangulares suportadas elasticamente. O trabalho finaliza com a aplicação do método a um estudo de caso envolvendo uma travessa particular. / One of the biggest challenges for hydraulic turbine manufacturers is to prevent vortex-induced vibration on the stay vanes and runner blades. Only regarding stay vanes, 28 cases of cracks or unusual noises attributed to such vibrations were reported in the past decades leading to huge costs due to repair, delays and lack of generation. The state of the art today is to use powerful and expensive commercial computational fluid dynamics software to address the required transient phenomena. The present work carries out a comprehensive survey on occurred events in stay vanes during the last 50 years. Then, an alternative approach, based only on free open-source tools, is proposed. From due justified simplifying assumptions, the problem is formulated two-dimensionally, in the stay vane cross section plane, taking the fluid-structure interaction into account. In such a strategy, the Navier-Stokes equations are solved using oomph-lib, an object-oriented, finite-element library. A special C++ computational code is developed to deal with the fluid-structure interaction problem, in which turbulence is considered through a special algorithm, based on the Baldwin-Lomax model. The proposed method is validated through comparisons with an aerodynamics benchmark and an experimental measurement of oscillating rectangular bars both available in the literature. The method is finally applied to a case study of a particular stay vane. Keywords: Hydraulic turbine. Fluid-structure interaction. Vortex-induced vibration.

Orientador: Renato Pavanello / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-10T06:16:16Z (GMT). No. of bitstreams: 1 RomeroCueto_Liliana_M.pdf: 2272152 bytes, checksum: 792e6fbefc070243b8ea9c5f015ace8a (MD5) Previous issue date: 2007 / Resumo: O comportamento de sistemas acoplados fluido-estrutura pode ser modificado pela presença de defeitos estruturais. O estudo desta modificações, usando metodos numericos de simulação, permite aos engenheiros, estimar o novo estado do sistema a partir a avaliação da variação da resposta dinâmica e vibroacustica. Desta forma, a analise de sensibilidade da resposta do sistema, permite verificar a viabilidade da aplicação das tecnicas apresentadas, em tecnicas de detecção de falhas, controle de ruido e otimização estrutural. Este trabalho apresenta um estudo do comportamento vibroacustico de paineis aeronauticos reforçados e a analise de sensibilidade da resposta, destes sistemas, na presença de falhas estruturais. È elaborado um modelo estrutural baseado na teoria de cascas e vigas que é acoplado a um meio fluido acustico interior bi e tridimensional. O metodo de Elementos Finitos é usado para a solução das equações que governam o comportamento acustico-estrutural do sistema acoplado, adotando a formulação não simetrica, com os deslocamentos e pressões como variaveis incognitas dos dominios estrutural e acustico respectivamente. São apresentadas tecnicas analiticas de analise de sensibilidade da resposta em frequencia, estrutural e acustica, em relação as modificações estruturais. A precisão numerica dos metodos de sensibilidade, propostas para o caso vibroacustico, é avaliada usando-se a tecnica da Diferenças Finitas como referencia. Com a finalidade de validar os modelos implementados, são apresentados alguns testes representativos. Na modelagem de paineis aeronauticos, dois modelos de fuselagem de aeronave são apresentados. São calculadas resposta estrutural e vibroacustica e analise de sensibilidade da resposta para modelos de fuselagem 2D e 3D na presença de modificações estruturais / Abstract: The behavior of coupled fluid-structural systems can be modified due to presence of located structural defects. The study of these modifications, by simulation, using numerical methods, allows engineers to predict the new state of the system from the evaluation of structural and acoustic responses. Sensitivity analysis of acoustic and structural responses, with respect to structural geometry modifications, allow to verify the viability to apply techniques to be presented in this work for damage detection, noise control and structural optimization. In this work a formulation for vibro-acoustic sensitivity analysis of aircraft stiffener panels under structural defects is presented. Shells and beams structural elements for two and three-dimensional models coupled with an interior fluid domain are formulated. The Finite Element Method is used to solve the structural-acoustic coupled governing equations, adopting structural displacements and fluid pressure as variables. Structural and vibro-acoustic sensitivity techniques, with respect to structural modifications are presented. The accuracy of sensitivity methods is estimated using the Finite Difference technique as a reference. Respect to aircraft panels models, two fuselage models are described. Dynamic and Vibro-acoustic response and response sensitivity with respect to structural defects are calculated. / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica

A interação entre as ondas vibratórias geradas pela resposta estrutural de uma placa e os modos acústicos de uma cavidade, que é uma das principais fontes de ruído no interior de cabines de automóveis, é estudada empregando técnicas computacionais e experimentais. Para investigar este fenômeno foi construído um protótipo que consiste de uma cavidade feita de acrílico que simula o interior de um veículo. O teto deste modelo pode ser substituído por uma placa de alumínio flexível, para que se possa gerar efeitos de acoplamento entre o campo acústico e a resposta estrutural. Inicialmente, foi feito um estudo analítico do modelo, empregando-se técnicas de elementos finitos com o objetivo de extrair as freqüências naturais e modos vibroacústicos em uma faixa de freqüência abaixo de 300 Hz, em dois casos diferentes: com e sem a presença de assentos dentro da cavidade. Depois disso, o modelo foi submetido a vários testes experimentais, para se obterem as funções de resposta em freqüência, as freqüências naturais vibroacústicas e os modos acústicos e estruturais, usando excitação acústica e estrutural e empregando métodos padrões de análise modal. Para fins de comparação, as FRF\'s vibroacústicas foram obtidas de duas maneiras diferentes: com excitação estrutural e resposta acústica e com excitação acústica e resposta estrutural. As freqüências naturais e os modos vibroacústicos foram medidos usando excitação acústica e estrutural. Neste caso, é mostrada a importância do controle da força excitadora. Alguns parâmetros experimentais foram variados, tais como posicionamento das excitações acústica e estrutural e o tipo de suspensão. Finalmente, é apresentada uma compilação das diferentes técnicas para a extração dos parâmetros modais vibroacústicos, propondo novas metodologias para a medição dos modos acústicos e estruturais em freqüências abaixo de 70 Hz, controlando a força excitadora e utilizando métodos de ajuste de curvas. É mostrado que o uso de excitação acústica e estrutural merece alguns cuidados, visto que a aplicação dos métodos traz restrições de acordo com a faixa de freqüência e a escolha da melhor forma de extração dos modos vibroacústicos depende da freqüência estudada. / The interaction between the vibration waves generated by the structural response of a plate and the acoustic modes of a cavity, which is one of the main noise sources in automobile passenger cabin, is addressed employing both computational and experimental techniques. In order to investigate this phenomenon a prototype which consists of a cavity made of acrylic that resembles a vehicle interior was built. The roof of this model can be replaced by a flexible aluminum plate in order to generate coupling effects between the acoustic field and the structural response. Initially an analytical study was performed on the model employing finite element techniques with the aim of extracting the vibroacoustic natural frequencies and mode shapes in a frequency range below 300 Hz in two different cases: with and without the presence of seats inside the cavity. After that, the model was submitted to several experimental tests in order to obtain the frequency response functions, vibroacoustic natural frequencies and mode shapes, using acoustic and structural excitation and employing standard modal testing methods. For comparison purposes, the vibroacoustic FRF\'s were obtained in two different manners: with the structural excitation and acoustic response and with acoustic excitation and structural response. The vibroacoustic natural frequencies and mode shapes were measured employing both acoustic and structural excitation. In this case, it is shown the importance of the control of the exciting force. Some experimental parameters were varied, such as positioning of the acoustic and structural sources and the suspension type. Finally, a compilation of different techniques for the vibroacoustic modal parameters extraction is presented, proposing new methodologies for the acoustic and structural modes measurements in frequencies below 70 Hz, using force control and curve fitting methods. It is shown that both acoustic and structural excitations have restrictions concerning the frequency range and the choice of the better method for the extraction of the vibroacoustic mode shapes and natural frequencies depends on it.

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.

In the presented thesis work, the meshfree method with distance fields was coupled with the lattice Boltzmann method to obtain solutions of fluid-structure interaction problems. The thesis work involved development and implementation of numerical algorithms, data structure, and software. Numerical and computational properties of the coupling algorithm combining the meshfree method with distance fields and the lattice Boltzmann method were investigated. Convergence and accuracy of the methodology was validated by analytical solutions. The research was focused on fluid-structure interaction solutions in complex, mesh-resistant domains as both the lattice Boltzmann method and the meshfree method with distance fields are particularly adept in these situations. Furthermore, the fluid solution provided by the lattice Boltzmann method is massively scalable, allowing extensive use of cutting edge parallel computing resources to accelerate this phase of the solution process. The meshfree method with distance fields allows for exact satisfaction of boundary conditions making it possible to exactly capture the effects of the fluid field on the solid structure.