Orientador: João Batista Campos Silva / Resumo: Nas últimas décadas, há uma crescente preocupação em mensurar os parâmetros dinâmicos do sangue, dadas as imensas perdas de vidas por doenças cardiovasculares sendo, muitas delas, por aneurismas. A formação e desenvolvimento de um aneurisma é, predominantemente, degenerativo e provém de uma complexa interação entre os efeitos biológicos da parede arterial, os estímulos de escoamento e tensões provenientes da hemodinâmica. A tensão cisalhante na parede e o cíclico campo de pressão são um dos principais fatores responsáveis pela formação e crescimento de um aneurisma. Logo, há a necessidade de se conhecer os campos de velocidades e pressão além das tensões cisalhantes e efetivas na parede. Uma análise numérica é mais promissora que uma experimental. Uma análise experimental in-vivo possui impasses éticos e morais, sem contar a necessidade de um grande investimento. Outra vantagem de um estudo numérico é a disponibilidade de softwares livres de extração de tomografias que permite a extração da geometria sem a necessidade de um método invasivo que ocorreria em estudo experimental in vivo. Portanto com o auxílio de simulações numéricas (Ansys®), considerando o efeito multi-física de interação fluido estrutura (FSI) pela metodologia de elementos e volumes finitos foi possível verificar o efeito dos fatores que levam a formação e crescimento de aneurisma na aorta abdominal. Os aneurismas estudados foram modelos geométricos e reais sendo um dos reais obtidos através de imagens DICOM... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: In recent decades, there is growing concern in measuring the dynamic parameters of the blood, given the immense loss of life from cardiovascular disease in human history, including, aneurysms. The formation and development of an aneurysm is predominantly degenerative and results from a complex interaction between the biological effects of the arterial wall and the flow and stress effects from hemodynamics. The stress in the wall and in the cyclic field of pressure is one of the main factors for the formation and growth of an aneurysm that degenerates until its rupture. Therefore, it is necessary to know the velocity and pressure fields as well as the shear stress and effective stress on the wall. A numerical analysis is more promising than an experimental one since an in-vivo experimental analysis has ethical and moral impasses, not counting the need for a large investment. Another advantage of a numerical study is the availability of free softwares for tomography analysis that allows the extraction of the geometry without the need for an invasive method that would occur in an in vivo experimental study. Therefore, with the aid of numerical simulations (Ansys®), considering the multi-physical effect of fluid structure interaction (FSI) by the methodology of finite elements and finite volumes it was possible to verify the effect of factors that lead to the formation and growth of abdominal aortic aneurysm. The studied aneurysms came from geometric models and from real examples... (Complete abstract click electronic access below) / Mestre

Une méthode innovante pour simuler des interactions fluide-structure complexes tout en gardant un bon compromis temps de calcul/précision est présenté.Pour réduire le temps de simulation des modèles d’ordre réduits sont utilisés au lieu des modèles complets aussi bien pour les modèles structuraux que pour les modèles aérodynamiques. Un des challenges de base était d'utiliser des modèles industrielles hautes fidélités. La technique de condensation dynamique est utilisée pour réduire la taille du modèle éléments finis structures et la décomposition aux valeurs propres est utilisé sur une base de données aérodynamiques construite à partir de simulations CFD.Les non-linéarités structurelles sont réintroduites à posteriori.Une comparaison poussée des méthodes classique d'interpolation comme des méthodes de spline, d’interpolation sur des Manifold de Grassmann avec des méthodes innovantes d'apprentissage statistiques a été amené.Afin de valider complètement la méthodologie développée, une maquette expérimentale visant à imiter le comportement du carénage au sol avant le décollage a été conçue.Ce cas a pu être assimilé à une plaque avec des raideurs de liaisons dans une couche de mélange.La validation de cette méthode est réalisée en comparant les résultats des simulations numériques avec les données enregistrées pendant des essaies en soufflerie. On pourra ainsi comparer aussi bien des champs que des mesures locales. L'ensemble des essais a permis d'améliorer la compréhension de ce phénomène vibratoire qui mène à des problèmes récurrents de fatigue dans cette sous structures.Cette méthode est enfin appliquée à une structure aéronautique: les carénages de volet hypersustentateur / An innovative method for numerical simulating complex problems of fluid structure interaction, such as non-linear transients, characterized by good performances and high precision is presented in this manuscript. To cut down the simulation time, reduced order models are used for both the aerodynamic and structural modules. High fidelity industrial models have been used. A technique of dynamic condensation is employed to reduce the size of the finite element model while the technique of Proper Orthogonal Decomposition is used on a database of aerodynamic pressures built from CFD simulations. Structural non-linearities are reintroduced a posteriori. Different interpolation techniques such as the classic spline interpolation, interpolation on a Grassmann Manifold with more innovative methods of statistical learning have been compared. In order to validate the developed methodology a test campaign has been designed to reproduce a simplified mechanism of interaction inspired by a flap track fairing in take-off configuration. A plate whose stiffness depends on the springs at its attachment to the wind tunnel test section floor is immersed in a mixing layer. In parallel to the test activities a numerical model of the test rig has been developed. The validation of the methodology of fluid structure interaction is done through direct comparison between test data and simulation results. The testing activities have granted a deeper comprehension of the vibratory phenomenon that has led to recurrent fatigue problems on the impacted structures. The methodology is ultimately applied to an industrial problem: the load prediction on flap track fairings excited by engine exhaust.

