O presente trabalho vem propor um estudo numérico qualitativo de um modelo hidro-elástico para simulação das vibrações induzidas por vórtices (VIV) em cabos submersos. O modelo é composto de um sistema acoplado de equações que descrevem tanto a estrutura quanto o fluido. A estrutura é modelada através de um oscilador elástico clássico e discretizada através de um método de massas concentradas. O fluido é considerado através de osciladores discretos baseados em modelos fenomenológicos com equação de Van der Pol ou Ginzburg-Landau. O acoplamento entre os dois ocorre na direção transversal ao escoamento local através da força de sustentação. O estudo foi feito para verificar quais tipos de comportamentos de resposta podem ser obtidos através do uso dos modelos fenomenológicos e, também, avaliar o potencial de tais modelos no ciclo de projeto na indústria atual, com foco na indústria petrolífera. Esta adequabilidade ao uso em projeto leva em conta, em primeiro lugar, a qualidade das respostas, mas também leva em conta a viabilidade da execução de simulações numéricas em tempo apropriado. Assim, um modelo numérico foi desenvolvido para permitir a integração, no domínio do tempo, deste sistema hidro- elástico. Esta integração é feita através de um método explícito de Euler e permite tratar não-linearidades dos osciladores elástico e fluido. Um conjunto de situações foi simulado numericamente, incluindo configurações verticais e também configurações em catenária. As simulações indicam alto grau de acoplamento entre os osciladores, o que leva a ressonância entre ambos na maior parte dos casos. Com relação ao comportamento do oscilador estrutural em si, foi possível observar fenômenos interessantes tais como lock-in e travelling waves, além do terceiro harmônico da vibração transversal detectado em alguns experimentos reportados em literatura. Com relação à adequabilidade do uso do modelo em projeto, considera-se que isto seja possível dada a relativa rapidez das simulações e aos resultados promissores que ainda podem ser melhorados. / The present work proposes a qualitative numerical study of a hydro-elastic model in order to simulate vortex-induced vibrations in submerged cables. The model is composed of a coupled system of equations describing the structure and the fluid. The structure is modeled through a classical elastic oscillator and discretized using a lumped mass approach. The fluid is modeled through discrete oscillators based on phenomenological models using Van der Pol or Ginzburg- Landau type equations. The coupling between these two oscillators is carried by the lift force. The study was carried in order to verify which response behavior can be obtained by using such phenomenological models and also, to evaluate their potential to be used in the offshore industry. In order to accomplish that, the model must provide reliable answers and must allow fast responses when it comes to simulation time. In this sense, a numerical model was developed to allow the time domain integration of this hydro-elastic model. The numerical integration is performed by a simple explicit Euler algorithm and allows dealing with non-linearities of both oscillators, elastic and fluid. A set of conditions was numerically simulated, including top tensioned risers and catenary. Simulations indicate a high level of coupling between the two oscillators and thus resonance is achieved in most part of the cases. With respect to the structural oscillator analysis, it was possible to observe interesting phenomena such as lock-in, travelling waves and also the third harmonic observed in some experiments and reported in the literature. Regarding the feasibility of the model usage in the daily offshore industry life, it seems possible once the model generally runs fast and some promising results were achieved.

