O presente trabalho tem por objetivo desenvolver uma formulação numérica baseada em Análise Isogeométrica para o estudo de problemas de interação fluido-estrutura (IFE) em aplicações envolvendo corpos rígidos submersos, onde escoamentos incompressíveis de fluidos Newtonianos com superfície livre são considerados. Propõe-se o emprego da Análise Isogeométrica por permitir a unificação entre os procedimentos de pré-processamento e análise, melhorando assim as condições de continuidade das funções de base empregadas tanto na discretização espacial do problema como na aproximação das variáveis do sistema de equações. O sistema de equações fundamentais do escoamento é formado pelas equações de Navier-Stokes e pela equação da conservação de massa, descrita segundo a hipótese de pseudo-compressibilidade, em uma formulação cinemática ALE (Arbitrary Lagrangean- Eulerian). A consideração da superfície livre no escoamento se dá tratando o fluido como um meio bifásico, através do método Level Set. O corpo rígido apresenta não linearidade na rotação e restrições representadas por vínculos elásticos e amortecedores viscosos, sendo a equação de equilíbrio dinâmico resolvida através do método de Newmark. O esquema de acoplamento sólido-fluido adotado é o particionado convencional, que impõe condições de compatibilidade cinemáticas e de equilíbrio sobre a interface sólido-fluido, analisando ambos os meios de maneira sequencial. A discretização das equações governantes é realizada através do esquema explícito de dois passos de Taylor-Galerkin, aplicado no contexto da Análise Isogeométrica. Por fim, são analisados alguns problemas da Dinâmica de Fluidos Computacional, de onde se concluiu que os resultados obtidos são bastante consistentes com os fenômenos envolvidos, com as ferramentas exclusivas da Análise Isogeométrica, como o refinamento k, melhorando a convergência dos resultados. Para escoamentos bifásicos, verificou-se que o método Level Set obteve resultados bastante promissores apresentando, entretanto, uma dissipação numérica excessiva. Propõe-se, para estudos futuros, a elaboração de esquemas numéricos que conservem melhor o volume da fase líquida do escoamento. / The present work aims to development of a numerical formulation based on Isogeometric Analysis for the study of Fluid-Structure Interaction problems in applications involving rigid bodies submerged, considering incompressible Newtonian flows with free surface. The use of the Isogeometric Analysis allows unification between the preprocessing and analysis steps, improving then the continuity of the base functions employed, both in the spatial discretization and approximation of the variables in the system of equations. The fundamental flow equations are formed by the Navier-Stokes and the mass conservation, described by de pseudo-compressibility hypothesis, in an ALE (Arbitrary Lagrangean-Eulerian) kinematic formulation. The free surface consideration of the flow is handled treating the fluid like a two- phase medium, using the Level Set method. The rigid body considers nonlinearity in rotation, and restrictions represented by elastic springs and viscous dampers, with the dynamic equilibrium equation being resolved using the Newmark’s method. The solid-fluid coupling scheme is the conventional partitioned, which imposes kinematics and equilibrium compatibility conditions on the solid-fluid interface, analyzing both mediums in a sequential manner. The governing equations are discretized using the explicit two step Taylor-Galerkin method, applied in an Isogeometric Analisys context. Finally, some Computational Fluid Dinamics problems are analysed, from which it was concluded that the results obtained are quite consistent with phenomena involved, with the unique tools of Isogeometric Analysis, such as k-refinement, improving the convergence of the results. For biphasic flows, it was verified that the Level Set method obtained very promising results, presenting, however, an excessive numerical dissipation. For future studies, it is proposed the elaboration of numerical schemes that better preserve the volume of the liquid phase of the flow.