Le travail engagé dans cette thèse porte sur l'étude numérique des Interactions Fluide-structure en hydrodynamique. Dans une première partie, une analyse détaillée des méthodes de couplage (schémas décalés) a été effectuée sur un cas académique. Il s'agit de la résolution de l'équation non-linéaire de Burgers dans un domaine mobile, dont I'interface mobile est représentée par un système de type masse ressort. Selon la discrétisation en temps et la linéarisation du problème couplé, on distingue quatre schémas de couplages différents : explicite, semi-implicite, implicite-externe et implicite-interne. Une étude comparative des performances en vitesse de convergence et en temps de calcul de ces schémas a été effectuée. Les performances varient suivant le schéma de couplage utilisé. Le schéma explicite permet un calcul rapide en comparaison des autres schémas. En revanche il n'assure pas la conservation de l'énergie mécanique à I'interface fluide-structure. D'où le problème de stabilité du schéma numérique. Ce problème ne se pose pas pour les algorithmes de couplage implicites, car dans ce cas la conservation de l'énergie à I'interface est assurée. Il s'agit en effet d'une condition de convergence du schéma implicite. Ce schéma requière plus de temps de calcul, mais il est nécessaire pour avoir plus de précision dans les résultats. Par ailleurs, I'analyse des déplacements de I'interface fluide-structure montre que l'écart entre la position de I'interface comme étant le bord mobile du fluide et la position de la structure, dépend principalement du schéma d'actualisation du maillage choisi.Dans une deuxième partie une extension de l'étude des algorithmes de couplage à un problème plus concret d'IFS est effectuée. Un hydrofoil en pilonnement et tangage est ainsi étudié. L'équation de la dynamique de I'hydrofoil est écrite en considérant un centre de rotation situé à une distance non nulle du centre de gravité.Ce qui rend l'équation non-linéaire et introduit un couplage des deux modes pilonnement et tangage) ainsi qu'un amortissement du tangage. La dynamique de I'hydrofoil est étudiée pour différentes configurations : en mouvement libre ou forcé, dans un fluide au repos ou en écoulement. On observe que le mouvement de I'hydrofoil est pseudo périodique amorti. L'évolution des charges hydrodynamiques suit également cette tendance et tend vers un point d'équilibre. L'étude vibratoire montre bien une modification des fréquences propres du système, qui varient suivant que le fluide est au repos ou en écoulement. Le problème est également couplé à l'équation de la position du centre de pression, qui dépend de la position de I'hydrofoil et de l'écoulement. Celle-ci présente une singularité lorsque la portance et la traînée s'annulent simultanément.Enfin Les équations prenant en compte la présence d'un fluide non-homogène à I'interface fluide-structure, du type des écoulements cavitants par poche stationnaire ou auto-oscillante, ont été développés. La méthode consiste à séparer les variables du fluide en écoulement autour d'un hydrofoil immobile d'une part et celles de l'écoulement généré par la vibration de I'hydrofoil d'autre part. Il en résulte un opérateur de masse ajoutée non symétrique en milieu non homogène et un opérateur d'amortissement ajouté dû au taux de variations de masse volumique à l’interface dans le cas auto-oscillant. L'ensemble se traduit par une modulation au cours du temps des fréquences propres et des amplitudes du système. / A numerical study of Fluid Structure Interaction (FSI) in hydrodynamic case is adressed in this thesis. Thirstly, the analysis of coupling methods (staggered schemes) was established to an academic case. It corresponds to the resolution of non linear Burgers equation in a moving domain where the moving interface is assimilated to a mass spring system. According to the time discretisation and linearization of the coupled problem, four coupling scheme can be defined : explicit, semi-implicit, implicit-outer and implicit-inner. A comparative performance study in convergence and computing time were performed. The performance depends on the coupling scheme used. The explicit scheme requires less time compared to the others schemes. However it does not allow the mechanical energy conservation at the interface, inducing the stability issue of the numerical scheme. This instabilities does not arise for the implicit coupling algorithms because the energy conservation at the interface is fulfilled. lndeed, a convergence condition is added for implicit schemes. Even though these schemes require more computing time, they are necessary to get better precision. Inter alia, the fluid-structure interface analysis shows that the gap between the interface taken as the moving boundary and the structure position mostly depends on the actualization scheme of the chosen mesh.In the second part, the coupling algorithm study is extended to physical problem of FSI. A hydrofoil in heave and pitch immersed in a fluid flow is then studied. The equation of hydrofoil movement takes account the distance between the rotation center and the center of gravity. This causes the equation to be nonlinear and introduces a coupling of the two movements (heave and pitch) and a damping of the heave movement. The hydrofoil dynamic is studied for different configurations : forced movements or not, immersed in a fluid at rest or a flowing one. It shows that the hydrofoil movement is pseudo-periodic followed by a damping movement. The hydrodynamic forces tend to follow the same evolution and converge to an equilibrium point. The vibration study clearly shows a frequency modification of the system that depends on the fluid flow (at rest or with an inflow). The problem is also coupled to center of pressure position's equation which depends on the hydrofoil position and the fluid flow. The trend of the position presents a singularity when the lift and drag coefficients vanishes at the same time.Last part, the equation that take into account the inhomogeneous characteristic of the fluid at the fluid-structure interface as well as sheet cavitation in steady or unsteady case, was developed. The method allows the separation of the fluid variables when flowing around the fixed hydrofoil on one hand and the flow generated by the hydrofoil vibration one the other. This introduces an asymmetric added mass operator and an added damping operation due to the variation of the density of the fluid at the interface in unsteady case.The whole system results in a natural frequencies and amplitudes modulation over time.