Estudou-se a aplicação de uma modelagem RANS, com a utilização do modelo de turbulência k-? na modelagem do escoamento no interior de um sistema de lavagem de gases, composto por um lavador Venturi e uma coluna de spray, por técnicas de fluidodinâmica computacional (CFD), com a fase líquida resolvida pela abordagem Lagrangiana e a atomização secundária dessas por diferentes modelos de quebra. Utilizou-se o pacote comercial de CFD da ANSYS CFX 15.0 com as simulações realizadas em estado transiente. Verificaram-se os resultados obtidos pelo simulador para um lavador Venturi em escala piloto com dados da literatura, obtendo-se boa adequação para a perda de pressão ao longo do equipamento. Comparou-se a influência do tamanho das gotas de líquido da injeção na distribuição de tamanho de gotas ao longo do equipamento utilizando o modelo de quebra CAB (Cascade Atomization and drop Breakup), obtendo-se uma melhor dispersão das gotas de líquido na secção transversal do lavador ao se utilizar na injeção, uma distribuição de tamanho de gotas ao invés de gotas com mesmo diâmetro. Simularam-se diferentes condições de vazão de líquido e gás em um lavador Venturi em escala industrial e compararam-se os resultados de perda de pressão obtidos com correlações da literatura, obtendo-se boa adequação para as regiões iniciais do equipamento. Utilizando o modelo de quebra CAB e a abordagem Euler-Lagrange, comparou-se a distribuição de tamanho de gotas ao longo do lavador Venturi industrial para diferentes vazões de gás. Com a simulação do lavador Venturi, utilizou-se os resultados obtidos neste e simulou-se o sistema de lavagem de gases completo (coluna de Spray e lavador Venturi). Comparou-se ao final a retirada de poluentes por meio da impactação inercial e da difusão, obtendo-se uma melhor probabilidade de retirada de poluentes com o aumento da vazão de líquido. / The flow field inside a wet scrubber was simulated with a RANS model using Computational Fluid Dynamics (CFD) in which the multiphase flow was solved by an Euler-Lagrange approach. Turbulence was modeled by the k-? model and droplet breakup was assumed to occur. The transient flow equations were solved using a commercial CFD package (ANSYS CFX 15.0). Simulations of a pilot plant Venturi scrubber were compared with literature data, in which a good agreement level was achieved for pressure loss through the scrubber. Droplet size distribution was evaluated throughout the scrubber by changing droplet diameter of injected liquid and using a Cascade Atomization and drop Breakup model (CAB). A better liquid dispersion inside the scrubber was achieved for a droplet size distribution injection when compared with a single diameter droplet liquid injection. After the results verification by comparing with data from a pilot scale Venturi, pressure loss for different flow conditions was compared between numerical simulations and literature correlations in an industrial size Venturi scrubber, achieving a good agreement for the entrance section of the scrubber. Droplet size distribution throughout the industrial Venturi scrubber for different gas mass flow conditions were also compared using the CAB model and Euler-Lagrange approach. Removal of contaminants was evaluated by two parameters, being inertial impaction of particulate and diffusion to droplet surface, by which a better contaminant removal was achieved for higher liquid flow rates.