Orientador: Prof. Dr. André Fenili / Dissertação (mestrado) - Universidade Federal do ABC, Programa de Pós-Graduação em Engenharia Mecânica, 2015. / No presente trabalho, é investigado o modelo matemático de uma viga não linear com seção retangular em rotação. Foi considerado, no modelo, curvatura linear e condições de contorno de viga livre-engastada. Para simulação, considera-se ainda ângulo de ataque constante e as não linearidades resultantes do acoplamento entre a velocidade angular do eixo de rotação e vibração transversal da viga. As forças externas de arrasto e sustentação ao longo do comprimento da viga foram modeladas utilizando teoria de faixas. Estas forças são funções não lineares da velocidade transversal da viga e da velocidade angular do eixo de rotação. As equações de governo para a viga foram discretizadas no espaço utilizando-se o método das diferenças finitas. Isso resulta em uma sistema de equações diferenciais para os valores nodais da discretização espacial da viga. Esse conjunto de equações diferenciais é integrado numericamente no domínio do tempo utilizando-se um método de Runge-Kutta de quarta ordem. / In this work, we investigate the mathematical model of a nonlinear beam with rectangular cross section in rotation. It was considered, in the model, linear curvature and boundary conditions of clamped-free. To simulation, still considered constant attack angle and the nonlinearities resulting from the coupling between the angular velocity of the rotation axis and transverse vibration of the beam. The external drag and lift forces along the length of the beam were modeled as turbulent drag and potential generated by the lift force. These forces are dependent on the speed of nonlinear excitations acting on the beam. The solution in space wasobtained from the finite difference method using 5 points, and the solution in the time from the Runge- Kutta method of 4th order.

Orientador: Prof. Dr. Juan Pablo Julca Avila / Dissertação (mestrado) - Universidade Federal do ABC, Programa de Pós-Graduação em Engenharia Mecânica, 2015. / Atualmente, os risers rígidos em catenária ou SCRs (Steel Catenary Risers), que são tubos longos de aço, apresentam-se como a melhor solução técnico-econômica na transferência de petróleo e gás desde o solo marinho até uma plataforma flutuante. Os SCRs são de fácil fabricação, resistem altas pressões internas e hidrostáticas e também resistem altas temperaturas. Porém, cuidado especial deve-se tomar no cálculo dos momentos fletores e força axial interna no ponto de contato com o solo marinho, sendo estes parâmetros cruciais no projeto. Por outro lado, devido a que os SCRs interagem com o solo marinho, a plataforma à qual está conectada, correntezas e com o escoamento interno, a teoria de SCRs é complexa e não tem sido totalmente desenvolvida, requerendo para seu estudo a teoria de vigas curvas, tópicos de mecânica dos sólidos e dos fluidos, dinâmica não-linear, mecânica de ondas e mecânica dos solos. Este trabalho tem como principal objetivo o desenvolvimento e implementação de uma ferramenta computacional para análise estática e dinâmica bidimensional de risers rígidos e flexíveis dispostos em catenária em contato com o solo marinho. A discretização espacial do riser é feita usando elementos finitos não lineares tipo de viga, incluindo grandes deslocamentos e rotações. A formulação co-rotacional é utilizada para o tratamento da não linearidade geométrica. O método iterativo-incremental de Newton-Raphson é usado para resolver as equações de equilíbrio estático e dinâmico. A integração no tempo das equações dinâmicas é feita usando o esquema implícito de Newmark. A fim de garantir a estabilidade do esquema numérico implementado quando são impostos deslocamentos no topo do riser pelo método de penalização, é introduzido nas equações dinâmicas um termo de amortecimento estrutural para a filtragem das frequências espúrias induzidas por este tipo de excitação. O solo marinho é modelado como uma fundação elástica-linear do tipo Winkler e o método de penalização é usado para a imposição da condição de não penetração. Simulações estáticas e dinâmicas de problemas geometricamente não lineares foram conduzidas para a avaliação do elemento de viga plana implementado neste trabalho. Os resultados obtidos foram comparados com resultados da literatura para a validação do código. A ferramenta computacional foi aplicada satisfatoriamente para resolver problemas estáticos e dinâmicos de risers rígidos e flexíveis. / Steel catenary risers (SCR) are slender steel pipes that hang free in the ocean, this represents the best technical and economical solution for the oil and gas transfer from the seabed to the floating platform. SCRs are of easy manufacturing, high internal and external pressure resistance and also high temperature resistance. Special care should be taken in the calculation of the bending and tension stresses at the touch down point (TDP) as this parameters are of main importance in the calculation of fatigue resistance. On the other side, as the riser interacts with many other elements as seawater currents, internal flow, floating platform and seabed SCRs theory is complex and is not yet well developed, requiring for its study deep knowledge of curved beam theory, solid and fluid mechanics, non-linear dynamics, wave theory and soil mechanics. The main objective of this work is the development and implementation of a computational tool for the static and dynamic two-dimensional analysis of steel catenary and flexible risers, special attention is given to the seabed contact phenomena, to this end, numerical methods for the solution of dynamic equations were implemented into a MATLAB code. The spacial discretization of the riser geometric domain was made by finite element procedures, the large deflections and rotations, inherent to risers geometric non linearity, were treated by means of the co-rotational formulation. The incremental-iterative Newton-Raphson scheme is used to solve the equations of static and dynamic equilibrium. Time domain integration is made using Newmarks implicit method. To guarantee the numerical stability of the implemented code when imposed a time-varying nodal displacement by the penalty method an structural damping is introduced. This damping filters spurious frequencies induced by the penalty method. The seabed is modeled as an elastic foundation of Winkler type, once again the penalty method is used to enforce the non-penetration condition. Static and dynamic simulations of beams with geometrical non linearity were conducted in order to test the stability and accuracy of the implemented code. These results were compared with those available in specialized literature in order to validate the code. This computational tool was successfully applied to the static and dynamic analysis of steel catenary risers.