Dans le contexte de lutte contre les nuisances sonores, cette thèse porte sur le développement de méthodes de résolution efficaces par éléments finis, pour des problèmes de vibroacoustique interne avec interfaces dissipatives, dans le domaine des basses fréquences. L’étude se limite à l’utilisation de solutions passives telles que l’intégration de matériaux poreux homogènes et isotropes, modélisés par une approche fondée sur la théorie de Biot-Allard. Ces modèles étant coûteux en terme de résolution, un des objectifs de cette thèse est de proposer une approche modale pour la réduction du problème poroélastique, bien que l’adéquation d’une telle approche avec le comportement dynamique des matériaux poreux soit à démontrer. Dans un premier temps, la résolution de problèmes couplés élasto-poro-acoustiques par sous-structuration dynamique des domaines acoustiques et poreux est établie. L’approche modale originale proposée pour les milieux poroélastiques, ainsi qu’une procédure de sélection des modes significatifs, sont validées sur des exemples 1D à 3D. Une deuxième partie présente une méthode combinant l’utilisation des modèles réduits précédemment établis avec une procédure d’approximation de solution par approximants de Padé. Il est montré qu’une telle combinaison offre la possibilité d’accroître les performances de la résolution (allocation mémoire et ressources en temps de calcul). Un chapitre dédié aux applications permet d’évaluer et comparer les approches sur un problème académique 3D, mettant en valeur leurs performances encourageantes. Afin d’améliorer les méthodes établies dans cette thèse, des perspectives à ces travaux de recherche sont apportées en conclusion. / In the context of interior noise reduction, the present work aims at proposing Finite Element (FE) solution strategies for interior structural-acoustic applications including 3D modelling of homogeneous and isotropic poroelastic materials, under timeharmonic excitations, and in the low frequency range. A model based on the Biot-Allard theory is used for the poroelastic materials, which is known to be very costly in terms of computational resources. Reduced models offer the possibility to enhance the resolution of such complex problems. However, their applicability to porous materials remained to be demonstrated.First, this thesis presents FE resolutions of poro-elasto-acoustic coupled problems using modal-based approaches both for the acoustic and porous domains. The original modal approach proposed for porous media, together with a dedicated mode selection and truncation procedure, are validated on 1D to 3D applications.In a second part, modal-reduced models are combined with a Padé approximants reconstruction scheme in order to further improve the efficiency.A concluding chapter presents a comparison and a combination of the proposed methods on a 3D academic application, showing promising performances. Conclusions are then drawn to provide indications for future research and tests to be conducted in order to further enhance the methodologies proposed in this thesis.