Nos últimos anos, devido à crescente preocupação com a sustentabilidade, foram despendidos mundialmente grandes esforços para a redução do consumo de energia pelos sistemas prediais. Em países tropicais como o Brasil, a ventilação natural é uma maneira efetiva e econômica para melhorar o conforto térmico dentro de habitações. Ela contribui para diminuir o uso do condicionamento de ar e renovar o ar da edificação, ajudando a reduzir as chances de se ter a Síndrome do Edifício Doente (Sick Building Syndrome SBS) e melhorando a qualidade do ar interno (Indoor Air Quality IAQ). Para tirar proveito destas vantagens da ventilação natural, o comportamento do fluxo de ar dentro da edificação deve ser analisado considerando o clima local. Existem diversos códigos computacionais baseados na Dinâmica de Fluidos Computacional (Computational Fluid Dynamics CFD) que podem ser utilizados para esta finalidade, no entanto, o CFD é um campo que requer conhecimentos altamente especializados e experiência prática para se obter bons resultados e este conhecimento geralmente está além da formação da maioria dos engenheiros e arquitetos. Com base nas dificuldades listadas e na necessidade de complementar a formação de engenheiros e arquitetos nesta área do conhecimento, este trabalho está focado em dois objetivos: o primeiro é implementar um simulador numérico computacional baseado no algoritmo Solution Algorithm for Transient Fluid Flows SOLA e as condições de contorno necessárias para a simulação da ventilação, sendo que a validação do simulador foi realizada por comparação com resultados numéricos e experimentais existentes na literatura. O segundo objetivo é propor uma ferramenta prática para a análise da ventilação natural na fase de projeto, com uma abordagem baseada na teoria de sistemas nebulosos, para identificar as melhores configurações de aberturas para um dado leiaute. Para isto, adotou-se a idéia utilizada por Givoni em seu estudo experimental: o espaço interno de uma sala é dividido em subdomínios onde a velocidade média do ar, sob diversas configurações de aberturas, é registrada. Como as velocidades médias refletem bem a eficácia da ventilação no subdomínio, elas formam a base para a definição espacial da função de pertinência para boa circulação de ar dentro da sala, considerando cada configuração de abertura. No entanto, ao invés de usar resultados experimentais, uma série de simulações computacionais baseadas em CFD, foram executadas para compor um banco de dados para avaliação das funções de pertinência. Por outro lado, temos o leiaute, que é produzido durante a concepção do projeto. Na medida em que o leiaute provê as informações para elaborar os requisitos do fluxo de ar, a função de pertinência relacionada ao fluxo de ar em cada subdomínio deve ser avaliada baseada no leiaute e nos requisitos do usuário. Acertando os requisitos providos pelo leiaute com a eficácia do fluxo de ar provido pela configuração de abertura, pode ser identificada a configuração que melhor se adapta ao leiaute. Nos casos analisados neste trabalho, o método mostrou-se promissor, indicando a configuração típica que melhor atende aos requisitos de projeto com uma boa conformidade com os resultados obtidos pela simulação da sala completa, incluindo a mobília. / In the last years, due to the concerns on sustainability, a great effort in energy saving of building systems is being carried out worldwide. In tropical countries such as Brazil, the natural ventilation is an effective and economical option to improve thermal comfort inside the dwellings, to avoid the use of costly HVAC systems, and to renew the indoor air, contributing to mitigate the Sick Building Syndrome (SBS) and to improve the Indoor Air Quality (IAQ). In order to take advantage of the natural ventilation, the behaviors of the airflow inside the buildings must be analyzed considering the local climate. There are many computer simulation codes based on Computational Fluid Dynamics (CFD) that may be adopted for this purpose, however, CFD is a field that requires highly specialized knowledge and experience to achieve good results and this expertise, which is needed to obtain reliable numerical results, is generally beyond the formation of the most part of architects and engineers. Owing to these difficulties and on the necessity to form engineers and architects in this area of knowledge, this work is focused in two main objectives: The first one is to implement a numeric computational simulation program based on Solution Algorithm for Transient Fluid Flows (SOLA) and the boundary conditions needed to simulate ventilation. The validation of the code is made by comparing the numerical results with results obtained using numerical or experimental methods published by other authors. The second objective is to propose a practical tool for the analysis of natural ventilation in the design of dwellings, with an approach based on the concepts of the Fuzzy Systems Theory to identify the best configurations of the openings for a given layout. For this, the idea used by Givonis experimental study is adopted: the inner space of a room is divided in sub-domains whose mean air velocities under different opening configurations are recorded. As the mean velocities reflect very well the effectiveness of the ventilation in the sub-domain, they form the basis for the definition of spatial distribution of the membership function for good air circulation inside the room concerning each opening configuration. However, instead of the experimental ones, a series of computer simulations were carried out to build a database for the assessment of the membership functions. On the other hand, we have the sketch of the layout, which is produced during the conceptual stage of design. As the layout provides the information about the requirements for the airflow, the membership function regarding the desirable air flow for each sub-domain might be assessed based on the layout and considering the user\'s requirements. By matching the requirements provided by the layout with the effectiveness of the airflow provided by the opening configurations, the opening configuration that best fits the layout can be identified. In the cases analyzed in this work, the method shows promising results. The typical configuration that best fits the design requirements with a good conformity with the results was obtained by full room simulation, including the furniture.

Flow machines are very important to industry, being widely used on various processes. Performance improvements are relevant factors and can be achieved by using optimization methods, such as topology optimization. Thus, this work aims to develop a method to design radial flow machine rotors operating on laminar regime, by implementing a topology optimization formulation based on density model. The design of a rotor involves firstly modelling the fluid flow by using the Navier-Stokes equations on a rotating reference frame and using the Finite Element Method for solving the differential equations. To determine the material distribution on the domain, a porous flow model based on the Darcy equation is employed by using an inverse permeability that interpolates between fluid and solid. In the optimization phase, it is defined a multi-objective function that aims to minimize the viscous energy dissipation, vorticity and power. The optimization problem is implemented using the FEniCS environment and the libraries dolfin-adjoint and pyIpopt. The optimized topologies are verified with the ANSYS software. The resulting topologies are post-processed and a CAD model is created. The rotors are manufactured by using a 3D printer, the complete prototype is built by coupling an electric brushless motor and an experimental characterization is performed by measuring fluid flow and pressure head given by the pumps. Experimental and computational results are compared and the improvement is verified. / Máquinas de fluxo são muito importantes para a indústria, sendo utilizadas em diversos processos. Assim, melhorias de desempenho são fatores relevantes e podem ser alcançadas com a utilização de métodos de otimização, como a otimização topológica. Este trabalho visa desenvolver uma metodologia para projetar rotores de máquinas de fluxo radiais que operam em escoamento laminar implementando-se a formulação de otimização topológica baseada no modelo de densidades. O projeto de rotores envolve, primeiramente, a modelagem do escoamento utilizando-se as equações de Navier-Stokes em um referencial rotativo e a utilização do Método de Elementos Finitos para a resolução das equações diferenciais. A distribuição de material no domínio é feita empregando-se um modelo de escoamento em meio poroso baseado nas equações de Darcy, utilizando-se a permeabilidade inversa que interpola o elemento entre sólido e fluido. Na fase de otimização é definida uma função multi-objetivo, que visa minimizar dissipação de energia viscosa, a vorticidade e a potência. O problema de otimização é implementado utilizando-se o ambiente FEniCS para a resolução do sistema de elementos finitos e as bibliotecas dolfin-adjoint e pyIpopt para o algorithmo de otimização. As topologias otimizadas são verificadas com o software ANSYS. As topologias resultantes são pós-processadas para a criação de um modelo CAD dos rotores. Os rotores são construídos utilizando-se a impressão 3D, o protótipo completo é montado acoplando-se um motor elétrico sem escovas e a caracterização experimental é feita medindo-se a vazão e o ganho de pressão dados pelas bombas. Por fim, os resultados experimentais e computacionais são comparados e uma melhoria de desempenho é observada.