Contrairement aux liquides ordinaires, les fluides complexes comme le sang exhibent des comportements étranges qui dépendent essentiellement des structures sous-jacentes qui les composent. La simulation des écoulements sanguins continue de poser un formidable défi pour les modélisations théoriques et numériques dont l'intérêt est de développer des méthodes et des outils de simulation pour la communauté médicale. Nous proposons dans cette thèse une contribution à ce projet qui sera majoritairement centré sur les aspects numériques et informatiques. Nous nous sommes particulièrement intéressés à l'interaction entre le sang et la paroi vasculaire, qui joue un rôle important dans les grandes artères comme l'aorte. Nous nous sommes aussi investis dans la simulation du transport des cellules sanguines dans le sang. Pour la résolution des équations aux dérivées partielles décrivant nos modèles d'hémodynamique, nous avons choisi d'utiliser des méthodes numériques dont la précision pourra être accrue de manière arbitraire. Dans ce but, les principaux ingrédients qui ont été mis en oeuvre sont (i) la méthode des éléments finis basée sur des approximations de Galerkin d'ordre arbitraire en espace et géométrie, (ii) la méthode ALE pour la prise en compte de la mobilité des domaines pour des déplacements d'ordre arbitraire, (iii) les couplages implicites et semi-implicites pour l'interaction fluide-structure. Nous proposons également une nouvelle formulation de la méthode de la frontière élargie visant à modéliser le transport de particules déformables immergées dans un fluide. Nos simulations numériques se sont appuyées sur la librairie de calcul Feel++, spécialisée dans la résolution d'EDP. Outre l'implémentation des modèles physiques, nous y avons développé diverses fonctionnalités nécessaires à la mise en oeuvre de nos méthodes : interpolation, méthode de Galerkin non standard, méthode ALE, environnement pour l'interaction fluide-structure. De plus, de par la taille des géométries et la complexité des modèles mis en jeu, le passage au calcul parallèle a été indispensable pour pouvoir réaliser nos simulations. Ainsi, nous avons décrit le développement qui a été effectué dans cette librairie pour permettre le déploiement de nos programmes sur des architectures parallèles. / Unlike ordinary liquids, complex fluids such as blood exhibit strange behavior mainly dependent underlying structures that compose them. Simulation of blood flow remains a formidable challenge for theoretical and numerical modeling whose interest is to develop methods and simulation tools for the medical community. We propose in this work a contribution to this project will be mainly focused on numerical and computational aspects. We are particularly interested in the interaction between the blood and the vascular wall, which plays an important role in the large arteries as the aorta. We also invested in transport simulation of blood cells in the blood. For solving partial differential equations describing our hemodynamic models we chose to use numerical methods whose accuracy can be increased arbitrarily. For this purpose, the main ingredients that have been used are textit {(i)} the finite element method based on Galerkin approximations of arbitrary order in space and geometry, (i) the ALE method for taking into account the mobility of areas for movements of arbitrary order, (ii) the implicit and semi-implicit coupling for fluid-structure interaction. We also propose a new formulation of the method of boundary extended to model the transport of deformable particles immersed in a fluid. Our numerical simulations were based on the library Feel++ specializing on the EDP resolution. In addition to implementing physical models, we have developed various features are necessary for the implementation our methods: interpolation, Galerkin non-standard method, ALE method, framework for fluid-structure interaction. In addition, the size and geometry of the complexity of the models involved, the transition to parallel computing was essential in order to achieve our simulations. Thus, we have described the development was carried out in the library to allow deployment of our programs on parallel architectures.