L'objet de cette thèse est l'étude de la stabilisation de modèles d'interaction fluide-structure par des contrôles de dimension finie agissant sur la frontière du domaine fluide. L'écoulement du fluide est décrit par les équations de Navier-Stokes incompressibles tandis que l'évolution de la structure, située à la frontière du domaine fluide, satisfait une équation d'Euler-Bernoulli avec amortissement. Dans le chapitre 1, nous étudions le cas où le contrôle est une condition aux limites de Dirichlet sur les équations du fluide (contrôle par soufflage/aspiration). Nous obtenons des résultats de stabilisation locale du système non-linéaire autour d'une solution stationnaire instable de ce système. Dans les chapitres 2 et 3, nous nous intéressons au cas où le contrôle est une force appliquée sur la structure (contrôle par déformation de paroi). Dans le chapitre 2, nous considérons un modèle simplifié, où l'équation d'Euler-Bernoulli pour la structure est remplacée par un système de dimension finie. Nous construisons des lois de contrôle pour les systèmes de dimension infinie, ou pour leurs approximations semi-discrètes, capables de stabiliser les systèmes linéarisés avec un taux de décroissance exponentielle prescrit, et localement les systèmes non-linéaires. Nous présenterons des résultats numériques permettant de vérifier l'efficacité de ces lois de contrôles. / The aim of this thesis is to study the stabilization of fluid-structure interaction models by finite dimensional controls acting at the boundary of the fluid domain. The fluid flow is described by the incompressible Navier-Stokes equations while the displacement of the structure, localized at the boundary of the fluid domain, satisfies a damped Euler-Bernoulli beam equation. First, we study the case where the control is a Dirichlet boundary condition in the fluid equations (control by suction/blowing). We obtain local feedback stabilization results around an unstable stationary solution of this system. Chapters 2 and 3 are devoted to the case where control is a force applied to the structure (control by boundary deformation). We consider, in chapter 2, a simplified model where the Euler-Bernoulli equation for the structure is replaced by a system of finite dimension. We construct feedback control laws for the infinite dimensional systems, or for their semi-discrete approximations, able to stabilize the linearized systems with a prescribed exponential decay rate, and locally the nonlinear systems. We present some numerical results showing the efficiency of the feedback laws.

Made available in DSpace on 2014-06-11T19:27:13Z (GMT). No. of bitstreams: 0 Previous issue date: 2010-02-26Bitstream added on 2014-06-13T20:16:16Z : No. of bitstreams: 1 souza_ap_me_ilha.pdf: 2793163 bytes, checksum: cac24da914ff0917c1d9438da0bc50f8 (MD5) / A Bioengenharia está cada vez mais presente em todo o mundo. Ela estuda problemas que tenham dificuldades em análise experimental em laboratório. Este trabalho determina a força exercida por um grampo cirúrgico quando o mesmo é aplicado em determinada região Arterial a fim de ocluir o fluxo sanguíneo. Esta força foi medida utilizado uma simulação computacional. Duas áreas foram consideradas: uma representando a parede Arterial e a outra, representando o sangue. Esta simulação realizou-se em regime de acoplamento entre os dois domínios utilizando interação fluido-estrutura. A modelagem foi feita considerando dois domínios distintos: estrutura e fluido. O fluido é considerado incompressível e Newtoniano e é governado pelas equações de Navier-Stokes. As paredes da estrutura são modeladas a partir da Lei de Hooke. A solução numérica calcula: os campos de pressão e velocidade do fluido, campo de deslocamento da estrutura e a força aplicada pelo grampo para que ocorra a obstrução do fluxo sanguíneo naquele local / The Bioengineering is increasingly present in the fields of scientific research throughout the world. It studies problems that have difficulty in experimental analysis in the laboratory. This work measures the force exerted by a surgical clip when it is applied in a given region Arterial to occlude blood flow. This force was measured using a computational simulation. It was modeled two cylinders, one representing the arterial wall and the other in this first, representing the blood. This simulation was carried out under the coupling between the two domains using fluid-structure interaction. Modeling was done in three-dimensions, considering two distinct areas: one as structure and other as fluid. The fluid was considered incompressible and Newtonian. It is governed by the Navier-Stokes equations. The walls of the structure were modeled from the Hooke's Law. The numerical solution calculate: pressure fields, and fluid velocity, displacement field of the structure and the force applied by the clip for the occurrence of obstruction of blood flow there