Neste trabalho, é apresentada uma análise de estabilidade linear para um sistema flowline-riser. O modelo utiliza as equações de continuidade para as fases líquidas e gasosas e uma equação de momento para a mistura onde são considerados os efeitos de atrito. A mistura água-óleo é modelada de forma homogênea. Para a determinação da fração de vazio, adota-se o modelo de fluxo de deriva, baseado em várias correlações drift. Se o escoamento é estratificado a fração de vazio é modelada através do modelo de equilíbrio local de Taitel e Dukler (1976). Para a caracterização dos fluidos é adotado um modelo de equilíbrio de fase black-oil onde a transferência de massa é considerada entre as correntes óleo e gás em função das condições locais de pressão e temperatura. É considerada a abordagem de parâmetros distribuídos, onde diversas geometrias com discretização variável podem ser consideradas tais como riser em catenária e riser lazy wave. Para realizar a análise de estabilidade linear, as equações do modelo são linearizadas em torno do estado estacionário e discretizadas pelo método das diferenças finitas, onde foi utilizado um programa escrito em Matlab. A partir do sistema linearizado, é avaliada a estabilidade do estado estacionário pelas raízes do polinômio característico, que são solução do problema de autovalores e autovetores. É avaliada a convergência numérica do modelo e mapas de estabilidade são apresentados para vários sistemas de produção de petróleo. Os resultados obtidos numericamente são comparados com pontos operacionais de sistemas reais de produção de petróleo. Observou-se que o modelo convergiu para todos os casos avaliados e apresentou ótima concordância com os dados de campo. / This work presents a linear stability analysis for a flowline-riser system. The model considers continuity equations for the liquid and gas phases, and a momentum equation for the mixture that accounts for friction effects. The water-oil mixture is modeled as being homogeneous. The void fraction is determined by using the drift-flux model based on several drift correlations. For stratified flow, the void fraction is expressed by the local equilibrium model of Taitel and Dukler (1976). Fluid characterization uses a black-oil model that considers mass transfer between the oil and gas flows as a function of the local pressure and temperature conditions. In the proposed approach with distributed parameters, several geometries with variable discretization can be considered, such as catenary and lazy wave risers. To perform the linear stability analysis, the equations of the model are linearized around the stationary state and discretized using the finite difference method, implemented using custom-written code in Matlab. From the linearized system, the stability of the steady state is evaluated by computing the roots of the characteristic polynomial equation of the eigenvalues and eigenvectors problem. Convergence of the numerical model is evaluated and stability maps for several oil production systems were presented. Numerical results were compared with the actual measurements of oil production systems. The model converged in all cases tested and presented an excellent agreement with field data.