Submitted by Kilvya Braga (kilvyabraga@hotmail.com) on 2018-04-27T11:35:02Z No. of bitstreams: 1 DIEGO PAES DE ANDRADE PEÑA - DISSERTAÇÃO (PPGEM) 2016.pdf: 5093848 bytes, checksum: c6e79e54502ec5a0ff9ff5b410ffd362 (MD5) / Made available in DSpace on 2018-04-27T11:35:02Z (GMT). No. of bitstreams: 1 DIEGO PAES DE ANDRADE PEÑA - DISSERTAÇÃO (PPGEM) 2016.pdf: 5093848 bytes, checksum: c6e79e54502ec5a0ff9ff5b410ffd362 (MD5) Previous issue date: 2016-09-02 / CNPq / A aeroelasticidade é o campo da ciência que estuda a correlação entre as forças aerodinâmicas, elásticas e de inércia. Tal ciência é de grande importância no campo aeronáutico uma vez que as estruturas alares são flexíveis, devem suportar os esforços aerodinâmicos e serem rígidas o suficiente para garantir que esteja livre de todos os problemas aeroelásticos característicos (divergência, eficiência de controle, flutter e buffeting) dentro da faixa operacional de velocidades desenvolvida pela aeronave. Realizou-se uma análise modal da estrutura a fim de se conhecer os modos naturais de vibração e as respectivas frequências naturais. Para tal, utilizou-se o ANSYS Structural e o método dos elementos finitos, além de um estudo de malha para verificar a convergência dos resultados. Estudou-se também a influência da posição do lastro na ponta da placa plana, que causa a diminuição da segunda frequência natural. Além disso, realizou-se uma análise bidimensional de um volume de controle do tipo C-Grid, uma vez que o tamanho do volume de controle em uma análise aerodinâmica computacional é um fator extremamente importante. Com um volume de controle grande, tem se mais elementos na malha, caso o mesmo seja pequeno, as condições de contorno juntamente com os tamanhos dos elementos podem interferir nos resultados dos campos de velocidade e pressão em torno da estrutura. Nesse contexto, utilizou-se do software ANSYS Fluent para a simulação aerodinâmica da placa plana inclinada e obtenção dos coeficientes aerodinâmicos de sustentação e arrasto CL e CD. Os resultados foram comparados com resultados experimentais em túnel de vento de Goudeseune (SELIG; ROBERT; WILLIAMSON, 2011). Através do cálculo do Grid Convergence Index (GCI) e da comparação dos resultados numéricos com os dados experimentais constatou-se a convergência e conseguiu-se determinar um tamanho de volume de controle com erro baixo e aceitável. A análise fluido-estrutura acoplada de duas vias foi realizada com o ANSYS Structural para analisar a dinâmica estrutural através do método dos elementos finitos e o ANSYS CFX para resolver o campo do escoamento mediante método dos volumes finitos. Obtiveram-se o comportamento oscilatório da estrutura, além do coeficiente de amortecimento e tensões de von Mises. Analisando o comportamento transiente da dinâmica estrutural mediante um fluxo aerodinâmico constante (velocidade fixa). As simulações representaram bem o fenômeno, já que com o aumento da velocidade, o escoamento induz maior amortecimento à estrutura quando comparado com baixas velocidades. / The aeroelasticity is the field of science that studies the relationship between the aerodynamic elastic and inertia forces. Such knowledge is of great importance in the aviation field since the wing structures are flexible, must withstand the aerodynamic loads and be rigid enough to ensure that it is free from all aeroelastic problems like divergence, control efficiency, flutter and buffeting within the operating speed range. We carried out a modal analysis of the structure in order to know the natural vibration modes and natural frequencies. To this end, we used the ANSYS Structural with finite element method, a mesh study to verify the convergence of the results. It is also studied the influence of the slender body position of the tip of the flat plate, which causes the decrease of the second natural frequency. Furthermore, there was a twodimensional analysis of a volume control type C-Grid, since the control volume aerodynamic size in a computational analysis is an extremely important factor. A large volume of control has more elements in the mesh if it is small, the boundary conditions together with the sizes of elements may affect the results of the velocity field and pressure around the structure. In this context, we used the ANSYS FLUENT for the aerodynamic simulation of the inclined flat plate, and obtaining the aerodynamic support, and drag coefficients CL and CD. The results were compared with experimental results of Goudeseune wind tunnel (SELIG; ROBERT; WILLIAMSON, 2011). By calculating the Grid Convergence Index (GCI) and comparing the numerical results with experimental data found the convergence and managed to determine a control volume size with low and acceptable error. The fluid-structure coupled two-way analysis was performed using ANSYS Structural to analyze the structural dynamics through the finite element method and ANSYS CFX to resolve the flow field by the finite volume method. It was possible to obtain the oscillatory behavior of the structure, besides the damping coefficient and von Mises stresses. Analyzing the transient behavior of structural dynamics by a constant aerodynamic flow (fixed speed), the simulations represented the phenomenon as well, since with the increase in speed, the flow induces cushioning structure as compared to low speed