Submitted by Jeane Batista de Carvalho null (carvalhojeanetiagina@yahoo.com.br) on 2017-05-23T15:08:27Z No. of bitstreams: 1 Dissertação_CarvalhoJeane.pdf: 8189158 bytes, checksum: ebb9bc23a7d48fbc4154db51680f263e (MD5) / Approved for entry into archive by Luiz Galeffi (luizgaleffi@gmail.com) on 2017-05-23T18:27:50Z (GMT) No. of bitstreams: 1 carvalho_jb_me_ilha.pdf: 8189158 bytes, checksum: ebb9bc23a7d48fbc4154db51680f263e (MD5) / Made available in DSpace on 2017-05-23T18:27:50Z (GMT). No. of bitstreams: 1 carvalho_jb_me_ilha.pdf: 8189158 bytes, checksum: ebb9bc23a7d48fbc4154db51680f263e (MD5) Previous issue date: 2017-03-24 / Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) / Nas últimas décadas, há uma crescente preocupação em mensurar os parâmetros dinâmicos do sangue, dadas as imensas perdas de vidas por doenças cardiovasculares sendo, muitas delas, por aneurismas. A formação e desenvolvimento de um aneurisma é, predominantemente, degenerativo e provém de uma complexa interação entre os efeitos biológicos da parede arterial, os estímulos de escoamento e tensões provenientes da hemodinâmica. A tensão cisalhante na parede e o cíclico campo de pressão são um dos principais fatores responsáveis pela formação e crescimento de um aneurisma. Logo, há a necessidade de se conhecer os campos de velocidades e pressão além das tensões cisalhantes e efetivas na parede. Uma análise numérica é mais promissora que uma experimental. Uma análise experimental in-vivo possui impasses éticos e morais, sem contar a necessidade de um grande investimento. Outra vantagem de um estudo numérico é a disponibilidade de softwares livres de extração de tomografias que permite a extração da geometria sem a necessidade de um método invasivo que ocorreria em estudo experimental in vivo. Portanto com o auxílio de simulações numéricas (Ansys®), considerando o efeito multi-física de interação fluido estrutura (FSI) pela metodologia de elementos e volumes finitos foi possível verificar o efeito dos fatores que levam a formação e crescimento de aneurisma na aorta abdominal. Os aneurismas estudados foram modelos geométricos e reais sendo um dos reais obtidos através de imagens DICOM. / In recent decades, there is growing concern in measuring the dynamic parameters of the blood, given the immense loss of life from cardiovascular disease in human history, including, aneurysms. The formation and development of an aneurysm is predominantly degenerative and results from a complex interaction between the biological effects of the arterial wall and the flow and stress effects from hemodynamics. The stress in the wall and in the cyclic field of pressure is one of the main factors for the formation and growth of an aneurysm that degenerates until its rupture. Therefore, it is necessary to know the velocity and pressure fields as well as the shear stress and effective stress on the wall. A numerical analysis is more promising than an experimental one since an in-vivo experimental analysis has ethical and moral impasses, not counting the need for a large investment. Another advantage of a numerical study is the availability of free softwares for tomography analysis that allows the extraction of the geometry without the need for an invasive method that would occur in an in vivo experimental study. Therefore, with the aid of numerical simulations (Ansys®), considering the multi-physical effect of fluid structure interaction (FSI) by the methodology of finite elements and finite volumes it was possible to verify the effect of factors that lead to the formation and growth of abdominal aortic aneurysm. The studied aneurysms came from geometric models and from real examples, being one of the real ones obtained through DICOM images. / CNPq: 131153/2015-3