Os métodos para se projetar um edifício com segurança adequada para a situação de incêndio são pouco estudados e aplicados no Brasil. Essa falta de estudo é ainda mais visível quando nos deparamos com o assunto compartimentação, tópico pouco desenvolvido e ainda não normatizado pela ABNT. A compartimentação, especialmente a vertical, é fundamental para a segurança à vida, pois minimiza a propagação do incêndio entre pavimentos de edifícios. Para validar os procedimentos normatizados de dimensionamento de estruturas em situação de incêndio é necessário considerar que as edificações são compartimentadas. Atualmente as exigências sobre compartimentação são apresentadas apenas em Instruções Técnicas do Corpo de Bombeiros, variando de estado para estado. Neste trabalho é apresentado um estudo sobre a compartimentação, uma análise de sua influência na segurança contra incêndio das edificações e a comparação da visão do assunto em países de quatro continentes distintos. Após um estudo comparativo entre a visão e normatização do Brasil, Portugal, Inglaterra, Hong Kong e Estados Unidos foram realizadas simulações computacionais com base na teoria dos volumes finitos e da fluidodinâmica computacional de forma a se tentar verificar de forma mais aprofundada a influência das exigências de compartimentação na real propagação de um incêndio. Por fim, o principal foco deste trabalho é demonstrar que é possível utilizar a Dinâmica dos Fluidos Computacionais (CFD) de maneira a auxiliar o estudo do incêndio, e possibilitar uma análise mais aprofundada do assunto. / Methods for designing a building with adequate fire safety are poorly studied and applied in Brazil. This lack of study is even more present in the compartmentation subject, much-undeveloped topic and not yet standardized by ABNT. The compartmentation, especially the vertical, is critical to life safety, since it minimizes the spread of fire between building floors, and to validate the standardized procedures for structures in fire situation, because an important hypothesis of design methods is the vertical compartmentation. Currently the requirements of compartmentation are only present in the Technical Instructions of the States Fire Department, varying from state to state. This paper presents a study on compartmentation, its influence on the fire safety of buildings and a comparison of the vision of some countries of four different continents on this subject. After a comparative study of the vision and standardization between countries such as Brazil, Portugal, England, Hong Kong and the United States, computer simulation based on the theory of the finite volumes and computational fluid dynamics were carried out in order to try to check thoroughly the influence of compartmentation requirements in a real fire spread. Ultimately, the main focus of this work is to demonstrate that the Computational Fluid Dynamics (CFD) can be used in order to assist the study of fire, in order to allow a more in-depth analysis of the subject.

Made available in DSpace on 2015-03-04T18:57:51Z (GMT). No. of bitstreams: 1 Thesis_Golbert_2013.pdf: 18646737 bytes, checksum: 7adc57aa9e06e2477c908c15f1f9afb8 (MD5) Previous issue date: 2013-04-29 / O objetivo deste trabalho é estudar a modelagem do escoamento de fluidos incompressíveis, visando modelar a hemodinâmica presente no sistema cardiovascular humano. Para tanto, vamos usar método de lattice Boltzmann (LBM), baseado na cinética mesoscópica, que permite simular o comportamento macro-contínuo da dinâmica de fluidos. Acoplado a este método, será usado um método de fronteira imersa para modelar as interações entre fluido e estrutura (como ocorre entre o sangue e a parede arterial). Ainda, visando a modelagem de condições fisiológicas realistas, será empregada uma abordagem de acoplamento de modelos dimensionalmente heterogêneos. Para isto serão desenvolvidos algoritmos para acoplar o LBM com um método de diferenças finitas, para a modelagem unidimensional do escoamento de sangue em vasos deformáveis. Serão detalhados diversos aspectos e desenvolvimentos teóricos dos métodos introduzidos e faremos estudos de caráter numérico, através de simulações de problemas estacionários e transientes, cujas características são similares às encontradas na modelagem do escoamento sanguíneo em artérias. Visando prover técnicas e orientações para a modelagem do escoamento sanguíneo no sistema arterial em regime fisiológico de forma acurada e computacionalmente eficiente por meio do uso do LBM.