Resumo: 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 / Abstract: 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 / Orientador: Aparecido Carlos Gonçalves / Coorientador: Marcio Antonio Bazani / Banca: Gilberto Pechoto de Melo / Banca: Edson Antonio Capello Sousa / Mestre

O presente trabalho tem por objetivo desenvolver uma formulação numérica baseada em Análise Isogeométrica para o estudo de problemas de interação fluido-estrutura (IFE) em aplicações envolvendo corpos rígidos submersos, onde escoamentos incompressíveis de fluidos Newtonianos com superfície livre são considerados. Propõe-se o emprego da Análise Isogeométrica por permitir a unificação entre os procedimentos de pré-processamento e análise, melhorando assim as condições de continuidade das funções de base empregadas tanto na discretização espacial do problema como na aproximação das variáveis do sistema de equações. O sistema de equações fundamentais do escoamento é formado pelas equações de Navier-Stokes e pela equação da conservação de massa, descrita segundo a hipótese de pseudo-compressibilidade, em uma formulação cinemática ALE (Arbitrary Lagrangean- Eulerian). A consideração da superfície livre no escoamento se dá tratando o fluido como um meio bifásico, através do método Level Set. O corpo rígido apresenta não linearidade na rotação e restrições representadas por vínculos elásticos e amortecedores viscosos, sendo a equação de equilíbrio dinâmico resolvida através do método de Newmark. O esquema de acoplamento sólido-fluido adotado é o particionado convencional, que impõe condições de compatibilidade cinemáticas e de equilíbrio sobre a interface sólido-fluido, analisando ambos os meios de maneira sequencial. A discretização das equações governantes é realizada através do esquema explícito de dois passos de Taylor-Galerkin, aplicado no contexto da Análise Isogeométrica. Por fim, são analisados alguns problemas da Dinâmica de Fluidos Computacional, de onde se concluiu que os resultados obtidos são bastante consistentes com os fenômenos envolvidos, com as ferramentas exclusivas da Análise Isogeométrica, como o refinamento k, melhorando a convergência dos resultados. Para escoamentos bifásicos, verificou-se que o método Level Set obteve resultados bastante promissores apresentando, entretanto, uma dissipação numérica excessiva. Propõe-se, para estudos futuros, a elaboração de esquemas numéricos que conservem melhor o volume da fase líquida do escoamento. / The present work aims to development of a numerical formulation based on Isogeometric Analysis for the study of Fluid-Structure Interaction problems in applications involving rigid bodies submerged, considering incompressible Newtonian flows with free surface. The use of the Isogeometric Analysis allows unification between the preprocessing and analysis steps, improving then the continuity of the base functions employed, both in the spatial discretization and approximation of the variables in the system of equations. The fundamental flow equations are formed by the Navier-Stokes and the mass conservation, described by de pseudo-compressibility hypothesis, in an ALE (Arbitrary Lagrangean-Eulerian) kinematic formulation. The free surface consideration of the flow is handled treating the fluid like a two- phase medium, using the Level Set method. The rigid body considers nonlinearity in rotation, and restrictions represented by elastic springs and viscous dampers, with the dynamic equilibrium equation being resolved using the Newmark’s method. The solid-fluid coupling scheme is the conventional partitioned, which imposes kinematics and equilibrium compatibility conditions on the solid-fluid interface, analyzing both mediums in a sequential manner. The governing equations are discretized using the explicit two step Taylor-Galerkin method, applied in an Isogeometric Analisys context. Finally, some Computational Fluid Dinamics problems are analysed, from which it was concluded that the results obtained are quite consistent with phenomena involved, with the unique tools of Isogeometric Analysis, such as k-refinement, improving the convergence of the results. For biphasic flows, it was verified that the Level Set method obtained very promising results, presenting, however, an excessive numerical dissipation. For future studies, it is proposed the elaboration of numerical schemes that better preserve the volume of the liquid phase of the flow.