Made available in DSpace on 2014-12-17T14:58:00Z (GMT). No. of bitstreams: 1 RairamFCA_DISSERT2.pdf: 4758176 bytes, checksum: bfb1653549a50848b4721bda9a78bd6e (MD5) Previous issue date: 2010-03-26 / The pumping through progressing cavities system has been more and more employed in the petroleum industry. This occurs because of its capacity of elevation of highly viscous oils or fluids with great concentration of sand or other solid particles. A Progressing Cavity Pump (PCP) consists, basically, of a rotor - a metallic device similar to an eccentric screw, and a stator - a steel tube internally covered by a double helix, which may be rigid or deformable/elastomeric. In general, it is submitted to a combination of well pressure with the pressure generated by the pumping process itself. In elastomeric PCPs, this combined effort compresses the stator and generates, or enlarges, the clearance existing between the rotor and the stator, thus reducing the closing effect between their cavities. Such opening of the sealing region produces what is known as fluid slip or slippage, reducing the efficiency of the PCP pumping system. Therefore, this research aims to develop a transient three-dimensional computational model that, based on single-lobe PCP kinematics, is able to simulate the fluid-structure interaction that occurs in the interior of metallic and elastomeric PCPs. The main goal is to evaluate the dynamic characteristics of PCP s efficiency based on detailed and instantaneous information of velocity, pressure and deformation fields in their interior. To reach these goals (development and use of the model), it was also necessary the development of a methodology for generation of dynamic, mobile and deformable, computational meshes representing fluid and structural regions of a PCP. This additional intermediary step has been characterized as the biggest challenge for the elaboration and running of the computational model due to the complex kinematic and critical geometry of this type of pump (different helix angles between rotor and stator as well as large length scale aspect ratios). The processes of dynamic generation of meshes and of simultaneous evaluation of the deformations suffered by the elastomer are fulfilled through subroutines written in Fortan 90 language that dynamically interact with the CFX/ANSYS fluid dynamic software. Since a structural elastic linear model is employed to evaluate elastomer deformations, it is not necessary to use any CAE package for structural analysis. However, an initial proposal for dynamic simulation using hyperelastic models through ANSYS software is also presented in this research. Validation of the results produced with the present methodology (mesh generation, flow simulation in metallic PCPs and simulation of fluid-structure interaction in elastomeric PCPs) is obtained through comparison with experimental results reported by the literature. It is expected that the development and application of such a computational model may provide better details of the dynamics of the flow within metallic and elastomeric PCPs, so that better control systems may be implemented in the artificial elevation area by PCP / O sistema de bombeamento por cavidades progressivas est? sendo cada vez mais empregado na ind?stria do petr?leo, devido ? sua capacidade de eleva??o de ?leos altamente viscosos ou de fluidos com grandes concentra??es de areia ou outras part?culas s?lidas. Uma Bomba de Cavidades Progressivas (BCP) ? composta, basicamente, por um rotor - uma pe?a met?lica de forma semelhante a um parafuso exc?ntrico, e um estator - um tubo de a?o revestido internamente por uma h?lice dupla, a qual pode ser r?gida ou deform?vel/elastom?rica. Em geral, uma BCP ? submetida a uma combina??o de press?o do po?o com press?o gerada pelo pr?prio processo de bombeio. Em BCPs elastom?ricas, essa combina??o de esfor?os comprime o estator, gerando ou aumentando a folga existente entre o rotor e o estator, reduzindo, portanto, o efeito de veda??o entre suas cavidades. Tal abertura da regi?o de selagem produz o que ? conhecido como escorregamento do fluido, diminuindo, com isso, a efici?ncia de sistema de bombeio por BCP. Dessa maneira, este trabalho se prop?e a desenvolver um modelo computacional tridimensional transiente do processo din?mico da intera??o fluido-estrutural (FSI) que ocorre no interior de BCPs met?licas e elastom?ricas. O objetivo principal ? avaliar, a partir do uso do modelo desenvolvido, as caracter?sticas din?micas de efici?ncia de bombeio por BCPs, em fun??o de informa??es locais e instant?neas detalhadas dos campos de velocidade, press?o e deforma??o no seu interior. Para o alcance de tais metas (desenvolvimento e uso do modelo), fez-se necess?rio o desenvolvimento de uma metodologia pr?pria para gera??o de malhas computacionais din?micas, m?veis e deform?veis, representando as regi?es fluida e estrutural de uma BCP. Tal procedimento caracterizou-se como o maior desafio para a elabora??o do modelo computacional, devido ? cinem?tica complexa e ? geometria cr?tica desse tipo de bomba (?ngulos de h?lice diferentes entre rotor e estator e grandes diferen?as de escala de comprimento). Os processos de gera??o din?mica das malhas e de avalia??o simult?nea das deforma??es sofridas pelo elast?mero s?o realizados atrav?s de sub-rotinas em linguagem Fortran 90, as quais interagem dinamicamente com o software de din?mica dos fluidos computacional CFX/ANSYS. Desde que o modelo linear el?stico ? empregado para avaliar as deforma??es elastom?ricas, n?o ? necess?rio usar nenhum software para an?lise estrutural. Entretanto, uma proposta inicial para simula??o din?mica no ANSYS empregando-se modelos constitutivos hiper-el?sticos para o elast?mero ? tamb?m apresentada no presente trabalho. A valida??o dos resultados produzidos com a presente metodologia (gera??o de malha, simula??o do escoamento em BCPs met?licas e simula??o da intera??o fluido-estrutural em BCPs elastom?ricas) ? obtida atrav?s da compara??o com resultados experimentais reportados pela literatura. Vislumbra-se que o desenvolvimento e aplica??o de tal ferramenta computacional poder?o fornecer maiores detalhes da din?mica do escoamento no interior de BCPs met?licas e elastom?ricas, de maneira que melhores sistemas de controle possam ser implementados na ?rea de eleva??o artificial por BCPs

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / In the present thesis it is proposed and evaluated, both numerically and experimentally, an inverse procedure for the indirect determination of axial loads applied to submersed pipe-like structures, based on their dynamic responses. The investigation is motivated by the existence of practical problems encountered in the oil industry. An experimental bench has been designed and built, consisting in a reservoir inside which a tubular stainless steel beam has been mounted and tested. Special fixtures have been designed in such a way to enable to apply controlled axial loads and represent different types of boundary conditions. In parallel, computational routines have been developed for the two-dimensional modeling of the structure accounting for the effects of axial loads, flexible supports and fluid-structure interaction, based on the finite element approach. Having in mind the difficulties which are expected to be encountered when the methodology be applied in real conditions, some special dynamic test procedures have been considered, including Operational Modal Analysis (OMA), which enables to identify modal parameters from output-only measurements. Numerous scenarios have been considered using either numerically simulated or experimentally measured responses. As for the resolution of the inverse problem, two strategies have been investigated: the first consists in the deterministic resolution of a constrained optimization problem based on evolutionary algorithms, and the second, which enables to account for the presence of uncertainties in the experimental data, is a stochastic approach based on Bayesian inference, combined with Markov chains and Metropolis-Hastings algorithm. The results obtained confirm the operational feasibility and satisfactory accuracy provided by the suggested identification approaches. / Na presente tese de doutorado é proposto e avaliado numérica e experimentalmente um procedimento inverso para determinação indireta de carregamentos axiais aplicados a estruturas tubulares submersas a partir de suas respostas dinâmicas. A investigação é motivada pela existência de problemas práticos evidenciados pelo setor de tecnologia submarina da indústria petrolífera. Nesta proposta, as cargas axiais, que na prática não podem ser medidas diretamente, são identificadas através da resolução de um problema inverso, formulado como um problema de otimização, a partir das respostas dinâmicas da estrutura. Uma bancada experimental foi projetada e construída, composta de um reservatório dentro do qual foi ensaiado um tubo metálico de seção circular. Mecanismos de fixação e aplicação de carga à estrutura foram especialmente projetados de modo a permitir consideração de dois tipos diferentes de condição de contorno. Paralelamente, rotinas computacionais foram desenvolvidas para a modelagem numérica bidimensional da estrutura incluindo os efeitos de interação fluido-estrutura e das cargas axiais, com base no Método de Elementos Finitos. Tendo em vista o objetivo da aplicação da metodologia proposta em situações práticas, as quais envolvem dificuldades de execução de ensaios em ambientes submarinos, foram investigados procedimentos de ensaios dinâmicos especialmente adaptados a estas condições. Com este intuito, foi analisado o emprego da técnica de análise modal experimental denominada OMA (Operational Modal Analysis), que permite obter os parâmetros modais sem conhecimento das forças de excitação da estrutura. Numerosos cenários de identificação foram estudados utilizando tanto respostas dinâmicas simuladas numericamente, quanto respostas medidas experimentalmente. Visando considerar a influência de incertezas nos dados experimentais, o problema de identificação da carga axial também foi tratado utilizando uma abordagem estocástica, com base em inferência bayesiana, a partir da simulação de cadeias de Markov, associada ao algoritmo Metropolis-Hastings. Os resultados obtidos atestam a viabilidade operacional e a precisão satisfatória do procedimento de identificação proposto. / Doutor em Engenharia Mecânica

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / In this work, the combined multi-direct forcing and immersed boundary method (IBM) were presented to simulate uid-structure interaction problems. The multi-direct forcing is used aim at satisfying the no-slip condition in the immersed boundary. For the numerical simulations was used a multi-purpose computer code that is being developed in the MFlab - Fluid Mechanics Laboratory of Federal University of Uberl^andia. Tests are made to validate the numerical schemes and routines were implemented to simulate uid-structures interaction problems. Furthermore, computational tools are developed to construct and manage and optimize the use of a Beowulf cluster where all the parallel simulations presented in this work were done. The Method of Manufactured Solutions has been used for order-of-accuracy verication in the computational uid dynamics code. Two uid-structure interaction problems were studied using this methodology. The rst is a ow over a sphere for some Reynolds numbers. The results were compared to empirical results, obtaining satisfactory approximations. The second one is a immersed simple pendulum. For this problem the results are in agreement with physics. Indeed, these are preliminar results. New tests must be done to make progress in the methodology. Improvements are proposed in the IBM, in the uid-structure model, in the turbulence model, in the method used to discretize the uid domain. It is also proposed to apply the methodology to real problems as risers and valves. / O presente trabalho tem como principal objetivo a aplicação do método multifoçagem (MMF) para solução numérica tridimensional de problemas de interação uidoestrutura, buscando-se garantir a condição de não-escorregamento na região da fronteira imersa. Para as simulações numéricas foi utilizado um código computacional multipropósito em desenvolvimento no MFlab - Laboratório de Mecânica dos Fluidos da Universidade Federal de Uberlândia. Foram feitas modificações nesse código para que se pudesse validá-lo para solução de problemas com fronteira imersa e foi implementada uma rotina para solução de um problema de interação uido-estrutura total. Além disso, foi desenvolvido um pacote de ferramentas computacionais que possibilitou instalar e melhorar o desempenho de um cluster do tipo Beowulf utilizado para o desenvolvimento das simulações num eriças em paralelo do presente trabalho. Utilizando o Método das Soluções Manufaturadas foram obtidas soluções sintetizadas para as equações de Navier-Stokes, o que possibilitou obter a ordem de convergência numérica do código computacional para problemas contínuos e a validação deste código para problemas envolvendo corpos imersos ao combinar a o método das soluções manufaturadas com a metodologia de fronteira imersa. Na sequência foi solucionado o problema de escoamento ao redor de uma esfera parada, cujos resultados foram comparados com referencias empíricas, obtendo-se boa aproximação. Ainda para esse caso foi feita a avalição da norma L2 para as soluções num eriças obtidas nos pontos lagrangianos verificando a garantia da condição de não-escorregamento e feita uma análise da inuência dos número de ciclos utilizados no método multi-forçagem. Foi vericado que a solução numérica obtida depende do número de ciclos o que faz com que seja necessário se estabelecer um critério de convergência para este método. Um segundo problema de interação uido-estrutura total foi estudado. Consiste em um pêndulo simples imerso em um uido que parte de uma dada posição angular inicial e oscila em torno da sua posição de equilíbrio, até parar. Para esse caso foram feitas análises quantitativas. Os resultados são preliminares mas coerentes com a física do problema, indicando que a metodologia é adequada para solução deste tipo de problema. / Mestre em Engenharia Mecânica