Uma metodologia para Determinação de Tração em Vôo é apresentada. Trata-se de uma variação do Método do Erro Residual, um dos métodos de Análise do Escoamento/Bocal apresentados no documento SAE AIR 1703A. Na medida em que o valor de tração não é uma medida direta, o método procura minimizar os desvios do resultado dos cálculos, que levam em considerações medições nos bocais do motor, em relação ao valor real. Para isso, são utilizados ensaios em solo, em que se comparam resultados calculados e medidos em bancada. O resultado dessa comparação é um escalar, o FPC, aplicado à razão de pressão no Fan de modo a minimizar os erros de vazão mássica e de tração líquida em relação aos valores calculados. Para cada razão de pressão no Fan (FPR) é obtido um FPC a partir da minimização de uma função que considera o erro residual de tração e vazão em solo. A metodologia apresentada propõe o uso de diferentes pesos para os desvios de tração e vazão. O procedimento é implementado através de rotinas em VISUAL BASIC, no ambiente do MICROSOFT EXCEL. Tais rotinas são usadas com dados de ensaios em vôo e em solo para dois motores de mesmo modelo usados na aviação comercial. Inicialmente, os dados de ensaios são utilizados para a determinação de dados de FPC e tração e vazão em vôo. É possível verificar a repetibilidade do método e a variação da tração líquida com parâmetros como rotação do Fan, altitude e Mach. Em seguida, analisa-se a variação dos pesos da tração e da vazão na determinação de FPC e nos valores finais de tração líquida e vazão mássica. Com esses dados, verifica-se a sensibilidade dos erros residuais para o ensaio em solo e a variação dos valores finais de tração e vazão em relação aos valores de pesos. Finalmente, são feitas considerações sobre o uso da metodologia proposta em uma campanha de ensaios em vôo.