Une méthode complètement Eulérienne reposant sur un modèle 1-fluide est présentée afinde résoudre les problèmes d’interaction fluide-structure élastique. L’interface entre le fluideet la structure élastique est représentée par une fonction level-set, transportée par le champde vitesse du fluide et résolue par un schéma d’ordre élevé WENO 5. Les déformationsélastiques sont calculées sur la grille eulérienne à l’aide des caractéristiques rétrogrades.Nous utilisons différents modèles d’hyperélasticité, afin de générer puis d’intégrer les forcesélastiques comme terme source des équations de Navier Stokes. Le couplage vitesse/pressionest résolu par une méthode de correction de pression et les équations sont discrétisées parla méthode des volumes finis sur la grille eulérienne. La principale difficulté réside dansles grands déplacements de fluide autour du solide, source d’instabilités numériques. Afind’éviter ces problèmes, nous effectuons périodiquement une redistanciation de la level-setet une extrapolation linéaire des caractéristiques rétrogrades. Dans un premier temps,nous effectuons la vérification et la validation de notre approche à l’aide de plusieurs castests comme celui proposé par Turek. Ensuite, nous appliquons notre méthode à l’étudedu phénomène d’atténuation des vagues par des structures élastiques. Il s’agit d’une desvoies possibles pour réduire l’impact des fortes houles sur notre littoral. De plus dans lalittérature et à notre connaissance, seules des structures élastiques rigides ou élastiquesmais monodimensionnelles ont été utilisées pour réaliser ces études. Nous proposons deplacer des structures élastiques sur les fonds marins et analysons leur capacité d’absorptionde l’énergie produite par les vagues. / A fully Eulerian method is developed to solve the problem of fluid-elastic structure interactionsbased on a 1-fluid method. The interface between the fluid and the elastic structureis captured by a level set function, advected by the fluid velocity and solved with a WENO5 scheme. The elastic deformations are computed in an Eulerian framework thanks to thebackward characteristics. We use the Neo Hookean or Mooney Rivlin hyperelastic modelsand the elastic forces are incorporated as a source term in the incompressible Navier-Stokesequations. The velocity/pressure coupling is solved with a pressure-correction methodand the equations are discretized by finite volume schemes on a Cartesian grid. The maindifficulty resides in that large deformations in the fluid cause numerical instabilities. Inorder to avoid these problems, we use a re-initialization process for the level set and linearextrapolation of the backward characteristics. First, we verify and validate our approachon several test cases, including the benchmark of FSI proposed by Turek. Next, we applythis method to study the wave damping phenomenon which is a mean to reduce thewaves impact on the coastline. So far, to our knowledge, only simulations with rigid orone dimensional elastic structure has been studied in the literature. We propose to placeelastic structures on the seabed and we analyse their capacity to absorb the wave energy

O objetivo deste trabalho é estudar e desenvolver procedimentos numéricos adequados para a análise de problemas da Engenharia do Vento Computacional (EVC). O escoamento é analisado a partir das equações de Navier-Stokes para um fluido Newtoniano e de uma equação de conservação de massa considerando a hipótese de pseudo-compressibilidade, ambas em um processo isotérmico. Na presença de escoamentos turbulentos emprega-se a Simulação de Grandes Escalas (“LES”) com os modelos clássico e dinâmico de Smagorinsky para as escalas inferiores à resolução da malha. Dois modelos numéricos de Taylor-Galerkin para a análise do escoamento são estudados: o esquema explícito de dois passos e o esquema explícito-iterativo. O Método dos Elementos Finitos (MEF) é empregado para a discretização do domínio espacial utilizando o elemento hexaédrico trilinear isoparamétrico com integração reduzida das matrizes em nível de elemento. Em problemas envolvendo efeitos de interação fluido-estrutura emprega-se um esquema de acoplamento particionado com características superiores de conservação, permitindo, inclusive, o uso de subciclos entre as análises do fluido e da estrutura e de malhas não compatíveis na interface. A estrutura é considerada como um corpo deformável constituído de um material elástico linear com a presença de nãolinearidade geométrica. O MEF é também usado para a discretização da estrutura, empregando-se para tanto o elemento hexaédrico trilinear isoparamétrico com integração reduzida e controle de modos espúrios. A equação de equilíbrio dinâmico é integrada no tempo utilizando o método implícito de Newmark no contexto do método de estabilização α- Generalizado. Na presença de estruturas deformáveis, o escoamento é descrito através de uma formulação arbitrária Lagrangeana-Euleriana (ALE). Ao final, comparações com exemplos numéricos e experimentais são apresentadas para demonstrar a viabilidade dos algoritmos desenvolvidos, seguindo-se com as conclusões do trabalho e as sugestões para trabalhos futuros. / Analysis and development of numerical tools to simulate Computational Wind Engineering (CWE) problems is the main goal of the present work. The isothermal flow is analyzed using the Navier-Stokes equations for viscous fluids and a mass conservation equation obtained according to the pseudo-compressibility assumption. Turbulent flows are simulated employing Large Eddy Simulation (LES) with the classical and dynamic Smagorinsky’s models for subgrid scales. Two Taylor-Galerkin models for the flow analysis are investigated: the explicit two-step scheme and the explicit-iterative scheme. The Finite Element Method (MEF) is employed for spatial discretizations using the eight-node hexahedrical isoparametric element with one-point quadrature. Fluid-structure interaction problems are analyzed with a coupling model based on a conservative partitioned scheme. The Finite Element Method (MEF) is employed for spatial discretizations using the eight-node hexahedrical isoparametric element with one-point quadrature. Fluid-structure interaction problems are analyzed with a coupling model based on a conservative partitioned scheme. Subcycling and nonmatching meshes for independent discretizations of the fluid and structure domains are also available. The structure is considered as a deformable body constituted by a linear elastic material with geometrically nonlinear effects. The FEM is used for the spatial discretization of the structure as well. Eight-node hexahedrical isoparametric elements with one-point quadrature and hourglass control are adopted in this process. The implicit Newmark algorithm within the framework of the α-Generalized method is employed for the numerical integration of the dynamic equilibrium equation. An arbitrary Lagrangean-Eulerian (ALE) description is adopted for the kinematic description of the flow when deformable structures are analyzed. Numerical and experimental examples are simulated in order to demonstrate the accuracy of the developed algorithms. Concluding remarks and suggestions for future works are pointed out in the last chapter of the present work.

Este trabalho pode ser dividido em três etapas principais. Inicialmente é proposta uma formulação estabilizada do método dos elementos finitos (MEF) para solução de problemas de escoamento incompressível governado pela equação de Navier-Stokes. Esta formulação foi implementada em um código computacional e testada através de diversos exemplos numéricos. Alguns elementos finitos com diferentes pares de função de interpolação da velocidade e pressão, consagrados na literatura, e também elementos finitos menos populares, foram investigados e seus resultados e performance comparados. A segunda etapa consiste na formulação do problema estrutural. Buscou-se por uma formulação dinâmica, não linear, capaz de simular movimentos complexos de estruturas sujeitas a grandes deslocamentos e grandes deformações durante longos intervalos de tempo. A etapa final deste trabalho é a proposição de um método para solução de problemas de Interação Fluido Estrutura (IFE) que utiliza o conceito de fronteiras imersas como alternativa a abordagens ALE (Arbitrary Lagrangian Eulerian) clássicas. Elementos Finitos Generalizados, juntamente com Multiplicadores de Lagrange, são utilizados para prover descontinuidade nos campos de velocidade e pressão do fluido ao longo da interface com a estrutura. O acoplamento dos dois problemas é realizado utilizando um método implícito e alternado (staggered scheme), que possui a vantagem de permitir, facilmente, a implementação de códigos computacionais desenvolvidos para resolver isoladamente o problema fluido e/ou estrutural. / This work is divided in three parts. Initially, it is presented a stabilized Finite Element Method formulation to solve fluid flow problems governed by the incompressible Navier-Stokes Equations. This formulation was implemented in a computer code and validated throughout several numeric simulations. Some well-known finite elements with different pairs of velocity/pressure approximations, as well as some other less popular elements, were investigated and their performance compared. The second part describes the Structural Problem formulation. This formulation is able to simulate nonlinear dynamic problems involving large displacements and finite strains during long period of time. In the final part of this work, it is proposed a Fluid-Structure Interaction method based on an immersed interface approach in opposition to classical ALE (Arbitrary Lagrangian Eulerian) approaches. Generalized Finite Elements, together with Lagrange Multipliers, are used to provide velocity and pressure discontinuities on the fluid domain across the immersed interface. To couple both fluid and structural problems, an implicit staggered scheme is adopted, which allows the easy implementation of already developed black box computer codes.