O crescimento da exploração do petróleo e gás em águas ultraprofundas fez com que a demanda por unidades flutuantes crescesse e, com isso, aumentasse a ocorrência do fenômeno de Movimentos Induzidos por Vórtices (ou VIM, de Vortex-Induced Motions) que age em plataformas com casco cilíndrico, como é o caso das monocolunas e spars. Efeito natural desse aumento nas ocorrências foi um maior interesse pelo fenômeno de VIM, simplificadamente investigado em cilindros lisos e curtos como os utilizados nesta pesquisa. A fim de contribuir com o conhecimento dos fundamentos fluido-dinâmicos nesta área, foram investigados quatro diferentes cilindros flutuantes de baixa razão de aspecto, todos caracterizados pelo valor típico de L/D = 2 (comprimento imerso por diâmetro, onde L= 250 mm e D= 125 mm), mas com diferentes condições de arredondamento da extremidade imersa, estas caracterizadas por quatro razões distintas entre o raio de adoçamento e o raio do cilindro, especificamente r/R= 0,00; 0,25; 0,50 e 1,00. Com este objetivo de observar a resposta dos cilindros curtos sob efeito dos vórtices gerados em diferentes formatos de extremidade imersa, experimentos foram realizados no tanque de provas do Instituto de Pesquisas Tecnológicas do Estado de S~ao Paulo (IPT), compreendendo velocidades de reboque entre 0,024 a 0,154 m/s, ou seja, uma faixa de números de Reynolds entre 3.300 e 19.200. Os resultados de amplitude adimensional de resposta observados na direção transversal mostram influência moderada das modificações na extremidade imersa, ligeiramente mais acentuadas no modelo com razão de arredondamento de r/R=0,25, sendo que para este caso, as amplitudes máximas ficaram abaixo daquelas exibidas pelos demais modelos, com uma tendência geral de queda para as velocidades de reboque mais altas. A variação da geometria da extremidade imersa também influenciou moderadamente os resultados de amplitude adimensional de resposta na direção longitudinal, onde o modelo com r/R=1,00 apresentou menores amplitudes nas velocidades mais altas. Por outro lado, a razão entre as frequências de resposta e as frequências naturais de referência mostraram comportamentos típicos, semelhantes aos encontrados em análises de VIV de cilindros longos montados em suportes elásticos com dois graus de liberdade. Além disso, a observação dos movimentos dos modelos no plano da superfície livre não indicaram mudanças significativas nas trajetórias em função do arredondamento da extremidade livre. No tocante aos coeficientes de força (arrasto e sustentação), as comparações entre os modelos foram semelhantes, mas os coeficientes para o cilindro sem chanfro (r/R=0) mostraram-se maiores para grande parte das velocidades ensaiadas. Para o modelo com r/R=1,00, as força (arrasto e sustentação não se mostraram com a mesma intensidade do que aquelas nos demais modelos, o que, desta forma, se refletiu nos menores valores de coeficiente encontrados. Sob uma ótica geral, os resultados levam a crer que o modelo fluido responsável pelas oscilações nos cilindros com L/D = 2 não tem relação com a existência de vórtices de ponta e aresta, mas com a emissão de vórtices em forma de arco. / The oil and gas exploration in deep and ultra-deep water has increased in the past years, and at the same time, increasing the oating vessel demand, thereby the vortex-induced motion phenomenon (VIM) has been present for cylindrical hull structures, for example, spar and monocolumns oating units. The impact of this, enhance the VIM phenomenon interest with simplify investigations on smooth cylinders, as was used in the present research. In order to expand the knowledge, experiments were made in four floating cylinders with low aspect of ratio, L/D = 2 (Length / Diameter) were tested with different free end corner shape types, namely by the relation between chamfer rounding radius (r) divided by the radius of cylinder (R) (r/R = 0.0, 0.25, 0.5 and 1.0). For the initial case, r/R =0.0 represents at tip and r/R = 1.0 the semi-hemispherical tip. The aims were to understand the effect of different free-end types on VIV behaviour of cylinders. The oating circular cylinders, i.e. unit mass ratio m* = 1(structural mass/displaced uid mass), were elastically supported by a set of linear springs to provide low structural damping on the system and allow six degrees of freedom. The aim is to understand the effects of the free end vortex in the cylinder movements, varying the free end shape. The experiments were carried out in a towing tank of IPT - Institute for Technological Research of the state of São Paulo, Brazil. The range of velocities tested is 0.024 to 0.154 m/s , and Reynolds number covered 3;300 <= Re <= 19;200. The nondimensional amplitudes results for transversal direction, show a moderated effect when the free end shape has changed, more significant in the r/R=0.25 case, where the maximum amplitudes present lower values compare with the others cylinders. For in-line nondimensional amplitudes, model with r/R= 1.00 shown lower values for high velocities. On the other hand, the frequencies ratios between in-line/transverse and natural sistem frequency, presented typical VIV response for long cylinders with two degree of freedom. In addition, the trajectory observations in the XY plane have not illustrated differences for cylinders with rounded edges. Force coeffcients outcome (drag and lift), the no-rounded edge case (r/R=0.0), presented higher coeffcients for the large proportion velocities. Furthermore, the model with a semi-spherical shape in the free end (r/R=1.00), demonstrated the both force coeffcients have not acted with the same severity as was observed for the others cylinders, and thus, the r/R=1.00 models has presented lowest values. In general, the present results suggest the answer for oscillation on L/D=2 models is caused by arc-type vortex shedding, not by the free end vortex effects.

O estudo do comportamento de fluidos é um antigo domínio das ciências da natureza. Ultimamente, fenômenos de engenharia que eram estudados empiricamente passaram a ser estudados com auxílio computacional. A Dinâmica de Fluidos Computacional (DFC) é a área da ciência da computação que estuda métodos computacionais para simulação de escoamento de fluidos, e muitas vezes é a forma mais prática, ou a única, de se observar fenômenos de interesse no escoamento. Este projeto de Mestrado procurou investigar, no âmbito da simulação de um escoamento bifásico, métodos computacionais para representar a interface entre dois fluidos imiscíveis. A separação dos fluidos por meio de uma interface é necessária para assegurar que, propriedades como viscosidade e densidade, específicas de cada fluido, sejam utilizadas corretamente para o cálculo do movimento de seus respectivos fluidos. Desenvolvemos um método lagrangeano sem a utilização de malhas com o objetivo de suprir algumas restrições de trabalhos prévios. Para representar a interface entre os dois fluidos, este método utiliza uma técnica de reconstrução de superfícies baseada em aproximações de superfícies algébricas de alta ordem. Os resultados numéricos reportados neste documento evidenciam o potencial da nossa abordagem / The study of the behaviour of fluids is an ancient field in natural sciences. Recently, engineering phenomena that were empirically studied started to be done with computacional aid. The Computational Fluid Dynamics (CFD) is the area of science that studies computational methods for computer simulation of fluid flow, and often is the most practical way, or the only, to observe phenomena of interest in flow. This Masters degree project sought to investigate, in the context of the simulation of biphasic flows, computational methods to represent the interface between two immiscible fluids. The separation of fluids by the means of an interface is required to ensure that, during the simulation, the physical properties of a fluid, like density and viscosity (specific of each fluid) are properly used in the calculus of the respective fluid motion. We developed a lagrangean method without the use of mesh with the goal of alleviating some of the previous works restrictions. To represent the interface between the two fluids, this method uses a surface reconstruction technique based on approximations of high order algebraic surfaces. The numerical results reported herein show the potential of our approach